| 查看: 703 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hilihuo 的 3 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hilihuo至尊木虫 (著名写手)
|
[交流]
【求助】用MATLAB解复杂的微分方程 (急!!!)
|
||
|
请问如何用MATLAB解复杂的微分方程 (急!!!)谢谢写出程序和方法!!谢谢!!! 要解的方程如图: [ Last edited by fegg7502 on 2009-12-7 at 22:53 ] |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有133人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 16796.9
- 散金: 67
- 红花: 6
- 帖子: 2784
- 在线: 717.1小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
3楼2009-12-01 16:59:25
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 16796.9
- 散金: 67
- 红花: 6
- 帖子: 2784
- 在线: 717.1小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
2楼2009-11-29 22:47:08
change0618
铁杆木虫 (著名写手)
方丈大师
- 计算强帖: 1
- 应助: 44 (小学生)
- 金币: 17724.5
- 红花: 17
- 帖子: 2413
- 在线: 546.7小时
- 虫号: 496517
- 注册: 2008-01-19
- 专业: 化学反应工程
★ ★ ★ ★
mingdong(金币+2,VIP+0):感谢交流！ 12-2 22:06
hilihuo(金币+2,VIP+0):谢谢热心帮助！！十分感激！！ 12-4 21:30
mingdong(金币+2,VIP+0):感谢交流！ 12-2 22:06
hilihuo(金币+2,VIP+0):谢谢热心帮助！！十分感激！！ 12-4 21:30
|
function test clear clc global k1 k2 k3 k7 lamda1 lamda2 k1 = 1; k2 = 1; k3 = 1; k7 = 1; lamda1 = 1; lamda2 = 1; [x,y]=ode23s(@myfun,[0,10],1); [x y] plot(x,y) function f=myfun(t,x) global k1 k2 k3 k7 lamda1 lamda2 f=-1/(k1*exp((lamda1+lamda2)*x))-(k3*exp(lamda1*x)+k7*exp(lamda2*x))/(k1*k2*exp((lamda1+lamda2)*x))/x; |
4楼2009-12-02 14:37:20
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 16796.9
- 散金: 67
- 红花: 6
- 帖子: 2784
- 在线: 717.1小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
谢谢。但仍然算不出来。请问你可以算出来么？
★
fegg7502(金币+1,VIP+0):鼓励交流2 12-19 07:48
fegg7502(金币+1,VIP+0):鼓励交流2 12-19 07:48
|
运行结果如下： >> function test clear clc global k1 k2 k3 k7 lamda1 lamda2 k1=1; k2=1; k3=1; k7=1; lamda1=1; lamda2=1; [x,y]=ode23s(@myfun,[0,10],1); [x y] plot(x,y) function f=myfun(t,x) global k1 k2 k3 k7 lamda1 lamda2 f=-1/(k1*exp((lamda1+lamda2)*x))-(k3*exp(lamda1*x)+k7*exp(lamda2*x))/(k1*k2*exp((lamda1+lamda2)*x))/x; ??? function test | Error: Function definitions are not permitted at the prompt or in scripts. |
5楼2009-12-02 20:13:37
