p110_3.m

来自「matlab仿真的文件包然后写出其具体功能(至少要20个字)。尽量不要让站长把时」· M 代码 · 共 10 行

wns=[0.1:0.1:1]; z=0.707; yy=[]; t=0:.1:12;
for i=1:length(wns)
   wn=wns(i); 
   y=wn*exp(-z*wn*t).*sin(wn*sqrt(1-z^2)*t)/sqrt(1-z^2);
   yy=[yy;y];
end
plot(t,yy)

mesh(wns,t,yy'); set(gca,'XDir','reverse'); axis([0,1,0,12,0,0.5])