图12.4.1由得到尺度函数.docVIP

图12.4.1 由得到尺度函数 很容易证明(12.4.4)式所决定的的傅里叶变换是 (12.4.8) 而对移位后的，其傅里叶变换为 (12.4.9) 如果为偶数，式中，若为奇数，则式中。 分析(12.4.6)式，我们发现 (12.4.10) 满足我们在第十章所讨论的二尺度差分方程。同时，可求出 (12.4.11) 是有界的。当时， (12.4.12) 同样也满足二尺度差分方程，同理可求出 (12.4.13) 因此，在时不同的可构成一个多分辨率分析。由1.7节关于正交基频域的性质，由于(12.4.11) 和(12.4.13)式的右边不等于1，因此的整数移位之间不构成正交基。由(9.8.40)式，我们可将“正交化”，即令 (12.4.14)对作反变换，得尺度函数，则，可形成一族正交基。再由第十章的方法可得到正交归一的小波函数。 在双正交的情况下，我们可不必对作(12.4.14)式的正交化，而直接用作适当移位后的作为尺度函数，如(12.4.5)～(12.4.7)式所示。这样选定后，Daubechies令(12.3.11)式中的等于零，并令。，从而得到了在双正交条件下样条小波分析滤波器和重建滤波器的系数，即 ， (12.4.15a) ， (12.4.15b) ， (12.4.16a) ， (12.4.16b) (12.4.15a)和 (12.4.15b)分别对应(12.3.5b)式和 (12.3.6b)式，而(12.4.16)式是(12.3.5a)式、(12.3.6a)式和(12.3.11)式的结合。 由(12.4.15)式可以看出，仅和有关，而和无关；由(12.4.16)式，不但和有关，而且还和有关，也即取决于和。给定不同的和，就可求出一对和。将(12.4.15)式和(12.4.8)及(12.4.9)式相比较可以看出，尺度函数的傅里叶变换的“阶次”和中的等价，也即即是得到时由卷积的次数，或称之为的“阶次”。 现给出不同和组合情况下、、、、和的系数。 情况1. 令，则必有，由(12.4.15b)式，有 所以， 即 令，则必有，由(12.4.16b)式，有 [注]：Daubechies在文献[5]中令，故和本书定义的有区别。 所以 即 在时的尺度函数即是Haar尺度函数，即，对，其余为零。又由于在时的，由(12.2.12b)式，必有。 易知在该情况下的小波函数即是Haar小波，即 我们知道， Haar小波属正交小波，即DB1，但因为它是对称的，故又属双正交小波。我们记在该情况下的，，和分别为，，及。前面的1代表，后面的1代表，以下均相同。 在的情况下，我们再令，则必有，由(12.4.16b)式，有 所以 因为仍为1，所以， 不变，，及可由上一节的公式推出。但MATLAB中的wavefun.m文件可用来产生这些函数，如图12.4.2所示。 情况2：令，则，由(12.4.5a)式，有 再令，则，由(12.4.16a)式，有 即 按此方法类推，读者不难得出在不同和的组合下的及，如表12.4.1所示。和的合适组合是： ， ， ， 由该表可以看出，和，特别是，若不考虑前面的，其分母的系数都是2的整次幂，因此有利于在计算机上快速实现。此外，在不同的值下，都是精确已知的，这些都是样条点双正交小波的优点。在不同组合下的，，及如图12.4.2~12.4.6所示。其中图12.4.2给出的是，和5时的及。图中左边标的“bior 1.3 phi-D”，即是指双正交小波的尺度函数，对应，，“D”代表分解，右边图标的“psi”指的是“”。以下各图的标法均相同。另外，图中的横坐标是MATLAB按正的坐标求出的，这和12.3节所给出的支撑范围有区别。 图12.4.2 ，时用于分解的和 图12.4.5给出的是，，和时的(图中标为phi-R，R代表重建)。显然，时的即是Haar尺度函数，它即是图12.4.1的(a)图。而时的即是图12.4.1(b)，时的即是图12.4.1(c)。它们分别是，和。显然，图中时的应是。注意，只和有关，而和无关。 图12.4.6给出的是， ；，；，和，时的。显然，不仅和有关，而且也和有关。 图12.4.3，，时用于分解的和 图12.4.4，，时用于分解的和 12.4.5，时用于重建的 图12.4.6，， ；，；，和，时用于重建的 表12.4.1给出了在时取不同值时，的系数，由于在，和8，，和9时的的系数过长，故表中没有列入，详细数据