1、2008信息与计算科学专业计算方法习题参考解答 江世宏编第一章 绪论姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。1 若误差限为,那么近似数0.003400有几位有效数字？（有效数字的计算）解：，故具有3位有效数字。2 具有4位有效数字的近似值是多少？（有效数字的计算）解：，欲使其近似值具有4位有效数字，必需，即3 已知，是经过四舍五入后得到的近似值，问，有几位有效数字？（有效数字的计算）解：，而，故至少具有2位有效数字。故至少具有2位有效数字。4 设，的相对误差为，求的误差和相对误差？（误差的计算）解：已知，则误差为 则相对误差为 5测得某圆柱体高
2、度的值为，底面半径的值为，已知，求圆柱体体积的绝对误差限与相对误差限。（误差限的计算）解：绝对误差限为相对误差限为6 设的相对误差为,求的相对误差。（函数误差的计算）解：，7计算球的体积，为了使体积的相对误差限为，问度量半径时允许的相对误差限为多大？（函数误差的计算）解：球体积为 ，欲使，必须 。8 设，求证：（1）（2）利用（1）中的公式正向递推计算时误差逐步增大；反向递推计算时误差逐步减小。（计算方法的比较选择）解：如果初始误差为，若是向前递推，有可见，初始误差的绝对值被逐步地扩大了。如果是向后递推，其误差为可见，初始误差的绝对值被逐步减少了。第二章 插值法姓名 学号 班级 习题主要考察点
3、：拉格朗日插值法的构造，均差的计算，牛顿插值和埃尔米特插值构造，插值余项的计算和应用。1 已知，求的拉氏插值多项式。（拉格朗日插值）解法一（待定系数法）：设，由插值条件，有解得：。故 。解法二（基函数法）：由插值条件，有2 已知，用线性插值求的近似值。（拉格朗日线性插值）解：由插值节点与被插函数，可知，其线性插值函数为的近似值为。3 若为互异节点，且有试证明。（拉格朗日插值基函数的性质）解：考虑辅助函数，其中，。是次数不超过的多项式，在节点（）处，有这表明，有n+1个互异实根。故，从而对于任意的均成立。4 已知，用抛物线插值计算的值并估计截断误差。（拉格朗日二次插值）解：由插值条件，其抛物线插
4、值函数为将代入，计算可得：。其余项为： 其中， 故误差的上界为：。5 用余弦函数在，三个节点处的值，写出二次拉格朗日插值多项式, 并近似计算及其绝对误差与相对误差，且与误差余项估计值比较。（拉格朗日二次插值）解：由插值条件，二次拉格朗日插值多项式为绝对误差为：相对误差为：余项为：，其中，其余项的上界为：比较可知，实际计算所得的绝对误差较余项公式所估计出的值要小一些。6 已知函数值，求函数的四阶均差和二阶均差。（均差的计算）解：采用列表法来计算各阶均差，有xy一阶均差二阶均差三阶均差四阶均差0611043461814/34823661/362126529/311/151/15从表中可查得：。xy
5、一阶均差二阶均差48211072/3346186故。其实，根据均差的对称性，该值在第一个表中就可以查到。7 设求之值，其中，而节点互异。（均差的计算）解：由均差可以表示成为函数值的线性组合，有而 ，故。8 如下函数值表012419233建立不超过三次的牛顿插值多项式。（牛顿插值多项式的构造）解：先构造均差表xf(x)一阶均差二阶均差三阶均差0119822314343-10-8-11/4故 。9求一个次数小于等于三次多项式，满足如下插值条件：，。（插值多项式的构造）解法一（待定系数法）：设，则，由插值条件，有解得：。故 解法二（带重节点的均差法）：据插值条件，造差商表xy一阶差商二阶差商三阶差商
6、122422431312852故 10 构造一个三次多项式，使它满足条件（埃尔米特插值）。解：设，利用插值条件，有解得：。11 设。(1)试求在上的三次埃尔米特插值多项式，使得，以升幂形式给出。(2)写出余项的表达式。（埃尔米特插值及其余项的计算）。解：，设，解得：，。故 。，其中，。12 若，试证明： （插值余项的应用）解：以为插值条件，作线性插值多项式，有其余项为故 。13 设求使；又设 ，则估计余项的大小。（插值误差的估计）解：由插值条件，有解得：从而 其余项为第三章 函数逼近姓名 学号 班级 习题主要考察点：最小二乘法，最佳平方逼近，正交多项式的构造。1 设，求于上的线性最佳平方逼近多
7、项式。（最佳平方逼近）解：，法方程组为解得：，线性最佳平方逼近多项式为：。2 令，且设,求使得为于 上的最佳平方逼近多项式。（最佳平方逼近）解：，法方程组为解得：，线性最佳平方逼近多项式为：。3证明：切比雪夫多项式序列在区间上带权正交。（正交多项式的证明）解：对于，有对于，有故，序列在-1，1上带权正交。4求矛盾方程组：的最小二乘解。（最小二乘法）解法一：求与，使得达到最小。于是，令即：，其最小二乘解为：。解法二：，记作，该矛盾方程组的最小二乘解，应满足以下方程组，即解之，得。5 已知一组试验数据22.53455.544.5688.59试用直线拟合这组数据. (计算过程保留3位小数)。（最小二
8、乘线性逼近）解：作矩阵，法方程为即解得：，。其直线拟合函数为。6 用最小二乘原理求一个形如的经验公式,使与下列数据相拟合.19253138 441932.34973.397.8（最小二乘二次逼近）解：等价于对数据表3616259611444 19361932.34973.397.8作线性拟合。其法方程组为：解得：，故经验公式为 。第四章 数值积分姓名 学号 班级 习题主要考察点：代数精度的计算，构造插值型求积公式（梯形，辛甫生公式），复化求积的计算，高斯公式的构造。1给定求积公式试确定使它的代数精度尽可能高。（代数精度的应用和计算）解：分别取，使上述数值积分公式准确成立，有；解得：。故求积公式
9、为。再取，左边=，右边=再取，左边=，右边=此求积公式的最高代数精度为3。2 求积公式，试确定系数,及，使该求积公式具有尽可能高的代数精确度，并给出代数精确度的次数。（代数精度的应用和计算）解：分别取，使求积公式准确成立，有解得：。求积公式为。再取，左边=右边故该求积公式的最高代数精度为2。3数值积分公式，是否为插值型求积公式，为什么？又该公式的代数精确度为多少？（插值型求积公式特征）解：令，故代数精度为1。由于求积节点个数为2，代数精度达到1次，故它是插值型的求积公式。4如果，证明用梯形公式计算积分所得到的结果比准确值大，并说明其几何意义。（梯形求积）解：梯形求积公式是由过点，的线性插值函数
10、在a,b上的定积分。注意到：在区间a,b上，而，有从而。其几何意义可作以下解释：在区间a,b上，故曲线下凹，直线位于曲线之上，因此，曲边梯形的面积小于梯形面积。5用的复化梯形公式计算积分，并估计误差。（复化梯形求积）解：，取求积节点为因，则误差大约为：。6设,则用复化辛甫生公式计算，若有常数使 ，则估计复化辛甫生公式的整体截断误差限。（复化辛甫生公式）解：7已知高斯求积公式 将区间0,1二等分，用复化高斯求积法求定积分的近似值。（高斯公式）解：对于作变量换，有对于作变量换，有8 试确定常数A，B，C和，使得数值积分公式有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少？它是否为高斯型的
11、？（代数精度的应用和计算，高斯点的特征）解：分别取，使上述数值积分公式准确成立，有；整理得：解得：。数值求积公式为再取，左边=，右边=再取，左边=，右边=可见，该数值求积公式的最高代数精度为5。由于该公式中的节点个数为3，其代数精度达到了次，故它是高斯型的。9设是0,1区间上带权的最高次幂项系数为1的正交多项式系（1）求。（2）构造如下的高斯型求积公式。（高斯求积）解（1）：采用施密特正交化方法，来构造带权且在0，1上正交的多项式序列取，设，且它与在0，1上带权正交，于是，故 。设，且它与、在0，1上带权正交，于是，解（2）：的零点为：。设 分别取，使上述求积公式准确成立，有，即解得：，。高斯
12、型求积公式为第五章 非线性方程求根姓名 学号 班级 习题主要考察点：二分法、迭代法、牛顿法和弦截法求根，迭代法求根的收敛性和收敛速度的讨论。1用二分法求方程的正根，要求误差小于0.05。（二分法）解：，在0，2连续，故0，2为函数的有根区间。（1）计算，故有根区间为1，2。（2）计算，故有根区间为。（3）计算，故有根区间为。（4）计算，故有根区间为。（5）计算，故有根区间为。（6）计算，故有根区间为。（7）计算，故有根区间为。（8）若取中点作为取根的近似值，其误差小于取近似根，可满足精度要求。2说明方程 在区间1，2内有惟一根，并选用适当的迭代法求（精确至3位有效数），并说明所用的迭代格式是收
13、敛的。（迭代法）解： ，故函数单调增加，因此，该方程在（1，2）之间存在着惟一的实根。取迭代函数 显然，且故迭代 （）对任意初始值收敛。对于初值，其迭代值分别为，由于，故作为近似值，已精确到了3位有效数字。3设有解方程的迭代法 (1)证明均有（为方程的根）。(2)此迭代法的收敛阶是多少，证明你的结论。 (3) 取用此迭代法求方程根的近似值，误差不超过，列出各次迭代值。（和收敛性讨论）解(1)：，（），故该迭代对任意初值均收敛于方程的根。解(2)：由，故有。，故该迭代的收敛速度是1阶的。解(3):取，代入迭代式，可计算出以下结果：，由于，取可满足精度要求。4设，试证明：由 ，得到的序列收敛于。（
14、收敛性证明）证明：由知，方程有根。由，当时，有，即序列收敛于。5 设方程在0,1内的根为，若采用迭代公式，试证明：均有为方程的根)；此迭代的收敛阶是多少，证明你的结论。（迭代法和收敛性讨论）解：迭代函数，当故迭代在区间上整体收敛。设，则，且故 故该迭代的收敛速度为1阶的。6 方程在附近有根，把方程写成3种不同的等价形式：(1) ，对应迭代格式：(2) ，对应迭代格式：(3) ，对应迭代格式：讨论这些迭代格式在时的收敛性。若迭代收敛，试估计其收敛速度，选一种收敛格式计算出附近的根到4位有效数字。（收敛速度的计算和比较）解：，故方程在上有根。，故方程在上有根。，故方程在上有根。对于迭代式（1）：，
15、 而，故该迭代局部收敛，且收敛速度为1阶的。对于迭代式（2）：在上，又，故该迭代在上整体收敛，且收敛速度为一阶的。对于迭代式（3）：在1，2上的值域为，该迭代式不收敛。取迭代式，进行计算，其结果如下：，取为近似值具有4位有效数字。7设 (1) 写出解 的牛顿迭代格式；(2) 证明此迭代格式是线性收敛的。（牛顿迭代的构造与收敛速度）解：牛顿迭代式为 ，方程的根为，因，故迭代局部收敛。又因，故迭代收敛速度为1阶。8 设计一个计算的牛顿迭代法，且不用除法（其中）。（牛顿迭代法）解：考虑方程，而，该迭代局部收敛。9 用牛顿法求的近似值，取或11为初始值，计算过程保留4位小数。（牛顿迭代的构造）解：考虑
16、方程，取为初始值，计算其迭代值如下：，取为初始值，计算其迭代值如下：，10设是非线性方程的m重根，试证明：迭代法具有至少2阶的收敛速度。（收敛速度证明）解：设是非线性方程的m重根，则，且及，其牛顿迭代函数为牛顿迭代式故该迭代的收敛速度至少是2阶的。11设是非线性方程的m重根，证明：用牛顿迭代法求只是线性收敛。（收敛速度证明）解：设是非线性方程的m重根，则，且及，其牛顿迭代函数为牛顿迭代式故收敛速度为1阶的。12设，在附近有直到阶的连续导数，且，试证：迭代法在附近是阶收敛的。（收敛速度证明）解：将在点附近作泰勒展式，有，其中，在与之间。于是： ，其中，在与之间。由于，故，从而。因此，迭代的收敛速
17、度为p。第六章 常微分方程数值解姓名 学号 班级 习题主要考察点：欧拉方法的构造，单步法的收敛性和稳定性的讨论，线性多步法中亚当姆斯方法的构造和讨论。1 用改进的欧拉公式，求以下微分方程的数值解（取步长），并与精确解作比较。（改进的尤拉公式的应用）解：原方程可转化为 ，令，有解此一阶线性微分方程，可得 。利用以下公式求在节点处的数值解，其中，初值为。MATLAB程序如下：x(1)=0;%初值节点y(1)=1;%初值fprintf(x(%d)=%f,y(%d)=%f,yy(%d)=%fn,1,x(1),1,y(1),1,y(1);for i=1:5 yp=y(i)+0.2*(y(i)-2*x(i
18、)/y(i);%预报值 yc=y(i)+0.2*(yp-2*x(i)/yp);%校正值 y(i+1)=(yp+yc)/2;%改进值 x(i+1)=x(i)+0.2;%节点值 yy(i+1)=sqrt(2*x(i+1)+1);%精确解fprintf(x(%d)=%f,y(%d)=%f,yy(%d)=%fn,i+1,x(i+1),i+1,y(i+1),i+1,yy(i+1);end程序运行的结果如下：x(1)=0.000000, y(1)=1.000000, yy(1)=1.000000x(2)=0.200000, y(2)=1.220000, yy(2)=1.183216x(3)=0.40000
19、0, y(3)=1.420452, yy(3)=1.341641x(4)=0.600000, y(4)=1.615113, yy(4)=1.483240x(5)=0.800000, y(5)=1.814224, yy(5)=1.612452x(6)=1.000000, y(6)=2.027550, yy(6)=1.7320512用四阶龙格库塔法求解初值问题，取, 求时的数值解. 要求写出由直接计算的迭代公式，计算过程保留3位小数。（龙格库塔方法的应用）解：四阶龙格-库塔经典公式为由于，在各点的斜率预报值分别为：四阶经典公式可改写成以下直接的形式：在处，有在处，有注：这两个近似值与精确解在这两点
20、的精确值十分接近。3 用梯形方法解初值问题证明其近似解为并证明当时，它收敛于原初值问题的准确解。解：显然，是原初值问题的准确解。求解一般微分方程初值问题的梯形公式的形式为对于该初值问题，其梯形公式的具体形式为，于是：亦即：注意到：，令，有从而 即：当时，收敛于原初值问题的准确解。4对于初值问题，证明当时，欧拉公式绝对稳定。（显式和隐式欧拉公式的稳定性讨论）证明：显式的欧拉公式为从而，由于，因此，显式欧拉公式绝对稳定。隐式的欧拉公式为，由于，因此，隐式的欧拉公式也是绝对稳定的。5证明：梯形公式无条件稳定。（梯形公式的稳定性讨论）解：对于微分方程初值问题其隐式的梯形公式的具体形式可表示为，从而由，
21、可知，故隐式的梯形公式无条件稳定。6设有常微分方程的初值问题，试用泰勒展开法，构造线性两步法数值计算公式，使其具有二阶精度，并推导其局部截断误差主项。（局部截断误差和主项的计算）解：假设，利用泰勒展式，有又 欲使其具有尽可能高的局部截断误差，必须，从而 ，于是数值计算公式为 。该数值计算公式的局部截断误差的主项为7已知初值问题取步长，利用阿当姆斯公式，求此微分方程在0，10上的数值解，求此公式的局部截断误差的首项。（阿当姆斯公式的应用）解：假设，利用泰勒展开，有，而该阿当姆斯两步公式具有2阶精度，其局部截断误差的主项为。取步长，节点（），注意到，其计算公式可改写为仅需取一个初值，可实现这一公式
22、的实际计算。其MATLAB下的程序如下：x0=0;%初值节点y0=0;%初值for n=0:99 y1=y0+0.02*n+0.01; x1=x0+0.1; fprintf(x(%3d)=%10.8f,y(%3d)=%10.8fn,n+1,x1,n+1,y1); x0=x1; y0=y1;end运行结果如下：x( 1)=0.10000000,y( 1)=0.01000000x( 2)=0.20000000,y( 2)=0.04000000x( 3)=0.30000000,y( 3)=0.09000000x( 4)=0.40000000,y( 4)=0.16000000x( 5)=0.50000
23、000,y( 5)=0.25000000x( 6)=0.60000000,y( 6)=0.36000000x( 7)=0.70000000,y( 7)=0.49000000x( 8)=0.80000000,y( 8)=0.64000000x( 9)=0.90000000,y( 9)=0.81000000x( 10)=1.00000000,y( 10)=1.00000000x( 11)=1.10000000,y( 11)=1.21000000x( 12)=1.20000000,y( 12)=1.44000000x( 13)=1.30000000,y( 13)=1.69000000x( 14)=1
24、.40000000,y( 14)=1.96000000x( 15)=1.50000000,y( 15)=2.25000000x( 16)=1.60000000,y( 16)=2.56000000x( 17)=1.70000000,y( 17)=2.89000000x( 18)=1.80000000,y( 18)=3.24000000x( 19)=1.90000000,y( 19)=3.61000000x( 20)=2.00000000,y( 20)=4.00000000x( 21)=2.10000000,y( 21)=4.41000000x( 22)=2.20000000,y( 22)=4.8
25、4000000x( 23)=2.30000000,y( 23)=5.29000000x( 24)=2.40000000,y( 24)=5.76000000x( 25)=2.50000000,y( 25)=6.25000000x( 26)=2.60000000,y( 26)=6.76000000x( 27)=2.70000000,y( 27)=7.29000000x( 28)=2.80000000,y( 28)=7.84000000x( 29)=2.90000000,y( 29)=8.41000000x( 30)=3.00000000,y( 30)=9.00000000x( 31)=3.1000
26、0000,y( 31)=9.61000000x( 32)=3.20000000,y( 32)=10.24000000x( 33)=3.30000000,y( 33)=10.89000000x( 34)=3.40000000,y( 34)=11.56000000x( 35)=3.50000000,y( 35)=12.25000000x( 36)=3.60000000,y( 36)=12.96000000x( 37)=3.70000000,y( 37)=13.69000000x( 38)=3.80000000,y( 38)=14.44000000x( 39)=3.90000000,y( 39)=1
27、5.21000000x( 40)=4.00000000,y( 40)=16.00000000x( 41)=4.10000000,y( 41)=16.81000000x( 42)=4.20000000,y( 42)=17.64000000x( 43)=4.30000000,y( 43)=18.49000000x( 44)=4.40000000,y( 44)=19.36000000x( 45)=4.50000000,y( 45)=20.25000000x( 46)=4.60000000,y( 46)=21.16000000x( 47)=4.70000000,y( 47)=22.09000000x(
28、 48)=4.80000000,y( 48)=23.04000000x( 49)=4.90000000,y( 49)=24.01000000x( 50)=5.00000000,y( 50)=25.00000000x( 51)=5.10000000,y( 51)=26.01000000x( 52)=5.20000000,y( 52)=27.04000000x( 53)=5.30000000,y( 53)=28.09000000x( 54)=5.40000000,y( 54)=29.16000000x( 55)=5.50000000,y( 55)=30.25000000x( 56)=5.60000
29、000,y( 56)=31.36000000x( 57)=5.70000000,y( 57)=32.49000000x( 58)=5.80000000,y( 58)=33.64000000x( 59)=5.90000000,y( 59)=34.81000000x( 60)=6.00000000,y( 60)=36.00000000x( 61)=6.10000000,y( 61)=37.21000000x( 62)=6.20000000,y( 62)=38.44000000x( 63)=6.30000000,y( 63)=39.69000000x( 64)=6.40000000,y( 64)=4
30、0.96000000x( 65)=6.50000000,y( 65)=42.25000000x( 66)=6.60000000,y( 66)=43.56000000x( 67)=6.70000000,y( 67)=44.89000000x( 68)=6.80000000,y( 68)=46.24000000x( 69)=6.90000000,y( 69)=47.61000000x( 70)=7.00000000,y( 70)=49.00000000x( 71)=7.10000000,y( 71)=50.41000000x( 72)=7.20000000,y( 72)=51.84000000x(
31、 73)=7.30000000,y( 73)=53.29000000x( 74)=7.40000000,y( 74)=54.76000000x( 75)=7.50000000,y( 75)=56.25000000x( 76)=7.60000000,y( 76)=57.76000000x( 77)=7.70000000,y( 77)=59.29000000x( 78)=7.80000000,y( 78)=60.84000000x( 79)=7.90000000,y( 79)=62.41000000x( 80)=8.00000000,y( 80)=64.00000000x( 81)=8.10000
32、000,y( 81)=65.61000000x( 82)=8.20000000,y( 82)=67.24000000x( 83)=8.30000000,y( 83)=68.89000000x( 84)=8.40000000,y( 84)=70.56000000x( 85)=8.50000000,y( 85)=72.25000000x( 86)=8.60000000,y( 86)=73.96000000x( 87)=8.70000000,y( 87)=75.69000000x( 88)=8.80000000,y( 88)=77.44000000x( 89)=8.90000000,y( 89)=7
33、9.21000000x( 90)=9.00000000,y( 90)=81.00000000x( 91)=9.10000000,y( 91)=82.81000000x( 92)=9.20000000,y( 92)=84.64000000x( 93)=9.30000000,y( 93)=86.49000000x( 94)=9.40000000,y( 94)=88.36000000x( 95)=9.50000000,y( 95)=90.25000000x( 96)=9.60000000,y( 96)=92.16000000x( 97)=9.70000000,y( 97)=94.09000000x(
34、 98)=9.80000000,y( 98)=96.04000000x( 99)=9.90000000,y( 99)=98.01000000x(100)=10.00000000,y(100)=100.00000000第七章 线性方程组的迭代解法姓名 学号 班级 习题主要考察点：雅可比、高斯-塞德尔迭代法解线性方程组，及其收敛性讨论。1证明：迭代格式收敛，其中。（迭代法收敛性判断）解：因，故迭代收敛。2若用雅可比迭代法求解方程组迭代收敛的充要条件是。（雅可比迭代法的收敛性）解：原线性方程组的等价方程组为其雅可比迭代式为其收敛的充要条件是，即。3 用雅可比、高斯-塞德尔迭代法，求解方程组是否收敛？
35、为什么？若将方程组改变成为再用上述两种迭代法求解是否收敛？为什么？（雅可比、高斯-塞德尔迭代法的收敛性）解：雅可比迭代式为其，故雅可比迭代发散。高斯-塞德尔迭代式为其，故高斯-塞德尔迭代发散。对于线性方程组，即，其雅可比迭代为，其，故雅可比迭代收敛。，其，故高斯-塞德尔迭代收敛。4证明解线性方程组的雅可比迭代收敛，其中。（雅可比迭代收敛性判断）解：雅可比迭代为，其，故雅可比迭代收敛。5已知方程组，其中，(1) 试讨论用雅可比迭代法和高斯-塞德尔迭代法求解此方程组的收敛性。(2) 若有迭代公式，试确定的取值范围，使该迭代公式收敛。（雅可比迭代法、高斯-塞德尔迭代法和一般迭代法的收敛性讨论）解：雅
36、可比迭代式为其，故雅可比迭代收敛。高斯-塞德尔迭代式为其，故高斯-塞德尔迭代收敛。对于以下迭代式故的特征值为，。当时，有，从而迭代收敛。6给出矩阵，（为实数），试分别求出的取值范围：(1) 使得用雅可比迭代法解方程组时收敛；(2) 使得用高斯-塞德尔迭代法解方程组时收敛。（雅可比、高斯-塞德尔迭代法及收敛性讨论）解：雅可比迭代为当时，使雅可比迭代收敛。高斯-塞德尔迭代为仍然是当时，使高斯-塞德尔迭代收敛。7设，(1) 设是由雅可比迭代求解方程组所产生的迭代向量，且，试写出计算的精确表达式。(2) 设是的精确解，写出误差的精确表达式。(3) 如构造如下的迭代公式解方程组，试确定的范围，使迭代收敛
37、。（雅可比迭代及其收敛判断）解：原线性方程组等价于，其雅可比迭代为将上述迭代式记作，从而而，若记，则，于是，当为偶数时，当为奇数时，总之，。的特征多项式为，故其特征值为1，3。对于迭代其迭代矩阵为，其特征值为，。当时，该迭代收敛。8对于给定的线性方程组（1）讨论雅可比迭代法与高斯-塞德尔迭代法的收敛性。（2）对收敛的方法，取初值，迭代两次，求出。（雅可比,高斯-塞德尔迭代法的计算和比较）解：，雅可比迭代式为：，取计算迭代阵的特征值故，雅可比迭代收敛。高斯-塞德尔迭代式，取计算迭代阵的特征值故，高斯-塞德尔迭代发散。对于雅可比迭代，取，可得，进一步的计算可知，。9 证明对称矩阵当为正定矩阵，且只
38、有当时，用雅可比迭代法求解方程组才收敛。（雅可比迭代法的收敛性）解：矩阵A的各级顺序主子式分别为，当时，上述各级顺序主子式均大于零，故A正定。其雅可比迭代式为其迭代矩阵的特征多项式为迭代矩阵的特征值为和。当时，雅可比迭代收敛。第八章 线性方程组的直接解法姓名 学号 班级 习题主要考察点：高斯消去法，LU分解法，平方根法和追赶法解线性方程组。1用高斯消去法解方程组。(高斯消去法的应用)解：用，依次左乘方程组两边（采用高斯消去法），有解得：，。2用LU分解法求解线性方程组。（LU分解法的应用）解：原线性方程组的系数矩阵，右端列向量分别为，其中，解，可得：，解，可得：。3设，求A的LU分解。（LU分
39、解法的应用）解：，4试用“追赶法”解方程组，其中：，（追赶法的应用）解：追的过程：由第1个方程，有代入第2个方程，解出再代入第3个方程，解出。赶的过程：将代入，得到。再将代入，得。故原线性方程组的解为。5设，求（条件数的计算）解：的特征值分别为2，7，故。从而，故。6求证：，（范数的性质）证明：设和为n阶方阵，对于任意给定n维向量，有故 。由有：，从而 。7求证：。（范数的性质）证明：，设的最大特征值为，其特征向量为，且从而。，故。8对矩阵，求，和。（范数，条件数的计算）解：行范数，列范数，且为实对称矩阵。其特征值分别为，。的特值值分别为，。故，。9方程组，其中，是对称的且非奇异。设有误差，则原方程组变化为，其中为解的误差向量，试证明：，其中和分别为的按模最大和最小的特征值。（范数的性质，误差的分析）证明：是对称的且非奇异，则，于是，所证明的不等式等价于，亦即由，有 ，从而，注意到，有于是10证明：若为严格对角占优矩阵，则非奇异。（严格对角占优矩阵的性质）解：据严格对角占优，有用反证法，假设为奇异阵，则存在着非零向量，使得。由于，不妨设，则，于是，这显然与条件相矛盾，故为非奇异阵。42
