怎么用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-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大（共振点）；正弦振动在任一瞬间只包含一种频率的振动，而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率，所以样品上的共振点会同时激发...
点击进入详情页
本回答由富港检测技术（东莞）有限公司_提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
