数值分析(第五版)计算实习题第四章作业
数值分析第四章习题
数值分析第四章习题第四章习题1. 采用数值计算方法，画出dt t t x y x ?=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
〖答案〗1.65412. 求函数x e x f 3sin )(=的数值积分?=π 0 )(dx x f s ，并请采用符号计算尝试复算。
〖答案〗s = 5.1354Warning: Explicit integral could not be found. > In sym.int at 58s =int(exp(sin(x)^3),x = 0 .. pi)3. 用quad 求取dx x e x sin 7.15?--ππ的数值积分，并保证积分的绝对精度为910-。
〖答案〗1.087849437547794. 求函数5.08.12cos 5.1)5(sin )(206.02++-=t t t et t f t 在区间]5,5[-中的最小值点。
〖答案〗最小值点是-1.28498111480531 相应目标值是-0.186048010065455. 设0)0(,1)0(,1)(2)(3)(22===+-dt dy y t y dt t dy dt t y d ，用数值法和符号法求5.0)(=t t y 。
〖答案〗数值解y_05 = 0.78958020790127符号解ys =1/2-1/2*exp(2*t)+exp(t)ys_05 =.789580356470605529168507052137806. 求矩阵b Ax =的解，A 为3阶魔方阵，b 是)13(?的全1列向量。
〖答案〗x =0.06670.06670.06677. 求矩阵b Ax =的解，A 为4阶魔方阵，b 是)14(?的全1列向量。
〖答案〗解不唯一x =-0.0074 -0.0809 0.1397 0.0662 0.0588 0.1176 -0.0588。
《数值分析》第四章答案
习题41. 给定x x f =)(在144,121,100=x 3点处的值，试以这3点建立)(x f 的2次(抛物)插值公式，利用插值公式115求的近似值并估计误差。
再给13169=建立3次插值公式，给出相应的结果。
解：x x f =)( 2121)(-='x x f ，2341)(--=''x x f ,2583)(-='''x x f ,27)4(1615)(--=x x f，72380529.10)115(=f1000=x ， 1211=x ， 1442=x ， 1693=x 100=y ， 111=y ， 122=y ， 133=y))(())(())(())(())(())(()(1202102210120*********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 ----+----+----= )121144)(100144()121115)(100115(12)144121)(100121()144115)(100115(11)144100)(121100()144115)(121115(10)115(2----⨯+----⨯+----⨯=L=2344)6(1512)23(21)29(1511)44)(21()29)(6(10⨯-⨯⨯+-⨯-⨯⨯+----⨯72276.1006719.190683.988312.1=-+=))()((!3)()()(2102x x x x x x f x L x f ---'''=-ξ ，144100<<ξ )44115()121115()100115()(max 61)115()115(1441002-⨯-⨯-⋅'''≤-≤≤x f L f x 296151083615⨯⨯⨯⨯⨯≤-001631.0101631.02=⨯=- 实际误差 22101045.0)115()115(-⨯=-L f))()(())()(())()(())()(()(312101320130201032103x x x x x x x x x x x x y x x x x x x x x x x x x y x L ------+------= ))()(())()(())()(())()((23130321033212023102x x x x x x x x x x x x y x x x x x x x x x x x x y ------+------+ )169100()144100()121100()169115()144115()121115(10)115(3-⨯-⨯--⨯-⨯-⨯=L )169121()144121()100121()169115()144115()100115(11-⨯-⨯--⨯-⨯-⨯+)169144()121144()100144()169115()121115()100115(12-⨯-⨯--⨯-⨯-⨯+)144169()121169()100169()144115()121115()100115(13-⨯-⨯--⨯-⨯-⨯+)48()23(21)54()29(1511)69()44()21()54()29()6(10-⨯-⨯-⨯-⨯⨯+-⨯-⨯--⨯-⨯-⨯= 254869)29()6(1513)25(2344)54()6(1512⨯⨯-⨯-⨯⨯+-⨯⨯-⨯-⨯⨯+ 723571.10409783.0305138.2145186.11473744.1=+-+= ))()()((!4)()()(3210)4(3x x x x x x x x f x L x f ----=-ξ，169100<<ξ)169115)(144115)(121115)(10115(101615241)115()115(73----⨯⨯⨯≤--L f )54()29()6(151016152417-⨯-⨯-⨯⨯⨯⨯=- 0005505.0105505.03=⨯=-实际误差 321023429.0)115()115(-⨯=-L f 2. 设j x 为互异节点),,1,0(n j =求证： (1)k nj j k j x x l x =∑=)(0),,1,0(n k =;(2)0)()(0=-∑=x l x x j knj j ),,1(n k =。
数值分析课程第五版课后习题答案(李庆扬等)1之欧阳育创编
第一章 绪论（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 n x n 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 n k 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 = 1.7863 × 10 − 2 。 55.982
8、当 N 充分大时，怎样求 ∫ [解]因为 ∫
N +1 N
1 dx ？ 1+ x2
1 dx = arctan( N + 1) − arctan N ，当 N 充分大时为两个相近数相 1+ x2
减，设 α = arctan( N + 1) ， β = arctan N ，则 N + 1 = tan α ， N = tan β ，从而 tan(α − β ) = 因此 ∫
5、计算球体积要使相对误差限为 1%，问度量半径 R 允许的相对误差是多少？ 4 ε * ( π (R* )3 ) 4 3 [解]由 1% = ε r* ( π ( R * ) 3 ) = 可知， 4 3 * 3 π (R ) 3 ′ 4 4 4 ε * ( π ( R * ) 3 ) = 1% × π ( R * ) 3 = π ( R * ) 3 ε * ( R * ) = 4π ( R * ) 2 × ε * ( R * ) ， 3 3 3
ε * ( y n ) = 10ε * ( y n −1 ) = 10 n ε * ( y 0 ) ，
1 1 从而 ε * ( y10 ) = 1010 ε * ( y 0 ) = 1010 × × 10 − 2 = × 10 8 ，因此计算过程不稳定。 2 2 12、计算 f = ( 2 − 1) 6 ，取 2 ≈ 1.4 ，利用下列公式计算，哪一个得到的结果最 好？ 1 ( 2 + 1)
* r
x= x
*
ε ( x * ) = n( x * ) n −1 2% x * = 2n% ⋅ x * ，
数值分析第五版计算实习题
弟二草插值法3.卜列数据点的插值可以得到平方根函数的近似，在区间064］上作图。
(1〉用这9个点做8次多项式插值Q x)。
(2)用三次样条(第一边界条件)程岸求S(X)。
从得到结果石在［0.64］ 1：・哪个插值更粘确：在区间［0,1］ I：•两种插值哪个更精确?(1) 8次多项式插值：(1)8次多项式插值:首先建立新的M-file:输入如卜代码(此为拉格朗口插值的功能函数)并保存function f=Language(x,y,x0)%求Li知数据点的拉格朗Fl插值多项式%己知数据点的x坐标向量：x%已知数据点的y坐标向量:y%插值的x坐标：x0%求得的拉格朗H插值多项式或在X0处的插值：fsyms t;ifi(lcngth(x)=length(y))n=length(x);elsedisp(*x和y的维数不相等!)；return;end %检错tbr(i=l:n)i=y(i);fbr(j=1:i-l)l=l*(t-x(j))/(x(i)-x(j))；end;for(j=i-M:n)end;for(j=i+l:n) l=l*(t-x(j))/(x(i)-x(j))； end;simplify(f);if(i==n) if|nargin=3)f=subs(C't\xO);else f=collcct(f);f=vpa(f,6);endendend再建立新的M-file：输入：clear;x=[0 1 49 16 25 36 49 64];y=[0：l：8]；%计算拉格朗口基丞数%计算拉格朗ri插值函数%化简%计算插值点的曲数值%将插值多项式展开%将插值多项式的系数化成6位精度的小数f=Uinguage(x,y) 运行得到f=1.32574*1-381410*t A2+.604294e-1 *t A3+.222972e-3 *t A5-.542921 e-5*t A6+.671268e・7T7・.328063e・9T8・.498071 e-2*t A4 这就是8次多项式插值L s(x)= 1.32574怜.381410*t A2+.604294e-1 *t A3+.222972e-3 *t A5-.542921 e-5*t A6+.671268e-7*t A7-.328063e-9*t A8-.498071 e-2*t A4. (2)三次样条插值：建立新的M-filc：输入：clear;x=[0 I 49 1625 36 4964];尸[0:8];t=[0:0.1:64];Y=t.A(0.5);O=Language(x,y)f= 1,32574*t-.381410*t.A2+.604294e-1 *t.A3+.222972e-3*t.A5-.542921 e・5*(. W+.671268e-7*t.A7-.328063e-9*t.A8-.498071 e-2 *t.A4;S=interp l(x,y,t.'spline,);plol(x,y,o；(・YY.lf.'b'」S'g:');grid;运行程序得到如下图：从结果屮很明显可以看出在[0.64].上.三次样条插值更精确，儿乎与原函数帀合。
最新数值分析第四章数值积分与数值微分习题答案
第四章 数值积分与数值微分1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度：10121012112120(1)()()(0)();(2)()()(0)();(3)()[(1)2()3()]/3;(4)()[(0)()]/2[(0)()];hhhh hf x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-⎰⎰⎰⎰解：求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。
（1）若101(1)()()(0)()hhf x dx A f h A f A f h --≈-++⎰令()1f x =，则1012h A A A -=++令()f x x =，则110A h A h -=-+令2()f x x =，则3221123h h A h A -=+ 从而解得011431313A h A h A h -⎧=⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =，则3()0hhhhf x dx x dx --==⎰⎰101()(0)()0A f h A f A f h --++=令4()f x x =，则4551012()52()(0)()3hhhhf x dx x dx h A f h A f A f h h ---==-++=⎰⎰故此时，101()()(0)()hhf x dx A f h A f A f h --≠-++⎰故101()()(0)()h hf x dx A f h A f A f h --≈-++⎰具有3次代数精度。
（2）若21012()()(0)()hhf x dx A f h A f A f h --≈-++⎰令()1f x =，则1014h A A A -=++令()f x x =，则110A h A h -=-+令2()f x x =，则32211163h h A h A -=+ 从而解得11438383A h A h A h -⎧=-⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =，则22322()0hhhhf x dx x dx --==⎰⎰101()(0)()0A f h A f A f h --++=令4()f x x =，则22452264()5hhhhf x dx x dx h --==⎰⎰510116()(0)()3A f h A f A f h h --++=故此时，21012()()(0)()hhf x dx A f h A f A f h --≠-++⎰因此，21012()()(0)()h hf x dx A f h A f A f h --≈-++⎰具有3次代数精度。
数值分析第五版课后习题答案
数值分析第五版课后习题答案数值分析是一门应用数学的分支学科，主要研究如何利用数值方法解决实际问题。
在学习这门课程的过程中，课后习题是不可或缺的一部分。
本文将对《数值分析第五版》的课后习题进行一些探讨和解答。
第一章是数值分析的导论，主要介绍了误差分析和计算方法的基本概念。
在课后习题中，有一道题目是关于误差传播的。
假设有一个函数f(x, y) = x^2 + y^2，其中x和y的测量误差分别为Δx和Δy，要求计算f(x, y)的误差。
解答：根据误差传播公式，可以得到f(x, y)的误差为Δf = √[(∂f/∂x)^2 *(Δx)^2 + (∂f/∂y)^2 * (Δy)^2]。
对于本题而言，∂f/∂x = 2x，∂f/∂y = 2y。
代入公式，得到Δf = √[(2x)^2 * (Δx)^2 + (2y)^2 * (Δy)^2] = 2√(x^2 * (Δx)^2+ y^2 * (Δy)^2)。
第二章是插值与多项式逼近的内容。
其中一道习题涉及到拉格朗日插值多项式。
给定n+1个数据点(x0, y0), (x1, y1), ..., (xn, yn)，要求构造一个n次多项式p(x)，使得p(xi) = yi (i = 0, 1, ..., n)。
解答：拉格朗日插值多项式的表达式为p(x) = Σ(yi * Li(x))，其中Li(x) = Π[(x - xj) / (xi - xj)]，j ≠ i。
将数据点代入表达式中，即可得到所求的多项式。
第三章是数值微积分的内容，其中一道习题是关于数值积分的。
给定一个函数f(x)，要求使用复化梯形公式计算定积分∫[a, b]f(x)dx。
解答：复化梯形公式的表达式为∫[a, b]f(x)dx ≈ h/2 * [f(a) + 2Σf(xi) + f(b)]，其中h = (b - a)/n，xi = a + i * h (i = 1, 2, ..., n-1)。
根据给定的函数f(x)，代入公式中的各个值，即可得到近似的定积分值。
数值分析计算实习题答案
数值分析计算实习题答案数值分析计算实习题答案数值分析是一门研究如何利用计算机对数学问题进行近似求解的学科。
在数值分析的学习过程中，实习题是一种重要的学习方式，通过实践来巩固理论知识，并培养解决实际问题的能力。
本文将为大家提供一些数值分析计算实习题的答案，希望能够帮助大家更好地理解和掌握数值分析的相关知识。
一、插值与拟合1. 已知一组数据点，要求通过这些数据点构造一个一次插值多项式，并求出在某一特定点的函数值。
答案：首先，我们可以根据给定的数据点构造一个一次插值多项式。
假设给定的数据点为(x0, y0), (x1, y1)，我们可以构造一个一次多项式p(x) = a0 + a1x，其中a0和a1为待定系数。
根据插值条件，我们有p(x0) = y0，p(x1) = y1。
将这两个条件代入多项式中，可以得到一个方程组，通过求解这个方程组，我们就可以确定a0和a1的值。
最后，将求得的多项式代入到某一特定点，就可以得到该点的函数值。
2. 已知一组数据点，要求通过这些数据点进行最小二乘拟合，并求出拟合曲线的表达式。
答案：最小二乘拟合是一种通过最小化误差平方和来找到最佳拟合曲线的方法。
假设给定的数据点为(x0, y0), (x1, y1)，我们可以构造一个拟合曲线的表达式y =a0 + a1x + a2x^2 + ... + anx^n，其中a0, a1, ..., an为待定系数。
根据最小二乘拟合原理，我们需要最小化误差平方和E = Σ(yi - f(xi))^2，其中yi为实际数据点的y值，f(xi)为拟合曲线在xi处的函数值。
通过求解这个最小化问题，我们就可以确定拟合曲线的表达式。
二、数值积分1. 已知一个函数的表达式，要求通过数值积分的方法计算函数在某一区间上的定积分值。
答案：数值积分是一种通过将定积分转化为数值求和来近似计算的方法。
假设给定的函数表达式为f(x)，我们可以将定积分∫[a, b]f(x)dx近似为Σwi * f(xi)，其中wi为权重系数，xi为待定节点。
数值分析第五版计算实习题
数值分析计算实习题第二章2-1程序：clear;clc;x1=[0.2 0.4 0.6 0.8 1.0];y1=[0.98 0.92 0.81 0.64 0.38];n=length(y1);c=y1(:);for j=2:n %求差商for i=n:-1:jc(i)=(c(i)-c(i-1))/(x1(i)-x1(i-j+1));endendsyms x df d;df(1)=1;d(1)=y1(1);for i=2:n %求牛顿差值多项式df(i)=df(i-1)*(x-x1(i-1));d(i)=c(i)*df(i);enddisp('4次牛顿插值多项式');P4=vpa(collect((sum(d))),5) %P4即为4次牛顿插值多项式,并保留小数点后5位数pp=csape(x1,y1, 'variational');%调用三次样条函数q=pp.coefs;disp('三次样条函数');for i=1:4S=q(i,:)*[(x-x1(i))^3;(x-x1(i))^2;(x-x1(i));1];S=vpa(collect(S),5)endx2=0.2:0.08:1.08;dot=[1 2 11 12];figureezplot(P4,[0.2,1.08]);hold ony2=fnval(pp,x2);x=x2(dot);y3=eval(P4);y4=fnval(pp,x2(dot));plot(x2,y2,'r',x2(dot),y3,'b*',x2(dot),y4,'co');title('4次牛顿插值及三次样条');结果如下：4次牛顿插值多项式P4 = - 0.52083*x^4 + 0.83333*x^3 - 1.1042*x^2 + 0.19167*x + 0.98三次样条函数x∈[0.2,0.4]时， S = - 1.3393*x^3 + 0.80357*x^2 - 0.40714*x + 1.04 x∈[0.4,0.6]时，S = 0.44643*x^3 - 1.3393*x^2 + 0.45*x + 0.92571x∈[0.6,0.8]时，S = - 1.6964*x^3 + 2.5179*x^2 - 1.8643*x + 1.3886 x∈[0.8,1.0]时，S =2.5893*x^3 - 7.7679*x^2 + 6.3643*x - 0.80571输出图如下2-3（1）程序：clear;clc;x1=[0 1 4 9 16 25 36 49 64];y1=[0 1 2 3 4 5 6 7 8];%插值点n=length(y1);a=ones(n,2);a(:,2)=-x1';c=1;for i=1:nc=conv(c,a(i,:));endq=zeros(n,n);r=zeros(n,n+1);for i=1:n[q(i,:),r(i,:)]=deconv(c,a(i,:));%wn+1/(x-xk) end 三次样条插值曲线4次牛顿插值曲线Dw=zeros(1,n);for i=1:nDw(i)=y1(i)/polyval(q(i,:),x1(i));%系数endp=Dw*q;syms x L8;for i=1:nL8(i)=p(n-i+1)*x^(i-1);enddisp('8次拉格朗日插值');L8=vpa(collect((sum(L8))),5)xi=0:64;yi=polyval(p,xi);figureplot(xi,yi,x1,y1,'r*');hold ontitle('8次拉格朗日插值');结果如下：8次拉格朗日插值L8 =- 3.2806e-10*x^8 + 6.7127e-8*x^7 - 5.4292e-6*x^6 + 0.00022297*x^5 - 0.0049807*x^4 + 0.060429*x^3 - 0.38141*x^2 + 1.3257*x输出图如下：第五章4-1（3）程序：clc;clear;y= (x)sqrt(x).*log(x);a=0;b=1;tol=1e-4;p=quad(y,a,b,tol);fprintf('采用自适应辛普森积分结果为： %d \n', p); 结果如下：采用自适应辛普森积分结果为： -4.439756e-01第九章9-1(a)程序：clc;clear;a=1;b=2;%定义域h=0.05;%步长n=(b-a)/h;y0=1;%初值f= (x,y) 1/x^2-y/x;%微分函数Xn=linspace(a,b,n+1);%将定义域分为n等份Yn=zeros(1,n);%结果矩阵Yn(1)=y0;%赋初值%以下根据改进欧拉公式求解for i=1:nxn=Xn(i);xnn=Xn(i+1);yn=Yn(i);yp=yn+h*f(xn,yn);yc=yn+h*f(xnn,yp);yn=(yp+yc)/2;Yn(i+1)=yn;endXn=Yn;%以下根据经典四阶R-K法公式求解for i=1:nxn=Xn(i);yn=Yn(i);k1=f(xn,yn);k2=f(xn+h/2,yn+h/2*k1);k3=f(xn+h/2,yn+h/2*k2);k4=f(xn+h,yn+h*k3);yn=yn+h/6*(k1+2*k2+2*k3+k4);Yn(i+1)=yn;enddisp(' 改进欧拉法四阶经典R-K法'); disp([Xn' Yn'])结果如下：改进欧拉法四阶经典R-K法1 10.99887 0.998850.99577 0.99780.99114 0.996940.98532 0.996340.97857 0.996030.97111 0.996060.96311 0.996450.9547 0.997230.94598 0.998410.93705 10.92798 1.0020.91883 1.00440.90964 1.00730.90045 1.01060.89129 1.01430.88218 1.01840.87315 1.02290.86421 1.02780.85538 1.03310.84665 1.0388（b）程序：clc;clear;a=0;b=1;%定义域H=[0.1 0.025 0.01];%步长y0=1/3;%初值f= (x,y) -50*y+50*x^2+2*x;%微分函数xi=linspace(a,b,11);Y=1/3*exp(-50*xi)+xi.^2;%准确解Ym=zeros(1,11);for j=1:3h=H(j);n=(b-a)/h;Xn=linspace(a,b,n+1);%将定义域分为n等份Yn=zeros(1,n);%结果矩阵Yn(1)=y0;%赋初值for i=1:nxn=Xn(i);yn=Yn(i);k1=f(xn,yn);k2=f(xn+h/2,yn+h/2*k1);k3=f(xn+h/2,yn+h/2*k2);k4=f(xn+h,yn+h*k3);yn=yn+h/6*(k1+2*k2+2*k3+k4);Yn(i+1)=yn;endfor k=1:11m=0.1/h;Ym(k)=Yn(1+(k-1)*m);enddelta=Ym-Y;fprintf('步长为： %d \n', h);disp(' 四阶经典R-K法准确解误差'); disp([Ym' Y' delta'])end结果如下：步长为： 1.000000e-01四阶经典R-K法准确解误差0.33333 0.33333 04.6055 0.012246 4.593263.062 0.040015 63.022864.05 0.09 863.9611844 0.16 118431.6235e+05 0.25 1.6235e+052.2256e+06 0.36 2.2256e+063.0509e+07 0.49 3.0509e+074.1823e+08 0.64 4.1823e+085.7333e+09 0.81 5.7333e+097.8594e+10 1 7.8594e+10步长为： 2.500000e-02四阶经典R-K法准确解误差0.33333 0.33333 00.013015 0.012246 0.000768940.040063 0.040015 4.82e-050.090037 0.09 3.6857e-050.16004 0.16 3.6723e-050.25004 0.25 3.6722e-050.36004 0.36 3.6722e-050.49004 0.49 3.6722e-050.64004 0.64 3.6722e-050.81004 0.81 3.6722e-051 1 3.6722e-05步长为： 1.000000e-02四阶经典R-K法准确解误差0.33333 0.33333 00.012256 0.012246 9.5673e-060.040016 0.040015 7.8252e-070.090001 0.09 6.6347e-070.16 0.16 6.6226e-070.25 0.25 6.6225e-070.36 0.36 6.6225e-070.49 0.49 6.6225e-070.64 0.64 6.6225e-070.81 0.81 6.6225e-071 1 6.6225e-07由结果可知，步长越小，结果越精确。
李庆扬-数值分析第五版第4章习题答案(20130714)
4h
A1
A0
A1
0 hA1 hA1
16 3
h3
h2 A1
h2 A1
4h A1 A0 A1
解出
8 A1 3 h
A0
4 3
h
A1
8 3
h
令 f (x) x3 ，
有 8h (h)3 4h 03 8h h3 8 h4 8 h4 0 h x3dx 0
有 h (h)4 4h 04 h h4 2 h5 0 h x4dx 1 x5 h 2 h5
3
3
33
h
5 h 5
因此，具有 3 次代数精度。
2h
2） 2h f (x)dx A1 f (h) A0 f (0) A1 f (h)
由于需要确定 3 个 未知量，因此，需要给定 3 个方程。
3
3
3
33
h
令 f (x) x4 ，
有 8h (h)4 4h 04 8h h4 16 h5 0 h x4dx 1 x5 2h 64 h5
3
3
3
3
h
5 2h 5
因此，具有 3 次代数精度。
1
3） 1 f (x)dx [ f (1) 2 f (x1) 3 f (x2)] / 3
T (k) m
4m 4m 1
T (k 1) m1
1 4m 1
T (k) m1
,
k
1, 2,3...
9、什么是高斯型求积公式？它的求积节点是如何确定的？它的代数精确度是多少？为何称 它是具有最高代数精确度的求积公式？ 答： 如果求积公式
an
n
f (x)(x) d x Ak f (xk )
数值分析第五版全答案chap4
第四章 数值积分与数值微分1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度： 1012101211212(1)()()(0)();(2)()()(0)();(3)()[(1)2()3()]/3;(4)()[(0)()]/2[(0)()];h h h h hf x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-⎰⎰⎰⎰解：求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。
（1）若101(1)()()(0)()h hf x dx A f h A f A f h --≈-++⎰令()1f x =，则1012h A A A -=++令()f x x =，则110A h A h -=-+令2()f x x =，则3221123h h A h A -=+从而解得 011431313A h A h A h -⎧=⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =，则3()0h h hhf x dx x dx --==⎰⎰101()(0)()0A f h A f A f h --++=故101()()(0)()h hf x dx A f h A f A f h --=-++⎰成立。
令4()f x x =，则4551012()52()(0)()3h h hhf x dx x dx hA f h A f A f h h---==-++=⎰⎰故此时，101()()(0)()h hf x dx A f h A f A f h --≠-++⎰故101()()(0)()h hf x dx A f h A f A f h --≈-++⎰具有3次代数精度。
数值分析第五版全答案chap4
第四章数值积分与数值微分1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度：h(1)f(x)dx A f(h)A f(0)A f(h);101h2h(2)f(x)dx A f(h)A f(0)A f(h);10 1 2h1(3)f(x)dx[f(1)2f(x)3f(x)]/3;121h2(4)f(x)dx h[f(0)f(h)]/2ah[f(0)f(h)];解：求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。
（1）若h(1)f(x)dx A f(h)A f(0)A f(h)101 h令f(x)1，则2h A A A101令f(x)x，则0A h A h11令2f(x)x，则2 3322h h A h A11从而解得4A h31A h131A h13令3f(x)x，则h h3f(x)dx x dx0 hhA1f(h)A0f(0)A1f(h)0 h故h f(x)dx A f(h)A f(0)A f(h)成立。
101令4f(x)x，则h h452f(x)dx x dx h hh52A f(h)A f(0)A f(h)h10135故此时，hh f(x)dx A f(h)A f(0)A f(h)101h故h f(x)dx A f(h)A f(0)A f(h)101具有3次代数精度。
2h （2）若2h f(x)dx A f(h)A f(0)A f(h)101令f(x)1，则4h A A A101令f(x)x，则0A h A h11令2f(x)x，则16 3322h h A h A11从而解得4A h38A h138A h13令3f(x)x，则2h2h3f(x)dx x dx0 2h2hA1f(h)A0f(0)A1f(h)0 2h故2h f(x)dx A f(h)A f(0)A f(h)成立。
数值分析课程第五版课后习题答案(李庆扬等)1
第一章 绪论（12）之阳早格格创做1、设0>x ，x 的相对付缺面为δ，供x ln 的缺面.[解]设0*>x 为x 的近似值，则有相对付缺面为δε=)(*x r ，千万于缺面为**)(x x δε=，进而xln 的缺面为δδεε=='=*****1)()(ln )(ln x xx x x ，相对付缺面为****ln ln )(ln )(ln xxx x rδεε==.2、设x 的相对付缺面为2%，供n x 的相对付缺面.[解]设*x 为x 的近似值，则有相对付缺面为%2)(*=x r ε，千万于缺面为**%2)(x x =ε，进而nx 的缺面为nn x x n x n 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 . [解]53232323*42*4*2*2*41***4*2*1088654.01021)430.56(461.561021)430.56(461.561021)430.56(031.01021430.561)()()(1)()/(-----=⨯≈⨯⨯=⨯⨯=⨯⨯+⨯⨯=+=⎪⎪⎭⎫⎝⎛∂∂=∑x x x x x x x f x x e n k k kεεε. 5、预计球体积要使相对付缺面限为1%，问度量半径R 允许的相对付缺面是几？ [解]由3*3**3**)(34))(34())(34(%1R R R r ππεπε==可知，)()(4)()(34)(34%1))(34(**2***3*3*3**R R R R R R επεπππε⨯='⎥⎦⎤⎢⎣⎡=⨯=， 进而***31%1)(R R ⨯=ε，故300131%1)()(*****=⨯==RR R r εε.6、设280=Y ，按递推公式),2,1(78310011 =-=-n Y Y n n 预计到100Y ，若与982.27783≈（五位灵验数字，）试问预计100Y 将有多大缺面？[解]令n Y 表示n Y 的近似值，n n n Y Y Y e -=)(*，则0)(0*=Y e ，而且由982.2710011⨯-=-n n Y Y ，78310011⨯-=-n n Y Y 可知， )783982.27(100111-⨯--=---n n n n Y Y Y Y ，即=-⨯-=-⨯-=--)783982.27(1002)()783982.27(1001)()(2*1**n n n Y e Y e Y e ，进而982.27783)783982.27()()(0*100*-=--=Y e Y e ，而31021982.27783-⨯≤-，所以3100*1021)(-⨯=Y ε. 7、供圆程01562=+-x x 的二个根，使它起码具备四位灵验数字（982.27783≈）[解]由78328±=x 与982.27783≈（五位灵验数字）可知，982.55982.2728783281=+=+=x （五位灵验数字）.而018.0982.2728783282=-=-=x ，惟有二位灵验数字，不切合题意.然而是22107863.1982.55178328178328-⨯==+=-=x .8、当N 充分大时，何如供⎰++1211N N dx x？ [解]果为N N dx xN Narctan )1arctan(1112-+=+⎰+，当N 充分大时为二个相近数相减，设)1arctan(+=N α，N arctan =β，则αtan 1=+N ，βtan =N ，进而11)1(1)1(tan tan 1tan tan )tan(2++=++-+=+-=-N N N N N N βαβαβα，果此11arctan 11212++=-=+⎰+N N dx x N Nβα. 9、正圆形的边少约莫为100cm ，应何如丈量才搞使其里积缺面不超出12cm ?[解]由)(2)(])[())((*****2*2**l l l l l εεε='=可知，若央供1))((2**=l ε，则2001100212))(()(*2****=⨯==l l l εε，即边少应谦脚2001100±=l .10、设221gt S =，假定g 是准确的，而对付t 的丈量有1.0±秒的缺面，道明当t 减少时S 的千万于缺面减少，而相对付缺面却缩小.[道明]果为******1.0)()()()(gt t gt t dtdS S ===εεε，***2******51)(2)(21)()()(t t t t g t gt S S S r====εεεε，所以得证.11、序列{}n y 谦脚递推闭系),2,1(1101 =-=-n y y n n ，若41.120≈=y （三位灵验数字），预计到10y 时缺面有多大？那个预计历程宁静吗？[解]设n y 为n y 的近似值，n n n y y y -=)(*ε，则由⎪⎩⎪⎨⎧-==-110210n ny y y 与 ⎩⎨⎧-==-11041.110n n y y y 可知，20*1021)(-⨯=y ε，)(1011---=-n n n n y y y y ，即 )(10)(10)(0*1**y y y n n n εεε==-，进而82100*1010*1021102110)(10)(⨯=⨯⨯==-y y εε，果此预计历程不宁静. 12、预计6)12(-=f，与4.12≈，利用下列公式预计，哪一个得到的截止最佳？6)12(1+，3)223(-，3)223(1+，27099-.[解]果为1*1021)(-⨯=f ε，所以对付于61)12(1+=f ，2417*11*10211054.61021)14.1(6)4.1()(---⨯<⨯=⨯⨯+='=e f f e ，有一位灵验数字； 对付于32)223(-=f ，1112*22*10211012.01021)4.123(6)4.1()(---⨯<⨯=⨯⨯⨯-='=e f f e ，不灵验数字； 对付于33)223(1+=f ，2314*33*10211065.21021)4.123(6)4.1()(---⨯<⨯=⨯⨯⨯+='=e f f e ，有一位灵验数字；对付于270994-=f ，111*44*10211035102170)4.1()(⨯<⨯=⨯⨯='=--e f f e ，不灵验数字. 13、)1ln()(2--=x x x f ，供)30(f 的值.若启仄圆用六位函数表，问供对付数时缺面有多大？若改用另一等价公式)1ln()1ln(22-+-=--x x x x 预计，供对付数时缺面有多大？[解]果为9833.298991302==-（六位灵验数字），4*1021)(-⨯=x ε，所以2442**11*102994.010219833.293011021)13030(1)()()(---⨯=⨯⨯-=⨯⨯---='=x e f f e ，6442**22*108336.010219833.29301102111)()()(---⨯=⨯⨯+=⨯⨯-+-='=x x x e f f e .14、试用消元法解圆程组⎩⎨⎧=+=+2101021102101x x x x ，假定惟有三位数预计，问截止是可稳当？[解]透彻解为110210,110101*********--=-=x x .当使用三位数运算时，得到1,121==x x ，截止稳当.15、已知三角形里积c ab s sin 21=，其中c 为弧度，20π<<c ，且丈量a ，b ，c 的缺面分别为c b a ∆∆∆,,，道明里积的缺面s ∆谦脚cc b b a a s s ∆+∆+∆≤∆. [解]果为c c ab b c a a c b x x f s nk k k ∆+∆+∆=∆∂∂=∆∑=cos 21sin 21sin 21)()(1， 所以cc b b c c c c b b c c c ab cc ab b c a a c b ss ∆+∆+∆≤∆+∆+∆=∆+∆+∆=∆tan sin 21cos 21sin 21sin 21. 第二章 插值法（40-42）1、根据（2.2）定义的范德受止列式，令⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=----nn n n n nn n x x xx xx x x x x x x x V21211020110111),,,,(，道明)(x V n 是n 次多项式，它的根是121,,,-n x x x ，且)())(,,,(),,,,(101101110------=n n n n n x x x x x x x V x x x x V .[道明]由∏∏∏∏-=---=-=-=--⋅=-⋅-=1110111010110)(),,,()()(),,,,(n j j n n n j j n i i j j i n n x x x x x V x x x x x x x x V 可得供证.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 的近似值.[解]若与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 .4、给出 900,cos ≤≤x x 的函数表，步少 )60/1(1='=h ，若函数具备5位灵验数字，钻研用线性插值供x cos 近似值时的总缺面界.[解]设插值节面为h x x x x +=<<010，对付应的x cos 值为10,y y ，函数表值为10,y y ，则由题意可知，5001021-⨯≤-y y ，5111021-⨯≤-y y ，近似线性插值多项式为01011011)(x x x x y x x x x y x L --+--=，所以总缺面为()100101110100100101110100101111,,)()())((2cos )()())((!2)()()()()()()()(x x x x x x y y x x x x y y x x x x x x x x y y x x x x y y x x x x f x L x L x L x f x L x f x R ∈---+---+---=---+---+--''=-+-=-=ξξξ，进而55555201051015100101110100101047.310211094.621102114400121102142110211021))((21))((cos 21)(-------⨯=⨯+⨯⨯=⨯+⨯=⨯+≤--⨯⨯+--⨯⨯+---≤---+---+--≤h x x x x x x x x x x x x x x x x y y x x x x y y x x x x x R ξ.5、设3,2,1,0=+=k kh x x k，供)(max 220x l x x x ≤≤.[解])3)()((max 21)()2()3)()((max))()(())()((max)(max 000300032120231023033030h x x h x x x x h h h h h x x h x x x x x x x x x x x x x x x x x l xx x xx x x x x x x x -----=------=------=≤≤≤≤≤≤≤≤.令)34()383()43()3)()(()(0220302020203000x h hx x x h h x x x h x x h x x h x x x x x f ++-++++-=-----=，则)383()43(23)(202002h h x x x h x x x f ++++-='，进而极值面大概为 hx h h x h h x x h x h x x 37437)43(6)383(12)43(4)43(2002020200±+=±+=++-+±+=，又果为30)20714(271375371374)374(h h h h h x f -=--⨯-⨯-=-+， 30)71420(271357371374)374(h h h h h x f +-=-⨯+⨯+=++， 隐然)374()374(00h x f h x f ++≤-+，所以277710)71420(27121)374(21)(max 3303230+=+=++=≤≤h h h x f h x l x x x . 6、设),,1,0(n j x j=为互同节面，供证：1）),,1,0()(0n k x x l x kn j j k j =≡∑=；2）),,2,1()()(0n k x x l x x knj j k j =≡-∑=；[解]1）果为左侧是k x 的n 阶推格朗日多项式，所以供证创制. 2）设k x y y f )()(-=，则左侧是k x y y f )()(-=的n 阶推格朗日多项式，令x y =，即得供证.7、设[]b a C x f ,)(2∈且0)()(==b f a f ，供证)(max )(81)(max 2x f a b x f b x a b x a ''-≤≤≤≤≤. [解]睹补充题3，其中与0)()(==b f a f 即得.8、正在44≤≤-x 上给出x e x f =)(的等距节面函数表，若用二次插值供x e 的近似值，要使截断缺面不超出610-，问使用函数表的步少h 应与几？[解]由题意可知，设x 使用节面h x x -=10，1x ，h x x +=12举止二次插值，则插值余项为()201112102,)],()[)](([6))()((!3)()(x x h x x x x h x x ex x x x x x f x R ∈+----=---'''=ξξξ，令)()3(3)]()[)](([)(2211221213111h x x x h x x x x h x x x x h x x x f -+-+-=+----=，则)3(63)(22112h x x x x x f -+-='，进而)(x f 的极值面为h x x 331±=，故3932)331()331(33)(max2h h h h x f xx x =-⋅+⋅=≤≤，而 343422739326)(max 6)(20h e h e x f e x R x x x =≤≤≤≤ξ，要使其不超出610-，则有63410273-≤h e ，即22226210472.010389.74863.310243---⨯=⨯≈⨯≤ee h . 9、若n n y 2=，供n y 4∆及n y 4δ.[解]nn n n n n nn n n n n n n n n j jn j j n j jn n y y y y y y j y E j y I E y 22282242322162242624244)1(34)1(24)1(14)1(04)1(4)1(4)1()(123441322314040440444=+⨯-⨯+⨯-⨯=+⨯-⨯+⨯-=⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=-=∆++++++++=-+=-∑∑.22221221122413211204024024021)4(2142121422282242322162242624244)1(34)1(24)1(14)1(04)1(4)1(4)1(4)1()(--------++--++=-+=-=---=+⨯-⨯+⨯-⨯=+⨯-⨯+⨯-=⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-+⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎭⎫ ⎝⎛-=-=∑∑∑n n n n n n n n n n n n n n n n j jn j j n j j j njj jn n y y y y y y j y E j y E Ej y E E y δ. 10、如果)(x f 是m 次多项式，记)()()(x f h x f x f -+=∆，道明)(x f 的k 阶好分)0()(m k x f k ≤≤∆是k m -次多项式，而且0)(=∆+x f l m （l 为正整数）.[道明]对付k 使用数教归纳法可证. 11、道明k k k k k k g f g f g f ∆+∆=∆+1)(. [道明]kk k k k k k k k k k k k k k k k k k k k k k k g f g f g g f g f f g f g f g f g f g f g f g f ∆+∆=-+-=-+-=-=∆++++++++++1111111111)()()(.12、道明∑∑-=+-=∆--=∆11001n k k k n n n k kk f g g f g f g f .[道明]果为01111111111011)()]()([)(g f g f g f f g f f g g g f f g g f f g g fn n n k k k k k n k k k k k k k n k k k k k n k k k n k k k-=-=-+-=∆+∆=∆+∆∑∑∑∑∑-=++-=+++-=+-=+-=，故得证.13、道明：0102y y y n n j j∆-∆=∆∑-=.[道明]01112)(y y y y y n n j j j n j j ∆-∆=∆-∆=∆∑∑-=+-=.14、若n n n n x a x a x a a x f ++++=--1110)( 有n 个分歧真根n x x x ,,,21 ，道明⎩⎨⎧-=-≤≤='-=∑1,20,0)(11n k a n k x f x n nj j k j. [道明]由题意可设∏=-=---=ni i n n n x x a x x x x x x a x f 121)()())(()( ，故∏≠=-='nji i i j n j x x a x f 1)()(，再由好商的本量1战3可知：)!1()(1],,[1)()()1(1111-==-='-=≠==∑∏∑n x a x x x a x x a x xf x n k n n k n nj nj i i i j n k jnj j k j，进而得证.15、道明n 阶均好有下列本量：1）若)()(x cf x F =，则],,,[],,,[1010n n x x x cf x x x F =； 2）若)()()(x g x f x F +=，则],,,[],,,[],,,[101010n n n x x x g x x x f x x x F +=.[道明]1）],,,[)()()()()()(],,,[1000000010n nj nji i i jj nj nji i i jj nj nji i i jj n x x x cf x xx f c x xx cf x xx F x x x F =-=-=-=∑∏∑∏∑∏=≠==≠==≠=.2）],,,[],,,[)()()()()()()()()(],,,[10100000000010n n nj nji i i jj nj nji i i jj nj nji i i jj j nj nji i i jj n x x x g x x x f x xx g x xx f x xx g x f x xx F x x x F +=-+-=-+=-=∑∏∑∏∑∏∑∏=≠==≠==≠==≠=.16、13)(47+++=x x x x f ，供]2,,2,2[71f ，0!80!8)(]2,,2,2[)8(81===ξf f . [解]1!7!7!7)(]2,,2,2[)7(71===ξf f ，]2,,2,2[810 f .17、道明二面三次埃我米特插值余项是()1212)4(3,,!4/)())(()(++∈--=k k k k x x x x x x f x R ξξ，并由此供出分段三次埃我米特插值的缺面限. [解]睹P30与P33，缺面限为k nk f h h '+≤≤0max 278)(ω. 18、XXXXXXXXXX ．19、供一个次数不下于4次的多项式)(x P ，使它谦脚0)0()0(='=P P ，1)1()1(='=P P ，1)2(=P .[解]设1223344)(a x a x a x a x a x P ++++=，则122334234)(a x a x a x a x P +++='，再由0)0()0(='=P P ，1)1()1(='=P P ，1)2(=P 可得：⎪⎪⎪⎩⎪⎪⎪⎨⎧++++==+++='=++++==='===012341234012*********)2(1234)(1)1(1)0(0)0(0a a a a a P a a a a x P a a a a a P a P a P 解得⎪⎪⎪⎪⎪⎪⎩⎪⎪⎪⎪⎪⎪⎨⎧==-===432141234900aa a a a .进而4)3()96(4492341)(2222234-=+-=+-=x x x x x x x x x P .20、设],[)(b a C x f ∈，把[]b a ,分为n 仄分，试构制一个台阶形的整次分段插值函数)(x n ϕ，并道明当∞→n 时，)(x n ϕ正在[]b a ,上普遍支敛到)(x f .[解]令n i x f x f x ii ii x x x x x xi ,,3,2,1,2)(inf)(sup )(11 =+=≤≤≤≤--ϕ.21、设)1/(1)(2x x f +=，正在55≤≤-x 上与10=n ，按等距节面供分段线性插值函数)(x I h ，预计各节面中面处的)(x I h 与)(x f 的值，并预计缺面.[解]由题意可知，1=h ，进而当[]1,+∈k k x x x 时，)(])1(1[1)()1(1)1(1111)(2121211211k k kk k k k k k k k k h x x k h x x k h x x x x k x x x x k l f l f x I -+++-+-=--+++--+=+=++++++.22、供2)(x x f =正在[]b a ,上的分段线性插值函数)(x I h ，并预计缺面.[解]设将[]b a ,区分为少度为h 的小区间b x x x a n =≤≤≤= 10，则当[]1,+∈k k x x x ，1,,2,1,0-=n k 时， 进而缺面为))(())((!2)()(112++--=--''=k k k k x x x x x x x x f x R ξ， 故4))(()(212h x x x x x R k k ≤--=+.23、供4)(x x f =正在[]b a ,上的分段埃我米特插值，并预计缺面. [解]设将[]b a ,区分为少度为h 的小区间b x x x a n =≤≤≤= 10，则当[]1,+∈k k x x x ，1,,2,1,0-=n k 时，)(4)(42121)()(121312113112141121141111++++++++++++++++-⎪⎪⎭⎫⎝⎛--+-⎪⎪⎭⎫ ⎝⎛--⎪⎪⎭⎫ ⎝⎛--+⎪⎪⎭⎫ ⎝⎛--+⎪⎪⎭⎫ ⎝⎛--+⎪⎪⎭⎫ ⎝⎛--='+'++=k k k kk k k k k k k k k k k kk k k kk kk k k k k k k k k k h x x x x x x x x x x x x x x x x x x xx x x x x x x x x x x x x x f f f f x I ββαα，进而缺面为212212)4(2)()()()(!4)()(++--=--=k k k k x x x x x x x x f x R ξ， 故16)()()(42122h x x x x x R k k ≤--=+.24、给定数据表如下：试供三次样条函数)(x S ，并谦脚条件： 1）6868.0)53.0(,0000.1)25.0(='='S S ； 2）0)53.0()25.0(=''=''S S .[解]由05.025.030.00=-=h ，09.030.039.01=-=h ，06.039.045.02=-=h ，08.045.053.03=-=h ，及（8.10）式)1,,1(,,111-=+=+=---n j h h h h h h jj j j jj j j μλ可知，14909.005.009.01011=+=+=h h h λ，5206.009.006.02122=+=+=h h h λ，7408.006.008.03233=+=+=h h h λ，14509.005.005.01001=+=+=h h h μ，5306.009.009.02112=+=+=h h h μ，7308.006.006.03223=+=+=h h h μ，由（8.11）式)1,1(]),[],[(311-=+=+-n j x x f x x f g j j j j j j jμλ可知，7541.2700019279)900768145500477149(3)30.039.05477.06245.014525.030.05000.05477.0149(3])()(145)()(149[3]),[],[(3121201012111011==⨯+⨯⨯=--⨯+--⨯⨯=--+--=+=x x x f x f x x x f x f x x f x x f g μλ.413.2100046332564)6004635390076852(3)39.045.06245.06708.05330.039.05477.06245.052(3])()(53)()(52[3]),[],[(3232312123222122=⨯+⨯=⨯+⨯⨯=--⨯+--⨯⨯=--+--=+=x x x f x f x x x f x f x x f x x f g μλ.0814.27001457140011894634)8004727360046374(3)45.053.06708.07280.07339.045.06245.06708.074(3])()(73)()(74[3]),[],[(3343423234333233==⨯+⨯=⨯+⨯⨯=--⨯+--⨯⨯=--+--=+=x x x f x f x x x f x f x x f x x f g μλ.进而1）矩阵形式为：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⨯-⨯-=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡7871.1413.21112.26868.0730814.2413.20000.11497541.227405325201452321m m m ，解得⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡6570.08278.09078.0321m m m ，进而∑=+=nj j j j j x m x y x S 0)]()([)(βα.2）此为自然鸿沟条件，故862.2500477325.030.05000.05477.03)()(3],[30101100=⨯=--⨯=--⨯==x x x f x f x x f g ；145.2800572345.053.06708.07280.03)()(3],[3111=⨯=--⨯=--⨯==---n n n n n n n x x x f x f x x f g ，矩阵形式为：⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡145.20814.2413.27541.2862.227400732740005325200014521490001243210m m m m m ，不妨解得⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡43210m m m m m ，进而∑=+=nj j j j j x m x y x S 0)]()([)(βα.25、若],[)(2b a C x f ∈，)(x S 是三次样条函数，道明 1）⎰⎰⎰⎰''-''''+''-''=''-''babababadx x S x f x S dx x S x f dx x S dx x f )]()()[(2)]()([)]([)]([222；2）若),,1,0()()(n i x S x f i i ==，式中ix 为插值节面，且b x x x a n =<<<= 10则)]()()[()]()()[()]()()[(a S a f a S b S b f b S dx x S x f x S ba '-'''-'-'''=''-''''⎰.[解]1）⎰⎰⎰⎰⎰⎰⎰⎰''-''=''-''=''-''''+''=''-''''+''-''=''-''''+''-''=''-''''+''-''bababab a ba b ababadxx S dx x f dxx S x f dx x S x f x S x f dxx S x f x S x S x f dxx S x f x S x S x f dxx S x f x S dx x S x f 222222)]([)]([)]([)]([)]()()][()([)]()()}[(2)]()({[)]()()[(2)]()([)]()()[(2)]()([.2）由题意可知，[]b a x A x S ,,)(∈='''，所以)]()()[()]()()[()]()([)]()()[()]()()[()]()([)]()()[()]()()[()()]()([)]}()()[({)]()()[(a S a f a S b S b f b S x S x f A a S a f a S b S b f b S dx x S x f A a S a f a S b S b f b S dxx S x S x f x S x f x S dx x S x f x S b a b ab ab aba'-'''-'-'''=--'-'''-'-'''='-'-'-'''-'-'''=''''-'-'-'''=''-''''⎰⎰⎰.补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并预计插值余项.[解]由1)(000===-e x y y ，111)(-==e x y y 可知，xe x e x x e x x x x x y x x x x y x L )1(1)1(0101011)(111010110101-+=+--=--⨯+--⨯=--+--=---，余项为()1,0),1(2))((!2)()(101∈-=--''=-ξξξx x e x x x x f x R ， 故8141121)1(max max 21)(10101=⨯⨯=-⨯⨯≤≤≤-≤≤x x e x R x ξξ. 2、设4)(x x f =，试利用推格朗日插值余项定理写出以2,1,0,1-为插值节面的三次插值多项式. [解]由插值余项定理，有x x x x x x x x x x x x x x x x x x x f x R 22)1)(2()2)(1()1(!4!4))()()((!4)()(234223210)4(3+--=--=--+=----=ξ，进而x x x x x x x x x R x f x L 22)22()()()(23234433-+=+---=-=.3、设)(x f 正在[]b a ,内有二阶连绝导数，供证：)(max )(81)]()()()([)(max 2x f a b a x a b a f b f a f x f b x a bx a ''-≤---+-≤≤≤≤.[证]果为)()()()(a x ab a f b f a f ---+是以a ，b 为插值节面的)(x f 的线性插值多项式，利用插值多项式的余项定理，得到：))()((21)]()()()([)(b x a x f a x a b a f b f a f x f --''=---+-ξ，进而)(max )(81)(41)(max 21))((max )(max 21)]()()()([)(max 22x f a b a b f b x a x f a x a b a f b f a f x f b x a b a b x a ba b x a ''-=-⋅''=--⋅''≤---+-≤≤≤≤≤≤≤≤≤≤ξξξξ.4、设15)(37++=x x x f ，供好商]2,2[10f ，]2,2,2[210f ，]2,,2,2[710 f 战]2,,2,2[810 f .[解]果为7)1()2(0==f f ，1691252)2()2(371=+⨯+==f f ，167051454)4()2(372=+⨯+==f f ，所以162716912)1()2(]2,2[10=-=--=f f f ，826821691670524)2()4(]2,2[21=-=--=f f f ，27023162826822]2,2[]2,2[]2,2,2[02102121=-=--=f f f ， 1!7!7!7)(]2,,2,2[)7(71===ξf f ，0!80!8)(]2,,2,2[)8(810===ξf f .5、给定数据表：5,4,3,2,1=i ，供4次牛顿插值多项式，并写出插值余项. [解]由好商表可得4次牛顿插值多项式为：)6)(4)(2)(1(1801)4)(2)(1(607)2)(1(65)1(34)6)(4)(2)(1(1801)4)(2)(1(607)2)(1(65)1(34)(4----+------+--=----+------+--=x x x x x x x x x x x x x x x x x x x x x N ，插值余项为()7,1),7)(6)(4)(2)(1(!5)()()5(4∈-----=ξξx x x x x f x R .6、如下表给定函数：4,3,2,1,0=i ，试预计出此列表函数的好分表，并利用牛顿背前插值公式给出它的插值多项式. [解]构制好分表：由好分表可得插值多项式为：32)1(3322)1(332)1()(2020004++=-++=⨯-++=+∆-+∆+=+t t t t t t t t f t t f t f th x N .第三章 函数迫近与预计（80-82）1、（a ）利用区间变更推出区间为[]b a ,的伯恩斯坦多项式；（b ）对付x x f sin )(=正在⎥⎦⎤⎢⎣⎡2,0π上供1次战3次伯恩斯坦多项式并绘出图形，并与相映的马克劳林级数部分战缺面搞出比较. [解]（a ）令t a b a x )(-+=，则[]1,0∈t ，进而伯恩斯坦多项式为∑=-=nk k n x P n k a b f x f B 0)())((),(，其中kn k k x a b x k n x P ---⎪⎪⎭⎫ ⎝⎛=)()(. （b ）令t x 2π=，则[]1,0∈t ，进而伯恩斯坦多项式为∑==nk k n x P n kf x f B 0)()2(),(π，其中k n k k x x k n x P --⎪⎪⎭⎫ ⎝⎛=)2()(π. xx x x x x x f x x f x P kf x f B k k =+⎪⎭⎫⎝⎛-⨯=⨯+⎪⎭⎫ ⎝⎛-⨯=⎪⎭⎫⎝⎛-⎪⎪⎭⎫ ⎝⎛+⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛==∑=202sin 20sin 211)2(201)0()()2(),(010101πππππππ；3223323223223223312213033)533(21)32(4383)2(233)4(23)2(233)2(232sin )2(33sin )2(36sin 20sin )2(33)2()2(23)3()2(13)6()2(03)0()()6(),(x x x x x x x x x x x x x x x x x x x x x x f x x f x x f x x f x P kf x f B k k ----=+-++-=+-+-=⨯+-⨯+-⨯+⎪⎭⎫⎝⎛-⨯=-⎪⎪⎭⎫ ⎝⎛+-⎪⎪⎭⎫ ⎝⎛+-⎪⎪⎭⎫ ⎝⎛+-⎪⎪⎭⎫ ⎝⎛==∑=πππππππππππππππππππππ.2、供证：（a ）当Mx f m ≤≤)(时，M x f B m n ≤≤),(；（b ）当x x f =)(时，x x f B n =),(.[道明]（a ）由∑==nk k n x P nk f x f B 0)()(),(及Mx f m ≤≤)(可知，∑∑∑∑====≤≤≤≤nk k nk k n n k k n k k x P M x MP x f B x mP x P m 0)()(),()()(，而1)]1([)1()(00=-+=-⎪⎪⎭⎫⎝⎛=∑∑=-=nnk k n k nk k x x x x k n x P ，进而得证. （b ）当x x f =)(时，xx x x x x k n k n x x xx k n k n x x k n k n n k x x k n n k f x P n k f x f B n n k k n k n k k n k nk kn k f nk kn k nk k n =-+=----=------=--⨯==-⎪⎪⎭⎫ ⎝⎛==--=--=----=-==-=∑∑∑∑∑110)1(1)1()1(110)0(00)]1([)1()!1(!)!1()1()]!1()1[()!1()!1()1()!(!!)1()()()(),(.3、正在次数不超出6的多项式中，供x x f 4sin )(=正在[]π2,0的最佳普遍迫近多项式.[解]由[]π2,0,4sin ∈x x 可知，14sin 1≤≤-x ，进而最小偏偏好为1，接错面为ππππππππ815,813,811,89,87,85,83,8，此即为6)(H x P ∈的切比雪妇接错面组，进而)(x P 是以那些面为插值节面的推格朗日多项式，可得0)(=x P .4、假设)(x f 正在[]b a ,上连绝，供)(x f 的整次最佳普遍迫近多项式.[解]令)(infx f m bx a ≤≤=，)(sup x f M bx a ≤≤=，则2)(mM x f +=正在[]b a ,上具备最小偏偏好2m M -，进而为整次最佳迫近一次多项式.5、采用常数a ，使得ax x x -≤≤310max 达到极小，又问那个解是可唯一？[解]果为ax x -3是奇函数，所以ax x ax x x x -=-≤≤-≤≤311310max max ，再由定理7可知，当)34(4141333x x T ax x -==-时，坐即43=a ，偏偏好最小.6、供x x f sin )(=正在⎥⎦⎤⎢⎣⎡2,0π上的最佳一次迫近多项式，并预计缺面.[解]由πππ22sin 2sincos )()()(221=--=='=--=x x f ab a f b f a 可得π2arccos2=x ，进而最佳一次迫近多项式为ππππππππππππ2arccos1242)2arccos 21(224)22arccos0(2)]2sin(arccos 0[sin 21)2()]()([2122212--+=-+-=+-++=+-++=x x x x a x a x f a f y 7、供x e x f =)(正在[]1,0上的最佳一次迫近多项式.[解]由101)()()(01212-=--=='=--=e e e e xf a b a f b f a x 可得)1ln(2-=e x ，进而最佳一次迫近多项式为)1ln(212)1()]1ln(21)[1(2)2)1ln(0)(1(][21)2()]()([21)1ln(0212---+-=---+=-+--++=+-++=-e e e x e e x e e e x e e e x a x a x f a f y e .8、怎么样采用r ，使r x x p +=2)(正在[]1,1-上与整偏偏好最小？r 是可唯一？[解]由r r x x p x x +=+=≤≤-≤≤-1)(max )(max 21111，r r x x p x x =+=≤≤-≤≤-)(min )(min 21111可知当与整偏偏好最小时，r r =+1，进而21-=r .另解：由定理7可知，正在[]1,1-上与整偏偏好最小的二次多项式为21)12(21)(21222-=-=x x x T ，进而21-=r .9、设13)(34-+=x x x f ，正在[]1,0上供三次最佳迫近多项式. [解]设所供三次多项式为)(3x P ，则由定理7可知81)188(81)(21)()(2424433+-=+-==-x x x x x T x P x f ，进而893)81()13()81()()(232434243-+=+---+=+--=x x x x x x x x x f x P .10、令[]1,0),12()(∈-=x x T x T n n ，供)(*0x T 、)(*1x T 、)(*2x T 、)(3x T . [解]由[]1,0),12()(∈-=x x T x T n n 可知，令[]1,1,211-∈+=t t x ，则[]1,1),()121(-∈=+t t T t T n n ，进而⎪⎪⎩⎪⎪⎨⎧⎥⎦⎤⎢⎣⎡∈+⎥⎦⎤⎢⎣⎡-∈=1,21),121(21,1),()(00*x x T x x T x T n . 11、试证{})(*x T n 是正在[]1,0上戴权21xx -=ρ的正接多项式.？12、正在[]1,1-上利用插值极小化供x x f arctan )(=的三次近似最佳迫近多项式.[解]由题意可知，插值节面为)3,2,1(,812cos =-k k π，即ππππ87cos ,85cos ,83cos ,81cos 4321====x x x x ，则可供得)(3x L .13、设x e x f =)(正在[]1,1-上的插值极小化近似最佳迫近多项式为)(x L n ，若∞-nL f 有界，道明对付所有1≥n ，存留常数n n βα,，使得)11()()()()(11≤≤-≤-≤++x x T x L x f x T n n n n n βα.[道明]由题意可知，[]1,1),()!1(2)()()(1)1(-∈+=-++ξξx T n f x L x f n n n n ，进而与)!1(2)(min )1(11+=+≤≤-n x f nn x n α，)!1(2)(max )1(11+=+≤≤-n x f n n x n β，则可得供证.14、设正在[]1,1-上543238401653841524381211)(x x x x x x -----=ϕ，试将)(x ϕ落矮到3次多项式并预计缺面.[解]果为x x T x 16545161355-+=，8181244-+=x T x ，所以323232307254510241234096199310241029)16545(3840165)81(3841524381211)(~x x x x x x x x x x ---=-------=ϕ，缺面为0056.040962381384016516138415)(~)(≈=+≤-x x ϕϕ.15、正在[]1,1-利用幂级数项数俭朴供x x f sin )(=的3次迫近多项式，使缺面不超出0.005.[解]果为 ++-+++-=+)!12()1(!5!3sin 1253n x x x x x n n ，与前三项，得到!5!3)(535x x x x L +-=，缺面为0002.0!71)(sin 5≈≤-x L x ，又果为 x x T x 16545161355-+=，所以3次迫近多项式为3333227384383)16545(!51!3sin x x x x x x x +-=-+-=，此时缺面为005.010986.71611201!714<⨯≈⨯+-. 16、)(x f 是[]a a ,-上的连绝奇（奇）函数，道明不管n 是奇数大概奇数，)(x f 的最佳迫近多项式n n H x F ∈)(*也是奇（奇）函数. [解])(x f 的最佳迫近多项式是由切比雪妇多项式得到的，再由切比雪妇多项式的本量4即得.17、供a 、b 使⎰-+202]sin [πdx x b ax 为最小，并与1题及6题的一次迫近多项式缺面做比较. [解]由2120ππ=⎰dx ，8220ππ=⎰dx x ，243202ππ=⎰dx x ，1sin 200==⎰πxdx d ，1cos |)cos (sin 2020201=---==⎰⎰πππxdx x x xdx x d ，可得⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡1124882322a b ππππ，解得⎪⎪⎩⎪⎪⎨⎧=-==-=1148.0)3(86644.0)4(2423ππππb a . 18、],[)(),(1b a C x g x f ∈，定义 （a ）()⎰''=aadx x g x f g f )()(,；（b ）())()()()(,a g a f dx x g x f g f aa+''=⎰. 问它们是可形成内积？[解]（a ）果为()0)(0)]([,0)(2='⇔='=⇒=⎰x f dx x f f f x f ba ，然而反之不可坐，所以不形成内积. （b ）形成内积.19、用许瓦兹不等式（4.5）预计⎰+161dx xx 的上界，并用积分中值定理预计共一积分的上下界，并比较其截止. [解]1961.026113121)131()11()()11(110131012612106≈==+-=+≤+⎰⎰⎰x x dxx dx xdx x x .果为[]1,0,12666∈≤+≤x x xx x ，所以7112141106106106=≤+≤=⎰⎰⎰dx x dx x x dx x . 20、采用a ，使下列积分与最小值：⎰-1022)(dx ax x ，⎰--112dx ax x .[解]481)45(51512131)2()(22142321022+-=+-=+-=-⎰⎰a a a dx x a ax x dx ax x ，进而45=a .当0=a 时，12121100111112=+=+-==-⎰⎰⎰⎰---xdx xdx dx x dx ax x ，当0≠a 时，由02=-ax x ，可得接面为ax 1=，若1>a ，则1323121316161)2131()2131()3121()2131()()()(222012310321123012102112112>+=++++-=-+-+-=-+-+-=----⎰⎰⎰⎰a a a aa a x ax ax x x ax dxx ax dx ax x dx x ax dx ax x a aa a，若01>≥a ，则1)2131()3121()()(012102112=----=-+-=-⎰⎰⎰--a a dx x ax dx ax x dx ax x .共理可知，当01<≤-a 时，1112=-⎰-dx ax x ，当1-<a 时，1112>-⎰-dx ax x ，进而当1≤a 时，积分博得最小.21、设{}x span ,11=ϕ，{}1011002,x x span =ϕ，分别正在21,ϕϕ上供一元素，使其为]1,0[2C x ∈的最佳仄圆迫近，并比较其截止.[解]由1110=⎰dx ，2110=⎰xdx ，31102=⎰dx x ，41103=⎰dx x 可知，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡41313121211b a ，解得⎪⎩⎪⎨⎧=-=161b a ，即正在1ϕ上为⎪⎭⎫ ⎝⎛-1,61. 由201110100100=⋅⎰dx x x ，202110101100=⋅⎰dx x x ，203110101101=⋅⎰dx x x ，1031102100=⋅⎰dx x x ，1041102101=⋅⎰dx x x 可知， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡104110312031202120212011b a ，解得⎪⎪⎩⎪⎪⎨⎧-≈⨯⨯⨯-=≈⨯⨯⨯=148.37510310420320298243.37510410320220199b a ，即正在2ϕ上为()148.375,243.375-.22、x x f =)(正在[]1,1-上，供正在{}421,,1x x span =ϕ上的最佳仄圆迫近.[解]由1100111=+-=⎰⎰⎰--xdx xdx dx x ，2113013112=+-=⎰⎰⎰--dx x dx x dx x x ，。
数值分析课程第五版课后习题答案(李庆扬等)(OCR)
根是ｘ，，２…，ｘ－，且Ｖ。ｘ，ｘ…·，ｘ）＝Ｖ，Ｃｘ６，ｘ…·）（ｘ－ｘ）…（ｘ－ｘ）。
Ｖ，（ｘｏ，ｘ，…ｘ－ｘ）＝１１】 －ｘ，）用ａ－ｘ，）
［证明］由
可得求证。
＝Ｖ，（Ｃｘ８，ｘ，…，ｘＸ））１１（ｘ－ｘ）
２、当ｘ＝１－１，２时，ｆ（ｘ）＝０，－３．４，求ｆ（ｘ）的二次插值多项式。
Ｌ，（ｘ）＝ｙ％（（ｘｘ６－－ｘｘ，）（（ｘｘ－２ｘ－ｘ２２））
ｙ＝ｆ（ｘ）＝ｆ０．５）＝－０．６９３１４７，ｙ２＝ｆ（ｘ）＝ｆ（０．６）＝－０．５１０８２６，则
Ｌ２（ｘ）＝ｙ。 （ｘ－ｘ）（ｘ－ｘ２）
（ｘ６－ｘ）ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ）
（ｘ－ｘ）（ｘ－ｘ２）
（ｘ－ｘｏ）（ｘ－ｘ） （ｘ２－ｘｏ）（ｘ２－ｘ）
＝－０．９１６２９１×．（０（．ｘ４－０－．０５．）５（）ｘ（－００．．６４）－０．６－．
３０—＋２—９．ｘ９５８３ｘ３１￣０２＇＝０．８３３６×１０４
１４、试用消元法解方程ｘ组１＋１０＂ｘ＝１００
ｘ＋ｘ２＝２
，假定只有三位数计算，问结果是否
可靠？
［解］精确解为ｘ１＝０１００－＊１ １０＂－２ ，当使用三位数运 算时，得到
ｘ ＝１，ｘ２＝１，结果可靠。
１５、已知三角形面积ｓ＝ｓ去＝ ａｂｓｉｎｃ，其中ｃ为弧度，０＜ｃ＜ 且测量ａ，ｂ，ｃ
位有效数字；ｘ＝５６．４３０有５位有效数字；ｘ＝７×１０有２位有效数字。 ４、利用公式（３．３）求下列各近似值的误差限，其中ｘ，ｘ；，ｘ，ｘ；均为第３题所给
的数。
（１）ｘ＋ｘ２＋ｘ：
ｅ（ｘ＋ｘ写＋ｘ）＝＞
［解］
Ｅ（ｘ）＝Ｅ（ｘ）＋Ｅ（ｘ）＋Ｅ（ｘ；）
３＋ｔｘ１０＝１．０５×１０３
（２）ｘｘｘ；
数值分析(第五版)计算实习题第四章作业
第四章：1、（1）：复合梯形建立m文件：function t=natrapz(fname,a,b,n)h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b); f=feval(fname,a+h:h:b-h+0.001*h); t=h*(0.5*(fa+fb)+sum(f));输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,10)输出：ans =-0.417062831779470输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,100)输出：ans =-0.443117908008157输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,1000)输出：ans =-0.444387538997162复合辛普森建立m文件：function t=comsimpson(fname,a,b,n)h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b);f1=feval(fname,a+h:h:b-h+0.001*h);f2=feval(fname,a+h/2:h:b-h+0.001*h);t=h/6*(fa+fb+2*sum(f1)+4*sum(f2));输入：>> syms x>> f=inline('sqrt(x).*log(x);');>> format long；>>comsimpson(f,eps,1,10)输出：ans =-0.435297890074689输入：>>syms x>>f=inline('sqrt(x).*log(x);');>>comsimpson(f,eps,1,100)输出：ans =-0.444161178415673输入：>>syms x>>f=inline('sqrt(x).*log(x);');>>comsimpson(f,eps,1,1000)输出：ans =-0.444434117614180（2）龙贝格建立m文件：function [RT,R,wugu,h]=Romberg(fun,a,b,wucha,m) %RT是龙贝格积分表%R是数值积分值%wugu是误差估计%h是最小步长%fun是被积函数%a b是积分下、上限%m是龙贝格积分表中行最大数目%wucha是两次相邻迭代值的绝对误差限n=1;h=b-a;wugu=1;x=a;k=0;RT=zeros(4,4);RT(1,1)=h*(feval(fun,a)+feval(fun,b))/2;while((wugu>wucha)&(k<m)|(k<4))k=k+1;h=h/2;s=0;for j=1:nx=a+h*(2*j-1);s=s+feval(fun,x);endRT(k+1,1)=RT(k,1)/2+h*s;n=2*n;for i=1:kRT(k+1,i+1)=((4^i)*RT(k+1,i)-RT(k,i))/(4^i-1);endwugu=abs(RT(k+1,k)-RT(k+1,k+1));endR=RT(k+1,k+1);输入：>>fun=inline('sqrt(x).*log(x)');>> [RT,R,wugu,h]=Romberg(fun,eps,1,1e-5,13)输出：RT =1 至5 列-0.000000268546145 0 0 0-0.245064670140209 -0.326752804004897 0 0-0.358104125949240 -0.395783944552250 -0.400386020588741 0 0-0.408090073087781 -0.424752055467295 -0.426683262861631 -0.427100679405645 0-0.429474601629505 -0.436602777810080 -0.437392825966266 -0.437562819031419 -0.437603847029951-0.438389494461832 -0.441361125405941 -0.441678348578999 -0.441746372747455 -0.4417627788404596 列-0.441766844267449R =-0.441766844267449wugu =4.065426989774412e-06h =0.031250000000000（3）自适应辛普森输入：>> f=inline('sqrt(x).*log(x)');>> q=quad(f,0,1,1e-4)输出：q =-0.4439755729517282.（1）复合辛普森建立m文件function q=combinesimpson2(F,x0,a,b,n)%复合Simpson多元求积公式%F—被积函数%x0—被积函数自变量%[a,b]积分区间%n—区间份数x=linspace(a,b,n+1);q=0;for k=1:nq=q+subs(F,x0,x(k))+4*subs(F,x0,(x(k)+x(k+1))/2)+subs(F,x0,x(k+1)); endq=q*(b-a)/n/6;输入：>> clear>> syms x y;>> F=exp(-x.*y);>> s=combinesimpson2(combinesimpson2(F,'x',0,1,4),'y',0,1,4)输出：s =exp(-1)/576 + exp(-1/2)/144 + exp(-1/4)/72 + exp(-3/4)/144 + exp(-1/8)/36 +exp(-3/8)/36 + exp(-5/8)/72 + exp(-7/8)/72 + (5*exp(-1/16))/144 + exp(-3/16)/24 + exp(-5/16)/36 + exp(-7/16)/36 + exp(-9/16)/144 + exp(-1/32)/36 + exp(-3/32)/18 + exp(-5/32)/36 + exp(-7/32)/36 + exp(-9/32)/36 + exp(-15/32)/36 + exp(-21/32)/36 + exp(-1/64)/36 + exp(-3/64)/18 + exp(-5/64)/18 + exp(-7/64)/18 + exp(-9/64)/36 + exp(-15/64)/18 + exp(-21/64)/18 + exp(-25/64)/36 + exp(-35/64)/18 + exp(-49/64)/36 + 47/576>> double(s)ans =0.796599967946203高斯求积公式function q=gaussquad(F,x0,a,b,n)%Gauss求积公式%F—被积函数%x0—被积函数自变量%[a,b]积分区间%n—节点个数syms t;F=subs(F,x0,(b-a)/2*t+(a+b)/2);[x,A]=gausspoints(n);q=(b-a)/2*sum(A.*subs(F,t,x));输入：>> clear>> syms x y;F=exp(-x.*y);>> s=gaussquad(gaussquad(F,x,0,1,4),y,0,1,4)输出：s =0.7966（2）复合辛普森输入：>> syms x y;>> f=exp(-x.*y);>> s=combinesimpson2(combinesimpson2(f,y,0,sqrt(1-x^2),4),x,0,1,4)输出：s =(3^(1/2)*(exp(-3^(1/2)/4) + 2*exp(-3^(1/2)/8) + 2*exp(-3^(1/2)/16) + 2*exp(-(3*3^(1/2))/16) + 4*exp(-3^(1/2)/32) + 4*exp(-(3*3^(1/2))/32) + 4*exp(-(5*3^(1/2))/32) + 4*exp(-(7*3^(1/2))/32) + 1))/576 + (7^(1/2)*(exp(-(3*7^(1/2))/16) + 2*exp(-(3*7^(1/2))/32) + 2*exp(-(3*7^(1/2))/64) + 2*exp(-(9*7^(1/2))/64) + 4*exp(-(3*7^(1/2))/128) + 4*exp(-(9*7^(1/2))/128) + 4*exp(-(15*7^(1/2))/128) + 4*exp(-(21*7^(1/2))/128) + 1))/1152 + (15^(1/2)*(exp(-15^(1/2)/16) + 2*exp(-15^(1/2)/32) + 2*exp(-15^(1/2)/64) + 2*exp(-(3*15^(1/2))/64) + 4*exp(-15^(1/2)/128) + 4*exp(-(3*15^(1/2))/128) + 4*exp(-(5*15^(1/2))/128) + 4*exp(-(7*15^(1/2))/128) + 1))/1152 + (15^(1/2)*(exp(-(7*15^(1/2))/64) + 2*exp(-(7*15^(1/2))/128) + 2*exp(-(7*15^(1/2))/256) + 2*exp(-(21*15^(1/2))/256) + 4*exp(-(7*15^(1/2))/512) + 4*exp(-(21*15^(1/2))/512) + 4*exp(-(35*15^(1/2))/512) + 4*exp(-(49*15^(1/2))/512) + 1))/1152 + (39^(1/2)*(exp(-(5*39^(1/2))/64) + 2*exp(-(5*39^(1/2))/128) + 2*exp(-(5*39^(1/2))/256) + 2*exp(-(15*39^(1/2))/256) + 4*exp(-(5*39^(1/2))/512) + 4*exp(-(15*39^(1/2))/512) + 4*exp(-(25*39^(1/2))/512) + 4*exp(-(35*39^(1/2))/512) + 1))/1152 + (55^(1/2)*(exp(-(3*55^(1/2))/64) + 2*exp(-(3*55^(1/2))/128) + 2*exp(-(3*55^(1/2))/256) + 2*exp(-(9*55^(1/2))/256) + 4*exp(-(3*55^(1/2))/512) + 4*exp(-(9*55^(1/2))/512) + 4*exp(-(15*55^(1/2))/512) + 4*exp(-(21*55^(1/2))/512) + 1))/1152 + (63^(1/2)*(exp(-63^(1/2)/64) + 2*exp(-63^(1/2)/128) + 2*exp(-63^(1/2)/256) + 2*exp(-(3*63^(1/2))/256) + 4*exp(-63^(1/2)/512) + 4*exp(-(3*63^(1/2))/512) + 4*exp(-(5*63^(1/2))/512) + 4*exp(-(7*63^(1/2))/512) + 1))/1152 + 1/24>> double(s)ans =0.670113633359095。
数值分析课程第五版课后习题答案(李庆扬等)
第一章 绪论（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 =ε，从而n x 的误差为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、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第四章：1、（1）：复合梯形建立m文件：function t=natrapz(fname,a,b,n)h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b); f=feval(fname,a+h:h:b-h+0.001*h); t=h*(0.5*(fa+fb)+sum(f));输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,10)输出：ans =-0.417062831779470输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,100)输出：ans =-0.443117908008157输入：>> syms x>> f=inline('sqrt(x).*log(x);'); >> natrapz(f,eps,1,1000)输出：ans =-0.444387538997162复合辛普森建立m文件：function t=comsimpson(fname,a,b,n)h=(b-a)/n;fa=feval(fname,a);fb=feval(fname,b);f1=feval(fname,a+h:h:b-h+0.001*h);f2=feval(fname,a+h/2:h:b-h+0.001*h);t=h/6*(fa+fb+2*sum(f1)+4*sum(f2));输入：>> syms x>> f=inline('sqrt(x).*log(x);');>> format long；>>comsimpson(f,eps,1,10)输出：ans =-0.435297890074689输入：>>syms x>>f=inline('sqrt(x).*log(x);');>>comsimpson(f,eps,1,100)输出：ans =-0.444161178415673输入：>>syms x>>f=inline('sqrt(x).*log(x);');>>comsimpson(f,eps,1,1000)输出：ans =-0.444434117614180（2）龙贝格建立m文件：function [RT,R,wugu,h]=Romberg(fun,a,b,wucha,m) %RT是龙贝格积分表%R是数值积分值%wugu是误差估计%h是最小步长%fun是被积函数%a b是积分下、上限%m是龙贝格积分表中行最大数目%wucha是两次相邻迭代值的绝对误差限n=1;h=b-a;wugu=1;x=a;k=0;RT=zeros(4,4);RT(1,1)=h*(feval(fun,a)+feval(fun,b))/2;while((wugu>wucha)&(k<m)|(k<4))k=k+1;h=h/2;s=0;for j=1:nx=a+h*(2*j-1);s=s+feval(fun,x);endRT(k+1,1)=RT(k,1)/2+h*s;n=2*n;for i=1:kRT(k+1,i+1)=((4^i)*RT(k+1,i)-RT(k,i))/(4^i-1);endwugu=abs(RT(k+1,k)-RT(k+1,k+1));endR=RT(k+1,k+1);输入：>>fun=inline('sqrt(x).*log(x)');>> [RT,R,wugu,h]=Romberg(fun,eps,1,1e-5,13)输出：RT =1 至5 列-0.000000268546145 0 0 0-0.245064670140209 -0.326752804004897 0 0-0.358104125949240 -0.395783944552250 -0.400386020588741 0 0-0.408090073087781 -0.424752055467295 -0.426683262861631 -0.427100679405645 0-0.429474601629505 -0.436602777810080 -0.437392825966266 -0.437562819031419 -0.437603847029951-0.438389494461832 -0.441361125405941 -0.441678348578999 -0.441746372747455 -0.4417627788404596 列-0.441766844267449R =-0.441766844267449wugu =4.065426989774412e-06h =0.031250000000000（3）自适应辛普森输入：>> f=inline('sqrt(x).*log(x)');>> q=quad(f,0,1,1e-4)输出：q =-0.4439755729517282.（1）复合辛普森建立m文件function q=combinesimpson2(F,x0,a,b,n)%复合Simpson多元求积公式%F—被积函数%x0—被积函数自变量%[a,b]积分区间%n—区间份数x=linspace(a,b,n+1);q=0;for k=1:nq=q+subs(F,x0,x(k))+4*subs(F,x0,(x(k)+x(k+1))/2)+subs(F,x0,x(k+1)); endq=q*(b-a)/n/6;输入：>> clear>> syms x y;>> F=exp(-x.*y);>> s=combinesimpson2(combinesimpson2(F,'x',0,1,4),'y',0,1,4)输出：s =exp(-1)/576 + exp(-1/2)/144 + exp(-1/4)/72 + exp(-3/4)/144 + exp(-1/8)/36 +exp(-3/8)/36 + exp(-5/8)/72 + exp(-7/8)/72 + (5*exp(-1/16))/144 + exp(-3/16)/24 + exp(-5/16)/36 + exp(-7/16)/36 + exp(-9/16)/144 + exp(-1/32)/36 + exp(-3/32)/18 + exp(-5/32)/36 + exp(-7/32)/36 + exp(-9/32)/36 + exp(-15/32)/36 + exp(-21/32)/36 + exp(-1/64)/36 + exp(-3/64)/18 + exp(-5/64)/18 + exp(-7/64)/18 + exp(-9/64)/36 + exp(-15/64)/18 + exp(-21/64)/18 + exp(-25/64)/36 + exp(-35/64)/18 + exp(-49/64)/36 + 47/576>> double(s)ans =0.796599967946203高斯求积公式function q=gaussquad(F,x0,a,b,n)%Gauss求积公式%F—被积函数%x0—被积函数自变量%[a,b]积分区间%n—节点个数syms t;F=subs(F,x0,(b-a)/2*t+(a+b)/2);[x,A]=gausspoints(n);q=(b-a)/2*sum(A.*subs(F,t,x));输入：>> clear>> syms x y;F=exp(-x.*y);>> s=gaussquad(gaussquad(F,x,0,1,4),y,0,1,4)输出：s =0.7966（2）复合辛普森输入：>> syms x y;>> f=exp(-x.*y);>> s=combinesimpson2(combinesimpson2(f,y,0,sqrt(1-x^2),4),x,0,1,4)输出：s =(3^(1/2)*(exp(-3^(1/2)/4) + 2*exp(-3^(1/2)/8) + 2*exp(-3^(1/2)/16) + 2*exp(-(3*3^(1/2))/16) + 4*exp(-3^(1/2)/32) + 4*exp(-(3*3^(1/2))/32) + 4*exp(-(5*3^(1/2))/32) + 4*exp(-(7*3^(1/2))/32) + 1))/576 + (7^(1/2)*(exp(-(3*7^(1/2))/16) + 2*exp(-(3*7^(1/2))/32) + 2*exp(-(3*7^(1/2))/64) + 2*exp(-(9*7^(1/2))/64) + 4*exp(-(3*7^(1/2))/128) + 4*exp(-(9*7^(1/2))/128) + 4*exp(-(15*7^(1/2))/128) + 4*exp(-(21*7^(1/2))/128) + 1))/1152 + (15^(1/2)*(exp(-15^(1/2)/16) + 2*exp(-15^(1/2)/32) + 2*exp(-15^(1/2)/64) + 2*exp(-(3*15^(1/2))/64) + 4*exp(-15^(1/2)/128) + 4*exp(-(3*15^(1/2))/128) + 4*exp(-(5*15^(1/2))/128) + 4*exp(-(7*15^(1/2))/128) + 1))/1152 + (15^(1/2)*(exp(-(7*15^(1/2))/64) + 2*exp(-(7*15^(1/2))/128) + 2*exp(-(7*15^(1/2))/256) + 2*exp(-(21*15^(1/2))/256) + 4*exp(-(7*15^(1/2))/512) + 4*exp(-(21*15^(1/2))/512) + 4*exp(-(35*15^(1/2))/512) + 4*exp(-(49*15^(1/2))/512) + 1))/1152 + (39^(1/2)*(exp(-(5*39^(1/2))/64) + 2*exp(-(5*39^(1/2))/128) + 2*exp(-(5*39^(1/2))/256) + 2*exp(-(15*39^(1/2))/256) + 4*exp(-(5*39^(1/2))/512) + 4*exp(-(15*39^(1/2))/512) + 4*exp(-(25*39^(1/2))/512) + 4*exp(-(35*39^(1/2))/512) + 1))/1152 + (55^(1/2)*(exp(-(3*55^(1/2))/64) + 2*exp(-(3*55^(1/2))/128) + 2*exp(-(3*55^(1/2))/256) + 2*exp(-(9*55^(1/2))/256) + 4*exp(-(3*55^(1/2))/512) + 4*exp(-(9*55^(1/2))/512) + 4*exp(-(15*55^(1/2))/512) + 4*exp(-(21*55^(1/2))/512) + 1))/1152 + (63^(1/2)*(exp(-63^(1/2)/64) + 2*exp(-63^(1/2)/128) + 2*exp(-63^(1/2)/256) + 2*exp(-(3*63^(1/2))/256) + 4*exp(-63^(1/2)/512) + 4*exp(-(3*63^(1/2))/512) + 4*exp(-(5*63^(1/2))/512) + 4*exp(-(7*63^(1/2))/512) + 1))/1152 + 1/24>> double(s)ans =0.670113633359095。