| 4 | 1/1 | 返回列表 |
| 查看: 1204 | 回復: 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為離子導電率,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ā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料求調(diào)劑 +15 | hyf hyf hyf 2026-02-28 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境調(diào)劑 +7 | 柒槿levana 2026-03-01 | 7/350 |
|
|
[考研] 070300化學 280 一志愿太原理工 求調(diào)劑 +7 | 拾玖壹 2026-03-04 | 7/350 |
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
[考研] 295求調(diào)劑 +4 | 小賽不吃香菜 2026-03-04 | 4/200 |
|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 材料專碩290求調(diào)劑 +3 | 杰尼龜aaa 2026-03-04 | 3/150 |
|
|
[考研] 歡迎采礦、地質(zhì)、巖土、計算機、人工智能等專業(yè)的同學報考 +9 | pin8023 2026-02-28 | 12/600 |
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
[考研] 085700資環(huán)求調(diào)劑,初始279,六級已過,英語能力強 +3 | 085700資環(huán)調(diào)劑 2026-03-03 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ErMiao1020 2026-03-02 | 9/450 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 化工335求調(diào)劑 +5 | 摸摸貓貓頭 2026-03-02 | 5/250 |
|
|
[考研] 考研復試調(diào)劑,過國家線的同學都可報名 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 303求調(diào)劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
[考研] 328求調(diào)劑 +3 | aaadim 2026-03-01 | 5/250 |
|
|
[論文投稿]
求助coordination chemistry reviews 的寫作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|