2基于数值积分法的控制系统仿真

《2基于数值积分法的控制系统仿真》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、2基于数值积分法的控制系统仿真实验目的李德和092220381理解并学握基于数值积分的控制系统仿真方法。2初步理解利用计算机仿真分析控制系统的方法.3进一步理解2阶（高阶）系统的时间响应特性。4锻炼编写调试MATLAB/C程序能力。实验内容1对下列传递函数进行matlab仿真。(1)G】(s)=1(5s+1)(10s+1)(20s+1)利用RK-4法，程序如卜:clear;clc;num=[1];den=conv(conv([51]A[101])f[201]);gl=tf(num,den);[A,B,C,D]=ssda(gl);xO=[0;0;0];%

2、3oOpxk=xO;u=1;h=1/50;ST=100;iST=round(ST/h);t=zeros(iST,1);y=zeros(iSTf1);for(ih=1:iST)t(ih+l)=t(ih)+h;KI=A*xk+B*u;xkp=xk+0・5*h*Kl;K2=A★xkp+B*u;xkp=xk+0・5*h*K2;K3=A★xkp+B*u;xkp=xk+h*K2;K4:=A*xkp+B*u;xkp=xk+h/6*(KI+2*K2+2*K3+K4);y(ih+1)=C*xkp;xk=xkp;endfigure();plot(tzyf1r1z1line

3、width1z3);ylabel(*yf;grid;holdon;gl=ss(A,BzCzD);[sy,st]=step(gl);plot(st,sy,*b*);legend(1RK-4'「'Step*);仿真结果：(2)G

4、(s)=S+1_(5s+1)(10s+1)(20s+1)利用RK-4法，程序同上，结果如下:(2)G](s)=s-1(5s+1)(10s+1)(20s+1)利用RK-4法，程序同上，结果如下:2欧拉法编程如卜:clear;num二⑴；den=conv(conv([51],[101]),[201]);g=tf(num,den);[A

5、,B,C,D]二ssdata(g);h=1/50;u=1;xO=[0;0;0];xk=xO;y(l)=C*xk+D*u;t(l)=0;%初始化for(ih=l:80/h)t(ih+l)=t(ih)+h;kl=A*xk+B*u;xkp=xk+h*kl;y(ih+l)=C*xkp+D*u;xk=xkp;endfigure();plot(t,y,'r','linewidth',3);ylabel('y');grid;holdon;[sy,st]=step(g);plot(st,sy,'b‘);legend('eu'z'Step');