| 5 | 1/1 | 返回列表 |
| 查看: 1205 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問題,作者唯有千尋將贈送您 100 個(gè)金幣 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
唯有千尋新蟲 (小有名氣)
|
[求助]
1stOpt多參數(shù)約束擬合代碼錯(cuò)誤 已有1人參與
|
||
求助1stopt多參數(shù)非線性代碼錯(cuò)誤,去掉保證擴(kuò)散系數(shù)大于0的兩行代碼可以運(yùn)行,得出的解和文獻(xiàn)結(jié)果相差甚遠(yuǎn),加入兩行約束條件程序報(bào)錯(cuò),請大神指點(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! |
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
鐵桿木蟲 (職業(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 |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] A區(qū)一本交叉課題組,低分調(diào)劑,招收機(jī)械電子信息通信等交叉方向 +39 | lisimayy 2026-03-04 | 52/2600 |
|
|---|---|---|---|---|
|
[考研] 中科大材料299求調(diào)劑 +8 | DAIjiayo 2026-03-05 | 13/650 |
|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué),0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | +5 | 橙子cyx 2026-03-06 | 6/300 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考研] 324求調(diào)劑 +5 | wxz2 2026-03-03 | 8/400 |
|
|
[考研] 0805總分292,求調(diào)劑 +13 | 幻想之殤 2026-03-01 | 13/650 |
|
|
[考研] 材料與化工328求調(diào)劑 +4 | 。,。,。,。i 2026-03-02 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
[考研] 083000,總分284,求調(diào)劑 +3 | 徐yr 2026-03-04 | 3/150 |
|
|
[考研] 274求調(diào)劑 +9 | 一個(gè)學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[基金申請] 請問大家,研究風(fēng)險(xiǎn)與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 268求調(diào)劑 +10 | 簡單點(diǎn)0 2026-03-02 | 14/700 |
|