ex524.m

详细介绍了Mathlab的使用方法与技巧

dnum=[3.4 -1.5];dden=[1 -1.6 0.8]; Ts=0.05;sysd=tf(dnum,dden,Ts);
num1=0.25;den1=[1 1];num2=1;den2=[1 2];
sys1=tf(num1,den1);sys2=tf(num2,den2);sysc2=sys1*sys2;
syscd=c2d(sysc2,Ts,'zoh');dsys=sysd*syscd;
dsysb=feedback(dsys,1);[dnum,dden]=tfdata(dsysb,'v');
%以上求出闭环Z域脉冲传递函数模型：dsysb；dnum,dden，Ts=0.05
dstepb(dsysb,Ts,0.05)                 %绘出单位阶跃响应图，得到系统响应的时域特性参数
n=80;t=0:Ts:5;                  %取采样次数为80；时间向量t：初值0，终值：5s
figure(2)
dstep(dnum,dden,n);            %单位阶跃响应图，注意横坐标以采样次数n为单位
figure(3);step(dsysb,t);     %单位阶跃响应图，注意横坐标以时间n为单位，t=n*Ts
figure(4);dstep(dnum,dden);  %绘出单位脉冲响应图，系统自动选择采样次数n
figure(5);lsim(dsysb,t);                    %绘制单位斜坡响应图
figure(6); plot((1/5)*sin(2*pi*t)); hold on;
dlsim(dnum,dden,sin(2*pi*t)); grid on;      %绘制 激励下的响应图