www.art-1co0m.co.kr
Company Logo
脉冲激励响应函数impulse() – 绘制响应曲线。
impulse(sys, t)
– 不画图,通过函数返回值得到响应的相关数据。 [y,x,t]=impulse(sys,t) [y,x]=impulse(sys,t) y=impulse(num,den,t)
Copyright © by ARTCOM PT All rights reserved.
www.art-co3m.co.kr
Company Logo
将系统的特征方程式写成: 胡尔维茨行列式可列写为:
在an>o的情况下,系统稳定的 充要条件为:上述行列式的各阶 主子式均大于零,即
Copyright © by ARTCOM PT All rights reserved.
0 1 0 0
2],系统的输入为单位阶跃函数,求有初始条件的零输入响应曲线和总的
输出响应曲线。
A=[-2 -2.5 -0.5;1 0 0;0 1 0]; B=[1;0;0]; C=[0 1.5 1]; D=0; G=ss(A,B,C,D); t=0:.2:100; step(G,t) hold on y1=step(G,t); y=initial(G,[1 0 2],t); plot(t,y) xlabel('t'); ylabel('y')
www.art-com.co.kr
1.2 MATLAB判断稳定性的方法
Company Logo
ⅱ直接判别。直接求根
1.2 MATLAB判断稳定性的方法 如果我们只对绝对稳定性感兴趣,可使用代数判据 ⅰ把胡尔维茨行列式作为一个矩阵a输入。 ⅱ用下列语句计算出各阶行列式的值: [m,n]=size(a) for i=1:n b=a(1:i,1:i) h(i)=det(b) end • 求闭环特征方程的根。roots() • 为零极点模型,看极点是否在s右半平面。 • 对状态空间形式(闭环),求A阵的特征值eig(A)。 • 求解Lyapunov方程,命令格式为x=lyap(A,B,C)。