下载文档原格式
习题一（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 迭代法收敛.。
习题一1、取3.14,3.15，722，113355作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
解：14.31=x312110211021--⨯=⨯≤-x π 所以，1x 有三位有效数字绝对误差：14.3-=πe ，相对误差：ππ14.3-=r e 绝对误差限：21021-⨯≤ε，相对误差限：213106110321-+-⨯=⨯⨯=r ε21122105.0105.01084074.000840174.015.315.3---⨯=⨯≤⨯==-=πx所以，2x 有两位有效数字绝对误差：15.3-=πe ，相对误差：ππ15.3-=r e 绝对误差限：11021-⨯=ε，相对误差限：11061-⨯=r ε31222105.0105.01012645.00012645.0722722---⨯=⨯≤⨯==-=πx所以，3x 有三位有效数字绝对误差：722-=πe ，相对误差：ππ722-=r e绝对误差限：21021-⨯=ε，相对误差限：21061-⨯=r ε1133551=x 7166105.0105.01032.000000032.0113355---⨯=⨯≤⨯==-π 所以，4x 有七位有效数字 绝对误差：113355-=πe ，相对误差：ππ113355-=r e绝对误差限：61021-⨯=ε，相对误差限：61061-⨯=r ε3、下列各数都是对准确数四舍五入后得到的近似数，试分别指出它们的绝对误差限和相对误差限，有效数字的位数。
5000,50.31,3015.0,0315.04321====x x x x解：0315.01=x m=-13141*10211021---⨯=⨯≤-x x 所以，n=3，1x 有三位有效数字绝对误差限：41021-⨯=ε，相对误差：2110611021-+-=⨯=n r a ε3015.02=x m=04042*10211021--⨯=⨯≤-x x所以，n=4，1x 有四位有效数字绝对误差限：41021-⨯=ε，相对误差：3110611021-+-=⨯=n r a ε50.313=x m=24223*10211021--⨯=⨯≤-x x所以，n=4，1x 有四位有效数字绝对误差限：21021-⨯=ε，相对误差：3110611021-+-=⨯=n r a ε50004=x m=44404*10211021-⨯=⨯≤-x x所以，n=4，1x 有四位有效数字绝对误差限：5.010210=⨯=ε，相对误差：23110105211021--+-=⨯=⨯=n r a ε 4、计算10的近似值，使其相对误差不超过%1.0。
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 是精确值ａ的近似值（1） ａ=π，x=3.1 （2） ａ=1/7，x=0.143 （3） ａ=π/1000，x=0.0031 （4） ａ=100/7，x=14.3 试估计x 的绝对误差和相对误差。
解：（1） e=∣3.1-π∣≈0.0416， δr = e/∣x ∣≈0.0143 （2） e=∣0.143-1/7∣≈0.0143 δr = e/∣x ∣≈0.1 （3） e=∣0.0031-π/1000∣≈0.0279 δr = e/∣x ∣≈0.9 (4) e=∣14.3-100/7∣≈0.0143 δr = e/∣x ∣≈0.0012. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。
解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［-x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.497073. 设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑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),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析课后习题部分参考答案Chapter 1（P10）5. 求2的近似值*x ，使其相对误差不超过%1.0。
解： 4.12=。
设*x 有n 位有效数字，则nx e -⨯⨯≤10105.0|)(|*。
从而，1105.0|)(|1*nr x e -⨯≤。
故，若%1.0105.01≤⨯-n，则满足要求。
解之得，4≥n 。
414.1*=x 。
（P10）7. 正方形的边长约cm 100，问测量边长时误差应多大，才能保证面积的误差不超过12cm 。
解：设边长为a ，则cm a 100≈。
设测量边长时的绝对误差为e ，由误差在数值计算的传播，这时得到的面积的绝对误差有如下估计:e ⨯⨯≈1002。
按测量要求，1|1002|≤⨯⨯e 解得，2105.0||-⨯≤e 。
Chapter 2（P47）5. 用三角分解法求下列矩阵的逆矩阵：⎪⎪⎪⎭⎫ ⎝⎛--=011012111A 。
解：设()γβα=-1A。
分别求如下线性方程组：⎪⎪⎪⎭⎫ ⎝⎛=001αA ，⎪⎪⎪⎭⎫ ⎝⎛=010βA ，⎪⎪⎪⎭⎫⎝⎛=100γA 。
先求A 的LU 分解（利用分解的紧凑格式），⎪⎪⎪⎭⎫⎝⎛-----3)0(2)1(1)1(2)0(1)1(2)2(1)1(1)1(1)1(。
即，⎪⎪⎪⎭⎫ ⎝⎛=121012001L ，⎪⎪⎪⎭⎫⎝⎛---=300210111U 。
经直接三角分解法的回代程，分别求解方程组，⎪⎪⎪⎭⎫ ⎝⎛=001Ly 和y U =α，得，⎪⎪⎪⎭⎫ ⎝⎛-=100α；⎪⎪⎪⎭⎫ ⎝⎛=010Ly 和y U =β，得，⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛=323131β；⎪⎪⎪⎭⎫ ⎝⎛=100Ly 和y U =γ，得，；⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--=313231γ。
所以，⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛---=-3132132310313101A 。
（P47）6. 分别用平方根法和改进平方根法求解方程组：⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----816211515311401505231214321x x x x 解: 平方根法：先求系数矩阵A 的Cholesky 分解（利用分解的紧凑格式），⎪⎪⎪⎪⎪⎭⎫ ⎝⎛----1)15(2)1(1)5(3)3(3)14(2)0(1)1(1)5(2)2(1)1(，即，⎪⎪⎪⎪⎪⎭⎫⎝⎛--=121332100120001L ，其中，TL L A ⨯=。
第一章 绪论3．下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位，试指出它们是几位有效数字：**1 1.1021x =,**20.031x =, **3385.6x =, **456.430x =,**57 1.0.x =´解：*1 1.1021x =是五位有效数字；是五位有效数字；*20.031x =是二位有效数字；是二位有效数字； *3385.6x =是四位有效数字；是四位有效数字；*456.430x =是五位有效数字；是五位有效数字； *57 1.0.x =´是二位有效数字。
是二位有效数字。
4．利用公式(2.3)求下列各近似值的误差限：(1) ***124x x x ++,(2) ***123x x x ,(3) **24/x x . 其中****1234,,,x x x x 均为第3题所给的数。
题所给的数。
解：解：*41*32*13*34*151()1021()1021()1021()1021()102x x x x x e e e e e -----=´=´=´=´=´ ***124***1244333(1)()()()()1111010102221.0510x x x x x x e e e e ----++=++=´+´+´=´***123*********123231132143(2)()()()()1111.10210.031100.031385.610 1.1021385.6102220.215x x x x x x x x x x x x e e e e ---=++=´´´+´´´+´´´»**24****24422*4335(3)(/)()()110.0311056.430102256.43056.43010x x x x x x x e e e ---+»´´+´´=´= 6．设028Y =，按递推公式11783100n n Y Y -=- （n=1,2,n=1,2,……）计算到100Y 。
数值分析参考答案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．略。
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x -==≈⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x --==≈-⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.15＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.141592653.1428571430.0012644930.00137e x π=-=-=-≈-相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯有效数字： 因为π＝3.14159265…=0.314159265…×10， 223.1428571430.3142857143107==⨯，m=1。
数值分析课后习题及答案第一章绪论（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使得达到最小。
对求导得到：，令，得到，这说明取平均值在最小二乘意义下误差达到最小。
东南大学数值分析上机题答案-CAL-FENGHAI.-(YICAI)-Company One1数值分析上机题第一章17.（上机题）舍入误差与有效数 设∑=-=Nj N j S 2211，其精确值为）111-23（21+-N N 。
（1）编制按从大到小的顺序1-1···1-311-21222N S N +++=，计算N S 的通用程序；（2）编制按从小到大的顺序121···1)1(111222-++--+-=N N S N ，计算N S 的通用程序；（3）按两种顺序分别计算210S ,410S ,610S ,并指出有效位数（编制程序时用单精度）；（4）通过本上机题，你明白了什么？ 解： 程序：（1）从大到小的顺序计算1-1···1-311-21222N S N +++=： function sn1=fromlarge(n) %从大到小计算sn1format long ; sn1=single(0); for m=2:1:nsn1=sn1+1/(m^2-1); end end（2）从小到大计算121···1)1(111222-++--+-=NN S Nfunction sn2=fromsmall(n) %从小到大计算sn2format long ; sn2=single(0); for m=n:-1:2sn2=sn2+1/(m^2-1); end end （3）总的编程程序为：function p203()clear allformat long;n=input('please enter a number as the n:') sn=1/2*(3/2-1/n-1/(n+1));%精确值为snfprintf('精确值为%f\n',sn);sn1=fromlarge(n);fprintf('从大到小计算的值为%f\n',sn1);sn2=fromsmall(n);fprintf('从小到大计算的值为%f\n',sn2);function sn1=fromlarge(n) %从大到小计算sn1format long;sn1=single(0);for m=2:1:nsn1=sn1+1/(m^2-1);endendfunction sn2=fromsmall(n) %从小到大计算sn2format long;sn2=single(0);for m=n:-1:2sn2=sn2+1/(m^2-1);endendend运行结果：从而可以得到N值真值顺序值有效位数2100.740050 从大到小0.740049 5从小到大0.740050 6 4100.749900 从大到小0.749852 3从小到大0.749900 6 6100.749999 从大到小0.749852 3从小到大0.749999 6 （4）感想：通过本上机题，我明白了，从小到大计算数值的精确位数比较高而且与真值较为接近，而从大到小计算数值的精确位数比较低。
