example6_39.m

来自「<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.」· M 代码 · 共 26 行

for alph=[-1:1]
    alph
    num=[1,alph];
    den=[2 7 24 19];
    [A,B,C,D]=tf2ss(num,den);
    cam=ctrb(A,B);
    rcam=rank(cam);
    ob=obsv(A,C);
    roam=rank(ob);
    CAM=ctrb(A,B);
    rcam=rank(CAM);
    n=3;
    if rcam==n
        disp('System is controlled')
        if rcam<n
            disp('System is no controlled')
        end
        ob=obsv(A,C);
        roam=rank(ob);
        if roam==n
            disp('System is observable')
        else
            disp('System is no observable')
        end
    end