 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1201 | 回復: 3 | ||
| 【懸賞金幣】回答本帖問題，作者唯有千尋將贈送您 100 個金幣 | ||
[求助]
1stOpt多參數(shù)約束擬合代碼錯誤 已有1人參與
|
||
 求助1stopt多參數(shù)非線性代碼錯誤，去掉保證擴散系數(shù)大于0的兩行代碼可以運行，得出的解和文獻結果相差甚遠，加入兩行約束條件程序報錯，請大神指點如何加上約束條件能正常運行代碼如下：Title "Nonlinear Fit for Defect diffusion coefficient"; Constant F=96485.334,R=8.3145,b3=0,d3=0;//法拉第常數(shù)F，R氣體常數(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ù)預指數(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! |
回復此樓» 猜你喜歡
德國亥姆霍茲Hereon中心招收兩位醫(yī)用鎂合金腐蝕與LPSO相變方向2026公派博士生
已經(jīng)有0人回復
-
推薦一款可以AI輔助寫作的Latex編輯器SmartLatexEditor，超級好用，推薦試試
已經(jīng)有11人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有100人回復
-
2026-CJ開始申報了
已經(jīng)有1人回復
-
西安電子科學大學杭州研究院劉麗香教授招收智能多模態(tài)傳感器和微型儲能器件方向博士
已經(jīng)有13人回復
-
重慶交通大學光子學微結構與器件課題組招收2026年碩士研究生信息
已經(jīng)有1人回復
-
一志愿鄭大材料學碩298分，求調(diào)劑
已經(jīng)有6人回復
-
尋合作：應力腐蝕多尺度模擬
已經(jīng)有2人回復
dingd
鐵桿木蟲 (職業(yè)作家)
- 計算強帖: 4
- 應助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
【答案】應助回帖
★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
獨孤神宇: 金幣+5, 鼓勵交流 2021-03-08 18:56:47
感謝參與，應助指數(shù) +1
獨孤神宇: 金幣+5, 鼓勵交流 2021-03-08 18:56:47
|
參考下： 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 |
» 本帖已獲得的紅花（最新10朵）
唯有千尋2樓2021-03-08 17:11:55
3樓2021-03-08 20:39:18
Cfdstudio2
禁蟲 (正式寫手)
|
本帖內(nèi)容被屏蔽 |
4樓2021-03-14 18:50:46
相關版塊跳轉(zhuǎn)
我要訂閱樓主 唯有千尋 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 中國林科院林化所（南京）2026年招收化學/材料/環(huán)境工程等背景碩士研究生3名 +3 | realstar2006 2026-02-27 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 化工專碩調(diào)劑 +4 | 利好利好. 2026-03-03 | 7/350 |
|
|
|
[考研] 085600材料與化工調(diào)劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè)，申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
[考研] 一志愿華中科技大學，化學專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 環(huán)境工程專碩307求調(diào)劑 +3 | ccc! 2026-03-03 | 3/150 |
|
|
|
[考研] 考研復試調(diào)劑，過國家線的同學都可報名 +7 | 黑！在干嘛 2026-02-28 | 8/400 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研] 化學，材料，環(huán)境類求調(diào)劑 +7 | 考研版棒棒 2026-03-02 | 7/350 |
|
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
|
[考研] 261求調(diào)劑 +3 | 陸lh 2026-03-01 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
|
[基金申請] 此成果不能導入原因：元數(shù)據(jù)必填信息不完整，可 進行補充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 281求調(diào)劑 +5 | 2026計算機_誠心 2026-03-01 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
|
[考研] 材料調(diào)劑 +6 | 愛擦汗的可樂冰 2026-02-28 | 7/350 |
|
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-02-28 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|