damp1.m

《Matlab控制系统计算机辅助设计》随书程序

y0=0.15;wn=sqrt(2);
zeta1=3/(2*sqrt(2));
zeta2=1/(2*sqrt(2));
t=[0:0.1:10];
%
t1=acos(zeta1)*ones(1,length(t));
t2=acos(zeta2)*ones(1,length(t));
c1=(y0/sqrt(1-zeta1^2));
c2=(y0/sqrt(1-zeta2^2));
y1=c1*exp(-zeta1*wn*t).*sin(wn*sqrt(1-zeta1^2)*t+t1);
y2=c2*exp(-zeta1*wn*t).*sin(wn*sqrt(1-zeta2^2)*t+t2);
%
bu=c2*exp(-zeta2*wn*t);bl=-bu;
%
plot(t,y1,'-',t,y2,'-',t,bu,'--',t,bl,'--'),grid;
xlabel('Time[sec]'),ylabel('y(t)Displacement[m]');
text(0.2,0.85,['overdamped zeta1=',num2str(zeta1),],'sc');
text(0.2,0.80,['underdamped zeta2=',num2str(zeta2),],'sc');