基于Sierpinski分形上的多分辨分析
A Multiresolution Analysis Based on Sierpinski Fractal
DOI: 10.12677/PM.2019.98124, PDF, HTML, XML, 下载: 927  浏览: 2,399 
作者: 李岱琪, 李万社:陕西师范大学,数学与信息科学学院,陕西 西安
关键词: Sierpinski垫分形多分辨分析框架Sierpinski Gasket Fractal Multiresolution Analysis Framework
摘要: 根据经典的由迭代仿射函数系统产生的分形的扩张和平移所生成的σ有限测度空间上的多分辨分析的概念,本文讨论了基于Sierpinski垫上的多分辨分析,通过定义酉扩张算子D和酉算子T,并结合分形几何中的测度知识证明了基于Sierpinski垫的膨胀分形集上的多分辨分析。
Abstract: According to the classic theory, the notion of multiresolution analysis based on σ-finite measure spaces built from dilations and translations on a fractal arising from an iterated affine function system. A multiresolution analysis based on Sierpinski gasket is discussed in this paper. By defining unitary dilation operator D and unitary operator T, and combining the measure knowledge in fractal geometry, the multiresolution analysis on the inflated fractal set based on Sierpinski gasket is proved.
文章引用:李岱琪, 李万社. 基于Sierpinski分形上的多分辨分析[J]. 理论数学, 2019, 9(8): 980-988. https://doi.org/10.12677/PM.2019.98124

1. 引言

随着科学研究的推进,小波理论体系不仅可以定义在实数域上,还可以定义在其他域上,比如阿贝尔群、分层树、一些分形集等 [1] [2] [3] [4]。

近年来,小波与分形的结合已经成为科学研究的热点问题。小波分析技术可以有效地揭示分形局域标度性质,小波分析还具有放大和移位的功能,它是从远处到近处观察形体,这与分形的本质是尺度变化相类似。Mallat和Meyer建立的多分辨分析(Multi-Resolution Analysis, MRA)为小波分析奠定了基础,多分辨分析与分形之间也有一定的联系。多分辨分析是从远到近观察形体,先观察其轮廓,再观察其线条,进一步观察物体的细节纹理,这体现了从低分辨到高分辨的思想,对具有自相似性质的分形的观察也是这样,通过从大到小的不同尺度变换,在越来越小的尺度上观察越来越丰富的细节,这也是从低分辨到高分辨的观察过程。

多分辨分析是由S. Mallat引入的。他在空间概念上形象的说明了小波的多分辨特性。1989年,Mallat在探究小波变换多分辨分析理论与图像处理的应用时,受到塔式算法的启发,提出了信号的塔式多分辨分析分解与重构的快速算法,即著名的Mallat算法。MRA形成了构造正交小波基的一个框架,比如常用的Daubechies紧支撑正交小波基可以看作是该框架下的产物 [5]。

D. Dutkay和P. Jorgensen介绍了由迭代仿射函数系统产生的分形的扩张和平移所生成的 σ 有限测度空间上的多分辨分析的概念 [6],但主要是基于一维分形集的情形,比如一维康托尔集上的多分辨分析。本文讨论了基于Sierpinski垫上的多分辨分析,证明了基于Sierpinski垫的膨胀分形集上的多分辨分析,丰富了多分辨分析的内涵。

2. 预备知识

在本文中,用R表示实数集合,Z表示所有整数集合, 表示复数空间。下面先讨论多分辨分析的定义。

定义1 [6] 空间 L 2 ( R ) 中一列闭子空间 { V j } j Z 称为 L 2 ( R ) 的一个依尺度函数 φ 的多分辨分析,如果该序列满足下列条件:

(1) 嵌套性: V j V j + 1 , j Z

(2) 逼近性: j Z V j = { 0 } j Z V j ¯ = L 2 (R)

(3) 伸缩性: f ( x ) V j f ( 2 x ) V j + 1 , j Z

(4) 平移不变性: f ( x ) V 0 f ( x k ) V 0 , k Z

(5) Riesz基存在性:

存在函数 φ V 0 ,使得 { φ ( x k ) } k Z 构成 V 0 的一个Riesz基,即函数序列 { φ ( x k ) } k Z 线性无关,且存在常数A和B,满足 0 < A B < + ,使得对任意的 f ( x ) V 0 ,总存在序列 { c k } k Z l 2 使得

f ( x ) = k = + c k φ ( x k ) A f 2 2 k = + | c k | 2 B f 2 2 ,则称 φ 为尺度函数,并称 φ 生成 L 2 ( R ) 的一个多

分辨分析 { V j } j Z

特别地,若 { φ ( x k ) } k Z 构成 V 0 的一个标准正交基,则称 φ 为正交尺度函数;相应地,称 φ 生成 L 2 ( R ) 的一个正交多分辨分析 { V j } j Z

由一维多分辨分析的张量积空间构造二维多分辨分析,我们给出如下定义:

定义2 [7] [8]

1. 用 L 2 ( R 2 ) 表示平面上的平方可积函数空间,即

f ( x , y ) L 2 ( R 2 ) + + | f ( x , y ) | 2 d x d y <

{ V j 1 } { V j 2 } 是由尺度函数 φ 1 ( x ) φ 2 ( y ) 生成的两个多分辨分析,则可以得到 V j 1 V j 2 的张量积

空间 V j = V j 1 V j 2 。由于 V j 1 的基底为 { 2 j / 2 φ 1 ( 2 j x k ) } V j 2 的基底为 { 2 j / 2 φ 2 ( 2 j y l ) } ,所以 V j 的基底为

{ 2 j φ 1 ( 2 j x k ) φ 2 ( 2 j y l ) } 。对于二元函数 ,引入记号 f j ; k , l ( x , y ) = 2 j f ( 2 j x k , 2 j y l ) ,记 φ ( x , y ) = φ 1 ( x ) φ 2 ( y ) ,则 { φ j ; k , l ( x , y ) : k , l Z } V j 的基底。

这样, { V j } 就形成 L 2 ( R 2 ) 中的一个多分辨分析, 就是相应的尺度函数。

2. 在 L 2 ( R 2 ) 函数空间的一串子空间序列集合 { V j } j Z 称为依尺度函数 φ ( x , y ) 的多分辨分析:

(1) 嵌套性: V j ( x , y ) V j + 1 ( x , y ) , j Z

(2) 逼近性: j Z V j ( x , y ) = { 0 } j Z V j ( x , y ) = L 2 (R2)

(3) 伸缩性: f ( x , y ) V j ( x , y ) f ( 2 x , 2 y ) V j + 1 ( x , y ) , j Z

(4) 平移不变性: f ( x , y ) V 0 ( x , y ) f ( x k , y l ) V 0 ( x , y ) , k , l Z

(5) 规范正交基: φ ( x , y ) = φ 1 ( x ) φ 2 ( y ) V 0 ( x , y ) V 0 ( x , y ) 的一组标准正交基。

W j V j V j + 1 中的正交补空间, { ψ j ; k , l ( x , y ) = 2 j ψ ( 2 j x k , 2 j y l ) } 构成 W j 的基底,则由多分辨分析可以得到两个结果:

V j + 1 = V j W j

L 2 ( R 2 ) = V 0 W 0 W 1 W 2

3. Sierpinski垫上的多分辨分析

给定平面 R 2 上的一个单位正三角形 S 0 ,在每边上取中点,然后两两连结,构成四个边长为1/2的正三角形,把中间的三角形去掉(但是保留其三条边),得到集合 S 1 (见图1),然后在 S 1 的三个三角形的每边取中点,分别连结成边长为1/22的小正三角形,又分别去掉中间一个(同样保留每个三角形的三边),剩下一共有32个三角形。如此不断继续下去,得到一个平面集列:可见 S k 是由 3 k 个边长为1/2k的正三角形组成。 { S k } 的极限集 S = k = 0 S k 称为Sierpinski (谢尔宾斯基)垫片,它是一个著名的分形集 [9] [10]。

S k 的总面积为 3 4 ( 3 4 ) k ,所以当 k 时,S的面积等于如下极限: lim k 3 4 ( 3 4 ) k = 0

我们将原三角形 S 0 视为 L 2 ( R 2 ) 空间中的子空间 V 0 ,把每次去掉的部分子空间记为 W k ,将余下的部分记为 V k (如图2)。

容易看出任意两个不同区域的交集是空集,表示它们的特征函数(该区域的元素集合) C V ( x ) C W ( x ) 彼此正交,简记为 C V ( x ) : = V , C W ( x ) : = W ,则有 V 1 W 1 = { 0 } V 2 W 2 W 1 = { 0 } V N W N W N 1 W 1 = 0

S0 S1S2S3S

Figure 1. Establishment of Sierpinski gasket

图1. Sierpinski垫的形成过程

Figure 2. MRA of Sierpinski gasket

图2. Sierpinski垫上的多分辨分析

我们发现 V k W k 互不相交,说明分别定义在 V k W k 上面的特征函数彼此正交,对于不同的子空间 W k ,它们之间彼此正交,即定义在它们上面的特征函数彼此正交。然而对于不同的子空间 V k V l ,并不是正交的,我们有 V 0 = V 1 W 1 W 1 V 1 V 0 中的正交补空间,依此类推, V 1 = V 2 W 2 V 2 = V 3 W 3 ,因此,

在V上定义尺度函数 φ ( x , y ) ,W上定义小波函数 ψ ( x , y ) 。将 V k 上的基函数 φ m , n ( x , y ) W k 上的基函数 ψ m , n ( x , y ) 结合起来展开空间中的信号 f ( x , y ) L 2 ( R 2 ) 空间的正交分解为 V 0 = V 1 W 1 = V 2 W 2 W 1 = = V N ( k = 1 N W k ) 。相应的信号 f ( x , y ) 的分解过程为

f 0 ( x , y ) = f 1 ( x , y ) + d 1 ( x , y ) = f 2 ( x , y ) + d 2 ( x , y ) + d 1 ( x , y ) = f 3 ( x , y ) + d 3 ( x , y ) + d 2 ( x , y ) + d 1 ( x , y ) = = f N ( x , y ) + k = 1 N d k ( x , y )

我们对 空间赋予多分辨分析的几何特征,将多分辨分析这一抽象的数学定义与实用信号结合起来。

4. 基于Sierpinski垫的膨胀分形集上的多分辨分析

4.1. 定义直角三角形Sierpinski垫

特别地,我们旨在于研究直角三角形Sierpinski垫。

R 2 上构建一个直角三角形 S 0 ,顶点分别是 ( 0 , 0 ) , ( 0 , 1 ) , ( 1 , 0 ) .设对角扩张矩阵 A = ( 2 0 0 2 ) ,令

S 1 = [ A 1 ( S 0 + τ 0 ) ] [ A 1 ( S 0 + τ 1 ) ] [ A 1 ( S 0 + τ 2 ) ]

其中, τ 0 = ( 0 , 0 ) , τ 1 = ( 1 , 0 ) , τ 2 = ( 0 , 1 ) ,则

S 2 = [ A 1 ( S 1 + τ 0 ) ] [ A 1 ( S 1 + τ 1 ) ] [ A 1 ( S 1 + τ 2 ) ]

继续归纳下去,有

S n + 1 = [ A 1 ( S n + τ 0 ) ] [ A 1 ( S n + τ 1 ) ] [ A 1 ( S n + τ 2 ) ]

因此,我们可以得到 R 2 上的紧子集的一个嵌套序列 { S n } n = 0

然后,定义Sierpinski垫分形 S = n = 0 S n (见图3),这个Sierpinski垫S满足自相似关系

A ( S ) = S [ S + ( 1 , 0 ) ] [ S + ( 0 , 1 ) ]

我们知道S的Hausdorff维数是 s = log 3 log 2 ,记 H s ( S ) 表示在S上的s维Hausdorff测度,简记为 H ,则 。更一般地,如果E是S的一个Borel子集,则 H ( A 1 ( E ) ) = 1 3 H ( E )

4.2. 基于Sierpinski垫的膨胀分形集上的多分辨分析

定义一个与Sierpinski垫S有关的一个膨胀分形集

R s = j = ( m , n ) Z 2 [ A j ( S + ( m , n ) ) ]

对于 R s 的每一个Borel子集E,测度满足

H ( A 1 ( E ) ) = 1 3 H (E)

.

在Hilbert空间 L 2 ( R s , H ) 上,构建酉扩张算子D和酉算子T,定义如下:

(1)

T ( m , n ) f ( s , t ) = f ( s m , t n ) . (2)

这些算子满足一个标准的交换关系:

定理1

设D和 { T ( m , n ) : ( m , n ) Z 2 } 是上述定义的酉算子,则

证明:由(1)式和(2)式可得:

T ( m , n ) D ( f ) ( s , t ) = 3 T ( m , n ) f ( 2 s , 2 t ) = 3 f ( 2 s 2 m , 2 t 2 n ) = 3 T ( 2 m , 2 n ) f ( 2 s , 2 t ) = D T ( 2 m , 2 n ) f ( s , t )

所以 T ( m , n ) D = D T ( 2 m , 2 n ) , ( m , n ) Z 2

现在构建一个 L 2 ( R s , H ) 上的多分辨分析。定义一个闭子空间 V 0 L 2 ( R s , H )

V 0 = s p a n ¯ { T ( m , n ) ( X S ) ( m , n ) Z 2 } ,

其中 X S 是Sierpinski垫三角形上的特征函数(对应标准多分辨分析中的尺度函数)。Sierpinski垫S满足自相似关系

A ( S ) = S [ S + ( 1 , 0 ) ] [ S + ( 0 , 1 ) ]

在测度为0的集合上,上述的并集是互相非交的并集,且它的特征函数 X S 满足扩张方程

X S ( A 1 ( s , t ) ) = X S ( s , t ) + X S ( s 1 , t ) + X S ( s , t 1 )

通过构造, V 0 在算子 { T ( m , n ) ( m , n ) Z 2 } 下是不变量。对于每个 j Z ,定义

V j = D j (V0)

可以看出

V 1 = D V 0 = D ( s p a n ¯ { T ( m , n ) ( X S ) ( m , n ) Z 2 } ) = s p a n ¯ { D T ( m , n ) ( X S ) ( m , n ) Z 2 } = s p a n ¯ { T ( m / 2 , n / 2 ) D ( X S ) ( m , n ) Z 2 }

因此, V 0 D ( V 0 ) = V 1 。由此看出,闭子空间 { V j } j = 形成了 L 2 ( R s , H ) 的闭空间的一个递增的嵌套序列。

下面证明子空间 { V j } j = 形成一个多分辨分析

定理2

{ V j } j = 是上述构造的 L 2 ( R s , H ) 的子空间,设D, 是定理1中的酉算子,则有

(1)

(2) T ( m , n ) ( X S ) , X S = δ ( m , n ) , ( 0 , 0 ) , ( m , n ) Z 2

(3) j = V j ¯ = L 2 ( R s , H )

(4) j = V j = { 0 }

证明:

(1) 由Sierpinski垫S的自相似关系和对应的扩张方程可知,因为

X S ( A 1 ( s , t ) ) = X S ( s , t ) + X S ( s 1 , t ) + X S ( s , t 1 ) = X S ( s , t ) + T ( 1 , 0 ) ( X S ( s , t ) ) + T ( 0 , 1 ) ( X S ( s , t ) )

又因为

D ( f ) ( s , t ) = 3 f ( 2 s , 2 t )

D 1 ( 3 f ( 2 s , 2 t ) ) = f ( s , t )

D 1 ( 3 X S ( 2 s , 2 t ) ) = X S ( s , t )

D 1 ( 3 X S ( s , t ) ) = X S ( 2 1 s , 2 1 t )

所以

D 1 ( X S ) = 1 3 [ X S + T ( 1 , 0 ) ( X S ) + T ( 0 , 1 ) ( X S ) ]

(2) 因为

T ( m , n ) ( X S ( s , t ) ) = X S ( s m , t n )

所以

T ( m , n ) ( X S ) , X S = X S ( s m , t n ) , X S = δ ( m , n ) , ( 0 , 0 )

(3) 我们将证明对于任意Hausdorff可测子集 E R s ,其中 H ( E ) < ,则 X E 含于的闭包中。

因为 R s = j = ( m , n ) Z 2 [ A j ( S + ( m , n ) ) ] ,我们可以写一个集合 E = E ( j , m , n ) ,其中

E ( j , m , n ) = E [ A j ( S + ( m , n ) ) ]

容易证明,每个集合 E ( j , m , n ) 的特征函数含于 span { D j T ( m , n ) ( X S ) | j , m , n Z } 的闭包。因此,通过扩张和平移,不失一般性,我们可以假定任意可测集 E S

V 是S中所有子三角形的左下顶点的集合,对于每一个 { v } V ,有一个递减的交集: { v } = n = 1 T n

其中 T n 是S的一个平移的n次扩张。我们可以得到 H ( { v } ) = lim n H ( T n ) = lim n 1 3 n = 0 。由测度的次

可数可加性得, H ( V ) = 0

S = S V ,则 H ( S ) = H ( S ) = 1 。因为S是一个度量空间,所以 S 也是,尽管 S 不再是闭的。

因此,不失一般性,我们假定 E S ,令

T 是由 S 的子Sierpinski垫组成的,它是 S 的子集的一个半代数,也就是说, T 中元素的有限非交的并形成了 S 子集的一个代数。我们用 A 表示这个代数,在 A 上作用一个Hausdorff测度,得到一个 A 上的集值函数,记为 u * ,它满足Caratheodory扩张定理的条件。因此, u * 可以延伸为 S 的所有子集的一个外测度, u * 与代数 A 上的Hausdorff测度 H 一致。因此, u * 决定了可测集上的一个 σ 代数 M ,这个可测集含有包含 A 的最小的 σ 代数 B 。如果将外测度 u * 限制到 M 上,则变成测度,记为u,显然 ( X , M , u ) 是一个完备的测度空间。而且,因为 是有限的,由Caratheodory扩张定理可知,将 A 上的 H 延伸成含有 A 的最小的 σ 代数是唯一的。因此,不论在 B 上还是 M 上,都有 u = H 。剩下的问题是,当构造子集上的外Hausdorff测度时所产生的 σ 代数是否大于利用Caradiodory扩张定理产生的Hausdorff测度 u = H σ 代数 M

注意,集合 T 族除了是 S 子集的半代数外,也是 S 的任意子集的一个Vitali覆盖,也就是说,对于每一个 x E 和每个 δ > 0 ,都有一个子集 T T ,其中 x T , 0 < H ( T ) δ 。利用Vitali覆盖上的Caratheodory扩张定理构造的外测度是度量外测度,因此,对于任意子集A和B,其中 d ( A , B ) > 0 ,满足 u * ( A B ) = u * ( A ) + u * ( B ) 。此外,如果 v * 是度量空间X上的度量外测度,则该外测度下的可测集的 σ 代数包含X的Borel集的 σ 代数。我们将这个结果应用于 ( S , M , u ) 可以推导出 u * 是度量外测度,且 M 包含 S 的Borel集。因此,由 u * 构造的 S 可测集的 σ 代数包含 S 的Borel集以及 S 的开集和闭集。因此,在S的Borel子集上, H 与u是一致的 [11]。

最后,如果利用Caratheodory扩张定理在集合X上用代数 A 构造有限测度u,则给定任意可测集 G X ,存在 A 的一个元R,使得 u ( G Δ R ) < ε

由上述结果以及 S 的可测子集E可知,存在一个 F σ F E ,其中 H ( E Δ F ) = H ( E F ) = 0 。因

为在 S 的Borel集的 σ 代数上, u = H H ( F Δ i = 1 n T i ) < ϵ ,则可推导出存在一个有限集Sierpinski子三角形 T i ,其中 u ( F Δ i = 1 n T i ) < ϵ 。但是,因为 H ( E Δ F ) = 0 ,所以 H ( E Δ i = 1 n T i ) < ϵ 。因此, X E 含于

span { D j T ( m , n ) ( X S ) | j , m , n Z } 的闭包 [12] [13] [14]。

(4) 对于 f j = V j ,如果 ( x , y ) support ( f ) ,则对于每一个 j Z ,存在一些 ( u j , v j ) Z 2 ,使得

( x , y ) A j ( S + ( u j , v j ) )

对于使 x 2 + y 2 < 2 2 j 1 成立的足够大的j,有 ( u j , v j ) { ( 0 , 0 ) , ( 0 , 1 ) , ( 1 , 0 ) } ,且对于这些j, ( u j , v j )

常量。因此, ( x , y ) 一定是 j = A j S j = A j ( S + ( 0 , 1 ) ) j = A j ( S + ( 1 , 0 ) ) 之一。因为 f V j

定是 A j ( S + ( u , v ) ) 上的常量,所以每个并是嵌套的意味着f在这些并上都是常量.因为每个并的测度是无限的,所以这些常量必须都是0 [15]。

因此, { V j } j = 构成了 L 2 ( R s , H ) 上的一个多分辨分析, { T ( m , n ) ( X S ) } 构成了 V 0 的标准正交基。

5. 结束语

由于小波分析可以同时在时域和频域内进行局部化信号,并对信号在不同尺度上进行分解和重构,所以研究分形集上的多分辨分析具有重要的现实意义,为分形集上的小波变换奠定基础。本文首先讨论了基于Sierpinski垫上的多分辨分析,然后通过构建酉扩张算子D和酉算子T,结合分形几何中的测度知识证明了基于Sierpinski垫的膨胀分形集上的多分辨分析。

参考文献

[1] Yang, S.Z. and Cheng, Z.X. (2002) Orthonormal Multi-Wavelets on the Interval [0,1] with Multiplicity. Acta Mathematica Sinica, 45, 789-796.
[2] Nakahira, K. and Miyamoto, A. (2018) Parseval Wavelets on Hierarchical Graphs. Applied and Computational Harmonic Analysis, 44, 414-445.
https://doi.org/10.1016/j.acha.2016.05.004
[3] 师东利, 李万社. 基于二维康托尔二元群的小波变换[J]. 纺织高校基础科学学报, 2017, 30(3): 311-317.
[4] 聂伟平, 师东利, 李万社. 二维康托尔集上的多分辨分析[J]. 西安文理学院学报, 2016, 19(5): 26-29.
[5] 成礼智. 小波的理论与应用[M]. 北京: 科学出版社, 2004: 75-77.
[6] Dutkay, D. and Jorgensen, P.E.T. (2006) Wavelets on Fractals. Revista Matematica Iberoamericana, 22, 131-180.
https://doi.org/10.4171/RMI/452
[7] 方清城. MATLAB R2016a小波分析22个算法实现[M]. 北京: 电子工业出版社, 2018: 338-342.
[8] 沙震, 阮火军. 分形与拟合[M]. 杭州: 浙江大学出版社, 2005: 16-18.
[9] Christensen, O. (2005) Wavelet Frames and Multiresolution Analysis. Pure and Applied Mathematics Series, 272, 73-103.
https://doi.org/10.1201/9781420026511.pt2
[10] Bohnstengel, J. and Kessebohmer, M. (2010) Wavelets for Iterated Function Systems. Functional Analysis, 259, 583-601.
https://doi.org/10.1016/j.jfa.2010.04.014
[11] David, G. and Semmes, S. (1997) Fractured Fractals and Broken Dreams. Clarendon Press, Oxford.
[12] D’Andrea, J., Merrill, K.D. and Packer (2008) Fractal Wavelets of Dutkay-Jorgensen Type for the Sierpinski Gasket Spaces. Communications in Contemporary Mathematics, 451, 69-88.
https://doi.org/10.1090/conm/451/08758
[13] Baggett, L.W., Medina, H.A. and Merrill, K.D. (1999) Generalized Multi-Resolution Analyses and a Construction Procedure for All Wavelet Sets in R^n. Fourier Analysis Applications, 5, 563-573.
https://doi.org/10.1007/BF01257191
[14] Strichartz, R.S. (1997) Piecewise Linear Wavelets on Sierpinski Gasket Type Fractals. The Journal of Fourier Analysis and Applications, 3, 387-415.
https://doi.org/10.1007/BF02649103
[15] Cao, S.P. and Qiu, H. (2017) Some Properties of the Derivatives on Sierpinski Gasket Type Fractals. Constructive Approximation, 46, 319-347.
https://doi.org/10.1007/s00365-017-9385-3