| 查看: 812 | 回復(fù): 11 | ||
[求助]
擬合問(wèn)題 已有3人參與
|
||
| 已知一個(gè)函數(shù)去擬合一系列點(diǎn);這個(gè)函數(shù)以e為基底的指數(shù)(并非簡(jiǎn)單的多項(xiàng)式),一些列點(diǎn)為(x1,x2,y),有10組,函數(shù)中有未知參數(shù)三個(gè),怎么編程求這些未知參數(shù)? |
木蟲(chóng) (正式寫(xiě)手)
專(zhuān)家顧問(wèn) (正式寫(xiě)手)

木蟲(chóng) (正式寫(xiě)手)
最小二乘法,對(duì)于指數(shù)函數(shù),先待定系數(shù),然后與實(shí)驗(yàn)對(duì)應(yīng)值求誤差,對(duì)所有求得的誤差,做平方和,得到關(guān)于待定參數(shù)的函數(shù),然后該函數(shù)分別對(duì)參數(shù)求偏導(dǎo)數(shù),令偏導(dǎo)數(shù)等于零,求得待定系數(shù),進(jìn)而知道函數(shù)表達(dá)式。樓主可以網(wǎng)上或者找個(gè)課本看看具體的原理。我只能幫你到這了 |
至尊木蟲(chóng) (文壇精英)
|
你這個(gè)其實(shí)是系數(shù)擬合,你可以用matlab自己編程,也可以使用1stopt,matlab使用nlinfit函數(shù),但是它是依賴(lài)于初值的選擇的,由于你沒(méi)有給出函數(shù)形式和數(shù)據(jù),下面給你一個(gè)例子,你自己參考吧 首先編寫(xiě)一個(gè)函數(shù),其對(duì)應(yīng)于你的已知函數(shù),不知你的函數(shù)有幾個(gè)變量幾個(gè)系數(shù),下面給出的是1個(gè)變量,2個(gè)系數(shù) function y=fxnh(beta,x) a=beta(1); b=beta(2); x1=x; y=beta(1).*(6/pi^0.5*x1.^0.5*beta(2)-3*x1*beta(2)^2);%此處為函數(shù)形式 然后再編寫(xiě)一個(gè)腳本進(jìn)行運(yùn)算 x=xlsread('2.xlsx');%變量數(shù)據(jù) y=xlsread('1.xlsx');%因變量數(shù)據(jù) beta0=[1,0.0001]';%兩個(gè)系數(shù)的初值 [b,r]=nlinfit(x,y,'fxnh',beta0); |

至尊木蟲(chóng) (文壇精英)
|
你可以自己參考一下我給的代碼,也可以使用1stopt,不過(guò)我估計(jì)你自編的可能性不大,你把程序修改為如下,一個(gè)函數(shù) function y=fxnh(beta,x) a=beta(1); b=beta(2); c=beta(3); x1=x(:,1); x2=x(:,2); y=beta(1).*exp(beta(2).*x1).*x2.^beta(3);%此處為函數(shù)形式 注意x是所有變量的總的數(shù)據(jù),不同變量按列區(qū)分 然后再編寫(xiě)一個(gè)腳本進(jìn)行運(yùn)算 x=xlsread('2.xlsx');%變量數(shù)據(jù) y=xlsread('1.xlsx');%因變量數(shù)據(jù) beta0=[1,1,1]';%三個(gè)系數(shù)的初值 [b,r]=nlinfit(x,y,'fxnh',beta0); 注意x,y都是列向量 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 材料專(zhuān)碩274一志愿陜西師范大學(xué)求調(diào)劑 +8 | 薛云鵬 2026-03-13 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +15 | 大米飯! 2026-03-15 | 18/900 |
|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 機(jī)械專(zhuān)碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化,一志愿:哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|