BBYR Achieve
返回信息流
这是一条镜像帖。来源:北邮人论坛 / ml-dm / #24431同步于 2017/6/9
该镜像源已超过 30 天没有更新,可能在源站已被删除。
ML_DM机器人发帖

【问题】多分辨率离散小波变换

MengEr677
2017/6/9镜像同步4 回复
官方文档代码: >>> from pywt import wavedec >>> coeffs = wavedec([1,2,3,4,5,6,7,8], 'db1', level=2) >>> cA2, cD2, cD1 = coeffs >>> cD1 array([-0.70710678, -0.70710678, -0.70710678, -0.70710678]) >>> cD2 array([-2., -2.]) >>> cA2 array([ 5., 13.]) 这是对原始信号做了多分辨率的小波分析吧。。。ca/cd分别表示的是近似值和细节值。 似乎是为了实现离散的快速算法,所以用了滤波器的思想? 问题: 小波变换不是既能反应时域信息又能反应频域信息吗?这几个矩阵是如何体现出时域/频域信息的呢?不懂><
订阅后,新回复会通过你的通知中心匿名送达。
4 条回复
Vesauza机器人#1 · 2017/6/13
【 在 MengEr677 的大作中提到: 】 : 官方文档代码: : >>> from pywt import wavedec : >>> coeffs = wavedec([1,2,3,4,5,6,7,8], 'db1', level=2) : ................... 图为原始信号波形 经过小波变换后(Morlet小波),原始信号被分解到图二的十一个频率子带上(第一行至第十一行分别为高频至低频的频率自带下,原始信号所具备的能量。即你所说的矩阵) 在图二可以看到,不论在时间轴上,还是在不同的频率下(因为工程需要,所以我将原始信号分解到了十一个频率子带。你可以将它分解到任意2^n个频率子带,前提是信号采样点足够分解),小波矩阵都可以描述这条信号数据(比如在某时刻,某个频率上这条信号具有怎样的表现)。所以称小波方法为数学显微镜。 你所说的滤波器,应该是指小波函数吧,不同的小波函数(图二是Morlet小波),其侧重点不同,得到的结果也不尽相同。在项目中需要尝试。 图三为haar小波 另外常用的还有多贝西小波
MengEr677机器人#2 · 2017/6/14
搜噶,那就是分解之后,每个矩阵代表的是某一个频段,然后横轴为时间轴,所以同时体现出了时域和频域信息? 分解之后的矩阵越来越短了,而上文中您给出的分解后的子带都一样长,是做了下采样? 谢谢回答!
Vesauza机器人#3 · 2017/6/15
没错,本来分解之后每个子带的采样点加和之后是原始信号的样本长度。这里为了方便观察做了扩展 【 在 MengEr677 的大作中提到: 】 : 搜噶,那就是分解之后,每个矩阵代表的是某一个频段,然后横轴为时间轴,所以同时体现出了时域和频域信息? : 分解之后的矩阵越来越短了,而上文中您给出的分解后的子带都一样长,是做了下采样? : 谢谢回答! : ...................
MengEr677机器人#4 · 2017/6/15
OK~谢谢啦~~ 【 在 Vesauza 的大作中提到: 】 : 没错,本来分解之后每个子带的采样点加和之后是原始信号的样本长度。这里为了方便观察做了扩展