| 查看: 1793 | 回复: 6 | ||
liuhaizhen木虫 (初入文坛)
|
[求助]
求mathematica编程
|
|
dy(1)/dt=(a+b+c)*y(1)+g*y(2)+b*y(3) dy(2)/dt=-(a+b)*y(3) dy(3)/dt=(a-b-c)*y(3)-g*y(1)+2*g*y(2)*y(1) 其中a,b,c,g为常量 |
回复此楼» 猜你喜欢
自荐读博
已经有7人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有6人回复
-
求个博导看看
已经有18人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- mathematica 编程求助
已经有9人回复
- 如果编程语言是汽车……
已经有16人回复
- 这个复杂式子如何计算出求导结果
已经有5人回复
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
2楼2011-11-17 00:30:03
【答案】应助回帖
★ ★
余泽成(金币+2): 谢谢参与应助！ 2011-11-19 14:28:38
liuhaizhen(金币+10): 感谢你的应助 2011-11-28 10:00:33
liuhaizhen(金币+5): 优秀 2011-12-07 10:16:41
余泽成(金币+2): 谢谢参与应助！ 2011-11-19 14:28:38
liuhaizhen(金币+10): 感谢你的应助 2011-11-28 10:00:33
liuhaizhen(金币+5): 优秀 2011-12-07 10:16:41
|
下面是程序代码，语法肯定没错，只要你再给定初始条件，就可以得到结果。否则解出来的是通解。 eq1 = y1'[t] == (a + b + c)*y1[t] + g*y2[t] + b*y3[t] eq2 = y2'[t] == -(a + b)*y3[t] eq3 = y3'[t] == (a - b - c)*y3[t] - g*y1[t] + 2*g*y2[t]*y1[t] DSolve[{eq1, eq2, eq3}, {y1[t], y2[t], y3[t]}, t] |
3楼2011-11-17 09:46:25
4楼2011-11-17 09:48:37
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
5楼2011-11-17 16:25:28
6楼2011-11-17 22:02:13
7楼2011-11-23 18:58:36