| 5 | 1/1 | 返回列表 |
| 查看: 1184 | 回復(fù): 3 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||||
forestsong12銀蟲 (小有名氣)
|
[求助]
1stopt軟件擬合非線性方程,計算結(jié)果全為0,這是為什么 已有2人參與
|
||||
|
Parameters T; Variable x,y; Function y= 3.741832*10^-16/(exp(0.01438786/(T*x*10^-3))-1)/(x*10^-3)^5; Data; 405.425 19197 424.055 14239 441.669 11853 463.348 6832.2 541.255 4912.8 566.321 3806.9 609.001 2420.1 711.297 1356.1 744.153 810.13 764.477 613.14 804.786 412.09 852.208 315.52 900.307 463.06 迭代數(shù): 33 計算用時(時:分:秒:微秒): 00:00:00:319 優(yōu)化算法: 麥夸特法(Levenberg-Marquardt) + 通用全局優(yōu)化法 計算結(jié)束原因: 達(dá)到收斂判斷標(biāo)準(zhǔn) 均方差(RMSE): 7875.40731737143 殘差平方和(SSE): 806286525.388597 相關(guān)系數(shù)(R): 0.985793202644828 相關(guān)系數(shù)之平方(R^2): 0.971788238380747 決定系數(shù)(DC): -0.75817360362487 卡方系數(shù)(Chi-Square): 67231.0399999996 F統(tǒng)計(F-Statistic): 378.908299540322 參數(shù) 最佳估算 ---------- ------------- T 27.2788094020072 ====== 結(jié)果輸出 ===== No 實(shí)測值y 計算值y 1 19197 0.0000000 2 14239 0.0000000 3 11853 0.0000000 4 6832.2 0.0000000 5 4912.8 0.0000000 6 3806.9 0.0000000 7 2420.1 0.0000000 8 1356.1 0.0000000 9 810.13 0.0000000 10 613.14 0.0000000 11 412.09 0.0000000 12 315.52 0.0000000 13 463.06 0.0000000 |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
公式本身不合適。 先來看你的數(shù)據(jù) x y 405.425 19197 424.055 14239 441.669 11853 463.348 6832.2 541.255 4912.8 566.321 3806.9 609.001 2420.1 711.297 1356.1 744.153 810.13 764.477 613.14 804.786 412.09 852.208 315.52 900.307 463.06 你給出方程簡化得: y=3.741832*0.1*(x^(-5)) / (exp(14.38786/T/x)-1) 注意x是幾百數(shù)量級的,那么式中3.741832*0.1*(x^(-5))將是非常小的數(shù),至少小于10^-11,4^5與9^5相差50倍;而y是數(shù)百到幾萬的數(shù)值,那么顯然 exp(14.38786/T/x)-1 這一項(xiàng)應(yīng)當(dāng)是非常接近于0的數(shù),也就是說 exp(14.38786/T/x)這一項(xiàng)應(yīng)當(dāng)非常接近于1,這要求T非常大;如前述,3.741832*0.1*(x^(-5))在x整個變化范圍中數(shù)值變化達(dá)50倍,而至于指數(shù)函數(shù)而言,其變化范圍是很窄的。 也就是,不可能存在合適的T值能滿足高擬合的質(zhì)量。 至于你說y的計算值為0的問題,在1stopt算法設(shè)置中,以簡面體爬山法計算,可得: 相關(guān)系數(shù)(R): 0.985799655457293 相關(guān)系數(shù)之平方(R^2): 0.971800960699717 決定系數(shù)(DC): -0.758172807118177 卡方系數(shù)(Chi-Square): 67230.9893136268 F統(tǒng)計(F-Statistic): 379.084211127347 參數(shù) 最佳估算 ---------- ------------- T 117315023909293 ====== 結(jié)果輸出 ===== No 實(shí)測值y 計算值y 1 19197 0.0038974 2 14239 0.0032564 3 11853 0.0027671 4 6832.2 0.0022845 5 4912.8 0.0012269 6 3806.9 0.0010237 7 2420.1 0.0007655 8 1356.1 0.0004113 9 810.13 0.0003434 10 613.14 0.0003083 11 412.09 0.0002510 12 315.52 0.0001996 13 463.06 0.0001603 從結(jié)果中可見,T的估算值非常大。 簡而概之,就是公式不合適。此外,在擬合的時候,公式輸入能化簡就化簡,有利于避免諸如過擬合之類的問題,從而得到高質(zhì)量的擬合結(jié)果。 |

鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
孅華 2026-03-22 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 11/550 |
|
|
[考研]
求調(diào)劑材料學(xué)碩080500,總分289分
5+3
|
@taotao 2026-03-19 | 21/1050 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +5 | 苑豆豆 2026-03-20 | 5/250 |
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項(xiàng)目,四級已過 +6 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
[考研] 尋找調(diào)劑 +4 | 倔強(qiáng)芒? 2026-03-21 | 4/200 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[基金申請] 山東省面上項(xiàng)目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|