%计算不同的ξ 值所对应的振型 for j=1:10; eta(j)=0.1*j; wd(j)=wn*sqrt(1-eta(j)^2); %求振幅 A
a=sqrt((wn*x0*eta(j)+v0)^2+(x0*wd(j))^2)/wd(j); %求相位角
phi=atan2(wd(j)*x0,v0+eta(j)*wn*x0); %设定自变量数组 t t=0:tf/1000:tf; %求过渡过程
x(j,:)=a*exp(-eta(j)*wn*t).*sin(wd(j)*t+phi); end
%在同一个图形窗口中绘制不同的ξ 值所对应的振型 plot(t,x(1,:),t,x(2,:),t,x(3,:),t,x(4,:),... t,x(5,:),t,x(6,:),t,x(7,:),t,x(8,:),... t,x(9,:),t,x(10,:)) grid on
%新建一个图形窗口，绘制三维网格图 figure mesh(x)
