Pure Mathematics
Vol. 10  No. 11 ( 2020 ), Article ID: 38902 , 9 pages
10.12677/PM.2020.1011132

基于多分辨分析的小波紧框架构造

李晨,李万社

陕西师范大学,数学与信息科学学院,陕西 西安

收稿日期:2020年11月1日;录用日期:2020年11月20日;发布日期:2020年11月27日

摘要

本文主要研究了离散序列基于多分辨分析(MRA)生成小波紧框架的构造方法。首先介绍了多分辨分析和基于多分辨分析的小波紧框架的相关知识,然后给出所定义的离散序列生成平方可积函数空间小波紧框架的充分条件,最后利用Matlab工具给出两个数值算例构造的细分函数和小波函数的函数图像。

关键词

小波紧框架,多分辨分析,细分函数

Wavelet Tight Frame Construction Based on Multi-Resolution Analysis

Chen Li, Wanshe Li

School of Mathematics and Information Science, Shaanxi Normal University, Xi’an Shaanxi

Received: Nov. 1st, 2020; accepted: Nov. 20th, 2020; published: Nov. 27th, 2020

ABSTRACT

This paper mainly studies the construction method of tight wavelet frame generated by discrete sequences based on multi-resolution analysis (MRA). First, we introduce the related knowledge of multi-resolution analysis and wavelet tight frames based on multi-resolution analysis. And then sufficient conditions for discrete sequences are given to generate wavelet tight frame with square integrable function space. Finally, the function images of refinable function and wavelet function constructed by two numerical examples are given by Matlab.

Keywords:Wavelet Tight Frame, Multi-Resolution Analysis, Refinable Function

Copyright © 2020 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

1. 引言

经过多年的研究和发展,框架理论已经形成了一套相对比较完整的理论体系,其中小波框架和Gabor框架是框架理论中的一个十分重要研究领域 [1] [2]。小波框架和Gabor框架的大量文献目前在信号图像处理、时频分析、量子力学、电子工程的滤波器设计等多个相关方面被广泛应用 [3] [4]。由于基于多分辨分析(MRA)的小波紧框架,为分段函数提供了稀疏逼近且保证了快速分解算法和重构算法的存在,故基于多分辨分析的小波紧框架在图像处理中得以许多应用 [5] [6] [7]。基于的多分辨分析小波紧框架可以看作为基于正交多分辨分析小波的推广,它与基于正交多分辨分析小波是小波分析研究中的两大热点,其区别主要在于冗余性质。正是由于基于小波多分辨分析紧框架具有的冗余性使得应用基于小波多分辨分析紧框架在恢复原始信号的数值计算中更加稳定,与原始信号误差更小 [8] [9] [10] [11]。

2. MRA和基于MRA小波紧框

在本节中,主要介绍多分辨分析(MRA)以及基于多分辨分析小波紧框架的相关知识。在本文中,我们先给出平方可积函数空间 L 2 ( R ) 上的相关定义:

f , g = R f ( x ) g ( x ) ¯ d x f , g L 2 ( R )

f 2 = f , f f L 2 ( R )

f ^ ( ξ ) = R f ( x ) e i ξ x d x f L 2 ( R )

a ^ ( ξ ) = k Ζ a ( k ) e i k ξ a l 2 ( Ζ )

1987年,Meyer和Mallat第一次提出了多分辨分析的概念,将多分辨分析思想应用到小波分析中,给出了构造 L 2 ( R ) 正交小波基的一般方法,并提出了快速小波变换分解和重构算法,即Mallat算法。1988年,Daubechies构造了基于多分辨分析的紧支撑正交小波基。在本文中,我们应用C. De Boor等人给出的更为一般的多分辨分析 [12]。

对于 φ L 2 ( R ) ,记 φ j , k = p j / 2 φ ( p j x k ) j , k Ζ ,此处 p Ζ + p > 1 。记 V j = Span { φ j , k : k Ζ } ¯ 为空间 L 2 ( R ) 中的由函数族 { φ j , k : k Ζ } 经过平移伸缩所张成的空间 L 2 ( R ) 中的线性闭子空间。下面给出多分辨分析的概念。

定义2.1:闭子空间序列 { V j } j Z 称为一个多分辨分析(MRA),如果 { V j } j Z 满足下列条件:

1) V j V j + 1 , j Ζ

2) j Ζ V j = { 0 }

3) j Ζ V j ¯ = L 2 ( R )

则称 { V j } j Z 为函数 φ 生成空间 L 2 ( R ) 的一个多分辨分析。

定义2.2:若函数 φ L 2 ( R ) 存在离散序列 a 0 l 2 ( Ζ ) 满足

φ ( x ) = p k Ζ a 0 ( k ) φ ( p x k ) , x R (1)

则称函数 φ ( x ) 为细分函数(或尺度函数),称序列 a 0 为细分函数 φ 的细分面具。对等式(1)两边同时做傅里叶变换,得到

φ ^ ( p ξ ) = a ^ 0 ( ξ ) φ ^ ( ξ ) , ξ R (2)

类似于多分辨分析中细分函数 φ 的定义,定义小波函数族 Ψ = { ψ l } l = 1 M 1

ψ l ( x ) = p k Ζ a l ( k ) φ ( p x k ) , x R (3)

其中称离散序列族 { a l } l = 1 M 1 l 2 ( Ζ ) 为函数族 { ψ l } l = 1 M 1 的小波面具,同理对等式(3)两边同时做傅里叶变换,得到

ψ ^ l ( p ξ ) = a ^ l ( ξ ) φ ^ ( ξ ) , x R (4)

从而可定义小波系 Χ ( Ψ ) ,即为

Χ ( Ψ ) = { ψ l , j , k = p j / 2 ψ l ( p j x k ) : 1 l M 1 , j , k Ζ } (5)

定理2.1:设函数 φ L 2 ( R ) 中的紧支撑细分函数满足 φ ^ ( 0 ) = 1 ,并定义由函数 φ 生成的闭子空间 V j = Span { φ j , k : k Ζ } ¯ ,则可得闭子空间序列 { V j } j Z 形成一个多分辨分析。

定义2.3:希尔伯特空间H上的函数集合 { f k : k J } 称为一个框架,如果对 f H ,均存在正常数A,B使得下式成立

A f 2 k J | f , f k | 2 B f 2

并称A,B分别为框架的下界和上界。

特别的,若有 A = B = 1 ,则称框架 { f k : k J } 为紧框架,即 f 2 = k J | f , f k | 2

定理2.2:(酉扩展定理)设 φ L 2 ( R ) 为一个细分函数且满足 φ ^ ( 0 ) = 1 ,其细分面具 a 0 l 2 ( Ζ ) 为有限支撑序列,根据(3)式给出由小波面具 { a l } l = 1 M 1 l 2 ( Ζ ) 生成的函数族 { ψ l } l = 1 M 1 ,如果面具 { a 0 , a 1 , , a M 1 } 满足

l = 0 M 1 a ^ l ( ξ ) a ^ l ( ξ + 2 π μ ) ¯ = δ μ , 0 a . e . ξ R , μ Ζ / p Ζ (6)

等价于 m Ζ Ω j = ( p Ζ + j ) supp ( a 0 ) j Ζ / p Ζ

l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = 1 p δ m , 0 (7)

则根据(5)式中所定义的小波系 Χ ( Ψ ) 形成 L 2 ( R ) 的一个小波紧框架,且由函数 φ 生成空间 L 2 ( R ) 中的多分辨分析,即闭子空间 { V j } j Z

3. p-带小波紧框架的构造

在本节中,我们主要根据酉扩展定理给出离散有限支撑序列 { a l } l = 1 M 1 l 2 ( Ζ ) 通过(5)式所定义的小波系统形成 L 2 ( R ) 小波紧框架的充分条件。

定理3.1:设 a l 2 ( Ζ ) 为一个有限支撑序列,由离散序列a来定义离散序列 a l ( m ) = a ( m ) e i 2 π l m / M l = 0 , 1 , , M 1 ,然后根据(1)式和(3)式分别定义细分函数 φ 和函数族 Ψ = { ψ l } l = 1 M 1 ,则若满足:

i) M N

ii) k Ω j | a ( k ) | 2 = 1 p M Ω j = ( p Ζ + j ) [ 0 , N 1 ] j Ζ / p Ζ

iii) k Ω j a ( k ) = 1 p

则由(1)式所定义的细分函数 φ 生成空间 L 2 ( R ) 的一个多分辨分析,且根据(3)式和(5)式定义的小波系 Χ ( Ψ ) 形成 L 2 ( R ) 的一个小波紧框架。

为了证明定理3.1给出以下的引理。

G ( ξ ) = [ a ^ 0 ( ξ ) a ^ 1 ( ξ ) a ^ M 1 ( ξ ) a ^ 0 ( ξ + 2 π p ) a ^ 1 ( ξ + 2 π p ) a ^ M 1 ( ξ + 2 π p ) a ^ 0 ( ξ + 2 π ( p 1 ) p ) a ^ 1 ( ξ + 2 π ( p 1 ) p ) a ^ M 1 ( ξ + 2 π ( p 1 ) p ) ]

引理3.1:设面具 { a 0 , a 1 , , a M 1 } 满足(6)式,则对 ξ R

k = 0 p 1 | a ^ l ( ξ + 2 π k p ) | 2 1 , l = 0 , 1 , , M 1 (8)

证明:面具 { a 0 , a 1 , , a M 1 } 满足(6)式,即 G ( ξ ) G ( ξ ) = Ι p ,其中 G ( ξ ) 表示矩阵 G ( ξ ) 的共轭转置。根据矩阵的相关性质,我们只需证明对 l = 0 时成立即可。

G 0 ( ξ ) = [ a ^ 1 ( ξ ) a ^ 2 ( ξ ) a ^ M 1 ( ξ ) a ^ 1 ( ξ + 2 π p ) a ^ 2 ( ξ + 2 π p ) a ^ M 1 ( ξ + 2 π p ) a ^ 1 ( ξ + 2 π ( p 1 ) p ) a ^ 2 ( ξ + 2 π ( p 1 ) p ) a ^ M 1 ( ξ + 2 π ( p 1 ) p ) ]

且记 α = ( a ^ 0 ( ξ ) , a ^ 0 ( ξ + 2 π p ) , , a ^ 0 ( ξ + 2 π ( p 1 ) p ) ) T ,则可得

G 0 ( ξ ) G 0 ( ξ ) = Ι p α α T = ( β 1 , β 2 , , β p )

其中 { β 1 = ( 1 | a ^ 0 ( ξ ) | 2 , a ^ 0 ( ξ + 2 π p ) a ^ 0 ( ξ ) ¯ , , a ^ 0 ( ξ + 2 π ( p 1 ) p ) a ^ 0 ( ξ ) ¯ ) T β 2 = ( a ^ 0 ( ξ ) a ^ 0 ( ξ + 2 π p ) ¯ , 1 | a ^ 0 ( ξ + 2 π p ) | 2 , , a ^ 0 ( ξ + 2 π ( p 1 ) p ) a ^ 0 ( ξ + 2 π p ) ¯ ) T β p = ( a ^ 0 ( ξ ) a ^ 0 ( ξ + 2 π ( p 1 ) p ) ¯ , a ^ 0 ( ξ + 2 π p ) a ^ 0 ( ξ + 2 π ( p 1 ) p ) ¯ , , 1 | a ^ 0 ( ξ + 2 π ( p 1 ) p ) | 2 ) T

计算可知,矩阵 G ( ξ ) G ( ξ ) 的p个特征值分别为 λ 1 = λ 2 = = λ p 1 = 1 , λ p = 1 k = 0 p 1 | a ^ l ( ξ + 2 π k p ) | 2

又因为矩阵 G ( ξ ) G ( ξ ) 为一个半正定矩阵,从而有 λ p 0 ,即 l = 0 时(8)式成立。

引理3.2:设 a l 2 ( Ζ ) 为一个有限支撑序列,即 a ( m ) = 0 m [ 0 , N 1 ] Ζ ,且可由序列 a l 2 ( Ζ ) 定义M个离散紧支撑序列,即 a l ( m ) = a ( m ) e i 2 π l m / M l = 0 , 1 , , M 1 且如果满足:

i) M N

ii) k Ω j | a ( k ) | 2 = 1 p M Ω j = ( p Ζ + j ) [ 0 , N 1 ] j Ζ / p Ζ

则有定义的离散序列 a l ( m ) l = 0 , 1 , , M 1 满足 l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = 1 p δ m , 0

证明: m [ 0 , N 1 ] Ζ ,由 a ( m ) = 0 可得

m [ 0 , N 1 ] Ζ a l ( m ) = a ( m ) e i 2 π l m / M = 0 l = 0 , 1 , , M 1 。则有 m [ 1 N , N 1 ]

l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = 0

1 m N 1 (或 1 N m 1 ),由(i) M N 可知

l = 0 M 1 e i 2 π l m / M = 0

则当 1 m N 1 (或 1 N m 1 )时,

l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = k Ω j l = 0 M 1 a l ( k ) a l ( k + m ) ¯ = k Ω j a l ( k ) a l ( k + m ) ¯ l = 0 M 1 e i 2 π l m / M = 0

下面考虑当 m = 0 时,

l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = k Ω j l = 0 M 1 a l ( k ) a l ( k + m ) ¯ = k Ω j a l ( k ) a l ( k + m ) ¯ l = 0 M 1 e i 2 π l m / M = M k Ω j | a ( k ) | 2

由(ii) k Ω j | a ( k ) | 2 = 1 p M 可得 l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = 1 p

即证 l = 0 M 1 k Ω j a l ( k ) a l ( k + m ) ¯ = 1 p δ m , 0 , l = 0 , 1 , , M 1

引理3.3:设序列 a l 2 ( Ζ ) 为一个有限支撑序列,对于 Ω j = ( p Ζ + j ) [ 0 , N 1 ] j Ζ / p Ζ ,满足 k Ω j a ( k ) = 1 p ,定义离散序列 a l ( m ) = a ( m ) e i 2 π l m / M l = 0 , 1 , , M 1 ,且满足(6)式,则根据定义的细分函数 φ L 2 ( R ) 且满足 φ ^ ( 0 ) = 1

证明:根据定义有 k Ω j a 0 ( k ) = k Ω j a ( k ) = 1 p ,则可知 k Ζ a 0 ( k ) = 1 ,又由 a 0 有限支撑,从而得 a ^ 0 ( 0 ) = 1 ,即存在唯一紧支撑细分函数 φ 满足 φ ^ ( 0 ) = 1 ,且其傅里叶变换为 φ ^ ( ξ ) = j = 0 a ^ 0 ( p j ξ )

下面定义函数列 { f ^ n } n N

f ^ n ( ξ ) = a ^ 0 ( p 1 ξ ) f ^ n 1 ( p 1 ξ ) = j = 1 n a ^ 0 ( p j ξ ) f ^ 0 ( p j ξ )

其中 f ^ 0 = 1 [ π , π ] ,且 { f ^ n } n N 的依点列收敛到 φ ^

f ^ n 2 = π p n π p n | j = 1 n a ^ 0 ( p j ξ ) | 2 d ξ = 0 2 π p n | j = 1 n a ^ 0 ( p j ξ ) | 2 d ξ = μ = 0 p 1 2 π μ p n 1 2 π ( μ + 1 ) p n | j = 1 n a ^ 0 ( p j ξ ) | 2 d ξ = μ = 0 p 1 0 2 π p n 1 | j = 1 n a ^ 0 ( p j ξ + p n j 1 2 π μ ) | 2 d ξ = 0 2 π p n 1 | j = 1 n a ^ 0 ( p j ξ ) | 2 ( μ = 0 p 1 | a ^ 0 ( p n ξ + 2 π μ p ) | 2 ) d ξ

根据引理3.1,当 l = 0 时,即有 k = 0 p 1 | a ^ 0 ( ξ + 2 π k p ) | 2 1

从而 f ^ n 2 0 2 π p n 1 | j = 1 n a ^ 0 ( p j ξ ) | 2 d ξ = f ^ n 1 2

根据归纳法,得

f ^ n 2 f ^ 0 2 = 2 π n 0

因为 { f ^ n } n N 依点列收敛到 φ ^ ,根据法图引理可得 φ ^ 2 lim inf n f ^ n 2 2 π ,即证 φ L 2 ( R )

利用以上三个引理和定理2.2易证定理3.1。

4. 3-带小波紧框架数值算例

例1. 当 M = 9 , p = 3 时,我们可得离散序列 a 0 = 1 9 ( 1 , 1 , , 1 9 ) ,从而可以分别给出相应的细分面具和小波面具 a l ( m ) l = 0 , 1 , , M 1 ,由 a l ( m ) 生成的细分函数和函数族 { ψ l } l = 1 M 1 的函数图像如图1

φ ψ 1 ψ 2 ψ 3 ψ 4 ψ 5 ψ 6 ψ 7 ψ 8

Figure 1. The refinable function and framelets if, M = 9 , p = 3 (the blue and red images correspond to the real and imaginary parts of the function respectively)

图1. M = 9 , p = 3 时细分函数和小波函数的函数图像(其中蓝色、红色图像分别对应函数的实部图像、虚部图像)

φ ψ 1 ψ 2 ψ 3 ψ 4 ψ 5 ψ 6 ψ 7 ψ 8 ψ 9 ψ 10 ψ 11 ψ 12 ψ 13 ψ 14 ψ 15

Figure 2. The refinable function and framelets if, M = 16 , p = 4 (the blue and red images correspond to the real and imaginary parts of the function respectively)

图2. M = 16 , p = 4 时细分函数和小波函数的函数图像(其中蓝色、红色图像分别对应函数的实部图像、虚部图像)

例2. 当 M = 16 , p = 4 时,我们可得离散序列 a = 1 16 ( 1 , 1 , , 1 16 ) ,从而可以分别给出相应的细分面具和小波面具 a l ( m ) , l = 0 , 1 , , M 1 ,由 a l ( m ) 生成的细分函数和函数族 { ψ l } l = 1 M 1 的函数图像如图2

文章引用

李 晨,李万社. 基于多分辨分析的小波紧框架构造
Wavelet Tight Frame Construction Based on Multi-Resolution Analysis[J]. 理论数学, 2020, 10(11): 1106-1114. https://doi.org/10.12677/PM.2020.1011132

参考文献

  1. 1. 黄永东, 程正兴. 带小波紧框架的显式构造方法[J]. 数学物理学报, 2007, 27(1): 7-18.

  2. 2. 黄焱. 基于MRA的三带小波紧框架的参数化研究[D]: [硕士学位论文]. 北京: 北京化工大学, 2012.

  3. 3. Mallat, S. (1989) A The-ory of Multiresolution Signal Decomposition: The Wavelet Transform. IEEE Transactions, 11, 674-693. https://doi.org/10.1109/34.192463

  4. 4. Zhang, Z. and Saito, N. (2009) Constructions of Periodic Wavelet Frames Using Extension Principles. Applied and Computation Harmonic Analysis, 27, 12-23. https://doi.org/10.1016/j.acha.2008.10.003

  5. 5. Dong, B. and Shen, Z. (2010) MRA-Based Wavelet Frames and Applications, IAS/Park City Mathematics Series: The Mathematics of Image Processing, Vol. 19, Park City Mathematics Institute, 7-185.

  6. 6. Fan, Z., Ji, H. and Shen, Z. (2016) Dual Gramian Analysis: Duality Principle and Unitary Extension Principle. Mathematics of Computation, 85, 239-270. https://doi.org/10.1090/mcom/2987

  7. 7. Fan, Z., Heinecke, A. and Shen, Z. (2016) Duality for Frames. Journal of Fourier Analysis and Applications, 22, 71-136. https://doi.org/10.1007/s00041-015-9415-0

  8. 8. Christensen, O. (2008) Frames and Bases: An Introductory Course. Birkhauser, Boston. https://doi.org/10.1007/978-0-8176-4678-3

  9. 9. Hur, Y. and Lubberts, Z. (2017) New Constructions of Nonseparable Tight Wavelet Frames. Linear Algebra and its Applications, 534, 13-35. https://doi.org/10.1016/j.laa.2017.08.002

  10. 10. 方清城. MATLAB R2016a 小波分析22个算法实现[M]. 北京: 电子工业出版社, 2018: 338-342.

  11. 11. 郝菁, 崔丽鸿. 构造 带对偶小波紧框架一个充分条件的证明[J]. 北京化工大学学报(自然科学版), 2011, 38(2): 134-138.

  12. 12. Daubechies, I., Han, B., Ron, A. and Shen, Z.W. (2003) Framelets: MRA-Based Constructions of Wavelet Frames. Applied and Computational Harmonic Analysis, 14, 1-46. https://doi.org/10.1016/S1063-5203(02)00511-0

期刊菜单