| 查看: 1645 | 回復(fù): 10 | ||
| 【懸賞金幣】回答本帖問題,作者天藍(lán)的將贈(zèng)送您 90 個(gè)金幣 | ||
[求助]
1stopt軟件非線性擬合求助 已有2人參與
|
||
![]() Parameter C[,100], E[5000,], F=88[5,];Variable x, y; Function y=C*(1-(165.289*(x/(F*exp(E/6428.74225)))^0.5)*(185.185*(x/(F*exp(E/6428.74225)))^0.5)*(140.252*(x/(F*exp(E/6428.74225)))^0.5)); Data; //x y 0 0 10 1.86 22 2.85 34 3.49 46 3.99 58 4.36 71.5 4.75 83 4.96 107 5.35 131 5.79 167 6.32 191 6.63 200 6.73 擬合結(jié)果 Iterations: 35 Elapsed Time (Hr:Min:Sec:Msec): 00:00:14:105 Algorithms: Universal Global Optimation(UGO1) Stop Reason: Convergence tolerance reached Root of Mean Square Error (RMSE): 1.89309282451297 Residual Sum of Squares (SSE): 46.5894057488923 Correlation Coef. (R): 0 R-Square: 0 Determination Coef. (DC): 1.52511654600143E-16 Chi-Square: 10.6288164821448 F-Statistic: 5 Parameters Best Estimate ---------- ------------- c 4.38331076909242 e 200415.247116993 f 847455.285339582 ====== Output Results ===== No Observed y Calculated y 1 0 4.38331076909242 2 1.86 4.38331076909242 3 2.85 4.38331076909242 4 3.49 4.38331076909242 5 3.99 4.38331076909242 6 4.36 4.38331076909242 7 4.75 4.38331076909242 8 4.96 4.38331076909242 9 5.35 4.38331076909242 10 5.79 4.38331076909242 11 6.32 4.38331076909242 12 6.63 4.38331076909242 13 6.73 4.38331076909242 目標(biāo)根據(jù)下面公式擬合出拋物線,但是擬合曲線是一條直線(紅色部分);一開始我沒有添加參數(shù)范圍,可是結(jié)果還是一條直線,后來加了大概的參數(shù)范圍,結(jié)果還是直線,希望1stopt軟件非線性擬合的大神指導(dǎo)一下,謝謝 目標(biāo)公式和擬合圖 見附件圖片 ![]() |
|
目標(biāo)公式和擬合圖片在下面分享百度網(wǎng)盤里面 鏈接: https://pan.baidu.com/s/1BlEsbyaUPfhmEyxJb2bsgQ 提取碼: my4u 復(fù)制這段內(nèi)容后打開百度網(wǎng)盤手機(jī)App,操作更方便哦 |
版主 (知名作家)

版主 (知名作家)
|
不約束范圍也是一樣的結(jié)果。 檢查一下模型是否輸入有誤 發(fā)自小木蟲Android客戶端 |

鐵桿木蟲 (職業(yè)作家)
|
樓主公式應(yīng)該不對(duì)吧: 原公式:y=C*(1-(165.289*(x/(F*exp(E/6428.74225)))^0.5)*(185.185*(x/(F*exp(E/6428.74225)))^0.5)*(140.252*(x/(F*exp(E/6428.74225)))^0.5)); 等價(jià)于:y=C*(1-(165.289*185.185*140.252*(x/(F*exp(E1)))^1.5)); “F*exp(E1)”可視為一個(gè)參數(shù)(否則無窮多組解),因此該公式又可以簡化為:y=C-a*x^1.5; 效果能好嗎?還不如一個(gè)簡單的如:y=a*x^b; |
|
謝謝,仔細(xì)檢查發(fā)現(xiàn),公式輸入錯(cuò)誤,正確應(yīng)該是y=C*(1-(1-165.289*Sqrt(x/F)*exp(-E/12857.4845))*(1-185.185*Sqrt(x/F)*exp(-E/12857.4845))*(1-140.252*Sqrt(x/F)*exp(-E/12857.4845)))。還有一個(gè)問題,現(xiàn)在擬合是拋物線了,但是單擊1次擬合的結(jié)果,和單擊2次到的擬合結(jié)果,為啥參數(shù)差別很大呀?相關(guān)系數(shù)多是0.999,這是說明有好幾組最優(yōu)解嗎?熱執(zhí)行和執(zhí)行2個(gè)菜單有啥區(qū)別呀? |
送紅花一朵 |
謝謝,我的模型參數(shù)確實(shí)輸入不正確,正常的話,應(yīng)該是y=C*(1-(1-165.289*Sqrt(x/F)*exp(-E/12857.4845))*(1-185.185*Sqrt(x/F)*exp(-E/12857.4845))*(1-140.252*Sqrt(x/F)*exp(-E/12857.4845))),這樣就不是你說的過擬合了。還有一個(gè)問題,現(xiàn)在擬合是拋物線了,但是單擊1次擬合的結(jié)果,和單擊2次到的擬合結(jié)果,為啥參數(shù)差別很大呀?相關(guān)系數(shù)多是0.999,這是說明有好幾組最優(yōu)解嗎?熱執(zhí)行和執(zhí)行2個(gè)菜單有啥區(qū)別呀? |
鐵桿木蟲 (職業(yè)作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 298 +8 | 小西笑嘻嘻 2026-03-03 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +5 | 1一11 2026-03-03 | 5/250 |
|
|
[考研] 江蘇省農(nóng)科院招調(diào)劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 一志愿東北大學(xué)化學(xué)314分求調(diào)劑 +5 | lr1212.. 2026-03-02 | 5/250 |
|
|
[考研] 化工專碩348,一志愿985求調(diào)劑 +8 | 弗格個(gè) 2026-02-28 | 11/550 |
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 考研復(fù)試調(diào)劑,過國家線的同學(xué)都可報(bào)名 +7 | 黑!在干嘛 2026-02-28 | 8/400 |
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 295求調(diào)劑。一志愿報(bào)考鄭州大學(xué)化學(xué)工藝學(xué)碩,總分295分 +8 | yl1 2026-03-02 | 9/450 |
|
|
[考研] 一志愿中石油(華東)本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
[考研] 一志愿華南理工大學(xué)材料與化工326分,求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
[考研] 0856化工專碩求調(diào)劑 +15 | 董boxing 2026-03-01 | 15/750 |
|
|
[考研] 292求調(diào)劑 +7 | yhk_819 2026-02-28 | 7/350 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
[考研] 313求調(diào)劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|