| 5 | 3/1 | 返回列表 |
| 查看: 23814 | 回復(fù): 73 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
dshuanger木蟲 (小有名氣)
|
[求助]
matlab非線性擬合模型參數(shù),大神你在哪!
|
|||
|
本人非計算模擬專業(yè),小白一枚,但最近需要模擬參數(shù),不知如何運用matlab求參數(shù),顧求助各位大神。 模型表達式: y=a/E1 +a/E2[1-exp(-x/t)]+a*x/E3 其中:x=[2,5,10,15,30,60,120,200,300,550,700]; y=[8.66667,10,10.66667,12.66667,13.33333,15,15.33333, 16.66667,17.33333,18.33333,20]; a=2 未知參數(shù):E1,E2,t,E3 不知道能否通過matlab擬合參數(shù) 請各位大神能鼎力相助,最好能給出代碼或者計算方法,謝謝大家! [ Last edited by dshuanger on 2013-11-5 at 11:06 ] |
|
我自己實驗數(shù)據(jù)的擬合代碼: Parameters b0,b1,b2,b3,b4,b5,b6,b7,b8,b9; Variable X1,X2,X3,Y; Function Y=b0+b1*X1+b2*X2+b3*X3+b4*X1*X2+b5*X1*X3+b6*X2*X3+b7*X1^2+b8*X2^2+b9*X3^2; Data; 0.0000 8.0000 7.0000 63.9 0.0000 24.0000 11.5000 56.9 350.0000 40.0000 5.5000 95.4 350.0000 0.0000 11.5000 54.3 700.0000 16.0000 5.5000 92.7 700.0000 32.0000 10.000 92.9000 1050.0000 48.0000 4.0 95.80 1050.0000 0.0000 10.000 87.8 1400.0000 16.0000 4.0000 92.6000 1400.0000 32.0000 8.5000 98.700 1750.0000 48.0000 2.5 71.2000 1750.0000 8.0000 8.5000 95.1000 2100.0000 24.0000 2.5000 64.8000 2100.0000 40.0000 7.0000 67.2000 擬合數(shù)據(jù)結(jié)果: 均方差(RMSE): 1.54658329060161 殘差平方和(SSE): 33.4868782467533 相關(guān)系數(shù)(R): 0.995247203455277 相關(guān)系數(shù)之平方(R^2): 0.990516995985549 決定系數(shù)(DC): 0.990516995985549 卡方系數(shù)(Chi-Square): 0.228909500020559 F統(tǒng)計(F-Statistic): 1253.42179901157 參數(shù) 最佳估算 ---------- ------------- b0 37.0459960300393 b1 0.0487926069483363 b2 1.50487196705523 b3 5.5454028094221 b4 -0.000558451632571275 b5 0.00153209146092655 b6 -0.0281388598195167 b7 -2.10419568056527E-5 b8 -0.0122539097071058 b9 -0.492832403832336 ====== 結(jié)果輸出 ===== No 實測值y 計算值y 1 63.9 61.3939773 2 56.9 56.9333928 3 95.4 96.6663312 4 54.3 56.3074838 5 92.7 94.5904221 6 92.9 91.8798377 7 95.8 96.2627759 8 87.8 87.3372240 9 92.6 93.6201624 10 98.7 96.8095779 11 71.2 69.1925163 12 95.1 93.8336688 13 64.8 64.7666071 14 67.2 69.7060227 和文獻里的不一致。文獻的二次型多項式見21樓圖1: |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +14 |
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
|
y=2/e1 +2/e2/(1-exp(-x/t))+a*x/3; // a=2 x=[2,5,10,15,30,60,120,200,300,550,700]; y=[8.66667,10,10.66667,12.66667,13.33333,15,15.33333, 16.66667,17.33333,18.33333,20]; 相關(guān)系數(shù)(R): 1 相關(guān)系數(shù)之平方(R^2): 1 決定系數(shù)(DC): 0.999999999996653 卡方系數(shù)(Chi-Square): 1.64243431014082E-11 參數(shù) 最佳估算 ---------- -------------- e1 =1.2090897803241 e2 =3.11591327892529 t =0.00209072546506921 ====== 輸出結(jié)果 ===== No. Observed y Calculated y 1 8.66667 8.66666365993411 2 10 9.99999365993411 3 10.66667 10.6666636599341 4 12.66667 12.6666636599341 5 13.33333 13.3333236599341 6 15 14.9999936599341 7 15.33333 15.3333236599341 8 16.66667 16.6666636599341 9 17.33333 17.3333236599341 10 18.33333 18.3333236599341 11 20 19.9999936599341 |
榮譽版主 (文壇精英)
![]() |
專家經(jīng)驗: +518 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 能源材料化學課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學上- 2026-03-14 | 9/450 |
|
|
[考研] 考研化學學碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 【0856】化學工程(085602)313 分,本科學科評估A類院校化學工程與工藝,誠求調(diào)劑 +7 | 小劉快快上岸 2026-03-11 | 8/400 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 285化工學碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 277材料科學與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
[考研] 工科,求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
[考研] 0703化學求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|