bass_pp.m

来自「包含大量遗传算法程序」· M 代码 · 共 5 行

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);