ex8_9.m

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

% ex8_9
% 状态反馈等效回路的计算及其他问题
G=tf(10,[1 6 5 0]);
ts=5.65;
Mp=0.043;
syms s k1 k2 k3 K;
k1=1;
lu=4/ts;
tp=-log(Mp)/lu;
wd=pi/tp;
beta=atan(wd/lu);
s1=-lu+beta*j;
s2=- lu-beta*j;
s3=floor(-lu*10)
S=sym2poly((s-s1)*(s-s2)*(s-s3))
d0=10*k1*K-S(4);
d1=5+10*K*(k2+k3)-S(3);
d2=6+2*k2*K-S(2),
digits(5);
[K,k2,k3]=solve(d0,d1,d2,'K','k2','k3');
K=vpa(K),k2=vpa(k2),k3=vpa(k3)