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

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

至尊木蟲 (文壇精英)
|
你可以自己參考一下我給的代碼,也可以使用1stopt,不過我估計你自編的可能性不大,你把程序修改為如下,一個函數(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ū)分 然后再編寫一個腳本進(jìn)行運算 x=xlsread('2.xlsx');%變量數(shù)據(jù) y=xlsread('1.xlsx');%因變量數(shù)據(jù) beta0=[1,1,1]';%三個系數(shù)的初值 [b,r]=nlinfit(x,y,'fxnh',beta0); 注意x,y都是列向量 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ,六級已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 321求調(diào)劑 +8 | 何潤采123 2026-03-18 | 10/500 |
|
|
[考研]
|
.6lL 2026-03-18 | 6/300 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 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)格式了,假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|