| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1212 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者唯有千尋將贈(zèng)送您 100 個(gè)金幣 | ||
| 當(dāng)前只顯示滿(mǎn)足指定條件的回帖,點(diǎn)擊這里查看本話(huà)題的所有回帖 | ||
唯有千尋新蟲(chóng) (小有名氣)
|
[求助]
1stOpt多參數(shù)約束擬合代碼錯(cuò)誤 已有1人參與
|
|
求助1stopt多參數(shù)非線(xiàn)性代碼錯(cuò)誤,去掉保證擴(kuò)散系數(shù)大于0的兩行代碼可以運(yùn)行,得出的解和文獻(xiàn)結(jié)果相差甚遠(yuǎn),加入兩行約束條件程序報(bào)錯(cuò),請(qǐng)大神指點(diǎn)如何加上約束條件能正常運(yùn)行 代碼如下:Title "Nonlinear Fit for Defect diffusion coefficient"; Constant F=96485.334,R=8.3145,b3=0,d3=0;//法拉第常數(shù)F,R氣體常數(shù),文獻(xiàn)中的兩個(gè)參數(shù)設(shè)置為0; VarConstant T=[773.15,973.15,1073.15]; Parameters a1[0,1E-4],b1,c1,d1,a2[0,1E-4],b2,c2,d2,a3[0,1E-4],c3;//a1,a2,a3為參數(shù)預(yù)指數(shù)因子; Variables y,x1,x2,x3,x4;//y為離子導(dǎo)電率,x1,x2,x3為各離子濃度 Function y=(F^2)*x1*(a1*(1+b1*x4))*exp(-(c1+d1*x4)/(R*T))/(R*T)+(F^2)*x2*(a2*(1+b2*x4))*exp(-(c2+d2*x4)/(R*T))/(R*T)+(4*F^2)*x3*(a3*(1+b3*x4))*exp(-(c3+d3*x4)/(R*T))/(R*T); a1*(1+b1*x4)>0;//保證擴(kuò)散系數(shù)大于0; a2*(1+b2*x4)>0;//保證擴(kuò)散系數(shù)大于0; RowData; 119685.18,124330.77,127828.002,130751.07,132786.77; 4403.185,4851.994,5364.097,5588.501,5715.09; 214.3,185.53,151.007,116.483,87.71; 789.697,501.999,312.118,185.531,151.007; 0.0274,0.0174,0.0108,0.0064,0.0052; RowData; 50940.98,64564.546,69523.315,80276.015,84869.4,94160.57,97083.63,103764.92,106061.615,106374.8; 983.51,1362.565,1621.014,2160.885,2448.05,2993.664,3246.37,3694.35,3883.88,4205.5; 2034.532,1460.2,1236.21,914.587,822.694,627.422,564.246,472.353,437.893,409.176; 3441.642,2930.49,2706.499,2321.697,2068.992,1712.91,1586.554,1299.39,1173.04,1012.223; 0.119,0.1015,0.09375,0.0804,0.0717,0.05934,0.055,0.045,0.0406,0.0351; RowData; 25155.38,37004.228,41649.812,53864.044,58509.628,70097.488,74429.887,83460.067,88105.65; 593.555,791.372,940.732,1350.516,1509.876,2141.63,2397.74,2966.884,3251.455; 4071.022,3092.095,2682.312,1919.659,1669.236,1225.304,1065.944,815.521,752.915; 4674.315,4071.022,3757.993,3188.85,2932.735,2460.346,2204.231,1828.596,1669.236; 0.162,0.141,0.13,0.1105,0.1016,0.085,0.0764,0.063,0.0578; 程序報(bào)錯(cuò): Incompatible types: 'Double' and 'Boolean' Incompatible types: 'Double' and 'Boolean' Compile failed, check your program codes please! |
新蟲(chóng) (小有名氣)
鐵桿木蟲(chóng) (職業(yè)作家)
|
參考下: Root of Mean Square Error (RMSE): 1512.6424091951 Sum of Squared Residual: 54914089.394293 Correlation Coef. (R): 0.974780606831772 R-Square: 0.950197231455318 Adjusted R-Square: 1.17033437542122 Determination Coef. (DC): 0.950332931884681 F-Statistic: -35.7233854532568 Parameter Best Estimate -------------------- ------------- a1 3.00320288098963E-5 b1 -6.17282059059677 c1 5688.58025090443 d1 -127407.822744204 a2 4.67944748234004E-6 b2 13.236089621059 c2 -7474.57974538798 d2 963214.59538852 a3 1.04462488285875E-6 c3 -2496.32779108848 |
禁蟲(chóng) (正式寫(xiě)手)
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 282分材料專(zhuān)業(yè)求調(diào)劑院校 +10 | 楓橋ZL 2026-03-09 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑,307分 +8 | 張泳銘1 2026-03-09 | 8/400 |
|
|
[考研] 320材料與化工,求調(diào)劑 +11 | 鶴遨予卿 2026-03-04 | 13/650 |
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 301求調(diào)劑 +11 | 朝天椒的雙馬尾 2026-03-05 | 11/550 |
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
[考研] 材料專(zhuān)碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
[考研] 化學(xué) 0703求調(diào)劑 總分293 一志愿211 +9 | 土土小蟲(chóng) 2026-03-03 | 11/550 |
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
[考研] 0703化學(xué)求調(diào)劑,總分320分,一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|
[考研]
|
旅行中的紫葡萄 2026-03-03 | 4/200 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研] 化學(xué)0703求調(diào)劑 學(xué)碩 理/工科均可 總分279 +3 | 1一11 2026-03-03 | 5/250 |
|