作者YoursEver (beyond the time!)
看板MATLAB
标题Re: [讨论] 小波(Wavelet)转换,以及FFT的问题
时间Sat Oct 19 08:19:11 2019
※ 引述《abcc50717 (鬼彻妖刀)》之铭言:
: 小弟目前接触到小波理论,
: 使用到了"离散小波转换" 使用的程式码大致如下:
: 右图是原始讯号的fft,得到在50跟950Hz的地方有值,表示转换是正确的。
: 再来看到左图跟中间的频谱图,在(500+50)Hz以及(500-50)Hz的部分竟然有值,
: 但这是原始讯号并不存在的频率。
: 照理来说小波转换,还原之後的高低频讯号做FFT之後,不应该出现讯号本身没有的频率吧?
: 还是我对小波转换有什麽误会呢?
: 附上我的Matlab Code, https://reurl.cc/24LKRv GOOGLE云端
: → YoursEver: 你对DSP有误会. 10/19 08:13
Fourier{ source * wfilter} = Fourier{ source } x Fourier{ wfilter }
^ convolution ^ multiplication
所以有其他的component在,才是正常的.
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.136.129.68 (台湾)
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1571444353.A.E88.html
1F:推 abcc50717: 对欸~小波本身是一个filter 所以与source择积会改变 10/21 18:09
2F:→ abcc50717: source。 10/21 18:10
3F:→ abcc50717: 不过比较好奇的是 wavelet分解之後,高频部分含有原始 10/21 18:15
4F:→ abcc50717: 讯号不存在的频率。 但重构讯号s1=a1+d1(参考.m CODE) 10/21 18:19
5F:→ abcc50717: 刚刚发现重构讯号s1没有完全等於原始讯号s 10/21 18:24
6F:→ abcc50717: s1与s 的差异在10^-16 order以下 10/21 18:25
你要重念perfect reconstruction那一段.
※ 编辑: YoursEver (114.32.63.235 台湾), 10/22/2019 01:22:10
7F:推 abcc50717: 好的 感谢你的关键字指引~ 来研究一下 10/22 16:56