| 查看: 1412 | 回復(fù): 11 | |||
[求助]
請(qǐng)高版本1stopt代跑 已有2人參與
|
|
這是我自己寫(xiě)的代碼,麻煩大佬代跑一下,如果寫(xiě)的有錯(cuò)還請(qǐng)指教!拜托 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 |
鐵桿木蟲(chóng) (職業(yè)作家)
|
就是一個(gè)式子,ca,cb,cc,cd,ce是變量,也就是數(shù)據(jù),ka,kb,kc,kd,ke,kf是參數(shù),需要求這幾個(gè)參數(shù) 發(fā)自小木蟲(chóng)Android客戶(hù)端 |
版主 (知名作家)
|
計(jì)算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 優(yōu)化算法: 通用全局優(yōu)化算法(UGO1) 目標(biāo)函數(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é)果輸出 ====== 目標(biāo) ce 計(jì)算 ce 0.588533232 0.585434820888929 0.56087973 0.563011233652204 0.531698367 0.538199947346761 0.514160096 0.508410001555901 ====== 計(jì)算結(jié)束 ====== |

版主 (知名作家)
|
你已經(jīng)將 r 定義為常數(shù)了,除非r也有一組數(shù)據(jù) 發(fā)自小木蟲(chóng)Android客戶(hù)端 |

送紅花一朵 |
您看這樣行嗎,是不是這個(gè)代碼就可以使我的目標(biāo)函數(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 我再問(wèn)您幾個(gè)小問(wèn)題,怎么選取我的目標(biāo)函數(shù)是哪個(gè)呢?這兩種不同目標(biāo)函數(shù)的選取對(duì)我最后得到的參數(shù)有沒(méi)有影響呢? 拜托了,搞懂這些就去花錢(qián)讓人代跑,最后一步啦!。。 |
版主 (知名作家)
|
這樣是可以計(jì)算,但是每一行數(shù)據(jù)都會(huì)計(jì)算出一個(gè)r值出來(lái)。。。 ****************** 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)計(jì)(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. 實(shí)測(cè)值 r 計(jì)算值 r 1 3.2707 3.2692790256648 2 3.2707 3.27221542633868 3 3.2707 3.27274334006356 4 3.2707 3.26855834103416 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 280分 +15 | yyqqhh 2026-03-03 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 一志愿蘇大材料工程專(zhuān)碩293求調(diào)劑 +5 | 鋼鐵大炮 2026-03-04 | 5/250 |
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿江南大學(xué),食品科學(xué)與工程,總分,320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
[考研] 290求調(diào)劑 +6 | Bananaiy 2026-03-04 | 6/300 |
|
|
[考研] 334求調(diào)劑 +8 | Trying] 2026-03-06 | 8/400 |
|
|
[考研] 085701環(huán)境工程專(zhuān)業(yè),初試305,均過(guò)國(guó)家A區(qū)線 +7 | 卡卡來(lái)了@ 2026-03-07 | 8/400 |
|
|
[考研] 材料工程專(zhuān)碩調(diào)劑自薦信,初試323 +4 | 李白26 2026-03-07 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
[考研] 求調(diào)劑,不管什么專(zhuān)業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
[考研] 一志愿南京大學(xué)資源與環(huán)境專(zhuān)業(yè)283分求調(diào)劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|
[考研] 一志愿中石油,化學(xué)工程,301求調(diào)劑 +4 | 豚饅枳枳 2026-03-04 | 4/200 |
|
|
[考研] 273求調(diào)劑 +5 | 星星111222 2026-03-02 | 7/350 |
|
|
[考研] 考研282分求調(diào)劑,接受跨專(zhuān)業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
[考研] 一志愿武理085601專(zhuān)碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
[考研] 334求調(diào)劑 +3 | Trying] 2026-03-03 | 4/200 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 286 +6 | ksncj 2026-03-04 | 6/300 |
|
|
[考研] 293求調(diào)劑 +4 | 是樂(lè)渝哇 2026-03-03 | 4/200 |
|