数值分析简明教程第二版课后习题答案(供参考)
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
0.1算法
1、 (p.11,题1)用二分法求方程013
=--x x 在[1,2]内的近似根,要求误差不
超过10-3.
【解】 由二分法的误差估计式31
1*102
1
2||-++=≤=-≤
-εk k k a b x x ,得到100021≥+k .两端取自然对数得96.812
ln 10
ln 3≈-≥
k ,因此取9=k ,即至少需
2、(p.11,题2) 证明方程210)(-+=x e x f x
在区间[0,1]内有唯一个实根;使用
二分法求这一实根,要求误差不超过2102
1
-⨯。
【解】 由于210)(-+=x e x f x
,则)(x f 在区间[0,1]上连续,且
012010)0(0<-=-⨯+=e f ,082110)1(1>+=-⨯+=e e f ,即0)1()0(<⋅f f ,由连续函数的介值定理知,)(x f 在区间[0,1]上至少有一个零点.
又010)('>+=x e x f ,即)(x f 在区间[0,1]上是单调的,故)(x f 在区间[0,1]内有唯一实根.
由二分法的误差估计式211*1021
2
12||-++⨯=≤=-≤-εk k k a b x x ,得到1002≥k .
两端取自然对数得6438.63219.322
ln 10
ln 2=⨯≈≥k ,因此取7=k ,即至少需二分
0.2误差
1.(p.12,题8)已知e=2.71828…,试问其近似值7.21=x ,71.22=x ,x 2=2.71,718.23=x 各有几位有效数字?并给出它们的相对误差限。
【解】有效数字:
因为111021
05.001828.0||-⨯=
<=- x e ,所以7.21=x 有两位有效数字; 因为1
2102105.000828.0||-⨯=<=- x e ,所以71.22=x 亦有两位有效数字;
因为3
3102
10005.000028.0||-⨯=<=- x e ,所以718.23=x 有四位有效数字;
%85.17.205
.0||111=<-=
x x e r ε; %85.171.205
.0||222=<-=
x x e r ε; %0184.0718
.20005
.0||333=<-=
x x e r ε。 评 (1)经四舍五入得到的近似数,其所有数字均为有效数字;
(2)近似数的所有数字并非都是有效数字.2.(p.12,题9)设72.21=x ,
71828.22=x ,0718.03=x 均为经过四舍五入得出的近似值,试指明它们的绝对误差(限)
与相对误差(限)。
【解】 005.01=ε,31
1
11084.172.2005
.0-⨯≈<
=
x r εε; 000005.02=ε,622
21084.171828
.2000005
.0-⨯≈<
=x r εε;
00005.03=ε,43
3
31096.60718
.000005
.0-⨯≈<
=
x r εε;
评 经四舍五入得到的近似数,其绝对误差限为其末位数字所在位的半个单位.
3.(p.12,题10)已知42.11=x ,0184.02-=x ,4
310184-⨯=x 的绝对误差限均为
2105.0-⨯,问它们各有几位有效数字?
【解】 由绝对误差限均为2105.0-⨯知有效数字应从小数点后两位算起,故42.11=x ,有
三位;0184.02-=x 有一位;而0184.0101844
3=⨯=-x ,也是有一位。
1.1泰勒插值和拉格朗日插值
1、(p.54,习题1)求作x x f sin )(=在节点00=x 的5次泰勒插值多项式)(5x p ,并计算
)3367.0(5p 和估计插值误差,最后将)5.0(5p 有效数值与精确解进行比较。
【解】由x x f sin )(=,求得x x f
cos )()
1(=;x x f sin )()2(-=;x x f
cos )()
3(-=;
x x f sin )()4(=;x x f
cos )()
5(=;x x f sin )()6(-=,所以
插值误差:)(5x R 66060)6(!
61
)(!6|)sin(|)(!6|)(|x x x x x f ≤-=-=
ξξ,若5.0=x ,则 )3367.0(5p 3303742887.0!
53367.0!33367.03367.05
3≈+-=,而
566
5105.01002.2!
63367.0)3367.0(--⨯<⨯≈≈R ,精度到小数点后5位,
故取33037.0)3367.0(5=p ,与精确值 330374191.0)3367.0sin()3367.0(==f 相比
较,在插值误差的精度内完全吻合!
2、(p.55,题12)给定节点4,3,1,13210===-=x x x x ,试分别对下列函数导出拉格朗日余项:
(1)234)(3
+-=x x x f ; (2)3
4
2)(x x x f -=
【解】依题意,3=n ,拉格朗日余项公式为 ∏=-=3
)
4(3)(!4)()(i i x x f
x R ξ (1)0)()
4(=x f
→ 0)(3=x R ;
(2)因为!4)()
4(=x f ,所以
3、(p.55,题13)依据下列数据表,试用线性插值和抛物线插值分别计算)3367.0sin(的近