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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助1stopt拟合动力学参数
51/1返回列表
查看: 2219  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chen1215

新虫 (初入文坛)

[求助] 求助1stopt拟合动力学参数 已有1人参与

最近做反应动力学模型，需要拟合微分方程组的几个参数，之前用matlab，ode45求解，得到的数值解套用非线性最小二乘拟合，可是总也拟合不好。不知道是初值不合适还是模型有问题，但模型之前有文献报道过，有问题的可能性较小。
想求助用1stopt拟合试试，看看能不能得到比较好的效果~
我对1stopt 不太熟，试着写了下程序，不知道有没有问题，具体如下：
Parameter a,b,c,d,e,f,g;
Variable t,x,y,z,l,m,n;
Function x’=-a*x*n*(0.000384369/(1+d*x+e*y+f*z+g*m))；
y’=(a*x*n-b*y*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
z’=(b*y*n-c*z*n)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
l’=c*z*n*(0.000384369/(1+d*x+e*y+f*z+g*m))；
m’=(a*x+b*y+c*z)*(0.000384369/(1+d*x+e*y+f*z+g*m))；
n’= -m’；
RowData；
t=0 5 10 20 30 40 60;
x=0.930682749        0.598029237        0.541233252        0.502127725        0.46364134        0.442517119        0.404660654;
y=0         0.192347906        0.219628094        0.23819892        0.251976191        0.252181432        0.264979117;
z=0         0.044174154        0.058131509        0.066640197        0.073187605        0.084552378        0.093852475;
l=0         0.096131452        0.111689894        0.123715906        0.141877613        0.15143182        0.167190503;
m=0          0.569090571        0.670960794        0.742627034        0.823984239        0.875581648        0.954255576;
n=4.591368228        4.022277657        3.920407434        3.848741194        3.767383988        3.71578658        3.637112651;
就是已知x,y,z,l,m,n 随时间变化的值，想拟合出模型中的参数 a,b,c,d,e,f,g(这些参数的值都是大于等于0）
之前有下过低版本的1stopt，结果点运行没反应，不知道是不是版本问题还是代码写的有问题。。

求各路牛人帮忙用高版本运行下，非常非常感谢~~
回复此楼

» 本帖@通知

@dingd

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-03-31 17:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
如果都要求大于0：

均方差(RMSE):0.291509496670823
残差平方和(SSE):3.05920031937395
相关系数(R): 0.41204372086458
相关系数之平方(R^2): 0.169780027903928
确定系数(DC): -32.2049580782807
F统计(F-Statistic): -0.0103799073026604

参数                  最佳估算
--------------------        -------------
a        2379421492.52982
b        5477871570.00215
c        4827348354.29146
d        0.000157698484986776
e        28.2125368960529
f        9.53786176032263E-8
g        333524832.481022

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

文件: 数据文件-1
No        t        目标 x        计算 x        目标 y        计算 y        目标 z        计算 z        目标 l        计算 l        目标 m        计算 m        目标 n        计算 n
1        5        0.598029237        0.461228912051106        0.192347906        0.212216074803243        0.044174154        0.154698496352178        0.096131452        0.102539265793472        0.569090571        0.18430565779089        4.022277657        4.40706257020911
2        10        0.541233252        0.352387819112766        0.219628094        0.194209366200932        0.058131509        0.186948110939243        0.111689894        0.19713745274706        0.670960794        0.259898922792036        3.920407434        4.33146930520796
3        20        0.502127725        0.241106037906148        0.23819892        0.153262848683966        0.066640197        0.191130610932665        0.123715906        0.34518325147722        0.742627034        0.355979862816147        3.848741194        4.23538836518385
4        30        0.46364134        0.179231653280274        0.251976191        0.121526115166401        0.073187605        0.174250824742148        0.141877613        0.455674155811178        0.823984239        0.419251043228081        3.767383988        4.17211718477191
5        40        0.442517119        0.138636909389627        0.252181432        0.0975440616385669        0.084552378        0.153342396668236        0.15143182        0.541159381303569        0.875581648        0.465201935622196        3.71578658        4.1261662923778
6        60        0.404660654        0.0883552573120677        0.264979117        0.0646893622260779        0.093852475        0.114305811509068        0.167190503        0.663332317952786        0.954255576        0.527617939775541        3.637112651        4.06375028822446
赞一下 回复此楼
3楼2014-04-03 14:01:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

终极爱好者

铜虫 (初入文坛)
想学 matlab 不知道 改怎么获取资料
赞一下 回复此楼
2楼2014-04-01 09:45:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chen1215

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by dingd at 2014-04-03 14:01:19
如果都要求大于0：

均方差(RMSE):0.291509496670823
残差平方和(SSE):3.05920031937395
相关系数(R): 0.41204372086458
相关系数之平方(R^2): 0.169780027903928
确定系数(DC): -32.2049580782807
F统计(F- ...

真是太感谢了。不过这个拟合效果确实挺差的，不知道是不是模型不太合适。。
你能帮我再试试这组数据么？谢谢~
t=0        2        5        10        15        20        30
x=0.97929461        0.848102002        0.791547773        0.742120992        0.665706681        0.627442016        0.567283612
y=0        0.111151416        0.154799598        0.178808192        0.203933465        0.230873341        0.255998614
z=0        0.029197641        0.031321393        0.058558032        0.078165633        0.1017788        0.122183142
l=0        0.000965082        0.003370178        0.008463798        0.012727551        0.019200453        0.030677511
m=0        0.142077349        0.222319921        0.295346436        0.454731222        0.492032299        0.601852621
n=1.929382737        1.787305387        1.707062816        1.634036301        1.474651515        1.437350438        1.327530115
赞一下 回复此楼
4楼2014-04-03 15:01:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蜗牛虫

银虫 (正式写手)
引用回帖:
3楼: Originally posted by dingd at 2014-04-03 14:01:19
如果都要求大于0：

均方差(RMSE):0.291509496670823
残差平方和(SSE):3.05920031937395
相关系数(R): 0.41204372086458
相关系数之平方(R^2): 0.169780027903928
确定系数(DC): -32.2049580782807
F统计(F- ...

您好，您是怎么做到的？

肯求授以渔：）


非常感谢
赞一下 回复此楼
多思考勤记录嗯加油
5楼2015-12-07 10:55:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] 我面上完蛋了 +9 且听虎啸 2026-08-20 10/500 2026-08-24 20:08 by 苦寒来香
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭