完整版計算機科學計算答案

1、1、上下三角矩陣的乘積、逆仍為上下三角矩陣2、AB 與 BA 跡相同 tr(AB)=tr(BA),如果A或者B可逆,那么AB與BA特征值相同1)、tr(AB)n mai ,kbk,ii 1 k 1m nbk,iai, kk 1 i 1tr (BA)2)、由A 1( I AB)AI BA,或者B I ABB 1 I BA ,兩邊取行列式并令其為令,即得到證實.3、 有上條性質(zhì)可知：不存在滿足AB-BA=I條件的方陣A、B由于:tr(AB-BA)=tr(AB)-tr(BA)=0 豐 tr(I)=n4、A和B都嚴格對角占優(yōu),但是 A 土 B未必嚴格對角占優(yōu)(例： B=-A或者B=A )5、A和B可逆

2、時1) K cond(A),由于 II I =I AA11 | v | A | | A1III I兩邊消去I I I 即得2) 由 1)得到：| A | 1/ IA-1 I , I A1 I 1/ IA|3) 與 2)比照有 p (A) 1/ P (A)k _k k 4) 如果I A| v 1時 akA收斂,那么 akA 必發(fā)散,而| A| 1時 akA發(fā)散,那么 akA k未必收斂 k 06、cond(AB) cond(A)cond(B),利用范數(shù)相容性立即可得由此引出的不等式：I A-1-B-1 I I A-1IIB-1|IA-BIVk、? II- 1- 1 IIII 1 1UII 1 1

3、IIII1 1U由于 A - B = | B - A = | B (I- BA ) = | B (A- B)A IB-1 IIA- BIIA-1I對應地有 IA-BI I A|IB|IA1- B -1 I7、 A非奇異,B奇異,那么對于算子范數(shù)有1/Cond(A) | A-B | / | A |由于B奇異,那么存在 y乒0,使得By=0,從而有x=y/ I y |乒0, | x | =1,并且Bx=0,-1 _-1 _-1-1-1A Bx=0, x-A Bx=x, A (A-B)x=x , 1= | x | = | A (A-B)x | max A (A-B)x = x1=I A-1 (A-B

4、) I I A-1 I I A-B I ,兩邊除 I A| 即可得證.由6) -7 )組合,還可以得到更多的不等式.8、正規(guī)陣同時又是三角陣,那么它一定是對角陣9、 酉陣同時又是三角陣,那么它一定是對角陣,并且對角元的模為110、對稱矩陣的奇異值是特征值的絕對值n階實對稱矩陣如果有 n個互異的奇異值,那么它有 n個互異的特征值11、相似變換、酉變換、正交變換不改變方陣的跡和行列式由于上述變換不改變特征值,從而不改變特征值的和和乘積,從而也不改變跡與行列式.12、 酉變換、正交變換不改變向量的2-范數(shù),從而不改變兩點之間的歐幾里得距離.13、 酉陣特征值的模為1,正交陣特征值的絕對值為1 14、

5、x和e-x在任何區(qū)間a,b上線性無關(guān)設(shè)cx+de-x=0,由于e-x永遠不為零,如果 c乒0,那么有x/e-x=-d/c,顯然x/e-x在任何區(qū) 間都不會是一個常數(shù), 從而必須c=0,這樣一來d只有為零,因此只有當c=d=0等式才成立, 線性無關(guān).15、如果 I A| V 1 那么 I (I-A)-1 1/(1- I A|由于p (A)Ak收斂到s=0kA,兩邊乘(I-A)不影響收斂性, k 0(I-A ) s=(I-A)kA =I-lim(Ak 0k)=I,所以 Ak =(I-A) -1 ,(I-A)Ak 117、AC Cn*n,對任意范數(shù)有,lim圳Ak|(A)*首先存在某種范數(shù)(Ak)

6、|A (Ak),而(Ak) k(A)所以k(A) |Ak|k(A) k(A)(1/ k(A),取 k(A)得到k(A) |Ak|2 k(A),對不等式同時取極限即得到pm寸時| (A).*.*再根據(jù)范數(shù)的等價性cA | IA | c2|A I對不等式同時取極限即得到對任意范數(shù)的結(jié)果lim k |Ak(A)18、A非奇異,對算子范數(shù)有由于A (maxA、y|1 minA1y.AA1y minA 1yAxmin min Axx 111 119、.,tr(AHA)20、uavh21、11| 1,Cond( A) Cond(A)22、x,y為向量,那么有平行四邊形關(guān)系22y 2 x y 2)X：23、

7、A Cn*n,并且為Hermite正定陣,JxH Ax為范數(shù)正性和齊次性好證實,只證三角不等式A為Hermite陣有A UDU H,并且A正定, A UD1/2 xHAy于是對角陣d. 0,(ud1/2)hHd / Od / O |_|Id / O |_|_|d / O |_|1/2r1/2H/r1/2H H/c1/2Hx UD (UD ) y (UD ) x) (UD ) y)1/2、H 、H1/2、H(xHAx)(yHAy) |x | y ,同理可得所以x y2(x y)H A(x2 y22x y (x24、ahatA225、AH A26、A 為 Hermite陣,27、A 為 Hermi

8、te陣,a2yH Axy) xH Ax yH Ay xH Ay yH Axiyi)2Cond(AH A)2 Cond(A)22,A A1 nA?AIAI28、A c Cn*n,那么lim AkI的充要條件為A=Ik由于 I lim Aklim AAk 1 A lim Ak 1 AI 所以 A=Ikkk 1同理如果limAk B,并且B可逆,那么B=29、證實：1) A為實斜對稱陣,那么 eA為正交陣(A=-At)2) A 為 Hermite 陣,那么 eiA 為酉陣(AH=A)證實：1) eA(eA)TeA(eAT) eAe A e0I(eA)TeAiA / iA、H iA / (iA)H、

9、iA iAH iA iA 0iA、H iA2) e (e ) e (e ) e e e e e I (e ) e30、A Cn*n,A 2V1,證實 Ln(I+A) | 2 A 2/(1-A 2)由于：Ln(1 x) x x2 /2 x3/3(1)n 1xn/n收斂半徑r=1, p (A) I A| 2v 1,所以Ln(I+A)收斂,A21 A2A%I A|：kkl|Ln(I A)?k1Mlk IA%oM231、A 3 2 求：eAt,sinA,A1002 332、xn+1= 4 (x n)的迭代收斂條件,1)映內(nèi)性 xn C a,b , (f) (xn)C a,b 2)壓縮性 I 4 I L

10、v 1停機判據(jù)xkL xkxk 1 ,收斂速度xk1 LLk1 Lxx33、Newton迭代法,單根為二階收斂-j Xk 1f ( ) j xkxk 1lim7 c limsk x 22f ( ) k X x 2xk乙k /xk 1xk 234、Newton迭代法,重根變線性收斂,如果知道重數(shù)m, xk 1xk m f(xk)仍二階收f (xk)35、Newton 迭代法中,xk 1xk始點f (xQf (xk)下山因子入V,那么無法下山,要另選初36、弦割法xk 1 xkf(xk)(xk xk 1)f(xk) f(xk 1)的收斂階為1.61837、分半法的收斂速度為(b-a) /2n-1(

11、xk 1 xk)238、Aitken 加速公式 xk (xk 1), xk 1(xk),xkxk 1 xk 12 xkXk 139、Jacobi、Gauss-Seidel和超松弛(SOR)法的分量形式和矩陣形式(k 1)(J)xiaiin(k) aij xjj 1,j iD 1(L U)xk D 1b(G-S) x(k 1jbii 1(kaij xjj 11)n(k)、aijxj ) , xi 1(DL) 1Uxk(D L) 1b(k 1)(SOR) xi (bi aiii 1(kaij xjj 11)n(k)、 (k aijxj ),為j i 11)(1)xi(k)x(k 1)40、41、4

12、2、43、(D1L) (1k)D Ux (DL)1迭代法xk 1 Bxk f中,(B) 1時收斂,|B|1更收斂矩陣A嚴格對角占優(yōu),Jacobi法和(B) 1時可以得到超松弛法中 0Legendre正交多項式,x 1,1,Gauss-Seidel 法收斂2,(X)|b| 1時也是如此1,Pn(x)1fx22nn! dxn1非標準正交Pm(x)Pn(x)dx12n 1-(3x2 2遞推公式 Pn1(x) (2n xPn(x) /n前幾項 P01,P(x) x, P2 (x)n 1Gauss-Legendre機械求積公式11f(x)dx0, x0, A02,2f(0)1, x , A.3A11,f

13、(x)dx2,x 0, J|,AA25人 8-,A1-99有5階代數(shù)精度1)1 ,3 c、(5x3x),.2有1階代數(shù)精度f( I11f(x)dx1f (土),有3階代數(shù)精度3538539f(此)8f(0) 丁(伊22(1 x2)般地 咒 f 八 ,k 0,1,23.它的一般求積系數(shù)和零點復雜,但是權(quán)函數(shù) (n 1)Pn i(Xk)2n為1,在積分中不出現(xiàn),A 2k 044、Chebyshev正交多項式,11,1, (X) 2 ,Tn(x) cos(narccosx).1 x22k 1奇點為 xk cos( ), (k0,1,2,3,.n)遞推公式 Tn 1(x) 2xTn (x) Tn1(x

14、)2(n 1)前幾項：T01,T1 x,T2 2x2 1,T3 4x3 3x,T4 8x4 8x2 1非標準正交112 Tm (x)Tn (x)dx mmGauss-Chebyshev機械求積公式n0, xcos()0, Ao21 f(x)dx2x1.1f(0)有1階代數(shù)精度2k 1n 1,xo,icos()4Ai1 f2 1 1 x2dx (f(211.2).2)有3階代數(shù)精度-一 3 .n 2, x0, , Ao2A2A21 f (x) ,dx - (f (3 1 1 x23 33項 f(0) 0有5階代數(shù)精度2k 1般地Ak,xkcos( ), (k 0,1,2,3,.n)因此它的求積系

15、數(shù)和零點簡單,n 12(n 1)n缺點是權(quán)函數(shù)不整潔,A,k 045、利用兩點Chebyshev求積公式和利用 Legendre機械求積公式11f(x)dx1 .1 x2f (x)1 11 f()2 . 211 f(x)dxf(1、工 ,一-尸)兩者都3次代數(shù)精度,對f (x)3ax3bx2 cx d 精確ib a t)dt2bb a 1f(x)dxa2af(b),如何選擇a,b使之有更高的代數(shù)精度?令 f(x)=1、x,2 12, xxdx b3 032 ； 1 23-x2-,b -50 55于是 o xf(x)dx-f(3),此時 L&x2dx3 50-x2-(3)2 ,因此該公式最多有3

16、 5次代數(shù)精度147、如果 x f (x)dx af (b)cf (1),如何選擇a,b,c使之有更高的代數(shù)精度?令 f(x)=1、x, x2,xdx0I1-, xxdx ab0 3 02x546、利用三點Legendre機械求積公式47、如果 0 Txf (x)dxa , b1537,c再令f (x)7 A35(1)1 L 222 n 12xx dx ab c -x2-,萬程聯(lián)立解出070 7?,丁是0衣皿火175f(35f(1)51575因此該公式最多有 2次代數(shù)精度,不夠理想.48、Jim | Ak |A|,但是 lim AkkkA,例如:但是pm Ak49、Householder 鏡像

17、變換 H1 有 HT H H 1,HtH H2 I1,由于H變換特征值為土 1,奇異值為1, |det(H)|Hx x,x 2Txx,如果 Tx 0(正交),1;如果 Tx 0(平行),(1) x 2 Tx,(1)Tx2 TTx2Tx 12,150、利用H變換可將任意向量x,變?yōu)榈乳L度的向量y, iixi2II y2x x_yy2為鏡面H的標準法向量51、利用復化梯形公式和 Simpson公式求 ；cos(x)dx,如果誤差小于0.01區(qū)間應該分幾份?h2, h2h2Tn11制(f(b)f (a)芒 sin質(zhì))sin(0)節(jié) 0.01, h 0.34, (b a)/h -/0.34 4.6至少

18、要分5份,計算6個函數(shù)值SnI1h840(2)4(f(b)f(a)h4sin(-) sin(0)2880228800.01,h 2.32, (b a)/h -/2.32 0.67只利用原區(qū)間即可,計算 3個函數(shù)值.并不代表實際面積誤差, 特別以上公式在估算面積近似值時,要理解為有向面積的代數(shù)和,當(f (b)f(a)或(f(b)f (a)等零或者很小時,不能簡單套用.52、復內(nèi)積(x, y) yHxxiyi,由此：(x, y) (y,x),(x, y) (x, y)般實內(nèi)積(x,y)Wixiyi,由此:(f ,g)nWif(xi)f(yi)(g, f),(f, f)i 153、函數(shù)內(nèi)積(f,g

19、)b(x)f(x)g(x)dx a(g,f),(f,f) f|：54、正定矩陣:x 0,有xT Ax 0如果A正定,aH 0,A的順序主子式正定,A的任意主子式正定3n 255、Gauss消去法：乘除次數(shù) 一 n333n no( 一),加減次數(shù)一335n6t s2n3、2加減乘除總次數(shù)o(U)存儲量n23Gauss選主元增加比較次數(shù)3o(一),Gauss 選主3元消去法計算量為o(n3)56、Cholesky分解法：A (對稱)A=LDL T, (A又還正定),A=LL T,如果指定 A各對 角元符號(例如都大于零),那么A=LL T唯一.1. c乘除次數(shù)一(n 9n62n)3,n o(一),加減次數(shù)61 , 36(n3