whqs8426212

銅蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 13.2
散金: 375
紅花: 7
帖子: 782
在線: 685.7小時(shí)
蟲號: 2572504
注冊: 2013-07-30
專業(yè): 數(shù)論

[求助] 求助 1stopt低版本方程的擬合（含若干未知參數(shù)）已有1人參與

公式如圖；
**********************************
******************************88
一組R隨時(shí)間變化的實(shí)驗(yàn)數(shù)據(jù)
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);
***************************
***************************
求未知參數(shù)：
k ，kHo, KO3, [HO]/[O3]

88.png

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有212人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有68人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:

1stopt擬合中出現(xiàn)的問題已經(jīng)有8人回復(fù)
隱函數(shù)曲線擬合已經(jīng)有8人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有6人回復(fù)
LM算法擬合問題已經(jīng)有9人回復(fù)
求助大神幫忙用1stopt擬合下模型參數(shù) 已經(jīng)有9人回復(fù)
1stopt自變量超多時(shí)是否有簡單語法？AND 自變量多數(shù)據(jù)量大的話1stopt都會(huì)算的特慢么已經(jīng)有7人回復(fù)
請教一下各位大神怎么求解 Runge–Kutta method的微分方程的各個(gè)系數(shù) 已經(jīng)有12人回復(fù)
求大神跑下1stOpt，一個(gè)sharemodel，謝謝！已經(jīng)有7人回復(fù)
求助--用1stopt擬合藍(lán)寶石襯底的GaN薄膜紅外反射光譜已經(jīng)有6人回復(fù)
1stopt解方程組已經(jīng)有12人回復(fù)
用1stOpt進(jìn)行復(fù)數(shù)曲線擬合時(shí)，結(jié)果出錯(cuò)。急求教！！！已經(jīng)有5人回復(fù)
求高版本1stopt擬合，已經(jīng)有11人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)
用1stopt擬合非線性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
求高版本1stopt，解決一個(gè)問題，微分方程組參數(shù)擬合已經(jīng)有3人回復(fù)
用1stopt或者origin非線性擬合磁導(dǎo)率已經(jīng)有5人回復(fù)
有1stOpt4.0版本的朋友幫忙擬合一下復(fù)數(shù)吧感激不盡吶已經(jīng)有4人回復(fù)
1stopt或matlab如何求解以下的非線性方程并擬合出相應(yīng)曲線？已經(jīng)有13人回復(fù)
用origin或1stopt擬合一條隱函數(shù)方程曲線已經(jīng)有4人回復(fù)

1樓 2014-10-05 19:55:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

先把公式用簡單的字母簡化下，別用化學(xué)符號，更容易理解些。

贊一下

回復(fù)此樓

2樓2014-10-05 20:16:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

whqs8426212

銅蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 13.2
散金: 375
紅花: 7
帖子: 782
在線: 685.7小時(shí)
蟲號: 2572504
注冊: 2013-07-30
專業(yè): 數(shù)論

引用回帖:

2樓: Originally posted by dingd at 2014-10-05 20:16:08
先把公式用簡單的字母簡化下，別用化學(xué)符號，更容易理解些。

****************實(shí)驗(yàn)數(shù)據(jù)*************
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時(shí)，C(1)，C(2)的數(shù)值
***********函數(shù)*************************

fuction：
dC(1)=-k(1)*t;
dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1）

贊一下

回復(fù)此樓

3樓2014-10-06 00:51:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

第二個(gè)公式：dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1）
完全可簡化等同于：dC(2)=-k(2)*C(2)*C(1）

否則會(huì)出現(xiàn)過擬合現(xiàn)象，參數(shù)組值不唯一。

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值數(shù)量上相差太多了，相對合理的結(jié)果參考如下；
均方差(RMSE):0.00797587120432916
殘差平方和(SSE):0.00063614521468047
相關(guān)系數(shù)(R): 0.960512267688022
相關(guān)系數(shù)之平方(R^2): 0.922583816379187
確定系數(shù)(DC): 0.776373194295518
F統(tǒng)計(jì)(F-Statistic): 6.71167688419456

參數(shù) 最佳估算
-------------------- -------------
k1 0.526046976304358
k2 243.193548940279

c119.jpg

c120.jpg

贊一下

回復(fù)此樓

4樓2014-10-06 20:36:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

whqs8426212

銅蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 13.2
散金: 375
紅花: 7
帖子: 782
在線: 685.7小時(shí)
蟲號: 2572504
注冊: 2013-07-30
專業(yè): 數(shù)論

引用回帖:

4樓: Originally posted by dingd at 2014-10-06 20:36:44
第二個(gè)公式：dC(2)=-(k(2)*k(4)+k(3))*C(2)*C(1）
完全可簡化等同于：dC(2)=-k(2)*C(2)*C(1）

否則會(huì)出現(xiàn)過擬合現(xiàn)象，參數(shù)組值不唯一。

1stOpt求解：

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

什么是過擬合現(xiàn)象，是把特別小的參數(shù)忽略嗎

[ 發(fā)自小木蟲客戶端 ]

贊一下

回復(fù)此樓

5樓2014-10-07 16:28:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

擬合如下方程：y=a+b*c*x
無法得到穩(wěn)定解，有無窮多組b和c值，原方程等同于：y=a+b*x

贊一下

回復(fù)此樓

6樓2014-10-07 20:24:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

whqs8426212

銅蟲 (正式寫手)

應(yīng)助: 5 (幼兒園)
金幣: 13.2
散金: 375
紅花: 7
帖子: 782
在線: 685.7小時(shí)
蟲號: 2572504
注冊: 2013-07-30
專業(yè): 數(shù)論

引用回帖:

6樓: Originally posted by dingd at 2014-10-07 20:24:37
擬合如下方程：y=a+b*c*x
無法得到穩(wěn)定解，有無窮多組b和c值，原方程等同于：y=a+b*x

我將Code放入版本的1stOPT，結(jié)果報(bào)錯(cuò)'數(shù)據(jù)長度不夠‘
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;

贊一下

回復(fù)此樓

7樓2014-10-08 12:36:29

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

【答案】應(yīng)助回帖

data后面少“;”號。還有如果用的是網(wǎng)上流行的破解1.5版的話，是不能進(jìn)行微分方程擬合計(jì)算的。

贊一下

回復(fù)此樓

8樓2014-10-08 12:59:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

水中的夢88

木蟲 (正式寫手)

應(yīng)助: 52 (初中生)
金幣: 2378.9
散金: 26
紅花: 5
帖子: 323
在線: 195.7小時(shí)
蟲號: 1308847
注冊: 2011-05-28
性別: MM
專業(yè): 結(jié)構(gòu)工程

引用回帖:

8樓: Originally posted by dingd at 2014-10-08 12:59:11
data后面少“;”號。還有如果用的是網(wǎng)上流行的破解1.5版的話，是不能進(jìn)行微分方程擬合計(jì)算的。

請問那個(gè)版本以上的可以做，謝謝！

贊一下

回復(fù)此樓

9樓2015-10-19 20:59:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dingd

鐵桿木蟲 (職業(yè)作家)

計(jì)算強(qiáng)帖: 4
應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號: 291104
注冊: 2006-10-28

引用回帖:

9樓: Originally posted by 水中的夢88 at 2015-10-19 20:59:37
請問那個(gè)版本以上的可以做，謝謝！...

5.0以上吧

回復(fù)此樓

10樓2015-10-19 21:15:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[基金申請] 山東省面上項(xiàng)目限額評審 +3	石瑞0426 2026-03-19	3/150	2026-03-21 18:59 by RRRKKK
[考研] 求調(diào)劑 +3	.m.. 2026-03-21	4/200	2026-03-21 16:25 by barlinike
[考研] 求調(diào)劑 +3	白QF 2026-03-21	3/150	2026-03-21 13:12 by zhukairuo
[考研] 085601調(diào)劑 358分 +3	zzzzggh 2026-03-20	4/200	2026-03-21 10:21 by luoyongfeng
[考研] 化學(xué)求調(diào)劑 +4	臨澤境llllll 2026-03-17	5/250	2026-03-21 02:23 by JourneyLucky
[考研] 265求調(diào)劑 +9	梁梁校校 2026-03-17	9/450	2026-03-21 02:17 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 材料專業(yè)求調(diào)劑 +6	hanamiko 2026-03-18	6/300	2026-03-21 00:24 by JourneyLucky
[考研] 294求調(diào)劑材料與化工專碩 +15	陌の森林 2026-03-18	15/750	2026-03-20 23:28 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿中南化學(xué)（0703）總分337求調(diào)劑 +8	niko- 2026-03-19	9/450	2026-03-20 21:57 by luoyongfeng
[考研] 295材料求調(diào)劑，一志愿武漢理工085601專碩 +5	Charlieyq 2026-03-19	5/250	2026-03-20 20:35 by JourneyLucky
[考研] 求調(diào)劑 +3	eation27 2026-03-20	3/150	2026-03-20 19:32 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[論文投稿] 申請回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu