当前位置：文档之家 › 《数值分析》杨大地 答案(第四章)

《数值分析》杨大地 答案(第四章)

  • 格式：docx
  • 大小：48.21 KB
  • 文档页数：2
6 2 1 1 1 1
⑴ 2 3 1 ，求按模最大特征值和对应的特征向量，精确到小数三位。
解：由幂法公式有:
1i n
//P67 页
(1) | a r | max | a i | ,其中 ai 是uk −1 = (a1 , a2 , . . . , an )T 的各分量； (2) y k 1
∴ i t 为������ − ������������ 的特征值。 令 为 A 的特征向量，则有： A i 又∵ A tI A tI i t i t ∴ 也为������ − ������������ 的特征向量； ∴ i t 是 A tI 的特征值，且 A 和������ − ������������特征向量相同。
T
1
当计算到第 9 次时，λ 1 的小数点前三位精度开始稳定，满足题目要求，所以此时 A 矩阵的按模最大特 征值 1 =7.288，对应的特征向量为（1.000,0.523,0.242）T 5．若 A 的特征值为 1, 2 ,, n , t 是一实数，证明： i t 是 A tI 的特征值，且特征向量不变． 证明： ∵ A 的特征值为 ∴｜ i I A ｜=0 假设������是 A tI 的特征值，则有：
k 0 1 2 3 4 5 6 7 8 9 ukT 1.000,1.000,1.000 9.000,6.000,3.000 7.6667,4.3333,2.000 7.3913,3.9565,1.8261 7.3177,3.8530,1.7824 7.2966,3.8231,1.7701 7.2906,3.8146,1.7666 7.2887,3.8119,1.7655 7.2882,3.9112,1.7652 7.2880,3.8109,1.7651 ykT 1.000,1.000,1.000 1.000,0.6667,0.3333 1.000,0.5653,0.2609 1.000,0.5353,0.2471 1.000,0.5265,0.2436 1.000,0.5240,0.2426 1.000,0.5232,0.2423 1.000,0.5230,0.2422 1.000,0.5229,0.2422 9.0000 7.6667 7.5913 7.3177 7.2966 7.2916 7.2887 7.2882 7.2880
12 ...n ，设 为 A 的特征值（i=1,2,...n） i
∴ห้องสมุดไป่ตู้
������������ − ������ − ������������ = 0 ������ + ������ ������ − ������ = 0 ������ + ������ = ������������ ������ = ������������ − ������
数字分析第 4 章
4.1 填空题： (1) 幂法主要用于求一般矩阵的按模最大特征值，Jacobi 旋转法用于求对称矩阵的全部特征值； (2) 古典的 Jacobi 法是选择绝对值最大的一对非对角元素将其消为零; (3) QR 方法用于求一般矩阵的全部特征值，反幂法加上原点平移用于一个近似特征值的精确化和求 出对应的特征向量． 4.2 用幂法求矩阵．
u k 1 ； | ar |
(3)������������ = ������������������−1 = ������1 ������������−1 ；
~ 是 yk-1 绝对值最大的分量，ak 是 uk 绝对值最大的分量； (4) t k ~ k , 其中 a k 1 a
k 1
a
(5) | t k t k 1 | ，令 1 t k ,x1=yk-1,退出运算，否则返回(1)重做以上步骤。 假设 U0 ={1.000,1.000,1.0000}，有

合集下载

数值分析习题解答

数值分析习题解答

6.（1）设(1,0,5,2)Tx =-，试求12,,x x x∞（2）设40004402A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦，试求12,,,F A A A A ∞ 解12128,5;6,8,FxxxA AAA∞∞=======；4．设05813622,10612422A b ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦， （1）试对A 进行PLU 分解：PA LU =； （2）根据PLU 分解求解Ax b =。

解 （1）162201011,102,00100.517100L U P ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪=== ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭（2）(1,1,1)Tx =8．分别用Householder 变换法和MGS 法对A 进行QR 分解⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=542112111A解 (1) Householder 法对A 进行QR 分解[]()()()123123,,,1,2,2,1,1,4,1,1,5T T TααααααA ===--=-令()11,2,2Tαα==，调用算法2.1有[]13,,42212Tu ρβ=-==，所以 []1122333100412210102422123330012212333T uu β---⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥--⎢⎥⎢⎥⎢⎥H =I -=-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦--⎢⎥⎢⎥⎣⎦, 故1333003033--⎡⎤⎢⎥H A =-⎢⎥⎢⎥-⎣⎦再令()0,3Tα'=-,调用算法2.1得20110H ⎡⎤'=⎢⎥⎣⎦,则 2100001010⎡⎤⎢⎥H =⎢⎥⎢⎥⎣⎦,21333033003--⎡⎤⎢⎥H H A =-⎢⎥⎢⎥-⎣⎦故121223331212,0333221003T TQ R -----⎡⎤⎡⎤⎢⎥⎢⎥=H H =--=-⎢⎥⎢⎥⎢⎥⎢⎥---⎣⎦⎣⎦. 10．设131112000,110001A b ⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦求Ax b =的最小二乘问题的全部解。

数值分析课后习题及答案

数值分析课后习题及答案

第一章 绪论（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 ，并估计插值余项。

《数值分析》第四章答案

《数值分析》第四章答案

习题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 =。

数值分析课后习题答案

数值分析课后习题答案


0 1
0 10 1 1 0 0 0 1
0 0 12 1 1 2 0 0 0

1 2
0 0 0 1 1 0
1 2

1 2


1 2
1
0 0 0 1 0

1 2

1 2


0
1 2

1 2
0
0
0
341 1 1
2-5.对矩阵A进行LDLT分解和GGT分解，并求解方程组
Ax=b，其中
16 4 8
1
A 4 5 4 , b 2
8 4 22
3
解
16 A 4
4 5
84
44 11
2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中
2 1 1 A1 3 2
4 ,b6
1 2 2
5
解
2 A 1
1 3
1 2


2 11
22
1
5 2
1

3 21来自，所以 A12
1
2 1 1



5 3
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
解
3 2 6 4 10 7 0 7 10 7 0 7

r1r2
消元

10 7 0 7 3 2 6 4 0 0.1 6 6.1
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623