ex7_23.m

来自「华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 」· M 代码 · 共 19 行

% ex7_23
% 计算模态矩阵
a=[-9 1 0;-26 0 1;-24 0 0];
b=[2 5 0]';
c=[1 2 -1];
d=0;
v=eig(a);
n=length(v);
syms lamda ad add; 
ad=(lamda*eye(n)-a);
add=inv(ad)*det(ad)
for i=1:n;
    ab=subs(add,'lamda',v(i)),
    T(:,i)=ab(:,n),
end;
A=inv(T)*a*T,
B=inv(T)*b,
C=c*T,
D=d