fig6_30.m

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

% fig6_30
% 速度反馈控制系统的设计
% 与原系统比较
figure('pos',[50,50,200,200],'color','w');
a=1;
G=tf([1],[1 a 0]);
Mp=0.2;
zeta=0:.1:0.9;
mp=exp(-pi*zeta./sqrt(1-zeta.*zeta)); 
zetan=spline(mp,zeta,Mp);
tp=1;
wd=pi/tp;
wn=wd/sqrt(1-zetan*zetan);
Ko=wn*wn;
Kf=(2*zetan*wn-a)/Ko;
Gc=tf([Kf 0],1),
GL=feedback(G*Ko,Gc),
Gcl=feedback(GL,1),
step(tf(Gcl.num{1},Gcl.den{1}),10);
hold on;
Gold=feedback(G,1);
step(Gold,10)