| 查看: 1095 | 回复: 6 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
毛毛虫飞了新虫 (著名写手)
|
[求助]
求帮忙解微分方程组 已有1人参与
|
|||
微分方程组如下，由于自己数学不行，也不会用matlab，所以恳请大神帮忙解一下下面的微分方程组。其中x y z m为因变量，t为自变量，k1 k2 k3 k4均为常量。 @月只蓝 @beefly 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有139人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有1人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 大家帮忙看看这个偏微分方程组如何求解析解和数值解？希望大家给个思路，尤其解析解
已经有2人回复
- 四阶龙格库特法解微分方程组
已经有6人回复
- matlab----解微分方程组，求大神帮忙！！
已经有6人回复
- 请教前辈：该二元二阶微分方程组中的两个变量的表达式如何求呢？
已经有16人回复
- 求大神指导用matlab解偏微分方程组
已经有1人回复
- 微分方程的数值解求法
已经有3人回复
- MATLAB求解偏微分方程组
已经有5人回复
- 求会解微分方程组的大神指点下——————
已经有0人回复
- matlab解含有定积分的微分方程组。方程组如图所示，求大神帮忙！！！
已经有2人回复
- 四阶非线性偏微分方程组求解
已经有6人回复
- 求助，用matlab解微分方程组，希望高手能给予指点
已经有17人回复
- 求组！求解一个微分方程组的符号解，并画出图像！！！
已经有7人回复
- 帮忙用Mathematica求解一个微分方程组
已经有6人回复
- 周期性微分方程组求解
已经有0人回复
- 有熟悉偏微分方程粘性解的吗？烦请看一下，求帮忙
已经有0人回复
- 求助matlab解一个含有微分方程的方程组
已经有1人回复
- [求助]跪求高手帮忙解一个微分方程组
已经有5人回复
- 【求助】求高人帮我求解常系数微分方程组
已经有3人回复
- 【求助】急求各位大侠，关于用ode45解微分方程
已经有8人回复
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
4楼2016-10-20 23:41:23
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
代码： [x,y,z,m] = dsolve('Dx=-k1*x','Dy=k1*x-k2*y',... 'Dz=k2*y-k3*z','Dm=k3*z-k4*m',... 'x(0)=1,y(0)=0,z(0)=0,m(0)=0','t'); x = simplify(x) y = simplify(y) z = simplify(z) m = simplify(m) 结果： x=(k1*k2*k3*exp(-k2*t))/((k1 - k2)*(k2 - k3)*(k2 - k4)) - (k1*k2*k3*exp(-k1*t))/((k1 - k2)*(k1 - k3)*(k1 - k4)) - (k1*k2*k3*exp(-k3*t))/((k1 - k3)*(k2 - k3)*(k3 - k4)) + (k1*k2*k3*exp(-k4*t))/((k1 - k4)*(k2 - k4)*(k3 - k4)) y=exp(-k1*t) z=(k1*exp(-k2*t))/(k1 - k2) - (k1*exp(-k1*t))/(k1 - k2) m=(k1*k2*exp(-k1*t))/((k1 - k2)*(k1 - k3)) - (k1*k2*exp(-k2*t))/((k1 - k2)*(k2 - k3)) + (k1*k2*exp(-k3*t))/((k1 - k3)*(k2 - k3)) |
2楼2016-10-20 22:45:01
512095518
专家顾问 (正式写手)
-
专家经验: +52 - 应助: 118 (高中生)
- 金币: 7665.4
- 红花: 78
- 帖子: 712
- 在线: 186.1小时
- 虫号: 3195576
- 注册: 2014-05-10
- 性别: GG
- 专业: 大气环境与全球气候变化
- 管辖: 计算模拟
【答案】应助回帖
|
我仔细看了下我给你的代码回复，有一个地方要注意，就是解的顺序。 也就是dsolve得到的4个解的顺序，matlab可能是按m,x,y,z这样排列的。所以正确的解应该是： m=(k1*k2*k3*exp(-k2*t))/((k1 - k2)*(k2 - k3)*(k2 - k4)) - (k1*k2*k3*exp(-k1*t))/((k1 - k2)*(k1 - k3)*(k1 - k4)) - (k1*k2*k3*exp(-k3*t))/((k1 - k3)*(k2 - k3)*(k3 - k4)) + (k1*k2*k3*exp(-k4*t))/((k1 - k4)*(k2 - k4)*(k3 - k4)) x=exp(-k1*t) y=(k1*exp(-k2*t))/(k1 - k2) - (k1*exp(-k1*t))/(k1 - k2) z=(k1*k2*exp(-k1*t))/((k1 - k2)*(k1 - k3)) - (k1*k2*exp(-k2*t))/((k1 - k2)*(k2 - k3)) + (k1*k2*exp(-k3*t))/((k1 - k3)*(k2 - k3)) |
3楼2016-10-20 23:32:01
毛毛虫飞了
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1589.3
- 散金: 673
- 红花: 3
- 帖子: 1361
- 在线: 210.5小时
- 虫号: 3051043
- 注册: 2014-03-14
- 专业: 能源化工
5楼2016-10-21 01:19:24