| 3 | 1/1 | 返回列表 |
| 查看: 640 | 回復(fù): 2 | ||
xxz903金蟲 (小有名氣)
|
[求助]
1stOpt求助!無法計(jì)算~說是fatal:line too long(more than 1023 characters)
|
|
求高手解答啊,代碼如下: Title "單電解質(zhì)參數(shù)回歸"; Parameters r1,q1,r2,q2,t12,t21; Variable m,y1,y2; Constant z=10; conststr x1=m/(m+h); conststr x2=1-x1; conststr o1=x1*r1/(x1*r1+x2*r2); conststr o2=x2*r2/(x1*r1+x2*r2); conststr b1=x1*q1/(x1*q1+x2*q2); conststr b2=x2*q2/(x1*q1+x2*q2); conststr l1=0.5*z*(r1-q1)-(r1-1); conststr l2=0.5*z*(r2-q2)-(r2-1); SharedModel; Function y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12)); Function y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21)); data; 0 0 0 0.1 -0.673344553 -0.003214304 0.2 -0.783071888 -0.005447496 0.3 -0.841647189 -0.008713882 0.4 -0.869884359 -0.010990339 0.5 -0.884307686 -0.014309599 0.6 -0.889162064 -0.018687569 0.7 -0.889162064 -0.02207116 0.8 -0.88673193 -0.025484701 0.9 -0.879476759 -0.029975034 1 -0.869884359 -0.038731848 1.2 -0.850971266 -0.042636996 1.4 -0.830113036 -0.051995295 1.6 -0.807436327 -0.061533872 1.8 -0.783071888 -0.071257081 2 -0.759286983 -0.081169452 2.5 -0.701179352 -0.107972916 3 -0.650087691 -0.136153773 3.5 -0.599656837 -0.165810242 4 -0.553385238 -0.197052419 4.5 -0.512493681 -0.227286842 5 -0.47320876 -0.260553522 |
鐵桿木蟲 (職業(yè)作家)
金蟲 (小有名氣)
|
謝謝了!上一個(gè)的值有點(diǎn)問題 不應(yīng)該是負(fù)的,我修改了一個(gè)目標(biāo)函數(shù) 能再幫我擬合一個(gè)嗎? Title "單電解質(zhì)參數(shù)回歸"; Parameters r1,q1,r2,q2,t12,t21; Variable m,y1,y2; Constant z=10; conststr x1=m/(m+h); conststr x2=1-x1; conststr o1=x1*r1/(x1*r1+x2*r2); conststr o2=x2*r2/(x1*r1+x2*r2); conststr b1=x1*q1/(x1*q1+x2*q2); conststr b2=x2*q2/(x1*q1+x2*q2); conststr l1=0.5*z*(r1-q1)-(r1-1); conststr l2=0.5*z*(r2-q2)-(r2-1); SharedModel; Function y1=ln(o1/x1)+0.5*z*q1*ln(o1/b1)+o2*(l1-r1*l2/r2)-q1*ln(b1+b2*t21)+b2*q1*(t21/(b1+b2*t21)-t12/(b2+b1*t12))-ln(r1/r2)-05*z*q1*ln(q1*r2/(q2*r1))-(l1-r1*l2/r2)+q1*ln(t21)-q1*(1-t12); Function y2=ln(o2/x2)+0.5*z*q2*ln(o2/b2)+o1*(l2-r2*l1/r1)-q2*ln(b2+b1*t12)+b1*q2*(t12/(b2+b1*t12)-t21/(b1+b2*t21))-ln(r2/r1)-05*z*q2*ln(q2*r1/(q1*r2))-(l2-r2*l1/r1)+q2*ln(t12)-q2*(1-t21); data; 0 0 0 0.1 -0.673344553 -0.003214304 0.2 -0.783071888 -0.005447496 0.3 -0.841647189 -0.008713882 0.4 -0.869884359 -0.010990339 0.5 -0.884307686 -0.014309599 0.6 -0.889162064 -0.018687569 0.7 -0.889162064 -0.02207116 0.8 -0.88673193 -0.025484701 0.9 -0.879476759 -0.029975034 1 -0.869884359 -0.038731848 1.2 -0.850971266 -0.042636996 1.4 -0.830113036 -0.051995295 1.6 -0.807436327 -0.061533872 1.8 -0.783071888 -0.071257081 2 -0.759286983 -0.081169452 2.5 -0.701179352 -0.107972916 3 -0.650087691 -0.136153773 3.5 -0.599656837 -0.165810242 4 -0.553385238 -0.197052419 4.5 -0.512493681 -0.227286842 5 -0.47320876 -0.260553522 |
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|