数值计算第三章答案
数值计算方法(山东联盟)智慧树知到答案章节测试2023年中国石油大学(华东)
第一章测试1.数值计算方法研究的误差有()A:截断误差;B:观测误差;C: 模型误差；D:舍入误差.答案:AD2.A:只有模型误差、截断误差与观测误差。
B: 只有舍入误差、截断误差与观测误差；C:只有模型误差、观测误差与舍入误差；D:只有模型误差、截断误差与舍入误差；答案:C3.A:4位B:5位C:3位D:2位答案:A4.对于下列表达式，用浮点数运算，精度较高是A:B:C:D:答案:A5.A:B:C:D:答案:B第二章测试1.A:0.5000B:0.6250C:0.5625D:0.6875答案:C2.A:B:C:D:答案:CD3.关于Steffensen(斯蒂芬森)迭代方法，下列命题中正确的是：A:Steffensen迭代法使得收敛的迭代格式加速收敛，发散的迭代格式更快发散。
B:Steffensen迭代法使得某些发散的迭代格式变为收敛。
C:Steffensen迭代法使得任何收敛的迭代格式加速收敛。
D:Steffensen迭代法使得某些收敛的迭代格式加速收敛。
答案:BD4.关于Newton迭代法，下列命题中正确的是：A:求解任一方程的Newton迭代法都是2阶收敛的。
B:Newton迭代格式若收敛，则一定是超线性收敛的。
C:D:Newton迭代格式可能收敛也可能发散。
答案:CD5.A:6B:3C:5D:4答案:A第三章测试1.A:若求解失败，则说明矩阵A奇异。
B:算法的计算量与近似成正比。
C:若A的对角线元素的绝对值都大于1，则求解结果的精度一定较高。
D:只要A非奇异，则求解结果的精度一定较高。
答案:B2.列主元Gauss消去法与Gauss顺序消元法相比，优点是：A:提高了稳定性，减少了误差的影响。
B:方程组的系数矩阵奇异时也可以求解。
C:能求出方程组的精确解。
D:减少了计算量。
答案:A3.A:平方根法与Gauss列主元消去法相比，提高了稳定性，但增加了计算量。
B:只要是对称正定矩阵，就可用平方根法求解。
数值方法简明教程作业集答案
数值计算方法简明教程第一章1 *1x =1.7； *2x =1.73； *3x =1.732 。
2．3． （1） ≤++)(*3*2*1x x x e r 0.00050； （注意：应该用相对误差的定义去求） （2） ≤)(*3*2*1x x x e r 0.50517； （3） ≤)/(*4*2x x e r 0.50002。
4．设6有n 位有效数字，由6≈2.4494……，知6的第一位有效数字1a =2。
令3)1()1(1*1021102211021)(-----⨯≤⨯⨯=⨯=n n r a x ε 可求得满足上述不等式的最小正整数n =4，即至少取四位有效数字，故满足精度要求可取6≈2.449。
5. 答：（1）*x (0>x )的相对误差约是*x 的相对误差的1/2倍；（2）n x )(* 的相对误差约是*x 的相对误差的n 倍。
6． 根据********************sin 21)(cos 21sin 21)(sin 21sin 21)(sin 21)(c b a c e c b a c b a b e c a c b a a e c b S e r ++≤ =******)()()(tgc c e b b e a a e ++ 注意当20*π<<c 时，0**>>c tgc ，即1*1*)()(--<c tgc 。
则有)()()()(****c e b e a e S e r r r r ++<7．设20=y ，41.1*=y ，δ=⨯≤--2*001021y y 由 δ1*001*111010--≤-=-y y y y ，δ2*111*221010--≤-=-y y y yδ10*991*10101010--≤-=-y y y y即当0y 有初始误差δ时，10y 的绝对误差的绝对值将减小1010-倍。
而11010<<-δ，故计算过程稳定。
数值计算方法课后习题答案
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
丁丽娟《数值计算方法》五章课后实验题答案（源程序很详细,且运行无误）
丁丽娟《数值计算方法》五章课后实验题答案（源程序都是自己写的，很详细，且保证运行无误）我做的五章数值实验作业题目如下：第二章：1、2、3、4题第三章：1、2题第四章：1、2题第六章：2、3题第八章：1、2题第二章1：(1) 对A进行列主元素三角分解：function [l u]=myfun(A) n=size(A); for k=1:n for i=k:n sum=0; m=k; for j=1:(k-1) sum=sum+A(i,j)*A(j,k); end s(i)=A(i,k)-sum; if abs(s(m))<abs(s(i)) m=i; end end for j=1:n c=A(m,j); A(m,j)=A(k,j); A(k,j)=c; end for j=k:n sum=0; for r=1:(k-1) sum=sum+A(k,r)*A(r,j); end u(k,j)=A(k,j)-sum; A(k,j)=u(k,j); end for i=1:n l(i,i)=1; end for i=(k+1):n sum=0; for r=1:(k-1) sum=sum+A(i,r)*u(r,k); end l(i,k)=(A(i,k)-sum)/u(k,k); A(i,k)=l(i,k); end end 的列主元素三角分解：求A的列主元素三角分解：>>A=[1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70]; >>[L,U]=myfun(A) 结果：L = 1.0000 0 0 0 0 1.0000 1.0000 0 0 0 1.0000 0.5000 1.0000 0 0 1.0000 0.7500 0.7500 1.0000 0 1.0000 0.2500 0.7500 -1.0000 1.0000 U = 1.0000 1.0000 1.0000 1.0000 1.0000 0 4.0000 14.0000 34.0000 69.0000 0 0 -2.0000 -8.0000 -20.5000 0 0 0 -0.5000 -2.3750 0 0 0 0 -0.2500 (2) 求矩阵的逆矩阵A -1: inv(A) 结果为：ans = 5 -10 10 -5 1 -10 30 -35 19 -4 10 -35 46 -27 6 -5 19 -27 17 -4 1 -4 6 -4 1 （3）检验结果：E=diag([1 1 1 1 1]) A\E ans = 5 -10 10 -5 1 -10 30 -35 19 -4 10 -35 46 -27 6 -5 19 -27 17 -4 1 -4 6 -4 1 2: 程序：程序：function d=myfun(a,b,c,d,n) for i=2:n l(i)=a(i)/b(i-1); a(i)=l(i); u(i)=b(i)-c(i-1)*a(i); b(i)=u(i); y(i)=d(i)-a(i)*d(i-1); d(i)=y(i); end x(n)=d(n)/b(n); d(n)=x(n); for i=(n-1):-1:1 x(i)=(d(i)-c(i)*d(i+1))/b(i); d(i)=x(i); end 求各段电流量程序：求各段电流量程序：for i=2:8 a(i)=-2; end b=[2 5 5 5 5 5 5 5]; c=[-2 -2 -2 -2 -2 -2 -2]; V=220; R=27; d=[V/R 0 0 0 0 0 0 0]; n=8; I=myfun(a,b,c,d,n) 运行程序得：运行程序得：I = 8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.0477 3：程序：（1）求矩阵A和向量b的matlab程序：function [A b]=myfun(n) for i=1:n X(i)=1+0.1*i; end for i=1:n for j=1:n A(i,j)=X(i)^(j-1); end end for i=1:n b(i)=sum(A(i,:)); end 求n=5时A1，b1及A1的2-条件数程序运行结果如下：条件数程序运行结果如下： n=5；[A1,b1]=myfun(n) A1 = 1.0000 1.1000 1.2100 1.3310 1.4641 1.0000 1.2000 1.4400 1.7280 2.0736 1.0000 1.3000 1.6900 2.1970 2.8561 1.0000 1.4000 1.9600 2.7440 3.8416 1.0000 1.5000 2.2500 3.3750 5.0625 b1 = 6.1051 7.4416 9.0431 10.9456 13.1875 cond2=cond(A1,2)cond2 = 5.3615e+005 条件数程序运行结果如下：求n=10时A2，b2及A2的2-条件数程序运行结果如下：n=10; [A2,b2]=myfun(n) A2 = 1.0000 1.1000 1.2100 1.3310 1.4641 1.6105 1.7716 1.9487 2.1436 2.3579 1.0000 1.2000 1.4400 1.7280 2.0736 2.4883 2.9860 3.5832 4.2998 5.1598 1.0000 1.3000 1.6900 2.1970 2.8561 3.7129 4.8268 6.2749 8.1573 10.6045 1.0000 1.4000 1.9600 2.7440 3.8416 5.3782 7.5295 10.5414 14.7579 20.6610 1.0000 1.5000 2.2500 3.3750 5.0625 7.5938 11.3906 17.0859 25.6289 38.4434 1.0000 1.6000 2.5600 4.0960 6.5536 10.4858 16.7772 26.8435 42.9497 68.7195 1.0000 1.7000 2.8900 4.9130 8.3521 14.1986 24.1376 41.0339 69.7576 118.5879 1.0000 1.8000 3.2400 5.8320 10.4976 18.8957 34.0122 61.2220 110.1996 198.3593 1.0000 1.9000 3.6100 6.8590 13.0321 24.7610 47.0459 89.3872 169.8356 322.6877 1.0000 2.0000 4.0000 8.0000 16.0000 32.0000 64.0000 128.0000 256.0000 512.0000 b2 = 1.0e+003 * 0.0159 0.0260 0.0426 0.0698 0.1133 0.1816 0.2866 0.4451 0.6801 1.0230 cond2=cond(A2,2) cond2 = 8.6823e+011 条件数程序运行结果如下：求n=20时A3，b3及A3的2-条件数程序运行结果如下：n=20; [A3,b3]=myfun(n) A3 = 1.0e+009 * Columns 1 through 10 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 Columns 11 through 20 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0005 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0013 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0007 0.0015 0.0032 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0014 0.0032 0.0075 0.0000 0.0000 0.0000 0.0001 0.0002 0.0005 0.0012 0.0029 0.0070 0.0167 0.0000 0.0000 0.0001 0.0001 0.0004 0.0009 0.0023 0.0058 0.0146 0.0364 0.0000 0.0000 0.0001 0.0002 0.0006 0.0017 0.0044 0.0113 0.0295 0.0766 0.0000 0.0001 0.0002 0.0004 0.0011 0.0030 0.0080 0.0215 0.0581 0.1570 0.0000 0.0001 0.0002 0.0007 0.0018 0.0051 0.0143 0.0400 0.1119 0.3133 0.0000 0.0001 0.0004 0.0010 0.0030 0.0086 0.0250 0.0726 0.2105 0.6103 0.0001 0.0002 0.0005 0.0016 0.0048 0.0143 0.0430 0.1291 0.3874 1.1623 b3 = 1.0e+009 * Columns 1 through 10 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0002 0.0004 0.0010Columns 11 through 20 0.0025 0.0059 0.0132 0.0287 0.0606 0.1246 0.2494 0.4874 0.9316 1.7434 cond2=cond(A3,2) cond2 =3.2395e+022 由上述运行结果可知：它们是病态的，而且随着n的增大，矩阵的病态变得严重。
数值分析课后习题及答案
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
数值计算方法第3章3-04范数
是收敛的，称 A 为矩阵序列 A(k) 的收敛极限。
矩阵的收敛
记矩阵序列 A(k) 是收敛于 A 为： lim A(k) A 。 k
Rnn 上 的 矩 阵序 列 A(k) 是 收 敛 于 A 的 充 要 条件 为
lim
k
a(k ij
)
aij
。
其中
a(k ij
矩阵范数的另一个定义 设A Rnn ，矩阵A
A sup Ax
x 1 xR n
的范数
4 常用的矩阵范数
设 A [aij ]nn常用的矩阵范数有行（无穷）范数和列（一）范数。
n
A max aij 1in j1
n
A 1 max aij 1 jn i 1
)
和 aij
分别表示
A( k )
和
A
的第 i 行第
j
列的元素。
定义 设 A Rnn ，如果存在 R 使
Ax x
则称 为A 的一个特征值。x 就是特征值 对应的特征向量。
谱半径
定义 6：对于 Rnn 上的矩阵 A ，设 A 的特
征值为 1, 2 , , n ，称 ( A) max{1, 2 , ,n} 为 矩 阵 A 的 谱 半
但在各种范数下，考虑向量序列收敛性时结论时一致的，一致的含义
是收敛都收敛，且有相同的极限。
提出各种范数是为解不同问题时用的，即对某一个问题可能是某一种
范数方便，而另一种范数不方便。
向量范数的等价定理 给定 x Rn ，对于Rn
，
，总存在与x 无关的正常数m
，M
对一切 x Rn 成立。
数值计算方法(宋岱才版)课后答案
第一章 绪论一 本章的学习要求（1）会求有效数字。
（2）会求函数的误差及误差限。
（3）能根据要求进行误差分析。
二 本章应掌握的重点公式（1）绝对误差：设x 为精确值，x *为x 的一个近似值，称e x x **=-为x *的绝对误差。
（2）相对误差：r e e x***=。
（3）绝对误差限：e x x ε***==-。
（4）相对误差限：r x x xxεε*****-==。
（5）一元函数的绝对误差限：设一元函数()()()0,df f x f x dx εε***⎛⎫==⋅ ⎪⎝⎭则。
（6）一元函数的相对误差限：()()1r df f x dx f εε****⎛⎫=⋅ ⎪⎝⎭。
（7）二元函数的绝对误差限：设一元函数()()(),0,f f x y f y y εε***⎛⎫∂==⋅ ⎪∂⎝⎭则。
（8）二元函数的相对误差限：()()()1r f f f x y x y f εεε******⎡⎤⎛⎫∂∂⎛⎫⎢⎥=⋅+⋅ ⎪ ⎪∂∂⎝⎭⎢⎥⎝⎭⎣⎦。
三 本章习题解析1. 下列各数都是经过四舍五入得到的近似值，（1）试指出它们有几位有效数字，（2）分别估计1123A X X X ***=及224X A X **=的相对误差限。
12341.1021,0.031,385.6,56.430x x x x ****====解：（1）1x *有5位有效数字，2x *有2位有效数字，3x *有4位有效数字，4x *有5位有效数字。
（2）1111123231312123,,,,A A AA x x x x x x x x x x x x ∂∂∂====∂∂∂由题可知：1A *为1A 的近似值，123,,x x x ***分别为123,,x x x 近似值。
所以()()111rA A Aεε***=()()()12311111123A A A x x x A X X X εεε*******⎡⎤⎢⎥=++⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫⎛⎫∂∂∂ ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭43123131212311111010100.215222x x x x x x x x x **-**-**-***⎡⎤=⨯⨯+⨯⨯+⨯⨯=⎢⎥⎣⎦()222222424441,,,X A Ax A X x x x x ∂∂===-∂∂则有同理有2A *为2A 的近似值，2x *，4x *为2x ，4x 的近似值，代入相对误差限公式：()()222rA A Aεε***=()()24212224A A X X A X X εε*****⎡⎤⎢⎥=+⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫∂∂ ⎪ ⎪∂∂⎝⎭⎝⎭()33542224411*********X X X X X **--***⎡⎤⎢⎥=⨯⨯+⨯⨯=⎢⎥⎣⎦2. 正方形的边长大约为100cm ，怎样测量才能使其面积误差不超过21cm ? 解：设正方形的边长为x ，则面积为2S x =，2dsx dx=，在这里设x *为边长的近似值，S *为面积的近似值：由题可知：()()1ds s x dx εε***=≤⎛⎫ ⎪⎝⎭即：()21x x ε**⋅≤ 推出：()10.005200xcm ε*≤=。
现代数值计算课后答案
现代数值计算课后答案【篇一：数值计算课后答案4】>1、设x0?0,x1?1，写出f(x)?e?x的一次插值多项式l1(x)，并估计插值误差。
设插值函数为l1(x)?ax?b，由插值条件，建立线性方程组为?a?0?b?1??1?a?1?b?e?a?e?1?1解之得??b?1则l1(x)?(e?1?1)x?1因为y?(x)??e?x,y??(x)?e?x 所以，插值余项为r(x)?f(x)?p(x)??1f(n?1)(?)?(x)(n?1)!1(2)f(?)?(x)2!1?f(2)(?)(x?x0)(x?x1)2!1?e??(x?0)(x?1)(??(0,1))2所以1r(x)?maxe??maxx(x?1)0?x?120???1。
111??e?0??2482选用合适的三次插值多项式来近似计算f(0.2)和f(0.8)。
解：设三次插值多项式为f(x)?a0?a1x?a2x2?a3x3，由插值条件，建立方程组为?a0?a1?(?0.1)?a2?(?0.1)2?a3?(?0.1)3?0.995?23?a0?a1?0.3?a2?0.3?a3?0.3?0.995?23?a0?a1?0.7?a2?0.7?a3?0.7?0.765?a?a?1.1?a?1.12?a?1.13?0.4 5423?01即?a0?0.1a1?0.01a2?0.001a3?0.995?a0?0.1a1?0.01a2?0.001a3?0 .995?a?0.3a?0.09a?0.027a?0.995?0.4a1?0.08a2?0.028a3?0?0?1 23???a?0.7a?0.49a?0.343a?0.7650.8a1?0.48a2?0.344a3?1.76123?0?? ?0.4a1?0.72a2?0.988a3??0.311?a0?1.1a1?1.21a2?1.331a3?0.45 4??a0?0.1a1?0.01a2?0.001a3?0.995?0.4a1?0.08a2?0.028a3?0???0.32a2?0.288a3?1.76???0.384a3??3.831?解之得 ?a0?0.41?a??6.29?1?a??3.48?2??a3?9.98则所求的三次多项式为f(x)?0.41?6.29x?3.48x2?9.98x3。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
3.1证明：如果求积公式（3.4）对函数f （x ）和g （x ）都准确成立，则它对于线性组合af(x)+bg(x) (a,b 均为常数)亦准确成立. 因此，求积公式（3.4）具有m 次代数精度的充分必要条件是：它对任一小于等于m 次的多项均能准确成立，但对某个m+1次多项式不能准确成立.()()不能成立对与题设矛盾多项式都能准确成立，次多，即对任意的线性组合亦准确成立也能准确成立，则对若对的线性组合亦准确成立对次的多项式准确成立对于任意小于等于不准确成立，对的线性组合亦准确成立对成立次的多项式于等于根据定义可知：对于小次代数精度机械求积公式具有机械求积公式也成立对于线性组合同理可得机械求积公式都成立对于证明：1m 13213213200000)1(,,,,,,1,,,,,1,,,,,1),1,0()(2)()()]()([)()()]()([)()()()()()()()()(),(1++++=======∴+⋅∴⇐∴==∴⇒+∴+=+≈+∴≈≈∴≈≈∴∑∑⎰∑⎰∑⎰∑⎰∑⎰∑x m x x x x x x x x x x m x x x x x m j x x f m m x bg x af x bg x af A x bg A x af A dx x bg x af xbg A dx x bg x af A dx x af x g A dx x g x f A dx x f x g x f m m m m m m j nk k k nk k k ban k k k bank kkbank k k bank k k bank k k 3.2直接验证中矩形公式具有一次代数精度，而Simpson 公式则具有3次代数精度。
右边左边）（时当右边左边）（时当）（公式：次代数精度中矩形公式具有右边左边时当右边左边时当右边左边时当中矩形公式：知：解：根据代数精度定义=-=+++--===-=+++--==+++-≈∴≠--+=+--===-=+--===-=+--==+-≈⎰⎰⎰⎰⎰⎰⎰;2)()2(4)(6)(;2)()(;)()2(4)(6)(;)(1)()()2(4)(6)()(1;4)2()(;3)()(;2)2()(;2)()(;)2()(;)(1)()2()()(222232233322222a b b f b a f a f a b a b dx x f x x f a b b f ba f a f ab a b dx x f x f b f ba f a f ab dx x f Simpson a b a ab b b a f a b a b dx x f x x f a b b a f a b a b dx x f x x f a b ba f ab a b dx x f x f ba f ab dx x f b a b a ba b a b a b a ba次代数精度公式具有右边，左边）（时当右边左边）（时当右边左边）（时当3;242255)()2(4)(6)(;5)()(;4)()2(4)(6)(;4)()(;3)()2(4)(6)(;3)()(234324555544444333332Simpon b a ab b a b a a b b f b a f a f a b a b dx x f x x f a b b f b a f a f a b a b dx x f x x f a b b f b a f a f a b a b dx x f x x f b a b a ba ∴≠--++-=+++--===-=+++--===-=+++--==⎰⎰⎰3.3 已知数据表试分别用Simpson 法与复合梯形法计算积分dx e x ⎰5.11.1.）（复合梯形法：由题意知法：解：4817.46693.320042.31.0)]1.1()21.15.1(2)1.1([221.15.1247754.1)4817.46693.340042.3(151)]1.1()21.15.1(4)1.1([6)1.15.1(5.11.15.11.1+⨯+⨯=+++-≈==+⨯+=+++-≈⎰⎰f f f dx e n f f f dx e Simpon x x3.4若,0)(''>x f 证明用梯形求积公式计算积分dx x f ba⎰)(所得结果比准确值大，并说明几何意义.线围成的曲面面积为该直线与被积函数曲的直线，，梯形插值函数为连接被积函数为严格凸函数几何意义：果比准确值大梯形求积公式得到的结即）（为梯形求积公式，则的准确值，为证明：设1''11''''311))(,()),(,(0)(00)(,0],[)(12--)(T I b f b a f a x f T I T I x f a b b a f a b T I T dx x f I ba -∴>∴<<-∴>>-∈-=⎰ ξξ3.5分别用复合梯形法和复合Simpson 法计算积分dx e x ⎰1，怎样取n 才能保证计算结算结果有6位有效数字.178284.1))(2)21(4)1()0((46117828.1))(2)1()0((170214000005.0|1)-(e 1611801||(a))f -(b)(f )2(1801-||S -I |170000005.0|1)-(e 121||(a))f -(b)(f 12-||T -I |30314169117043'34n 2''22n 10=++++⨯==++⨯=≥≤=-≈≥≤-=-≈∑∑∑⎰===k k k k k k xx f x f f f S x f f f T n nn a b n nn a b dx e I 解得解得）（的准确值，为解：设3.6设⎪⎩⎪⎨⎧≤≤-+-+-+≤≤-+-+-+≤≤+=.3.02.0,)2.0(2)2.0(9.0)2.0(15.0009.1;2.01.0,)1.0(2)1.0(3.0)1.0(3.0001.11.00,1)(32323x x x x x x x x x x x f 分别用复合梯形法（n=6）和复合Simpson 法（n=3）计算积分dx x f ⎰3.00)(，并估计误差.)]()(2)([2516b f x f a f hT k k ++=∑=解：0]612[)63.0(1801)]()([)2(180100008125.0]039.0[36123.0)]()([12302425.0]035.1)009.1001.1(2)019.10035.1000125.1(41[183.0]3.0())2.0()1.0((2))25.0()15.0()05.0((4)0([183.0)]()21(2)21(4)([630250625.0]035.1)019.1009.10035.1001.1000125.1(21[123.0)]3.0()25.0()2.0()15.0()1.0()05.0((2)0([123.0433462''2621203≈--=--≈-=-⨯-=--≈-=++++++=++++++=+++++==++++++⨯=++++++⨯=∑∑==a f b f h S I a f b f h T I f f f f f f f b f x f x f a f h S f f f f f f f k k k k3.7导出中矩形公式的余项.3''2'''2'''))((241)2)((21)2()2(])2)((21)2()2([)]2()([)2()()2()()()2()()(a b f dxb a x f dx b a x b a f dx b a x f b a x b a f dxba f x f dxb a f dx x f ba f ab dx x f R ba f ab dx x f b a b a b a b a b a ba ba-=+-++-+=+-++-+=+-=+-=+--=+-≈⎰⎰⎰⎰⎰⎰⎰ηηη解：中矩形公式： 3.8.（略）3.9设（3.6）是Gauss 公式，证明它的求积系数恒大于零，求积系数之和等于2.12auss );()(1111k 11-∴--∏=≈∑⎰⎰=-≠=n G dx x x x x A x f A dx x f nk jk j nkj j k k 上式的代数精度为公式是其中证明：2x d 11)(111===∴∑⎰=-nk k A x x f 即准确成立上式对于3.10 1)证明)53(95)0(98)53(95)(11f f f dx x f ++-≈⎰-是Gauss 求积公式. 2)用3点Gauss 求积公式计算积分⎰-12dx ex .公式为且节点的代数精度为右边左边时当右边左边时当右边左边时当右边左边时当右边左边时当右边左边时当右边左边时）：当证明：Gauss f f f dx x f n f f f dx x f f f f dx x f x x f f f f dx x f x x f f f f dx x f x x f f f f dx x f x x f f f f dx x f x x f f f f dx x f x x f f f f dx x f x f )53(95)0(98)53(95)(132535)53(95)0(98)53(95)(;256)53(95)0(98)53(95;72)()(;0)53(95)0(98)53(95;0)()(;52)53(95)0(98)53(95;52)()(;0)53(95)0(98)53(95;0)()(;32)53(95)0(98)53(95;32)()(;0)53(95)0(98)53(95;0)()(;2)53(95)0(98)53(95;2)(1)(111111161151141131121111++-≈∴-⨯==++-≈∴≠=++-====++-====++-====++-====++-====++-====++-==⎰⎰⎰⎰⎰⎰⎰⎰⎰--------- 746814.0]959895[212122222]215321[41]21)53(21[11)2121(10≈++≈=+⨯--+-⨯--+--⎰⎰e e e dt edx e t x ）：3.11考虑求积公式),1()()(10Bf x Af dx x f +≈⎰选取求积系数A ，B 和求积节点0x ，使得求积公式具有尽可能高的代数精度，并指出达到的最高代数精度的次数.2;3610)1(41)31(43;41)()1(41)31(43)(;31;41;43321)3(21)()()2(21)()()1(1)(1)(10331020102011代数精度为右边左边时当）组成的方程组的）（）（解（时当时当时解：当∴≠=+==+≈∴===+===+===+===⎰⎰⎰⎰⎰f f dx x x x f f f dx x f C B A B Ax dx x f x x f B Ax dx x f x x f B A dx x f x f3.12设已给出2)1(1)(x x f +=的函数表试用三点公式计算f ’(x)在x=1.0,1.1,1.2的值，并估计误差.0025.001.0375.0|3)(||)2.1()2.1(|00125.001.0675.0|6)(||)1.1()1.1(|0025.001.0375.0|3)(||)0.1()0.1(|],[75.0)()1(24)(187.0)2066.032268.042500.0(2.01))2.1(3)1.1(4)0.1((21)2.1(217.0)2066.02500.0(2.01))2.1()0.1((21)1.1(247.0)2066.02268.042500.03(2.01))2.1()1.1(4)0.1(3(21)0.1(1.00.11.123'2'23'2'23'2'20353'''01=⨯≤=-=⨯≤-=-=⨯≤=-∈≤+-=-=⨯+⨯-=+-≈-=+-=+-≈-=-⨯+⨯-=-+-≈=-=-=-h f P f h f P f h f P f x x f x x f f f f h f f f h f f f f h f x x h ξξξξξ（注：文档可能无法思考全面，请浏览后下载，供参考。