小波函数与逆函数,及其离散形式表达
由于小波变换作为一种可以将原始信号分解为不同频带的数学工具,而且在上一章中我们已经简要的介绍了它与傅里叶函数的异同,和一些其他特点,所以在这一章中我们再继续探讨一些关于小波在信号分解和合成方面的数学原理。
首先,还是我们的小波函数,它长下面这个样子:
W a b x ( t ) = 1 a ∫ − ∞ ∞ x ( τ ) ψ ( τ − b a ) d τ W_a^bx(t)=\frac{1}{\sqrt{a}}\int_{-\infty}^{\infty}x(\tau)\psi\left(\frac{\tau-b}{a}\right)d\tau Wabx(t)=a1∫−∞∞x(τ)ψ(aτ−b)dτ
其中, ψ ( τ − b a ) \psi \left(\frac{\tau-b}{a}\right) ψ(aτ−b) 是小波函数,比如墨西哥帽小波,Haar小波等,这里的 a a a 和 b b b 是尺度和平移参数,又或者你也能在某些教科书中见到这样的写法 ψ ( t ) \psi(t) ψ(t) 。而 x ( τ ) x(\tau) x(τ) 是原始信号。
对于上式,其对应的逆变换为:
x ( t ) = 1 C ψ ∫ 0 ∞ ∫ − ∞ ∞ W a b x ( t ) ψ a , b ( t ) d a d b x(t)=\frac{1}{C_\psi}\int_{0}^{\infty}\int_{-\infty}^{\infty}W_a^bx(t)\psi_{a,b}(t)dadb x(t)=Cψ1∫0∞∫−∞∞Wabx(t)ψa,b(t)dadb
其中, ψ a , b ( t ) = 1 a ψ ( t − b a ) \psi_{a,b}(t)=\frac{1}{\sqrt{a}}\psi\left(\frac{t-b}{a}\right) ψa,b(t)=a1ψ(at−b), C ψ C_\psi Cψ是小波函数的常数。这个公式表示说,执行连续小波函数的逆运算,其实等价于小波函数与所得系数的连续卷积,即可以简要的表达为:
x ( t ) = ∫ − ∞ ∞ W ( t − τ ) ∗ ψ ( τ ) d τ x(t) = \int_{-\infty}^{\infty} W(t - \tau) * \psi (\tau) \ d\tau x(t)=∫−∞∞W(t−τ)∗ψ(τ) dτ
这里的小波函数 ψ ( τ ) \psi(\tau) ψ(τ)是一个满足正交条件的函数,即:
∫ − ∞ ∞ ψ ( t ) ψ ( t − τ ) d t = δ ( τ ) \int_{-\infty}^{\infty}\psi(t)\psi(t-\tau)dt=\delta(\tau) ∫−∞∞ψ(t)ψ(t−τ)dt=δ(τ)
其中, δ ( τ ) \delta(\tau) δ(τ)是单位冲击函数。
所以,我们可以立即得到一个和上面相似的,但是更简单点的方法来描述这组公式,也就是小波的离散形式,即对于离散信号 x [ n ] x[n] x[n],其小波公式为:
c i [ n ] = ∑ k x [ k ] ∗ h i [ n − k ] c_i[n] = \sum_{k} x[k] * h_{i}[n - k] ci[n]=k∑x[k]∗hi[n−k]
其中 c i [ n ] c_i[n] ci[n] 为小波系数, h i [ n ] h_i[n] hi[n] 为小波函数(Wavelet Function), ∗ * ∗ 表示卷积运算。而相对的,如果我们获得了一组小波系数 c i [ n ] c_i[n] ci[n],那么逆小波变换的公式为:
x [ n ] = ∑ i = 1 N h i [ n ] ∗ c i [ n ] x[n] = \sum_{i=1}^{N} h_{i}[n] * c_{i}[n] x[n]=i=1∑Nhi[n]∗ci[n]
其中 h i [ n ] h_i[n] hi[n] 为小波函数(Wavelet Function), ∗ * ∗ 表示卷积运算。
常见小波核函数
Haar 小波
Haar小波是最简单的小波,其公式为
h ( t ) = { 1 t ∈ [ 0 , 1 2 ] − 1 t ∈ [ 1 2 , 1 ] 0 otherwise h(t) = \begin{cases} 1 & t \in [0, \frac{1}{2}] \\ -1 & t \in [\frac{1}{2}, 1] \\ 0 & \text{otherwise} \end{cases} h(t)=⎩ ⎨ ⎧1−10t∈[0,21]t∈[21,1]otherwise
Haar 小波是最简单的小波核函数,具有较低的计算复杂度。它可以用于简单的图像压缩和去噪。
Daubechies 小波
Daubechies 小波是一类常用的小波核函数,其中最常用的是 D 4 D4 D4 小波,公式为:
h ( t ) = 1 + 3 4 ( 1 + t ) e − t 2 u ( t ) − 1 − 3 4 ( 1 − t ) e − t 2 u ( − t ) h(t) = \frac{1 + \sqrt{3}}{4} (1 + t) e^{-\frac{t}{2}} u(t) - \frac{1 - \sqrt{3}}{4} (1 - t) e^{-\frac{t}{2}} u(-t) h(t)=41+3(1+t)e−2tu(t)−41−3(1−t)e−2tu(−t)
其中 u ( t ) u(t) u(t) 为单位函数。Daubechies 小波是一类常用的小波核函数,具有较好的时频分解性能。它常用于图像压缩、信号处理、医学图像处理等应用中。
Symlets 小波
Symlets 小波是一类常用的小波核函数,其中最常用的是 S 8 S8 S8 小波,公式为:
h ( t ) = − 1 8 t 3 e − t 2 u ( t ) + 1 4 t e − t 2 u ( − t ) h(t) = -\frac{1}{8} t^3 e^{-\frac{t}{2}} u(t) + \frac{1}{4} t e^{-\frac{t}{2}} u(-t) h(t)=−81t3e−2tu(t)+41te−2tu(−t)
Symlets 小波是一类常用的小波核函数,具有较高的计算效率。它常用于图像压缩、信号处理、音频处理等应用中。
Coiflets 小波
Coiflets 小波是一类常用的小波核函数,其中最常用的是 C 6 C6 C6 小波,公式为:
h ( t ) = 1 + 6 32 t 5 e − t 2 u ( t ) − 1 + 6 8 t 3 e − t 2 u ( − t ) + 1 − 6 32 t e − t 2 u ( t ) h(t) = \frac{1 + \sqrt{6}}{32} t^5 e^{-\frac{t}{2}} u(t) - \frac{1 + \sqrt{6}}{8} t^3 e^{-\frac{t}{2}} u(-t) + \frac{1 - \sqrt{6}}{32} t e^{-\frac{t}{2}} u(t) h(t)=321+6t5e−2tu(t)−81+6t3e−2tu(−t)+321−6te−2tu(t)
Coiflets 小波是一类常用的小波核函数,具有较高的时域精度和较低的频域精度。它常用于图像压缩、信号处理、医学图像处理等应用中。
Marr 小波
h ( t ) = 1 2 [ 1 2 e − t 2 u ( t ) − 1 2 e t 2 u ( − t ) ] h(t) = \frac{1}{\sqrt{2}} \left[ \frac{1}{\sqrt{2}} e^{-\frac{t}{2}} u(t) - \frac{1}{\sqrt{2}} e^{\frac{t}{2}} u(-t) \right] h(t)=21[21e−2tu(t)−21e2tu(−t)]
其中 u ( t ) u(t) u(t) 为单位函数。Marr wavelet 具有窄带频率响应,可以对图像进行高通滤波。它在图像增强中常用于去除噪声或提取图像的纹理信息。
小波函数的关键特性
我们不能无脑的使用小波去进行信号或图像处理,所以必须对小波的重要特性做一定的了解。一般说来,小波具备如下关键特性:
- 小波具有时频均衡的特性,即在时域和频域上的分解精度相似。
- 小波具有多阶平移不变的特性,即平移小波函数的位置,可以得到完全相同的分解结果。
- 小波具有稳定性的特性,即经过小波变换后的信号对于原信号的变化非常敏感。
- 小波具有稀疏性的特性,即大多数信号在小波变换后的系数都非常小,只有很少的系数具有较大的值。
在下一章节中,我们将看看一些具体的应用,和如何人工实现一个小波用于信号处理的。