MATLAB的ode45函数解常微分方程的用法！！

Question

这个怎么用MATLAB的ode函数来弄呀，我看其他教程都是只有一个方程，；两个方程怎么弄呀，请教一下！

Answer 1

你好，两个方程一样求解：
代码供参考，该注释的都注释了

% 设置物块质量和弹簧刚度
m1 = 1;
m2 = 1;
k1 = 1;
k2 = 1;
k3 = 1;
% 设x(1)对应物块1的位移;x(2)对应物块1的速度
% 设x(3)对应物块2的位移;x(4)对应物块2的速度
odefun = @(t, x) [ x(2); % 物块1的速度对时间积分得到位移
    (k2*x(3) - (k1+k2)*x(1)) / m1; % 物块1的加速度对时间积分得到速度
    x(4); % 物块2的速度对时间积分得到位移
    ( -(k2+k3)*x(3) + k2*x(1)) / m2]; % 物块2的加速度对时间积分得到速度
x0 = [0.1; 0; 0; 1]; % 物块1初始位移0.1，初始速度0，物块2初始位移0，初始速度1
tspan = [0, 10]; % 运行10s
[t,x] = ode45(odefun, tspan, x0);%ode45求解
% 画图
plot(t,x(:,1),'r-');hold on
plot(t,x(:,2),'b-')
plot(t,x(:,3),'r--')
hold on
plot(t,x(:,4),'b--')
legend('物块1位移','物块1速度','物块2位移','物块2速度')

有帮助望题主给个采纳支持一下