24小时热门版块排行榜    

查看: 2011  |  回复: 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的回帖

水中的梦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的回帖
查看全部 11 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,0703化学310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-21 22:21 by peike
[考研] 354求调剂 +6 Tyoumou 2026-03-18 9/450 2026-03-21 20:47 by lbsjt
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 0703化学调剂 +4 妮妮ninicgb 2026-03-21 4/200 2026-03-21 18:39 by 学员8dgXkO
[考研] 生物学调剂 +3 Surekei 2026-03-21 3/150 2026-03-21 18:31 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 298求调剂 +4 上岸6666@ 2026-03-20 4/200 2026-03-21 17:14 by 学员8dgXkO
[基金申请] 学校已经提交到NSFC,还能修改吗? 40+4 babangida 2026-03-19 9/450 2026-03-21 16:12 by babangida
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 机械专硕299求调剂至材料 +3 kkcoco25 2026-03-16 4/200 2026-03-21 03:52 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 一志愿武汉理工材料工程专硕调剂 +9 Doleres 2026-03-19 9/450 2026-03-20 22:36 by JourneyLucky
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
信息提示
请填处理意见