| 查看: 888 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hilihuo 的 3 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hilihuo至尊木虫 (著名写手)
|
[交流]
【求助】用MATLAB解复杂的微分方程 (急!!!)
|
||
|
请问如何用MATLAB解复杂的微分方程 (急!!!)谢谢写出程序和方法!!谢谢!!! 要解的方程如图: [ Last edited by fegg7502 on 2009-12-7 at 22:53 ] |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有233人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
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 (幼儿园)
- 金币: 18155.9
- 散金: 67
- 红花: 6
- 帖子: 2958
- 在线: 717.3小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
2楼2009-11-29 22:47:08
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 18155.9
- 散金: 67
- 红花: 6
- 帖子: 2958
- 在线: 717.3小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
3楼2009-12-01 16:59:25
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 18155.9
- 散金: 67
- 红花: 6
- 帖子: 2958
- 在线: 717.3小时
- 虫号: 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
