| 查看: 2041 | 回复: 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为常量 |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有24人回复
-
号外号外
已经有8人回复
-
filecode，4个jtjc了
已经有17人回复
-
人气不行了
已经有9人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有5人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- mathematica 编程求助
已经有9人回复
- 如果编程语言是汽车……
已经有16人回复
- 这个复杂式子如何计算出求导结果
已经有5人回复
7楼2011-11-23 18:58:36
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
20