当前位置:高等教育资讯网  >  中国高校课件下载中心  >  大学文库  >  浏览文档

《小波分析讲义》第三章 连续小波变换

资源类别:文库,文档格式:PDF,文档页数:11,文件大小:190.48KB,团购合买
本章讨论连续小波变换。在第一节将引入小波的概念,小波变换与重构公式。第二节将介绍对尺度进行离散化后的二进小波变换。第三节将用小波变换刻划函数的正则性。
点击下载完整版文档(PDF)

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 第三章连续小波变换 本章讨论连续小波变换。在第一节将引入小波的概念,小波变换与重构公式。第二节将介绍 对尺度进行离散化后的二进小波变换。第三节将用小波变换刻划函数的正则性, 1连续小波变换 设∈L1(R)∩L2(R).我们引入平移,伸缩和调制的概念 定义:对s∈R\{0},t∈R,定义vsb(t)为 (t) v(t)由三种运算构成:i)平移b,位置参数,1)伸缩s,尺度(或1/频率),i)调制左,对峰 值的调整 例1:考虑Haar小波 v()={-1,1/2<t≤ 0,其它 经伸缩与平移后的函数,如20(t),1/20(t)及to1(t)如下:图形 例2设()=/为Gas函数。置o(t)=-0()=h(1-2)e-/2.它的 Fourier变 换为v() /2 从图形可以看出:当s增加时,小波函数vo(t)的“宽度”增加,因而在时域上的分辨率 变差,然而ψ(ω)的宽度变窄,故频率分辨率增强。反之亦然 定理1:经伸缩与平移后,υsb(t)的模保持不变,即 ‖次sb=‖(这部分地解释为什么引入因子 证:作替换u=(t-b)/s,则有 ‖2=∫b(t)2dt ∫(=2d ∫(u)|2 定义:函数ψ∈L2(R)(∩L1(R)称为母小波( mother wavelet),若

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 1 第三章 连续小波变换 本章讨论连续小波变换。在第一节将引入小波的概念,小波变换与重构公式。第二节将介绍 对尺度进行离散化后的二进小波变换。第三节将用小波变换刻划函数的正则性。 1 连续小波变换 设ψ ∈ L 1 (R) ∩ L 2 (R). 我们引入平移,伸缩和调制的概念。 定义: 对s ∈ R\ {0}, t ∈ R,定义ψs,b(t)为 ψs,b(t) = 1 p |s| ψ( t − b s ). (1) ψs,b(t)由三种运算构成:i) 平移b, 位置参数,ii) 伸缩s,尺度(或1/频率), iii) 调制√ 1 |s| , 对峰 值的调整. 例1: 考虑Haar小波 ψ(t) =    1, 0 ≤ t < 1/2, −1, 1/2 < t ≤ 1, 0, 其它. 经伸缩与平移后的函数,如ψ2,0(t),ψ1/2,0(t)及ψ0,1(t) 如下:图形 例2: 设θ(t) = √ 1 2π e −t 2/2为Gauss函数。置ψ(t) = −θ 00(t) = √ 1 2π (1 − t 2 )e −t 2/2。它的Fourier变 换为ψˆ(ω) = ω 2 e −ω 2/2 . 从图形可以看出:当|s|增加时,小波函数ψs,0(t)的“宽度”增加,因而在时域上的分辨率 变差,然而ψds,0(ω)的宽度变窄,故频率分辨率增强。反之亦然。 定理1:经伸缩与平移后,ψs,b(t)的模保持不变,即 kψs,bk = kψk ( 这部分地解释为什么引入因子 1 p |s| ). 证:作替换u = (t − b)/s,则有 kψs,bk 2 = R R |ψs,b(t)| 2 dt = R R 1 |s| ¯ ¯ψ( t−b s ) ¯ ¯ 2 dt = R R |ψ(u)| 2 du = kψk 2 . 定义:函数ψ ∈ L 2 (R) (∩L 1 (R))称为母小波(mother wavelet),若

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 2 (2) 公式(2)称为小波函数的可容许性条件( admissible condition 由于如∈L1故它的 Fourier变换v(u)连续。由(2)式,推得 (u)2 要此积分收敛,必须有:v(0)=0.也即 v(t)dt=0 上式表明,v的正的部分与负的部分与t轴所围的面积是相同的,它的图象是振荡的。类似 于正弦波。我们将对v的光滑性以及局部性作出进一步的限制,因而得到局部化的波形,故 有“小波”这一名称 例:1)高斯函数的差(DOG) v(1)=a-2e-1122)-e-/2,a6时,第二项是数值可忽略的,因而它变为 /2 v(u)=e(-)2. 3)样条函数 3t2-2t,0≤t≤1/2 v(t) (t-1)2,1/2≤t v(t)为奇函数 我们考察从到b后,它的时一频( Heisenberg)窗口的变化。按定义有(不妨设s> 0)。对平均位置°有

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 2 Cψ = Z R ¯ ¯ ¯ψˆ (ω) ¯ ¯ ¯ 2 |ω| dω 0)。对平均位置u s有

Lecture Notes on Wavelets, Chapter 3. by D Q. Dai, 2003 广t|(+)dt b, 即平移了b后再伸缩了s倍。而在频率上,有 广ub()d 1 s 从而它们的平均宽度分别满足 (o2)2 ∫(t-)2losb(t)2dt ∫(t-(su+b)21(=2)dt (o5)2 5)2 sw(sw 由于01.08=01,0,所以 Heisenberg盒子的面积保持不变。但它们的形状随着s的变化而 变化。这是与窗口 Fourier变换不相同的地方 我们现在考虑小波变换 定义:设v为母小波.∫∈L2(R),定义它的小波变换为 (Wf)(s, b)=(, s, b) ∫f(t)sb(t)dt 利用 Plancherel公式,小波变换也可以在频率域内写为 (Wf)(s, b) 我们一般要求小波函数在时域与频域都有良好的局部性,例如,紧支集、指数衰减等。 我们考察在尺度s,位置b的小波变换,在时域与频域从函数∫所抽取得信息

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 3 u s = 1 kψs,bk 2 R +∞ −∞ t|ψs,b(t)| 2 dt = 1 kψk 2 R +∞ −∞ t 1 s ¯ ¯ψ( t−b s ) ¯ ¯ 2 dt = su + b, 即平移了b后再伸缩了s倍。而在频率上,有 ξ s = 1 2πkψs,bk R +∞ −∞ ω ¯ ¯ ¯ψˆ s,b(ω) ¯ ¯ ¯ 2 dω = 1 2πkψk 2 R +∞ −∞ ωs ¯ ¯ ¯ψˆ(sω) ¯ ¯ ¯ 2 dω = 1 s 1 2πkψk 2 R +∞ −∞ ω ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 dω = ξ s 从而它们的平均宽度分别满足 (σ s t ) 2 = 1 kψs,bk 2 R (t − u s ) 2 |ψs,b(t)| 2 dt = 1 kψk 2 R (t − (su + b))2 1 s ¯ ¯ψ( t−b s ) ¯ ¯ 2 dt = s 2 (σt) 2 和 (σ s ω ) 2 = 1 kψs,bk 2 R (ω − ξ s ) 2 ¯ ¯ ¯ψˆ s,b(ω) ¯ ¯ ¯ 2 dω = 1 kψk 2 R (ω − ξ s ) 2 s ¯ ¯ ¯ψˆ(sω) ¯ ¯ ¯ 2 dω = (σω) 2 s 2 . 由于σ s t .σs ω = σt .σω,所以Heisenberg盒子的面积保持不变。但它们的形状随着s的变化而 变化。这是与窗口Fourier变换不相同的地方。 我们现在考虑小波变换。 定义:设ψ为母小波. f ∈ L 2 (R),定义它的小波变换为 (W f) (s, b) = hf, ψs,bi = R R f(t)ψs,b (t)dt = R R f(t)√ 1 |s| ψ ¡ t−b s ¢ dt. 利用Plancherel公式,小波变换也可以在频率域内写为 (W f) (s, b) = 1 2π Z R ˆf(ω) p |s|ψˆ (sω)e ibωdω. (4) 我们一般要求小波函数在时域与频域都有良好的局部性,例如,紧支集、指数衰减等。 我们考察在尺度s,位置b的小波变换,在时域与频域从函数f所抽取得信息

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 由于(Wf)(s,b)=∫f(t)sb(t)d,在时域上所利用的f的信息集中在时刻b,宽度为o st的位置上。它在时域上随着参数b的改变而作平移。当尺度s变小时,o也变小,从 而υb的分辨率增强,(W∫)(s,b)提取愈来愈局部的细节。当尺度增加时,σ增大,v。b的分 辨率减弱,(Wf)(s,b)获取的是∫的整体的性质。尺度s类似于地图绘制中所用的“比例尺” 高比例尺可获得整体的视图,如世界地图;低比例尺对应于象城市一类的局部的地图。这也 形象的解释位置参数的作用 在频域上,由于(Wf)(s,b)=1/(2)∫f(u)vb(u)du,所利用的f的频率宽度为os σ/s,与尺度s成反比。所以当s变小时,小波变换的频域分辨率变差,而当s增大时,它得 到增强。由于小波变换不是 Fourier变换,直观上尺度s与频率成反比关系。因而,我们可以 看出,小波变换在低频时,有较好的分辨率,这是以牺牲在时域的分辨率为代价而得到的。 在高频的时候,(对应于小尺度s),小波变换的分辨率变差。 调制参数1/√冈是由于归一化要求从数学上所导出的。它本身也有物理意义。比如显微 镜,当想知道细节时,你不能看到一个小窗口以外的。在此窗口内细节的强度都增强了,这 对应于|b(t)的幅值,即增加1/√s从这种意义上,小波变换可看作成“数学显微镜 定理2:设∫∈L2(R).则对固定的s≠0,有 wf(s,)2≤√‖1f2 证:设(1)=(-2).则有 =√问 这时(W)(,b)=(f*)(b)。由Youg不等式得 (Wf)(s,)‖ ‖f2=√s‖l1·‖f 定理3:( Calderon(1964), Grossman, Morlet(1984):设f,9∈L2().则有 wf)(s, b)(Wg)(s, b)db 且在L2意义下,有 f()=C//w(s,b)v() dsdb 证:由(4),(Wf)(s,b)关于b的 Fourier变换Wf(s,)=f()√/lsv(s).将 Parseval等式应用 于关于b的积分,得

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 4 由于(W f)(s, b) = R R f(t)ψs,b(t)dt ,在时域上所利用的f的信息集中在时刻b,宽度为σ s t = sσt 的位置上。它在时域上随着参数b的改变而作平移。当尺度s变小时,σ s t也变小,从 而ψs,b的分辨率增强,(W f)(s, b)提取愈来愈局部的细节。当尺度增加时,σ s t增大,ψs,b的分 辨率减弱,(W f)(s, b)获取的是f的整体的性质。尺度s类似于地图绘制中所用的“比例尺”。 高比例尺可获得整体的视图,如世界地图;低比例尺对应于象城市一类的局部的地图。这也 形象的解释位置参数的作用。 在频域上,由于(W f)(s, b) = 1/(2π) R R ˆf(ω)ψˆ s,b (ω)dω,所利用的 ˆf的频率宽度为σ s ω = σω/s,与尺度s成反比。所以当s变小时,小波变换的频域分辨率变差,而当s增大时,它得 到增强。由于小波变换不是Fourier变换,直观上尺度s与频率成反比关系。因而,我们可以 看出,小波变换在低频时,有较好的分辨率,这是以牺牲在时域的分辨率为代价而得到的。 在高频的时候,(对应于小尺度s),小波变换的分辨率变差。 调制参数1/ p |s|是由于归一化要求从数学上所导出的。它本身也有物理意义。比如显微 镜,当想知道细节时,你不能看到一个小窗口以外的。在此窗口内细节的强度都增强了,这 对应于|ψs,b(t)| 的幅值,即增加1/ p |s|。从这种意义上,小波变换可看作成“数学显微镜”。 定理2:设f ∈ L 2 (R). 则对固定的s 6= 0,有 kW f(s, ·)k2 ≤ p |s| kψk1 · kfk2 (5) 证:设ψ ∼ (t) = √ 1 |s| ψ ¡ − t s ¢ . 则有 ° ° ° ° ψ ∼ ° ° ° ° 1 = p |s| kψk1 这时(W f) (s, b) = µ f ∗ ψ ∼ ¶ (b)。由Young不等式得 k(W f) (s, ·)k2 ≤ ° ° ° ° ψ ∼ ° ° ° ° 1 kfk2 = p |s| kψk1 · kfk2 . 定理3: (Calderon(1964), Grossman, Morlet(1984)):设f, g ∈ L 2 (R). 则有 Z R Z R (W f) (s, b)(W g) (s, b)dbds s 2 = Cψ hf, gi (6) 且在L 2意义下,有 f (t) = C −1 ψ Z R Z R (W f) (s, b) ψs,b (t) dsdb s 2 . (7) 证: 由(4),(W f)(s, b)关于b的Fourier变换W f d(s, ξ) = ˆf(ξ) p |s|ψˆ(sξ). 将Parseval等式应用 于关于b的积分,得

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 (.(m9)b学=(W(,HW)空 ∫∫f()Vs(s)0(5)√sv(s)dk R R ∫∫(9(5)()d 2∫f()0(5)( k()-ds)ds g∫f()0()d=C, 这就证明了(6)式 为证明(7)式,考虑积分 c// (Wf)(s, b)vs, b(t) dsdb s1≤|s≤s2 ≤B 当81→0,s2,B→∞时的极限。记△1={(s,b):S1≤|s≤s2,|b≤B},△2=R2-△1我们 有 ∫∫△,(Wf)(s,b)kb()s2 ≤spc-/2wb94 ≤C(△fsb2)2.(sbPe ≤scJ∫△2(f(s,b)2≠)(∫Wgs,b)2e)2 ≤(∫J△2(W)(s.6)24)12 上式收敛于零,因为由(6)式,积分C()(s,b)24收敛 (8)式中的积分当与g∈L2作内积后 W(6)<>=≤△W9>学 ≤∫∫,‖f·b·‖b· =4B(-是)f·|l·啡‖ 从而(8)式属于L2(R) (7)式表明,从函数∫的小波变换系数Wf,可以通过小波v。b重构函数∫。分解时用的小 波可与重构时所用的小波不同,如下面的定理所述的 定理4:设的,如2∈D(,满足C1的=∫当巴d≠,2可微且∈D2(R),2()∈ L1(R),v1(0)=2(0)=0。若f∈L1()且有界,则在f的连续点t,有

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 5 R R R R (W f)(s, b)(wg)(s, b)db ds s 2 = 1 2π R R R R (W fˆ)(s, ξ)W gˆ (s, ξ)dξ ds s 2 = 1 2π R R R R ˆf(ξ) p |s|ψˆ(sξ)gˆ(ξ) p |s|ψˆ(sξ)dξ ds s 2 = 1 2π R R R R ˆf(ξ)gˆ(ξ) ¯ ¯ ¯ψˆ(sξ) ¯ ¯ ¯ 2 dξ ds |s| = 1 2π R R ˆf(ξ)gˆ(ξ)(R R |ψˆ(sξ)| 2 |ξ| ds)dξ = Cψ 2π R R ˆf(ξ)gˆ(ξ)dξ = Cψ , 这就证明了(6)式。 为证明(7)式,考虑积分 1 Cψ Z Z s1 ≤ |s| ≤ s2 |b| ≤ B (W f)(s, b)ψs,b(t) dsdb s 2 (8) 当s1 → 0, s2, B → ∞时的极限。记∆1 = {(s, b) : s1 ≤ |s| ≤ s2, |b| ≤ B}, ∆2 = R 2 − ∆1 我们 有 ° ° °f − 1 Cψ R R ∆1 (W f)(s, b)ψs,b(t) dsdb s 2 ° ° ° = sup kgk≤1 ¯ ¯ ¯ ¯ ¯ ¯ ≤ sup kgk≤1 ¯ ¯ ¯ Cψ −1 R R ∆2 (W f)(s, b)(W g)(s, b) dsdb s 2 ¯ ¯ ¯ ≤ sup kgk≤1 Cψ −1 ( R R ∆2 ¯ ¯ ¯W f(s, b2 dsdb s 2 ) ¯ ¯ ¯) 1/2 · ( R R ∆2 |W g(s, b)| 2 dsdb s 2 ) 1/2 ≤ sup kgk≤1 C −1 ψ ( R R ∆2 |(W f)(s, b)| 2 dsdb s 2 ) 1/2 ( R R R2 |W g(s, b)| 2 dsdb s 2 ) 1/2 ≤ (C −1 ψ R R ∆2 |(W f)(s, b)| 2 dsdb s 2 ) 1/2 上式收敛于零,因为由(6)式,积分Cψ −1 RR R2 |(W f)(s, b)| 2 dsdb s 2 收敛。 (8)式中的积分当与g ∈ L 2 作内积后。 ¯ ¯ ¯ R R ∆1 (W f)(s, b) dsdb s 2 ¯ ¯ ¯ ≤ R R ∆1 |W f| · || dsdb s 2 ≤ R R ∆1 kfk · kψs,bk · kψs,bk · kgk dsdb s 2 = 4B( 1 A1 − 1 A2 ) kfk · kgk · kψk 2 . 从而(8)式属于L 2 (R)。 (7)式表明,从函数f的小波变换系数W f,可以通过小波ψs,b重构函数f。分解时用的小 波可与重构时所用的小波不同,如下面的定理所述的: 定理4:设ψ1,ψ2 ∈ L 1 (R),满足Cψ1,ψ2 = R ψˆ 1(ω)ψˆ 2(ω) |ω| dω 6= 0, ψ2可微且ψ 0 2 ∈ L 2 (R), ψ2(t) ∈ L 1 (R),ψˆ 1(0) = ψˆ 2(0) = 0。若f ∈ L 1 (R)且有界,则在f的连续点t,有

Lecture Notes on Wavelets, Chapter 3. by D Q. Dai, 2003 ds f(t)=Cel tz lim s1≤|sl≤s2- 证:见 DAubechies,pp.28-30 2二进小波变换 在定理3和4中,利用了所有的尺度-∞0,对函数v将有严 格一些的要求。即有 du<+∞o (当是实值函数时,由于v(-u)=v(u),故第二个等式自动成立。) 我们有类似于定理3的重构公式 =C/=/w0 若进一步离散尺度s,它只取二进数值s=2(∈2),自然对于函数会有更多的限制 定义:设函数如∈L2(R),若存在正常数A,B(0<A≤B<+∞),使得 ∑l(2-u)|≤B a. e 则称为二进小波( dyadic wavelet) 定理5:设ψ为二进小波,若(9)式成立,则有 Aln 2< dw< B 2 而且,如果A=B,则 dw=2Aln 2 证明:由于 (2-Ja

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 6 f(t) = C −1 ψ1,ψ2 lims1→0 s2→∞ Z ds s 2 s1≤|s|≤s2 Z∞ −∞ ψ2,s,b(t)db 证:见I.Daubechies,pp. 28-30。 2 二进小波变换 在定理3和4中,利用了所有的尺度−∞ 0,对函数ψ将有严 格一些的要求。即有 Cψ = Z∞ 0 ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ |ω| dω = Z 0 −∞ ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ |ω| dω < +∞. (当ψ是实值函数时,由于ψˆ(−ω) = ψˆ(ω),故第二个等式自动成立。) 我们有类似于定理3的重构公式: f = C −1 ψ Z∞ 0 ds s 2 Z +∞ −∞ (W f)(s, b)ψs,b(t)db. 若进一步离散尺度s,它只取二进数值s = 2j (j ∈ Z),自然对于函数ψ会有更多的限制。 定义:设函数ψ ∈ L 2 (R), 若存在正常数A, B(0 < A ≤ B < +∞),使得 A ≤ X j ¯ ¯ ¯ψˆ(2−jω) ¯ ¯ ¯ 2 ≤ B, a.e. (9) 则称ψ为二进小波(dyadic wavelet) 定理5:设ψ为二进小波,若(9)式成立,则有 A ln 2 ≤ Z∞ 0 ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 ω dω Z 0 −∞ ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 ω dω ≤ B ln 2. 而且,如果A = B,则 Cψ = Z +∞ −∞ ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 |ω| dω = 2A ln 2 证明:由于 Z 2 1 ¯ ¯ ¯ψˆ(2−jω) ¯ ¯ ¯ 2 ω dω = 2−j+1 Z 2−j ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 ω dω

Lecture Notes on Wavelets, Chapter 3. by D Q. Dai, 2003 对(9)式,每项用ω相除,并在区间(1,2)上积分,得 j+1 所以 Aln 22(2s≤2=‖ 即 As∑厂 igito2-v(2-w dsb, gEL(R) 对任凶∈R,取函数g,使得正 ( Gauss函数),并令a→0+,我们 得(9) 定理6解释了,为什么条件(9)被称为稳定性条件

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 7 对(9)式,每项用ω相除,并在区间(1, 2)上积分,得 A Z 2 1 dω ω ≤ X j Z 2−j+1 2−j ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 ω dω ≤ B Z 2 1 dω ω 所以 A ln 2 ≤ Z ∞ 0 ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 ω dω ≤ B ln 2. 类似地,用−ω除,并在区间(−2, −1)上积分,得 A ln 2 ≤ Z 0 −∞ ¯ ¯ ¯ψˆ(ω) ¯ ¯ ¯ 2 −ω dω ≤ B ln 2. 对二进小波ψ,引入二进小波变换如下: (Wjf)(b) = ­ f(t), 2 jψ(2j (t − b))® 注意,在小波前的常数为2 j。 容易验证,Wjf关于b的Fourier变换为 Wdjf(ξ) = ˆf(ξ)ψˆ(2−j ξ). (10) 定理6: 稳定性条件(9)等价于 A kfk 2 ≤ X j kwjfk 2 ≤ B kfk 2 , ∀f ∈ L 2 (R). (11) 证明:在(9)式两边乘以 ˆf(ξ),利用(10)式和Plancherel公式,我们推得(11)。反之,若(11)式 成立,由(10)式,我们有 A 2π ° ° ° ˆf ° ° ° 2 ≤ X 1 2π ° ° ° ˆf(ξ)ψ(2−j ξ) ° ° ° 2 ≤ B 2π ° ° ° ˆf ° ° ° 2 , 即 A ≤ X j Z +∞ −∞ |g(ω)| 2 kgk 2 ¯ ¯ ¯ψˆ(2−jω) ¯ ¯ ¯ 2 dω ≤ B, ∀g ∈ L 2 (R). 对任ω0 ∈ R,取函数g,使得|g(ω)| 2 kgk 2 = 1 2 √ πα e − (ω−ω0) 2 4α (Gauss函数),并令α → 0 +,我们 得(9)。 定理6解释了,为什么条件(9)被称为稳定性条件

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 定义:函数∈L2(R)称为二进小波v的二进对偶,若 f()=∑/(w0)20(2(t-b)dbf∈D2(R) 二进对偶ψ也许不是唯一的,就象对偶框架的情形一样。对于稳定的二进小波,我们总 可以构造它的对偶二进小波 定理7:设函数v(t)是由它的 Fourier变换v(u)所定义 则v是稳定的二进小波,即有 ∑|(2-)≤ 而且,U是v的二进对偶 证明:注意到,(12)式的分母是2伸缩不变的,因而我们有 ∑(2-)P ∑;v(2-)P2 ∑kv(2-k)|2)2∑kv(2-b)2 由(9)推得(13) 下面证明是v的二进对偶。我们有 ∑,」W)(b)·20(2(t-b) ∑,1/(2x)f()(2-1)v(2-)e“d t 其中,我们利用了公式h1*h2=F-1(h1h2),并将第一式写为卷积的形式 3 Holder正则性刻划 连续小波变换可以用来刻划函数的光滑性 定义:函数∫在R上是指数a(0|≤Cls 证明:由于∫v(1)dt=0,我们有

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 8 定义:函数ψ˜ ∈ L 2 (R)称为二进小波ψ的二进对偶,若 f(t) = X j Z +∞ −∞ (Wjf)(b) · 2 jψ˜(2j (t − b))db, ∀f ∈ L 2 (R). 二进对偶ψ˜也许不是唯一的,就象对偶框架的情形一样。对于稳定的二进小波,我们总 可以构造它的对偶二进小波。 定理7: 设函数ψ˜(t)是由它的Fourier变换ψ ˆ˜(ω)所定义 ψ ˆ˜(ω) = ψˆ(ω) P k ¯ ¯ ¯ψˆ(2−kω) ¯ ¯ ¯ 2 . (12) 则ψ˜是稳定的二进小波,即有 1 B ≤ X j ¯ ¯ ¯ψ ˆ˜(2−jω) ¯ ¯ ¯ ≤ 1 A , (13) 而且,ψ˜是ψ的二进对偶。 证明: 注意到,(12)式的分母是2 j伸缩不变的,因而我们有 X j |ψ ˆ˜(2−jω)| 2 = P j |ψˆ(2−jω)| 2 ( P k |ψˆ(2−kω)| 2 ) 2 = 1 P k |ψˆ(2−kω)| 2 . 由(9)推得(13). 下面证明ψ˜是ψ的二进对偶。我们有 P j R ∞ −∞(Wjf)(b) · 2 jψ˜(2j (t − b))db = P j 1/(2π) R ∞ −∞ ˆf(ξ) ˆ ψ(2−j )ψ ˆ˜(2−j ξ)e itξdξ = 1/(2π) R ∞ −∞ ˆf(ξ)e itξdξ = f(t), 其中,我们利用了公式h1 ∗ h2 = F −1 (hˆ 1hˆ 2), 并将第一式写为卷积的形式。 3 H¨older正则性刻划 连续小波变换可以用来刻划函数的光滑性。 定义: 函数f在R上是指数α(0 | ≤ C|s| α+1/2 . 证明:由于R ψ (t) dt = 0,我们有

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 9 f,vb)=1-/v ((t)-f(b))dt 所以 《f,vsb ≤|s∫|(t32)|·k|t-bdt ≤k|s+∫lv(t) tl dt 我们也有如下的反向结果 定理9:设是紧支集的小波函数。假定∫∈L(R),连续且有界。若对某α∈(O,1),它的小 波变换满足 ≤k|s+ 则∫是a- Holder连续的。 证:取连续可微,紧支的函数v2,使得v2(0)=0.C2=1则由定理4有 f(t) 广厂1 s -oo db -oo dyl f(y)lv(y-b)s)lv2((t+h-b/s)-v2((t-b)/s) 由于|2(2+1)-2(2)≤团,且对某个R>0,supp,supp2c[R,F],因而上式不超过 CIhIJisp21lsr-4ds jio-tlklslR+1 db y-bkIsR If(y)ldy ≤C"hlJp21s=as y-t≤2|sR+1 f(yldy ≤C"hf1s-3(4s+2)/ds ≤C"|h

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 9 hf, ψs,bi = |s| − 1 2 Z ψ µ t − b s ¶ (f(t) − f(b))dt 所以 |hf, ψs,bi| ≤ |s| − 1 2 R ¯¯ψ ¡ t−b s ¢¯¯ · k |t − b| α dt ≤ k |s| α+ 1 2 R |ψ(t)| |t| α dt ≤ k 0 |s| α+ 1 2 . 我们也有如下的反向结果: 定理9:设ψ是紧支集的小波函数。假定f ∈ L 2 (R),连续且有界。若对某α ∈ (0, 1),它的小 波变换满足 || ≤ k |s| a+ 1 2 则f是α-H¨older连续的。 证:取连续可微,紧支的函数ψ2,使得ψˆ 2(0) = 0, Cψ,ψ2 = 1则由定理4有 f(t) = Z +∞ −∞ ds s 2 Z +∞ −∞ ψ2,s,b(t)db. 我们将关于s的积分分为两部分,并分别记为fL(t)( 大尺度s, |s| ≥ 1)和fs(t)( 小尺 度s, |s| ≤ 1)。 我们将证明fL(t)总是正则的, 1)fL(t)关于t一致有界, |fL(t)| ≤ R |s|≥1 ds s 2 R +∞ −∞ kfk · kψs,bk · |ψ2,s,b(t)|db ≤ C R |s|≥1 ds s 2 R +∞ −∞ |s| − 1 2 |ψ2 ¡ t−b s ¢ |db ≤ C R |s|≥1 |s| − 3 2 ||ψ2||L1 ds ≤ C 0 0, suppϕ,suppψ2 ⊂ [−R, R], 因而上式不超过 ≤ C|h| R |s|≥1 |s| −4ds R |b−t|≤|s|R+1 db R |y−b|≤|s|R |f(y)|dy ≤ C 0 |h| R |s|≥1 |s| −3ds R |y−t|≤2|s|R+1 |f(y)|dy ≤ C 00|h|kfk R |s|≥1 |s| −3 (4|s|R + 2)1/2ds ≤ C 000|h|

Lecture Notes on Wavelets, Chapter 3. by D.Q. Dai, 2005 由第一步的结果,当>1式,我们显然有L(t+h)-f(t≤Ch.至此,我们还未用到 关于小波系数的条件 3)fs(t)关于t一致有界 Is(t)l s Jssi s5-oo Illv 2,., (t)ldb ≤Ks卖/p +1/21c|-1/ l2((t-b/s)ldb ≤Kl2+snls-tod 4)当|-1(t+h-b)/s)-2(t-b/)|b ≤s2厂-dse(a2(t+h-b)/s)+12(t-b)/s) h≤|sk≤1 b≤|sR+|h s°C|hb/sdb s Cllvall2 sksh Is1-1sllsds+C"| fIf(t)-f(O)1lsl-1/2l(t-b)s)ldt ≤C∫+s-1(t-b)/s)t ≤Cl+12∫y+b/s|(y)dy ≤C"ls+2∫(y+|b/s9)(y)dy Cs|1/()a+ bg) 反之,我们有 定理1:设小波函数ψ有紧支集。函数∫∈L2(R)连续有界。若对某个>0和a∈(0,1,有 ≤C|s+关于b一致 和 (,)≤C中( log bll

Lecture Notes on Wavelets, Chapter 3, by D.Q. Dai,2003 10 由第一步的结果,当|h| > 1式,我们显然有|fL(t + h) − fL(t)| ≤ C|h|. 至此,我们还未用到 关于小波系数的条件。 3)fs(t)关于t一致有界. |fs(t)| ≤ R |s|≤1 ds s 2 R ∞ −∞ | ||ψ2,s,b(t)|db ≤ K R |s|≤1 ds s 2 R ∞ −∞ |s| α+1/2 |s| −1/2 |ψ2((t − b)/s)|db ≤ Kkψ2k1 R |s|≤1 |s| −1+αds ≤ C ||s| −1/2 |ψ2((t + h − b)/s) − ψ2((t − b)/s)|db ≤ R |s|≤|h| ds s 2 R ∞ −∞ db|s| α (|ψ2((t + h − b)/s)| + |ψ2((t − b)/s)|) + R |h|≤|s|≤1 ds s 2 R |t−b|≤|s|R+|h| |s| αC|h/s|db ≤ C 0kψ2k2 R |s|≤h |s| −2 |s| α |s|ds + C 00|h| R |h|≤|s|≤1 |s| −3+α (|s|R + |h|)ds ≤ C 000|h| α . 定理8和定理9刻划了整体H¨older连续性。对于局部连续性,我们有如下的结果。它们首 先被Holschneider和Tchmitchian用来研究Riemann函数 P∞ n=1 n −2 sin(n 2πt)的可微性。 定理10: 设函数ψ满足 R (1 + |t|)|ψ(t)|dt | ≤ C|s| 1/2 (|s| α + |b| α ). 证:通过平移,我们可以置t0 = 0. 由 R ψ(t)dt = 0, 得 | | = R |f(t) − f(0)||s| −1/2 |ψ((t − b)/s)|dt ≤ C R |t| α |s| −1/2 |ψ((t − b)/s)|dt ≤ C|s| α+1/2 R |y + b/s| α |ψ(y)|dy ≤ C 0 |s| α+1/2 R (|y| α + |b/s| α )|ψ(y)|dy ≤ C 00|s| 1/2 (|s| α + |b| α ). 反之,我们有 定理11: 设小波函数ψ有紧支集。函数f ∈ L 2 (R)连续有界。若对某个γ > 0和α ∈ (0, 1], 有 |hf, ψs,bi| ≤ C |s| r+ 1 2 关于b 一致 和 |hf, ψs,b+t0 i| ≤ C |s| 1 2 µ |s| α + |b| α |log |b||¶

点击下载完整版文档(PDF)VIP每日下载上限内不扣除下载券和下载次数;
按次数下载不扣除下载券;
24小时内重复下载只扣除一次;
顺序:VIP每日次数-->可用次数-->下载券;
共11页,试读已结束,阅读完整版请下载
相关文档

关于我们|帮助中心|下载说明|相关软件|意见反馈|联系我们

Copyright © 2008-现在 cucdc.com 高等教育资讯网 版权所有