2010-2011学年《数值分析》答案

合集下载

数值分析课后部分习题答案

习题一（P.14）1. 下列各近似值均有4个有效数字,300.2,521.13,001428.0***===z y x ,试指出它们的绝对误差和相对误差限.解*20.001428=0.142810x -=⨯有4个有效数，即4n =，2m =-由有效数字与绝对误差的关系得绝对误差限为611101022m n --⨯=⨯，由有效数字与相对误差的关系得相对误差限为(1)3111101022n a ---⨯=⨯； *213.521=0.1352110y =⨯有4个有效数，即4n =，2m =由有效数字与绝对误差的关系得绝对误差限为211101022m n --⨯=⨯，由有效数字与相对误差的关系得相对误差限为(1)3111101022n a ---⨯=⨯； *12.300=0.230010z =⨯有4个有效数，即4n =，1m =由有效数字与绝对误差的关系得绝对误差限为311101022m n --⨯=⨯，由有效数字与相对误差的关系得相对误差限为(1)3111101024n a ---⨯=⨯.2．下列各近似值的绝对误差限都是31021-⨯,试指出它们各有几位有效数字.***2.00021,0.032,0.00052x y z ===解*12.000210.20002110x ==⨯，即1m =由有效数字与绝对误差的关系得 311101022m n --⨯=⨯，即3m n -=-，所以，2n =；*10.0320.3210y ==⨯，即1m =由有效数字与绝对误差的关系得 311101022m n --⨯=⨯，即3m n -=-，所以，4n =；*30.000520.5210z -==⨯，即3m =-由有效数字与绝对误差的关系得 311101022m n --⨯=⨯，即3m n -=-，所以，0n =.4.设有近似数35.2,84.1,41.2***===z y x 且都有3位有效数字，试计算***z y x S +=，问S 有几位有效数字.解方法一因*1*1*12.41=0.24110, 1.840.18410, 2.350.23510x y z =⨯==⨯==⨯都有3位有效数字，即3n =，1m =，则211|(*)|101022m n e x --≤⨯=⨯，211|(*)|101022m n e y --≤⨯=⨯，211|(*)|101022m n e z --≤⨯=⨯，|(**)||*(*)*(*)|*|(*)|*|(*)|e y z z e y y e z z e y y e z ≈+≤+222112.3510 1.8410 2.0951022---≤⨯⨯+⨯⨯=⨯，221|(***)||(*)(**)|10 2.095102e x y z e x e y z --+≈+≤⨯+⨯1110.259510102--=⨯≤⨯，又1***=2.41 1.84 2.350.673410x y z ++⨯=⨯，此时1m =，1m n -=-，从而得2n =.方法一因*1*1*12.41=0.24110, 1.840.18410, 2.350.23510x y z =⨯==⨯==⨯都有3位有效数字，即3n =，1m =，则211|(*)|101022m n e x --≤⨯=⨯，2110(*)2|(*)|=||* 2.41r e x e x x -⨯≤， 211|(*)|101022m n e y --≤⨯=⨯，2110(*)2|(*)|=||* 1.84r e y e y y -⨯≤，211|(*)|101022m n e z --≤⨯=⨯，2110(*)2|(*)|=||* 2.35r e z e z z -⨯≤|(**)||(*)(*)|r r r e y z e y e z ≈+，***|(***)||(*)(**)|******r r rx y z e x y z e x e y z x y z x y z +≈+++2.41 1.84 2.35|(*)||(*)+(*)|2.41 1.84 2.35 2.41 1.84 2.35r rr e x e y e z ⨯≤++⨯+⨯22211110 1.8410 2.35102222.41 1.84 2.35 2.41 1.84 2.35 2.41 1.84 2.35---⨯⨯⨯⨯⨯≤+++⨯+⨯+⨯20.385410-<⨯21102-<⨯，由有效数字与绝对误差的关系得2n =.5.序列{}n y 有递推公式),2,1(,1101 =-=-n y y n n若41.120≈=y （三位有效数字），问计算10y 的误差有多大，这个计算公式稳定吗？解用0ε表示0y 的误差，由41.120≈=y ，得0=0.0042ε，由递推公式),2,1(,1101 =-=-n y y n n ，知计算10y 的误差为810=0.4210ε⨯，因为初始误差在计算的过程中被逐渐的放大，这个计算公式不稳定.习题2 ( P.84)3.证明()1nkk lx ==∑，对所有的x其中()k l x 为Lagrange 插值奇函数. 证明令()1f x =，则()1i f x =，从而 0()()()()nnn k k k k k L x l x f x l x ====∑∑，又(1)1()()()0(1)!n n n f R x x n ξω++==+，可得 ()()1n l x f x ==，从而()1nkk lx ==∑.4. 求出在=012x ，，和3处函数2()1f x x =+的插值多项式. 解方法一因为给出的节点个数为4，而2()1f x x =+从而余项(4)34()()()04!f R x x ξω==，于是233()()()()=+1L x f x R x f x x =-=（n 次插值多项式对次数小于或等于的多项式精确成立）.方法二因为(0)1(1)2(2)5(3)10f f f f ====，，，，而0(1)(2)(3)1()=-(1)(2)(3)(01)(02)(03)6x x x l x x x x ---=------，1(2)(3)1()=(2)(3)(10)(12)(13)2x x x l x x x x --=-----，2(1)(3)1()=-(1)(3)(20)(21)(23)2x x x l x x x x --=-----，3(1)(2)1()=(1)(2)(30)(31)(32)6x x x l x x x x --=-----，从而30123()()(0)()(1)()(2)()(3)L x l x f l x f l x f l x f =+++2=+1x .5. 设2()[,]f x C a b ∈且()()0f a f b ==，求证21max |()|()max |()|8a x ba xb f x b a f x ≤≤≤≤''≤-.证明因()()0f a f b ==，则1()0L x =，从而1()()()()()2!f f x R x x a x b ξ''==--，由极值知识得 21max |()|()max |()|8a x ba xb f x b a f x ≤≤≤≤''≤-6. 证明 (()())()()()(+)f x g x f x g x f x g x h ∆=⋅∆+∆⋅. 证明由差分的定义(()())(+)()()()f xg x f xh g x h f x g x ∆=+-[(+)()()(+)][()()()()]f x h g x h f x g x h f x g x h f x g x =+-++-()()()(+)f x g x f x g x h =⋅∆+∆⋅或着 (()())(+)()()()f x g x f x h g x h f x g x ∆=+-[(+)()()()][()()()()]f x hg xh f x h g x f x h g x f x g x =+-+++- ()()()()f x h g x f x g x =+⋅∆+∆⋅7. 证明 n 阶差商有下列性质(a ) 如果()()F x cf x =，则0101[,,,][,,,]n n F x x x cf x x x =. (b ) 如果()()()F x f x g x =+，则010101[,,,][,,,][,,,]n n n F x x x f x x x g x x x =+.证明由差商的定义 (a ) 如果()()F x cf x =，则12011010[,,,]-[,,,][,,,]n n n n F x x x F x x x F x x x x x -=-120110[,,,]-[,,,]n n n cf x x x cf x x x x x -=-120110[,,,]-[,,,]n n n f x x x f x x x c x x -=⋅-01[,,,]n cf x x x =.(b ) 如果()()()F x f x g x =+，则12011010[,,,]-[,,,][,,,]n n n n F x x x F x x x F x x x x x -=-12120110110[[,,,][,,,]]-[[,,,][,,,]]n n n n n f x x x g x x x f x x x g x x x x x --++=-12011120110,,,]-[,,,][,,,][,,,]+n n n n n n f x x x f x x x g x x x g x x x x x x x ---=--[ 0101[,,,][,,,]n n f x x x g x x x =+8. 设74()3431f x x x x =+++，求0172,2,,2]f [，0182,2,,2]f [.解由P.35定理7的结论(2)，得7阶差商0172,2,,2]=3f [(()f x 的最高次方项的系数)，8阶差商0182,2,,2]=0f [(8阶以上的差商均等与0).9. 求一个次数不超过4次的多项式()P x ，使它满足：(0)(0)0P P '==，(1)(1)1P P '==，(2)1P =.解方法一先求满足插值条件(0)0P =，(1)=1P ，(2)1P =的二次插值多项式2()P x 213=22x -+(L-插值基函数或待定系数法)，设()P x 22=()(1)(2)(1)(2)P x Ax x x Bx x x +--+--213=22x x -+2+(1)(2)(1)(2)Ax x x Bx x x --+-- 从而()P x '323=4B +(39)(641)(2)2x A B x A B x A -+-+-++，再由插值条件(0)0P '=，(1)1P '=，得3=,4A -1=,4B所以 ()P x 213=22x x -+231(1)(2)(1)(2)44x x x x x x ---+--，即 ()P x 41=4x 332x -29+4x .方法二设()P x 23401234=a a x a x a x a x ++++，则 ()P x '231234=234a a x a x a x +++由插值条件(0)(0)0P P '==，(1)(1)1P P '==，(2)1P =，得010********0123400++++1+2+3+41+2+4+8+161a a a a a a a a a a a a a a a a =⎧⎪=⎪⎪=⎨⎪=⎪=⎪⎩ 解得 234931=,=-,=424a a a ，从而()P x 41=4x 332x -29+4x . 方法三利用埃尔米特插值基函数方法构造. 10. 下述函数()S x 在[1,3]上是3次样条函数吗？3232321,12()=92217,23x x x x S x x x x x ⎧-++≤≤⎨-+-+≤≤⎩ 解因为22362,12()=31822,23x x x S x x x x ⎧-+≤≤'⎨-+-≤≤⎩， 66,12()=618,23x x S x x x -≤≤⎧''⎨-+≤≤⎩而12(2)=1=(2)S S ，12(2)=2=(2)S S ''，12(2)=6=(2)S S ''''，又()S x 是三次函数，所以函数()S x 在[1,3]上是3次样条函数.补设f (x )=x 4，试利用L-余项定理写出以-1,0,1,2为插值节点的三次插值多项式.解因为 (4)34()()()(+1)(1)(2)4!f R x x x x x x ξω==--，从而3233()()()22L x f x R x x x x =-=+-习题3 ( P.159)1．设n k k x 0)}({=ϕ为],[b a 上具有权函数0)(≥x ω的正交多项式组且)(x k ϕ为首项系数为1的k 次的多项式，则n k k x 0)}({=ϕ于],[b a 线性无关.解方法一因为n k k x 0)}({=ϕ为],[b a 上具有权函数0)(≥x ω的正交多项式组,则其Gram 行列式不等于零，采用反证法：若{}n ϕϕϕ,,,10 于],[b a 线性相关，于是，存在不全为零,,,,10n c c c 使0011()()()0,[,]n n c x c x c x x a b ϕϕϕ+++=∈上式两边与i ϕ作内积得到0011(,)(,)(,)0(0,1,,)i i n i n c c c i n ϕϕϕϕϕϕ+++==，由于{}i c 不全为零，说明以上的齐次方程组有非零解),,,,(10n c c c 故系数矩阵的行列式为零，即{}0,,,10=n G ϕϕϕ 与假设矛盾.方法二因为n k k x 0)}({=ϕ为],[b a 上具有权函数0)(≥x ω的正交多项式组,则其Gram 行列式不等于零，由( P.95)定理2得n k k x 0)}({=ϕ于],[b a 线性无关.2．选择α，使下述积分取得最小值1221()[],a x x dx α--⎰120()()x b e x dx α-⎰解1221()[]a x x dx αα-∂-∂⎰1221=[]x x dx αα-∂-∂⎰1221=2[]()x x x dx α--⋅-⎰5112=5x α-4=5α，令1221[]=0x x dx αα-∂-∂⎰，得=0α. 12()()x b e x dx αα∂-∂⎰120=()xe x dx αα∂-∂⎰1=2()()x e x x dx α-⋅-⎰2=23α- 令120()=0x e x dx αα∂-∂⎰，得=3α.3．设],3,1[,1)(∈=x xx f 试用},1{1x H 求)(x f 一次最佳平方逼近多项式.解取权函数为()x x ω=(为了计算简便)，则32311(1,1)42x xdx ===⎰,33321126(1,)(,1)33x x x x dx ====⎰, 343311(,)204x x x x dx ===⎰,33111((),1)2f x xdx x x=⋅==⎰，3232111((),)42x f x x x dx x =⋅==⎰，得法方程0126423264203a a ⎡⎤⎢⎥⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎢⎥⎣⎦，解得011211311a a ⎧=⎪⎪⎨⎪=-⎪⎩，所以)(x f 的一次最佳平方逼近多项式1123()1111P x x =-. 8．什么常数C 能使得以下表达式最小？ ∑=-ni x i iCe x f 12))((解21(())i n x i i f x Ce C =∂-∂∑1=2(())()i i nx x i i f x Ce e =-⋅-∑，令21(())=0i nx i i f x Ce C =∂-∂∑，得121()(),iinx x ii nx xx i f x ef x e C e e e=-=⋅==∑∑（）（，）. 14．用最小二乘法求解矛盾方程组2+314921x y x y x y =⎧⎪-=-⎨⎪-=-⎩. 解方法一方程组可变形为31+22491122x y x y x y ⎧=⎪⎪-=-⎨⎪⎪-=-⎩，原问题转化成在已知三组离散数据3142211()922t f t ----下求一次最小二乘逼近函数1()P x x yt =+(x 与y 为一次函数的系数，t 为自变量)，取1H 基{}1,t ,求解法方程331133321113()()i i i i i i i i i i i x t f x t t t f x y =====⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦∑∑∑∑∑, 即3-3-93737-32x y ⎡⎤⎡⎤⎡⎤⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦，得到矛盾方程组的解为37=-3156=31x y ⎧⎪⎪⎨⎪⎪⎩. 方法二方程组可变形为31+22491122x y x y x y ⎧=⎪⎪-=-⎨⎪⎪-=-⎩，令(,)I x y 2223111=+-+4+9++2222x y x y x y --（）（）（）(,)I x y x ∂∂3111=2+-+24+9+2+2222x y x y x y ⨯⨯-⨯-（）（）（）=6618x y -+，(,)I x y y ∂∂331111=+44+9+222222x y x y x y ⨯--⨯--⨯-（）（）（） 37=3372x y -+- 令(,)0(,)0I x y x I x y y∂⎧=⎪∂⎪⎨∂⎪=⎪∂⎩，得3373372x y x y -=-⎧⎪⎨-+-⎪⎩，解之得矛盾方程组的解为37315631x y ⎧=-⎪⎪⎨⎪=⎪⎩. 习题47. 对列表函数 124810()0152127x f x求(5)(5).f f '''，解一阶微商用两点公式(中点公式),得(8)(2)10(5),63f f f -'≈= 二阶微商用三点公式(中点公式),首先用插值法求(5)f , 由(4)5,(8)21,f f ==得一次插值函数1()411,L x x =-从而 1(5)(5)9f L ≈=,于是,2(2)2(5)(8)4(5).39f f f f -+''≈= 8. 导出数值数分公式)]23()2(3)2(3)23([1)(3)3(h x f h x f h x f h x f h x f ---++-+≈并给出余项级数展开的主部.解由二阶微商的三点公式(中点公式),得213()[()2()()]2222h h h f x f x f x f x h h ''-≈+--+-,213()[()2()()]2222h h h hf x f x f x f x h ''+≈+-++-从而 (3)()()22()h h f x f x f x h''''+--≈3133=[()3()3()()]2222h h f x h f x f x f x h h +-++--- 将33()()()()2222h h f x h f x f x f x h ++--，，，分别在x 处展开,得2(3)3(4)4(5)55331313()=()()()()()()222!23!21313()()()()+()(1)4!25!2f x h f x f x h f x h f x h f x h f x h O h '''++⋅+⋅+⋅+⋅+⋅2(3)3(4)4(5)5511()=()()()()()()222!23!211()()()()()(2)4!25!2h h h h f x f x f x f x f x h h f x f x O h '''++⋅+⋅+⋅+⋅+⋅+2(3)3(4)4(5)5511()=()()()()()()()222!23!211()()()()()(3)4!25!2h h h h f x f x f x f x f x h h f x f x O h '''-+⋅-+⋅-+⋅-+⋅-+⋅-+2(3)3(4)4(5)55331313()=()()()()()()()222!23!21313()()()()()(4)4!25!2f x h f x f x h f x h f x h f x h f x h O h '''-+⋅-+⋅-+⋅-+⋅-+⋅-+(1)－(2)×3 +(3)×3－(4), 得(5)222131()[()2()()]()()22228h h h f x f x f x f x h f x h O h h ''--+--+-=-+,即余项主部为(5)21()8f x h -习题 5 (P. 299)3. 设n n R A ⨯∈为对称矩阵，且011≠a ，经高斯消去法一步后，A约化为11120T a a A ⎡⎤⎢⎥⎣⎦，试证明2A 亦是对称矩阵. 证明设1111()=T ij aa A a A α⎛⎫= ⎪⎝⎭，其中 21311=n a a a α⎛⎫⎪ ⎪ ⎪ ⎪⎝⎭，121311=n a a a a ⎛⎫ ⎪ ⎪ ⎪ ⎪⎝⎭，22232123=n n n nn a a a a a a a ⎛⎫⎪⎪ ⎪⎝⎭，则经高斯消去法一步后，A 约化为111111110TT a a A a a α⎡⎤⎢⎥⎢⎥-⎢⎥⎣⎦，因而211111T A A a a α=-，若n n R A ⨯∈为对称矩阵，则1A 为对称矩阵，且1=a α，易知211111T A A a a α=-为对称矩阵. 13. 设⎥⎦⎤⎢⎣⎡=989999100A(1) 计算2||||,||||A A ∞； (2) 计算∞)(A Cond ，及2)(A Cond . 解 (1) 计算||||=199A ∞,⎥⎦⎤⎢⎣⎡=989999100A，其特征值为1,299λ=，又⎥⎦⎤⎢⎣⎡=989999100A 为对称矩阵，则2=T A A A 的特征值为221,2(99λ=±，因此2||||99A ===+；(2)1989999100A --⎡⎤=-⎢⎥-⎣⎦，1||||=199A -∞，所以1()=||||||||=9801Cond A A A -∞∞∞⋅，1989999100A --⎡⎤=-⎢⎥-⎣⎦为对称矩阵，其特征值为1,299λ=-± 则1112()=()T A A A ---的特征值为221,2(99λ=，因此12||||99A -===+所以1222()=||||||||Cond A A A -⋅2(99=+15. 设,n n n A R x R ⨯∈∈，求证（1）1xx n x ∞∞≤≤； (2)∞∞≤≤An A A n11.证明 (2) 由（1）1x x n x∞∞≤≤，得1AxAx n Ax∞∞≤≤，则 11Ax Ax n Ax n x xx∞∞∞∞≤≤，从而11max max max nnnx Rx Rx RAxAx n Ax n xxx∞∞∀∈∀∈∀∈∞∞≤≤，由算子范数的定义max nx RAx Ax∞∞∀∈∞=，111max nx RAx A x∀∈=，得∞∞≤≤An A A n11.17. 设n n R W ⨯∈为非奇异阵，又设x为n R 上一向量范数，定义WxWx=，求证：Wx是nR 上向量的一种范数（称为向量的W 一范数）.证明 ①正定性，因Wx为一向量，0WxWx =≥，下证=0=0Wxx ⇔，⇒“”若=0Wx 即=0Wx ，由向量范数的正定性得=0Wx ，n n R W ⨯∈为非奇异阵，所以=0x ；⇐“”若=0x ，则=0Wx ，由向量范数的正定性得=0Wx 即=0Wx.②齐次性，任意实数α有=Wx W x Wxααα=，由向量范数的齐次性，得=WWxW x Wx Wx xααααα===；③ 三角不等式，任意实数,n n x R y R ∈∈，有+(+)=+Wx yW x y Wx Wy=，再由向量范数的三角不等式，得+(+)=+WWWx yW x y Wx Wy Wx Wy xy=≤+=+.习题 6 (P.347)1. 设有方程组(b )1231231232211221x x x x x x x x x +-=⎧⎪++=⎨⎪++=⎩，考查用Jacobi迭代法,G-S 迭代法解此方程组的收敛性.解系数矩阵分裂如下,122111221A -⎛⎫⎪= ⎪ ⎪⎝⎭D L U =--10022110112200-⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪=---- ⎪ ⎪⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭ Jacobi迭代矩阵为1()J D L U -=+=02211220-⎛⎫ ⎪-- ⎪ ⎪--⎝⎭， J 的特征方程为2211022λλλ-=，展开得 30λ=，即01λ=<，所以用Jacobi 迭代法解此方程组是收敛的.G-S 迭代矩阵为1()G D L U -=-11022=11012210--⎛⎫⎛⎫⎪ ⎪⋅- ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭100022=110010210-⎛⎫⎛⎫ ⎪ ⎪-⋅- ⎪⎪ ⎪ ⎪-⎝⎭⎝⎭122=023002-⎛⎫ ⎪- ⎪ ⎪⎝⎭， G 的特征方程为12221002λλλ---=-，展开得 (1)(2)(2)0λλλ---=，即1λ=或2λ=，由迭代基本定理得用G-S 迭代法解此方程组是不收敛的.4. 设有方程组Ax b =，其中A 为对称正定阵，且有迭代公式(1)()()()k k k x x b Ax ω+=+- (0,1,k =)，试证明当20ωβ<<时，上述迭代法收敛(其中A 的特征值满足0()A αλβ<≤≤).证明 A 为对称正定阵， A 的特征值满足0()A αλβ<≤≤，且20ωβ<<，则0()2A ωλ<<又迭代公式可变形为(1)()()k k x I A x bωω+=-+ (0,1,k =)，从而迭代矩阵 B I A ω=-，迭代矩阵的特征值为1()A ωλ-，且满足11()1A ωλ-<-<，即 |()|1B λ<，由迭代基本定理得该迭代法是收敛的.5. 设111a a A aa a a⎛⎫ ⎪= ⎪ ⎪⎝⎭，其中a 为实数，试确定a 满足什么条件时，解Ax b =的Jacobi 迭代法收敛.解系数矩阵分裂如下,111a a A aa a a ⎛⎫ ⎪= ⎪ ⎪⎝⎭D L U =--1001100a a aa aa--⎛⎫⎛⎫⎛⎫⎪ ⎪ ⎪=---- ⎪ ⎪⎪ ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭Jacobi迭代矩阵为1()J D L U -=+=000a a aa a a--⎛⎫⎪-- ⎪ ⎪--⎝⎭，J 的特征方程为0a aa a aaλλλ=，展开得 323320a a λλ--=，即a λ=-或2a λ=-，()max{||,|2|}J a a ρ=--()1J ρ<当且仅当1122a -<<，所以当1122a -<<时，解Ax b=的Jacobi 迭代法收敛.。

福州大学2010-2011年数值分析考题及答案1

得分评卷人
1、若向量 x (4, 2,3) ，则
T
x 2 =___ 29 _________
＝____ 6 ____，A 的
2、
1 1 A ，则 A 的谱半径 -5 1
=____6____
3、确定求积公式尽量高,则 A0=_

1
1
f ( x)dx A0 f (1) A1 f (0) A2 f '(1) 中的待定参数，使其代数精度
0 2 0 5、设 B 2 1 2 ,试用平面旋转矩阵对矩阵 A 进行 QR 分解，其中 Q 为正交 0 2 1
矩阵，R 为上三角阵（8 分）
4
记A1 A, 先将A的第一列变得与e1平行 cos = 0 2 0,sin = 1 04 04 0 1 0 0 1 0 0 P A 2 P A1 1 12 12 0 0 0 1
3、
h 用二步法 yn1 yn [ f ( xn , yn ) f ( xn1 , yn1 )] 求解一阶常微分方程初值问题 2
y f ( x, y ) 问：如何选择参数 , 的值，才使该方法的阶数尽可能地高？写出 y ( x0 ) y0
此时的局部截断误差主项，并说明该方法是几阶的。证明:局部截断误差为：
( x x )l ( x) 等于
i 0 i i
4
（ a ） 1 (c) 2 (d) 4
(a)
0
(b)
3、设 f ( x) 3x5 4 x 4 x 2 1 和节点 xk k / 2, k 0,1 则差商 f [ x0 , x1 x5 ] (a) 4 (b) 2 (c) 3 (d) 1 （（ c ） c ）

《数值分析》第一章答案

习题11．以下各表示的近似数，问具有几位有效数字？并将它舍入成有效数。

（1）*1x ＝451.023， 1x ＝451.01；（2）*2x ＝－0.045 113， 2x ＝－0.045 18；（3）*3x ＝23.421 3， 3x ＝23.460 4；（4）*4x ＝31，4x ＝0.333 3；（5）*5x ＝23.496， 5x ＝23.494；（6）*6x ＝96×510， 6x ＝96.1×510；（7）*7x ＝0.000 96， 7x ＝0.96×310-；（8）*8x ＝－8 700， 8x ＝－8 700.3。

解：(1) =*1x 451.023 =1x 451.01=-1*1x x 0.01311021-⨯≤，1x 具有4位有效数字。

→1x 451.0(2) -=*2x 0.045 113 -=2x 0.045 18=-<⨯-2*241021x x 0.045 18045113.0-=0.000 06731021-⨯<2x 具有2位有效数字，045.02-→x(3)=*3x 23.4213 =3x 23.4604=-3*3x x =-4604.234213.23=-4213.234604.23110210391.0-⨯≤3x 具有3位有效数字，4.233→x (不能写为23.5)(4) =*4x 31，=4x 0.3333=-4*4x x 41021000033.0-⨯< ,4x 具有4位有效数字，=4x 0.3333 (5) =*5x 23.496，=5x 23.494=-5*5x x =-494.23496.2321021002.0-⨯<5x 具有4位有效数字， →5x 23.50 (不能写为23.49)(6) =*6x 51096⨯71096.0⨯= =6x 5101.96⨯710961.0⨯==-6*6x x 710001.0-⨯72101021--⨯⨯≤6x 具有2位有效数字，57610961096.0⨯=⨯=x(7) =*7x 0.00096 371096.0-⨯=x 3*71096.0-⨯=x =-7*7x x 0 7x 精确 (8) 8700*8-=x 8x 3.8700-=8*8x x -010213.0⨯≤= 8x 具有4位有效数字，8x 8700-=精确 2.以下各数均为有效数字： (1) 0.1062 + 0.947; (3)2.747⨯6.83;(2)23.46―12.753; (4)1.473 / 0.064 。

《数值分析》参考答案

参考答案第1章一、选择题1. D2. C3. A4. B5. B二、填空题1. 函数题头 H1行帮助信息函数体注释部分函数题头2. nargin varargin3. A=rand(4)4. 单引号三、解答题1. for 语句和while 语句均可以实现循环执行的功能。

二者的区别在于，for 循环语句一般适用于已知道循环次数，而不知道循环运算的目标的问题，而while 循环语句则相反，一般适用于已知循环目标，而循环次数未知的问题。

2. 程序如下：function [highavg,weightavg]=avg_high_weight(varargin) n=length(varargin); highsum=0; weightsum=0; for i=1:n highsum=highsum+varargin{i}(1);weightsum=weightsum+varargin{i}(2);endhighavg=highsum/n; weightavg=weightsum/n;第2章一、选择题1. A2. B3. A4. C5. D二、填空题1. 1.7 1.73 1.7322. 3 13. 5％4. 3三、解答题1. 解：1*()()nn x nxx x ε-≈-1***()()n nr nxx x x x x nnxxε---≈=()0.02r ne x n ==2数值分析2. 解：*1 1.1021x =有五位有效数字；*20.031x =有两位有效数字；*3385.6x =有四位有效数字；*47 1.0x =⨯有一位有效数字。

3. 解：（1）*******124124()()()()x x x x x x εεεε++≤++433111101010222---=⨯+⨯+⨯3*1.0510ε-=⨯=（2）*********123231113()()x x x x x x x x x ε⋅⋅≈⋅-+⋅****221233()()x x x x x x -+⋅-*0.197ε≈=（3）******2242244**2441(/) |()()|()x x x x x x x xx ε≤---****2224**44|()()|r r x x x x xxεε=-***224*4||[|()||()|]r r x e x e x x≤+331110100.0312256.4800.03156.480--⎡⎤⨯⨯⎢⎥=+⎢⎥⎢⎥⎣⎦5*10ε-≤=4. 解：33**34433()43r R RV Rππεπ-=*2**2R R R R R RRR R-++=⋅*223R R R RR-≈⋅*3R R R-=⋅1%=故*1300R R R-=5. 解：设Y =*27.983Y =，*31102Y Y δ-=-≤⨯，028Y =，*028Y =，*0000Y Y δ=-=*111282827.983100Y Y ⎛⎛⎫-=---⨯ ⎪⎝⎝⎭1100δ≤，**22111127.983100100Y Y Y Y ⎛⎛⎫-=-⨯--⨯ ⎪⎝⎝⎭**111()()100Y Y Y Y =---112100100100δδδ≤+=仿此可得：*100n n n Y Y δ-≤则*31001001001101002Y Y δδ--≤==⨯即计算100Y 的误差界不超过31102-⨯参考答案 36. 解：解方程25610x x -+=得：28282x =±±由第5题知27.983具有五位有效数字，故可取：1282827.98355.983x =++=21280.0178655.983x =-≈=7. 解：设正方形的边长为x ，则其面积为2y x =。

数值分析部分答案

计算，解
Q f(x) ln(x Jx21),f(30)In(30 s/899)设u ^y899, y f (30)则u*
yu
u
1*
g u
0.0167
3
若改用等价公式
ln(x•.厂1)In (x1)
贝卩f(30)In(30x899)
此时
* *
yr u
u
1*
u
59.9833
7
第二章插值法
2
X
0.4
0.5
(y2*)10 (y「)
2
(y2*)10 (y°*)
S*)1010(yo*)
101011022
(x1)7
6* *
7y x
(x 1)
* *
y x
*2*
(32x)g x
6* *
*y g x
3 2x
* *
y x
(3 2.2)3计算y值，则
1
(3 2x )4
1*
7y x
(3 2x )7'
* *
y x
(3 2 <2)
(3)(x2/x4)
0.031 385.6
1.1021 385.6
x；
*ቤተ መጻሕፍቲ ባይዱ
(X4)
X4(X2)
* 2
X4
131
1056.43010
2 2
56.430 56.430
5
解：球体体积为V 4R
3
则何种函数的条件数为
2
Rgl R
1 V丨
43
-R3
3
3
r(V*) Cpgr(R*)3r(R*)
Cp
又Qr(V*)1

数值分析课后习题答案

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.求以下矩阵的特征值和特征向量。

《数值分析》第一章答案

《数值分析》第⼀章答案习题11．以下各表⽰的近似数，问具有⼏位有效数字？并将它舍⼊成有效数。

（1）*1x ＝451.023， 1x ＝451.01；（2）*2x ＝－0.045 113， 2x ＝－0.045 18；（3）*3x ＝23.421 3， 3x ＝23.460 4；（4）*4x ＝31， 4x ＝0.333 3；（5）*5x ＝23.496， 5x ＝23.494；（6）*6x ＝96×510， 6x ＝96.1×510；（7）*7x ＝0.000 96， 7x ＝0.96×310-；（8）*8x ＝－8 700， 8x ＝－8 700.3。

解：(1) =*1x 451.023 =1x 451.01=-1*1x x 0.01311021-?≤，1x 具有4位有效数字。

→1x 451.0(2) -=*2x 0.045 113 -=2x 0.045 18=-241021x x 0.045 18045113.0-=0.000 06731021-?<2x 具有2位有效数字，045.02-→x(3)=*3x x =-4604.234213.23=-4213.234604.231 10210391.0-?≤3x 具有3位有效数字，4.233→x (不能写为23.5) (4) =*4x 31 ，=4x 0.3333=-4*4x x 41021000033.0-?<,4x 具有4位有效数字，=4x 0.3333(5) =*5x 23.496，=5x 23.494=-5*5x x =-494.23496.2321021002.0-?<5x具有4位有效数字，→5x 23.50 (不能写为23.49)(6) =*6x 51096?710961.0?==-6*6x x 710001.0-?72101021--??≤6x 具有2位有效数字，57610961096.0?=?=x(7) =*7x 0.00096 371096.0-?=x3*71096.0-?=x =-7*7x x 0 7x 精确(8) 8700*8-=x 8x 3.8700-=8*8x x -010213.0?≤=8x 具有4位有效数字，8x 8700-=精确2.以下各数均为有效数字： (1) 0.1062 + 0.947; (3)2.747?6.83; (2)23.46―12.753; (4)1.473 / 0.064 。

中国石油大学《数值分析》2011年考试试题A卷及答案

f (4)(x)
1 2880
1 n
4
6
1 2
104
,
仅要 n 4 1 101 2.54 ,取 n 3 即对将[1,2] 作 6 等分,则有 240
(8 分)
2
1 ln xdx
1 [0 4(ln 7 ln 3 ln 11) 2(ln 4 ln 5) ln 2] 0.38628716327880 .
0.000040074
（ 4 分）
七、（10 分）(1)牛顿迭代格式
x(k 1)
x(k)
f f
(x(k ) ) '(x(k) )
x(k)
x(k) 1 (2
(x(k) )2 )(x(k) )1
1
(1 (2
)(
x( )(
)k ) 2 x(k ) )1
(2)
x(k 1)
lim
k
x(k)
1 1
fgdx
，取( x) ax bx3 ， f ( x) sin x ，则法方程为
(0 ,0 )
(1
,
0
)
(0 ,1) (1 , 1 )
a b
( (
f f
,0 ,1
) )
（ 4 分）
其中 0,0
1
x xdx
1
2， 3
0 ,1
(1 )(x(k) )2
lim
k
1
(2
)(x(k ) )1
c0
2
c 1
(5 分) (5 分)
1
x(k) 2
x(k) 3
1
x(k) 1
x(k) 3
/2
x3( k
1)

数值分析参考答案(第二章)doc资料

（2）若，则
证明：
（1）
得证。
+
得证。
14．求及。
解：
若
则
15．证明两点三次埃尔米特插值余项是
解：
若，且插值多项式满足条件
插值余项为
由插值条件可知
且
可写成
其中是关于的待定函数，
现把看成上的一个固定点，作函数
根据余项性质，有
由罗尔定理可知，存在和，使
即在上有四个互异零点。
根据罗尔定理，在的两个零点间至少有一个零点，
数值分析参考答案(第二章)
第二章插值法
1．当时， ,求的二次插值多项式。
解：
则二次拉格朗日插值多项式为
2．给出的数值表
X
0.4
0.5
0.6
0.7
0.8
lnx
-0.916291
-0.693147
-0.510826
-0.356675
-0.223144
用线性插值及二次插值计算的近似值。
解：由表格知，
若采用线性插值法计算即，
则
若采用二次插值法计算时，
3．给全的函数表，步长若函数表具有5位有效数字，研究用线性插值求近似值时的总误差界。
解：求解近似值时，误差可以分为两个部分，一方面，x是近似值，具有5位有效数字，在此后的计算过程中产生一定的误差传播；另一方面，利用插值法求函数的近似值时，采用的线性插值法插值余项不为0，也会有一定的误差。因此，总误差界的计算应综合以上两方面的因素。
解：函数的展式为
其中
又是次数为的多项式
为阶多项式
为阶多项式
依此过程递推，得是次多项式

《数值分析》所有参考答案

习题11 -以下各表示的近似数，问具有几位有效数字？并将它舍入成有效数(1)% = 451.023(2)x；=-0.045 113(3)x3 = 23.421 3,* 1(4)x4=3(5)x5 = 23.496,* /-(6)x6= 96x 10 ,(7)x；= 0.000 96,(8)x8 =-8 700, 解：(1) x；=451.023x1= 451.01;x2=—0.045 18;x3= 23.460 4;x4= 0.333 3;x5= 23.494;x6= 96.1 x 105;x；= 0.96X 10 'x8= —8 7003 x^ 451.01* 1 _1 一#x1—= 0.013兰一汇10 —, x1具有4 位有效数字。

％t451.02(2) x；二-0.045 113 x2二-0.045 18=0.045 1 8- 0.045113 =0.000 067 - 10 _32X2具有2位有效数字，x^ -0.045⑶x3 =23.4213 x3= 23.4604*X3— X3 = 23.4213 - 23.4604 = 23.4604 — 23 .4213 = 0.0391 X3具有3位有效数字，X3 > 23 .4 （不能写为23.5）* 1⑷ x4二,x4二0.3333 J 10_1 23二 23 .496 - 23.494 二 0.002X 6具有2位有效数字,75x 6 =0.9610= 96 102•以下各数均为有效数字:(1) 0.1062 + 0.947;(2)23.46— 12.753;(4) 1.473 / 0.064。

问经过上述运算后，准确结果所在的最小区间分别是什么? 解：(1) X i =0.1062， X 2 =0.947， X i +X 2 =1.05321e( )+ e(x 2 )兰 e( )+ e(x 2)兰一汉 10*X 4=0.000033：：： -10 一4 2,X 4具有4位有效数字，X 4 二0.3333(5) x 5 = 23.496, x 5 = 23.494X5具有4位有效数字,x 5 > 23.50 不能写为 23.49)(6)*57X 6 = 96100.96 10 57X 6=96.1 10 =0.96110*X6=0.001 10 _7< -10 ° 10 一7 2X 7 = 0.00096X 7 -0.9610° *X7-0.96 10’*X7=0X 7精确(8)二 -8700 x8二 -8 7 0.3*X8-X 81 = 0.3102X 8具有4位有效数字，X 8二-8700 精确e(xd| 兰丄。

数值分析课后习题答案

x2 6.6667x2 8.205
再解
1
15 56
x31.785,7得 x35.769
1 25069x4 0.47847x4 1.4872
1 x5 5.3718 x5 5.3718
2-10.证明下列不等式:
(1)x-yx-z+z-y; (2)|x-y|x-y;
证明 (1)x-y=(x-z)+(z-y)x-z+z-y
b.用Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=0.
102 x Байду номын сангаасy 1
100y 100
再用列主元Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=1.
x y
y 1
2
2-8.用追赶法求解方程组:
4 1
x1 100
1 4 1
x2 0
3-8.判定求解下列方程组的SOR方法的收敛性.
2 1 0 0 x1 1
1
0 0
2 1 0
1 2 1
0 12
x2 x3 x4
0 00
解直接可验证系数矩阵A是负定矩阵,所以-A是对称
1-3.为了使101/2的相对误差小于0.01%,试问应取几位有效数字?
解因为101/2=3.162…=0.3162…10,若具有n位有效数字,则其绝对误差限为0.5 101-n ,于是有
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
1 2
0
12 1,
1 2
1 2
0
12

2011数值分析试题及答案

122446
由于f(x)二si nx的4阶导数在[0,二]上的最大值为：M4=1,所以
5
误差为：|I-S2|：：——44=0.006641
2880x24
6.求解初值问题」y=sin(x+2y),0兰x兰2的改进Euler方法是否收敛？为什
.y(0) = 1
么？
解：由于|sin(x 2y)-sin(x 2y)|二| 2cos(x 2 )(y-y) 2 | y-y |
5.设f(x) = 4x33x-5，求差商f[0,1], f[1,2,3,4]和f[1,2,3,4,5]。
f(D…f(0)
解：f[0,1]==2-(-5) = 7
1-0
f [1,2,3,4^4，f[1,2,3,4,5]=0
3.解线性方程组丿X1-2忑=2的Jacobi迭代法是否收敛，为什么?
+9x2=3
即，函数f(x, y)二sin(x•2y)连续，且关于变量y满足Lipschitz条件，所以，改进Euler方法收敛。
所以，a=0, b=5/6，拟合曲线为：y=5/6x2
3.求满足条件f(0)=1,f(1)=2,f(2) =0,f(1)=0的三次插值多项式Ha(x)
的表达式。
解：设H3(x)二(^2)(ax2bx c)，则有:
1213
所以，H3(x) (x-2)(x2x 1) (x-3x-2)。
22
11
4.确定求积公式Jf(x)dx痒三f(-1)+Af(0)+A2f(1)中的待定系数，使其代数精度尽可能高，并问此公式是不是插值型求积公式.
解：令公式对f(x) = 1,x都精确成立，得：A,・A2= 3/2, A2= 1/2,
o
• • •

江理数值分析(计算方法)期末试卷3及参考答案

5 2
2xk3 5 3xk2 2
将 x0 1.5 代入上式，得
x1 1.34286 ， x2 1.37012 ， x3 1.32920 , x4 1.32827 , x5 1.32826
x5 x4 0.00001 104
所以，方程的近似根
x5 1.32826
5．解，Jacobi 迭代公式是：
(2) 求出用 Jacobi 迭代法和 Gauss-Seidel 迭代法求解该方程组的迭代
序列，同时证明该迭代法是收敛的。
矩阵.并判断用 Gauss-Seidel 迭代法求解该方程组的收敛性。（15 分）
第 3 页︵共 3 页︶
姓名
江西理工大学教务处
专业、班级
2007-2008-2 数值分析 A 参考答案
一．填空(每空 3 分，共 30 分)
1. 截断误差
2. x(x 2) ， x(x 1) ，
10
3. 1
2
4. xk1
xk
xk2 2xk
f (xk ) f (xk )
5. 6，5， 26 ，9
二．计算
1. 构造重节点的差商表：
n
x
y
一阶二阶三阶
0
1
0
1
2
-5
-5
2
3
-6
-1
2
3
4
3
9
5
72
3.解：改进的 Euler 公式是：
具体到本题中，求解的公式是：
代入求解得： y1 1.4 , y1 1.54
y2 2.276, y2 2.4832
yn1 yn hf (xn , yn )
yn1

(完整版)数值分析课后习题答案

第一章绪论习题一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),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。

数值分析课后题答案

数值分析第二章2．当1,1,2x =-时，()0,3,4f x =-,求()f x 的二次插值多项式。

解：0120121200102021101201220211,1,2,()0,()3,()4;()()1()(1)(2)()()2()()1()(1)(2)()()6()()1()(1)(1)()()3x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------==-+--则二次拉格朗日插值多项式为220()()k k k L x y l x ==∑0223()4()14(1)(2)(1)(1)23537623l x l x x x x x x x =-+=---+-+=+- 6．设,0,1,,j x j n =为互异节点，求证：（1）0()nkkj j j x l x x=≡∑ (0,1,,);k n =（2）0()()0nk jj j xx l x =-≡∑ (0,1,,);k n =证明（1）令()kf x x = 若插值节点为,0,1,,j x j n =，则函数()f x 的n 次插值多项式为0()()nk n j j j L x x l x ==∑。

插值余项为(1)1()()()()()(1)!n n n n f R x f x L x x n ξω++=-=+又,k n ≤(1)()0()0n n f R x ξ+∴=∴=0()nk kj j j x l x x =∴=∑ (0,1,,);k n =000(2)()()(())()()(())nk j j j n nj i k i k j j j i nnik ii kj j i j x x l x C x x l x C x x l x =-==-==-=-=-∑∑∑∑∑0i n ≤≤又由上题结论可知()nk ij jj x l x x ==∑()()0ni k i ik i k C x x x x -=∴=-=-=∑原式∴得证。

《数值分析》第二章答案

习题21. 分析下列方程各存在几个根，并找出每个根的含根区间：(1) 0cos =+x x ； (2) 0cos 3=-x x ； (3) 0sin =--x e x ； (4) 02=--x e x 。

解：(1) 0cos =+x x (A) x x x f cos )(+= ，0sin1)(≥-='x x f ,),(∞-∞∈x10cos 0)0(=+=f ,01cos 1)1cos(1)1(<+-=-+-=-f ∴ 方程(A) 有唯一根 ]0,1[*-∈x (2) 0cos 3=-x x (B) x x x f c o s 3)(-=,0sin 3)(>+='x x f , ),(∞-∞∈x 时010c o s03)0(<-=-⨯=f ，01cos 31cos 13)1(>-=-⨯=f ∴ 方程(B) 有唯一根 ]1,0[*∈x (3)sin =--xex (C)xex -=sinx x f sin )(1=, xex f -=)(2方程(C)有无穷个正根，无负根在[22,2πππ+k k ] 内有一根 )(1k x ，且0]2[lim )(1=-∞→πk x k k在[ππππ++k k 2,22]内有一根)(2k x ，且0])12([lim )(2=+-∞→πk x k k (示图如下) 3,2,1,0=k)(2x f x(4)02=--xex(D) xex-=2,)(21x x f = xex f -=)(2方程(D) 有唯一根 ]1,0[*∈x 当 0<x 时 (D)与方程2x ex -=- (E) 同解当 0<x 时 (E)无根 2. 给定方程 012=--x x ； (1)(2)若在[0 , 2]上用二分法求根，要使精确度达到6位有效数，需二分几次？解：012=--x x1) 01)(2=--=x x x f 1)1(-=f , 025.0)5.1(<-=f ,1)2(=f]2,5.1[*∈x, 618034.1251*=+=x)(5.1- 1.75(+) 2(+) )(5.1- 1.625(+) 1.75(+) )(5.1-1.5625(+) 1.625(+))(5625.1- )(59375.1-1.625(+)1102103125.02)5625.1625.1(-⨯<=-6.159375.1*≈≈x2位有效近似值为 1.6 2)00==a a ， 20==b b)(21k k k b a c +=kk k a b c x 2121*=-≤-+5102121-⨯≤k，51102≥-k60.162ln 10ln 51=≥-k∴ 只要2等分18次3. 为求0353=--x x 的正根，试构造3种简单迭代格式，判断它们是否收敛，且选择一种较快的迭代格式求出具有3位有效数的近似根。

数值分析参考答案

数值分析参考答案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．略。

2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准

2010-2011学年第一学期末数值分析考试试题A参考答案与评分标准D第2页，共5页第3页，共5页第4页，共5页()()()1121212,,n n n n n n h y y k k k f t y k f t h y hk +⎧=++⎪⎪=⎨⎪=++⎪⎩，它是 2 阶方法。

二、解答下列各题(每小题12分，共24分) 1. 用LU分解法求解线性方程组1234230156715112133721182x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪=⎪ ⎪ ⎪- ⎪ ⎪ ⎪--⎝⎭⎝⎭⎝⎭；解：() 2 3 0 1 -53 -2 1 2 41 2 1 -2 41 2 -1 2 1 -11 -2 2 -3015671511|2 1 133 7201182LU A b ⎛⎫-⎛⎫ ⎪⎪ ⎪-⎪ ⎪=−−−→ ⎪- ⎪ ⎪ ⎪--⎝⎭ ⎪⎝⎭分解1 0 0 023 0 13 1 0 0 0 -2 1 21 2 1 0 0 0 1 -21 2 -1 1 0 0 0 154,,41L U y -⎛⎫⎛⎫⎛ ⎪ ⎪⎪ ⎪ === ⎪ ⎪ ⎪ ⎪ -⎝⎭⎝⎭⎝ 1-2 2-1,x ⎫⎛⎫⎪ ⎪⎪ ⎪=⎪ ⎪⎪ ⎪⎭⎝⎭.. LU 矩阵(或对应元素每算对两个给1分) 2. x19 25 31 38 44 y 19.0 32.3 49.0 73.3 97.8 求形如2y a bx =+的拟合函数及误差平方和，并估计变量y 在21x =处的值。

第5页，共5页解：令a cb ⎛⎫= ⎪⎝⎭，1361162519611144411936A ⎛⎫ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪⎝⎭， ........................... 3分对应的正规方程组TTA Ac A y =为55327271.453277277699369321.5a b ⎛⎫⎛⎫⎛⎫= ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭....................... 7分解此方程得0.972578660.05003512a b ⎛⎫⎛⎫= ⎪ ⎪⎝⎭⎝⎭...................................... 9分即最终的拟合函数 (2)0.972578660.05003512y x =+ 10分误差平方和为 ................................ ()()4220.1502kkk y x y σ==-=∑ 11分2(21)0.972578660.050035122123.0380683y =+⨯= ....... 12分三(本题12分)（1）Gauss 求积公式(()11585()d 0.600.6999f x x f f f -≈-++⎰的代数精度是多少？（2）利用上述公式计算定积分111d 3x x-+⎰，并与真值进行比较；0.60.7745967≈）（3）将区间[1,1]-等分成两个小区间，对每个小区间分别利用（1）中的Gauss 求积公式进行计算，并将最后得到的111d 3x x-+⎰的近似值与（2）进行比较。