7-1 小波分析之傅里叶分析

7-1 小波分析之傅里叶分析

1.线性时不变滤波

\qquad 一个算子$L$的时移不变性是指：如果输入的信号$f(t)$延迟$\tau$个单位成为$f_{\tau }(t)=f(t-\tau )$则输出信号也延时$\tau$个单位，即$$g(t)=Lf(t)⟹g(t-\tau )=L{f}_{\tau }(t)$$

2.Fourier积分

2.1$L^1(R)$上的傅里叶变换

• $L^1(R)$上的傅里叶变换
  \qquad 傅里叶积分定义为$$\hat{f}(\omega )={\int }_{-\infty }^{\infty }f(t)e^{-i\omega t}dt$$它的意义是度量$f$在频率$\omega$处的震荡信息。
• $L^1(R)$上的傅里叶逆变换
  \qquad 如果$f\in L^1(R)$,$\hat{f}\in L^1(R)$则$$f(t)=\frac{1}{2\pi }{\int }_{-\infty }^{\infty }\hat{f}(\omega )e^{i\omega t}d\omega$$
• 卷积
  $$f\ast g(t)={\int }_{-\infty }^{+\infty }f(u)g(t-u)du$$
  性质：

1. 可交换性：$f\ast g=g\ast f$
2. 微分性：$\frac{d}{dt}(f\ast g)(t)=\frac{df}{dt}\ast g(t)=f(t)\ast \frac{dg}{dt}$
3. Dirac卷积：$f\ast {\delta }_{\tau }(t)=f(t-\tau )$

• Fourier变换与卷积
  \qquad 若$f\in L^1(R),h\in L^1(R)$,则$g=h\ast f\in L^1(R)$且$$\hat{g}(\omega )=\hat{h}(\omega )\hat{f}(\omega )$$.
  证明：$\hat{g}(\omega )={\int }_{-\infty }^{+\infty }g(t)e^{-i\omega t}dt={\int }_{-\infty }^{+\infty }({\int }_{-\infty }^{+\infty }f(t-u)h(u)du)e^{-i\omega t}dt$又因为$f,h$均可积，故积分可换续，令$u=u,v=t-u$得${\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }{e}^{-i(u+v)\omega }f(v)h(u)dudv=({\int }_{-\infty }^{\infty }f(v)e^{-i\omega v}dv)({\int }_{-\infty }^{\infty }h(u)e^{-i\omega u}du)$
  证毕
• Fourier变换的性质
  

2.2 $L^2(R)$上的傅里叶变换

• $L^2(R)$上特征函数的傅里叶变换
  $[-1,1]$上的特征函数的傅里叶变换为$$\hat{f}(\omega )={\int }_{-1}^1{e}^{-i\omega t}dt=\frac{2sin\omega }{\omega }$$
  证明：$${\int }_{-\infty }^{+\infty }{e}^{-i\omega t}dt=-\frac{1}{i\omega }{e}^{-i\omega t}{|}_{-1}^1=-\frac{1}{i\omega }{e}^{-i\omega }+\frac{1}{i\omega }{e}^{i\omega }=\frac{2sin\omega }{\omega }$$证毕
• $Parsaeval$等式及$Plancherel$等式
  \qquad 若$f,h\in L_1(R)\cap L_2(R)$，则$${\int }_{-\infty }^{+\infty }f(t)g^{\ast }(t)dt=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }\hat{f}(\omega ){\hat{g}}^{\ast }(\omega )d\omega$$,对于$h=f$则有$${\int }_{-\infty }^{+\infty }|f(t){|}^{2}dt=\frac{1}{2\pi }{\int }_{-\infty }^{+\infty }|\hat{f}(\omega ){|}^{2}d\omega$$
• 脉冲响应
  \qquad 理想的低通滤波器有传递函数$\hat{\varphi }=1_{[-\xi ,\xi ]}$,它截取低频$[-\xi ,\xi ]$的信息。由逆傅里叶积分可得$${\varphi }_{\xi }(t)=\frac{1}{2\pi }{\int }_{-\xi }^{+\xi }{e}^{-i\omega t}dt=\frac{sin(\xi t)}{\pi t}\text{\hspace{1em}(1)}$$
• 高斯函数的傅里叶变换
  \qquad 高斯函数的傅里叶变换也是高斯函数，$f(t)=e^{-t^2},\hat{f}(\omega )=\sqrt{\pi }{e}^{-\frac{{\omega }^2}{4}}$

2.3 性质

• 正则性与衰减性
  \qquad 如果$f$满足$${\int }_{-\infty }^{+\infty }|\hat{f}(\omega )|(1+|\omega {|}^p)d\omega <+\infty$$则函数$f$是有界函数，且$p$次连续可微、导数有界。
  \qquad 这个结果也证明，若存在常数$K$及$ϵ>0$使得$$|\hat{f}(\omega )|<\frac{K}{1+|\omega {|}^{p+1+ϵ}}$$则$f\in C^p$.
• 测不准原理
  \qquad 我们知道只有当$f$在时域上变化缓慢时，$|\hat{f}(\omega )|$才会在高频快速衰减。因此当$\hat{f}$的能量分散在所有频率上时$f$也必须分散在相对宽的时间域上。
  Heisenberg测不准原理：$f\in L^2(R)$的时频变差满足如下不等式$${\sigma }_t^2{\sigma }_{\omega }^2\ge \frac{1}{4}$$其中等号成立当且仅当存在$(u,\xi ,a,b)\in R^2\times C^2$,使得$$f(t)=a{e}^{i\xi t-b(t-u{)}^2}$$.
  注：${\sigma }_t^2=\frac{1}{||f|{|}^2}{\int }_{-\infty }^{+\infty }(t-u{)}^2|f(t){|}^{2}dt$,${\sigma }_{\omega }^2=\frac{1}{2\pi ||f|{|}^2}{\int }_{-\infty }^{+\infty }(\omega -\xi {)}^2|\hat{f}(\omega ){|}^{2}d\omega$
  证明：
  
  
• Gibbs震荡
  \qquad 对任意$\xi >0,$有$$u\ast {\varphi }_{\xi }(t)={\int }_{-\infty }^{\xi t}\frac{sinx}{\pi x}dx$$
  证明： 由$(1)$可得理想低通脉冲响应为${\varphi }_{\xi }(t)=\frac{sin\xi t}{\pi t}$，因此$$\begin{gathered} u\ast {\varphi }_{\xi }(t)={\int }_{-\infty }^{+\infty }u(\tau )\frac{sin\xi (t-\tau )}{ \\ \pi (t-\tau )}d\tau ={\int }_0^{+\infty }\frac{sin\xi (t-\tau )}{ \\ \pi (t-\tau )}d\tau \text{\hspace{1em}(2)} \end{gathered}$$
  其中$$u(t)=\{\begin{array}{ll}1,& \text{if t≥0}\\ 0,& \text{else}\end{array}$$
  令$x=\xi (t-\tau )$则$$u\ast {\varphi }_{\xi }(t)={\int }_{-\infty }^{\xi t}\frac{\sin x}{\pi x}dx$$证毕.
• 图像的全变差

1. 全变差：若$f$可微，它的全变差定义为$$\Vert f{\Vert }_V={\int }_{-\infty }^{+\infty }|f^{\prime }(t)|dt$$如果 { x p } p \{x_p\}_p {xp​}p​满足$f^{\prime }(x_p)=0$即 { x p } p \{x_p\}_p {xp​}p​为$f$的局部极值点集合，则有$$\Vert f{\Vert }_V=\sum _p|(f(x_{p+1})-f(x_p))|$$它测量了$f$的整体震荡幅度。
2. 图像全变差：假设$f(x_1,x_2)$可微，它的全变为$$\Vert f{\Vert }_V=\int \int |\overrightarrow{\nabla }f(x_1,x_2)|d{x}_{1}d{x}_2$$此处梯度向量的模定义为$$|\overrightarrow{\nabla }f(x_1,x_2)|=(|\frac{\partial f(x_1,x_2)}{\partial x_1}{|}^2+|\frac{\partial f(x_1,x_2)}{\partial x_2}{|}^2{)}^{\frac{1}{2}}$$.

2.4 二维傅里叶变换

$R^n$上的傅里叶变换是一维情况下的推广。为了图像处理的应用需要，我们做简单的回顾：二维傅里叶变换的积分为$$\hat{f}({\omega }_1,{\omega }_2)={\int }_{-\infty }^{+\infty }{\int }_{-\infty }^{+\infty }f(x_1,x_2)e^{-i{\omega }_1{x}_1-i{\omega }_2{x}_2}d{x}_{1}d{x}_2$$在极坐标下，$e^{i{\omega }_1{x}_1+i{\omega }_2{x}_2}$可表示为$e^{i\xi (x_{1}cos\theta +x_{2}sin\theta )}$其中$\xi =\sqrt{{\omega }_1^2+{\omega }_2^2}$.可见$e^{i\xi (x_{1}cos\theta +x_{2}sin\theta )}$是沿着$\theta$传播、以$\xi$频率震荡的平面波。

• 二维性质
  令$x=(x_1,x_2),\omega =({\omega }_1,{\omega }_2),f(x)=f(x_1,x_2)$

1. 若$f,\hat{f}\in L^1(R^2)$则$$f(x)=\frac{1}{4{\pi }^2}\int \int \hat{f}(\omega )e^{-i\omega \cdot x}d\omega$$.
2. 卷积：$f,h\in L^1(R^2)$则$$g(x)=f\ast h(x)={\int }_{-\infty }^{+\infty }f(u)\cdot g(x-u)du$$的傅里叶变换满足$$\hat{g}(\omega )=\hat{f}(\omega )\hat{h}(\omega )$$
3. parseval
   $$\int \int f(x)g^{\ast }(x)dx=\frac{1}{4{\pi }^2}\int \int \hat{f}(\omega ){\hat{g}}^{\ast }(\omega )d\omega$$若$f=g$则上式为：$$\int \int |f(x){|}^{2}dx=\frac{1}{4{\pi }^2}\int \int |\hat{f}(\omega ){|}^{2}d\omega$$
4. 如果$f(x)$旋转$\theta$,$$f_{\theta }(x_1,x_2)=f(x_{1}cos\theta -x_{2}sin\theta ,x_{1}sin\theta +x_{2}cos\theta )$$则$f_{\theta }$的傅里叶变换也旋转$\theta$,$${\hat{f}}_{\theta }(\omega )=\hat{f}({\omega }_{1}cos\theta -{\omega }_{2}sin\theta ,{\omega }_{1}sin\theta +{\omega }_{2}cos\theta )$$
   .