7-6 小波分析之小波变换

7-6 小波分析之小波变换

1.小波变换

1.小波变换

\qquad 小波变换是在伸缩和平移后的小波上分解信号。小波$\psi$是一个均值为零$L^2(R)$的函数$${\int }_{-\infty }^{+\infty }\psi (x)dx=0$$它有规范化的范数$\Vert \psi \Vert =1$，且能量集中在$t=0$的邻域内。对$\psi$伸缩$s$倍尺度，平移$u$个单位就得到一个时-频原子字典 D = { ψ u , s = 1 s ψ ( t − u s ) } u ∈ R , s ∈ R + \mathcal D=\{\psi_{u,s}=\frac{1}{\sqrt s}\psi(\frac{t-u}{s})\}_{u\in\mathbb R,s\in\mathbb R^+} D={ψu,s​=s ​1​ψ(st−u​)}u∈R,s∈R+​这些原子仍然有规范化的范数$\Vert {\psi }_{u,s}\Vert =1$.$f\in L^2(R)$在时间$u$和尺度$s$处的小波变换为$$Wf(u,s)=<f,{\psi }_{u,s}>={\int }_{-\infty }^{+\infty }f(t)\frac{1}{\sqrt{s}}{\psi }^{\ast }(\frac{t-u}{s})dt$$.

2.信号重构

信号重构： 设$\psi \in L^2(R)$,是实函数且满足$${\mathcal{C}}_{\psi }={\int }_0^{+\infty }\frac{\hat{\psi }|(\omega ){|}^2}{\omega }d\omega <+\infty$$则对任意的$f\in L^2(R)$有$$f(t)=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }{\int }_{-\infty }^{+\infty }Wf(u,s)\frac{1}{\sqrt{s}}\psi (\frac{t-u}{s})du\frac{ds}{s^2}$$及$${\int }_{-\infty }^{+\infty }|f(t){|}^{2}dt=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }{\int }_{-\infty }^{+\infty }|Wf(u,s){|}^{2}du\frac{ds}{s^2}$$.
证明：记${\psi }_s(t)=s^{-1/2}\psi (t/s)$,带入$Wf(u,s)$得$$Wf(u,s)={\int }_{-\infty }^{+\infty }f(t)\frac{1}{\sqrt{s}}{\psi }^{\ast }(\frac{t-u}{s})dt=f\ast {\overline{\psi }}_s(u)$$其中，${\overline{\psi }}_s(t)=s^{-1/2}\psi (-t/s)$,令$$b(t)=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }{\int }_{-\infty }^{+\infty }Wf(u,s)\frac{1}{\sqrt{s}}\psi (\frac{t-u}{s})du\frac{ds}{s^2}$$$$=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }Wf(\cdot ,t)\ast {\psi }_s(t)\frac{ds}{s^2}$$$$=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }f\ast {\overline{\psi }}_s\ast {\psi }_s(t)\frac{ds}{s^2}$$对$b(t)$做傅里叶变换得$$\hat{b}(\omega )=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }\hat{f}\sqrt{s}\hat{\psi }(s\omega {)}^{\ast }\sqrt{s}\hat{\psi }(s\omega )\frac{ds}{s^2}=\frac{\hat{f}(\omega )}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }|\hat{\psi }(s\omega ){|}^2\frac{ds}{s}$$,由于$\psi (t)$是实函数，由7-1节性质可知$|\hat{\psi }(-s){|}^2=|\hat{\psi }(s)|$.做变量替换$\xi =s\omega$可得$$\hat{b}(\omega )=\frac{\hat{f}(\omega )}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }\frac{\hat{\psi }|(\xi ){|}^2}{\xi }d\xi$$,又因为$${\mathcal{C}}_{\psi }={\int }_0^{+\infty }\frac{\hat{\psi }|(\omega ){|}^2}{\omega }d\omega <+\infty$$,所以$\hat{b}=\hat{f}$,于是有$b(t)=f(t)$得证。
下证第二个式子：由$Wf(u,s)=f\ast {\overline{\psi }}_s(u)$可得$$\widehat{Wf(u,s)}=\hat{f}\sqrt{s}\hat{\psi }(s\omega )$$于是由Plancherel公式有$$\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }{\int }_{-\infty }^{+\infty }|Wf(u,s){|}^{2}du\frac{ds}{s^2}=\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{f}(\omega ){|}^2|\sqrt{s}\hat{\psi }(s\omega ){|}^{2}du\frac{ds}{s^2}$$由Fubini定理可得$$\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }{\int }_{-\infty }^{+\infty }|Wf(u,s){|}^{2}du\frac{ds}{s^2}=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{f}(\omega ){|}^2(\frac{1}{{\mathcal{C}}_{\psi }}{\int }_0^{+\infty }|\sqrt{s}\hat{\psi }(s\omega ){|}^2\frac{ds}{s^2})d\omega =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|f(t){|}^{2}dt$$得证。

3.解析小波

解析小波： 如果一个函数$f_a\in L^2(R)$的傅里叶变换在负频率处的值为0，则称这个函数是解析的，即$\hat{f}(\omega )=0,\omega <0.$,解析函数必是复的，但完全由其实部刻画。事实上，它实部的傅里叶变换为$$\hat{f}(\omega )=\frac{{\hat{f}}_a(\omega )+{\hat{f}}_a^{\ast }(-\omega )}{2}$$,这个关系可被转换为$$\hat{f}(\omega )=\{\begin{array}{ll}2\hat{f}(\omega )& \text{若}\omega \ge 0\\ 0& \text{若}\omega <0\end{array}$$.

2.小波基

\qquad 我们可以构造小波$\psi$,使得伸缩平移函数族 { ψ j , n ( t ) = 1 2 j ψ ( t − 2 j n 2 j ) } ( j , n ) ∈ Z 2 \{\psi_{j,n}(t)=\frac{1}{\sqrt {2^j}}\psi(\frac{t-2^jn}{2^j})\}_{(j,n)\in Z^2} {ψj,n​(t)=2j ​1​ψ(2jt−2jn​)}(j,n)∈Z2​是$L^2(R)$的一组标准正交基。