ex8_1.m

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

% ex8_1
% 采用状态反馈配置闭环极点
a=[-1 3;0 -2];
b=[1 1]';
c=[1 0];
d=0;
G=ss(a,b,c,d);
p=[-5 -6];
syms k1 k2 s;k=[k1 k2]; 
eq=simple(det(s*eye(size(a))-a+b*k))
er=1;
for i=1:2,
    er=simple(er*(s-p(i))),
end;
er=er-eq;
[k1 k2]=solve(jacobian(er,'s'),subs(er,'s',0))