利用MATLAB进行系统时域分析,matlab控制系统时域分析教程

7.1 零输入响应

1、[调用格式] initial (sys,x0) initial (sys,x0,t)

initial (sys1,sys2….sysn,x0) initial (sys1,sys2….sysn,x0,t)

initial (sys1,’plotstyle1’,sys2, ’plotstyle2’,….,x0)

[y t x]=initial(sys,x0)

Sys—输入系统 x0—初始状态下状态变量

t—采样时间矢量 plotstyle—绘制零输入响应曲线的特性

[y t x]—表示输出矢量 采用时间 状态变量

重要提示：initial命令适用于状态方程模型sys

2、例：根据传递函数求解单一系统的零输入相应

初始状态下状态变量为

num=[0 0 1 1];den=[1 9 26 24];

[a b c d]=tf2ss(num,den);

sys=ss(a,b,c,d);

x0=[1 1 0]';

t=0:0.01:4;

initial(sys,x0,t);

7.2 阶跃响应

[调用格式] step(sys) step(sys ,t)

Step(sys1,sys2,….sysn) Step(sys1,sys2,….sysn,t)

step(sys1,’plotstyle1’,sys2, ’plotstyle2’,….,x0)

[y t x]=step(sys,x0)

例：根据状态方程输出多系统的阶跃响应

系统一：

系统二：

a1=[-4 -3;1 0];b1=[1 0]';c1=[3 0];d1=[2];

sys1=ss(a1,b1,c1,d1);

a2=[-5 -3;1 0];b2=[1 0]';c2=[1 1];d2=[0];

sys2=ss(a2,b2,c2,d2);

step(sys1,'r-',sys2,'b-.');

7.3 脉冲响应

[调用格式] impulse(sys) impulse(sys ,t)

impulse(sys1,sys2,….sysn) impulse(sys1,sys2,….sysn,t)

impulse(sys1,’plotstyle1’,sys2, ’plotstyle2’,….,x0)

[y t x]=impulse(sys,x0)

例：根据状态方程输出单一系统的脉冲响应

a=[-11 -44 -74 -40;1 0 0 0 ;0 1 0 0;0 0 1 0];b=[1 0 0 0]';

c=[0 0 20 0];d=[0];

sys