24小时热门版块排行榜    

查看: 1539  |  回复: 6

amenboss

铁虫 (初入文坛)

[交流] 求高手相助!matlab求反应动力学模型参数 已有2人参与

我通过实验有了x-t的实验数据,我现在想用matlab的函数求解模型中的参数,我的思路大概步骤是:首先假设参数的初值,并给其他参数赋值(已知),求出dx/dt,然后调用ODE函数求出Xsim的值,最后求所有参数点的Xsim和Xexp的差平方,以这个差平方和最小为优化目标求两个参数。
模型大致如图所示
求大神帮忙啊

求高手相助!matlab求反应动力学模型参数
1%D`%00ZJ3MRNB{N2)[N%HD.png
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
数据一起放上来看看。
2楼2014-12-14 13:54:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
2楼的大神是我们版块求解你提出问题的最强者,给出数据吧。
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
3楼2014-12-14 14:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amenboss

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by dingd at 2014-12-14 13:54:09
数据一起放上来看看。

求大神帮忙!
我的数据是如下:
t        x
0        2054.794521
300        1102.499315
600        695.659589
900        406.4431507
1200        198.1123288
1500        121.0650685
2100        72.21986301
2700        60.43630137
3300        44.2869863
3900        31.69109589
4500        26.82808219
5100        22.52260274


m,n,t,s的初值大致为
314.5833333,  0.328315976,  110.2308769,  2064.86589

初值可能会有出入,我也是大致估算的,没有准确数据。毕业论文遇到这个不熟悉的领域真是太难受了,希望大神能帮忙,感激涕零啊!!
4楼2014-12-15 16:42:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
公式中t是积分区间不是参数吧?
CODE:
Variable t,x;
ODEFunction x'=-(m+n*x)/(t+1/s/x^2);
Data;
t        x
0        2054.794521
300        1102.499315
600        695.659589
900        406.4431507
1200        198.1123288
1500        121.0650685
2100        72.21986301
2700        60.43630137
3300        44.2869863
3900        31.69109589
4500        26.82808219
5100        22.52260274

均方差(RMSE):13.8458076348104
残差平方和(SSE):2108.77027966193
相关系数(R): 0.999300935915745
相关系数之平方(R^2): 0.998602360522083
确定系数(DC): 0.998282949539454
F统计(F-Statistic): 2329.57626576002

参数                  最佳估算
--------------------        -------------
m        1228.02989306005
n        -0.52012529093965
s        1.39532314291721E-8

====== 结果输出 ======

文件: 数据文件-1
No        t        目标 x        计算 x
1        300        1102.499315        1113.41240644706
2        600        695.659589        680.535799025786
3        900        406.4431507        396.054156426377
4        1200        198.1123288        222.256063923777
5        1500        121.0650685        131.838213629664
6        2100        72.21986301        63.2589003819102
7        2700        60.43630137        39.8381649979752
8        3300        44.2869863        28.7549628853483
9        3900        31.69109589        22.4061726126003
10        4500        26.82808219        18.3200101730951
11        5100        22.52260274        15.4790969702342
求高手相助!matlab求反应动力学模型参数-1
c1.jpg

5楼2014-12-15 16:58:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

amenboss

铁虫 (初入文坛)

引用回帖:
5楼: Originally posted by dingd at 2014-12-15 16:58:26
公式中t是积分区间不是参数吧?

Variable t,x;
ODEFunction x'=-(m+n*x)/(t+1/s/x^2);
Data;
t        x
0        2054.794521
300        1102.499315
600        695.659589
900        406.4431507
1 ...

谢谢大神!
忘了说了,t不是变量,是参数。
这个拟合是用什么软件做的啊?
6楼2014-12-15 23:03:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你的定义:Variable t,x;

t怎么又成了待求参数?如果是参数你的公式就过拟合了,可以化简去掉一个。
7楼2014-12-16 09:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 amenboss 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[基金申请] 微信指数没变化,科研之友没阅读 +15 wangze12014 2026-07-28 19/950 2026-08-02 20:05 by 蔡棒棒菂
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
信息提示
请填处理意见