数值分析习题集和答案解析
数值分析课后习题及答案

第一章 绪论(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 ,并估计插值余项。
数值分析习题和答案解析(最新整理)

(1)
要使
应满足().
(2) 已知方程组
,则解此方程组的
Jacobi 迭代法是否收敛().它的渐近收敛速度 R(B)=
公式(6.13)直接计算即可。
对
,取 n=8,在分点处计算 f(x)的值构造函数表。
按式(6.11)求出
,按式(6.13)求得
,
积分
2. 用 Simpson 公式求积分 ,并估计误差 解:直接用 Simpson 公式(6.7)得
由(6.8)式估计误差,因
,故
3. 确定下列求积公式中的待定参数,使其代数精确度尽量 高,并指明求积公式所具有的代数精确度.
专业知识分享
WORD 格式
编辑整理
11. 填空题
(1) 满 足 条 件
的插值多项式
p(x)=( ).
(2)
,则 f[1,2,3,4]=( ),f[1,2,3,4,5]
=( ).
(3) 设
为互异节点, 为对应的四次插值基函
数,则
=( ),
=( ).
(4) 设
是区间[0,1]上权函数为 ρ(x)=x 的最
WORD 格式
编辑整理
误差估计由公式(5.19)得
这里 仍为 0.565 8. 求 一 个 次 数 不 高 于 四 次 的 多 项 式 p(x),使 它 满 足
解:这种题目可以有很多方法去做,但应以简单为宜。此处
可先造 使它满足
,显然 p(x)=x2(2-x)+Ax2(x-1)2
,再令
由 p(2)=1 求出 A= ,于是
5.计算
取 ,利用 :
式计算误差最小。
四个选项: 第二、三章 插值与函数逼近
习题二、三
数值分析题库答案(含详细解题步骤)

第 1 页/共 22 页1. 正方形的边长大约为100cm ,应怎样测量才干使面积误差不超过1cm 22. 已测得某场地长l 的值为110=*l m ,宽d 的值为80=*d m ,已知 2.0≤-*l l m,1.0≤-*d d m, 试求面积ld s =的绝对误差限与相对误差限.3.为使π的相对误差小于0.001%,至少应取几位有效数字?4.设x的相对误差界为δ,求n x的相对误差界.5.设有3个近似数a=2.31,b=1.93,c=2.24,它们都有3位有效数字,试计算p=a+bc的误差界和相对误差界,并问p的计算结果能有几位有效数字?第 3 页/共 22 页6. 已知333487.034.0sin ,314567.032.0sin ==,请用线性插值计算3367.0sin 的值,并预计截断误差.7. 已知sin0.32=0.314567, sin0.34=0.333487, sin0.36= 0.352274,用抛物插值计算sin0.3367的值, 并预计误差.8. 已知16243sin ,sin πππ===请用抛物插值求sin50的值,并预计误差9. . .6,8,7,4,1)(,5,4,3,2,1求四次牛顿插值多项式时设当==i i x f x第 5 页/共 22 页10. 已知4)2(,3)1(,0)1(=-=-=f f f , 求函数)(x f 过这3点的2次牛顿插 值多项式.11. 设x x f =)(,并已知483240.1)2.2(,449138.1)1.2(,414214.1)0.2(===f f f ,试用二次牛顿插值多项式计算(2.15)f 的近似值,并研究其误差12. 设],[)(b a x f 在上有四阶延续导数,试求满意条件)2,1,0()()(==i x f x P i i 及)()(11x f x P '='的插值多项式及其余项表达式.13. 给定3201219(),,1,,44f x x x x x ====试求()f x 在1944⎡⎤⎢⎥⎣⎦,上的三次埃尔米特插值多项式()P x ,使它满意11()()(0,1,2),()(),i i P x f x i P x f x ''===并写出余项第 7 页/共 22 页表达式.14. 设],1,0[,23)(2∈++=x x x x f 试求)(x f 在]1,0[上关于,,1{,1)(x span x =Φ=ρ}2x 的最佳平方逼近多项式15.已知实验数据如下:用最小二乘法求形如y=a+bx2的拟合曲线,并计算均方误差.16.已知数据表如下第 9 页/共 22 页x i 1 2 3 4 5 y iωi4 4.56 8 8.5 2 1 3 1 1试用最小二乘法求多项式曲线与此数据组拟合17. .1)(},1{span ,1]41[)(的最佳平方逼近多项式中的关于上的在在求==Φ=x x x x f ρ18. 决定求积公式⎰++≈10110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A , 使其代数精度尽量高,并指出所决定的求积公式的代数精度.19. 用复化辛普森公式计算积分⎰=10dx e I x , 问区间[0,1]应分多少等分才干使截断误差不超过?10215-⨯第 11 页/共 22 页20. 利用下表中给出的数据,分离用复化梯形公式和复化辛甫生公式计算定积分dx x I ln 21⎰=的近似值(要求结果保留到小数点后六位)21. 用复化梯形公式和复化辛甫生公式计算积分⎰=6.28.1)(dx x f I ,函数)(x f 在某些节点上的值如下图:(本题共14分)22. 决定公式⎰+≈101100)()()(x f A x f A dx x f x 的系数1010,,,x x A A ,使其具有最高代数精度23. 决定求积公式⎰++≈1110)1()(32)0()(f A x f f A dx x f 中的待定参数110,,A x A ,使其代数精度尽量高,并指出所决定的求积公式的代数精度第 13 页/共 22 页24.用LU 分解法求解以下方程组 (10分)123123142521831520x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭25.用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎪⎪⎪⎭⎫ ⎝⎛8892121514131615141321x x x26. 用LU 分解法求解以下方程组⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛⎪⎪⎪⎭⎫⎝⎛542631531321321x x x27. 设方程组b Ax =,其中⎪⎪⎪⎭⎫⎝⎛-=220122101A ,Tb ⎪⎭⎫ ⎝⎛-=32,31,21, 已知它有解Tx ⎪⎭⎫⎝⎛-=0,31,21,若右端有小扰动61021-∞⨯=bδ,试预计由此引起的解的相对误差.第 15 页/共 22 页28. 设方程组b Ax =,其中212 1.0001A -⎛⎫= ⎪-⎝⎭,11.0001b -⎛⎫= ⎪⎝⎭,当右端向量b 有误差00.0001δ⎛⎫= ⎪⎝⎭b 时,试预计由此引起的解的相对误差(用∞范数计算)29. 给定b Ax =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=111a a a a a a A 证实:(1) 当121<<-a 时,A 对称正定,从而GS 法收敛. (2) 惟独当2121<<-a 时,J 法收敛.30. 对于线性方程组⎪⎩⎪⎨⎧-=+-=-+=+1242043 16343232121x x x x x x x ,列出求解此方程组的Jacobi 迭代格式,并判断是否收敛。
数值分析练习题加答案(一)

数值分析期末考试一、 设80~=x ,若要确保其近似数的相对误差限为0.1%,则它的近似数x 至少取几位有效数字?(4分)解:设x 有n 位有效数字。
因为98180648=<<=,所以可得x 的第一位有效数字为8(1分) 又因为21101011000110821--⨯=<⨯⨯≤n ε,令321=⇒-=-n n ,可知x 至少具有3位有效数字(3分)。
二、求矩阵A 的条件数1)(A Cond (4分)。
其中⎥⎦⎤⎢⎣⎡=4231A 解:⎥⎦⎤⎢⎣⎡--=-5.05.1121A (1分) 1A =7(1分) 2711=-A (1分)249)(1=A Cond (1分)三、用列主元Gauss 消元法法求解以下方程组(6分)942822032321321321=++-=++--=+-x x x x x x x x x解:→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----5.245.2405.35.230914220321821191429142821120321 ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---→⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---8175835005,245.24091425.33.2305.245.2409142(4分) 等价三角方程组为:⎪⎪⎩⎪⎪⎨⎧-=-=+-=++,8175835,5.245.24,942332321x x x x x x (1分)回代得1,3,5123==-=x x x (1分)四、设.0,2,3,1,103)(3210234=-===-+-=x x x x x x x x f 1)求以3210,,,x x x x 为节3次Lagrange 多项式;(6分) 2)求以3210,,,x x x x 为节3次Newton 多项式;(6分)3)给出以上插值多项式的插值余项的表达式(3分)解:由0,2,3,13210=-===x x x x 可得10)(,34)(,1)(,11)(3210-==-=-=x f x f x f x f即得: +------+------=))()(())()(()())()(())()(()()(312101320130201032103x x x x x x x x x x x x x f x x x x x x x x x x x x x f x L=------+------))()(())()(()())()(())()(()(23130321033212023102x x x x x x x x x x x x x f x x x x x x x x x x x x x f+-+--+-⨯-+-+--+-⨯-)03)(23)(13()0)(2)(1()1()01)(21)(31()0)(2)(3(11x x x x x x326610.)20)(30)(10()2)(3)(1()10()02)(32)(12()0)(3)(1(34x x x x x x x x x -+--=+--+--⨯-+---------⨯2)计算差商表如下:i x )(i x f 一阶差商 二阶差商 三阶差商1 -11 3 -1 5 -2 34 -7 4 0-10-225-1则=+-----+-+-=)2)(3)(1()3)(1(4)1(511)(3x x x x x x x N326610x x x -+--3))2)(3)(1())()()((!4)()(3210)4(3+--=----=x x x x x x x x x x x x f x R ξ五、给定方程组b Ax =,其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=100131w w w w A 。
《数值分析》练习题及答案解析

《数值分析》练习题及答案解析一、单选题1. 以下误差公式不正确的是( D )A .()1212x x x x ∆-≈∆-∆B .()1212x x x x ∆+≈∆+∆C .()122112x x x x x x ∆≈∆+∆D .1122()x x x x ∆≈∆-∆ 2. 已知等距节点的插值型求积公式()()352kkk f x dx A f x =≈∑⎰,那么3kk A==∑( C )A .1 B. 2 C.3 D. 4 3.辛卜生公式的余项为( c )A .()()32880b a f η-''-B .()()312b a f η-''-C .()()()542880b a f η--D .()()()452880b a f η--4. 用紧凑格式对矩阵4222222312A -⎡⎤⎢⎥=-⎢⎥⎢⎥--⎣⎦进行的三角分解,则22r =( A ) A .1 B .12C .–1D .–25. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足( D ) A .()00l x =0,()110l x = B . ()00l x =0,()111l x = C .()00l x =1,()111l x = D . ()00l x =1,()111l x =6. 用二分法求方程()0f x =在区间[],a b 上的根,若给定误差限ε,则计算二分次数的公式是n ≥( D )A .ln()ln 1ln 2b a ε-++ B. ln()ln 1ln 2b a ε-+-C.ln()ln 1ln 2b a ε--+ D. ln()ln 1ln 2b a ε--- 7.若用列主元消去法求解下列线性方程组,其主元必定在系数矩阵主对角线上的方程组是( B )A .123123123104025261x x x x x x x x x -+=⎧⎪-+=⎨⎪-+=-⎩ B 。
数值分析(课后习题答案详解).ppt

x x 41 2 0 . 25 0 . 5451 1 1 再解 3 x 0 . 875 ,得 x 1 . 2916 2 2 2 0 3 1 . 7083 . 5694 x x 3 3
4 41 2 T 故得 GG 分解: A 1 2 3 2 2 3 3 3 1 1 16 11 4 2 T 3 1 LDL 分解为: A 1 4 4 1 2 3 1 1 9 1 2 2
一.习题1(第10页)
1-1.下列各数都是经过四舍五入得到的近似值 ,试分 别指出它们的绝对误差限,相对误差限和有效数字的位数.
x1=5.420,x2=0.5420,x3=0.00542,x4=6000,x5=0.6105.
解 绝对误差限分别为: 1=0.510-3,2=0.510-4, 3=0.510-5,4=0.5,5=0.5104 . 相对误差限分别为: r1=0.510-3/5.420=0.00923%, r2=0.00923%,r3=0.0923%,4=0.0083%,5=8.3%. 有效数位分别为: 4位,4位,3位,4位,1位. 1-2.下列近似值的绝对误差限都是0.005,试问它们有
2 11 2 1 2 故得 Crout 分解: A 4 3 13 6 12 1 1
1 2 11 2 1 2 LDM 分解为: A 21 13 3 3 4 1 1 1
几位有效数字. a=-1.00031,b=0.042,c=-0.00032
《数值分析》练习题及答案解析

《数值分析》练习题及答案解析第一章 绪论主要考查点:有效数字,相对误差、绝对误差定义及关系;误差分类;误差控制的基本原则;。
1. 3.142和3.141分别作为π的近似数具有( )和( )位有效数字.A .4和3B .3和2C .3和4D .4和4 答案:A2. 设 2.3149541...x *=,取5位有效数字,则所得的近似值x=___________ .答案:2.31503.若近似数2*103400.0-⨯=x 的绝对误差限为5105.0-⨯,那么近似数有几位有效数字 解:2*103400.0-⨯=x ,325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
4 . 14159.3=π具有4位有效数字的近似值是多少?解:10314159.0⨯= π,欲使其近似值*π具有4位有效数字,必需!41*1021-⨯≤-ππ,3*310211021--⨯+≤≤⨯-πππ,即14209.314109.3*≤≤π即取( , )之间的任意数,都具有4位有效数字。
第二章 非线性方程求根 主要考查点:二分法N 步后根所在的区间,及给定精度下二分的次数计算;非线性方程一般迭代格式的构造,(局部)收敛性的判断,迭代次数计算; 牛顿迭代格式构造;求收敛阶;1.用二分法求方程012=--x x 的正根,要求误差小于0.05。
(二分法)解:1)(2--=x x x f ,01)0(<-=f ,01)2(>=f ,)(x f 在[0,2]连续,故[0,2]为函数的有根区间。
"(1)计算01)1(<-=f ,故有根区间为[1,2]。
(2)计算041123)23()23(2<-=--=f ,故有根区间为]2,23[。
(3)计算0165147)47()47(2>=--=f ,故有根区间为]47,23[。
(4)计算06411813)813()813(2>=--=f ,故有根区间为]813,23[。
数值分析课后习题答案

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
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第一章绪论习题一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),使它满足解:这种题目可以有很多方法去做,但应以简单为宜。
此处可先造使它满足,显然,再令p(x)=x2(2-x)+Ax2(x-1)2由p(2)=1求出A=,于是9. 令称为第二类Chebyshev多项式,试求的表达式,并证明是[-1,1]上带权的正交多项式序列。
解:因10. 用最小二乘法求一个形如的经验公式,使它拟合下列数据,并计算均方误差.解:本题给出拟合曲线,即,故法方程系数法方程为解得最小二乘拟合曲线为均方程为11. 填空题(1) 满足条件的插值多项式p(x)=( ).(2) ,则f[1,2,3,4]=( ),f[1,2,3,4,5]=( ).(3) 设为互异节点,为对应的四次插值基函数,则=( ),=( ).(4) 设是区间[0,1]上权函数为ρ(x)=x的最高项系数为1的正交多项式序列,其中,则=( ),=( )答:(1)(2)(3)(4)第4章数值积分与数值微分习题41. 分别用复合梯形公式及复合Simpson公式计算下列积分.解本题只要根据复合梯形公式(6.11)及复合Simpson 公式(6.13)直接计算即可。
对,取n=8,在分点处计算f(x)的值构造函数表。
按式(6.11)求出,按式(6.13)求得,积分2. 用Simpson公式求积分,并估计误差解:直接用Simpson公式(6.7)得由(6.8)式估计误差,因,故3. 确定下列求积公式中的待定参数,使其代数精确度尽量高,并指明求积公式所具有的代数精确度.(1)(2)(3)解:本题直接利用求积公式精确度定义,则可突出求积公式的参数。
(1)令代入公式两端并使其相等,得解此方程组得,于是有再令,得故求积公式具有3次代数精确度。
(2)令代入公式两端使其相等,得解出得而对不准确成立,故求积公式具有3次代数精确度。
(3)令代入公式精确成立,得解得,得求积公式对故求积公式具有2次代数精确度。
4. 计算积分,若用复合Simpson公式要使误差不超过,问区间要分为多少等分?若改用复合梯形公式达到同样精确度,区间应分为多少等分?解:由Simpson公式余项及得即,取n=6,即区间分为12等分可使误差不超过对梯形公式同样,由余项公式得即取n=255才更使复合梯形公式误差不超过5. 用Romberg求积算法求积分,取解:本题只要对积分使用Romberg算法(6.20),计算到K=3,结果如下表所示。
于是积分,积分准确值为0.713272 6.用三点Gauss-Legendre求积公式计算积分.解:本题直接应用三点Gauss公式计算即可。
由于区间为,所以先做变换于是本题精确值7.用三点Gauss-Chebyshev求积公式计算积分解:本题直接用Gauss-Chebyshev求积公式计算即于是,因n=2,即为三点公式,于是,即故8. 试确定常数A,B,C,及α,使求积公式有尽可能高的代数精确度,并指出所得求积公式的代数精确度是多少.它是否为Gauss型的求积公式?解:本题仍可根据代数精确度定义确定参数满足的方程,令对公式精确成立,得到由(2)(4)得A=C,这两个方程不独立。
故可令,得(5)由(3)(5)解得,代入(1)得则有求积公式令公式精确成立,故求积公式具有5次代数精确度。
三点求积公式最高代数精确度为5次,故它是Gauss型的。
第五章解线性方程组的直接法习题五1. 用Gauss消去法求解下列方程组.解本题是Gauss消去法解具体方程组,只要直接用消元公式及回代公式直接计算即可。
故2. 用列主元消去法求解方程组并求出系数矩阵A的行列式detA的值解:先选列主元,2行与1行交换得消元3行与2行交换消元回代得解行列式得3. 用Doolittle分解法求的解.解:由矩阵乘法得再由求得由解得4. 下述矩阵能否作Doolittle分解,若能分解,分解式是否唯一?解:A中,若A能分解,一步分解后,,相互矛盾,故A不能分解,但,若A中1行与2行交换,则可分解为LU对B,显然,但它仍可分解为分解不唯一,为一任意常数,且U奇异。
C可分解,且唯一。
5. 用追赶法解三对角方程组Ax=b,其中解:用解对三角方程组的追赶法公式(3.1.2)和(3.1.3)计算得6. 用平方根法解方程组解:用分解直接算得由及求得7. 设,证明解:即,另一方面故8.设计算A的行范数,列范数及F-范数和2范数解:故9.设为上任一种范数,是非奇异的,定义,证明证明:根据矩阵算子定义和定义,得令,因P非奇异,故x与y为一对一,于是10. 求下面两个方程组的解,并利用矩阵的条件数估计.,即,即解:记则的解,而的解故而由(3.12)的误差估计得表明估计略大,是符合实际的。
11.是非题(若"是"在末尾()填+,"不是"填-):题目中(1)若A对称正定,,则是上的一种向量范数()(2)定义是一种范数矩阵()(3)定义是一种范数矩阵()(4)只要,则A总可分解为A=LU,其中L为单位下三角阵,U为非奇上三角阵()(5)只要,则总可用列主元消去法求得方程组的解()(6)若A对称正定,则A可分解为,其中L为对角元素为正的下三角阵()(7)对任何都有()(8)若A为正交矩阵,则()答案:(1)(+)(2)(-)(3)(+)(4)(-)(5)(+)(6)(+)(7)(-)(8)(+)第六章解线性方程组的迭代法习题六1.证明对于任意的矩阵A,序列收敛于零矩阵解:由于而故2. 方程组(1) 考查用Jacobi法和GS法解此方程组的收敛性.(2) 写出用J法及GS法解此方程组的迭代公式并以计算到为止解:因为具有严格对角占优,故J法与GS法均收敛。
(2)J法得迭代公式是取,迭代到18次有GS迭代法计算公式为取3. 设方程组证明解此方程的Jacobi迭代法与Gauss-Seidel迭代法同时收敛或发散解:Jacobi迭代为其迭代矩阵,谱半径为,而Gauss-Seide迭代法为其迭代矩阵,其谱半径为由于,故Jacobi迭代法与Gauss-Seidel法同时收敛或同时发散。
4. 下列两个方程组Ax=b,若分别用J法及GS法求解,是否收敛?解:Jacobi法的迭代矩阵是即,故,J法收敛、GS法的迭代矩阵为故,解此方程组的GS法不收敛。
5. 设,detA≠0,用,b表示解方程组Ax=f的J法及GS法收敛的充分必要条件.解J法迭代矩阵为,故J法收敛的充要条件是。
GS法迭代矩阵为由得GS法收敛得充要条件是6. 用SOR方法解方程组(分别取ω=1.03,ω=1,ω=1.1)精确解,要求当时迭代终止,并对每一个ω值确定迭代次数解:用SOR方法解此方程组的迭代公式为取,当时,迭代5次达到要求若取,迭代6次得7. 对上题求出SOR迭代法的最优松弛因子及渐近收敛速度,并求J法与GS法的渐近收敛速度.若要使那么J法GS法和SOR法各需迭代多少次?解:J法的迭代矩阵为,故,因A为对称正定三对角阵,最优松弛因子J法收敛速度由于,故若要求,于是迭代次数对于J法,取K=15对于GS法,取K=8对于SOR法,取K=58. 填空题(1)要使应满足().(2) 已知方程组,则解此方程组的Jacobi迭代法是否收敛().它的渐近收敛速度R(B)=().(3) 设方程组Ax=b,其中其J法的迭代矩阵是().GS法的迭代矩阵是().(4) 用GS法解方程组,其中a为实数,方法收敛的充要条件是a满足().(5) 给定方程组,a为实数.当a满足(),且0<ω<2时SOR迭代法收敛.答:(1)(2)J法是收敛的,(3)J法迭代矩阵是,GS法迭代矩阵(4)满足(5)满足第七章非线性方程求根习题七1.用二分法求方程的正根,使误差小于0.05 解使用二分法先要确定有根区间。
本题f(x)=x2-x-1=0,因f(1)=-1,f(2)=1,故区间[1,2]为有根区间。
另一根在[-1,0]内,故正根在[1,2]内。
用二分法计算各次迭代值如表。
其误差2. 求方程在=1.5附近的一个根,将方程改写成下列等价形式,并建立相应迭代公式.(1) ,迭代公式.(2) ,迭代公式.(3),迭代公式.试分析每种迭代公式的收敛性,并选取一种收敛最快的方法求具有4位有效数字的近似根解:(1)取区间且,在且,在中,则L<1,满足收敛定理条件,故迭代收敛。