matlab中wavedec2,Matlab中wavedec使用學習及詳解 | 学步园

matlab中wavedec2,Matlab中wavedec使用學習及詳解 | 学步园最近在使用對ECG訊號進行小波係數的提取,所以用到了wavedec這個matlab中内置函數。然後盡心學習,發現還挺高深,不是我學到的衹是一點皮毛,因爲要對一段波形進行特徵提取,所以我要看哪個小波係數最能代表波形特徵,然後進行保留。我使用的是:[C,L]=WAVEDEC(X,N,’wname’)這個會返回信號X在levelN(小波分解級數)的小波分解係數。輸出爲兩部分,第一部分C爲小波分解係…

大家好,又见面了,我是你们的朋友全栈君。

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

我使用的是:[C,L] = WAVEDEC(X,N,’wname’)

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

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

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

具體分解如圖所示:

3235440.html

c96dee36d1de517a8966c519c0890a70.png

生成的矩陣C爲:

d054306bcc6b918329d3033edca2491d.png

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

盜一張圖:

3235440.html

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

bed006aeb25f54bcaef4c2844839c1a7.png

更具體的

da7fb3911d186ed012629d88a0b2ddc4.png

3235440.html

首先訊號進來,進行高低頻過濾,主要是通過卷積的方式,然後在進行降采樣,大概就是將上一步所生成樣本進行篩選。通過降采樣就可以得到近似係數(approximation coefficient)和細節係數(detail coefficient)。

下圖是我的信號進行一步步的每次分解后圖形的顯示,爲了方便直觀感受,我將每一步生成的小波係數分別進行顯示輸出,分別是N=1,2,3,4的時候進行輸出:

d30391a971ab640fe0a1bd077791b10d.png

3235440.html

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

3235440.html

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

a0be33d1ef922f9393eeabdc2fad924b.png

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

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

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/163787.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • [GBA ROM列表]不断补完中……

    [GBA ROM列表]不断补完中……出处:http://www.ezflash.cn/vbb/showthread.php?s=3a42901b8c2aacb56a02e61f690b8a05&threadid=14&perpage=15&pagenumber=10001-F-Zero(J)(最高时速)(J)RCG320002-SuperMarioAdvance(J)(超级马里奥A)(J)ACT32000

    2022年7月11日
    77
  • 织梦dedecms如何去除版权中的Power by DedeCms

    织梦dedecms如何去除版权中的Power by DedeCms

    2021年9月19日
    47
  • pytest重试_arcmap重分类失败

    pytest重试_arcmap重分类失败安装:pip3installpytest-rerunfailures重新运行所有失败用例要重新运行所有测试失败的用例,请使用–reruns命令行选项,并指定要运行测试的最大次数:$py

    2022年7月28日
    4
  • ZIP 压缩解压命令

    ZIP 压缩解压命令压缩一个文件,打开cmd窗口 zip f:\wenjian.zip f:\wenjian.txt压缩一种格式的文件 zip f:\wenjian.zip f:\*.txt压缩文件夹及里面的内容 -r表示压缩文件夹wenjian及里面的内容 zip-rf:\wenjian.zip wenjian  移除压缩包里的一

    2022年5月16日
    53
  • js字符串操作方法(js对象转字符串)

    一、关于字符串分割1、slice(start,end);关于这个方法,一定要搞懂四个关键点:(1)截取字符串时不包括下标为end的元素。(2)end是可选参数,没有时,默认从start到结束的所有字符串。(3)String.slice与Array.slice区别。(4)参数为负数时,是如何处理的。其中第3点其实就是在JavaScript中字符串和数组都具有这个方法,它们…

    2022年4月18日
    51
  • DDPG存疑

    DDPG存疑存疑的点:value函数的参数是(s,a)。离散状态下s,a作为横纵坐标来获取Q值。q-learningDQN情况下,输出的是最大Q值。PG里面输出的是action的概率分布。DDPG里,actor网络输出的是Action。critic网络里输出的是Q值,依然和s,a有关,由s,a确定。只不过这里确定的方式是①将s,a分别经过一个输出维度为30的网络后 得到的值 相加(两个网络分别为:(s_dim,30)(a_dim,30)。②将s,a维度拼接再经过网络计算(这个网络的维度是(s_d

    2022年6月16日
    25

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号