Matlab中wavedec使用學習及詳解

from http://m.blog.csdn.net/blog/u010060391/42709317

2015-1-14阅读412 评论0

最近在使用對ECG訊號進行小波係數的提取，所以用到了wavedec這個matlab中内置函數。然後盡心學習，發現還挺高深，不是我學到的衹是一點皮毛，因爲要對一段波形進行特徵提取，所以我要看哪個小波係數最能代表波形特徵，然後進行保留。

我使用的是：[C,L] = WAVEDEC(X,N,\'wname\') 

這個會返回信號X在level N（小波分解級數）的小波分解係數。

輸出爲兩部分，第一部分C爲小波分解係數（是一個矩陣形式），而L則是表示相應小波係數的個數（矩陣内係數的個數）。

我使用的level爲4 ” [C,L] = wavedec(t_signal, 4, \'db8\')“

具體分解如圖所示：

生成的矩陣C爲：

而生成的L則相應的cA4,cD4.....cD1中係數的個數。

盜一張圖：

from:http://www.mathworks.com/help/wavelet/ref/wavedec.html

更具體的算法如圖所示：

首先訊號進來，進行高低頻過濾，主要是通過卷積的方式，然後在進行降采樣，大概就是將上一步所生成樣本進行篩選。通過降采樣就可以得到近似係數(approximation coefficient)和細節係數(detail coefficient)。下圖是我的信號進行一步步的每次分解后圖形的顯示，爲了方便直觀感受，我將每一步生成的小波係數分別進行顯示輸出，分別是N=1,2,3,4的時候進行輸出：

因爲這些生成的圖形的縱坐標不同，所以我將cD1,cD2,cD3,cD4,cA1輸出在同一張坐標軸上：

顯然cA1, cD3, cD4的係數更爲明顯，更爲代表一個波形，所以要對這三個小波係數進行保留。

下一步要把這個計算挪到C上面進行計算，正在尋找相應的源代碼。

有問題，大家可以一起交流。