数值分析第五版课后答案
- 格式:doc
- 大小:2.22 MB
- 文档页数:35
第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
第一章 绪论1.设,的相对误差为,求的误差。
0x >x δln x 解:近似值的相对误差为*x *****r e x x e x x δ-===而的误差为ln x ()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设的相对误差为2%,求的相对误差。
x n x 解:设,则函数的条件数为()n f x x ='()||()p xf x C f x =又, 1'()n f x nx -= 1||n p x nx C n n-⋅∴==又((*))(*)r p r x n C x εε≈⋅ 且为2(*)r e x ((*))0.02n r x nε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:,, , ,*1 1.1021x =*20.031x =*3385.6x =*456.430x =*57 1.0.x =⨯解:是五位有效数字;*1 1.1021x =是二位有效数字;*20.031x =是四位有效数字;*3385.6x =是五位有效数字;*456.430x =是二位有效数字。
*57 1.0.x =⨯4.利用公式(2.3)求下列各近似值的误差限:(1) ,(2) ,(3) .***124x x x ++***123x x x **24/x x 其中均为第3题所给的数。
****1234,,,x x x x 解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x εεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少?解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===A A (*)(*)3(*)r p r r V C R R εεε∴≈=A 又%1(*)1r V ε=故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设,按递推公式 (n=1,2,…)028Y =1n n Y Y -=-计算到(5位有效数字),试问计算将有多大误差?100Y 27.982≈100Y解: 1n n Y Y -=10099Y Y ∴=9998Y Y =9897Y Y =-……10Y Y =-依次代入后,有1000100Y Y =-即,1000Y Y =-, 27.982≈100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯的误差限为。
数值分析第五版课后答案2篇数值分析第五版课后答案(一)第一章1.1 机器精度的数值为2^-52 ≈2.22 × 10^-16。
1.2 Example 1.2设f(x) = (1 - cosx)/sinx,则f(0)的分母为0,无法进行数值计算。
1.3 Example 1.3设f(x) = (1 - cosx)/sinx,则f(0)的分子为0,因此有f(0) = 0。
1.4 Example 1.4(a) 将x的值从1.8改为1.799,则f(x)的值由-0.000000000000159为0.000000000000313,差值为0.000000000000472。
(b) 我们有f'(x) = sinx/(1 - cosx) - 1/sin^2x。
将x的值从1.8改为1.799,利用f(x)和f'(x)的值可以得到下面的近似式:f(x + Δx) ≈ f(x) + f'(x)Δx = -0.000000000000159 + 0.449787416887455×0.001 = -0.000000000000137。
与(a)中的结果相近。
1.5 Example 1.5(a) 当x很接近于0时,函数值的符号取决于cosx的符号,其中cosx接近于1。
因此,函数值为正。
(b) 当x很接近于π时,函数值的大小趋于无穷大,因为分母趋向于0,而分子不为0。
1.6 Example 1.6(a) 因为函数在x = 0处是奇函数,所以它的导数为偶函数。
(b) 首先,我们有f''(0) = -2,因此x = 0是最大值。
其次,我们有f''(x) = -2 - 8sin^2x。
由于-f''(x)在x = 0处是正的,我们有当x越接近0时,f''(x)越小,也就意味着函数在x = 0处是严格的最大值。
1.7 Example 1.7(a) 我们有f(x) = x^3 - 2x^2 - 5x + 6,f'(x) =3x^2 - 4x - 5和f''(x) = 6x - 4。
第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
数值分析课程第五版课后习题答案课后习题一:a) 求解非线性方程f(x) = x^3 - 2x - 5的根。
解答:可使用牛顿迭代法来求解非线性方程的根。
牛顿迭代法的迭代公式为:x_(n+1) = x_n - f(x_n)/f'(x_n),其中x_n为第n次迭代的近似解。
对于给定的方程f(x) = x^3 - 2x - 5,计算f'(x)的导数为f'(x) = 3x^2 - 2。
选择一个初始近似解x_0,并进行迭代。
迭代的终止条件可以选择两次迭代间的解的差值小于某个预设的精度。
b) 计算矩阵加法和乘法的运算结果。
解答:设A和B为两个矩阵,A = [a_ij],B = [b_ij],则A和B的加法定义为C = A + B,其中C的元素为c_ij = a_ij + b_ij。
矩阵乘法定义为C = A * B,其中C的元素为c_ij = ∑(a_ik * b_kj),k的取值范围为1到矩阵的列数。
c) 使用插值方法求解函数的近似值。
解答:插值方法可用于求解函数在一组给定点处的近似值。
其中,拉格朗日插值法是一种常用的方法。
对于给定的函数f(x)和一组插值节点x_i,i的取值范围为1到n,利用拉格朗日插值多项式可以构建近似函数P(x),P(x) = ∑(f(x_i) * l_i(x)),其中l_i(x)为拉格朗日基函数,具体表达式为l_i(x) = ∏(x - x_j)/(x_i - x_j),j的取值范围为1到n并且j ≠ i。
课后习题二:a) 解决数值积分问题。
解答:数值积分是求解定积分的数值近似值的方法。
常用的数值积分方法包括矩形法、梯形法和辛普森法。
矩形法采用矩形面积的和来近似曲边梯形的面积,梯形法采用等距离子区间上梯形面积的和来近似曲边梯形的面积,而辛普森法则利用等距离子区间上梯形和抛物线面积的加权和来近似曲边梯形的面积。
b) 使用迭代方法求解线性方程组。
解答:线性方程组的求解可以通过迭代方法来进行。
第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。
第一章 绪论1.设0x >,x 的相对误差为δ,求ln x 的误差。
解:近似值*x 的相对误差为*****r e x xe x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈进而有(ln *)x εδ≈2.设x 的相对误差为2%,求n x 的相对误差。
解:设()nf x x =,则函数的条件数为'()||()p xf x C f x = 又1'()n f x nx-=, 1||n p x nx C n n-⋅∴== 又((*))(*)r p r x n C x εε≈⋅且(*)r e x 为2((*))0.02n r x n ε∴≈3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*57 1.0.x =⨯解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =⨯是二位有效数字。
4.利用公式(2.3)求下列各近似值的误差限:(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x .其中****1234,,,x x x x 均为第3题所给的数。
解:*41*32*13*34*151()1021()1021()1021()1021()102x x x x x εεεεε-----=⨯=⨯=⨯=⨯=⨯***124***1244333(1)()()()()1111010102221.0510x x x x x x εεεε----++=++=⨯+⨯+⨯=⨯ ***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x εεεε---=++=⨯⨯⨯+⨯⨯⨯+⨯⨯⨯≈**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x xεεε---+≈⨯⨯+⨯⨯=⨯=5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π=则何种函数的条件数为23'4343p R V R R C V R ππ===(*)(*)3(*)r p r r V C R R εεε∴≈=又(*)1r V ε=%1故度量半径R 时允许的相对误差限为εr (V ∗)=13∗1%=13006.设028Y =,按递推公式1n n Y Y -= (n=1,2,…)计算到100Y 27.982≈(5位有效数字),试问计算100Y 将有多大误差?解:1n n Y Y -=-10099Y Y ∴=9998Y Y =9897Y Y =……10Y Y =依次代入后,有1000100Y Y =-即1000Y Y =27.982≈, 100027.982Y Y ∴=-*310001()()(27.982)102Y Y εεε-∴=+=⨯100Y ∴的误差限为31102-⨯。