下载文档原格式
习题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 =。
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
第一章题12给定节点01x =−，11x =，23x =，34x =，试分别对下列函数导出拉格朗日插值余项：（1）（1）3()432f x x x =−+（2）（2）43()2f x x x =−解（1）(4)()0f x =，由拉格朗日插值余项得(4)0123()()()()()()()04!f f x p x x x x x x x x x ξ−=−−−−=；（2）(4)()4!f x =由拉格朗日插值余项得01234!()()()()()()4!f x p x x x x x x x x x −=−−−−(1)(1)(3)(4)x x x x =+−−−.题15证明：对于()f x 以0x ，1x 为节点的一次插值多项式()p x ，插值误差01210()()()max ()8x x x x x f x p x f x ≤≤−′′−≤.证由拉格朗日插值余项得01()()()()()2!f f x p x x x x x ξ′′−=−−，其中01x x ξ≤≤，010101max ()()()()()()()()2!2!x x x f x f f x p x x x x x x x x x ξ≤≤′′′′−=−−≤−−01210()max ()8x x x x x f x ≤≤−′′≤.题22采用下列方法构造满足条件(0)(0)0p p ′==，(1)(1)1p p ′==的插值多项式()p x ：（1）（1）用待定系数法；（2）（2）利用承袭性，先考察插值条件(0)(0)0p p ′==，(1)1p =的插值多项式()p x .解（1）有四个插值条件，故设230123()p x a a x a x a x =+++，2123()23p x a a x a x ′=++，代入得方程组001231123010231a a a a a a a a a =⎧⎪+++=⎪⎨=⎪⎪++=⎩解之，得01230021a a a a =⎧⎪=⎪⎨=⎪⎪=−⎩23()2p x x x ∴=−；（2）先求满足插值条件(0)(0)0p p ′==，(1)1p =的插值多项式()p x ，由0为二重零点，可设2()p x ax =，代入(1)1p =，得1a =，2()p x x ∴=；再求满足插值条件(0)(0)0p p ′==，(1)(1)1p p ′==的插值多项式()p x ，可设22()(1)p x x bx x =+−，2()22(1)p x x bx x bx ′=+−+∵，代入(1)1p ′=，得1b =−，2223()(1)2p x x x x x x ∴=−−=−.题33设分段多项式323201()2112x x x S x x bx cx x ⎧+≤≤=⎨++−≤≤⎩是以0,1,2为节点的三次样条函数，试确定系数,b c 的值.解由(1)2S =得212b c ++−=，1b c ∴+=；223201()6212x x x S x x bx c x ⎧+<<′=⎨++<<⎩，由(1)5S ′=得625b c ++=，21b c ∴+=−；联立两方程，得2,3b c =−=，且此时6201()12212x x S x x b x +<<⎧′′=⎨+<<⎩，(1)8(1)S S −+′′′′==，()S x 是以0,1,2为节点的三次样条函数.题35用最小二乘法解下列超定方程组：24113532627x y x y x y x y +=⎧⎪−=⎪⎨+=⎪⎪+=⎩.解记残差的平方和为2222(,)(2411)(353)(26)(27)f x y x y x y x y x y =+−+−−++−++−令00f x f y ∂⎧=⎪∂⎪⎨∂⎪=∂⎪⎩，得3661020692960x y x y −−=⎧⎨−+−=⎩，解之得83027311391x y ⎧=⎪⎪⎨⎪=⎪⎩.题37用最小二乘法求形如2y a bx =+的多项式，使与下列数据相拟合：x1925313844y19.032.349.073.397.8解拟合曲线中的基函数为0()1x ϕ=，20()x x ϕ=，其法方程组为0001010001(,)(,)(,)(,)(,)(,)f a f b ϕϕϕϕϕϕϕϕϕϕ⎛⎞⎛⎞⎛⎞=⎜⎟⎜⎟⎜⎟⎝⎠⎝⎠⎝⎠，其中00(,)5ϕϕ=，0110(,)(,)5327ϕϕϕϕ==，11(,)7277699ϕϕ=，0(,)271.4f ϕ=，1(,)369321.5f ϕ=，解之得5320.97265472850.055696a b ⎧==⎪⎪⎨⎪==⎪⎩，20.97260.05y x ∴=+.第二章题3确定下列求积公式中的待定参数，使其代数精度尽量地高，并指明求积公式所具有的代数精度：（2）10120113()(()()424f x dx A f A f A f ≈++∫（2）从结论“在机械求积公式中，代数精度最高的是插值型的求积公式”出发，11000013()(224()11133()()4244x x A l x dx dx −−===−−∫∫，11110013()()144()11133()()2424x x A l x dx dx −−===−−−∫∫，11220011()242()31313()4442x x A l x dx dx −−===−−∫∫，10211123()()()(343234f x dx f f f ∴≈−+∫，当3()f x x =时，有左边＝113001()d d 4f x x x x ==∫∫，右边＝3332111232111231()()()()()()3432343432344f f f −+=⋅−⋅+⋅=，左边＝右边，当4()f x x =时，有左边＝114001()d d 5f x x x x ==∫∫，右边＝44421112321112337()()()()()()343234343234192f f f −+=⋅−⋅+⋅=，左边≠右边，所以该求积公式的代数精度为3.题8已知数据表x 1.11.3 1.5xe3.00423.66934.4817试分别用辛甫生法与复化梯形法计算积分 1.51.1x e dx∫.解辛甫生法1.51.1xe dx ∫()1.5 1.13.00424 3.66934.4817 1.477546−≈+×+=；复化梯形法1.51.1xe dx ∫()0.23.00422 3.66934.4817 1.482452≈+×+=.题17用三点高斯公式求下列积分值12041dxx π=+∫.解先做变量代换，设)(1+21=t x ，则1204d 1x x +∫＝112112418d d 124(1)1(1)4t t t t −−⋅=++++∫∫()2225888589994014141≈×+×+×++⎛⎞⎞++⎜⎟⎟⎝⎠⎠3.141068=.第三章用欧拉方法求解初值问题y ax b ′=+，(0)0y =：（1）试导出近似解n y的显式表达式；解（1）其显示的Euler 格式为：11111(,)()n n n n n n y y hf x y y h ax b −−−−−=+=+⋅+故122()n n n y y h ax b −−−=+⋅+⋯⋯100()y y h ax b =+⋅+将上组式子左右累加，得0021()n n n y y ah x x x nhb−−=+++++⋯(02(2)(1))ah h h n h n h nhb =+++−+−+⋯2(1)/2ah n n nhb=−+题10选取参数p 、q ，使下列差分格式具有二阶精度：1111(,)n n n n y y hK K f x ph y qhK +=+⎧⎨=++⎩.解将1K 在点(,)n n x y 处作一次泰勒展开，得11(,)n n K f x ph y qhK =++21(,)(,)(,)()n n x n n y n n f x y phf x y qhK f x y O h =+++()221(,)(,)(,)(,)(,)()(,)()n n x n n n n x n n y n n y n n f x y phf x y qh f x y phf x y qhK f x y O h f x y O h =++++++2(,)(,)(,)(,)()n n x n n n n y n n f x y phf x y qhf x y f x y O h =+++代入，得()21(,)(,)(,)(,)()n n n n x n n n n y n n y y h f x y phf x y qhf x y f x y O h +=++++2231(,)(,)(,)(,)()n n n n x n n n n y n n y y hf x y ph f x y qh f x y f x y O h +=++++而231()()()()()()2n n n n n h y x y x h y x hy x y x O h +′′′=+=+++23()(,())(,())(,())(,())()2n n n x n n n n y n n h y x hf x y x f x y x f x y x f x y x O h ⎡⎤=++++⎣⎦考虑其局部截断误差，设()n n y y x =，比较上两式，当12p =，12q =时，311()()n n y x y O h ++−=.第四章题2证明方程1cos 2x x=有且仅有一实根；试确定这样的区间[,]a b ，使迭代过程11cos 2k kx x +=对一切0[,]x a b ∈均收敛.解设1()cos 2f x x x=−，则()f x 在区间(,)−∞+∞上连续，且11(0)cos 0022f =−=−<，1(cos 022222f ππππ=−=>，所以()f x 在[0,]2π上至少有一根；又1()1sin 02f x x ′=+>，所以()f x 单调递增，故()f x 在[0,]2π上仅有一根.迭代过程11cos 2k k x x +=，其迭代函数为1()cos 2g x x=，[0,]2x π∀∈，110()cos 222g x x π≤=≤≤，()[0,]2g x π∴∈；1()sin 2g x x ′=−，1()12g x ′≤<，由压缩映像原理知0[0,2x π∀∈，11cos 2k kx x +=均收敛.注这里取[,]a b 为区间[0,]2π，也可取[,]a b 为区间(,)−∞+∞等.题5考察求解方程1232cos 0x x −+=的迭代法124cos 3k kx x +=+（１）（１）证明它对于任意初值0x 均收敛；（２）证明它具有线性收敛性；证（1）迭代函数为2()4cos 3g x x=+，(,)x ∀∈−∞+∞，()(,)g x ∈−∞+∞；又22()sin 133g x x ′=−≤<，由压缩映像原理知0x ∀，124cos 3k k x x +=+均收敛；（2）***1*2lim ()sin 03k k k x x g x x x x +→∞−′==−≠−（否则，若*sin 0x =，则*,x m m Z π=∈，不满足方程），所以迭代124cos 3k kx x +=+具有线性收敛速度；题7求方程3210x x −−=在0 1.5x =附近的一个根，证明下列两种迭代过程在区间[1.3,1.6]上均收敛：（１）（１）改写方程为211x x =+，相应的迭代公式为1211k k x x +=+；（２）（２）改写方程为321x x =+，相应的迭代公式为1k x +=解（1）3232211011x x x x x x −−=⇔=+⇔=+，迭代公式为1211k k x x +=+，其迭代函数为21()1g x x =+[1.3,1.6]x ∀∈，2221111.3 1.3906111 1.5917 1.61.6 1.3x ≤≈+≤+≤+≈<，()[1.3,1.6]g x ∴∈；又32()g x x ′=−，333222-0.9103==-0.48831.3 1.6x −−−≤≤，()0.91031g x ′≤<，由大范围收敛定理知0[1.3,1.6]x ∀∈，1211k k x x +=+均收敛；（2）3232101x x x x x −−=⇔=+⇔=1k x +=其迭代函数为()g x =[1.3,1.6]x ∀∈，1.3 1.3908 1.5269 1.6≤≈≤≤≈<，()[1.3,1.6]g x ∴∈；又()g x ′=，00.4912≤≤≤=，()0.49121g x ′≤<，由大范围收敛定理知0[1.3,1.6]x ∀∈，1k x +=均收敛.题5分别用雅可比迭代与高斯-塞德尔迭代求解下列方程组：1231231235325242511x x x x x x x x x +−=⎧⎪−+=⎨⎪+−=−⎩（2）其雅可比迭代格式为(1)()()123(1)()()213(1)()()312253512221121555k k k k k k k k k x x x x x x x x x +++⎧⎪=−+⎪⎪=−++⎨⎪⎪=++⎪⎩，取初始向量(0)000x ⎛⎞⎜⎟=⎜⎟⎜⎟⎝⎠，迭代发散；其高斯-塞德尔迭代格式为(1)()()123(1)(1)()213(1)(1)(1)312253512221121555k k k k k k k k k x x x x x x x x x ++++++⎧⎪=−+⎪⎪=−++⎨⎪⎪=++⎪⎩，取初始向量(0)000x ⎛⎞⎜⎟=⎜⎟⎜⎟⎝⎠，迭代发散.第六章题2用主元消去法解下列方程组）12312312323553476335x x x x x x x x x ++=⎧⎪++=⎨⎪++=⎩解（2）对其增广矩阵进行列主元消元得23553476347634763476235501/31/3105/32/331335133505/32/3301/31/31⎛⎞⎛⎞⎛⎞⎛⎞⎜⎟⎜⎟⎜⎟⎜⎟→→→⎜⎟⎜⎟⎜⎟⎜⎟⎜⎟⎜⎟⎜⎟⎜⎟⎝⎠⎝⎠⎝⎠⎝⎠347605/32/33001/52/5⎛⎞⎜⎟→⎜⎟⎜⎟⎝⎠回代求解上三角方程组1232333476523331255x x x x x x ⎧⎪++=⎪⎪+=⎨⎪⎪=⎪⎩得321214x x x =⎧⎪=⎨⎪=−⎩，所以412x −⎛⎞⎜⎟=⎜⎟⎜⎟⎝⎠.。
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 是精确值 a 的近似值。
3.14,7/100)4(143.0,7/1)2(0031.0,1000/)3(1.3,)1(========x a x a x a x a ππ试估计x 的绝对误差和相对误差。
解：（1）0132.00416.01.3≈=≈−=−=aee x a e r π （2）0011.00143.0143.07/1≈=≈−=−=a ee x a e r （3）0127.000004.00031.01000/≈=≈−=−=aee x a e r π （4）001.00143.03.147/100≈=≈−=−=aee x a e r2、已知四个数：001.0,25.134,0250.0,3.264321====x x x x 。
试估计各近似数的有效位数和误差限，并估计运算3211x x x =μ和1431/x x x =μ的相对误差限。
解：21111121101901.0,1021,3,10263.06.23−−⨯≈=⨯==⨯==x x x x n x r δδδ22214212102.0,1021,3,10250.00250.0−−−⨯≈=⨯==⨯==x x x x n x r δδδ 43332333103724.0,1021,5,1013425.025.134−−⨯≈=⨯==⨯==x x x x n x r δδδ 5.0,1021,1,101.0001.04443424==⨯==⨯==−−x x x x n x r δδδ 由相对误差限公式：i r ini n in ni i ir x x fx x f x x x f x x f u δδδ∂∂=∂∂=∑∑==1111),,(),,()(所以有：232123113211103938.0)(1)(−⨯≈++=x x x x x x x x x r δδδμμδ4971.0)(1)(4133141214311≈++−=x x x x x x x x x x r δδδμμδ 3、设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。
