(2021年整理)ode45求解多自由度動(dòng)力學(xué)方程實(shí)例

1、ode45求解多自由度動(dòng)力學(xué)方程實(shí)例ode45求解多自由度動(dòng)力學(xué)方程實(shí)例 編輯整理：尊敬的讀者朋友們：這里是精品文檔編輯中心，本文檔內(nèi)容是由我和我的同事精心編輯整理后發(fā)布的，發(fā)布之前我們對(duì)文中內(nèi)容進(jìn)行仔細(xì)校對(duì)，但是難免會(huì)有疏漏的地方，但是任然希望（ode45求解多自由度動(dòng)力學(xué)方程實(shí)例）的內(nèi)容能夠給您的工作和學(xué)習(xí)帶來便利。同時(shí)也真誠的希望收到您的建議和反饋，這將是我們進(jìn)步的源泉，前進(jìn)的動(dòng)力。本文可編輯可修改，如果覺得對(duì)您有幫助請(qǐng)收藏以便隨時(shí)查閱，最后祝您生活愉快 業(yè)績進(jìn)步，以下為ode45求解多自由度動(dòng)力學(xué)方程實(shí)例的全部內(nèi)容。ode45函數(shù)調(diào)用形式如下：t,y=ode45（odefun,tspa

2、n，y0）相關(guān)參數(shù)介紹如下:參數(shù)名稱參數(shù)說明odefun用于存放待求解的方程的m文件名，方程必須用y=f（t，y）的形式存放tspan指定自變量范圍的向量,通常用t0,tf指定y0函數(shù)的邊界條件，即y0=y(t0），對(duì)于方程組，y0也可以是向量例：若一三自由度多體動(dòng)力學(xué)系統(tǒng)方程如下：初始條件：由于方程必須用y=f(t,y）的形式存放，因此需要對(duì)方程組進(jìn)行降階處理.令則方程組可化為：因此建立m函數(shù)文件來定義此方程組如下：function dy=func(t,y)dy=zeros(6，1)；dy（1)=y（2）;dy（2)=0。5*（-3y(2)+1.5*y(4)100*y（1）+50*y(3）+

3、2.0sin（3。754*t); dy(3)=y(4）;dy（4)=0.5（1.5y(2）-3*y(4)+1。5y(6)+50*y（1)100*y（3)+50*y（5)-2。0*cos（2。2t）);dy（5）=y（6)； dy（6)=0.5（1。5y(4)3y（6）+50*y(3）-100*y（5）+1.0sin(2。8t);end在matlab命令窗口里輸入一下命令:y0=1 1 1 1 1 1;tspan=0 30;t，y=ode45(func,tspan，y0）；figure(1）plot(t，y（：,1)，t,y(：,3）,t，y(:，5)；legend（x1,x2，x3）；xlabel（時(shí)間（s)，fontsize,10)；ylabel(振動(dòng)位移曲線,fontsize,10);figure（2)plot（t，y（:,2)，t,y（:,4）,t,y(:,6)）；legend(v1,v2,v3);xlabe