ode45解动力学方程

%参数定义m1=1;m2=2;m3=3;k1=4;k2=5;k3=6;k4=7;P0=8;w=9;%微分方程A=[k1+k2-k10;-k2k2+k3-k3;0-k3k3+k4];dx=@(t,x)[x(4:6);(-A*x(1:3)+[P0*sin(w*t);0;0])./[m1;m2;m3]];%初始条件（x4~x6为x1'~x3'）x0=zeros(6,1);%求解，直接绘图ode45(dx,[05],x0)lstr=arrayfun(@(i){sprintf('x%i',i)},1:6);legend(lstr{:})legend(lstr{:},'Location','S','Orientation','horizontal')