Shannon-Gabor小波的快速计算

(9)
令 ,则Biblioteka (10)当 时，2
采用插值小波变换理论计算小波系数、构造插值算子，在提高精度的同时，也带来了计算量的问题。下面对二维偏微分方程的插值小波算法中的参数进行优化，以实现算法的快速计算。
二维偏微分方程的求解首先对所求解区域进行空间离散，形成关于时间的常微分方程组，然后再对常微分方程组的数值进行求解，对空间离散的方法主要包括小波有限元法和小波配置法[99-102]。小波配置法要求基函数具有插值特性，在相同计算精度下，计算速度比有限元方法快；且具有自适应性，在方程解的奇异点位置自动加密配点，在解平滑处配置点稀疏，可同时满足计算效率、精度。1996年S.Bertoluzza[99]选择Daubechies小波的自相关函数为基函数构造了一种单层小波配置法，该方法虽然对边界条件具有较强的适应性，但在整个求解区间上配置点是均匀选取的，计算量较大。OlegV.Vasilyev[100]构造的多层小波配置法则避免了上述算法的缺陷，但多层小波配置法中插值基函数的构造复杂，计算量大。文献[41]利用插值小波理论构造了一种求解偏微分方程的自适应quasi-Shannon小波配置法，该方法充分利用了quasi-Shannon小波的插值特性，相对于文献[99]的方法，多层插值函数的构造简单且计算量小，但只适用于一维的偏微分方程求解。下面通过构造二维多层自适应插值小波算子、插值小波基函数的选取、插值小波快速计算方法等将自适应小波配置法推广至二维偏微分方程的求解。
(4-4)
其中的插值算子定义如下：
(4-5)
为了得到统一的多层插值小波算子，需要将插值小波系数 、 、 表达成J层上所有配置点的权重和，因此定义限制算子为：
(4-6)
限制算子表示了多层间对应值相等的配点。利用限制算子，将公式（4-4）中的 映射到第J层上。
(4-7)
因此由公式（4-4）得出小波系数表达式：
(4-27)
其中, 。为了计算方便，将 ， 映射到第J层上，对应J层中的编号为 ， ，则得出以下公式：
(4-28)
其中，
由公式（4-27）、（4-28）可得：
(4-29)
公式(4-29)中的
是周期函数，根据函数的周期性和奇偶性，对公式进行简化，得出对应的递推公式：
(4-30)
根据正弦函数的周期性，只需要计算 在[0,1/2]范围内的值即可，应该说 越大，这种方法效率越高。
分析可得小波系数通过 直接计算的时间复杂度为O(n4.log2n)。
公式（4-6）中的第1项
可以记作：
利用限制算子将 映射到J层上的所有点处
(4-12)
公式（4-3）中的第1项可记作：
(4-13)
将公式（4-9）代入(4-3)中的第2项，第2项可写为：
(4-14)
因此，由公式（4-13）和（4-14）可将公式（4-3）写为：
Shannon-Gabor
1
插值小波基函数应是具有插值特性的紧支撑或指数衰减函数，如插值样条小波、Daubechies小波的自相关函数等。但插值样条小波不具备正交性，而Daubechies小波没有解析表达式，Haar小波同时具有紧支撑性和正交性，但不连续，Faber-Schauder小波具有插值特性和紧支撑性，但不光滑，而Shannon小波不具备紧支撑性，因此选择具有近似紧支撑性和插值特性的Shannon-Gabor小波函数。
------------------------------快速计算方法-------------------
1.Shannon-Gabor小波函数的快速计算
在小波自适应插值算法中，经常需要大量计算公式（4-16）中的Shannon-Gabor小波函数 ，当用j1层上的k点计算j2层上的配点n时，公式如下：
（3）
将 、 代入式(3)，可得出
(4)
其中， ， ，即 ，
公式（4）可被简化为
(5)
下面分别对Shannon-Gabor小波函数求一阶、二阶导数：
①Shannon-Gabor小波函数的一阶导数
(6)
将 代入公式，可得
(7)
当 时，
②Shannon-Gabor小波函数的二阶导数
( )=
(8)
将 代入公式，可得
(4-1)
(4-2)
显见，第j层的小波函数等于第j+1层的尺度函数，第j层配置点对应的第J层的三个小波点分别为 、 、 。多层小波配置法需要同时考虑不同离散栅格大小下的插值算子，根据插值小波变换理论，函数 的逼近表达式为
(4-3)
其中，
由插值小波变换理论可写出插值小波系数 ，小波系数表示经过插值得到的配点函数值与真实值之间的逼近误差。
（4-15）
因此插值算子 可记作：
(4-16)
由于在不同时刻配置点的数量和位置均不相同，所以自适应算法中，在不同的时间段，必须重新构造插值算子,因此插值算子的构造速度对整个算法的计算效率有极其重要的影响。本文构造出的插值算子具有显式表达式且在形式上将单层插值小波算子和多层插值小波算子统一起来，同传统方法相比，具有较小的计算量。
公式(4-29)中的
(4-31)
对公式（4-31）进行递推计算，只需要筛选计算有限个必要点，其他配点可以映射过来，减少了数据的存储量，提高了配点查找速度。
3图像处理结果
Denoised image
Amount of the wavelet collocation points is23488
(a)t= 0.0005
4.1.1
考虑二维图像函数 ，二维小波函数对其进行逼近，配置点选取如图4-1所示。由图4-1可看出，不同层上的配点关系如下：
1j层每行包含( )个元素。
2在第j层下标为k的元素是第k+1个配点，所在的行、列坐标为 ，k mod ( +1)。
3j层下标为k的配点，对应到第J层配点的下标为 ，则
根据张量积的定义，可得出二维小波函数 和二维尺度函数 。
Denoised image
Amount of the wavelet collocation points is19413
(b)t= 0.001
(4-8)
其中，
将公式（4-7）代入到（4-8），用保留的已知点坐标值直接求小波系数，可得以下公式：
(4-9)
公式（4-9）中 表示第J层的离散点， 为一常量矩阵。
所以，当 时，
(4-10)
其中， 且 =0，1，2，…，2j-1
,
显然，当 时，
(4-11)
由（4-10）和（4-11）可递归分别求出 。
Shannon-Gabor尺度函数定义如下：
(1)
其中 是离散点间距， （ 是任意参数）是窗口大小参数。
为取Shannon-Gabor尺度函数作为基函数，按多尺度分析理论对函数 在定义域[a,b]内进行均匀离散，取离散点个数为2j+1,( )，则变量x的离散点定义为：
则基函数为
(2)
令
， ,则若取具有插值特性的Shannon-Gabor小波作为试函数，即