7-5 小波分析之窗口傅里叶变换

7-5 小波分析之窗口傅里叶变换

1.定义

\qquad 从直观上理解，假设信号是一个声音，我们能听到不同时刻的不同声音。所以，一个信号是由时域和频域刻画的。
\qquad 设$g(t)=g(-t)$是一个实的对称窗函数，对它平移$u$各单位，以及调制（7-1小波分析之傅里叶分析有介绍平移与调制等性质）得到$$g_{u,\xi }(t)=e^{(i\xi t)}g(t-u)$$将$g$规范化为$\Vert g\Vert =1$，使对任意的$(u,\xi )\in {\mathbb{R}}^2$有$\Vert g_{u,\xi }\Vert =1。f\in L^2(\mathbb{R})$的窗口傅里叶变换为$$Sf(u,\xi )=<f,g_{(u,\xi )}>={\int }_{-\infty }^{+\infty }f(t)g(t-u)e^{-i\xi t}dt.$$

2.Heisenberg盒

\qquad 由于$g$是偶函数，故$g_{u,\xi }(t)=e^{(i\xi t)}g(t-u)$以$u$为中心，且它以$u$为中心的时间跨度不依赖于$u$和$\xi$，即$${\sigma }_t^2={\int }_{-\infty }^{+\infty }(t-u{)}^2|g_{u,\xi }(t){|}^{2}dt={\int }_{-\infty }^{+\infty }{t}^2|g(t){|}^{2}dt$$(注：$|e^{i\xi t}|=1$),因为$g$实对称，故其傅里叶变换也实对称。$g_{u,\xi }$的傅里叶变换为$${\hat{g}}_{u,\xi }(\omega )=\hat{g}(\omega -\xi )e^{-iu(\omega -\xi )}$$(注：傅里叶变换的调制与平移也参照第一讲内容)它是频率窗$\hat{g}$的$\xi$个平移，所以其中心频率为$\xi$.它以$\xi$为中心的频率跨度为$${\sigma }_{\omega }^2=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }(\omega -\xi {)}^2|{\hat{g}}_{u,\xi }(\omega ){|}^{2}d\omega ={\int }_{-\infty }^{+\infty }{\omega }^2|\hat{g}(\omega ){|}^{2}d\omega$$.它不依赖于$u,\xi$,所以$g_{u,\xi }$对应以$(u,\xi )$为中心，面积为${\sigma }_t^2{\sigma }_{\omega }^2$的Heisenberg盒。它的大小不依赖于$(u,\xi )$所以意味着窗口傅里叶变换在时频平面上有不变的分辨率。

3.信号的重构

\qquad 当时频指标$(u,\xi )$在${\mathbb{R}}^2$上移动时，原子$g_{u,\xi }$的Heisenberg盒覆盖了整个时-频空间。可以预见，$f$可由其窗口傅里叶变换$Sf(u,\xi )$来恢复。
窗口傅里叶变换重构公式： 若$f\in L^2(\mathbb{R})$,则$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }Sf(u,\xi )g(t-u)e^{i\xi t}d\xi du$$且$${\int }_{-\infty }^{+\infty }|f(t){|}^{2}dt=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }|Sf(u,\xi ){|}^{2}d\xi du$$
证明：令$f_{\xi }=Sf(u,\xi )$,应用傅里叶逆变换和Passval等式可证。
$Sf(u,\xi )={\int }_{-\infty }^{+\infty }f(t)g(t-u)e^{-i\xi t}dt=e^{-i\xi u}{\int }_{-\infty }^{+\infty }f(t)g(t-u)e^{i\xi (u-t)}dt=e^{-i\xi u}f\ast g_{\xi }(u)$,其中$g(t)=g(-t),g_{\xi }(u)=g(u)e^{-i\xi u}$，故$f_{\xi }$的傅里叶变换为$${\hat{f}}_{\xi }(\omega )=\hat{f}(\omega +\xi ){\hat{g}}_{\xi }(\omega +\xi )=\hat{f}(\omega +\xi )\hat{g}(\omega )$$,第一个等式的原因参照7-1节性质（卷积的傅里叶变换以及调制的傅里叶变换）。$g(t-u)$关于$u$的傅里叶变换为$\hat{g}(u)e^{-i\xi t}$,于是$$\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }Sf(u,\xi )g(t-u)e^{i\xi t}d\xi du=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\frac{1}{2\pi }{\hat{f}}_{\xi }(\omega )\ast \hat{g}(\omega )e^{-it\xi }{e}^{i\xi t}d\omega \text{\hspace{1em}(1)}$$,由Fubini定理，因为$f\in L^2$,所以$Sf\in L^2$,由$g$的定义可知$g\in L^2$,于是它们乘积的傅里叶变换也属于$L^2$,所以上式可积分换续，可得$$(1)=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\hat{f}(\omega )\hat{g}(\omega )e^{i\xi t}d\xi d\omega =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\hat{f}(\xi )e^{i\xi t}(\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{g}(\omega ){|}^{2}d\omega )d\xi =f(t)$$
定理的第二等式可由Plancherel公式导出：$$\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }|Sf(u,\xi ){|}^{2}d\xi du=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{S}f(\omega ,\xi ){|}^{2}d\xi d\omega =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{f}(\omega +\xi )\hat{g}(\omega ){|}^{2}d\xi d\omega$$再由Fubini定理可得$$\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }|Sf(u,\xi ){|}^{2}d\xi du=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }(\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{f}(\omega +\xi ){|}^{2}d\xi )|\hat{g}(\omega ){|}^{2}d\omega =\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\Vert f\Vert |\hat{g}(\omega ){|}^{2}d\omega =\Vert f\Vert$$