摘要-NTOU-海洋大学力学声响振动实验室.docVIP

河工系 雅鈞 題目:以力學觀點來看高斯消去法的矩陣操作 摘要 本計劃旨在培育大學生數學與工程的能力，使其以力學觀點看數學操作，以數學來建構物理模型。對其矩陣的數學操作之幾何與結構力學意義做一跨領域串聯。讓大學部學生除了學習基礎數學知識外，對其結構力學中次結構、子結構或超元素的相關物理概念，亦能有所體會。此外，再利用Mathematica軟體的動畫功能將各方法所推算出的理論結果以動畫呈現建構理論模型 在此，使用能量法來建構彈簧勁度矩陣。首先，求出各彈簧的相對位移、，可知 （1） （2） 其中、、各為a、b、c點的位移量。再利用系統的總勢能 （3） （4） （5） 整理後得 （6） 以及系統趨向最低能量 （7） （8） 將（7）式及（8）式整理後寫成矩陣形式 （9） 再將（9）式經由高斯消去法操作後得 （10） 由（10）式發現其化簡後之矩陣，內含等效彈簧常數，可將原串聯彈簧組，視為一等效彈簧，如圖二所示。 2.探討多自由度運動系統，並利用Mathematica來模擬其運動狀況 (1)兩端自由 以下例子皆為無阻尼之振動系統，且無外力作用。利用運動方程式求各振動頻率及振動模態。 已知 , （11） 先將各質點所受的力做分解，畫出自由體圖，如圖四所示。 藉由力平衡可知 質點 a： （12） 質點 b： （13） 質點 c： （14） 在此使用到運動方程式 （15） 我們將(12)-(14)式聯立可得 （16） 當質量與彈簧勁度皆取1單位時，可得質量矩陣 （17） 與勁度矩陣 （18） 利用數學解特徵值的方法 （19） 來求解振動頻率，三組振動頻率分別為0、1與，而與這三個振頻率所對應的振動模態依序為 （20） (2)左端固定，右端自由 此系統如圖五所示，為三個自由度且無外力作用之運動系統，探討系統中之振動頻率及振動模態 作法同上個範例，先對各個質點取自由體圖，討論其受力情形，如圖五所示。 已知 （21） 圖六 自由體圖 藉由力平衡可知 質點a： （22） 質點b： （23） 質點c： （24） 在此使用到運動方程式 （25） 為了列出質量與勁度矩陣，所以將各質點所受的力整理並得到以運動方程式的形式表現 （26） 當質量與彈簧勁度皆為1單位時，可得質量矩陣 （27） 與勁度矩陣 （28） 經由計算後，得知此三組振動頻率分別為0.445 、1.247 與1.803，而三組振動模態為 （29） 3.當樞紐元素為零時可藉由擴大矩陣法來解決 在高斯消去法的操作過程中若遇到某個樞紐元素為零時，利用擴大矩陣法可解決，再將矩陣分解成LDU的形式可讓之後運算更為簡易。 （30） 其中，為下三角矩陣；為對角矩陣；為上三角矩陣。若一聯立方程式為 （31） 令其中 （32） 經高斯消去法後，發現第二行第二列的樞紐元素為零，所以將擴充為，再將其矩陣作LDU分解，得 （33） 4.高斯消去法 一般高斯消去法是求解n元一次聯立方程組，利用消去法與代入得其解。當直接利用高斯消去法在進行運算時，發現有以下兩個問題 (1)樞紐元素為零時，無法直接繼續操作。 （34） 可利用行列互換的操作，解決此問題。第一步：將第一行元素取絕對值，則把絕對值最大所對應的那列元素與第一列元素交換 。 （35） 第二步依此逐層操作，如此才能使高斯消去法順利進行。 (2)捨位誤差： 以下這個例題，先取到小數點後三位時 （36） 第一步，化成上三角矩陣 （37） 第二步，反代回去求解 （38） 若再取到小數點後四位時，同理可得，先化成上三角矩陣 （39） 再反代回去求其解，得 （40） 因此，發現取的小數位數越多，得到的值越精準。 5. 直接利用反矩陣求解 （41） （42） 若A為病態矩陣，則此方法會求解不準。以下將介紹兩種方法來處理病態矩陣之問題。 (1)利用SVD法來求解 已知 （43） 我們可設 （44） A可分解為 （45） 其中，左酉矩陣為 （46） 對角矩陣為 （47） 右酉矩陣為 （48） 因此， （49） 令 （50） （51） 整理，得 （52） （53） 將條件代入，得 （54） （55） （56） (2)利用特徵值及特徵向量 利用上面的例子，求解特徵值及特徵向量。A可分解為 （57） 其中，D為對角矩陣 （58） V為特徵向量所建構而成的矩陣 （59） 令 （60） （61） 將條件代入整理，可得 （62） （63） （64） （65） 藉由上述的兩種方法我們可以解決病態