现在的位置: 首页 > 综合 > 正文

小波变换的概念、操作

2018年10月27日 ⁄ 综合 ⁄ 共 1488字 ⁄ 字号 评论关闭

from:http://bbs.matwav.com/?fromuid=500164

gjsdgjsd  小波变换版主

1。连续小波的概念。就是把一个可以称作小波的函数(从负无穷到正无穷积分为零)在某个尺度下与待处理信号卷积。改变小波函数的尺度,也就改变了滤波器的带通范围,相应每一尺度下的小波系数也就反映了对应通带的信息。本质上,连续小波也就是一组可控制通带范围的多尺度滤波器。

2。连续小波是尺度可连续取值的小波,里面的a一般取整数,而不像二进小波a取2的整数幂。从连续小波到二进小波再到正交离散小波,其实就是a、b都连续,a不连续、b连续,a、b都不连续的过程。操作他们的快速算法也就是卷积(快速傅里叶),多孔(a trous),MALLAT。在MATLAB里,也就是CWT,SWT,DWT。SWT称平稳小波变换、二进小波变换、或者非抽取小波变换。

3。从冗余性上:CWT>SWT>DWT,前面两个都冗余,后面的离散小波变换不冗余。

4。从应用上:CWT适合相似性检测、奇异性分析;SWT适合消噪,模极大值分析;DWT适合压缩

5。操作。就是在某个尺度上得到小波的离散值和原信号卷积,再改变尺度重新得到小波的离散值和原信号卷积。每一个尺度得到一个行向量存储这个尺度下的小波系数,多个尺度就是一个矩阵,这个矩阵就是我们要显示的时间-尺度图。

6。显示。“不要认为工程很简单”。我的一个老师说过的话。小波系数的显示还是有技巧的。很多人画出的图形“一片乌黑”就是个例子。第一步,一般将所有尺度下的小波系数取模;第二步,将每个尺度下的小波系数范围作映射,映射到你指定MAP的范围,比如如果是GRAY,你就映射到0-255;第三步,用IMAGE命令画图;第四步,设置时间和尺度坐标。MATLAB是个很专业的软件,它把这些做的很好,但也就使我们懒惰和糊涂,我是个好奇心重的人就研究了下。里面有个巧妙的函数把我说的(1,2)两个步骤封装在了一起,就是WCODEMAT,有兴趣的同学可以看看。

希望大家深入研究小波。

这里,还有要说的是,小波目前理论的热点:
1。不可分的小波或者具有可分性质的方向性小波;
2。XLET: CONTOURLET, WEDGELET, SHEARLET, BANDELET, RIDGELET, CURVELET; PLATELET.
3。多分辨率分析+多尺度几何分析的结合,才真正是我们所需要的。比如小波域的WEDGELET等等。

最后,几点建议:
1。理论研究和实际应用不同,工程上很多问题小波并不是最好的,在做项目的时候大家要实际情况,实际对待。
2。做研究,有些话题太成熟了,对于小波本身几乎很难再做下去,要在方向性和几何性上下功夫。对此,我也做的工作很少,毕竟研究方向转了。
3。看看DAUBECHIES,MALLAT,Minh N. Do的personal web,特别是preprint的文章,会有些启发。
4。最终大家还是最好自己提一些新的XLET或算法出来,这样才是有挑战性的工作,否则我们永远是落后的。
5。交叉领域和结合是值得做的,最近R2007也就增加了PCA+WAVELET等内容,还有用小波在数值分析领域求解偏微分等,多尺度现在是个很热的东西,SIAM为此都开了专门的期刊,这是数值分析最前沿的期刊可以看看。还有ELSEVIER的Applied and Computational Harmonic Analysis,这些虽然数学多了一点,但毕竟才有可能在工程上是“新”的。

大家可以到我的网站下载到.m文件
http://cem.ahu.edu.cn/Freecode/freecode.htm

抱歉!评论已关闭.