example6_35a.m

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

A=[-2 2 -2;0 -1 0;2 -6 1];
b=[0;1;2];
S=ctrb(A,b);
if det(S)~=0
    s1=inv(S);
end
P=[s1(3,:);s1(3,:)*A;s1(3,:)*A*A]
P1=inv(P)
A1=P*A*P1
b1=P*b