| 查看: 2838 | 回復(fù): 10 | ||
water11金蟲 (正式寫手)
|
[求助]
MATLAB曲線擬合 已有4人參與
|
|
|
請(qǐng)教大家一個(gè)用MATLAB進(jìn)行曲線擬合的問題。 現(xiàn)有100行數(shù)據(jù),每行3000個(gè)數(shù)據(jù),其實(shí)就是一個(gè)3000*100的矩陣,希望找到一個(gè)指數(shù)函數(shù)與此相擬合。我用MATLAB畫出了這100條曲線,可實(shí)在想不出對(duì)這些曲線擬合,找到一個(gè)與這些曲線相符的指數(shù)函數(shù)。請(qǐng)大家?guī)蛶兔,萬分感謝! |
銅蟲 (小有名氣)
|
我是以y=a*exp(x)+b作個(gè)例子,LZ具體的可以按自己的函數(shù)寫。只要在定義函數(shù)那里修改就好。 如果LZ的不知道函數(shù)是什么形式,我的思路是多作幾種,然后把100條線中,擬合效果最好的那個(gè)函數(shù)形式找出來 。 只作一咱函數(shù)形式的例子如下。 分三部分,生成數(shù)據(jù),擬合,輸出。 LZ是100條線,每條3000個(gè)點(diǎn)。 我的是5條線,每條10個(gè)點(diǎn)(省時(shí))。 具體的程序如下,不明白matlab的函數(shù)看幫助。 我沒說清楚的思路可以跟貼。 |

木蟲 (知名作家)
|
首先我覺得,用一條指數(shù)曲線擬合100條曲線,肯定是需要先將100行數(shù)據(jù)先求平均再進(jìn)行處理的。這點(diǎn)同意七樓觀點(diǎn)。其次,你怎么知道單一指數(shù)能擬合好?你必須對(duì)曲線進(jìn)行半對(duì)數(shù)作圖(semilogy)。如果半對(duì)數(shù)作圖的結(jié)果顯示接近一條直線才能認(rèn)為可用單指數(shù)回歸。否則,需要對(duì)數(shù)據(jù)進(jìn)行拉普拉斯逆變換(ILT)。實(shí)際可用非負(fù)約束的吉洪諾夫正則化算法來分析指數(shù)成分。樓主方便的話可以把數(shù)據(jù)發(fā)上來,我?guī)湍憧纯础?br />
[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
木蟲 (正式寫手)

木蟲 (著名寫手)

金蟲 (正式寫手)
金蟲 (正式寫手)
銅蟲 (小有名氣)
|
我的思路是將所有的數(shù)據(jù)進(jìn)行平均,然后 作一次曲線擬合。 具體的程序如下。 |

|
最好確定一個(gè)有限長(zhǎng)度的函數(shù): y=a(n)x^n+a(n-1)x^(n-1)+...a(0)x,n是有限的 然后采用最小二乘矩陣估計(jì),求出各個(gè)系數(shù),你簡(jiǎn)單計(jì)算一下矩陣方程,就很容易求出各個(gè)系數(shù)。。然后再matlab中通過矩陣求。。希望對(duì)你有幫助。。x^n級(jí)數(shù)是幾乎可以通用的,當(dāng)然你也可以擬合成傅里葉級(jí)數(shù)。。希望這個(gè)思路對(duì)你有幫助。。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 家佳佳佳佳佳 2026-03-29 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 290求調(diào)劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +9 | 景隨風(fēng) 2026-03-29 | 14/700 |
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +6 | 深V宿舍吧 2026-03-29 | 6/300 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 一志愿太原理工安全工程300分,求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
[考研] 壓國(guó)家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 291求調(diào)劑 +6 | HanBeiNingZC 2026-03-24 | 6/300 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
[考研] 0854AI CV方向招收調(diào)劑 +4 | 章小魚567 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|