24小时热门版块排行榜    

查看: 2241  |  回复: 10

whqs8426212

铜虫 (正式写手)

[求助] 求助 1stopt低版本 方程的拟合(含若干未知参数) 已有1人参与

公式如图;
**********************************
******************************88
一组R随时间变化的实验数据
t        0     0.1    0.2   0.3   0.4      0.5
[R]    0.1   0.09 0.08 0.07 0.069    0.66
[R0]=0.1
[O30]=6.25*10^(-5);   
***************************
***************************
求未知参数:
k   ,kHo,  KO3,  [HO]/[O3]

求助 1stopt低版本  方程的拟合(含若干未知参数)
88.png
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
先把公式用简单的字母简化下,别用化学符号,更容易理解些。
2楼2014-10-05 20:16:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqs8426212

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by dingd at 2014-10-05 20:16:08
先把公式用简单的字母简化下,别用化学符号,更容易理解些。

****************实验数据*************
t        0     0.1    0.2   0.3   0.4      0.5
C(1)    0.1 ,0.09 ,0.08, 0.06,0.05, 0.047

t        0     0.1    0.2   0.3   0.4      0.5
C(2) 0.01,0.0009,0.00008,0.000007,0.000000069,0.00000000066

初值就是t=0时,C(1),C(2)的数值
***********函数*************************

fuction:
dC(1)=-k(1)*t;
dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1)
3楼2014-10-06 00:51:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

第二个公式:dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1)
完全可简化等同于:dC(2)=-k(2)*C(2)*C(1)

否则会出现过拟合现象,参数组值不唯一。

1stOpt求解:
CODE:
Variable t,C1,C2;
ODEFunction C1'=-k1*t;
                 C2'=-k2*C2*C1;
Data;
t=0,0.1,0.2,0.3,0.4,0.5;
C1=0.1,0.09,0.08,0.06,0.05,0.047;
C2=0.01,0.0009,0.00008,0.000007,0.000000069,0.00000000066;

C1和C2值数量上相差太多了,相对合理的结果参考如下;
均方差(RMSE):0.00797587120432916
残差平方和(SSE):0.00063614521468047
相关系数(R): 0.960512267688022
相关系数之平方(R^2): 0.922583816379187
确定系数(DC): 0.776373194295518
F统计(F-Statistic): 6.71167688419456

参数                  最佳估算
--------------------        -------------
k1        0.526046976304358
k2        243.193548940279
求助 1stopt低版本  方程的拟合(含若干未知参数)-1
c119.jpg


求助 1stopt低版本  方程的拟合(含若干未知参数)-2
c120.jpg

4楼2014-10-06 20:36:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqs8426212

铜虫 (正式写手)

引用回帖:
4楼: Originally posted by dingd at 2014-10-06 20:36:44
第二个公式:dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1)
完全可简化等同于:dC(2)=-k(2)*C(2)*C(1)

否则会出现过拟合现象,参数组值不唯一。

1stOpt求解:

Variable t,C1,C2;
ODEFunction C1'=-k1*t;
       ...

什么是过拟合现象,是把特别小的参数忽略吗

[ 发自小木虫客户端 ]
5楼2014-10-07 16:28:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

拟合如下方程:y=a+b*c*x
无法得到稳定解,有无穷多组b和c值,原方程等同于:y=a+b*x
6楼2014-10-07 20:24:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

whqs8426212

铜虫 (正式写手)

引用回帖:
6楼: Originally posted by dingd at 2014-10-07 20:24:37
拟合如下方程:y=a+b*c*x
无法得到稳定解,有无穷多组b和c值,原方程等同于:y=a+b*x

我将Code放入版本的1stOPT,结果报错'数据长度不够‘
code:*******************************************

Parameters    k1,k2;
Variable     t,C1,C2;
ODEFunction   C1'=-k1*t;  C2'=-k2*C2*C1;
data t=0,0.1,0.2,0.3,0.4,0.5;
C1=0.1,0.09,0.08,0.06,0.05,0.047;
C2=0.01,0.0009,0.00008,0.000007,0.000000069,0.00000000066;
7楼2014-10-08 12:36:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

data后面少“;”号。还有如果用的是网上流行的破解1.5版的话,是不能进行微分方程拟合计算的。
8楼2014-10-08 12:59:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水中的梦88

木虫 (正式写手)

引用回帖:
8楼: Originally posted by dingd at 2014-10-08 12:59:11
data后面少“;”号。还有如果用的是网上流行的破解1.5版的话,是不能进行微分方程拟合计算的。

请问那个版本以上的可以做,谢谢!
9楼2015-10-19 20:59:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

引用回帖:
9楼: Originally posted by 水中的梦88 at 2015-10-19 20:59:37
请问那个版本以上的可以做,谢谢!...

5.0以上吧
10楼2015-10-19 21:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 whqs8426212 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 跳槽后在研项目怎么办? +4 简单化xn 2026-08-22 8/400 2026-08-23 07:20 by ilovexiaobin
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 QTy3jDtz1uLt 2026-08-21 11/550 2026-08-23 02:55 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见