MWMP简介.pptVIP

{PAGE} 下面是增加d8小波与balanced8小波，这些小波已经含在我所修改的coef.m文件中。并请注意d8是单小波。 下面是增加bio小波与ort4小波，这些小波已经含在我所修改的coef.m文件中。注意bio小波有些问题。 这是我增加的ort5,ort6平衡多小波。 这是我增加的jort6多小波。 平衡多小波不需要预滤波与后处理过程。对于其它多小波，进行小波变换前，都需要进行预滤波，以便将传统数据化成向量形式，能够进行小波变换。对小波变换后的数据进行反变换后，还要进行后处理，使向量形式的数据再变成传统数据。预小波有两种最常用的方法：重复预滤波与插值预滤波。 多小波预滤波文件的修改方法： 如果你增加的不是平衡多小波，则必须在文件coef_prep.m中针对此多小波，给出其预滤波方法PR与后处理方法PO，一般来说，它们都是可逆矩阵并且这两个矩阵的乘积为单位阵。 对每个非平衡多小波，其预小波的方法很多。重复预滤波对每个非平衡多小波都适用，但是，增加了计算量。插值预滤波效果较好，但是不保持正交性。还有其它各种各样的预滤波方法。 对于多小波来说，最大的计算难点是，每个多小波都有特定的预滤波与后处理方法。 下面是我编写的几个多小波的预滤波与后处理在文件coef_prep.m的实现。 在程序coef_prep.m的后面加入以下内容： MWMP中多小波的预滤波器的名称： ghmap对ghm小波，进行双正交插值预滤波。 ‘ghmorap’ ‘ghm’小波的正交预滤波器。 ‘clap’ ‘cl’小波的双正交插值预滤波器。 bih5ap bih52s , bih54n双正交插值预滤波器。 bih3apbih32s的双正交插值预滤波器。 ‘sa4ap’ ‘sa4’小波的正交插值预滤波器。 bighm2apbighm2的双正交插值预滤波器。 id恒等滤波器，即直接将信号按顺序化为向量形式。 ‘cl3ap’, ‘bioap’, ‘jort6ap’此是我增加的多小波的预滤波器。 以上的预滤波器中，为什么没有重复预滤波器，这是因为，重复预滤波用一个特定的函数来实现，因为对于各种多小波，其重复预滤波方法是一样的。 对给定的一维信号，prep1D_rr进行重复信号预滤波，即将每个单独的信号变为重复行向量的形式，此函数调用形式： rXn维预滤波后的数组=prep1D_rr(1Xn的一维信号数组，滤波器名称) 注意，此处的n一定要为integer*2^maxlevel的形式。例子： fp=prep1D_rr([1 2 3 4 5 6 7 8],ghm) fp = 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 7.0000 8.0000 0.7071 1.4142 2.1213 2.8284 3.5355 4.2426 4.9497 5.6569 [1 2 3 4 5 6 7 8]/sqrt(2) ans = 0.7071 1.4142 2.1213 2.8284 3.5355 4.2426 4.9497 5.6569 postp1D_rr是一维小波重复预滤波后滤波（后处理）函数。如果对给定的一维信号，使用prep1D_rr进行重复信号预滤波，然后使用dec1D_pe进行一维小波变换，将变换后的结果使用rec1D_pe进行小波逆变换，则逆变换后，是一个rXn的数组，设它为fhatp，它正好是一维信号经过预滤波后的形式，则可使用fhat=postp1D_rr(fhatp,flt)重构原来的一维信号，其中flt是coef.m中的滤波器。postp1D_rr和prep1D_rr是一对互逆的预滤波和后处理函数，注意，若数据完全重构，则小波正变换和逆变换的滤波器应该相同或者是对偶小波。例子： fp=prep1D_rr([1 2 3 4 5 6 7 8],ghm);transf=dec1D_pe(fp,ghm,2); fhatp=rec1D_pe(transf,ghm,2);fhat=postp1D_rr(fhatp,ghm) fhat = 1.0000 2.0000 3.0000 4.0000 5.0000 6.0000 7.0000 8.0000 prep2D_rr是二维多小波变换的预滤波函数。对于给定的二维信号，prep2D_rr使用重复预滤波。调用形式为fp=prep2D_rr(f,flt),其中f是一个nXn的原始信号矩阵,矩阵中的每个元素都是实数,n具有r*integer*2