matlab教程第八章simulink交互式仿真集成環(huán)境(參考版)
2024-08-25 13:34本頁面
【正文】 圖 37 圖 （ 6） 圖 38 圖 。angzzy=[0。 （ 4） 圖 36 圖 圖 （ 5） clear dampzzy=。 % mdlOutputs function sys=mdlOutputs(t,x,u) sys = x(2)。 sys =dx 。 % mdlDerivatives function sys=mdlDerivatives(t,x,u,dampzzy,gravzzy) dx(1)=dampzzy*x(1)gravzzy*sin(x(2))+u。 str = [ ]。 sys = simsizes(sizes)。 = 0。 = 1。 = 2。,num2str(flag)])。 otherwise error([39。 case 3, sys=mdlOutputs(t,x,u)。 case 1, sys=mdlDerivatives(t,x,u,dampzzy,gravzzy)。 (algebraic variable) 有代數(shù)環(huán) 無代數(shù)環(huán) 帶記憶 塊 34 0 1 2 3 4 5012with loopwithout loopmemory 圖 S 函數(shù)模塊 概述 SIMULINK 的運(yùn)行機(jī)理 用 M 文件表述的 S 函數(shù) S 函數(shù)的模板程序 S 函數(shù)模塊的創(chuàng)建和使用 【例 】本例演示：（ 1）用 S函數(shù)模塊為圖 型；（ 2）利用 SIMULINK 研究該單擺擺角 ? 的運(yùn)動(dòng)曲線；（ 3）用 S 函數(shù)動(dòng)畫模塊表現(xiàn)單擺的運(yùn)動(dòng)。 39。 39。 39。 contains 1 algebraic loop(s). Found algebraic loop containing block(s): 39。,4) Warning: Block diagram 39。,39。,39。),hold off legend(39。,2) plot(t3,y3,39。,39。,6),hold on plot(t2,y2,39。,39。]) disp([T1,T2,T3]) plot(t1,y1,39。,blanks(3),39。,blanks(3),39。 disp([blanks(4),39。,5)。[t3,x3,y3]=sim(39。T2=toc。exm080732_139。 tic。,5)。[t1,x1,y1]=sim(39。（ C）與前兩個(gè)模型做性能比較。演示：（ A）利用“記憶”模塊中斷代數(shù)環(huán)。演示：通過重組模型，直接 消除代數(shù)環(huán)，建立等價(jià)模型。 計(jì)算容差的選擇 代數(shù)環(huán)問題 代數(shù)環(huán)的形成 【例 】研究方程組 xkxkyykux ?? 213 , ???? 的解算問題。,0) ns1=length(xx1) ns2=length(xx2) ns1 = 15072 ns2 = 101 32 246 247x 104SymbolicODE45 ODE15S 圖 積分步長和容差 積分步長的選擇 【例 】以例 。,39。,39。) axis([246 247 ]) legend(39。,tt2,xx2(:,2),39。,tt1,xx1(:,2),39。 plot(t,x,39。exm080714_139。)。,39。 opts=simset(39。exm080714_139。)。,39。 opts=simset(39。x=eval(char(dxsym))。t39。t39。x(0)=1,Dx(0)=039。D2x+100*Dx+*x=039。本例演示：對(duì)于 Stiff方程，如果解算方法選擇不當(dāng)將產(chǎn)生嚴(yán)重后果。New39。Old39。LineWidth39。r39。LineWidth39。b39。 [ynew,tnew]=step(new,50)。 new=tf(1,coeff)。 a=fminunc(itae,a0,options)。off39。LargeScale39。 Tspan=(0:500)/10。 % 7 ss=s(end)。smodel39。,)。 % 5 opts=simset(39。Tspan39。base39。 （ 1） （ 2） 圖 （ 3） [] function ss=itae(aa) % global a ss a=aa。（ B） SIMULINK方塊模型、目標(biāo)函數(shù)和優(yōu)化程序之間的協(xié)調(diào)和參數(shù)傳遞。eig_Axb=(eig(Axb))39。exm080634_139。eig_Axa=(eig(Axa))39。exm080634_139。,[1,2]39。) xb=trim(39。exm080634_139。x239。x139。x239。x139。 [X1,X2]=portraitzzy(x1,x2,h)。x2=x1。 （ 5） h=。DX2(jj,ii)=dx2/L。 Z(jj,ii)=L。 dx2=x(2,2)x2(jj)。,h,opts)。% 12 [t,x,y]=sim(39。initialstate39。 for jj=1:n。X2=X1。 % 7 n=length(x1)。fixedstep39。ode539。solver39。x239。)。 end xlabel(39。,10,opts)。 27 [t,x,y]=sim(39。initialstate39。 nxx=size(xx,1)。1,1。1,1。1,1。 （ 1） （ 2） [] % clf。（ C）“一步仿真”計(jì)算方法。本例綜合演示：（ A） SIMULINK模型和 MATLAB指令的配合使用。,eA1=eig(A1)39。,[1,])。A A = 0 0 （ 4） [A1,B1,C1,D1]=linmod(39。exm080634_139。,4) 0 20 40 60 80 100020406080100120內(nèi)初值外初值 圖 模型的線性化問題 線性化的數(shù)學(xué)描述 連續(xù)系統(tǒng)的線性化模型 離散系統(tǒng)的線性化模型 模型線性化的算例
法律信息相關(guān)推薦
鄂ICP備17016276號(hào)-1