| 5 | 1/1 | 返回列表 |
| 查看: 2543 | 回復(fù): 7 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
alpshaode銀蟲 (小有名氣)
|
[求助]
求助一個(gè)最小二乘法matlab擬合程序。
|
||
| 求助一個(gè)用最小二乘法matlab擬合程序。函數(shù)是這樣的,y=exp(-x*1.22/(x+(1-x)*a)+(1-x)*1.60/(x*b+(1-x))-7.91), a,b為所要求得的擬合參數(shù),以及求出平均平方誤差。太感謝了,請給出具體的matlab程序。 |
taotie |

|
x=[0.31748 0.63496 0.95244 1.26992 1.5874 1.90488 2.22236 2.53984 2.85732 3.1748 3.9685 4.7622 5.5559 6.3496]; y=[72.70662369 73.31157463 73.9514314 74.62277256 75.32217667 76.0462223 76.79148803 77.55455241 78.33199401 79.1203914 81.11686331 83.10696505 85.03723673 86.85421846]; P0=[1 1 1 1 1 1]'; [P,r,J] = nlinfit(x',y',@myfun,P0); %得到P參數(shù),r殘差residual,JJacobian矩陣 [YY,delta] = nlpredci('myfun',x',P,r,J); plot(x,y,'k+',x,YY,'r') function yd = myfun(P,x) % myfun=@(x,P)(-P(1).*log(x)-P(1)*(1.023*x.^0.5)./(1+x.^0.5)+P(2)*x+P(3)*x.^1.5+P(4)*x.^2+P(5)*x.^2.5+P(6)*x.^3+72.14); yd=-P(1).*log(x)-P(1)*(1.023*x.^0.5)./(1+x.^0.5)+P(2)*x+P(3)*x.^1.5+P(4)*x.^2+P(5)*x.^2.5+P(6)*x.^3+72.14; end 程序解釋:非線性擬合,下面的myfun是要擬合的公式,待估計(jì)的參數(shù)為P,有幾個(gè)就p(幾)。P0是初始值,先隨意設(shè)置,估計(jì)一次后,用估計(jì)值在估計(jì)一遍。 x,y為數(shù)據(jù)。 多少給點(diǎn)分吧 |

銀蟲 (小有名氣)

金蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿浙江大學(xué)工科動(dòng)力工程370,數(shù)一121,專業(yè)課135,現(xiàn)在能去哪里 +3 | 080700調(diào)劑 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +3 | Yam. 2026-03-30 | 4/200 |
|
|
[考研] 08工科求調(diào)劑286 +5 | tgs_001 2026-03-28 | 5/250 |
|
|
[考研] 抱歉 +4 | 田洪有 2026-03-30 | 4/200 |
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
[考研] 各位老師好,我的一志愿為北京科技大學(xué)085601材料專碩 +10 | Koxui 2026-03-28 | 10/500 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
[考研] 283求調(diào)劑(080500) +14 | A child 2026-03-27 | 14/700 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 318一志愿吉林大學(xué)生物與醫(yī)藥 求調(diào)劑 +5 | 篤行致遠(yuǎn). 2026-03-28 | 5/250 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +5 | hanamiko 2026-03-29 | 5/250 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。。 +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4 | 小魚愛有機(jī) 2026-03-25 | 4/200 |
|
|
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3 | 江往賣魚 2026-03-26 | 3/150 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|