大家好,又见面了,我是你们的朋友全栈君。
最近在使用對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上面進行計算,正在尋找相應的源代碼。
有問題,大家可以一起交流。
发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/163787.html原文链接:https://javaforall.cn
【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛
【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...