应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 请教一下各位大神怎么求解 Runge–Kutta method的微分方程的各个系数
51/1返回列表
查看: 2140  |  回复: 12
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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


-1
222.jpg


-2
333.jpg


-3
444.jpg
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : Product_sampling_during_transient_continuou.pdf
    • 2014-07-12 22:08:43, 2.11 M

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
加油
1楼2014-07-12 22:05:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by qinzhong6138 at 2014-07-13 11:41:48

先谢谢版主了，我不懂matlab，所以对matlab代码也不是很懂。了解到1stOpt对解微分方程也很好，ps:dingd是这方面的行家的，多请教了。昨晚花了几个小时，大致了解一下这个软件，因此，自己编写了一些1stOpt ...

1stopt软件要高版本的才能做常微分方程的拟合，而且这款软件的拟合效果常比普通编写的MATLAB程序强。你可以把完整的1stopt的程序写好，让dingd专家跑一下。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2014-07-13 13:39:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-07-12 23:13:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinzhong6138

至尊木虫 (著名写手)

灰常灰常笨的小笨蛋
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-07-12 23:13:57
你是要拟合常微分方程组的参数吧。
方便的话，把数据也一起发上来。
如果只是要程序代码的话，参见：http://muchong.com/bbs/viewthread.php?tid=7575773&authorid=1122189


先谢谢版主了，我不懂matlab，所以对matlab代码也不是很懂。了解到1stOpt对解微分方程也很好，ps:dingd是这方面的行家的，多请教了。昨晚花了几个小时，大致了解一下这个软件，因此，自己编写了一些1stOpt格式的程序和数据，不知道您对这个了解吗？以后我会经常遇到类似的非线性拟合的数据，如果1stOpt真的很方便的话，我打算买一套，支持国产，于己于人都方便，不过这个软件价格有些肉疼啊，年轻老师，收入低。版主，您方便的话，能麻烦您用matlb和1stOpt都能帮我算一下吗，数据马上上传。@月只蓝，@ dingd
赞一下 回复此楼
加油
3楼2014-07-13 11:41:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qinzhong6138

至尊木虫 (著名写手)

灰常灰常笨的小笨蛋
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-07-12 23:13:57
你是要拟合常微分方程组的参数吧。
方便的话，把数据也一起发上来。
如果只是要程序代码的话，参见：http://muchong.com/bbs/viewthread.php?tid=7575773&authorid=1122189

不懂怎么传附件，，只好用网盘上传了，谢谢月版主，D专家有空的话，也麻烦帮助一下啊。@月只蓝，@ dingd
链接: http://pan.baidu.com/s/1dD4Bf7V 密码: 6mcn
赞一下 回复此楼
加油
4楼2014-07-13 11:47:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭