bass_pp.m

来自「反馈控制系统设计与分析（matlab程序）」· M 代码 · 共 6 行

function K=bass_plc(A,B,p)
n=length(B); a1=poly(p);
alpha=[a1(n:-1:2),1];
a=poly(A); aa=[a(n:-1:2),1];
L=hankel(aa); C=ctrb(A,B);
K=(a1(n+1:-1:2)-a(n+1:-1:2))*inv(L)*inv(C);