| 4 | 1/1 | 返回列表 |
| 查看: 1174 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +4 | 呼呼呼。。。。 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +7 | z1z2z3879 2026-03-16 | 9/450 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +9 | 負(fù)心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
[考研] 【0856】化學(xué)工程(085602)313 分,本科學(xué)科評估A類院;瘜W(xué)工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 10/500 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 285 求調(diào)劑 資源與環(huán)境 一志愿北京化工大學(xué) +3 | 未名考生 2026-03-10 | 3/150 |
|
|
[考研] 求調(diào)劑(材料與化工327) +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
[考研] 329求調(diào)劑 +3 | miaodesi 2026-03-12 | 4/200 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|