怎么用MATLAB求解二阶微分方程组
方程如下m1a1=-c1(v1-v0)-k1(z1-z0)+c2(v2-v1)+k2(z2-z1)-Fsin(wt);m2a2=-c2(v2-v1)-k2(z2-z1)+...
方程如下
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的数值解？ 展开
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（财富值+成长值）
1个回答
展开全部
首先z1z2z3的初值应该给出的， 否则无法求解，你少定解条件！！！我自己给了一个初始条件令他们都为0，取100个时间点为区间，然后进行了求解，结果如下：
你这个方程组好像结构动力学中的mck方程，结果类似正弦波，程序你可以私信我，建议你去再确定一下初始条件
追问
我编出来的程序有问题，是不是不能用dslove来解？
追答
是的，dslove应该是求解解析解的，求数值解要用ode系列的求解器
本回答被提问者和网友采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
系科仪器
2024-08-02 广告
科仪器致力于为微纳薄膜领域提供精益级测量及控制仪器，包括各种光谱椭偏、激光椭偏、反射式光谱等，从性能参数、使用体验、价格、产品可靠性及工艺拓展性等多个维度综合考量，助客户提高研发和生产效率，以及带给客户更好的使用体验。...
点击进入详情页
本回答由系科仪器提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
