xianxingxitong.m

来自「智能控制bp算法」· M 代码 · 共 27 行

clear;
a=[1 2 -1 2;2 6 3 0;4 7 -8 -5 ;7 2 1 6]; 
b=[-1 0 0 1]'; 
c=[-2 5 6 1]; 
d=7; 
[z,p,k]=ss2zp(a,b,c,d) 
ii=find(real(z)>0) % 检验零点的实部；求取零点实部大于零的个数 
n1=length(ii); 
jj=find(real(p)>0) % 检验极点的实部；求取极点实部大于零的个数 
n2=length(jj); 
% 判断系统是否稳定 
if(n2>0) 
disp(' 此系统不稳定极点是： ') 
disp(p(jj)) 
else 
disp(' 系统是稳定的！ ') 
end 
% 判断系统是否为最小相位系统 
if(n1>0) 
disp(' 此系统不是最小相位系统！ ') 
else 
disp(' 此系统是最小相位系统！ ') 
end 
% 绘制零极点图 
pzmap(p,z)