example8_3a.m

来自「这是一个MATLAB6。5 辅助控制系统设计与仿真例题的所有源码」· M 代码 · 共 15 行

% 判断系统可控性
n1=10;d1=conv(conv([1,0],[1,1]),[1,2]);
[a,b,c,d]=tf2ss(n1,d1);n=3;
CAM=ctrb(a,b);
if det(CAM)~=0;
    rcam=rank(CAM);
if rcam==n
    disp(系统可控')
elseif rcam<n
    disp('系统不可控')
end 
elseif det(CAM)==0
    disp('系统不可控')
end