计算方法引论课后答案
计算方法_课后习题答案
(4.5)(0.01172)
0.00879
（2）采用 Newton 插值多项式 y x N2(x) 根据题意作差商表：
i
xi
0
4
1
6.25
f (xi ) 2 2.5
一阶差商 2 9
2
9
3
2 11
二阶差商 4 495
N2 (7) 2 29 (7 4) ( 4 495) (7 4) (7 6.25) 2.6484848
1
e2
则根据二次Lagrange插值公式得：
L2 (x)
(x ( x0
x1)(x x2 ) x1)(x0 x2 )
y0
(x ( x1
x0 )(x x2 ) x0 )(x1 x2 )
y1
(x ( x2
x0 )(x x1) x0 )(x2 x1)
y2
2(x 1)(x 0.5) 2x(x 0.5)e1 4x(x 1)e0.5
8. 求作 f x xn1 关于节点 xi i 0,1, , n 的 Lagrange 插值多项式，并利用
插值余项定理证明
n
n
xin1li 0 1n xi
i0
i0
式中 li x 为关于节点 xi i 0,1, , n 的 Lagrange 插值基函数。
2 02 12 4 23 4 04 14 2 3
1 x2 3x 2 x 4 3x x2 6x 8 23 x x2 5x 4 1 x x2 3x 2
8
4
8
《计算方法引论》-徐翠微主编
《计算方法引论》-徐翠微主编2009 ~ 2010学年第一学期计算方法教案计0701-0703 4h第二章插值法知识点:拉格朗日插值法，牛顿插值法，余项，分段插值。
实际问题中，时常不能给出f(x)的解析表达式或f(x)解析表达式过于复杂而难于计算，能采集的只是一些f(x)的离散点值{xi,f(xi)}(i=0,1,2,…n)。
因之，考虑近似方法成为自然之选。
定义:设f(x)为定义在区间[a，b]上的函数，x0,x1,…,xn为[a，b]上的互异点，yi=f(xi)。
若存在一个简单函数,(x)，满足(插值条件),(xi)=f(xi)，i=0，1，…，n。
则称 ,(x)为f(x)插值函数，f(x)为被插函数,点x0,x1,…,xn为插值节点,点{xi,f(xi)},i=0,1,2,…n为插值点。
于是计算f(x)的问题就转换为计算 ,(x)。
构造插值函数需要解决:插值函数是否存在唯一;插值函数如何构造(L插值);插值函数与被插函数的误差估计和收敛性。
对插值函数 ,(x)类型有多种不同的选择，代数多项式常被选作插值函数。
P23(2.18)和(2.19)指出，存在唯一的满足插值条件的n次插值多项式p(x)。
但是需要计算范德蒙行列式，构造插值多n项式工作量过大，简单表达式不易得到，实际中不采用这类方法。
p(x)?f(x) n插值法是一种古老的数学方法，拉格朗日(Lagrange)、牛顿(Newton)等分别给出了不同的解决方法。
拉格朗日插值拉格朗日(Lagrange)插值的基本思想:把插值多项式p(x)的构造问题转化为n+1个插值基函数l(x)(i=0,1,…,n)的ni构造。
(1)线性插值?构造插值函数已知函数y=f(x)的两个插值点(x,y)，(x,y)，构造多项式y=p(x)，使p(x)=y,p(x)=y。
001111001111 《计算方法引论》、徐翠薇，高等教育出版社 2008年4月第三版第二章Lagrange插值法2009 ~ 2010学年第一学期计算方法教案计0701-0703 4h由直线两点式可知，通过A，B的直线方程为, y y 1 0 , , , y y ，， x x p ( x ) + 0 0 1 , x x 1 0变形为 x-x0 x-x1 y 1, ， p(x) y 10 x1-x0 x0-x1记 x-x0 x-x1 , l(x) , l(x) 10 x1-x0 x0-x1则p(x)=l(x)y+l(x)y10011插值完毕～注意性质:l(x)=l(x)=1，l(x)=l(x)=0，p(x)=y，p(x)=y。
计算方法各章习题及答案
第二章数值分析4^92.1 已知多项式通过下列点:1 3答案：q(x) = p(x) -r(x) X5X4X3-3X 1 .2 22.2观测得到二次多项式2的值：表中p2(x)的某一个函数值有错误，试找出并校正它.答案：函数值表中p2(-1)错误，应有p2(-1) = 0 .2.3利用差分的性质证明12■ 22■川，n2=n(n ■ 1)(2n ■ 1)/6.2.4当用等距节点的分段二次插值多项式在区间[-1,1]近似函数e x时，使用多少个节点能够保证误差不超过丄10-6.2答案：需要143个插值节点.2.5 设被插值函数f (x) • C4[a,b] , H3h)(x)是f (x)关于等距节点b — aa ^Xo :::捲：：：川：::x n=b的分段三次艾尔米特插值多项式，步长h .试估计n ||f(x)-H3h)(x)||：：.答案：||住)-出5)仪川：乞令人4.384第三章函数逼近3.1求f(x)二sinx,x，[0,0.1]在空间门=span{1,x, x2}上最佳平方逼近多项式，并给出平方误差.答案：f (x) =sin X的二次最佳平方逼近多项式为sin x p2(x) = -0.832 440 7 105 1.000 999 1x - 0.024 985 1x2,二次最佳平方逼近的平方误差为20.12 12■ = 0 (sinx) - P 2(x))2dx =0.989 310 7 10•3.2确定参数a,b 和c ,使得积分1 ---------------------------2 1 I (a,b,c)[ax 2 bx c -1 -x 2]dx 取最小值.J 1 — x 2810答案：a, b = 0, c =3 二3 二3.3 求多项式f (x) =2x 4 x 3 5x 2 1在[-1,1]上的3次最佳一致逼近多项式p(x)-答案：f (x)的最佳一致逼近多项式为p(x) = X ’ 7x2 3.43.4用幕级数缩合方法，求 f(x)=e x (―1兰XW1)上的3次近似多项式 p 6,3(x)，并估计 || f(X )-P 6,3(X )II ：：.答案：p5,3(x) =0.994 574 65 + 0.997 395 83x+0.542 968 75x 2 十 0.177 083 33x 3, || f (x) - p 6,3 (x) |^<0.006 572 327 71 一3.5 求f (x) -e x ( -1乞x 乞1)上的关于权函数「(X )-的三次最佳平方逼近小-x 2多项式 Q(x)，并估计误差 || f(x)-$(x)||2 和 || f(x)-S 3(x) ||：：.答案：§3(x) =0.994 571 0.997 308x 0.542 991x 20.177 347x 3,|| f (x) -S 3(x) ||2 = 0.006 894 83, || f (x) - §3(x)||严 0.006 442 575.第四章数值积分与数值微分14.1用梯形公式、辛浦生公式和柯特斯公式分别计算积分 X n dx (n -1,2,3,4)，并与精确值比较.答案：计算结果如下表所示I 2 0. 5 0. 333 333 0. 250 000 0. 208 333 I 30. 5 0. 333 333 0. 250 000 0. 200 000 精确值0. 50. 333 3330. 250 0000. 200 0004.2 确定下列求积公式中的待定参数，使得求积公式的代数精度尽量高，并指明所确定的求积公式具有的代数精度.h(])仁 f (x)dx 止 A_i f (-h) + A f (0) + A f (h)11 (2)J(x)dx： 3【f(-1) 2f(X i ) 3f(X 2)]hh2⑴ of(x)dxVf(O) f(h)「h[f g f(h)]答案：(1)具有三次代数精确度 (2)具有二次代数精确度 (3)具有三次代数精确度. 4.3 设h = % - X 0，确定求积公式r (x - x o ) f (x)dx = h 1 2[ Af (x o ) + Bf (x i )] + h 3[C 「(x o ) + Df^)] + R[ f ]xo中的待定参数 A, B,C, D ，使得该求积公式的代数精确度尽量高，并给出余项表达式.37 1 if 4)(叮)6答案：A = —, B— ,C —, D — , R[f]=— _) h ，其中 (x o ,xi).202030 20 14404.4设P 2(x)是以0,h,2h 为插值点的f(x)的二次插值多项式，用F 2(x)导出计算积分3h3 4 5If (x)dx 的数值积分公式I h ，并用台劳展开法证明：I - l h h f (0) O(h ).力83h3答案：I h P 2(x)dx h[ f(0) 3f (2h)].0 4(3)取7个节点处的函数值.1sin x4.6用变步长的复化梯形公式和变步长的复化辛浦生公式计算积分Idx .要x1o 1«求用事后误差估计法时，截断误不超过10和 10 .1(1) 运用复化梯形公式计算上述积分值，使其截断误差不超过丄10」. 2(2) 取同样的求积节点，改用复化辛浦生公式计算时，截断误差是多少？(3) 要求的截断误差不超过10“ ,若用复化辛浦生公式，应取多少个节点处的函数值? 答案：(1)只需n — 7.5，取9个节点，I : 0.9464.5 给定积分I 二1sin xdx|R n [f]耳一孟宀皿盂日中0.271估2 2答案：使用复化梯形公式时，I T^ 0.946满足精度要求；使用复化辛浦生公式时，I s4 =0.946 083满足精度要求.4.7 ( 1 )利用埃尔米特插值公式推导带有导数值的求积公式1 323 1 3>5.2用矩阵的直接三角分解法解方程组 广1 0 2 0、「5、0 10 1 X 2312 4 3X 3仃10 1 0 3丿 g<7；答案： &=2 , x 3 = 2 , x 2 = 1, X| = 1 .ba f(x)dx 二 其中余项为b —a(b 「a)2[f(a)f(b)] — ' 丿[f (b)-f (a)] R[f]， 2 12R[f]=U 54!30 f ( 4()),(a,b).其中(2)利用上述公式推导带修正项的复化梯形求积公式h 2 f(x)dx :T^—[ f (X N ) - f (x 。
计算方法-刘师少版第三章课后习题完整答案
⎡2 1 − 1⎤ ⎡ 1 ⎢4 − 1 3 ⎥ = ⎢l ⎢ ⎥ ⎢ 21 ⎢ ⎥ − 6 9 1 ⎣ ⎦ ⎢ ⎣l 31
1 l32
⎤ ⎡u11 ⎥⎢ ⎥⎢ 1⎥ ⎦⎢ ⎣
u12 u 22
u13 ⎤ u 23 ⎥ ⎥ u 33 ⎥ ⎦
将右端两矩阵相乘后比较两端，可得
l 2 = b2 − a 2 u1 = 4 −
1 15 = 4 4
u 2 = c2 / l2 = −
由
4 , 15
l3 = b3 − a3u 2 = 4 −
4 56 = 15 15
⎤ ⎡ ⎢4 0 0 ⎥ ⎡ y1 ⎤ ⎡ 2 ⎤ ⎥ ⎢ 15 y2 ⎥ =⎢ 0 ⎥⎢ 4⎥ ⎢− 1 ⎥ ⎢ ⎥ ⎢ 4 ⎥ ⎢ y ⎥ ⎢10⎥ ⎢ 56 ⎥⎣ 3 ⎦ ⎣ ⎦ ⎢ −1 ⎢ 15 ⎥ ⎦ ⎣
−1
消元
0 .5 1 .5 2 .5 0 .5 2 .5 1 .5 0 1 0 0 1 0 0 0 1 1 3 2
1 2 3 0 0 1 0 0 1 1
0 1 0
0 .5
− 0 .5 − 0 .5
0.5
0 ⎤ 0 ⎥ ⎥ 1 ⎥ ⎦ 0 ⎤ 1 ⎥ ⎥ 0 ⎥ ⎦
列选主
− 0.5 − 0.5
6
消元
0 .4 1 .2 0 .2
⎡ ⎢ ⎢ ⎢ ⎢ ⎣
得
3 2 3 3
0 2 3 − 6
0 ⎤ ⎡ y ⎤ ⎡5⎤ ⎥ 1 0 ⎥⎢ 3⎥ y2 ⎥ =⎢ ⎢ ⎥ ⎢ ⎥ ⎥ ⎢ ⎢ ⎥ 7 y ⎦ ⎣ 3⎦ ⎣ ⎥ 3⎥ ⎦ 1 6 1 3
y1 =
再由
5 3
,
计算方法引论课后答案.doc
解:
dx n nx n1d x nd x n % . xn xn x
12. 计算球的体积,为了使相对误差限为 1%, 问度量半径 R 时允许的相对误差限如何? 解: 已知 V
4 3 dR R ,设 dr R a ,则要使得 3 R
dr V
1 dV dlnV d ln R3 3 d ln R 3 d ln R 3dr R 3a 1% ,则 a 1% . 3 V
1 1 104 , db 103 , 2 2 又 a b 0.20532 10 , 1 1 1 d a b da db da db 104 103 0.55 103 102 , 2 2 2 所以 a b 有三位有效数字; 因为 a b 0.10475714 10 , 1 1 1 d a b b da a db 0.947 104 1.1062 103 0.60045 103 102 2 2 2 所以 a b 有三位有效数字.
n
x l x x , k 0,1,..., n
j 0 k j j k
证明: 由于
1 i j; l j xi ij 0 i j.
且
x l x 和 x
j 0 k j j
n
k
都为 k 次多项式,而且在 k+1 个不同的节点处的函数值都相同
L1 x
所以
100 121
x 121
10
121 100
x 100
11
115 L1 115 10.7143 .
利用后两个节点建立线性插值可得:
L1 x
计算理论导引习题答案[第2版]CHAP5new
![计算理论导引习题答案[第2版]CHAP5new](https://img.taocdn.com/s3/m/14a23d6565ce0508773213a8.png)
5.1 证明EQ CFG 是不可判定的。
解：只须证明ALL CFG ≤m EQ CFG 即可。
构造CFG G 1，使L(G 1)=∑*。
设计从ALL CFG 到EQ CFG 的归约函数如下： F=“对于输入＜G ＞，其中G 是CFG ：1）输出＜G ,G 1＞。
”若＜G ＞ALL CFG ，则<G ,G 1>EQ CFG 。
若＜G ＞ALL CFG ，则<G , G 1>EQ CFG 。
F 将ALL CFG 归约到EQ CFG 即ALL CFG ≤m EQ CFG∵ALL CFG 是不可判定的，∴EQ CFG 是不可判定的。
5.2证明EQ CFG 是补图灵可识别的。
证明：注意到A CFG ={<G,w>|G 是能派生串w 的CFG}是可判定的。
构造如下TM ： F=“输入<G ,H>，其中G ,H 是CFG ，1) 对于字符串S 1, S 2,，重复如下步骤。
2) 检测S i 是否可以由G 和H 派生。
3) 若G 和H 中有一个能派生w ，而另一个不能，则接受。
”F 识别EQ CFG 的补。
5.3 略。
5.4 如果A m B 且B 是正则语言，这是否蕴涵着A 也是正则语言？为什么？ 解：否。
例如：对非正则语言A={0n 1n |n 0}和正则语言B={0}，可以构造一个可计算函数f 使得：f(w)=⎩⎨⎧≠=n n nn 10w 1,10w 0, 于是w A f(w)B,故A m B 。
5.5 证明A TM 不可映射规约到E TM 。
证明：反证法假设A TM m E TM ， 则有TM m TM E A ≤。
而A TM 的补不是图灵可识别的，从而可知E TM 的补也不是图灵可识别的。
下面构造一个识别E TM 的补的图灵机S ：S=“输入<M>，M 是TM,1) 对i=1,2,…重复下一步。
2) 对S 1,S 2,…,S i 模拟M 运行i 步，若有接受，则接受。
计算方法——第一章——课后习题答案_刘师少
1.1 设3.14, 3.1415, 3.1416分别作为π的近似值时所具有的有效数字位数解 近似值x =3.14＝0.314×101,即m =1，它的绝对误差是 －0.001 592 6…，有31105.06592001.0-*⨯≤=- x x . 即n =3，故x =3.14有3位有效数字. x =3.14准确到小数点后第2位. 又近似值x =3.1416，它的绝对误差是0.0000074…，有 5-1*10⨯50≤00000740=-.. x x即m =1,n ＝5，x =3.1416有5位有效数字. 而近似值x =3.1415，它的绝对误差是0.0000926…，有 4-1*10⨯50≤00009260=-.. x x即m =1,n ＝4，x =3.1415有4位有效数字.这就是说某数有s 位数，若末位数字是四舍五入得到的，那么该数有s 位有效数字1.2 指出下列各数具有几位有效数字，及其绝对误差限和相对误差限：2.0004 －0.00200 9000 9000.00解 （1）∵ 2.0004＝0.20004×101, m=1 绝对误差限：4105.0000049.020004.0-*⨯≤≤-=-x x xm -n =-4,m =1则n =5,故x =2.0004有5位有效数字1x =2，相对误差限000025.010221102151)1(1=⨯⨯=⨯⨯=---n r x ε （2）∵ －0.00200= -0.2×10-2, m =-2 5105.00000049.0)00200.0(-*⨯≤≤--=-x x xm -n =-5, m =-2则n =3,故x =-0.00200有3位有效数字1x =2，相对误差限3110221-⨯⨯=r ε=0.0025 (3) ∵ 9000=0.9000×104, m =4， 0105.049.09000⨯<≤-=-*x x xm -n =0, m =4则n =4,故x =9000有4位有效数字4110921-⨯⨯=r ε＝0.000056 (4) ∵9000.00=0.900000×104, m =4，2105.00049.000.9000-*⨯<≤-=-x x xm -n =-2, m =4则n =6,故x =9000.00有6位有效数字 相对误差限为6110921-⨯⨯=rε＝0.000 00056由(3)与(4)可以看到小数点之后的0，不是可有可无的，它是有实际意义的.10-的近似值是多少？1.3 ln2=0.69314718…，精确到310-＝0.001，即绝对误差限是ε＝0.0005,解精确到3故至少要保留小数点后三位才可以.ln2≈0.6931.1 设3.14, 3.1415, 3.1416分别作为π的近似值时所具有的有效数字位数解近似值x=3.14＝0.314×101,即m=1，它的绝对误差是－0.001 592 6…，有.即n=3，故x=3.14有3位有效数字.x=3.14准确到小数点后第2位.又近似值x=3.1416，它的绝对误差是0.0000074…，有即m=1,n＝5，x=3.1416有5位有效数字.而近似值x=3.1415，它的绝对误差是0.0000926…，有即m=1,n＝4，x=3.1415有4位有效数字.这就是说某数有s位数，若末位数字是四舍五入得到的，那么该数有s位有效数字1.2指出下列各数具有几位有效数字，及其绝对误差限和相对误差限：2.0004 －0.00200 9000 9000.00 解（1）∵ 2.0004＝0.20004×101, m=1绝对误差限：m-n=-4,m=1则n=5,故x=2.0004有5位有效数字=2，相对误差限（2）∵ －0.00200= -0.2×10-2, m=-2m-n=-5,m=-2则n=3,故x=-0.00200有3位有效数字=2，相对误差限 =0.0025(3)∵ 9000=0.9000×104, m=4，m-n=0,m=4则n=4,故x=9000有4位有效数字＝0.000056(4) ∵9000.00=0.900000×104, m=4，m-n=-2,m=4则n=6,故x=9000.00有6位有效数字相对误差限为＝0.000 00056由(3)与(4)可以看到小数点之后的0，不是可有可无的，它是有实际意义的.1.3 ln2=0.69314718…，精确到的近似值是多少？解精确到＝0.001，即绝对误差限是e＝0.0005,故至少要保留小数点后三位才可以.ln2»0.6932.1 用二分法求方程在[1, 2]的近似根，要求误差不超过至少要二分多少？解：给定误差限e＝0.5×10－3，使用二分法时，误差限为只要取k满足即可，亦即只要取n＝10.2.3 证明方程1 -x–sin x＝0 在区间[0, 1]内有一个根，使用二分法求误差不超过0.5×10-4的根要二分多少次？证明令f(x)＝1－x－sin x，∵f(0)=1>0，f(1)=－sin1<0∴f(x)=1－x－sin x=0在[0，1]有根.又f¢(x)=-1－c os x<0 (xÎ[0.1]),故f(x) 在[0，1]单调减少，所以f(x) 在区间[0，1]内有唯一实根.给定误差限e＝0.5×10－4，使用二分法时，误差限为只要取k满足即可，亦即只要取n＝14.。
计算方法引论课后答案
计算方法引论课后答案第一章误差1.什么是模型误差，什么是方法误差？例如，将地球近似看为一个标准球体，利用公式 $A=4\pi r$ 计算其表面积，这个近似看为球体的过程产生的误差即为模型误差。
在计算过程中，要用到 $\pi$，我们利用无穷乘积公式计算 $\pi$ 的值：pi=2\cdot\frac{2}{1}\cdot\frac{2}{3}\cdot\frac{4}{3}\cdot\f rac{4}{5}\cdot\frac{6}{5}\cdot\frac{6}{7}\cdot\frac{8}{7}\cdot\ frac{8}{9}\cdot\cdots我们取前9项的乘积作为 $\pi$ 的近似值，得$\pi\approx3.xxxxxxxx5$。
这个去掉 $\pi$ 的无穷乘积公式中第9项后的部分产生的误差就是方法误差，也称为截断误差。
2.按照四舍五入的原则，将下列各数舍成五位有效数字：816.956，76.000，.322，501.235，.182，130.015，236.23.解：816.96，76.000，.501.24，.130.02，236.23.3.下列各数是按照四舍五入原则得到的近似数，它们各有几位有效数字？81.897，0.008，136.320，050.180.解：五位，三位，六位，四位。
4.若 $1/4$ 用 0.25 表示，问有多少位有效数字？解：两位。
5.若 $a=1.1062$，$b=0.947$，是经过舍入后得到的近似值，问：$a+b$，$a\times b$ 各有几位有效数字？已知 $da<\frac{1}{2}\cdot10^{-4}$，$db<\frac{1}{2}\cdot10^{-3}$，又 $a+b=0.\times10$。
begin{aligned}d(a+b)&=da+db\leq da+db=\frac{1}{2}\cdot10^{-4}+\frac{1}{2}\cdot10^{-3}=0.55\times10^{-3}<\frac{1}{2}\cdot10^{-2}end{aligned}所以 $a+b$ 有三位有效数字；因为 $a\timesb=0.xxxxxxxx\times10$。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
162
27
x 3 x x 10 x 3 x x 3
.
27
81/ 5
牛顿插值: 首先计算差商
3 1 02 1 3 2 1.333 6 10 4
0.3889 0.8889 0.1420
N3 x 1 x 3 0.3889 x 3 x 1.1420 x 3 x x 3.
也可以利用等距节点构造,首先计算差分
3.141587725... 这个去掉 的无穷乘积公式中第 9 项后的部分产生的误差就是方法误差,也成为截断误差.
2. 按照四舍五入的原则,将下列各数舍成五位有效数字: 816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236
可得前插公式 和后插公式
3 1 02 3 3 2 4 7 6 10 12 16 23
N3
x0
th
1
3t
7 2
t
t
1
23 6
t
t
1
t
2
;
N3
x3
th
10
12t
16 2
t
t
1
23 6
t
t
1
t
2
.
6. 确定一次数不高于 4 的多项式 x ,使 0 0,0 0, 1 1 1, 2 1.
解: 利用重节点计算差商
2
2
2
所以 a b 有三位有效数字.
6.
设 y1
0.9863, y2
0.0062 ,是经过舍入后作为 x1, x2 的近似值.求
1, y1
1 y2
的 计算值 与真
值的相对误差限及 y1 y2 与真值的相对误差限.
解:
已知 x1
y1 dx1, x2
y2
dx2
,
dx1 =
1 2
10-4
,dx2
3
R
dr V dV dlnV d ln R3 3 d ln R 3 d ln R 3dr R 3a 1% ,则 a 1 1% .
V
3
第二章 插值法与数值微分
1. 设 y x ,在 x 100,121,144 三处的值是很容易求得的,试以这三个点建立 y x 的
二次插值多项式,并用此多项式计算 115 的近似值,且给出误差估计.用其中的任意两点,构
00 000 11 11 1 1 1 0 1 2 1 0 1 1/ 2 1/ 4
则可构造 Hermite 插值函数满足题设条件:
H4 x 0 0 x 0 1 x 0 x 0 1 x 0 x 0 x 1 1 x 0 x 0 x 1 x 1
4 1 x4 3 x3 9 x2.
解: 已知 dx ,所以 ln x 的绝对误差 d ln x dx .
x
x
11. 设 x 的相对误差为 % ,求 xn 的相对误差.
解:
dxn xn
nxn1dx xn
nd x x
n % .
12. 计算球的体积,为了使相对误差限为 1%,问度量半径 R 时允许的相对误差限如何?
解: 已知V 4 R3 ,设 dr R dR a ,则要使得
10
x 121
100 100
x 144 121 144
11
x 121 x 100 144 121144 10012
所以 115 L2 115 10.7228 .
误差 R2 x
f
3!
x
x0
x
x1
x
x2
,
x0
,
x1,
x2
,
,所以
0.0006555 R2 115 0.001631
利用前两个节点建立线性插值函数可得:
第一章 误差
1. 试举例,说明什么是模型误差,什么是方法误差.
解: 例如,把地球近似看为一个标准球体,利用公式 A 4 r 2 计算其表面积,这个近似看为球体的过程产生
的误差即为模型误差.
在计算过程中,要用到 ,我们利用无穷乘积公式计算 的值:
2 2 2 ... q1 q2
其中
q1 2, qn1 2 qn , n 2,3,... 我们取前 9 项的乘积作为 的近似值,得
x 7 5 x 7 5 ,7 6
其中,
l0
x
7 0
6
x
x x
7
5 ,7 75
6;
,l17
x 6
.
则可由已知条件,可得
所以可得
hn,i
x
1
2ln,i
xi
x
xi
l2
n,i
x
;
hn,i
x
x
xi
l2
n,i
x
.
n
H2n1 x
1
2ln,i
xi
x
xi
l2
n,i
x
yi
x
xi
l2 n,i
x
yi
i0
8. 过 0,1 两点构造一个三次 Hermite 插值多项式,满足条件:
f 0 1, f 0 1 , f 1 2, f 1 1
lj x
x x0 ... x xj1
x j x0 ... x j x j1
x xj1 ... x xn
x j x j1 ... x j xn
证明
证明: 由于
n
x
k j
l
j
x
xk
,
k
0,1,...,
n
j0
lj
xi
ij
1 0
i j; i j.
n
且 xkj l j x 和 xk 都为 k 次多项式,而且在 k+1 个不同的节点处的函数值都相同 k 0,1,..., n , 所以 j0
L2
x
x 144 100 144
10
x 100 144 100
12
所以 115 L1 115 10.6818 .
与 115 的真实值比较,二次插值比线性插值效果好,利用前两个节点的线性插值比其他两个线性插值
效果好.此说明,二次插值比线性插值效果好,内插比外插效果好.
2. 利用(2.9)式证明 证明: 由(2.9)式
424
7. 寻找过 n 1个点 x0,x1,..., xn 的 2n 1 次多项式 H2n1 x ,满足条件:
H2n1 x0 f x0 , H2n1 x1 f x1 ,..., H2n1 xn f xn , H 2n1 x0 f x0 , H 2n1 x1 f x1 ,.., H 2n1 xn f xn .
Rx
max x0 x x1
f
x
x1
x0 8
2
,
x0
x
x1
Rx
f
2!
x
x0
x
x1
,
x0
x1
当 x0 x x1 时,
f
max x0 x x1
f
x
,
max
x0 xx1
x
x0
x
x1
1 4
x1
x0
2
所以
Rx
max x0 x x1
f
x
x1
x0 8
2
,
x0
x
x1
3. 若 xj 0,1,..., n 为互异节点,且有
L1
x
x 121 100 121
10
x 100 121 100
11
所以 115 L1 115 10.7143 .
利用后两个节点建立线性插值可得:
L1
x
x 144 121 144
11
x 121 144 121
12
所以 115 L1 115 10.7391.
利用前后两个节点建立线性插值可得:
2a 200
8. 用观测恒星的方法求得某地维度为 45 02 (读到秒),试问:计算 sin 将有多大误差?
解:
d sin cos d cos
45 0 2
1 2
.
9 . 真空中自由落体运动距离 s 与时间的关系由公式 s 1 gt2 确定,g 是重力加速度.现在假 2
设 g 是准确的,而对 t 的测量有 0.1s 的误差,证明 t 增加时,距离的绝对误差增加而相对误差
(2) 取第二类边界条件,作三次样条插值多项式.
80 3.153
(3) 用两种插值函数分别计算 x 75.5, x 78.3的函数值.
解: (1)做分段线性插值函数可得:
I5 x 2.768l0 x 2.833l1 x 2.903l2 x 2.979l3 x 3.062l4 x 3.153l5 x
cos
3!
x xi1 x xi x xi1
cos 1;令 g x x xi1 x xi x xi1 ,设 x xi1 th ,得
g xi1 th h3t t 1t 2,t 0, 2
又
t
t
t
1
t
2
,
t
0,
2
的最大值为
1
3 3
0.3849
,所以有
造线性插值函数,用得到的三个线性插值函数,计算 115 的近似值,并分析其结果不同的原
因.
解: 已知 x0 100, x1 121, x2 144; y0 10, y1 11, y2 12 ,
建立二次 Lagrange 插值函数可得:
L2
x
x 100
121 121