ex3_20.m

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

% ex3_20
% 系统前向通道传递函数G(s)=10/(s2+s)；单位负反馈，计算动态误差系数
den=[1 1 0];
num=[10];
G=tf(num,den);
ES=feedback(1,G);
nc=ES.num;
dc=ES.den;
syms s nums dens
nums=poly2sym(nc{1});
dens=poly2sym(dc{1});
y=sym2poly(taylor(nums/dens));
ly=length(y);
for i=1:ly,
    yy(i)=y(ly-i+1);
    if yy(i)==0,
        k(i)=inf;
    else,
        k(i)=1/yy(i);
    end;
end; 
k