【正文】 ndtn ??????? ??? ? ? （ 214） 參數(shù)設 定與求解 阻尼比 ? 分別取；應用 Matlab 對式 (211)和 (212),（ 214） 求解。需要對該力學模型進行深入的分析（借助 MATLAB 軟件進行仿真分析）。振動特性是振動篩非常重要的有別于非振動機械的一個本質特點，卻往往被設計者和制造者簡單化。因此，用高級語言編寫的程序一般代碼段較長，需要調用的子程序較多，整個程序的通讀性較差。應用 MATLAB的各種工具箱可以在很大程度上減小用戶編程時的復雜度，因此 MATLAB 在很廣的領域內得到了應用，其典型應用有；自動控制、圖像信號處理、生物醫(yī)學工程、語音處理、雷達工程、信號分析、振動理論、時序分析與建模、化學、統(tǒng)計學、經(jīng)濟學等 MATLAB 語言的優(yōu)點 MATLAB 作為一個以矩陣和數(shù)組為核心計算的軟件，對矩陣迭代法中的矩陣迭代計算尤其適合 [10]。這些年來，許多使用 MATLAB 的科學家、工程師和技術人員已經(jīng)開發(fā)出相當多的不同領域的應用程序。 (2)界面友好、編程效率高， MATLAB 的指令表達方式與標準教科書的數(shù)學表達式非常相近，用戶不需要有較高的計算機編程基礎，只要按照計算要求輸入表達式， MATLAB 將為用戶計算出結果。MATLAB 提供了良好的用戶界面，許多函數(shù)本身會自動繪制出圖形，而且會自動選取坐標刻度，繪制出直角坐標、極坐標、對數(shù)坐標下的二維和三維圖形，以及條形圖、直方圖、等高線圖、餅形圖、離散數(shù)據(jù)圖和瀑布圖等專用圖形。 MATLAB 的主要功能 :(1)數(shù)值計算功能，一條 MATLAB語句相當于幾十條 C 語言或 Fortran 語言的語句。 MATLAB 軟件的發(fā)展狀況及應用現(xiàn)狀 MATLAB 軟件概述： MATLAB 的名稱源自 Matrix Laboratory，是一門計算語言口 [7]。自 70 年代以來主要有： Ibrahim時域法（簡稱 LTD 法）、最小二乘復指數(shù)法 (LSCE 法 )、多參考點復指數(shù)法 (PRCE法 )、特征系統(tǒng)實現(xiàn)算法 (ERA)。模態(tài)參數(shù)辨識的頻域方法有：分量分析法、導納圓辨識方法、正交多項式曲線擬合、非線性優(yōu)化辨識方法等。模態(tài)測試技術主要采用同時測量輸入及輸出的方法，對一個振動系統(tǒng)來說，可以表示成圖 l1 所示的框圖 圖 模態(tài)分析框圖 Fig. Modal Analysis Diagram 通過測量激勵和響應，進行模念分析可以確定系統(tǒng)。模態(tài)分析是結構動力學中的～種 “逆問題 ”分析方法，它與傳統(tǒng)的 “正問題 ”方法（主要是指有限元方法）不同，是建立在實驗（或實測）的基礎上，采用實驗與理論相結合的方法來處理工程中的振動問題。線性振動的理論在發(fā)展過程中產(chǎn)生了一個重要分支，那就是模態(tài)分析理論。特別是多自由度系統(tǒng)的振動的理論，可以說既是振動力學的核心又是應用得最廣泛的振動理論。但是，建立振動系統(tǒng)的非線性力學模型難度大，求解困難，有些問題甚至無解可求。線性振動的一個重要特性是線性疊加原理成立。振動力學的物理基礎在 17 世紀已經(jīng)奠定，到了 18 世紀，振動力學已從物理學中獨立出來。對于低自由度的振動系統(tǒng)，方法一容易、快捷。固有頻率是決定系統(tǒng)振動特性的重 要物理量，它既是防止系統(tǒng)共振的依據(jù)，又是多自由度系統(tǒng)解耦分析（模態(tài)分析）的前提，因此研究某系統(tǒng)振動時，首先要求出系統(tǒng)的固有頻率。它在美國等發(fā)達國家的大學里已經(jīng)成為一種必須掌握的基本編程語言，而在國外的研究設計單位和工業(yè)部門，更是早己成為研究和解決工程計算問題的一種標準軟件。機械振動在許多情況下是有害的，人們想方設法避免它：另一方面，人們利用機械振動原理制造了各種機械或儀表來為人類服務。 many degrees of freedom main matrix iterative solution, our analysis based on abstract theory, while MATLAB programming The last iteration of data can be the desired data, so our calculations easier Using MATLAB programming and verify the correctness of the the process of operation, can quickly obtain multiple degrees of freedom vibration system and the main vibration mode natural frequency for the design to prevent resonance provide the theoretical basis for the preliminary analysis of the vibration of each ponent, and laid the decoupling of system response basis. Key words:vibrating system。單自由度頻率、阻尼、振型的分析，我們可以建立數(shù)學模型，最后通過利用 MATLAB 編程實現(xiàn)數(shù)據(jù)圖形；多自由度主要研究矩陣的迭代求解，我們在 分析抽象的理論的同時根據(jù) MATLAB 編程實現(xiàn)數(shù)據(jù)的迭代最后可以得到所要的數(shù)據(jù)，使我們的計算更加簡便。 I 基于 MATLAB 的振動模態(tài)分析 摘要 振動系統(tǒng)是研究機械振動的運動學和動力學，研究單自由系統(tǒng)的振動有著實際意義，因為工程上有許多問題通過簡化，用單自由度系統(tǒng)的振動理論就能得到滿意的結果。 利用 MATLAB 編程并驗證程序的正確性。 Single Degree of Freedom 。振動機械是 20 世紀后半期得到迅速發(fā)展的一類機械，它是利用振動原理來完成各種工藝過程的機械設備。在國內也有越來越多的科學技術工作者參加到學習和倡導這種語言的行列中來。主振型則為初步分析各構件的振動情況以及解耦分析奠定了基礎。但是在實際工程中，大多數(shù)振動 系統(tǒng)都是自由度較多，用特征矩陣方程式與特征方程式求解系統(tǒng)固有頻率與主振型這種傳纜的計算方法雖然從原則 上可行，但當自由度增加時，慣性、剛度陣的階數(shù)增高，計算量也急劇加大，這顯然很不方便。最主要的成就為線性振動理論的形成，它是與數(shù)學中的常 微分方程和偏微分方程同步發(fā)展的。非線性振動 [3]:描述其運動的方程為非絨性微分方程，相應的系統(tǒng)稱為非線性系統(tǒng)。在實際的工程應用中，很多情況下在誤差允許的范圍之內用線性的方法解決復雜的近線性問題 。線性振動在當今不僅是作為基礎科學的力學的一個重要組成部分，而且正走上向工程科學發(fā)展的道路，它在航空、機械、船舶、車輛、建筑、水利等工業(yè)技術部門中占有愈來愈重要的地位。在對選礦用振動篩進行分析時，需要通過實驗來驗證理論的正確性，振動實驗則需要用到模態(tài)分析技術。目前這一技術已發(fā)展成為解決工程中振動問題的重要手段，在機械、航空、航天、土木、建筑、造船、化工等遼寧工程技術大學畢業(yè)設計（論文） 3 工程領域被廣泛應用 [5]。自從 FFT 問世以來，目前廣泛采用寬頻帶激振技術。模態(tài)參數(shù)辨識的時域方法與模態(tài)參數(shù)辨識的頻域方法不同，它無需將所測得的響應與激勵的時間歷程信號變換到頻域中去，而是直接在時域中進行參數(shù)辨識。模態(tài)分析技術在動態(tài)載荷識別、模型修正與結構動力修改中有廣泛的應用，結構動態(tài)特征靈敏度分析是非常重要的方法之一。在工程計算領域，計算機技術的應用正逐步將科技人員從繁重的計算工作中解放出來。 (2)符號計算功能，利用 MATLAB 的符號計算功能可以清晰地獲得解的表達式，對于避免出錯和提高程序的可讀性均有很大的幫助。 (4)文字處理功能。同時使用 MATLAB 語言設計的程序，其編譯和執(zhí)行速度都超過了傳統(tǒng) c 和 Fortran 語言設計的程序，在工程計算方面的編程效率也高于其它編程語言。 MATLAB 的應用范圍： MATLAB 由主包和各種工具箱組成。就所查的資料看，以前的學者和研究人員迭代求解系統(tǒng)固有頻率與主振型時，大部分都是用 Visiual Basic 或 Fortran 語言來編寫程序 [11]。相反， MATLAB 則有簡潔、可讀性強等優(yōu)點?？陀^的說，一般的振動都是非線性的，但在許多情況下可以近似看作線性來處理。 本文主要利用 MATLAB 對振動系統(tǒng)進行模擬分析對于虛擬抽象的理論圖像化，處理單自由度振動的 3 個阻尼和強迫單自由度阻尼振動，多自由度系統(tǒng)振動矩陣迭代求解。程序如下： 王超：基于 MATLAB 的振動系統(tǒng)編程分析 8 clear,format pact。 k=1。r39。 w0=1。 figure(1),plot(t,y,39。t=0::18。 figure(1),plot(t,y,39。 1?? 時，振蕩系統(tǒng)等同于兩個一階系統(tǒng)串聯(lián)。可以利用下列 MATLAB 程序畫出 st? 在 0~ 范圍內 n? 和 n? 的變換曲線： % g=。 end plot(t,w)。 hold on。)。 title(39。圖所示的彈簧質量系統(tǒng)中，質量塊上作用有簡諧激振力 tPtP ?sin)( 0? (215) 遼寧工程技術大學畢業(yè)設計（論文） 11 其中 0P 為激振力幅， ? 為激振頻率。 將方程（ 215）的兩端同 除以質量 m ，并且令 22 nmc ??? （ 218） 其中 ? 為相對阻尼系數(shù)， n? 為相應的無阻尼系統(tǒng)的固有頻率，則方程（ 215）成為 xm? P(t) kx xc? c k x m P(t) 王超：基于 MATLAB 的振動系統(tǒng)編程分析 12 tmpxxx nn ???? s in2 02 ??? ??? （ 219） 上述方程特解可以通過 )sin( ?? ?? tBx 或者 tBtAx ?? s inc o s ?? 來求得，這里介紹用復數(shù)方法求式（ 219）的特解。 無阻尼系統(tǒng)對簡諧激勵的穩(wěn)態(tài)響應可以從式（ 226）得出。 beta=1./(sqrt((1lamda.^2).^2+(2*kesai*lamda).^2))。 plot(lamda,beta) hold on end axis([0 5 0 3])。 王超：基于 MATLAB 的振動系統(tǒng)編程分析 16 0 0 . 5 1 1 . 5 2 2 . 5 3 3 . 5 4 4 . 5 500 . 511 . 522 . 53頻率比振幅放大因子0 . 0 50 . 1 00 . 1 50 . 2 50 . 3 7 50 . 5 01 . 01 . 1 4 1 算例利用 MATLAB，繪制彈簧 質量系統(tǒng)在簡諧力作用下的響應曲線。 m=5。 f_0=F0/m。 xlabel(39。x(t)39。) 0 0 . 2 0 . 4 0 . 6 0 . 8 1 1 . 2 1 . 4 1 . 6 1 . 8 2 0 . 2 0 . 1 5 0 . 1 0 . 0 500 . 0 50 . 10 . 1 5tx(t)E x 3 . 1 1 本章小 結 基于 MATLAB 對單自由度自由振動繪制振動圖像，進行粘性阻尼，強迫振動振幅放大因子繪圖進行數(shù)據(jù)分析，使振動數(shù)據(jù)更加明顯。 式 （ 35）的分 母，叫做系統(tǒng)的特征方程。 將 i???? ?? 和i 帶入 公式 中，得 iii ??? ?A （ 39） 若將上式左端看作新列陣，上式表示：對于精確的主振型。為防止迭代過程中迭代列陣的元素變得過大或過小，每次迭代后需要使列陣歸一化，例如使它最后一個元素成為 1。 clear。 %*********迭代第 n 階主陣型 ******************** %n為計數(shù)器 % 遼寧工程技術大學畢業(yè)設計（論文） 23 n=1。 end %定義初始迭代向量 X(1) if n==1 X(:,1,n)=[1 1 1]39。 Y(:,1,n)=A(:,:,n)*X(:,1,n)。abs(X(2,i,n)X(2,i+1,n)) Y(:,i+1,n)=A(:,:,n)*X(:,i+1,n)。 end % X(:,:,n) f(:,n)=X(:,i,n)。 王超：基于 MATLAB 的振動系統(tǒng)編程分析 24 n=n+1。 X(:,:,1) disp(39。第三階主陣型的迭代結果 39。)。 本章小結 在工程振動中，確定系統(tǒng)固有頻率與主振型時是非常重要的。對微長度 dx ，有 dxxPdP ??? （ 42） x???? ??? tansin （ 43） 和 dxxxdd 22)t an ()s in ( ???????