ex6_14.m

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

% ex6_14
% 速度反馈控制系统的设计
% 希望系统最大超调量Mp<0.2
% 峰值时间Tp为1秒
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}));