版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

[求助] 求助大神幫忙用1stopt擬合下模型參數(shù) 已有1人參與

在下小碩
想請(qǐng)大神幫忙用1stopt跑下模型參數(shù)
數(shù)據(jù)
X=3091  3281  3508  3904  4108
Y=2874  3054  3271  3634  3832
Z=0.336  0.360  0.382  0.385  0.397
模型
dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*x*y
dy/dt=b1*x+b2*y-b3*z-b4*y*z
dz/dt=-c1*x+c2*y+c3*z-c4*z*z
參數(shù)
a1 a2 a3 a4 a5
b1 b2 b3 b4
c1 c2 c3 c4
如果可能的化，希望誤差可以<10^(-5)
感謝大神的幫助！

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有228人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有53人回復(fù)

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

1stopt數(shù)據(jù)擬合求助已經(jīng)有8人回復(fù)
1stopt軟件擬合非線(xiàn)性方程，計(jì)算結(jié)果全為0，這是為什么已經(jīng)有3人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有10人回復(fù)
求大神跑下1stOpt，一個(gè)sharemodel，謝謝！已經(jīng)有7人回復(fù)
呼叫版主，在線(xiàn)緊急求助，關(guān)于matlab中微分方程組參數(shù)擬合得問(wèn)題！已經(jīng)有12人回復(fù)
跪求高手幫忙看看這組數(shù)據(jù)如何進(jìn)行擬合已經(jīng)有7人回復(fù)
求大神幫你擬合一組數(shù)據(jù)~萬(wàn)分感激。。。已經(jīng)有6人回復(fù)
求大神指導(dǎo)1stopt 非線(xiàn)性最優(yōu)化問(wèn)題已經(jīng)有3人回復(fù)
求助--用1stopt擬合藍(lán)寶石襯底的GaN薄膜紅外反射光譜已經(jīng)有6人回復(fù)
求助1stopt擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有4人回復(fù)
求大神幫忙擬合一個(gè)非線(xiàn)性方程，求出模型參數(shù) 已經(jīng)有15人回復(fù)
拜求！matlab擬合動(dòng)力學(xué)參數(shù) 已經(jīng)有14人回復(fù)
1stopt擬合問(wèn)題求助已經(jīng)有7人回復(fù)
求大神幫忙用1stOpt 擬合復(fù)數(shù)數(shù)據(jù) 已經(jīng)有13人回復(fù)
1stopt中使用sharedmodel擬合曲線(xiàn)出現(xiàn)問(wèn)題，求指示��！已經(jīng)有6人回復(fù)
用1stopt擬合非線(xiàn)性方程結(jié)果與其他軟件擬合結(jié)果差異大已經(jīng)有5人回復(fù)
大家在使用1stOpt進(jìn)行數(shù)據(jù)擬合的時(shí)候有沒(méi)有遇到過(guò)這種情況？已經(jīng)有8人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】多變量非線(xiàn)性擬合【已完結(jié)】已經(jīng)有8人回復(fù)
1stopt擬合編程已經(jīng)有1人回復(fù)

1樓 2014-08-05 07:11:54

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

dingd

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

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

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
ziliaohan: 金幣+5, ★★★★★最佳答案 2014-08-05 16:58:04

給你一組參考：

均方差(RMSE):0.898138086625654
殘差平方和(SSE):9.67982427177108
相關(guān)系數(shù)(R): 0.936001558277226
相關(guān)系數(shù)之平方(R^2): 0.876098917097396
確定系數(shù)(DC): -32.6666746162752
F統(tǒng)計(jì)(F-Statistic): -78192.6674901789

參數(shù) 最佳估算
-------------------- -------------
a1 0.41660768767696
a2 -0.596899303711546
a3 -1672.37950977148
a4 0.057929293839314
a5 -2.54346768480167E-6
b1 0.391376424549718
b2 -0.295297695034362
b3 1025.02219900667
b4 0.0520438936287644
c1 -0.0915214045518645
c2 -0.098248199763037
c3 -2.47400188768122
c4 -10.6412526763723

數(shù)據(jù)應(yīng)該有些問(wèn)題：x、y與z大小存在數(shù)量級(jí)的差別，會(huì)嚴(yán)重影響數(shù)量級(jí)�。ū热鐉）的擬合精度。

贊一下(1人)

回復(fù)此樓

2樓2014-08-05 09:46:34

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

引用回帖:

2樓: Originally posted by dingd at 2014-08-05 09:46:34
給你一組參考：

均方差(RMSE):0.898138086625654
殘差平方和(SSE):9.67982427177108
相關(guān)系數(shù)(R): 0.936001558277226
相關(guān)系數(shù)之平方(R^2): 0.876098917097396
確定系數(shù)(DC): -32.6666746162752
F統(tǒng)計(jì)(F-St ...

那我把數(shù)據(jù)歸一化處理下，應(yīng)該會(huì)好一些

贊一下

回復(fù)此樓

3樓2014-08-05 16:41:59

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

內(nèi)容已刪除

回復(fù)此樓

4樓2014-08-05 16:55:50

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

dingd

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

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

x,y,z的初值都為0的話(huà)，后面的計(jì)算值都全為0；數(shù)據(jù)歸一的有問(wèn)題。

贊一下(2人)

回復(fù)此樓

5樓2014-08-05 20:06:40

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

引用回帖:

5樓: Originally posted by dingd at 2014-08-05 20:06:40
x,y,z的初值都為0的話(huà)，后面的計(jì)算值都全為0；數(shù)據(jù)歸一的有問(wèn)題。

感謝大神，耐心的幫助
我重新歸一處理下數(shù)據(jù)
還請(qǐng)大神再跑一次
X = 0.1728 0.1834 0.1961 0.2182 0.2296
Y = 0.1725 0.1833 0.1963 0.2181 0.2299
Z = 0.1806 0.1935 0.2054 0.2070 0.2134

贊一下

回復(fù)此樓

6樓2014-08-06 11:19:13

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

dingd

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

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

計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:30:03:00
均方差(RMSE):0.000238749203430912
殘差平方和(SSE):6.84014185666743E-7
相關(guān)系數(shù)(R): 0.999852515745996
相關(guān)系數(shù)之平方(R^2): 0.999705053243598
確定系數(shù)(DC): 0.999659468900762
F統(tǒng)計(jì)(F-Statistic): -6383.89529824175

參數(shù) 最佳估算
-------------------- -------------
a1 29.5516082698601
a2 -37.0355063680216
a3 -4.33534796194097
a4 3.37837345497585
a5 -18.2156770829613
b1 -31.1552076553908
b2 32.0949018446089
b3 -0.525118791751513
b4 6.63207723836531
c1 -23.5149166671026
c2 -23.5224081410208
c3 -0.976591458659144
c4 -4.870173216137

====== 結(jié)果輸出 ======

文件: 數(shù)據(jù)文件-1
No t 目標(biāo) x 計(jì)算 x 目標(biāo) y 計(jì)算 y 目標(biāo) z 計(jì)算 z
1 1 0.1834 0.183057385327924 0.1833 0.183837598389866 0.1935 0.193615350874645
2 2 0.1961 0.196269323838988 0.1963 0.195960460145724 0.2054 0.205140005956856
3 3 0.2182 0.218285645083849 0.2181 0.218090044028013 0.207 0.207201021981877
4 4 0.2296 0.229627009142356 0.2299 0.229855468869106 0.2134 0.21335309384729

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

ziliaohan

7樓2014-08-06 16:49:56

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

送紅花一朵

內(nèi)容已刪除

回復(fù)此樓

8樓2014-08-07 07:18:33

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

剛才把結(jié)果給老板看了下
比較滿(mǎn)意,和預(yù)想的有些出入
建議把模型中的最后一項(xiàng)修改成a5*y*y，會(huì)更好
即，模型為
dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*y*y
dy/dt=b1*x+b2*y-b3*z-b4*y*z
dz/dt=-c1*x+c2*y+c3*z-c4*z*z
參數(shù)還是
X = 0.1728  0.1834  0.1961  0.2182  0.2296
Y = 0.1725  0.1833  0.1963  0.2181  0.2299
Z = 0.1806  0.1935  0.2054  0.2070  0.2134
能不能還請(qǐng)跑一下
要求有點(diǎn)過(guò)分
還請(qǐng)大神原諒
先謝了

贊一下

回復(fù)此樓

9樓2014-08-07 08:24:59

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

ziliaohan

銅蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 47
散金: 10
帖子: 36
在線(xiàn): 33.6小時(shí)
蟲(chóng)號(hào): 2231633
注冊(cè): 2013-01-09
專(zhuān)業(yè): 偏微分方程

引用回帖:

7樓: Originally posted by dingd at 2014-08-06 16:49:56
計(jì)算用時(shí)(時(shí):分:秒:微秒): 00:30:03:00
均方差(RMSE):0.000238749203430912
殘差平方和(SSE):6.84014185666743E-7
相關(guān)系數(shù)(R): 0.999852515745996
相關(guān)系數(shù)之平方(R^2): 0.999705053243598
確定系數(shù)(DC): 0.9 ...

還想請(qǐng)大神再幫忙跑一下
主要是模型中的a5x*y項(xiàng)需要修改成a5y*y，
這樣模型才符合理論推導(dǎo)
原先那個(gè)項(xiàng)推導(dǎo)錯(cuò)了
老板要求重來(lái)一次
dx/dt=-a1*x-a2*y+a3*z+a4*x*z+a5*y*y
dy/dt=b1*x+b2*y-b3*z-b4*y*z
dz/dt=-c1*x+c2*y+c3*z-c4*z*z
參數(shù)還是
X = 0.1728　.1834　.1961　.2182　.2296
Y = 0.1725　.1833　.1963　.2181　.2299
Z = 0.1806　.1935　.2054　.2070　.2134
感謝大神

贊一下

回復(fù)此樓

10樓2014-08-08 07:01:33

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 ziliaohan 的主題更新

返回列表

普通表情龍兔虎貓

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

[考研] 321求調(diào)劑 +3	何潤(rùn)采123 2026-03-18	3/150	2026-03-18 21:27 by li123456789.
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 0703化學(xué)調(diào)劑 +7	妮妮ninicgb 2026-03-15	11/550	2026-03-18 19:48 by macy2011
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶(hù) 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專(zhuān)家開(kāi)始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黃鳥(niǎo)于飛Chao
[考研] 308求調(diào)劑 +4	是Lupa啊 2026-03-16	4/200	2026-03-17 17:12 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程專(zhuān)碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6	iwuli 2026-03-12	7/350	2026-03-16 21:18 by sculhf
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 0856專(zhuān)碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 中科大材料專(zhuān)碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 308求調(diào)劑 +3	是Lupa啊 2026-03-12	3/150	2026-03-13 14:30 by 求調(diào)劑zz

<rt id="wdwop"></rt>