Question

怎么用MATLAB求解二阶微分方程组

方程如下
m1a1=-c1(v1-v0)-k1(z1-z0)+c2(v2-v1)+k2(z2-z1)-Fsin(wt);
m2a2=-c2(v2-v1)-k2(z2-z1)+c3(v3-v2)+k3(z3-z2);
m3a3=-c3(v3-v2)-k1(z3-z2);
数据： m1=15;m2=8.8;m3=43.4;
k1=31000;k2=18000;k3=44130;
c1=830;c2=200;c3=1485;
z0=0;v0=0;z1z2z3初始一阶导数为0，v是z的一阶导数，a是z的二阶导数。

f=1000,w=2
怎么得到z1z2z3的数值解？

• 你的回答被采纳后将获得：
• 系统奖励15（财富值+成长值）+难题奖励30（财富值+成长值）

Answer 1

首先z1z2z3的初值应该给出的， 否则无法求解，你少定解条件！！！我自己给了一个初始条件令他们都为0，取100个时间点为区间，然后进行了求解，结果如下：

你这个方程组好像结构动力学中的mck方程，结果类似正弦波，程序你可以私信我，建议你去再确定一下初始条件

追问

我编出来的程序有问题，是不是不能用dslove来解？

追答

是的，dslove应该是求解解析解的，求数值解要用ode系列的求解器