2022年2022年离散小波变换 .pdf

《2022年2022年离散小波变换 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年离散小波变换 .pdf（11页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章 离散小波变换3.1 尺度与位移的离散化方法减小小波变换系数冗余度的做法是将小波基函数atata1)(,的,a限定在一些离散点上取值。1. 尺度离散化：一种最通常的离散方法就是将尺度按幂级数进行离散化，即取mmaa0（ m为整数，10a，一般取20a） 。如果采用对数坐标，则尺度a的离散取值如图 3.1 所示。图 3.1 尺度与位移离散方法2. 位移的离散化：当120a时，tta)(,。（1）通常对进行均匀离散取值，以覆盖整个时间轴。（2）要求采样间隔满足Nyquist采样定理，即采样频率大于该尺度下频率通带的 2 倍。3. )(,ta=？当m 增加 1 时，尺度增加一倍， 对应的频带减

2、小一半（ 见图 2.2） ，可见采样频率可以降低一半，即采样间隔可以增大一倍。因此，如果尺度0m时的间隔为sT ，则在尺度为m2时，间隔可取smT2。此时)(,ta可表示为);(2212221,tTntTntnmsmmmsmm记作Znm,为简化起见，往往把 t轴用sT 归一化，这样上式就变为名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 11 页 - - - - - - - - - nttmmnm22)(2,（3.1）4. 任意函数)(tf的离散小波变换为Rnmfdttt

3、fnmWT)()(),(,（3.2）DWT 与 CWT 不同，在尺度位移相平面上，它对应一些如图3.1 所示的离散的点，因此称之为离散小波变换。将小波变换的连续相平面离散化，显然引出两个问题：（1） 离散小波变换)(),(),(,ttfnmWTnmf是否完全表征函数)(tf的全部信息，或者说，能否从函数的离散小波变换系数重建原函数)(tf。（ 2）是 否任 意函 数)(tf都可以 表示 为以)(,tnm为基本 单元的 加权和ZnmnmnmtCtf,)()(？如果可以，系数nmC,如何求？上述两个问题可以归结为一个。假设条件（1）满足，可合理的选择，并对,a进行适当的离散（即适当的选择sTa ,

4、0） ，那么一定存在与小波序列nm,对应的nm,序列，使得问题（ 1）的重建简单地表示为Znmnmnmftf,)(（3.3）nm,称为nm,的对偶，它可以由一个基本小波)(t通过位移和伸缩取得：nttmmnm22)(2,由上式，若存在)()(2RLtg，则有nmnmnmgfgffg,=nmnmnmgf,),(=nmnmnmfg,=nmnmnmfg,也即nmnmnmgg,故问题（ 2）也成立，其中nmnmgC,名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 11 页 - -

5、 - - - - - - - 由于问题（ 1）和问题（ 2）是统一的，我们首先来看问题（1） ，该问题的数学语言描述如下：若小波系数nmf,表征)(tf的全部信息，则应有当21ff时，;,2,1nmnmffZnm,或当0f时，nmf,=0；Znm,当1f 和2f 很接近时，Znmnmf,1,和Znmnmf,2,也必然很接近。用范数的概念来描述，即当21ff为一个很小的数时，2,2,1,nmnmnmff也必然为一个很小的数，用数学公式来描述：2212,2,1,ffBffnmnmnm，RB也即22,fBfnmnm（3.4a）若要小波系数nmf,稳定的重建f，则必须有：当序列Znmnmf,1,和Zn

6、mnmf,2,很接近时，函数1f 和2f 也很接近，即,2,2nmnmffARA（3.4b）把（3.4a）和（ 3.4b）合到一起。我们便得到一个合理的离散小波变换，该小波变换对所有)()(2RLtf必须满足下述条件：;,22,2fBffAnmnmRBA,（3.4c）满足式（ 3.4c）的离散函数序列Znmnm,;,在数学上称为“框架” 。3.2 小波框架与离散小波变换的逆变换3.2.1 小波框架（1）小波框架的定义当由基本小波)(t经伸缩和位移引出的函数族名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - -

7、 - - - - 第 3 页，共 11 页 - - - - - - - - - sjjkjkTtaat020,)(；Zkj,（3.5）具有下述性质时：;,22,2fBffAjkkjBA0（3.6）便称Zkjkjt,)(构成了一个小波框架， 称上式为小波框架条件， 其频域表示为Zjj,)2(20（3.7）（2）小波框架的性质1）满足小波框架条件的)(,tkj，其基本小波)(t必定满足容许性条件。但是并不是满足容许性条件的小波，在任意离散间隔sT 及尺度基数0a 下都满足小波框架的条件。2）小波函数的对偶函数kttjjkj22)(2,也构成一个框架，其框架的上、下界是)(,tkj框架上、下界的倒数

8、：22,21,1fBffAjkkj（3.8）3）离散小波变换具有非伸缩和时移共变性。4）离散小波变换仍然具有冗余度。3.2.2 离散小波变换的逆变换与重建核问题1. 离散小波变换的逆变换如离散小波序列Zkjkjt,)(，构成一个框架，其上、下界分别为A和B，则当BA时（紧框架），由框架概念可知离散小波变换的逆变换为)(),(1)()(,)(,1tkjWTAttfAtfkjkjfkjjkj（3.9）当BA，而A，B比较接近时，作为一阶逼近，可取)(2)(,tBAtkjkj（3.10）则重建公式近似为)(),(2)()(,)(,tkjWTBAttftfkjkjfkjjkj（3.11）逼近误差的范数

9、为fBABAfRRf名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 11 页 - - - - - - - - - 由上式可见，A与B愈接近，逼近误差就愈小。为了保证kj ,能构成一个重建误差较小的框架就必须对基本小波在,a轴上的采样间隔提出更高要求：0a 不一定等于 2，sT 也不一定等于 1，以便于使A和B接近于相等，可以想像，当尺度间隔愈密，位移间隔愈小。离散栅格愈接近于覆盖整个 a半平面，AB /就愈接近于 1. 关于BA、与、0a，以及）（间的关系的部分结论如下：

10、如Znmnm,是一个框架，则框架的上界A、下界B满足下面的不等式：BdaA20)(log（3.12）特别对紧框架有：daA20)(log（3.13）举例：将 Marr 小波离散化为小波框架。Marr 小波是常用的一种连续小波形式。若将Marr 小波的尺度及位移分别离散化为ktaatjjkj020,)(则可证明，)(,tkj构成了一个)(2RL空间的小波框架， 其框架的上界A、下界B同、0a之间的关系如表3.1 表示。表 3.1 Marr 小波框架上、下界同0a和之间的关系0aABAB2 0.25 13.091 14.183 1.083 2 0.50 6.546 7.092 1.083 2 0.

11、75 4.364 4.728 1.083 2 1.00 3.223 3.596 1.161 2 1.25 2.001 3.454 1.726 2 1.50 0.325 4.221 12.984 20.25 27.273 27.278 1.0002 20.50 13.673 13.639 1.0002 21.00 6.768 6.870 1.015 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 11 页 - - - - - - - - - 21.50 2.609 6.48

12、3 2.485 3120.50 20.457 20.457 1.0000 3121.00 10.178 10.279 1.010 3121.50 4.629 9.009 1.947 4120.50 27.276 27.276 1.0000 4121.00 13.586 13.690 1.007 4121.50 6.594 11.590 1.758 由表 3.1 可知：1）当20a时，取;75.020a时，取; 1302a时，取1或402a时，取; 1均可使BA，可近似为紧框架。此时采用重建公式（ 3.9）可较精确地重构原函数。2）0a 一定时，AB /的值随增大而增大。3）给定一个0a 值，只

13、要足够小，总可以得到一个近似紧的小波框架。4）20a，1时，BA，不是紧框架。2. 重建核公式（1）正交性：只有当1BA时，框架)(,tkj变为正交基，此时经框架变换后的信息无任何冗余。但在其他情况下，框架)(,tkj并不正交，具有一定的相关性。因此经框架处理后所含的信息是有冗余的。（2）紧框架情况下的小波变换系数的相关性：将离散小波变换的逆变换公式（3.9）重写如下：)(),(1)(,tkjWTAtfkjkjf（3.14）其中RkjfdtttfkjWT)()(),(,（3.15a）则名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - -

14、- 名师精心整理 - - - - - - - 第 6 页，共 11 页 - - - - - - - - - RkjfdtttfkjWT)()(),(00,00（3.15b）将式（ 3.14）代入式（ 3.15b）得RkjjkjkffdtttkjWTAkjWT)( )(),(1),(00,00jkfkjjRkjkfkjWTkjkjKAdtttkjWTA),(),;,(1)()(),(100,00（3.16）其中)(),()()(),;,(0000,00ttdtttkjkjKkjkjkjRkj（3.17）分析说明：(1)与连续情况一样，式（ 3.16）给出任意一点),(00kj处小波变换之值与栅格

15、上其他各点小波变换系数之间的内在联系，称它为重建核方程，称K 为重建核，由小波框架本身决定。(2)并不是相平面上的任意离散函数),(kjF都可看作是某一函数的离散小波变换，只有它们之间满足（3.16）时才可以被看作为某一函数的离散小波变换序列。(3)无论将 Marr 小波如何离散，都不能使1BA，也即它不可能构成)(2RL空间的正交基。（Morlet 小波和 DOG 也是如此）3.3 二进小波变换对于尺度及位移均离散化的小波序列，若取离散栅格的20a，0，即相当于连续小波只在尺度上进行了二进制离散，而位移仍取连续变化， 我们称这类小波为二进小波，表示为kktatk2)(20,2（3.18）二进

16、小波介于连续小波和离散小波之间，它只是对尺度参量进行了离散化，而在时间域上的平移量仍保持连续变化，故二进小波具有时移共变性， 在奇异性检测、图像处理方面十分有用。 在讨论二进小波变换及逆变换公式时，我们仍借用离散小波框架理论对其进行分析。3.3.1 二进小波变换及其逆变换设小波函数为)(t，其傅里叶变换为)(，若存在二常数BA0，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 11 页 - - - - - - - - - 使得BAZkk2)2(（3.19）此时式（3.18

17、）定义的二进小波才是有意义的二进小波，即其逆变换存在。称式（3.19）为二进小波的稳定性条件；若BA，则称最稳定条件。若定义函数)(2RLf的二进小波变换系数为RkkdtttfttfWTkk2)(2)()()(2,22（3.20）由卷积定理，设)(2kWT的傅里叶变换为)(2kWT，则)2(2)()(22kjkeFWTk因此，稳定性条件（ 3.19）等价于对任意)(2RLf都有2222)(fBWTfAZkj（3.21）式（3.21）说明：1）二进小波)(,2tk构成了)(2RL的一个框架。2）二进小波)(,2tk的小波变换公式（式（ 3.20） ）及其逆变换公式存在。二进小波变换的重建公式为Z

18、kRdtWTtfkk)()()(,22（3.22）其中，)(,2tk为)(,2tk的对偶框架，其上、下界分别为11, AB。同离散小波框架相似，当BA时，)(1)(,2,2tAtkk（3.23）当BA时，)(,2tk的一阶近似为)(2,2tBAk（3.24）当AB接近于 1 时，其重构误差减小。当1AB采用高阶近似或递推的方法就可求得更精确的解。3.3.2 二进小波变换（1）与离散小波相同，二进小波也一定是一个允许小波，且有2ln)(2ln02BdA名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - -

19、 - - 第 8 页，共 11 页 - - - - - - - - - 2ln)(2ln02BdA特别是，当BA时，2ln)(02AdC（3.25）（2）二进小波变换时冗余的由框架理论可知， 当不满足1BA是，框架是冗余的， 也即二进变换系数之间具有一定的相关性， 它们之间的关系满足重建核方程。紧框架情况下的重建核方程如下：紧框架（BA）时，由（ 3.22）和（ 3.23）可知，重建公式为ZkRdtWTAtfkk)()(1)(,22（3.26）由于RkkdtttfttfWTkk2)(2)()()(2,22当尺度为 m，平移为0时，小波变换系数为RmmdtttfWTm2)(2)(022=RmZk

20、RmdttdtWTAkk2)()(210,222=RRmkZkRkmddtttWTAk22)(2210222=RZkRdkmKWTAk),;,()(102其中dtttKRkmkm22221022（3.27）此即为二进小波变换紧框架下的重建核方程。说明：由重建核方程可知，并不是任意函数序列Zkkg)(2都可以作为某一函数的二进小波变换， 而只有当它们满足重建核方程时，才可以看作是某一函数的二进小波变换。（3）二进小波变换具有平移不变性（时域平移不变性 ） ，即若)()(00tftf名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名

21、师精心整理 - - - - - - - 第 9 页，共 11 页 - - - - - - - - - 设)(tf的二进小波变换为)(2kWT，)(0tf的二进小波变换为)(2kWT，则有)()(022kkWTWT（3.28）证明略（习题）表 3.1 Marr 小波框架上、下界同0a和之间的关系0aABAB2 0.25 13.091 14.183 1.083 2 0.50 6.546 7.092 1.083 2 0.75 4.364 4.728 1.083 2 1.00 3.223 3.596 1.161 2 1.25 2.001 3.454 1.726 2 1.50 0.325 4.221 1

22、2.984 20.25 27.273 27.278 1.0002 20.50 13.673 13.639 1.0002 21.00 6.768 6.870 1.015 21.50 2.609 6.483 2.485 3120.50 20.457 20.457 1.0000 3121.00 10.178 10.279 1.010 3121.50 4.629 9.009 1.947 4120.50 27.276 27.276 1.0000 4121.00 13.586 13.690 1.007 4121.50 6.594 11.590 1.758 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 11 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 11 页 - - - - - - - - -