bass_pp.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function K=bass_pp(A,B,p)
a1=poly(p); a=poly(A);  % 求出原系统和闭环系统的特征多项式
L=hankel(a(end-1:-1:1)); C=ctrb(A,B);
K=(a1(end:-1:2)-a(end:-1:2))*inv(L)*inv(C);