資源描述:
《矩陣特征問題的求解簡》由會員上傳分享,免費在線閱讀,更多相關內容在教育資源-天天文庫。
1、§1引言矩陣特征問題的求解物理、力學及工程技術領域中的許多問題在數(shù)學上都可歸結為求方陣特征值和特征向量的問題,如振動問題:橋梁的振動、機械振動、電磁振蕩、地震引起的建筑物的振動等等;物理學中的臨界點、臨界值的確定;微分系統(tǒng)中的穩(wěn)定性研究。這些常見的問題都與方陣的特征值和特征向量有關。求方陣的特征值和特征向量是數(shù)學、物理學、力學、電磁電工學以及工程技術所面臨的一個共同問題。定義:矩陣A?Rn?n,若有數(shù)?和非零向量x?Rn使得Ax=?x,則稱?為A的特征值,x稱為?對應的特征向量。在線性代數(shù)中的求法:解特征多項式
2、?E-A
3、=0.利用線性代數(shù)中的上述方法計算特征值和特征
4、向量是十分困難的;我們的目的是尋求一種適合計算機運行的近似解法,且簡單、可行、有效。復習相關理論定義設矩陣A,B?Rn?n,若有可逆陣P,使定理1若矩陣A,B?Rn?n且相似,則B=P-1AP則稱A與B相似。(1)A與B的特征值完全相同;(2)若x是B的特征向量,則Px便為A的特征向量。定理2:設A?Rn?n具有完全的特征向量系,即存在n個線性無關其中?i為A的特征值,P的各列為相應于?i的特征向量。的特征向量構成Rn的一組基底,則經相似變換可化A為對角陣。即矩陣A與對角陣?相似的充要條件是A有n個線性無關的特征向量即有可逆陣P,使定理3:A?Rn?n,?1,…,?n
5、為A的特征值,則(2)A的行列式值等于全體特征值之積,即(1)A的跡數(shù)等于特征值之和,即定義:設A?Rn?n為對稱矩陣,對任意x?Rn,x≠0,則稱比值矩陣A關于向量x的瑞雷(Rayleigh)商瑞雷(L.Rayleigh)——英國數(shù)學家,1842-1919定理4設A?Rn?n為對稱矩陣,其特征值?1≥?2≥…≥?n,則1)對任意x?Rn,x≠0,2)3)定理5(Gerschgorin圓盤定理)設A?Rn?n,則表示以aii為中心,以半徑為的復平面上的n個圓盤。(2)如果矩陣A的m個圓盤組成的并集S(連通的)與其余(1)A的每一個特征值必屬于下述某個圓盤之中,n–m個
6、圓盤不連接,則S內恰包含m個A的特征值。1乘冪法定理6設A?Rn?n有n個線性無關的特征向量,若?1,?2,…,?n為A的n個特征值且滿足對任取初始向量x(0)?Rn,(x(0)?0)對乘冪公式確定的迭代序列{x(k)},有下述結論:乘冪法與反冪法乘冪法是一種求矩陣的按模最大的特征值及其特征向量的方法。(1)當時,對i=1,2,…,n收斂速度取決于的程度,r<<1收斂快,r?1收斂慢,且x(k)(當k充分大時)為相應于?1的特征向量的近似值。(2)當時a)若?1=?2,則主特征值?1及相應特征向量的求法同(1);收斂速度取決于程度。向量,分別為主特征值?1、?2相應的
7、特征向量的近似值。b)若?1=-?2,對i=1,2,…,nc)若,則連續(xù)迭代兩次,計算出x(k+1),x(k+2),然后對j=1,2,…,n解方程求出、后,由公式解出主特征值?1、?2。此時收斂速度取決于的程度。向量、分別為相應于?1,?2的特征向量的近似值。乘冪法方法步驟設為n維向量,令r=max(x)表示向量x分量中絕對值最大者。1.任意給定初始向量(非零)v(0)=u(0)?0,取r0=max(v(0));2.產生迭代序列:為求A矩陣的按模最大的特征值?1和其相應的特征向量?1,v(1)=Au(0),取r1=max(v(1)),…………v(2)=Au(1),取r
8、2=max(v(2)),v(k)=Au(k-1),取rk=max(v(k)),3.循環(huán)次數(shù)控制:當
9、rk-rk-1
10、
11、T;r1=max(v(1))=4;u(1)=1/4(2,4,1)T=(0.5,1,0.25)T;v(2)=Au(1)=(4.5,9,7.75)T;r2=max(v(2))=9;…………結果如下:r[1]=4.000000,u(1)=(0.500000,1.000000,0.250000)T
12、r1-r0
13、=4.000000,r[2]=9.000000,u(2)=(0.500000,1.000000,0.861111)T
14、r2-r1
15、=5.000000,r[3]=11.444445,u(3)=(0.500000,1.000000,0.730582)T
16、r