还剩51页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第四版数值分析习题第一章绪论
1.设x0x的相对误差为δ求的误差.
1.设x的相对误差为2%求的相对误差.
1.下列各数都是经过四舍五入得到的近似数即误差限不超过最后一位的半个单位试指出它们是几位有效数字:
1.利用公式
3.3求下列各近似值的误差限:其中均为第3题所给的数.
1.计算球体积要使相对误差限为1%问度量半径R时允许的相对误差限是多少
1.设按递推公式n=12…计算到.若取≈
27.982五位有效数字试问计算将有多大误差
1.求方程的两个根使它至少具有四位有效数字≈
27.
982.
1.当N充分大时怎样求
1.正方形的边长大约为100㎝应怎样测量才能使其面积误差不超过1㎝
1.设假定g是准确的而对t的测量有±
0.1秒的误差证明当t增加时S的绝对误差增加而相对误差却减小.
1.序列满足递推关系n=12…若三位有效数字计算到时误差有多大这个计算过程稳定吗
1.计算取利用下列等式计算哪一个得到的结果最好
1.求f30的值.若开平方用六位函数表问求对数时误差有多大若改用另一等价公式计算求对数时误差有多大
1.试用消元法解方程组假定只用三位数计算问结果是否可靠
1.已知三角形面积其中c为弧度且测量abc的误差分别为证明面积的误差满足第二章插值法
1.根据
2.2定义的范德蒙行列式令证明是n次多项式它的根是且.
1.当x=1-12时fx=0-34求fx的二次插值多项式.
1.给出fx=lnx的数值表用线性插值及二次插值计算ln
0.54的近似值.x
0.
40.
50.
60.
70.8lnx-
0.916291-
0.693147-
0.510826-
0.357765-
0.
2231441.给出cosx0°≤x≤90°的函数表步长h=1′=1/60°若函数表具有5位有效数字研究用线性插值求cosx近似值时的总误差界.
1.设k=0123求.
1.设为互异节点j=01…n求证:ii
1.设且求证
1.在上给出的等距节点函数表若用二次插值求的近似值要使截断误差不超过问使用函数表的步长应取多少
1.若求及.
1.如果是次多项式记证明的阶差分是次多项式并且为正整数.
1.证明.
1.证明
1.证明
1.若有个不同实根证明
1.证明阶均差有下列性质:i若则;i若则.
1.求及.
1.证明两点三次埃尔米特插值余项是并由此求出分段三次埃尔米特插值的误差限.
1.求一个次数不高于4次的多项式使它满足并由此求出分段三次埃尔米特插值的误差限.
1.试求出一个最高次数不高于4次的函数多项式以便使它能够满足以下边界条件.
1.设把分为等分试构造一个台阶形的零次分段插值函数并证明当时在上一致收敛到.
1.设在上取按等距节点求分段线性插值函数计算各节点间中点处的与的值并估计误差.
1.求在上的分段线性插值函数并估计误差.
1.求在上的分段埃尔米特插值并估计误差.
1.给定数据表如下:
0.
250.
300.
390.
450.
530.
50000.
54770.
62450.
67080.7280试求三次样条插值并满足条件ii
1.若是三次样条函数证明i;i若式中为插值节点且则.
1.编出计算三次样条函数系数及其在插值节点中点的值的程序框图可用
8.7式的表达式.第三章函数逼近与计算
1.a利用区间变换推出区间为的伯恩斯坦多项式.b对在上求1次和三次伯恩斯坦多项式并画出图形并与相应的马克劳林级数部分和误差做比较.
1.求证:a当时.b当时.
1.在次数不超过6的多项式中求在的最佳一致逼近多项式.
1.假设在上连续求的零次最佳一致逼近多项式.
1.选取常数使达到极小又问这个解是否唯一
1.求在上的最佳一次逼近多项式并估计误差.
1.求在上的最佳一次逼近多项式.
1.如何选取使在上与零偏差最小是否唯一
1.设在上求三次最佳逼近多项式.
1.令求.
1.试证是在上带权的正交多项式.
1.在上利用插值极小化求1的三次近似最佳逼近多项式.
1.设在上的插值极小化近似最佳逼近多项式为若有界证明对任何存在常数、使
1.设在上试将降低到3次多项式并估计误差.
1.在上利用幂级数项数求的3次逼近多项式使误差不超过
0.
005.
1.是上的连续奇偶函数证明不管是奇数或偶数的最佳逼近多项式也是奇偶函数.
1.求、使为最小.并与1题及6题的一次逼近多项式误差作比较.
1.、定义问它们是否构成内积
1.用许瓦兹不等式
4.5估计的上界并用积分中值定理估计同一积分的上下界并比较其结果.
1.选择使下列积分取得最小值:.
1.设空间分别在、上求出一个元素使得其为的最佳平方逼近并比较其结果.
1.在上求在上的最佳平方逼近.
1.是第二类切比雪夫多项式证明它有递推关系.
1.将在上按勒让德多项式及切比雪夫多项式展开求三次最佳平方逼近多项式并画出误差图形再计算均方误差.
1.把在上展成切比雪夫级数.
1.用最小二乘法求一个形如的经验公式使它与下列数据拟合并求均方误差.
192531384419.
032.
349.
073.
397.
81.观测物体的直线运动得出以下数据:时间秒
00.
91.
93.
03.
95.0距离米010305080110求运动方程.
1.在某化学反应里根据实验所得分解物的浓度与时间关系如下:时间0510152025303540455055浓度
01.
272.
162.
863.
443.
874.
154.
374.
514.
584.
624.64用最小二乘拟合求.
1.编出用正交多项式做最小二乘拟合的程序框图.
1.编出改进FFT算法的程序框图.
1.现给出一张记录试用改进FFT算法求出序列的离散频谱第四章数值积分与数值微分
1.确定下列求积公式中的待定参数使其代数精度尽量高并指明所构造出的求积公式所具有的代数精度:1;2;3;
4.
1.分别用梯形公式和辛普森公式计算下列积分:1;2;3;
4.
1.直接验证柯特斯公式
2.4具有5次代数精度.
1.用辛普森公式求积分并计算误差.
1.推导下列三种矩形求积公式:1;2;
3.
1.证明梯形公式
2.9和辛普森公式
2.11当时收敛到积分.
1.用复化梯形公式求积分问要将积分区间分成多少等分才能保证误差不超过设不计舍入误差
1.用龙贝格方法计算积分要求误差不超过.
1.卫星轨道是一个椭圆椭圆周长的计算公式是这里是椭圆的半长轴是地球中心与轨道中心椭圆中心的距离记为近地点距离为远地点距离公里为地球半径则.我国第一颗人造卫星近地点距离公里远地点距离公里试求卫星轨道的周长.
1.证明等式试依据的值用外推算法求的近似值.
1.用下列方法计算积分并比较结果.1龙贝格方法;1三点及五点高斯公式;1将积分区间分为四等分用复化两点高斯公式.
1.用三点公式和五点公式分别求在
1.
01.1和
1.2处的导数值并估计误差.的值由下表给出:
1.
01.
11.
21.
31.
40.
25000.
22680.
20660.
18900.1736第五章常微分方程数值解法
1.就初值问题分别导出尤拉方法和改进的尤拉方法的近似解的表达式,并与准确解相比较
2.用改进的尤拉方法解初值问题取步长h=
0.1计算,并与准确解相比较
3.用改进的尤拉方法解取步长h=
0.1计算,并与准确解相比较
4.用梯形方法解初值问题证明其近似解为并证明当时,它原初值问题的准确解
5.利用尤拉方法计算积分在点的近似值
6.取h=
0.2,用四阶经典的龙格-库塔方法求解下列初值问题1)2)
7.证明对任意参数t,下列龙格-库塔公式是二阶的
8.证明下列两种龙格-库塔方法是三阶的1)2)
9.分别用二阶显式亚当姆斯方法和二阶隐式亚当姆斯方法解下列初值问题取计算并与准确解相比较
10.证明解的下列差分公式是二阶的,并求出截断误差的首项
11.导出具有下列形式的三阶方法
12.将下列方程化为一阶方程组1)2)3)
13.取h=
0.25,用差分方法解边值问题
14.对方程可建立差分公式试用这一公式求解初值问题验证计算解恒等于准确解
15.取h=
0.2用差分方法解边值问题第六章方程求根
1.用二分法求方程的正根,要求误差
0.
052.用比例求根法求在区间
[01]内的一个根,直到近似根满足精度时终止计算
3.为求方程在附近的一个根,设将方程改写成下列等价形式,并建立相应的迭代公式1),迭代公式;2),迭代公式;3),迭代公式试分析每种迭代公式的收敛性,并选取一种公式求出具有四位有效数字的近似根
4.比较求的根到三位小数所需的计算量;1)在区间
[01]内用二分法;2用迭代法,取初值
5.给定函数,设对一切存在且,证明对于范围内的任意定数λ,迭代过程均收敛于的根
6.已知在区间[ab]内只有一根,而当axb时,,试问如何将化为适于迭代的形式?将化为适于迭代的形式,并求x=
4.5(弧度)附近的根
7.用下列方法求在附近的根根的准确值=
1.87938524…,要求计算结果准确到四位有效数字1用牛顿法;2)用弦截法,取;3)用抛物线法,取
8.用二分法和牛顿法求的最小正根
9.研究求的牛顿公式证明对一切且序列是递减的
10.对于的牛顿公式,证明收敛到,这里为的根
11.试就下列函数讨论牛顿法的收敛性和收敛速度
1212.应用牛顿法于方程,导出求立方根的迭代公式,并讨论其收敛性
13.应用牛顿法于方程,导出求的迭代公式,并用此公式求的值
14.应用牛顿法于方程和,分别导出求的迭代公式,并求
15.证明迭代公式是计算的三阶方法假定初值充分靠近根,求第七章解线性方程组的直接方法
1.考虑方程组a用高斯消去法解此方程组(用四位小数计算),a用列主元消去法解上述方程组并且与a比较结果
2.a设A是对称阵且,经过高斯消去法一步后,A约化为证明A2是对称矩阵b用高斯消去法解对称方程组
4.设A为n阶非奇异矩阵且有分解式A=LU,其中L为单位下三角阵,U为上三角阵,求证A的所有顺序主子式均不为零
5.由高斯消去法说明当时,则A=LU,其中L为单位下三角阵,U为上三角阵
6.设A为n阶矩阵,如果称A为对角优势阵证明若A是对角优势阵,经过高斯消去法一步后,A具有形式
7.设A是对称正定矩阵,经过高斯消去法一步后,A约化为,其中证明
(1)A的对角元素
(2)A2是对称正定矩阵;
(3)
(4)A的绝对值最大的元素必在对角线上;
(5)
(6)从
(2),
(3),
(5)推出,如果,则对所有k
8.设为指标为k的初等下三角阵,即(除第k列对角元下元素外,和单位阵I相同)求证当时,也是一个指标为k的初等下三角阵,其中为初等排列阵
9.试推导矩阵A的Crout分解A=LU的计算公式,其中L为下三角阵,U为单位上三角阵
10.设,其中U为三角矩阵a就U为上及下三角矩阵推导一般的求解公式,病写出算法b计算解三角形方程组的乘除法次数c设U为非奇异阵,试推导求的计算公式
11.证明(a)如果A是对称正定阵,则也是正定阵;(b)如果A是对称正定阵,则A可唯一写成其中L是具有正对角元的下三角阵
12.用高斯-约当方法求A的逆阵
13.用追赶法解三对角方程组,其中
14.用改进的平方根法解方程组
15.下述矩阵能否分解为LU(其中L为单位下三角阵,U为上三角阵)?若能分解,那么分解是否唯一?
16.试划出部分选主元素三角分解法框图,并且用此法解方程组.
17.如果方阵A有,则称A为带宽2t+1的带状矩阵,设A满足三角分解条件,试推导的计算公式,对1);2).
18.设,计算A的行范数,列范数,2-范数及F-范数
19.求证a,b
20.设且非奇异,又设为上一向量范数,定义试证明是上的一种向量范数
21.设为对称正定阵,定义,试证明为上向量的一种范数
22.设,求证
23.证明当且尽当x和y线性相关且时,才有
24.分别描述中(画图)
25.令是(或)上的任意一种范数,而P是任意非奇异实(或复)矩阵,定义范数,证明
26.设为上任意两种矩阵算子范数,证明存在常数,使对一切满足
27.设,求证与特征值相等,即求证
28.设A为非奇异矩阵,求证
29.设A为非奇异矩阵,且,求证存在且有估计
30.矩阵第一行乘以一数,成为证明当时,有最小值
31.设A为对称正定矩阵,且其分解为,其中,求证ab
32.设计算A的条件数
33.证明如果A是正交阵,则
34.设且为上矩阵的算子范数,证明第八章解方程组的迭代法
1.设方程组a考察用雅可比迭代法高斯-塞德尔迭代法解此方程组的收敛性;a用雅可比迭代法高斯-塞德尔迭代法解此方程组要求当时迭代终止.
2.设证明:即使级数也收敛.
3.证明对于任意选择的A序列收敛于零.4.设方程组迭代公式为求证:由上述迭代公式产生的向量序列收敛的充要条件是
5.设方程组ab试考察解此方程组的雅可比迭代法及高斯-塞德尔迭代法的收敛性
6.求证的充要条件是对任何向量x,都有
7.设,其中A对称正定,问解此方程组的雅可比迭代法是否一定收敛?试考察习题5a方程组
8.设方程组a求解此方程组的雅可比迭代法的迭代矩阵的谱半径;a求解此方程组的高斯-塞德尔迭代法的迭代矩阵的谱半径;a考察解此方程组的雅可比迭代法及高斯-塞德尔迭代法的收敛性
9.用SOR方法解方程组(分别取松弛因子)精确解要求当时迭代终止,并且对每一个值确定迭代次数
10.用SOR方法解方程组(取=
0.9)要求当时迭代终止
11.设有方程组,其中A为对称正定阵,迭代公式试证明当时上述迭代法收敛(其中)
12.用高斯-塞德尔方法解,用记的第i个分量,且a证明;a如果,其中是方程组的精确解,求证其中a设A是对称的,二次型证明a由此推出,如果A是具有正对角元素的非奇异矩阵,且高斯-塞德尔方法对任意初始向量是收敛的,则A是正定阵
13.设A与B为n阶矩阵,A为非奇异,考虑解方程组其中a找出下列迭代方法收敛的充要条件a找出下列迭代方法收敛的充要条件比较两个方法的收敛速度
14.证明矩阵对于是正定的,而雅可比迭代只对是收敛的
15.设,试说明A为可约矩阵
16.给定迭代过程,,其中,试证明如果C的特征值,则迭代过程最多迭代n次收敛于方程组的解
17.画出SOR迭代法的框图
18.设A为不可约弱对角优势阵且,求证解的SOR方法收敛
19.设,其中A为非奇异阵a求证为对称正定阵;b求证第九章矩阵的特征值与特征向量计算
1.用幂法计算下列矩阵的主特征值及对应的特征向量ab当特征值有3位小数稳定时迭代终止
2.方阵T分块形式为其中为方阵,T称为块上三角阵,如果对角块的阶数至多不超过2,则称T为准三角形形式,用记矩阵T的特征值集合,证明
3.利用反幂法求矩阵的最接近于6的特征值及对应的特征向量
4.求矩阵与特征值4对应的特征向量
5.用雅可比方法计算的全部特征值及特征向量,用此计算结果给出例3的关于p的最优值
6.a设A是对称矩阵,λ和是A的一个特征值及相应的特征向量,又设P为一个正交阵,使证明的第一行和第一列除了λ外其余元素均为零b对于矩阵λ=9是其特征值,是相应于9的特征向量,试求一初等反射阵P,使,并计算
7.利用初等反射阵将正交相似约化为对称三对角阵
8.设,且不全为零,为使的平面旋转阵,试推导计算第行,第j行元素公式及第i列,第j列元素的计算公式
9.设是由豪斯荷尔德方法得到的矩阵,又设y是的一个特征向量a证明矩阵A对应的特征向量是;b对于给出的y应如何计算x?
10.用带位移的QR方法计算ab全部特征值
11.试用初等反射阵A分解为QR,其中Q为正交阵,R为上三角阵,数值分析习题答案第一章绪论习题参考答案1.ε(lnx)≈1.1.有5位有效数字,有2位有效数字,有4位有效数字,有5位有效数字,有2位有效数字1.1.1.1.,1.1.1.,,故t增加时S的绝对误差增加,相对误差减小1.,计算过程不稳定1.,如果令,则,,,,,的结果最好1.,开平方时用六位函数表计算所得的误差为,分别代入等价公式中计算可得,1.方程组的真解为,而无论用方程一还是方程二代入消元均解得,结果十分可靠1.第二章插值法习题参考答案
1.;.
2..
3.线性插值取,则;二次插值取,则=-
0.
616707.
4.,其中.所以总误差界.
5.当时,取得最大值.
6.i对在处进行n次拉格朗日插值,则有由于,故有.ii构造函数在处进行n次拉格朗日插值,有.插值余项为,由于故有令即得.
7.以ab两点为插值节点作的一次插值多项式,据余项定理,,由于故
8.截断误差其中则时取得最大值.由题意,所以,
9.则可得,,则可得
10.数学归纳法证当时,为m-1次多项式;假设是m-k次多项式,设为,则为m-k+1次多项式,得证
11.右左
12.
13..
14.由于是的n个互异的零点,所以对求导得,则,记则由以上两式得
15.i.ii证明同上
16.
17.即均为的二重零点因而有形式作辅助函数则由罗尔定理,存在使得类似再用三次罗尔定理,存在使得又可得即
18.采用牛顿插值,作均差表一阶均差二阶均差01201110-1/2又由得所以
19.记则因为,所以在上一致连续当时,,此时有由定义知当时,在上一致收敛于
20.在每个小区间上表示为计算各值的C程序如下#includestdio.h#includemath.hfloatffloatx{return1/1+x*x;}floatIfloatxfloatafloatb{returnx-b/a-b*fa+x-a/b-a*fb;}voidmain{inti;floatx
[11]xcxx;x
[0]=-5;printfx
[0]=%f\nx
[0];fori=1;i=10;i++{x[i]=x[i-1]+1;printfx[%d]=%f\nix[i];}fori=0;i10;i++{xc=x[i]+x[i+1]/2;Ixcx[i]x[i+1];printfI[%d]=%f\ni+1Ixcx[i]x[i+1];}fori=0;i10;i++{xx=x[i]+x[i+1]/2;fxx;printff[%d]=%f\ni+1fxx;}}
21.在每个小区间上为
22.则在每个小区间上表示为
23.则三次样条插值函数表达式为i由,得,关于的方程组为
24.i因为所以右==左ii由于为三次函数,故为常数,又,则,所以第三章函数逼近与计算习题参考答案
1.a区间变换公式为,代入原公式可得新区间里的伯恩斯坦多项式为;b,相应的麦克劳林级数分别为,部分和误差则为,,大于伯恩斯坦多项式的误差
1.,故,当时,
1.,对任意不超过6次的多项式,在时,若有,则在上至少有7个零点,这与不超过6次矛盾,所以,就是所求最佳一致逼近多项式
1.设所求为,,由47页定理4可知在上至少有两个正负交错的偏差点,恰好分别为的最大值和最小值处,故由可以解得即为所求
1.原函数与零的偏差极大值点分别为,故,解方程可得出唯一解
1.,故,得,,故所求最佳一次逼近多项式为,又因为两个偏差点必在区间端点,故误差限为
1.,故由可以解得,,则有,故所求最佳一次逼近多项式为
1.切比雪夫多项式在上对零偏差最小,所求函数必为切比雪夫多项式的常数倍,,解得唯一解
1.作变换代入得,则在上的三次最佳逼近多项式为,作逆变换代入,则在上的三次最佳逼近多项式为
1.,,,,其中
1.,故正交
1.用的4个零点做插值节点可求得三次近似最佳逼近多项式为
1.,则有,其中由拉格朗日插值的余项表达公式可得出,令,则待证不等式成立,得证
1.由泰勒级数项数节约,在上有,即其中误差限为
1.,取为的近似,误差限为,再对幂级数的项数进行节约就可以得到原函数的3次逼近多项式,其误差限为,即为所求
1.当为上的奇函数时,设为原函数的最佳逼近多项式,则,对有,所以也是最佳逼近多项式,由最佳逼近多项式的唯一性,,即是奇函数同理可证,当为上的偶函数时,最佳逼近多项式也是偶函数
1.,为使均方误差最小,则有,解得
1.a,,c为常数,,,但当时,不满足定义,所以不构成内积(b),,,且当且仅当时,满足定义,所以构成内积
1.,,其中,则,由此可知用积分中值定理估计比许瓦兹不等式估计更精确
1.,时最小在时,值为,时,值为1,时,值为,时最小
1.要使最小,由拉格朗日乘子法可解得,误差为,要使最小,由拉格朗日乘子法可解得,误差为,前者误差小
1.上均为偶函数,也为偶函数,则最小,由拉格朗日乘子法可解得
1.,和差化积得证
1.由积分区间的对称性及勒让德多项式的奇偶性可知,,将原函数在此积分区间上按勒让德多项式三次展开就可以求得,,代入可得,均方误差为
1.,其中
1.,,解方程得,均方误差
1.经验公式为,最小二乘法解得,运动方程为
1.经验公式为,最小二乘法解得,浓度与时间的函数关系为
1.
1.
1.,,,,,,,,,第四章数值积分与数值微分习题参考答案
1.1公式可对均准确成立,即解得,具有3次代数精度2,具有3次代数精度3或具有2次代数精度4,具有3次代数精度
2.1==
0.1114=
0.
11162343.柯特斯公式为.其中.验证对于,均成立,但时不成立
4.=
0.63233,所以
5.1此差值型求积公式的余项为由于在上恒为正,故在上存在一点,使所以有
236.梯形公式和辛甫森公式的余项分别为其中,所以当时,,即两公式均收敛到积分,且分别为二阶和四阶收敛
7.设将积分区间分成n等分则应有其中,解得
8.首先算出,然后逐次应用3个加速公式计算结果如下表k
01230.
683940.
645240.
635410.
632940.
632340.
632130.
632120.
632130.
632120.63212所以,积分
9.,,所以=4×
7782.5×
1.56529=48728(可任选一种数值积分方法,如柯特斯公式)
10.由泰勒展开式有由于,用外推算法,令,则,,,即的近似值为
3.
1415911.1计算结果如下表k
01231.
333331.
166671.
116671.
103211.
111121.
100001.
098721.
099261.
098631.09862即积分I=
1.098622,令三点高斯公式五点高斯公式=
1.098623对每个积分用高斯公式,得I=
1.09854此积分精确值为
12.三点公式,,的误差的误差的误差五点公式误差分别为,,第五章常微分方程数值解法习题参考答案1.尤拉法表达式,误差,改进尤拉法表达式,无误差1.近似解准确解近似解准确解
0.
11.
111.
110340.
62.
040862.
044240.
21.
242051.
242810.
72.
323152.
327510.
31.
398471.
399720.
82.
645582.
651080.
41.
581811.
583650.
93.
012373.
019210.
51.
794901.
797441.
03.
428173.436561.近似解准确解
0.
10.
00550.
005162580.
20.
02192750.
02126920.
30.
05014440.
04918180.
40.
09093070.
08968000.
50.
1449920.1434691.,即,又由,则有当时,1.取步长h=
0.5,,f
0.5=
0.500000,f1=
1.14201,f
1.5=
2.50115,f2=
7.245021.1近似解2近似解
0.
21.
242800.
21.
727630.
41.
583640.
42.
743020.
62.
044210.
64.
094240.
82.
651030.
85.
829271.
03.
436551.
07.996011.,,,则1.1令,泰勒展开可得,,,同理有,代入龙格-库塔公式可得2类似1展开可得,,,同理有,代入龙格-库塔公式可得1.二阶显式公式为,代入得,二阶隐式公式为,代入得,真解为1.,,,,代入得,截断误差首项为1.,,,,,代入待定系数的公式中可得系数之间的关系式为,,,1.1,其中2,其中3,其中1.用差商逼近导数的方法把原边值问题转化为等价差分法方程组可得,解此方程组可得1.,初值条件等于准确解,由数学归纳法代入差分公式中可得,即差分法求出的解恒等于准确解1.差分方程,,代入得,第六章方程求根
1.令,则符号0021-
1121.5-
21.
521.75+
31.
51.
751.625+
41.
51.
6251.5625-
51.
56251.
6251.5938-
2.
3.1,在附近,,迭代公式收敛2,在附近,迭代公式收敛,迭代得近似值
1.4663,,,迭代公式发散
4.1二分14次得
0.0905456;2迭代5次得
0.
09052465.迭代函数,,由已知,有,所以即迭代过程收敛
6.将转化为,此时在附近,,所以迭代格式为,迭代三次得
4.
49347.1牛顿法迭代格式,迭代三次得
1.8792弦截法迭代格式,迭代三次得
1.8793抛物线法,故,则,迭代三次得
1.
8798.最小正根为
4.
49349.,即,即,序列单调递减
10.迭代函数为,且有,,.其中介于与之间将上式两边除以,并将处泰勒展开得,其中介于与之间将上式两边取极限,及,,得
11.1,迭代格式发散2,迭代格式收敛,且收敛到要使,则,为一阶收敛
12.令,迭代公式为,则,所以,又,所以,因此迭代格式为线性收敛
13.,取,迭代三次得
14.求的迭代公式分别为,设迭代函数为,则,.
15.记迭代函数,则,由上
①两边求导得则可得对
①式两边求二阶导数得则可得对
①式两边求二阶导数得则可得所以迭代公式是三阶方法,且.第七章解线性方程组的直接方法习题参考答案1.a高斯消去法解得;b列主元消去法解得1.a,故对称b高斯消去法解得1.a,b由及a的结论可得1.因为非奇异,的对角元不为零,又分解等价于高斯消去法,,由引理可知,矩阵的顺序主子式均不为零1.高斯消去法第步等价于左乘单位下三角矩阵,而顺序主子式均不为零保证所得矩阵对角元不为零,可进行第步消元,,1.,则是对角优势阵,故高斯消去法与部分选主元高斯消去法对于对称的对角优势阵每一步均选取同样的主元,得出的是同样的结果1.1,2,又有当时,故是对称正定矩阵,3,4若,令,由于和也是对称正定矩阵,代入得,矛盾,故的绝对值最大的元素必在对角线上,5,6对所有均有对称正定,1.,其中与位置互换1.对施行初等列变换,,进行次初等列变换后,令即为所求1.a若为阶可逆下三角矩阵,,则当时,而当时,,算法即从第一行开始顺序循环,同理可知若为阶可逆上三角矩阵,则当时,而当时,,算法即从最后一行开始逆序循环,b第k步循环进行k次乘除法,共进行次乘除法,c1.a,由此可知也是对称矩阵,,由此可知也是对称正定矩阵,b,得出唯一正对角元的下三角阵使得1.1.1.1.按高斯消去法,无法进行第二次消去,换行后可以分解,第二次消去可乘任意系数,分解不唯一,可唯一分解1.,解得1.高斯消去法公式中去掉即可推出该公式1.1.ab1.,,,故是上的向量范数1.,故,故是上的向量范数1.1.充分性若有和线性相关且,即,代入得;唯一性若有,由于,两边同时平方可得出,消去共同项可得,当且仅当和线性相关时等号成立1.以上图像分别为,,1.1.由向量范数的相容性可知存在常数,使得,于是令0,0,则对任意,均有不等式1.若,则就有,可推出即,同理可以推出,综合这两点即可得1.1.,则,故存在,1.,当时,,当时,,当时,有最小值71.a,b,1.,1.1.第八章解线性方程组的迭代法习题参考答案
1.aJacobi迭代矩阵特征方程为特征根均小于1,Jacobi迭代法收敛Gauss-Seidel迭代矩阵特征方程为特征根均小于1,Gauss-Seidel迭代法收敛bJacobi迭代格式为其中B如上,,迭代18次得,Gauss-Seidel迭代格式为其中G如上,,迭代8次得
2.证,则故,因此,即级数收敛
3.证设,一方面,,另一方面,因此,即序列收敛于零
4.证由已知迭代公式得迭代矩阵则特征多项式为解得,向量序列收敛的充要条件是,即
5.a谱半径,Jacobi迭代法不收敛;矩阵A对称正定,故Gauss-Seidel迭代法收敛b谱半径,Jacobi迭代法收敛;谱半径,Gauss-Seidel迭代法不收敛;
6.证必要性,则,对任意向量,有因而有,即充分性因对任何向量,都有,令,则即当时,的任一列向量的极限为A的对应的列向量,因而有
7.A对称正定,Jacobi迭代法不一定收敛,如题5a
7.aJacobi迭代矩阵的谱半径;bGauss-Seidel迭代矩阵的谱半径;c两种方法的谱半径均小于1,所以两种方法均收敛事实上,对于方程组,矩阵A为严格对角占优则Jacobi和Gauss-Seidel迭代法均收敛
7.取,迭代公式为使当时迭代终止,取时,迭代5次达到;取时,迭代6次达到;取时,迭代6次达到
7.迭代公式为取,,迭代8次达到精度要求
7.证所给迭代公式的迭代矩阵为,其n个特征值分别为,当时,有,因而,迭代法收敛
7.证a即为Gauss-Seidel迭代格式b由及,可得;其中,cd
13.a由已知,有,及,则,即由到的迭代矩阵为,所以由到的迭代矩阵为,则迭代方法收敛的充要条件为b由已知可推得,所以迭代矩阵为,则迭代方法收敛的充要条件为由迭代矩阵可以看出,b迭代法的收敛速度是a的2倍
14.证由于,当时,,所以A正定Jacobi迭代矩阵谱半径为,所以只对收敛
15.取排列阵,则A为可约矩阵
16.证迭代矩阵的特征方程为,若,则,所以,即对任给向量,迭代n次后,,其中,则即最多迭代n次收敛于方程组的解
17.用SOR方法解方程组,其中A对称正定,数组x用来存放解向量,用控制迭代终止,k表示迭代次数
18.证方程组的SOR迭代矩阵为,特征方程,即,记只要当时,,则的根均满足A不可约则G也不可约,又A为弱对角优势阵,则当且时,即时,G为不可约弱对角占优,于是有,故,SOR方法收敛
19.证a,设,则,为对称正定阵b因为为对称阵,所以左右左
20.证A为严格对角占优,则存在第九章矩阵的特征值与特征向量计算习题参考答案1.a取初始值111得
11.
0000000.
7500000.
0000008.
00000031.
0000000.617564-
0.
3711059.
54054051.
0000000.606413-
0.
3930959.
60407471.
0000000.605660-
0.
3941459.
60392181.
0000000.605583-
0.
3943779.605270b取初始值111得
10.
2857140.
7142861.
0000007.0000005-
0.
4893031.
0000000.
2497438.5614479-
0.
5947201.
0000000.
1591198.85623713-
0.
6035811.
0000000.
1514908.86782914-
0.
6039461.
0000000.
1511768.8687942.,使得,即,一定存在使得,则,,反之,故3.,由幂法得,原矩阵最接近6的特征值为,对应的特征向量为4.设特征向量为,则有,解得对应的特征向量为5.雅可比迭代进行五步可得,对应的特征向量分别为,,,最优值6.a,正交,则第一列,,又是对称矩阵,的第一行和第一列除外均为零b为反射阵,,解得,7.由豪斯荷尔德方法得,8.,解得,代入得
9.a,b由可求出初等反射阵,依次类推,10.a令,带位移QR方法计算可得,b令,带位移QR方法计算可得11.,故有输入初始节点,权函数及正交多项式次数n,计算判断计算是否令输入初始数组,等分点数,计算判断计算,,,计算,,,令否是判断判断否是否是k=0i=1i=n|P||P0|P0=P;|P0|ε输出xk;是否否是。