單擺運(yùn)動(dòng)過程的matlabsimulink建模與仿真-閱讀頁(yè)
2025-04-24 11:31本頁(yè)面
【正文】 （5） （6）仍然使用歐拉算法求解：將dv(t) = v(t + dt) v(t)和d(t) = (t + dt) (t)代入式（5）及式（6）中，并以仿真步進(jìn)量Δ作為dt 的近似，得到基于時(shí)間的遞推方程：v (t +Δ)= v (t)+(gsin( t) )Δ (t +Δ)=(t)Δ據(jù)此編寫仿真程序：在MALAB命令窗口輸入以下命令：subplot(2,1,1)dt=。 %仿真時(shí)間長(zhǎng)度t=0:dt:T。L=。k=3。 %初始擺角設(shè)置,不能超過π/2v0=0。 %程序存儲(chǔ)變量預(yù)先初始化，可提高執(zhí)行速度th=zeros(size(t))。th(1)=th0。th(i+1)=th(i)1./L.*v(i).*dt。plot39。set(B1,39。,39。)。LineStyle39。:39。set(get(AX(1),39。),39。,39。)。Ylabel39。String39。角位移\th(t)/rad39。xlabel(39。)。線速度v(t)39。legend(B2,39。,1)。 %仿真步進(jìn)T=16。%仿真計(jì)算時(shí)間序列g(shù)=。m=8。 %空氣阻力比例系數(shù)th0=pi。 %初始擺速設(shè)置v=zeros(size(t))。v(1)=v0。for i=1:length(t) %仿真求解開始v(i+1)=v(i)+(g*sin(th(i))k./m.*v(i)).*dt。endplot(t,v(1:length(t)))xlabel(39。)。線速度v(t)39。仿真圖象如下所示:由上圖可以看出，單擺圍繞懸點(diǎn)做圓周運(yùn)動(dòng)，速度在最高點(diǎn)的時(shí)候最低，在最低點(diǎn)的時(shí)候最大，現(xiàn)在驗(yàn)證最大速度值:帶入仿真初始值：v0=10 g= L=,與仿真結(jié)果明顯一致。2. 仿真的關(guān)鍵在于充分了解系統(tǒng)的物理化學(xué)等特性，并在此基礎(chǔ)上建立系統(tǒng)的數(shù)學(xué)模型。4. MATLAB作為一種輔助工具，功能十分強(qiáng)大，但是使用起來很簡(jiǎn)單，需要深
公司管理相關(guān)推薦
鄂ICP備17016276號(hào)-1