ex8_6.m

华东理工大学自动化系《控制系统分析、设计和应用》教材配套程序实例,有关说明： 1. 安装MATLAB软件。安装时

% ex8_6
% 确定状态反馈矩阵K
a=[-1 0 2;0 -2 1;1 0 3];
b=[2;0;1];
p=[-2+2*j -2-2*j -5];
syms s er;n=length(a);
co=ctrb(a,b);
er=1;
for i=1:n,
    er=simple(er*(s-p(i)));
end;
S=sym2poly(er),
delta=zeros(n,n);
for i=1:n;
    delta=delta+S(i)*a^(n-i+1);
end;
delta=delta+S(n+1)*eye(n);
z=zeros(1,n);
z(1,n)=1;
K=z*inv(co)*delta