| 4 | 1/1 | 返回列表 |
| 查看: 1202 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者唯有千尋將贈送您 100 個金幣 | ||
唯有千尋新蟲 (小有名氣)
|
[求助]
1stOpt多參數(shù)約束擬合代碼錯誤 已有1人參與
|
|
求助1stopt多參數(shù)非線性代碼錯誤,去掉保證擴散系數(shù)大于0的兩行代碼可以運行,得出的解和文獻結(jié)果相差甚遠,加入兩行約束條件程序報錯,請大神指點如何加上約束條件能正常運行 代碼如下:Title "Nonlinear Fit for Defect diffusion coefficient"; Constant F=96485.334,R=8.3145,b3=0,d3=0;//法拉第常數(shù)F,R氣體常數(shù),文獻中的兩個參數(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;//保證擴散系數(shù)大于0; a2*(1+b2*x4)>0;//保證擴散系數(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; 程序報錯: Incompatible types: 'Double' and 'Boolean' Incompatible types: 'Double' and 'Boolean' Compile failed, check your program codes please! |
鐵桿木蟲 (職業(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 |
新蟲 (小有名氣)
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
085600 英一數(shù)二272求調(diào)劑
5+6
|
vida_a 2026-03-01 | 47/2350 |
|
|---|---|---|---|---|
|
[考研] 能動297求調(diào)劑,本科川大 +3 | 邵11 2026-03-04 | 3/150 |
|
|
[考研] 295求調(diào)劑 +5 | 等春來, 2026-03-04 | 5/250 |
|
|
[考研] 0856材料與化工,270求調(diào)劑 +12 | YXCT 2026-03-01 | 14/700 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 283求調(diào)劑 +9 | 鹿沫笙 2026-03-02 | 9/450 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[論文投稿]
通訊作者寫誰,問題是你意想不到的問題
15+3
|
阿爾法啊 2026-03-01 | 3/150 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
[考研] 材料復(fù)試調(diào)劑 +5 | 學(xué)材料的點 2026-03-01 | 6/300 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考博] 誠招農(nóng)業(yè)博士 +3 | 心欣向榮 2026-02-28 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
[基金申請]
剛錄用,沒有期刊號,但是在線可看的論文可以放為代表作嗎
10+3
|
arang1 2026-03-01 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|