瑞aaa

新蟲 (正式寫手)

應助: 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

回復此樓

» 猜你喜歡

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

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

瑞aaa

新蟲 (正式寫手)

應助: 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

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

瑞aaa

新蟲 (正式寫手)

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

所有參數都大于0，再補充一下

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

贊一下

回復此樓

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

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

dingd

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

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

【答案】應助回帖

感謝參與，應助指數 +1

哪個是因變量？CE?

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

贊一下

回復此樓

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

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

瑞aaa

新蟲 (正式寫手)

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

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

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

贊一下

回復此樓

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

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

獨孤神宇

版主 (知名作家)

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

【答案】應助回帖

感謝參與，應助指數 +1

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

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

====== 結果輸出 ======

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

====== 計算結束 ======

贊一下

回復此樓

數值計算

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

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

瑞aaa

新蟲 (正式寫手)

應助: 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
計算結束原因: 達到收斂判斷標準
優(yōu)化算法: 通用全局優(yōu)化算法(UGO1)
目標函數值(最小): 8.94776269027934E-5
均方差(RMSE): 0.00472962980459415
殘差平方和(SSR): 8.94775923540212E-5
相關系數(R): 0.986304 ...

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

贊一下

回復此樓

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

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

獨孤神宇

版主 (知名作家)

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

【答案】應助回帖

引用回帖:

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

你已經將 r 定義為常數了，除非r也有一組數據

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

贊一下

回復此樓

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

瑞aaa

數值計算

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

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

瑞aaa

新蟲 (正式寫手)

應助: 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
你已經將 r 定義為常數了，除非r也有一組數據
...

您看這樣行嗎，是不是這個代碼就可以使我的目標函數變成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
我再問您幾個小問題，怎么選取我的目標函數是哪個呢？這兩種不同目標函數的選取對我最后得到的參數有沒有影響呢？
拜托了，搞懂這些就去花錢讓人代跑，最后一步啦�。。。�

贊一下

回復此樓

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

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

獨孤神宇

版主 (知名作家)

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

【答案】應助回帖

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

引用回帖:

9樓: Originally posted by 瑞aaa at 2019-04-07 20:53:13
您看這樣行嗎，是不是這個代碼就可以使我的目標函數變成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）^ ...

這樣是可以計算，但是每一行數據都會計算出一個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
相關系數(R): 0
相關系數之平方(R^2): 0
修正R平方(Adj. R^2): 2
確定系數(DC): 0
卡方系數(Chi-Square): 1.99922523760395E-6
F統(tǒng)計(F-Statistic): -0.4

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

====== 結果輸出 ======

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

贊一下

回復此樓

數值計算

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

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

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

[考研] 材料調劑 +3	ounce. 2026-03-05	6/300	2026-03-06 00:21 by wutongshun
[考博] 2026年博士名額撿漏 +4	科研ya 2026-03-04	5/250	2026-03-05 23:49 by 愛學學啊
[考研] 考研282分求調劑，接受跨專業(yè) +4	劉淄博 2026-03-04	7/350	2026-03-05 22:56 by 劉淄博
[考研] 313求調劑 +3	水流年lc 2026-02-28	4/200	2026-03-05 14:58 by 水流年lc
[考研] 347求調劑 +6	啊歐歐歐 2026-03-03	8/400	2026-03-05 11:40 by 0202liuyan
[考研] 一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學材料，有六級，有工作經驗 +3	H14528 2026-03-04	3/150	2026-03-05 10:35 by ms629
[考研] 085600材料調劑總分330 +5	池池丶 2026-03-03	5/250	2026-03-05 09:41 by oxidpl
[考研] 化學 0703求調劑總分293 一志愿211 +9	土土小蟲 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 070300化學 280 一志愿太原理工求調劑 +8	拾玖壹 2026-03-04	8/400	2026-03-05 08:52 by anny19840123
[考研] 材料學碩080500復試調劑294 +3	四葉zjz 2026-03-04	3/150	2026-03-05 07:09 by kunm555
[考研] 一志愿西交大材料學碩 346 求調劑 +3	zju51 2026-03-04	3/150	2026-03-04 13:01 by zhukairuo
[考研] 主題 +3	realstar2006 2026-02-27	3/150	2026-03-03 21:39 by AC.
[考研] 298求調劑一志愿中海洋 +3	lour. 2026-03-03	3/150	2026-03-03 20:41 by gxg2025
[考研] 江蘇省農科院招調劑1名 +5	Qwertyuop 2026-03-01	5/250	2026-03-03 17:50 by 梅子鴨鴨
[考研] 一志愿中科大能動297求調劑，本科川大 +3	邵11 2026-03-03	3/150	2026-03-03 14:50 by 新能源達人
[考研] 284求調劑 +6	天下熯 2026-03-02	6/300	2026-03-03 13:46 by Iveryant
[考研] 調劑 +5	13853210211 2026-03-02	7/350	2026-03-03 11:42 by 13853210211
[考研] 化工專碩342，一志愿大連理工大學，求調劑 +6	kyf化工 2026-02-28	7/350	2026-03-02 10:56 by 無際的草原
[考研] 275求調劑 +3	L-xin? 2026-03-01	6/300	2026-03-02 10:22 by 熱情沙漠
[論文投稿] 求助coordination chemistry reviews 的寫作模板 10+3	ljplijiapeng 2026-02-27	4/200	2026-03-01 09:07 by babero