p107.m

反馈控制系统分析与设计,虽然有点抽象

wn=1; zetas=[0:0.1:1,2,3,5]; yy=[]; t=0:.1:12;
for i=1:length(zetas)
   z=zetas(i);
   if z==0, y=1-cos(wn*t);
   elseif (z>0 & z<1), 
      wd=wn*sqrt(1-z^2); th=atan(sqrt(1-z^2)/z);
      y=1-exp(-z*wn*t).*sin(wd*t+th)/sqrt(1-z^2);
   elseif z==1, y=1-(1+wn*t).*exp(-wn*t);
   elseif z>1, 
      lam1=-z-sqrt(z^2-1); lam2=-z+sqrt(z^2-1); 
      y=1-0.5*wn*(exp(lam1*t)/lam1-exp(lam2*t)/lam2)/sqrt(z^2-1);
   end
   yy=[yy;y];
end
plot(t,yy)