ode45解动力学方程

来源：二三娱乐

% 参数定义

m1 = 1; m2 = 2; m3 = 3;

k1 = 4; k2 = 5; k3 = 6; k4 = 7;

P0 = 8; w = 9;

% 微分方程

A = [k1+k2 -k1 0; -k2 k2+k3 -k3; 0 -k3 k3+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,[0 5],x0)

lstr = arrayfun(@(i){sprintf('x%i',i)},1:6);

legend(lstr{:})

legend(lstr{:},'Location','S','Orientation','horizontal')

因篇幅问题不能全部显示，请点此查看更多更全内容

查看全文

违法及侵权请联系：TEL:199 1889 7713 E-MAIL:2724546146@qq.com

本站由北京市万商天勤律师事务所王兴未律师提供法律服务