fig3_37b.m

来自「包含了控制理论的各种仿真程序」· M 代码 · 共 33 行

%  Figure 3.37b      Feedback Control of Dynamic Systems, 4e
%                        Franklin, Powell, Emami
%

clf;
num=1;
a=5;

zeta =1;
den2=[1/(zeta*a) 1];
den1=[1 2*zeta 1];
den=conv(den1,den2);
t=0:.1:8;
y1=step(num,den,t);

zeta =.7;
den2=[1/(zeta*a) 1];
den1=[1 2*zeta 1];
den=conv(den1,den2);
y2=step(num,den,t);

zeta =.5;
den2=[1/(zeta*a) 1];
den1=[1 2*zeta 1];
den=conv(den1,den2);
y3=step(num,den,t);

axis([0 5 .1 .9])
plot(t,y1,'-',t,y2,'--',t,y3,'-.'),grid
title('Fig. 3.37b Step response with extra pole, \alpha= 5')
xlabel('\omegan t  ')
ylabel('y(t)')