| 4 | 1/1 | 返回列表 |
| 查看: 1173 | 回復(fù): 3 | ||||
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 實測值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 |
鐵桿木蟲 (職業(yè)作家)
銀蟲 (小有名氣)
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +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 這一項應(yīng)當(dāng)是非常接近于0的數(shù),也就是說 exp(14.38786/T/x)這一項應(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 實測值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é)果。 |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評估A類院;瘜W(xué)工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[基金申請] 被我言中:新模板不強調(diào)格式了,假專家開始管格式了 +3 | beefly 2026-03-14 | 3/150 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[文學(xué)芳草園] 伙伴們,祝我生日快樂吧 +17 | myrtle 2026-03-10 | 26/1300 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 255求調(diào)劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | Dizzy123@ 2026-03-10 | 3/150 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +5 | 瑪卡巴卡啊哈 2026-03-11 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|