| 5 | 1/1 | 返回列表 |
| 查看: 1624 | 回復: 5 | |||
| 【懸賞金幣】回答本帖問題,作者rollling將贈送您 5 個金幣 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
rollling新蟲 (初入文壇)
|
[求助]
1stOpt多參數(shù)約束擬合代碼錯誤 已有1人參與
|
||
|
大家好!關(guān)于1stOpt多參數(shù)約束擬合問題求助大家,其中參數(shù)a1,b1,c1,d1,a2,b2,c2,d2,a3,c3是待求參數(shù),變量y,x1,x2,x3,x4為實驗數(shù)據(jù),錯誤原因應(yīng)該是在Function下面三個約束條件。希望能夠有大神解決!代碼如下: Title "Defect diffusion coefficient"; Constant F=96485.334,R=8.3145,b3=0,d3=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; Variables y,x1,x2,x3,x4; 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))*exp(-(c1+d1*x4)/(R*T))>=0; a2*(1+b2*x4))*exp(-(c2+d2*x4)/(R*T))>=0; a3*(1+b3*x4))*exp(-(c3+d3*x4)/(R*T))>=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; 運行錯誤提示: error: missing operator or semicolon error: missing operator or semicolon error: missing operator or semicolon error: incompatible types: 'double' and 'boolean' D:\Program Files\1stOpt 8. Incompatible types: 'Double' and 'Boolean' |
版主 (知名作家)
|
你指的是約束y的計算范圍么? 這個8.0版本可以做到,通過循環(huán)計算約束每一組計算結(jié)果的范圍 可以看看使用手冊,里面有個這方面的例子。 發(fā)自小木蟲Android客戶端 |

版主 (知名作家)
|
將參數(shù)定義修改一下,a1=[0,1E-4] 發(fā)自小木蟲Android客戶端 |

新蟲 (初入文壇)
版主 (知名作家)
|
你值的是約束y的計算范圍么? 這個8.0版本可以做到,通過循環(huán)計算約束每一組計算結(jié)果的范圍 可以看看使用手冊,里面有個這方面的例子。 發(fā)自小木蟲Android客戶端 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +7 | zchqwer 2026-03-10 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 工科278分求調(diào)劑 +3 | 周慢熱啊 2026-03-12 | 5/250 |
|
|
[考研] 304求調(diào)劑(085602一志愿985) +9 | 化工人999 2026-03-09 | 9/450 |
|
|
[考研] 271求調(diào)劑 +7 | 生如夏花… 2026-03-11 | 7/350 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +5 | KRSLSR 2026-03-10 | 5/250 |
|
|
[考研] 工科0856專碩化學工程269能調(diào)劑嗎 +9 | 我想讀研11 2026-03-10 | 9/450 |
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽- 2026-03-09 | 18/900 |
|
|
[考研] 081200計算機科學與技術(shù)數(shù)一英一306求調(diào)劑 +3 | intankt 2026-03-08 | 3/150 |
|
|
[考研] 材料調(diào)劑,307分 +12 | 張泳銘1 2026-03-09 | 13/650 |
|
|
[考研] 327求調(diào)劑 +3 | Ffff03 2026-03-10 | 3/150 |
|
|
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4 | Tolkien 2026-03-07 | 8/400 |
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
[碩博家園] 木蟲好像不熱鬧了,是不是? +4 | 偏振片 2026-03-10 | 4/200 |
|
|
[考研] 求調(diào)劑,數(shù)一英一274分 +4 | 小菲會努力 2026-03-08 | 4/200 |
|
|
[考研] 求調(diào)劑,一志愿江南大學,食品科學與工程,總分,320 +3 | yyyyyukino 2026-03-07 | 3/150 |
|
|
[考研] 0701-322 求調(diào)劑 +3 | jiliuxian 2026-03-06 | 8/400 |
|
|
[考研] 347求調(diào)劑 +4 | 浮云滿足 2026-03-07 | 4/200 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 求調(diào)劑 +3 | snw石 2026-03-07 | 3/150 |
|
|
[考博] 2026申博自薦 六級440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|