ode45用法案例

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=…