當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

明雨天地

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 65
紅花: 2
帖子: 46
在線: 26小時(shí)
蟲號: 2211405
注冊: 2012-12-28
性別: GG
專業(yè): 環(huán)境污染化學(xué)

[求助] 求大哥大姐幫跑一下1stopt~新手啊，版本低啊已有1人參與

自己的軟件版本太低，

盼好心人、高版本1stopt出現(xiàn)~~~
非常感謝，

俺也知道大家的時(shí)間很寶貴，所以盡量多給金幣了~~~
一共有兩批數(shù)據(jù)，都是同樣的擬合過程

Constant a1=0.0001,a2=0.001,a3=0.01,a4=0.1,a5=1,a6=10,a7=100,a8=1000;
Parameter g(1:8);
Variables y1,y2,x;
SharedModel;
Function y1=g1*a1^2*x^2/(1+a1^2*x^2)+g2*a2^2*x^2/(1+a2^2*x^2)+g3*a3^2*x^2/(1+a3^2*x^2)+g4*a4^2*x^2/(1+a4^2*x^2)+g5*a5^2*x^2/(1+a5^2*x^2)+g6*a6^2*x^2/(1+a6^2*x^2)+g7*a7^2*x^2/(1+a7^2*x^2)+g8*a8^2*x^2/(1+a8^2*x^2);
y2=g1*a1*x/(1+a1^2*x^2)+g2*a2*x/(1+a2^2*x^2)+g3*a3*x/(1+a3^2*x^2)+g4*a4*x/(1+a4^2*x^2)+g5*a5*x/(1+a5^2*x^2)+g6*a6*x/(1+a6^2*x^2)+g7*a7*x/(1+a7^2*x^2)+g8*a8*x/(1+a8^2*x^2);
Data;
第一組數(shù)據(jù)：
714.954 218.475 0.0159155
790.545 205.354 0.0252244
843.715 191.499 0.0399779
880.668 172.557 0.0633607
916.241 179.116 0.10042
956.146 175.205 0.159155
995.362 179.39 0.252244
1039.92 184.729 0.39978
1082.48 192.036 0.633604
1127.35 199.892 1.0042
1173.93 208.977 1.59156
1222.63 219.333 2.52245
1273.37 231.457 3.99778
1325.44 245.833 6.33604
1376.67 263.674 10.042
1421.09 284.422 15.9155
第二組數(shù)據(jù)：
41.6235 69.5311 0.0159155
42.1067 72.1866 0.0252244
42.9638 77.6567 0.0399779
45.8138 84.8211 0.0633607
49.2256 94.4187 0.10042
54.8413 107.288 0.159155
63.4547 123.353 0.252244
76.7593 142.65 0.39978
97.0701 162.804 0.633604
124.786 178.922 1.0042
155.685 187.642 1.59156
183.444 189.678 2.52245
205.171 189.985 3.99778
220.368 192.472 6.33604
227.868 199.387 10.042
229.082 207.257 15.9155

回復(fù)此樓

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

明雨天地

» 猜你喜歡

請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有32人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有265人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有78人回復(fù)
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復(fù)
中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)

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

有沒有好心人告訴個管理員賬號或者明天幫忙查詢一下啊已經(jīng)有15人回復(fù)
求大神跑下1stOpt，一個sharemodel，謝謝！已經(jīng)有7人回復(fù)
請問一下，關(guān)于微分方程收斂性的問題已經(jīng)有4人回復(fù)
求大神利用1stopt幫我跑下程序已經(jīng)有5人回復(fù)
1stopt擬合出參數(shù)為0，為什么？已經(jīng)有6人回復(fù)
求vasp5.2.8版本已經(jīng)有6人回復(fù)
求高版本1stopt跑一個程序已經(jīng)有7人回復(fù)
有沒有人會用1stopt或者matlab擬合這樣數(shù)據(jù)，急求啊已經(jīng)有12人回復(fù)
求1stopt的安裝包已經(jīng)有8人回復(fù)
一區(qū)英語46求各位大哥大姐幫忙已經(jīng)有6人回復(fù)
1stopt解方程組已經(jīng)有12人回復(fù)
用1stopt軟件做微分方程組的參數(shù)估計(jì) 已經(jīng)有13人回復(fù)
1stopt求最小值已經(jīng)有37人回復(fù)
誰能幫我用高版本的 1stopt擬合一下，需要涉及復(fù)數(shù) 已經(jīng)有7人回復(fù)
求高版本1stopt擬合，已經(jīng)有11人回復(fù)
求大哥大姐幫跑一下1stopt~新手啊，版本低啊已經(jīng)有9人回復(fù)

1樓 2014-08-28 22:01:40

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

明雨天地

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 65
紅花: 2
帖子: 46
在線: 26小時(shí)
蟲號: 2211405
注冊: 2012-12-28
性別: GG
專業(yè): 環(huán)境污染化學(xué)

引用回帖:

2樓: Originally posted by dingd at 2014-08-28 23:02:17
第一組：
均方差(RMSE):13.893467804642
殘差平方和(SSE):6176.91032443595
相關(guān)系數(shù)(R): 0.981690409944436
相關(guān)系數(shù)之平方(R^2): 0.963716060976875
確定系數(shù)(DC): 0.868312941968537
F統(tǒng)計(jì)(F-Statistic): ...

謝謝親及時(shí)準(zhǔn)確的回復(fù)，太高興了。還有我想請教一下，您用的算法是通用全局優(yōu)化法(UGO1)？因?yàn)閿M合出來的g值要>0，所以如果不滿足要求的話，可能Constant a的值需要調(diào)整，Constant a（1:8）在0.0001~1000范圍內(nèi)由小到大的順序可以隨便調(diào)整，麻煩諸位大蝦了

贊一下

回復(fù)此樓

3樓2014-08-29 10:30:03

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

查看全部 7 個回答

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
明雨天地: 金幣+15, ★★★★★最佳答案 2014-08-29 10:34:25

第一組：
均方差(RMSE):13.893467804642
殘差平方和(SSE):6176.91032443595
相關(guān)系數(shù)(R): 0.981690409944436
相關(guān)系數(shù)之平方(R^2): 0.963716060976875
確定系數(shù)(DC): 0.868312941968537
F統(tǒng)計(jì)(F-Statistic): 409.372045865983

參數(shù) 最佳估算
-------------------- -------------
g1 21898585.462287
g2 -2373710.33510328
g3 19593.7601667467
g4 306.259092939804
g5 250.101741740911
g6 187.751871369503
g7 383.077779809732
g8 436.806292860088

第二組：
均方差(RMSE):33.521726288192
殘差平方和(SSE):35958.5962668948
相關(guān)系數(shù)(R): 0.933337714249498
相關(guān)系數(shù)之平方(R^2): 0.871119288840478
確定系數(shù)(DC): 0.620325659761535
F統(tǒng)計(jì)(F-Statistic): 8.99815686416928

參數(shù) 最佳估算
-------------------- -------------
g1 90511441.7006428
g2 -9939268.36476249
g3 92069.0559634856
g4 -29.7624522133981
g5 208.227698210748
g6 10.030656082362
g7 180.739905721451
g8 -147.54790424532

贊一下

回復(fù)此樓

2樓2014-08-28 23:02:17

已閱回復(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

引用回帖:

3樓: Originally posted by 明雨天地 at 2014-08-29 10:30:03

謝謝親及時(shí)準(zhǔn)確的回復(fù)，太高興了。還有我想請教一下，您用的算法是通用全局優(yōu)化法(UGO1)？因?yàn)閿M合出來的g值要>0，所以如果不滿足要求的話，可能Constant a的值需要調(diào)整，Constant a（1:8）在0.0001~1000 ...

參數(shù)范圍有限制條件，你不在開始就說明，誰能知道？

贊一下

回復(fù)此樓

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

明雨天地

4樓2014-08-29 11:06:58

已閱回復(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ún)?shù)都大于0，將“Parameter g(1:8);”改為“Parameter g(1:8)=[0,];”就行了：
第一組：
均方差(RMSE):17.068661604803
殘差平方和(SSE):9322.85468733688
相關(guān)系數(shù)(R): 0.971809999263162
相關(guān)系數(shù)之平方(R^2): 0.944414674667868
確定系數(shù)(DC): 0.787179155123501
F統(tǒng)計(jì)(F-Statistic): 322.766185328756

參數(shù) 最佳估算
-------------------- -------------
g1 110512.408867717
g2 5.44139883889035E-10
g3 5.46272427927806E-13
g4 255.175994425224
g5 257.777798483887
g6 182.247746525854
g7 388.040093025756
g8 432.334590079376

第二組：
均方差(RMSE):44.6901637863839
殘差平方和(SSE):63910.7436561223
相關(guān)系數(shù)(R): 0.944164598601182
相關(guān)系數(shù)之平方(R^2): 0.891446789251732
確定系數(shù)(DC): 0.268820694782057
F統(tǒng)計(jì)(F-Statistic): 9.15928197962095

參數(shù) 最佳估算
-------------------- -------------
g1 152778.991526402
g2 3.7122240049481E-11
g3 5.42311048362141E-12
g4 42.2422780669557
g5 192.74240634343
g6 35.6287907076214
g7 11.9918051835019
g8 6.21823318030801E-14

贊一下

回復(fù)此樓

6樓2014-08-29 11:36:10

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

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

[考研] 材料調(diào)劑 +3	匹克i 2026-03-23	3/150	2026-03-23 23:18 by peike
[考研] 291求調(diào)劑 +8	hhhhxn.. 2026-03-23	8/400	2026-03-23 23:15 by peike
[考研] 0854 考研調(diào)劑招生了！AI 方向 +4	pk3725069 2026-03-19	16/800	2026-03-23 23:09 by 汪��？！
[考研] 一志愿211 初試270分求調(diào)劑 +4	谷雨上岸 2026-03-23	5/250	2026-03-23 21:18 by 不惑可樂
[考研] 336求調(diào)劑 +4	收到VS 2026-03-20	4/200	2026-03-23 19:02 by macy2011
[考研] 考研化學(xué)308分求調(diào)劑 +7	你好明天你好 2026-03-23	8/400	2026-03-23 18:39 by macy2011
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +7	生物工程調(diào)劑 2026-03-17	12/600	2026-03-23 18:18 by YMU施老師
[考研] 352求調(diào)劑 +3	大米飯！ 2026-03-22	3/150	2026-03-22 23:28 by king123！
[考研] 289求調(diào)劑 +7	懷瑾握瑜l 2026-03-20	7/350	2026-03-22 15:57 by ColorlessPI
[考研] 299求調(diào)劑 +5	shxchem 2026-03-20	7/350	2026-03-21 17:09 by ColorlessPI
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 330求調(diào)劑 +4	小材化本科 2026-03-18	4/200	2026-03-20 23:13 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 290求調(diào)劑 +7	^O^乜 2026-03-19	7/350	2026-03-20 21:43 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 0856調(diào)劑，是學(xué)校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +5	ioodiiij 2026-03-17	5/250	2026-03-19 18:22 by zcl123
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴