最新第六章习题答案-数值分析

第六章习题解答1、设函数01(),(),,()n x x x φφφ 在[,]a b 上带权()x ρ正交，试证明{}()nj j x φ=是线性无关组。

证明：设0()nj jj l x φ==∑，两端与01()(,,,)kx k n φ= 作内积，由()jx φ的正交性可知，200(),()((),())((),())()()n n b k j j j k j k k k k k a j j x l x l x x l x x l x x dx φφφφφφρφ==⎛⎫==== ⎪⎝⎭∑∑⎰， 于是有001(,,,)k l k n == ，即{}()nj j x φ=是线性无关组。

2、试确定系数,a b 的值使22(()cos )ax b x dx π+-⎰达到最小。

解：定义02,[,]f g C π∈上的内积为20fgdx π⎰，取011(),()x x x ϕϕ==，()s x ax b =+，()cos f x x =，则法方程为0001010111(,)(,)(,)(,)(,)(,)f a f b ϕϕϕϕϕϕϕϕϕϕ⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 其中()2000112,dx ππϕϕ=⨯=⎰，()2201018,xdx ππϕϕ=⨯=⎰，()3211024,x xdx ππϕϕ=⨯=⎰，()2001,cos f xdx πϕ==⎰，()21012,cos f x xdx ππϕ==-⎰，于是方程组为22312812824a b πππππ⎛⎫⎛⎫ ⎪⎛⎫ ⎪ ⎪= ⎪ ⎪- ⎪ ⎪⎝⎭⎝⎭ ⎪⎝⎭，解之得1158506644.,.a b ==-。

3、已知函数11()(,)f x x =∈-，试用二类Chebyshev 多项式()n U x 构造此函数的二次最佳平方逼近元。

解：法一、取20121(),(),(),x x x x x ϕϕϕ===()()()00112222235,,,,,ϕϕϕϕϕϕ===，()()()011202203,,,,ϕϕϕϕϕϕ===，同时由二类Chebyshev 多项式的性质知 ()()()11101211028,,,,,f f f x ππϕϕϕ---======⎰⎰⎰于是可得法方程为0122203220003220835c c c ππ⎛⎫⎛⎫⎪ ⎪⎛⎫ ⎪ ⎪ ⎪⎪= ⎪ ⎪ ⎪⎪ ⎪ ⎪⎝⎭ ⎪ ⎪⎝⎭ ⎪⎝⎭，解之得0121.0308,0,0.7363c c c ===-， 于是()f x 的二次最佳逼近元是2001122() 1.03080.7363x c c c x ϕϕϕϕ=++=-法一、二类Chebyshev 多项式2012()1,()2,()41U x U x x U x x ===-，取内积权函数()()x f x ρ==，于是11200114(,)(1)3f U fU dx x dx ρ--==-=⎰⎰，1121111(,)2(1)0f U fU dx x x dx ρ--==-=⎰⎰，112222114(,)(41)(1)15f U fU dx x x dx ρ--==--=-⎰⎰ 由()n U x 正交性及(,)2n n U U π=可得0000(,)8(,)3f U c U U π==，1111(,)0(,)f U c U U ==，2222(,)8(,)15f U c U U π==-， 于是()f x 的二次最佳逼近元为001122()x c U c U c U ϕ=++=21632515x ππ- 4、设012{(),(),()}L x L x L x 是定义于[0,)+∞上关于权函数()xx eρ-=的首项系数为1的正交多项式组，若已知01()1,()1L x L x x ==-，试求出二次多项式2()L x 。

1.01
1.01
1
0.66
0.995
0.66
1.17
2
0.67
1.17
0.553333
1.223333
3
0.553333
1.165
0.517778
1.241111
4
0.556667
1.223333
0.505926
1.247037
5
0.517778
1.221667
0.501975
1.249012
6
0.518889
1-4.求方程x2-56x+1=0的两个根,使它们至少具有四 位有效数字 ( 7832.798)2.
解 x1=28+27.982=55.982,x2=1/x1=0.017863
精选课件
2
二.习题2 (第50页)
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
2.11
3-4.用J迭代法和G-S迭代法求解方程组Ax=b,其中
A 1 1
问取何值时这两种迭代法是收敛的? 解 J迭代法和G-S迭代法的迭代矩阵分别为
B 0 0
G 0 0
2
易得:(B)=||,(G)=2.故当||<1时两种方法都收敛.
3-7.给定方程组
精选课件
16
(1 ) 3 x 1 x 1 2 2 x x 22 3 4
x2 0
1 4 1
x3
0
1 4 1x4 0
1精选课4件x5 200
8
解
4 1 1 4 1
44
111 444
11

第一章 绪论姓名 学号 班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。

1 若误差限为5105.0-⨯,那么近似数0.003400有几位有效数字？（有效数字的计算） 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。

2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算） 解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（3.14109 , 3.14209）之间的任意数，都具有4位有效数字。

3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字？（有效数字的计算）解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。

2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。

4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算） 解：已知δ=-**xx x ，则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。

第六章课后习题解答(1)()()123(1)()213(1)()()312(01.21125551154213351010(1,1,1),17( 4.0000186,2.99999k k k k k k k k k Tx x x x x x x x x x x+++ìïï=---ïïïïïï=-+íïïïïï=-++ïïïî==-(17)解：(a )因系数矩阵按行严格对角占优，故雅可比法与高斯-塞德尔均收敛。

（b )雅可比法的迭代格式为取迭代到次达到精度要求(1)()()123(1)(1)()213(1)(1)(1)312(0)(8)15,2.0000012)21125551154213351010(1,1,1),8( 4.0000186,2.9999915,2.0000012)Tk k k k k k k k k TTx x x x x x x x x x++++++-ìïï=---ïïïïïï=-+íïïïïï=-++ïïïî==-高斯塞德尔法的迭代格式为x 取迭代到次达到精度要求1212:00.40.4.0.400.80.40.80||(0.8)(0.80.32)()1.09282031,00.40.4()00.160.6400.0320.672DL U I BD L U l l l l--骣--÷ç÷ç÷ç÷ç÷=+=--ç÷ç÷÷ç÷ç÷--÷ç桫-=-+-=>-æ--çççç=-=-ççççèlJJJS解(a )雅可比法的迭代矩阵B（）BB故雅可比迭代法不收敛高斯塞德尔法迭代矩阵131()||||0.81022101220||022023002SJBDL U I BD L Ul l¥--ö÷÷÷÷÷÷÷÷÷÷ç÷ø?<骣-÷ç÷ç÷ç÷ç÷=+=--ç÷ç÷÷ç÷ç÷--ç÷桫-=骣-÷ç÷ç÷ç÷ç÷=-=-ç÷ç÷÷ç÷ç÷ç桫llSJJ SB故高斯-塞德尔迭代法收敛。

7、计算的近似值，取。

利用以下四种计算格式，试问哪一种算法误差最小。

〔1〕〔2〕〔3〕〔4〕解：计算各项的条件数由计算知，第一种算法误差最小。

解：在计算机上计算该级数的是一个收敛的级数。

因为随着的增大，会出现大数吃小数的现象。

9、通过分析浮点数集合F=〔10，3，-2，2〕在数轴上的分布讨论一般浮点数集的分布情况。

10、试导出计算积分的递推计算公式，用此递推公式计算积分的近似值并分析计算误差，计算取三位有效数字。

解：此算法是数值稳定的。

第二章习题解答1.〔1〕 R n×n中的子集“上三角阵〞和“正交矩阵〞对矩阵乘法是封闭的。

〔2〕R n×n中的子集“正交矩阵〞，“非奇异的对称阵〞和“单位上〔下〕三角阵〞对矩阵求逆是封闭的。

设A是ｎ×ｎ的正交矩阵。

证明A-1也是ｎ×ｎ的正交矩阵。

证明：〔2〕A是ｎ×ｎ的正交矩阵∴A A-1 =A-1A=E 故〔A-1〕-1=A∴A-1〔A-1〕-1=〔A-1〕-1A-1 =E 故A-1也是ｎ×ｎ的正交矩阵。

设A是非奇异的对称阵，证A-1也是非奇异的对称阵。

A非奇异∴A可逆且A-1非奇异又A T=A ∴〔A-1〕T=〔A T〕-1=A-1故A-1也是非奇异的对称阵设A是单位上〔下〕三角阵。

证A-1也是单位上〔下〕三角阵。

证明：A是单位上三角阵，故|A|=1，∴A可逆，即A-1存在，记为〔b ij〕n×n由A A-1 =E，那么〔其中 j＞i时，〕故b nn=1, b ni=0 (n≠j)类似可得，b ii=1 (j=1…n) b jk=0 (k＞j)即A-1是单位上三角阵综上所述可得。

R n×n中的子集“正交矩阵〞，“非奇异的对称阵〞和“单位上〔下〕三角阵〞对矩阵求逆是封闭的。

2、试求齐次线行方程组Ax=0的根底解系。

A=解：A=～～～故齐次线行方程组Ax=0的根底解系为，3.求以下矩阵的特征值和特征向量。

1、解：将)(x V n 按最后一行展开，即知)(x V n 是n 次多项式。

由于ni i inn n n n i n x x x x x x x x x x V ...1...1 (1))(21110200---=，.1,...,1,0-=n i故知0)(=i n x V ，即110,...,,-n x x x 是)(x V n 的根。

又)(x V n 的最高次幂nx 的系数为)(...1...1..................1),...,,(101121112222102001101j n i j i n n n n n n n n n n n x x x x x x x x x x x x x x V -==∏-≤<≤-----------。

故知).)...()()(,...,,()(1101101------=n n n n x x x x x x x x x V x V6、解：（1）设.)(k x x f =当n k ,...,1,0=时，有.0)()1(=+x f n对)(x f 构造Lagrange 插值多项式，),()(0x l x x L j nj k j n ∑==其0)()!1()()()()(1)1(=+=-=++x w n f x L x F x R n n n n ξ，ξ介于j x 之间，.,...,1,0n j =故),()(x L x f n =即.,...,1,0,)(0n k xx l x kjnj k j ==∑=特别地，当0=k 时，10)(=∑=nj x j l。

（2）0)()1(1)()1()()(0000=-=⎪⎪⎭⎫ ⎝⎛-⎪⎪⎭⎫ ⎝⎛-=--=-===∑∑∑∑k j j i j i k j ki i j ii k j nj ki i j knj j x x x x i k x l x x i k x l x x ）利用（。

7、证明：以b a ,为节点进行线性插值，得)()()(1b f ab ax a f b a b x x P --+--=因0)()(==b f a f ，故0)(1=x P 。

第六章习题解答2、利用梯形公式和Simpson 公式求积分21ln xdx ⎰的近似值，并估计两种方法计算值的最大误差限。

解：①由梯形公式：21ln 2()[()()][ln1ln 2]0.3466222b a T f f a f b --=+=+=≈ 最大误差限3''2()111()()0.0833********T b a R f f ηη-=-=≤=≈ 其中，(1,2)η∈ ②由梯形公式：13()[()4()()][ln14ln()ln 2]0.38586262b a b a S f f a f f b -+=++=++≈ 最大误差限5(4)4()66()()0.0021288028802880S b a R f f ηη-=-=≤≈，其中，(1,2)η∈。

4、推导中点求积公式3''()()()()()()224baa b b a f x dx b a f f a b ξξ+-=-+<<⎰证明：构造一次函数P （x ），使'',()()2222a b a b a b a b P f P f ++++⎛⎫⎛⎫== ⎪ ⎪⎝⎭⎝⎭则，易求得'()()()()222a b a b a bP x f x f +++=-+ 且'()()()()222bbaa a ba b a b P x dx f x f dx +++⎡⎤=-+⎢⎥⎣⎦⎰⎰0()()()22ba ab a bf dx b a f ++=+=-⎰，令()b a P x dx Z =⎰现分析截断误差：令'()()()()()()-()222a b a b a b r x f x P x f x f x f +++=-=-- 由'''()()()2a b r x f x f +=-易知2a b x +=为()r x 的二重零点，所以可令2()()()2a b r x x x ϕ+=-，构造辅助函数2()()()()()2a b K t f t P t x t ϕ+=---，则易知： ()02a b K x K +⎛⎫== ⎪⎝⎭其中2a b t +=为二重根()K t ∴有三个零点 ∴由罗尔定理，存在''''''()(,)()0()2()0()2f a b K f K x K x ηηηη∈=-=∴=使即从而可知''2()()()()()22f a b r x f x P x x η+=-=- ∴截断误差[]''2()()()()()()()22b bb baaa af a b R f f x dx Z f x P x dx r x dx x dx η+=-=-==-⎰⎰⎰⎰ 2()2a b x +-Q 在(a,b)区间上不变号，且连续可积，由第二积分中值定理 ''''322''()()()()()()()(,)222224b b aa f ab f a b b a R f x dx x dx f a b ηξξξ++-=-=-=∈⎰⎰综上所述3''()()()()()()224baa b b a f x dx Z R f b a f f ξ+-=+=-+⎰证毕6、计算积分1x e dx ⎰，若分别用复化梯形公式和复化Simpson 公式，问应将积分区间至少剖分多少等分才能保证有六位有效数字？解：①由复化梯形公式的误差限32''522()1()()101212122T b a b a e R f h f e n n η---=-≤=≤⨯可解得：212.85n ≥即至少剖分213等分。

第一章答案
第二章答案
第三章答案
0 0.5 0.5 1 1 2.5000
5.0000 5.5000
第四章答案
2 10.5000 19.0000 19.5000
3 42.5000 91.0000 91.5000
4 170.5000 315.0000 315.5000
5 682.5000 1467.0000 1467.5000
第八章答案
练习： 第一章
答案
练习二 A 的哪个特征向量？ 若 A 的按模最大的特征值是单根，用幂法求此特征 值的收敛速度由什么量来决定？怎样改进幂法的收敛速度？
2、 反幂法收敛到矩阵的哪个特征向量？ 在幂法或者反幂法中，为什么每步都要将迭代向量规范化？
1.32
1.68
2.08
2.52
3.00
解答下列问题 （1）试列出相应的差分表； （2）写出牛顿向前插值公式； （3）用二次牛顿前插公式计算 f（0.225）；
例3已知当 x=-1，0，2，3时，对应的函数值为
，
，
，
，
，求 的四次 Newton 插值多项式。
例4 设 对 n=1,2,3时
，证明：
例5 设 （1）
第一章答案第二章答案第三章答案第四章答案050525000500005500010500019000019500021000000000000000380001950004250009100009150001700000000000000018199999999999999166363636363636371705000315000031550001623809523809523716578947368421051161794871794871796825000146700001467500016058823529411764161208791208791201603825136612021827305000505100005051500016014662756598241160349206349206351601109350237717910922500023483000023483500016003663003663004160074982958418521600238500851788743690500080827000080827500016000915583226515160021777865769151600069286350589则开根号得400011444626607140002722140595534000086607000640对应的特征向量为第五章答案第六章答案2727930204331053600038939418364475947673代入数据得132解

代过程收敛 ．
１６畅 考虑下述修正的 New ton 公式（Steffenson 方法）
xn ＋ １ ＝ xn － f（ xn ）／ Dn ，
Dn ＝
f（ xn ＋
f（ xn ）） － f（ xn ）
f（ xn ） ，
n≥ ０ ．
假定 f′（ x 倡 ） ≠ ０ ，证明它对单根是一个二阶方法 ．
１７畅 设在区间［ a ，b］上函数 y ＝ f （ x）可逆 ，即存在足够光滑的 函数 Q（ y） ，使 x ＝ Q（ y） ．再设 x 倡 ∈ ［ a ，b］为 f （ x） ＝ ０ 的根 ，试把 x 倡 ＝ Q（０） ＝ Q（ y － y）用 T aylor 公式展开 ，从而导出下列公式 ：
x０ ＝ ０畅 ６ ，计算精确到 １０ － ５ ．
２０畅 试确定常数 p 、q 、r ，使迭代公式
xk ＋ １
＝
p xk
＋
q
a x２
k
＋
r
a２ x５
k
产生的序列｛ xk ｝收敛到３ a ，并使其收敛阶尽可能高 ．
２１畅 用弦截法 求 f （ x） ＝ x３ ＋ ２ x２ ＋ １０ x － ２０ ＝ ０ 的 根 ，要 求
x 倡 － xn
＝１＋ ｛２［ f″（ x倡 ）］２ （ x倡
－ ［ f″（ x倡 ）］２ （ x倡 － xn ）３
－
xn ）２
－
１ ２
［
f″（
x倡
）］２ （ x倡
－
xn ）２ ｝（ x倡
－
xn ）
＝１－
２ ３
＝
１ ３
．
可知 ，在重根附近 ，（６畅 ２０）式的敛速是线性的 ．

第六章作业3（1）,02223=--+x x x 解：二分法m 文件：function [c,err,k]=bisect(f,a,b,delta) ya=feval(f,a);yb=feval(f,b); if ya*yb>0error('f(a)*f(b)>0'); endflag=1;k=0; while flag==1k=k+1;c=(a+b)/2;yc=feval(f,c); if yc==0a=c;b=c; elseif yb*yc>0 b=c;yb=yc; elsea=c;ya=yc; endif abs(b-a)<=2*delta flag=0; end endc=(a+b)/2;err=abs(b-a)/2; return 主程序：f=inline('x^3+2*x^2-x-2'); a=0;b=3;delta=0.5e-2; [c,err,k]=bisect(f,a,b,delta) 输出结果：c = 0.9990；err = 0.0029；k =9试位法m 文件function [c,err,k]=shiweufa(f,a,b,delta) ya=feval(f,a);yb=feval(f,b); if ya*yb>0error('f(a)*f(b)>0'); endflag=1;k=0; while flag==1k=k+1;c=b-yb*(b-a)/(yb-ya);yc=feval(f,c);if yc==0a=c;b=c; elseif ya*yc<0b=c;yb=yc; elsea=c;ya=yc; endif abs(yc)<=delta flag=0; end endc=b-yb*(b-a)/(yb-ya);err=abs(yc); return试位法主程序：f=inline('x^3+2*x^2-x-2'); a=0;b=3;delta=0.5e-2;[c,err,k]=shiweufa(f,a,b,delta) 输出结果：c =0.9995；err =0.0045；k =234.用迭代方法求解下列方程在给定0x 附近的根，要求误差不超过10-4. （3）.3,05ln 202==--x x x 解：迭代法m 文件function [c,err,k]=iteration(f,x0,delta) ya=feval(f,x0); k=1;while abs(ya-x0)>1.0e-4 x0=ya;ya=feval(f,ya); k=k+1; end c=ya;err=abs(ya-x0); 主程序：f=inline('sqrt(2*logx+5)'); x0=3;delta=1.0e-4;[c,err,k]=iteration(f,x0,delta)输出结果：c =3.4495；err =5.6447e-005；k = 810.分别用（1）牛顿法，取;20-=x （2）弦截法，取;1.2,210-=-=x x （3）抛物线法，取,2.2,1.2,2210-=-=-=x x x 求方程0433=+-x x 在20-=x 附近的根，并比较各算法的数值表现。

第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。

解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。

解：直接根据定义和式(1.2.2)(1.2.3)则得有5位有效数字，其误差限，相对误差限有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？（1）（2）解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。

（1）（2）4.近似数x*=0.0310,是 3 位有数数字。

5.计算取，利用：式计算误差最小。

四个选项：第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。

线性插值时，用0.5及0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值误差限，故2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h 应取多少?解：用误差估计式（5.8），令因得3. 若，求和.解：由均差与导数关系于是4. 若互异，求的值，这里p≤n+1.解：，由均差对称性可知当有而当P＝n＋1时于是得5. 求证.解：解：只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式（5.17）得其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。

对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电通，力1根保过据护管生高线产中0不工资仅艺料可高试以中卷解资配决料置吊试技顶卷术层要是配求指置，机不对组规电在范气进高设行中备继资进电料行保试空护卷载高问与中题带资2负料2，荷试而下卷且高总可中体保资配障料置2试时32卷，3各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并25工且52作尽22下可护都能1关可地于以缩管正小路常故高工障中作高资；中料对资试于料卷继试连电卷接保破管护坏口进范处行围理整，高核或中对者资定对料值某试，些卷审异弯核常扁与高度校中固对资定图料盒纸试位，卷置编工.写况保复进护杂行层设自防备动腐与处跨装理接置，地高尤线中其弯资要曲料避半试免径卷错标调误高试高等方中，案资要，料求编试技5写、卷术重电保交要气护底设设装。备备置管4高调、动线中试电作敷资高气，设料中课并技3试资件且、术卷料中拒管试试调绝路包验卷试动敷含方技作设线案术，技槽以来术、及避管系免架统不等启必多动要项方高方案中式；资，对料为整试解套卷决启突高动然中过停语程机文中。电高因气中此课资，件料电中试力管卷高壁电中薄气资、设料接备试口进卷不行保严调护等试装问工置题作调，并试合且技理进术利行，用过要管关求线运电敷行力设高保技中护术资装。料置线试做缆卷到敷技准设术确原指灵则导活：。。在对对分于于线调差盒试动处过保，程护当中装不高置同中高电资中压料资回试料路卷试交技卷叉术调时问试，题技应，术采作是用为指金调发属试电隔人机板员一进，变行需压隔要器开在组处事在理前发；掌生同握内一图部线纸故槽资障内料时，、，强设需电备要回制进路造行须厂外同家部时出电切具源断高高习中中题资资电料料源试试，卷卷线试切缆验除敷报从设告而完与采毕相用，关高要技中进术资行资料检料试查，卷和并主检且要测了保处解护理现装。场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。

第六章非线性方程的数值解法习题解答填空题：1. 求方程()x f x =根的牛顿迭代格式是__________________。

Ans ：1()1()n n n n n x f x x x f x +-=-'-2.求解方程在(1， 2)内根的下列迭代法中，（1) (2）(3)（4）收敛的迭代法是（A ）.A ．(1)和（2)B 。

（2)和（3) C. （3）和(4) D. (4）和（1)3.若0)()(<b f a f ，则0)(=x f 在),(b a 内一定有根. ( )4．用二分法求方程01)(3=-+=x x x f 在区间［0，1］内的根,进行一步后根的所在区间为 ,进行两步后根的所在区间为 . （答案［0。

5,1]， ［0.5,0。

75］）计算题:1、已知方程3210x x --=在 1.5x =附近有根，将方程写成以下三种不同的等价形式： ①211x x =+;②321x x =+11x x =-试判断以上三种格式迭代函数的收敛性,并选出一种较好的格式。

解:①令121()1x x ϕ=+，则'132()x x ϕ=-，'132(1.5)0.592611.5ϕ=≈<,故迭代收敛； ②令322()1x x ϕ=+则2'2322()(1)3x x x ϕ-=+，'2(1.5)0.45581ϕ≈<,故迭代收敛； ③令31()1x x ϕ=-,则'33()2(1)x x ϕ=-'3(1.5) 1.41421ϕ≈>，故迭代发散。

以上三中以第二种迭代格式较好。

2、设方程()0f x =有根,且'0()m f x M <≤≤。

试证明由迭代格式1()k k k x x f x λ+=-(0,1,2,)k =产生的迭代序列{}0k k x ∞=对任意的初值0(,)x ∈-∞+∞，当20Mλ<<时，均收敛于方程的根。

数值分析参考答案1.4 习题解答或提示1、解：（1）>> a=[1 2 3 ;4 5 6 ]'a =1 42 53 6（2）>> b=[9;7;5;3;1]b =97531（3）>> c=b(2:4)c =753（4）>> d=b(4:-1:1)d =3579（5）>> e=sort(b)e =13579（6）>> f=[3:b']f =3 4 5 6 7 8 92、解：>> x=[7 4 3 ];y=[-1 -2 -3];（1）>> u=[y,x]u =-1 -2 -3 7 4 3 （2）>> u=[x,y]u =7 4 3 -1 -2 -33、解：sum=0;a=[4 -1 2 -8 4 5 -3 -1 6 -7]; for i=1 : length(a)if a(i)>0, sum=sum+a(i); endendsumsum =214、解：m=input('input an array:')input an array:[1 2 5;3 1 2;4 1 3]m =1 2 53 1 24 1 35、解：sum(m)ans =8 4 10>> max(m)ans =4 2 5>> min(m)ans =1 1 26、解：function y=fun_es(x)y=0.5.*exp(x./3)-x.^2.*sin(x);>> fun_es(3)ans =0.0891>> fun_es([1 2 3])ans =-0.1437 -2.6633 0.08917、提示：本题主要考查的是随机数生成函数rand的使用方法，以及选取种子数的方法之一：使用clock命令。

可以参照课本的例1.5来编写函数。

8、解：function y=fun_xa()x=input('input the value of x：');n=input('input the value of n：');y=1;for i=1:1:ny=y+x^i/factorial(i); end>> fun_xa()input the value of x ：1 input the value of n ：4ans =2.70832.4 习题解答1 解：E(lnx)=(ln ’E(x)=)(1x E x =xδ=Er(x) 2. 解 Er(x 2)＝)(22x Er x xx ⨯=4% 3. 解：123451.1021,0.031,385.6,56.430,7 1.0x x x x x *****=====⨯分别有5 位，2位，4位，5位，2位有效数字4 解 4*1105.0)(-⨯=x E3*2105.0)(-⨯=x E1*3105.0)(-⨯=x E3*4105.0)(-⨯=x E=++)(*4*2*1x x x E +)(*1x E +)(*2x E )(*4x E =0.00105))()((*4*2x E x E E =)()()(*42*4*2*4*2x E x x x x E -5. 解 V=334r π Er(v)=)(//x Er V x dx dV ⨯⨯=3Er(x)%1)(3≤x Er%33.0)(≤x Er6. 解 7830100-=Y Y)783()(100E Y E =＝0.00057．解 x 1，2＝24561122-±＝56783±21,2105.0)x (-⨯=E 2105.0)783(-⨯=E98.27783≈x 1，2＝83.98 或 28.02 8．略。

* (5) x5 = 23.496， x5 = 23.494
1 2
1 * x5 − x5 = 23.496 − 23.494 = 0.002 < × 10 − 2 2
x5 具有 4 位有效数字， x5 → 23.50 (不能写为 23.49) x6 = 96.1 × 10 5 = 0.961 × 10 7
e( A1 ) 10 −2 =1 er ( A1 ) = ≤ 0.01 A1
不能肯定所得结果具有一位有效数字。
2 ) A* = 0.01 ( 2.01 + 2.00 ) ，
A2 = 0.01 (1.42 + 1.41) = 0.01 2.83 = 0.00353356 Λ
e( A2 ) = e(0.01
(2)23.46―12.753; (4)1.473 / 0.064 。
问经过上述运算后，准确结果所在的最小区间分别是什么？ 解：(1) x1 =0.1062， x2 =0.947， x1 + x2 =1.0532
e( x1 ) ≤
1 1 × 10 − 4 ， e( x2 ) ≤ × 10 − 3 2 2 1 1 × 10 − 4 + × 10 − 3 2 2
习题 1
1． 以下各表示的近似数，问具有几位有效数字？并将它舍入成有效数。
（1） x1 ＝451.023， （2） x 2 ＝－0.045 113， （3） x3 ＝23.421 3， （4） x 4 ＝ ， （5） x5 ＝23.496， （6） x6 ＝96×10 5 ， （7） x7 ＝0.000 96， （8） x8 ＝－8 700， 解：(1) x1 = 451.023
e( A2 ) 0.12486 × 10 −4 = 0.3533547 × 10 − 2 er ( A2 ) = ≤ 0.00353356 A2

数值分析课后习题及答案第一章绪论（12）第二章插值法（40-42）2、当时，，求的二次插值多项式。

[解]。

3、给出的数值表用线性插值及二次插值计算的近似值。

X 0.4 0.5 0.6 0.7 0.8 -0.916291 -0.693147 -0.510826 -0.357765 -0.223144 [解]若取，，则，，则，从而。

若取，，，则，，，则，从而补充题：1、令，，写出的一次插值多项式，并估计插值余项。

[解]由，可知，，余项为，故。

2、设，试利用拉格朗日插值余项定理写出以为插值节点的三次插值多项式。

[解]由插值余项定理，有，从而。

5、给定数据表：，1 2 4 6 7 4 1 0 1 1 求4次牛顿插值多项式，并写出插值余项。

[解]一阶差商二阶差商三阶差商四阶差商 1 42 1 -34 0 6 17 1 0 由差商表可得4次牛顿插值多项式为：，插值余项为。

第三章函数逼近与计算（80-82）26、用最小二乘法求一个形如的经验公式，使它与下列数据相拟合，并求均方误差。

19 25 31 38 44 19.0 32.3 49.0 73.3 97.8[解]由。

又，，，故法方程为，解得。

均方误差为。

27、观测物体的直线运动，得出以下数据：时间t（秒）0 0.9 1.9 3.0 3.9 5.0 距离s（米）0 10 30 5080 110 [解]设直线运动为二次多项式，则由。

，。

又，，，故法方程为，解得。

故直线运动为。

补充题：1、现测得通过某电阻R的电流I及其两端的电压U如下表：I ……U ……试用最小二乘原理确定电阻R的大小。

[解]电流、电阻与电压之间满足如下关系：。

应用最小二乘原理，求R使得达到最小。

对求导得到：。

令，得到电阻R为。

2、对于某个长度测量了n次，得到n个近似值，通常取平均值作为所求长度，请说明理由。

[解]令，求x使得达到最小。

对求导得到：，令，得到，这说明取平均值在最小二乘意义下误差达到最小。

数值分析智慧树知到课后章节答案2023年下长安大学长安大学第一章测试1.算法效率就是指算法的快慢（）答案:错2.数值分析的任务就是：根据要求解的数学问题去设计算法（）答案:错3.用3.14近似π的有效数字位数是( )答案:34.真值经‘四舍五入’得到的近似数一定是有效数 ( )答案:对5.自然底数e=2.718281828的近似数2.7，2.71，2.718，2.7182中，有效数有（）个答案:2第二章测试1.n+1个互异节点，能够构造多少个拉格朗日插值基函数？（）答案:n+12.插值条件越多，拉格朗日插值多项式和原函数之间的误差越小（）答案:错3.通过牛顿插值法构造插值多项式时，首先需要建立什么？（）答案:差商表4.相同插值条件下，牛顿插值多项式和拉格朗日插值多项式的次数是一样的（）答案:对5.埃尔米特插值相比于拉格朗日插值的区别在哪？（）答案:多了一些与导数相关的插值条件6.分段插值主要解决了什么问题？（）答案:插值次数过高7.分段插值的本质就是在多个区间上做了多次的拉格朗日插值（）答案:对第三章测试1.在C[a, b]中，是范数（）答案:对2.在C[a,b]中，内积诱导范数与函数的2-范数相等（）答案:对3.可以利用相邻三项的关系确定一个正交多项式系，且结果唯一（）答案:错4.连续函数最佳平方逼近法中涉及的范数是连续函数空间中哪种范数（）答案:2-范数5.内积空间中线性无关元素确定的Gram矩阵是实对称正定矩阵（）答案:对6.函数，则 ( )答案:7.连续函数最佳平方逼近法中，平方误差一定是一个（）答案:非负数8.在对一组离散数据进行函数近似时，可以选用的方法有（）答案:曲线拟合;Lagrange插值;Newton插值9.在离散数据最小二乘曲线拟合问题中，所涉及的范数是（）答案:实向量空间2-范数10.线性矛盾方程组的最小二乘解是存在且唯一的 ( )答案:错11.下面是Newton-Cotes公式中Cotes系数特点的是（）答案:全为正（时）;和为1;对称性第四章测试1.含有n+1个互异求积节点，代数精确度至少为n的数值求积公式是（）答案:存在且唯一2.具有n次代数精确度的数值求积公式是插值型求积公式 ( )答案:对3.含有n+1个求积节点的插值型求积公式至少具有n次代数精确度 ( )答案:对4.Simpson公式的代数精确度为（）答案:3次5.积分区间为[a, b]，Simpson公式的Cotes系数为（）答案:1/6 4/6 1/66.用相同的求积节点对同一定积分进行近似求解，通常复化Simpson比复化梯形公式更准确 ( )答案:对7.被积函数足够光滑，复化Simpson公式的收敛阶数是 ( )答案:错8.被积函数足够光滑，复化梯形公式的收敛阶数是 ( )答案:对9.被积函数足够光滑，复化Simpson公式的收敛阶数是（）答案:第五章测试1.为什么在消元前要选择主元？（）答案:增强算法稳定性2.当线性方程组Ax=b的系数矩阵A是（）时，可用回代法求解.答案:上三角矩阵3.用高斯顺序消去法解线性方程组时，消元能进行到底的充分必要条件是（）.答案:系数矩阵A的前n-1阶顺序主子式非零4.矩阵A的哪种分解对应着高斯顺序消去法？（）答案:Doolittle5.n阶三对角矩阵A能够进行三角分解的充要条件为（）.答案:A的前n-1阶顺序主子式都非零第六章测试1.向量序列还是矩阵序列，也不管是定义中的按范数收敛还是按分量收敛，不可转化为数列的收敛。