| 查看: 724 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 hilihuo 的 3 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
hilihuo至尊木虫 (著名写手)
|
[交流]
【求助】用MATLAB解复杂的微分方程 (急!!!)
|
||
|
请问如何用MATLAB解复杂的微分方程 (急!!!)谢谢写出程序和方法!!谢谢!!! 要解的方程如图: [ Last edited by fegg7502 on 2009-12-7 at 22:53 ] |
回复此楼» 猜你喜欢
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有298人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有22人回复
-
考博自荐
已经有4人回复
-
山东大学第二批博士研究生招生
已经有0人回复
-
中国科学院东莞材料科学与技术研究所-2026年博士招生-吴昊研究员-磁学与自旋电子学
已经有0人回复
-
《电磁学》教材推荐
已经有1人回复
-
【急招】合肥工大核聚变材料计算方向2026级工程博士生
已经有4人回复
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 17360.9
- 散金: 67
- 红花: 6
- 帖子: 2816
- 在线: 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
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 17360.9
- 散金: 67
- 红花: 6
- 帖子: 2816
- 在线: 717.1小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
2楼2009-11-29 22:47:08
hilihuo
至尊木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 17360.9
- 散金: 67
- 红花: 6
- 帖子: 2816
- 在线: 717.1小时
- 虫号: 746329
- 注册: 2009-04-11
- 性别: GG
- 专业: 污染控制化学
3楼2009-12-01 16:59:25
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