| 查看: 2311 | 回复: 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
» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有2人回复
-
070200求调剂，一志愿某211，288分
已经有18人回复
-
物理学I论文润色/翻译怎么收费?
已经有171人回复
-
法国博士后职位
已经有0人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
-
**
已经有1人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有0人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有2人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教一个不可压流动投影法的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
- 专业: 功能与智能高分子
9楼2014-07-13 14:28:13
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.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
- 专业: 功能与智能高分子
3楼2014-07-13 11:41:48
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
