| 查看: 2035 | 回复: 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为常量 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
能否退出参与的面上项目解除限项
已经有13人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
-
科研孤儿太难了
已经有19人回复
-
明天应该可查了！？
已经有4人回复
-
什么时候开奖？
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- mathematica 编程求助
已经有9人回复
- 如果编程语言是汽车……
已经有16人回复
- 这个复杂式子如何计算出求导结果
已经有5人回复
【答案】应助回帖
★ ★
余泽成(金币+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
mshwangg
至尊木虫 (正式写手)
- 程序强帖: 5
- 应助: 206 (大学生)
- 金币: 10702.8
- 散金: 100
- 红花: 19
- 帖子: 597
- 在线: 195.4小时
- 虫号: 576702
- 注册: 2008-06-21
- 专业: 物理学I
2楼2011-11-17 00:30:03
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
20