ex8_2.m

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

% ex8_2
% 已知系统的a，b和所需闭环极点p，确定状态反馈矩阵K
a=[0 1;-1 -2];
b=[1;-1];
p=[-3 -4];
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;
def=er-eq;
[k1 k2]=solve(jacobian(def,'s'),subs(def,'s',0))