ex621.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 18 行

A=[-10 -35 -50 -24;1 0 0 0;0 1 0 0;0 0 1 0];
B=[1;0;0;0];C=[1 7 24 24];D=[0];
[z,p,k]=ss2zp(A,B,C,D,1);                       %求零极点增益模型
Flagz=0;                                     %设标志
n=length(A);
for i=1:n 
    if real(p(i))>0
        Flagz=1;
    end
end
disp('系统的零极点模型为');z,p,k
if Flagz==1
    disp('系统不稳定');
else
    disp('系统是稳定的');
end
step(A,B,C,D);                                %求系统状态方程的阶跃响应