zeta = [0.01,0.02,0.05,0.5,1,1.01,1.05,1.5,2.0];
omega = 2*pi;
for i=1:9
[t,x] = ode45(@(t,x) ff1(t,x,i),[0 10],[0.2;0]);
plot(t,x(:,1),'DisplayName',['\zeta=' num2str(zeta(i))]);
hold on
legend("Location","east")
ylabel('x','Interpreter','latex','FontSize',12)
xlabel('t','Interpreter','latex','FontSize',12);
end
function dxdt = ff1(t,x,i)
zeta = [0.01,0.02,0.05,0.5,1,1.01,1.05,1.5,2.0];
omega = 2*pi;
dxdt = [x(2); -2*zeta(i)*omega*x(2)-omega*omega*x(1)];
end
二阶变一阶:dot_x1 = x2; dot_x2=…