瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

[求助] 請高版本1stopt代跑已有2人參與

這是我自己寫的代碼，麻煩大佬代跑一下，如果寫的有錯還請指教！拜托
Parameters KA,KB,KC,KD,KE,KF；
Variable CA,CB,CC,CD,CE,CF;
Constant r=3.2707;
Function r=KF(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE）^2;
Data；
3.914545455 2.146060606 4.249745353 2.737478991 0.588533232
3.828484848 2.232121212 4.308152457 2.706725388 0.56087973
3.767272727 2.293333333 4.340183215 2.703875994 0.531698367
3.711515152 2.349090909 4.37840252 2.683194959 0.514160096

回復此樓

» 猜你喜歡

重慶交通大學光子學微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息已經(jīng)有2人回復
一志愿鄭大材料學碩298分，求調(diào)劑已經(jīng)有6人回復
物理學I論文潤色/翻譯怎么收費? 已經(jīng)有263人回復
尋合作：應(yīng)力腐蝕多尺度模擬已經(jīng)有5人回復
考研交流已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有0人回復
幫我的英語口語老師找學生已經(jīng)有2人回復

1樓 2019-04-07 11:09:18

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

瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

不好意思著急寫錯了，沒有cf
Variable CA,CB,CC,CD,CE；

贊一下

回復此樓

2樓2019-04-07 11:10:13

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

瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

所有參數(shù)都大于0，再補充一下

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

贊一下

回復此樓

3樓2019-04-07 15:12:47

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

dingd

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

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

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

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

哪個是因變量？CE?

[ 發(fā)自手機版 http://m.gaoyang168.com/3g ]

贊一下

回復此樓

4樓2019-04-07 17:49:31

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

瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

就是一個式子，ca，cb，cc，cd，ce是變量，也就是數(shù)據(jù)，ka，kb，kc，kd，ke，kf是參數(shù)，需要求這幾個參數(shù)

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

贊一下

回復此樓

5樓2019-04-07 18:10:54

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

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

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

計算結(jié)束原因: 達到收斂判斷標準
優(yōu)化算法: 通用全局優(yōu)化算法(UGO1)
目標函數(shù)值(最小): 8.94776269027934E-5
均方差(RMSE): 0.00472962980459415
殘差平方和(SSR): 8.94775923540212E-5
相關(guān)系數(shù)(R): 0.986304513331934
確定系數(shù)(DC): 0.972186335183512

參數(shù) 最佳估算
---------- -------------
ka 11.8201249107019
kb 23.5726553856025
kc 0.843649223024818
kd 0.229359471088254
ke 6.45725464510755
kf 15.6674603816671

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

目標 ce 計算 ce
0.588533232 0.585434820888929
0.56087973 0.563011233652204
0.531698367 0.538199947346761
0.514160096 0.508410001555901

====== 計算結(jié)束 ======

贊一下

回復此樓

數(shù)值計算

6樓2019-04-07 18:39:55

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

瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

引用回帖:

6樓: Originally posted by 獨孤神宇 at 2019-04-07 18:39:55
計算結(jié)束原因: 達到收斂判斷標準
優(yōu)化算法: 通用全局優(yōu)化算法(UGO1)
目標函數(shù)值(最小): 8.94776269027934E-5
均方差(RMSE): 0.00472962980459415
殘差平方和(SSR): 8.94775923540212E-5
相關(guān)系數(shù)(R): 0.986304 ...

宇哥，我這個代碼怎么寫才能使目標與r對比呢？

贊一下

回復此樓

7樓2019-04-07 19:24:13

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

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

引用回帖:

7樓: Originally posted by 瑞aaa at 2019-04-07 19:24:13
宇哥，我這個代碼怎么寫才能使目標與r對比呢？...

你已經(jīng)將 r 定義為常數(shù)了，除非r也有一組數(shù)據(jù)

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

贊一下

回復此樓

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

瑞aaa

數(shù)值計算

8樓2019-04-07 20:07:44

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

瑞aaa

新蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 704.5
紅花: 3
沙發(fā): 1
帖子: 466
在線: 74小時
蟲號: 6018792
注冊: 2017-03-16
性別: GG
專業(yè): 化學工程及工業(yè)化學

送紅花一朵

引用回帖:

8樓: Originally posted by 獨孤神宇 at 2019-04-07 20:07:44
你已經(jīng)將 r 定義為常數(shù)了，除非r也有一組數(shù)據(jù)
...

您看這樣行嗎，是不是這個代碼就可以使我的目標函數(shù)變成r了嗎？
Parameters KA[0,],KB[0,],KC[0,],KD[0,],KE[0,],KF[0,]；
Variable CA,CB,CC,CD,CE;
Function r=KF(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE）^2;
Data；
3.2707 3.914545455       2.146060606       4.249745353       2.737478991       0.588533232
3.2707 3.828484848       2.232121212       4.308152457       2.706725388       0.56087973
3.2707 3.767272727       2.293333333       4.340183215       2.703875994       0.531698367
3.2707 3.711515152       2.349090909       4.37840252       2.683194959       0.514160096
我再問您幾個小問題，怎么選取我的目標函數(shù)是哪個呢？這兩種不同目標函數(shù)的選取對我最后得到的參數(shù)有沒有影響呢？
拜托了，搞懂這些就去花錢讓人代跑，最后一步啦�。。�！

贊一下

回復此樓

9樓2019-04-07 20:53:13

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

獨孤神宇

版主 (知名作家)

應(yīng)助: 490 (碩士)
貴賓: 0.008
金幣: 31016.3
散金: 802
紅花: 122
沙發(fā): 1
帖子: 5600
在線: 856.5小時
蟲號: 3522474
注冊: 2014-11-06
性別: GG
專業(yè): 機械動力學
管轄: 計算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
瑞aaa: 金幣+10, ★★★★★最佳答案 2019-04-13 07:07:08

引用回帖:

9樓: Originally posted by 瑞aaa at 2019-04-07 20:53:13
您看這樣行嗎，是不是這個代碼就可以使我的目標函數(shù)變成r了嗎？
Parameters KA,KB,KC,KD,KE,KF；
Variable CA,CB,CC,CD,CE;
Function r=KF(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE）^ ...

這樣是可以計算，但是每一行數(shù)據(jù)都會計算出一個r值出來。。。

******************

Parameters KA,KB,KC,KD,KE,KF;
ParameterDomain = [0, ];
Variable CA,CB,CC,CD,CE,r;
Function r=KF*(KA*KB*CA*CB-KC*KD*CC*CD/1.34)/(1+KA*CA+KB*CB+KC*CC+KD*CD+KE*CE)^2;
Data；
//CA,CB,CC,CD,CE,r;
3.914545455 2.146060606 4.249745353 2.737478991 0.588533232 3.2707
3.828484848 2.232121212 4.308152457 2.706725388 0.56087973 3.2707
3.767272727 2.293333333 4.340183215 2.703875994 0.531698367 3.2707
3.711515152 2.349090909 4.37840252 2.683194959 0.514160096 3.2707

*************************************************

均方差(RMSE): 0.00180815007608335
殘差平方和(SSR): 1.30776267905609E-5
相關(guān)系數(shù)(R): 0
相關(guān)系數(shù)之平方(R^2): 0
修正R平方(Adj. R^2): 2
確定系數(shù)(DC): 0
卡方系數(shù)(Chi-Square): 1.99922523760395E-6
F統(tǒng)計(F-Statistic): -0.4

參數(shù) 最佳估算
---------- -------------
ka 7.26870923240464
kb 76.8033683805592
kc 1.81885884331089E-8
kd 34.1212240739148
ke 80.0470405155196
kf 78.1355455967981

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

No. 實測值 r 計算值 r
1 3.2707 3.2692790256648
2 3.2707 3.27221542633868
3 3.2707 3.27274334006356
4 3.2707 3.26855834103416

贊一下

回復此樓

數(shù)值計算

10樓2019-04-07 21:33:19

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

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

[考研] 一志愿東北大學080500，英一數(shù)二343分求調(diào)劑 +5	Seele_v5 2026-03-04	5/250	2026-03-05 13:26 by Iveryant
[考研] 材料085600 303求調(diào)劑 +7	1bygone 2026-03-04	7/350	2026-03-05 09:39 by oxidpl
[考研] 0703 東華大學理學化學均過 a 線，總分 281 求調(diào)劑 +4	香香 Lu 2026-03-03	7/350	2026-03-05 09:38 by oxidpl
[考研] 320材料一志愿西工大專碩英二數(shù)二有兩年光伏工作經(jīng)驗 +5	錘某人 2026-03-04	8/400	2026-03-05 09:24 by oxidpl
[考研] 復試調(diào)劑 +3	呼呼？~+123456 2026-03-05	5/250	2026-03-05 08:19 by ms629
[考研] 環(huán)境調(diào)劑 +7	柒槿levana 2026-03-01	7/350	2026-03-04 22:00 by gwjxiaolang
[考研] 化學工程求調(diào)劑 +10	化工人999 2026-03-04	10/500	2026-03-04 21:21 by zhukairuo
[考研] 一志愿314求調(diào)劑 +7	202111120625 2026-03-03	7/350	2026-03-04 15:56 by zhukairuo
[考研] 281求調(diào)劑 +3	我是小小蔥蔥 2026-03-03	5/250	2026-03-04 14:23 by kakakapanpan
[考研] 一志愿西交大材料學碩 346 求調(diào)劑 +3	zju51 2026-03-04	3/150	2026-03-04 13:01 by zhukairuo
[考研] 085601 材料305分求助 +4	泡泡郵件 2026-03-03	6/300	2026-03-04 13:00 by 無際的草原
[考研] 復試調(diào)劑 +4	杰尼龜aaa 2026-03-03	4/200	2026-03-03 23:06 by zhukairuo
[考研] 298求調(diào)劑 +3	人間唯你是清歡 2026-03-03	4/200	2026-03-03 18:09 by ms629
[考研] 材料工程求調(diào)劑 +3	1431251 2026-03-03	3/150	2026-03-03 11:58 by EBSD
[考研] 0854復試調(diào)劑 276 +5	wmm9 2026-03-01	7/350	2026-03-03 02:49 by xiadaiyang
[考研] 285求調(diào)劑 +9	滿頭大汗的學生 2026-02-28	9/450	2026-03-02 20:29 by hypershenger
[考研] 275求調(diào)劑 +7	明遠求學 2026-03-01	7/350	2026-03-02 19:22 by zhukairuo
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學校有行學上就 +6	禮堂丁真258 2026-02-28	9/450	2026-03-02 12:04 by 52hz~~
[基金申請] 此成果不能導入原因：元數(shù)據(jù)必填信息不完整，可進行補充。 +4	Kittylucky 2026-03-02	5/250	2026-03-02 11:07 by jurkat.1640
[論文投稿] Optics letters投稿被拒求助 30+3	luckyry 2026-02-26	4/200	2026-03-01 09:06 by babero