| 查看: 2329 | 回复: 12 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
qinzhong6138至尊木虫 (著名写手)
灰常灰常笨的小笨蛋
|
[求助]
请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数 已有1人参与
|
|||
最近在做一点研究，但是被几个微分方程难住了，已经头疼好几天了 。现在发帖子求助一下各位牛人，也请伸出援手，帮一下我啊。我的目的是用matlab求解得到k1--k8的值，一共是6个微分方程的8个系数，如图所示。请大家指点一下啊！我问别人要了一段代码，但是这段代码还不够完整，现将其贴上来（本人不会写  ）和文献上面的公式贴上来，原始的文献也在附件。我在这里先谢谢大家了啊！！！ 600金币表示一下心意，不够的话，可以加啊function KineticsEst1_int clear all clc global CAm CBm CCm CDm CEm t = [ ]; CAm = []’; CBm = []'; CCm = []'; CDm=[]'; CEm=[]'; beta0 = [ ]; tspan = [ ]; CA0 = 0.11; CB0 = 0; CC0 = 0; CD0= 0.03; CE0= 0.106; C0=[0.31 0 0 0.02 0.10]; [beta,resnorm,resid,exitflag,output,lambda,jacobian] = ... lsqnonlin(@OptObjFunc,beta0,[0 0 0 0 0 0],[],[],tspan,[CA0 CB0 CC0 CD0 CE0]) ci = nlparci(beta,resid,jacobian) [t4plot CA4plot] = ode45(@KineticsEqs,[tspan(1) tspan(end)],[CA0 CB0 CC0 CD0 CE0],[],beta); plot(tspan,CAm,'bo',t4plot,CA4plot(:,1),'k-') hold on plot(tspan,CBm,'bo',t4plot,CA4plot(:,2),'k-') hold on plot(tspan,CCm,'bo',t4plot,CA4plot(:,3),'k-') legend('Exp','Model') xlabel('time, h') ylabel('concentration, mol/L') legend('Exp','Model') fprintf('Estimated Parameters:\n') fprintf('\tk1 = %.20f ¡À %.20f\n', beta(1),ci(1,2)-beta(1)) fprintf('\tk2 = %.20f ¡À %.20f\n', beta(2),ci(2,2)-beta(2)) fprintf('\tk3 = %.20f ¡À %.20f\n',beta(3),ci(3,2)-beta(3)) fprintf('\tk4 = %.20f ¡À %.20f\n',beta(4),ci(4,2)-beta(4)) fprintf('\tk5 = %.20f ¡À %.20f\n',beta(5),ci(5,2)-beta(5)) fprintf('\tk6 = %.20f ¡À %.20f\n',beta(6),ci(6,2)-beta(6)) （不会写了） % ------------------------------------------------------------------  111.jpg  222.jpg  333.jpg  444.jpg |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Product_sampling_during_transient_continuou.pdf
2014-07-12 22:08:43, 2.11 M
» 猜你喜欢
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有261人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有17人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
此贴已删除
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教一个不可压流动投影法的3阶Runge-Kutta显式时间推进的问题
已经有1人回复
- 请问Runge-Kutta算法的步长
已经有3人回复
- 【资源】Runge-Kutta解ODE的一些程序及应用（from MathWorks）
已经有3人回复
- 【求助】如何用Runge-Kutta迭代求解二阶常微分方程组【已解决】
已经有9人回复
qinzhong6138
至尊木虫 (著名写手)
灰常灰常笨的小笨蛋
- 应助: 9 (幼儿园)
- 金币: 22749.4
- 散金: 3540
- 红花: 14
- 帖子: 2194
- 在线: 2261.1小时
- 虫号: 585846
- 注册: 2008-08-08
- 专业: 功能与智能高分子
3楼2014-07-13 11:41:48
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68131.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
qinzhong6138: 金币+20, ★★★很有帮助, 先谢谢版主，等待你施以援手！ 2014-07-13 11:45:16
qinzhong6138: 金币+80, ★★★很有帮助, 谢谢版主的帮助，我从别的地方找到相应解决办法了，该帖子关了吧，谢谢啦 2014-07-16 14:52:15
感谢参与，应助指数 +1
qinzhong6138: 金币+20, ★★★很有帮助, 先谢谢版主，等待你施以援手！ 2014-07-13 11:45:16
qinzhong6138: 金币+80, ★★★很有帮助, 谢谢版主的帮助，我从别的地方找到相应解决办法了，该帖子关了吧，谢谢啦 2014-07-16 14:52:15
|
你是要拟合常微分方程组的参数吧。 方便的话，把数据也一起发上来。 如果只是要程序代码的话，参见：http://muchong.com/bbs/viewthread.php?tid=7575773&authorid=1122189 |
2楼2014-07-12 23:13:57
qinzhong6138
至尊木虫 (著名写手)
灰常灰常笨的小笨蛋
- 应助: 9 (幼儿园)
- 金币: 22749.4
- 散金: 3540
- 红花: 14
- 帖子: 2194
- 在线: 2261.1小时
- 虫号: 585846
- 注册: 2008-08-08
- 专业: 功能与智能高分子
|
不懂怎么传附件，  ，只好用网盘上传了，谢谢月版主，D专家有空的话，也麻烦帮助一下啊。@月只蓝，@ dingd链接: http://pan.baidu.com/s/1dD4Bf7V 密码: 6mcn |
4楼2014-07-13 11:47:26
qinzhong6138
至尊木虫 (著名写手)
灰常灰常笨的小笨蛋
- 应助: 9 (幼儿园)
- 金币: 22749.4
- 散金: 3540
- 红花: 14
- 帖子: 2194
- 在线: 2261.1小时
- 虫号: 585846
- 注册: 2008-08-08
- 专业: 功能与智能高分子
5楼2014-07-13 11:50:57
50