飞行器固有特性分析课件.ppt

?(4)特徵值的移軸性質式（6-21）兩邊分別減去，則有另一等價形式：或寫為式中式（6-21）和式（6-32）有相同的特徵向量，但特徵值相差μ，即：（6-31）（6-32）（6-33）6.3結構動態特性分析?(5)特徵值的分隔性質作移軸，並將作三角分解。則對角矩陣D中有i個負元素。如果有6.3結構動態特性分析?(6)位移展開定理以上討論的是廣義特徵值問題的一些基本特性，深入理解這些性質，對於求解特徵值問題很有幫助。對於n維空間中的任意向量｛x｝都可以按模態矩陣Φ展開：係數q可按下式確定：（6-34）（6-35）6.3結構動態特性分析6.3.2迭代法任意選取適當的初始向量{x1}，按迭代格式則向量序列{x1}{x2}將收斂於相應的特徵向量。向量迭代法又稱冪法，它既可用於標準特徵值問題，也可用於廣義特徵值問題。它不僅適合於對稱矩陣，也適合於非對稱矩陣。（6-36）6.3結構動態特性分析對任意向量則：按迭代格式式（6-19）有則：由於當k增大時，｛xk+1｝可能會變得很大或很小，因此，在迭代過程中，需要將迭代向量規一化（6-37）（6-38）（6-39）（6-40）（6-41）6.3結構動態特性分析?迭代法的迭代過程是自校正的。迭代向量中的誤差只能延遲收斂，而不會破壞收斂性。?根據特徵值的移軸性質可以構造帶移軸的向量迭代方法。只要選取合適的移軸量，就可以既使迭代收斂到所需要的特徵對，又可以加快收斂速度。6.3結構動態特性分析在廣義特徵值問題中，品質矩陣M是對稱正定，則一定存在非奇異矩陣6.3.3變換法?1.廣義特徵值問題化為標準特徵值問題可得到標準特徵值問題所以在上式中，前乘S-1，並令則有式中（6-45）（6-46）（6-47）6.3結構動態特性分析（6-48）?2.標準特徵值問題的變換法標準特徵值問題常用的變換法有雅可比方法（Jacobi）、Givens方法、Householder方法，在一般的矩陣代數教材中均有詳細敘述，在此我們只對雅可比方法簡要介紹。考察標準特徵值問題在經過k次變換後，有6.3結構動態特性分析（6-49）雅可比方法的思想是經過多次旋轉正交變換使矩陣對角化，對應的正交矩陣為6.3結構動態特性分析（6-50）?當則?當則，符號取決於的符號。?當時，矩陣A趨向於對角陣。6.3結構動態特性分析（6-51）?如果對廣義特徵值問題中的剛度矩陣K和品質矩陣M同時用雅可比方法作變換，得到矩陣M、矩陣K共同的主軸，則這種方法稱為廣義雅可比方法。?Givens方法與雅可比方法類似，也是進行座標旋轉變換，但它不是把實對稱矩陣A對角化，而只是三對角化。?Householdes方法也是一種將實對稱矩陣化為三對角陣的方法。6.3結構動態特性分析6.3.4Sturm序列二分法對於給定的對稱三對角矩陣把矩陣三對角化後，還需要求解三對角矩陣的特徵值問題。三對角矩陣求特徵值比一般矩陣要容易得多，常用的方法有QR方法和Sturm序列二分法。對於標準特徵值問題的三對角矩陣，則常採用Sturm序列的二分法。通常，用Householder法和Lanczos法將標準特徵值問題轉化為三對角矩陣。6.3結構動態特性分析（6-57）其特徵值的行列式為定義其零階主子式6.3結構動態特性分析（6-58）（6-59a）以後的各階主子式為p0,p1,p2,…，pn構成一個多項式序列，它是Sturm序列。6.2結構動態特性分析（6-59b）6.3.5大、中型特徵值問題的求解方法?1.Rayleigh-Ritz分析在實際問題中，人們最關心的不是全部特徵對，而只是其中的一小部分，例如，最低的前q階（qn）特徵對。在這種情況下，可以用一種近似的有效方法，將n階廣義特徵值問題化為q階廣義特徵值問題，這就是Rayleigh-RitZ分析法，或稱為Ritz變換法。6.2結構動態特性分析若要求系統的前p階特徵對，則先選取q≥p個線性無關的向量｛yi｝，i＝1，2，…，q，令｛x｝為這些向量的線性組合，有由式（6-10）得由Rayleigh商性質，從式（6-13），有6.2結構動態特性分析