对具有频率间断点的信号使用五阶消失矩的symlet 进行连续小波变换 在数学中，连续小波变换 （Continuous Wavelet Transform, CWT）是一种时频分析工具，通过让小波函数（Wavelet）的平移参数和尺度参数连续变化，提供信号的过完备表示。
函数
x
(
t
)
{\displaystyle x(t)}
a
∈
R
+
∗
{\displaystyle a\in \mathbb {R^{+*}} }
b
∈
R
{\displaystyle b\in \mathbb {R} }
X
w
(
a
,
b
)
=
1
|
a
|
1
/
2
∫
−
∞
∞
x
(
t
)
ψ
¯
(
t
−
b
a
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{|a|^{1/2}}}\int _{-\infty }^{\infty }x(t){\overline {\psi }}\left({\frac {t-b}{a}}\right)\,dt}
其中
ψ
(
t
)
{\displaystyle \psi (t)}
称为母小波，是在时间和频率域中的连续函数，
⋅
¯
{\displaystyle {\bar {\cdot }}}
表示
复共轭 。母小波的主要目的是为生成子小波（即母小波的尺度伸缩和平移）提供原函数。逆连续小波变换（Inverse continuous wavelet transform）可以用于恢复信号
x
(
t
)
{\displaystyle x(t)}
:
x
(
t
)
=
C
ψ
−
1
∫
0
∞
∫
−
∞
∞
X
w
(
a
,
b
)
1
|
a
|
1
/
2
ψ
~
(
t
−
b
a
)
d
b
d
a
a
2
{\displaystyle x(t)=C_{\psi }^{-1}\int _{0}^{\infty }\int _{-\infty }^{\infty }X_{w}(a,b){\frac {1}{|a|^{1/2}}}{\tilde {\psi }}\left({\frac {t-b}{a}}\right)\,db\ {\frac {da}{a^{2}}}}
ψ
~
(
t
)
{\displaystyle {\tilde {\psi }}(t)}
ψ
(
t
)
{\displaystyle \psi (t)}
对偶函数 ，且：
C
ψ
=
∫
−
∞
∞
ψ
^
¯
(
ω
)
ψ
~
^
(
ω
)
|
ω
|
d
ω
{\displaystyle C_{\psi }=\int _{-\infty }^{\infty }{\frac {{\overline {\hat {\psi }}}(\omega ){\hat {\tilde {\psi }}}(\omega )}{|\omega |}}\,d\omega }
为容许性常数，式中
⋅
^
{\displaystyle {\hat {\cdot }}}
ψ
~
(
t
)
=
ψ
(
t
)
{\displaystyle {\tilde {\psi }}(t)=\psi (t)}
C
ψ
=
∫
−
∞
+
∞
|
ψ
^
(
ω
)
|
2
|
ω
|
d
ω
{\displaystyle C_{\psi }=\int _{-\infty }^{+\infty }{\frac {\left|{\hat {\psi }}(\omega )\right|^{2}}{\left|\omega \right|}}\,d\omega }
一般来说，该常数为小波的容许性常数（Admissible constant）。满足：
0
<
C
ψ
<
∞
{\displaystyle 0<C_{\psi }<\infty }
的小波称为容许小波（Admissible wavelet）。容许小波满足
ψ
^
(
0
)
=
0
{\displaystyle {\hat {\psi }}(0)=0}
x
(
t
)
{\displaystyle x(t)}
x
(
t
)
=
1
2
π
ψ
^
¯
(
1
)
∫
0
∞
∫
−
∞
∞
1
a
2
X
w
(
a
,
b
)
exp
(
i
t
−
b
a
)
d
b
d
a
{\displaystyle x(t)={\frac {1}{2\pi {\overline {\hat {\psi }}}(1)}}\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{a^{2}}}X_{w}(a,b)\exp \left(i{\frac {t-b}{a}}\right)\,db\ da}
该式说明，小波可以定义为：
ψ
(
t
)
=
w
(
t
)
exp
(
i
t
)
{\displaystyle \psi (t)=w(t)\exp(it)}
其中
w
(
t
)
{\displaystyle w(t)}
小波转换(Wavelet Transform)的类型 [ 编辑 ] 小波转换 ( Wavelet Transform) 可依照输入与输出为连续或是离散(discrete)分成三种类型，
第一种，输入为连续，输出为连续，则称之为连续小波转换(Continuous Wavelet Transform)
第二种，输入为连续，输出为离散，则称之为连续离散系数小波转换(Continuous wavelet transform with discrete coefficients)
第三种，输入为离散，输出为离散，则称之为离散小波转换(Discrete Wavelet Transform)
并没有第四种，输入为离散输出为连续的小波转换，在应用中并不会将简单的讯号转换成更复杂的讯号 傅里叶转换 (Fourier Transform)与小波转换 比较共有四种类型
第一种，输入为连续，输出为连续，傅里叶转换(Fourier Transform)
第二种，输入为连续，输出为离散，傅里叶级数(Fourier Series)
第三种，输入为离散，输出为离散，离散傅里叶转换(Discrete Fourier Transform)
第四种，输入为离散，输出为连续，离散(时间)傅里叶转换(Discrete-time Fourier Transform) 连续小波转换介绍 [ 编辑 ] 连续小波转换 (Continuous Wavelet Transform)是一种用来分解一个连续时间函数，使它变成数个小波 (wavelet)。跟傅里叶变换 (Fourier Transform)不一样的是，连续小波转换可以建构一个具有良好时域 和频域 局部化的时频讯号。以数学来说，一个有连续时间性质且可积分的函数
x
(
t
)
{\displaystyle x(t)}
X
w
(
a
,
b
)
=
1
|
(
b
)
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {|(b)|}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt}
ψ
(
t
)
{\displaystyle \psi (t)}
a
{\displaystyle a}
b
{\displaystyle b}
a
{\displaystyle a}
(
−
∞
,
+
∞
)
{\displaystyle \left(-\infty ,+\infty \right)}
b
{\displaystyle b}
(
0
,
+
∞
)
{\displaystyle \left(0,+\infty \right)}
以时频分析的角度分析，当
b
{\displaystyle b}
以时频分析的角度分析，当
b
{\displaystyle b}
小波母函数的用途在于提供一个可以产生子波(Daughter Wavelet)的根源函数，而子波是小波母函数平移过或缩放过(或两者都有)的版本。如果要将已知且存在的讯号
x
(
t
)
{\displaystyle x(t)}
x
(
t
)
=
∫
0
∞
∫
−
∞
∞
1
b
2
X
w
(
a
,
b
)
1
|
(
b
)
|
ψ
1
(
t
−
a
b
)
d
a
d
b
{\displaystyle x(t)=\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{b^{2}}}X_{w}(a,b){\frac {1}{\sqrt {|(b)|}}}\psi _{1}({\frac {t-a}{b}})\,dadb}
ψ
1
(
t
)
{\displaystyle \psi _{1}(t)}
ψ
(
t
)
{\displaystyle \psi (t)}
∫
0
∞
∫
−
∞
∞
1
|
b
3
|
ψ
(
t
1
−
a
b
)
ψ
1
(
t
−
a
b
)
d
a
d
b
=
δ
(
t
−
t
1
)
{\displaystyle \int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{|b^{3}|}}\psi ({\frac {t_{1}-a}{b}})\psi _{1}({\frac {t-a}{b}})\,dadb={\boldsymbol {\delta }}(t-t_{1})}
有时
ψ
1
(
t
)
=
C
−
1
ψ
(
t
)
{\displaystyle \psi _{1}(t)=C^{-1}\psi (t)}
连续小波转换的问题 [ 编辑 ] 在实务上，相较于连续小波转换，通常较多会使用离散小波变换 、连续离散系数小波转换
连续小波转换组成要素 [ 编辑 ] 小波母函数(Mother Wavelet) [ 编辑 ] 举例，2个例子来说明小波母函数(Mother Wavelet):
Mother wavelet 的五大限制 [ 编辑 ] 紧支撑
支撑:函数的范围没有收敛区间
紧支撑: 函数的范围有收敛
实函数
奇对称或偶对称
消失矩
容许性条件
通常来说，我们会倾向选一个可以连续微分的小波母函数且拥有紧凑支撑(Compact Support)的尺度函数(Scaling Function)和高阶的消失矩(Vanishing Moment)。一个小波母函数是以这两个函数所组成:小波函数
ψ
(
t
)
{\displaystyle \psi (t)}
φ
(
t
)
{\displaystyle \varphi (t)}
m
k
=
∫
t
k
ψ
(
t
)
d
t
{\displaystyle m_{k}=\int t^{k}\psi (t)\,dt}
如果
m
0
=
m
1
=
m
2
=
.
.
.
.
.
=
m
p
−
1
=
0
{\displaystyle m_{0}=m_{1}=m_{2}=.....=m_{p-1}=0}
ψ
(
t
)
{\displaystyle \psi (t)}
p
{\displaystyle p}
a
−
i
{\displaystyle a^{-i}}
a
L
{\displaystyle a^{L}}
L
{\displaystyle L}
∫
t
k
ψ
(
t
)
d
t
=
∫
Ψ
∗
(
f
)
(
j
2
π
)
2
d
k
d
f
k
δ
(
f
)
d
f
d
t
{\displaystyle \int t^{k}\psi (t)\,dt=\int {\boldsymbol {\Psi ^{*}(f)}}({\frac {j}{2\pi }})^{2}{\frac {d^{k}}{df^{k}}}\delta (f)df\,dt}
⇒
∫
t
k
ψ
(
t
)
d
t
=
0
⟺
(
j
2
π
)
2
d
k
d
f
k
Ψ
∗
(
f
)
|
f
=
0
=
0
⟺
d
k
d
f
k
Ψ
∗
(
f
)
|
f
=
0
=
0
⟺
d
k
d
f
k
Ψ
(
f
)
|
f
=
0
=
0
{\displaystyle \Rightarrow \int t^{k}\psi (t)\,dt=0\iff ({\frac {j}{2\pi }})^{2}{\frac {d^{k}}{df^{k}}}{\boldsymbol {\Psi ^{*}(f)}}|_{f=0}=0\iff {\frac {d^{k}}{df^{k}}}{\boldsymbol {\Psi ^{*}(f)}}|_{f=0}=0\iff {\frac {d^{k}}{df^{k}}}{\boldsymbol {\Psi (f)}}|_{f=0}=0}
其中
Ψ
(
f
)
{\displaystyle {\boldsymbol {\Psi (f)}}}
ψ
(
t
)
{\displaystyle \psi (t)}
尺度函数 [ 编辑 ] 目的是简化反小波转换(Inverse wavelet transform) 过程 在连续小波转换中，定义尺度函数
φ
(
t
)
{\displaystyle \varphi (t)}
φ
(
t
)
=
∫
−
∞
∞
Φ
(
f
)
e
j
2
π
f
t
d
f
{\displaystyle \varphi (t)=\int _{-\infty }^{\infty }{\boldsymbol {\Phi (f)}}e^{j2\pi ft}df}
其中
Φ
(
f
)
{\displaystyle {\boldsymbol {\Phi (f)}}}
φ
(
t
)
{\displaystyle \varphi (t)}
|
Φ
(
f
)
|
2
=
∫
f
∞
|
Ψ
(
f
1
)
|
2
|
f
1
|
2
d
f
1
{\displaystyle {\left|{\boldsymbol {\Phi (f)}}\right|}^{2}=\int _{f}^{\infty }{\frac {{\left|{\boldsymbol {\Psi (f_{1})}}\right|}^{2}}{{\left|f_{1}\right|}^{2}}}d{f_{1}}}
小波函数
ϕ
(
t
)
{\displaystyle \phi (t)}
φ
(
t
)
{\displaystyle \varphi (t)}
缩放因子 [ 编辑 ] 缩放因子可以压缩或拉长一个讯号。当缩放因子的值相对低时，讯号会比较紧缩，也就是会造成一个更细致的图像。但是低缩放因子的缺点是它的效果无法覆盖一个讯号的持续期间。另一方面，当缩放因子的值相对高时，讯号会比较被拉长造成一个比较粗糙的图像，但是它的效果会持续整个讯号的期间。
连续小波转换性质 [ 编辑 ] 1. 输入与输出的性质关系
基本上，连续小波转换是输入资料序列和一组由小波母函数所产生函数的卷积 。这个卷积可以用快速傅里叶变换来计算。除非小波母函数为虚数函数，在正常的情况下，输出信息
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
‖
X
w
(
a
,
b
)
|
2
{\displaystyle \|X_{w}(a,b)|^{2}}
2. 时间轴上的位移
输入函数
x
(
t
)
{\displaystyle x(t)}
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
若
x
(
t
)
{\displaystyle x(t)}
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
则
x
(
t
−
τ
)
{\displaystyle x(t-\tau )}
X
w
(
a
−
τ
,
b
)
{\displaystyle X_{w}(a-\tau ,b)}
3. 时间轴上的缩放
当进行小波转换的函数在时间轴上拉长或压缩时，输出函数也会有相对应的变化： 若
x
(
t
)
{\displaystyle x(t)}
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
则
x
(
t
/
σ
)
{\displaystyle x(t/\sigma )}
σ
X
w
(
a
/
σ
,
b
/
σ
)
{\displaystyle {\sqrt {\sigma }}X_{w}(a/\sigma ,b/\sigma )}
4. 帕瑟伐定理(Parseval's Theory)
在
φ
(
t
)
=
C
−
1
ψ
(
t
)
{\displaystyle \varphi (t)=C^{-}1\psi (t)}
∫
|
x
(
t
)
|
2
d
t
=
1
C
∫
0
∞
∫
−
∞
∞
1
b
2
|
X
w
(
a
,
b
)
|
2
d
a
d
b
{\displaystyle \int {|x(t)|}^{2}\,dt={\frac {1}{C}}\,\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{b^{2}}}{|X_{w}(a,b)|}^{2}\,da\,db}
相较于加伯变换 [ 编辑 ] 加伯变换 在处理讯号时不管是高频或是低频，尺度皆是相同的
而小波转换则会根据不同的频率改变其本身的的尺度
小波变换的分辨率在a-axis中，不因a值的改变而改变，但延著不同的b值改变以得到较好的结果。
修正型连续小波转换 [ 编辑 ] 利用上述尺度函数的定义，我们可以定义修正型的连续小波转换。将原本的连续小波转换定义为
X
w
(
a
,
b
)
=
1
b
′
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
′
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {b^{\prime }}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b^{\prime }}})\,dt}
在
a
{\displaystyle a}
0
<
b
′
<
b
0
{\displaystyle 0<b^{\prime }<b_{0}}
L
X
w
(
a
,
b
0
)
=
1
b
0
∫
−
∞
∞
x
(
t
)
φ
(
t
−
a
b
0
)
d
t
{\displaystyle {LX}_{w}(a,b_{0})={\frac {1}{\sqrt {b_{0}}}}\int _{-\infty }^{\infty }x(t)\varphi ({\frac {t-a}{b_{0}}})\,dt}
借由新定义的函数
L
X
w
(
a
,
b
0
)
{\displaystyle {LX}_{w}(a,b_{0})}
x
(
t
)
=
1
C
ψ
[
∫
0
b
0
∫
−
∞
∞
1
b
′
5
/
2
X
w
(
a
,
b
′
)
ψ
(
t
−
a
b
′
)
d
a
d
b
+
∫
−
∞
∞
1
b
0
3
/
2
L
X
w
(
a
,
b
0
)
φ
(
t
−
a
b
0
)
d
a
]
{\displaystyle x(t)={\frac {1}{C_{\psi }}}\left[\int _{0}^{b_{0}}\int _{-\infty }^{\infty }{\frac {1}{{b^{\prime }}^{5/2}}}X_{w}(a,b^{\prime })\psi ({\frac {t-a}{b^{\prime }}})\,da\,db+\int _{-\infty }^{\infty }{\frac {1}{{b_{0}}^{3/2}}}{LX}_{w}(a,b_{0})\varphi ({\frac {t-a}{b_{0}}})\,da\right]}
此建构
x
(
t
)
{\displaystyle x(t)}
C
ψ
=
∫
0
∞
|
Ψ
(
f
)
|
2
|
f
|
d
f
{\displaystyle C_{\psi }=\int _{0}^{\infty }{\frac {|{\Psi (f)|}^{2}}{|f|}}\,df}
此式中
Ψ
(
f
)
{\displaystyle \Psi (f)}
ψ
(
f
)
{\displaystyle \psi (f)}
通常在设计母小波函数时，会要求
C
ψ
<
∞
{\displaystyle C_{\psi }\ <\ \infty }
小波量值图(Scalogram) [ 编辑 ] 在许多文献资料中，常用小波量值图(Scalogram)来表示连续小波转换后的结果。其定义如下
S
c
x
(
a
,
b
)
=
|
X
w
(
a
,
b
)
|
2
=
1
|
b
|
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
|
2
{\displaystyle Sc_{x}(a,b)={\left|X_{w}(a,b)\right|}^{2}={\frac {1}{|b|}}{\left|\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt\right|}^{2}}
此处的定义即为连续小波转换结果的绝对值平方，用以视觉化连续小波转换的结果。
小波量值图之于小波转换函数的意义和频谱图 之于短时距时频平分析的意义相似。
在实际应用上通常以三个轴来显示，分别代表时间、频率与小波量值图的振幅。若是在二维图片则是利用颜色深浅来表示小波量值图的强度。
连续小波转换应用 [ 编辑 ] 小波转换最热门的一个应用为图像压缩。用小波转换式的编码在图像压缩可以提供显著的图像品质改善且给予更高的压缩比率。因为小波转换可以分解一个复杂的讯息或图案成基本型式，它在音乐档案和图型辨识上被广泛的使用。此外，我们还可以在以下的科学研究领域见到小波转换的应用: 边缘检测，解偏微分方程，脑电瞬态信号检测，滤波器设计，心电图分析，衣料分析和商业资讯分析。
离散变数连续小波转换 [ 编辑 ] 在离散变数连续小波转换(Continuous Wavelet Transform with Discrete Coefficients)中，原本
X
w
(
a
,
b
)
=
1
|
(
b
)
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {|(b)|}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt}
中的
a
,
b
{\displaystyle a,\,b}
a
=
n
2
−
m
,
b
=
2
−
m
{\displaystyle a=n2^{-m},\ b=2^{-m}}
则离散变数连续小波转换则重新表示为
X
w
(
a
,
b
)
=
2
m
/
2
∫
−
∞
∞
x
(
t
)
ψ
(
2
m
−
n
)
d
t
{\displaystyle X_{w}(a,b)=2^{m/2}\int _{-\infty }^{\infty }x(t)\psi (2^{m}-n)\,dt}
其中
n
∈
Z
,
n
∈
(
−
∞
,
∞
)
{\displaystyle n\in Z,\ n\in (-\infty ,\infty )}
n
∈
Z
,
n
∈
(
−
∞
,
∞
)
{\displaystyle n\in Z,\ n\in (-\infty ,\infty )}
选择离散变数连续小波转换的主要目的，在于简化连续小波转换在实作上的复杂性，并且利用快速算法增加应用价值。实际上离散变数连续小波转换算是连续小波转换的一种特例，部分文献将其当作离散小波转换讨论。
此处的选择
a
{\displaystyle a}
b
{\displaystyle b}
X
w
(
n
,
m
−
1
)
{\displaystyle X_{w}(n,m-1)}
X
w
(
n
,
m
)
{\displaystyle X_{w}(n,m)}
离散变数反转连续小波转换 [ 编辑 ] 在
a
=
n
2
−
m
,
b
=
2
−
m
{\displaystyle a=n2^{-m},\ b=2^{-m}}
x
(
t
)
=
∑
m
=
−
∞
∞
∑
n
=
−
∞
∞
2
m
/
2
ψ
1
(
2
m
−
n
)
X
w
(
n
,
m
)
{\displaystyle x(t)=\sum _{m=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }2^{m/2}\psi _{1}(2^{m}-n)X_{w}(n,m)}
ψ
1
(
t
)
{\displaystyle \psi _{1}(t)}
ψ
(
t
)
{\displaystyle \psi (t)}
∑
m
=
−
∞
∞
∑
n
=
−
∞
∞
2
m
ψ
1
(
2
m
−
n
)
ψ
(
2
m
t
1
−
n
)
=
δ
t
−
t
1
{\displaystyle \sum _{m=-\infty }^{\infty }\sum _{n=-\infty }^{\infty }2^{m}\psi _{1}(2^{m}-n)\psi (2^{m}t_{1}-n)=\delta {t-t_{1}}}
又此条件可表示为
∫
−
∞
∞
2
m
ψ
1
(
2
1
m
−
n
1
)
ψ
(
2
m
t
1
−
n
)
d
t
=
δ
m
−
m
1
δ
n
−
n
1
{\displaystyle \int _{-\infty }^{\infty }2^{m}\psi _{1}(2_{1}^{m}-n_{1})\psi (2^{m}t_{1}-n)\,dt=\delta {m-m_{1}}\delta {n-n_{1}}}
离散变数连续小波优缺点 [ 编辑 ] 优点:
快速算法
正交性质（Orthogonal ）
非均匀频率分析（Non-uniform frequency analysis） 缺点:
小波转换应用 [ 编辑 ] 小波转换的应用有以下两项特点:
信号的频率分布，会随着不同的时间 (或地点 )有较大变化
多尺度的分析扮演重要的角色
大采样间隔
⟶
{\displaystyle \longrightarrow }
小采样间隔
⟶
{\displaystyle \longrightarrow }
应用:
影像压缩，例如JPEG /JPEG2000
边缘角落侦测
特征辨识
强调前景压缩背景
滤波器设计
声音讯号
指纹辨识
金融
气象分析 相关条目 [ 编辑 ] 参考文献 [ 编辑 ] 
![{\displaystyle x(t)={\frac {1}{C_{\psi }}}\left[\int _{0}^{b_{0}}\int _{-\infty }^{\infty }{\frac {1}{{b^{\prime }}^{5/2}}}X_{w}(a,b^{\prime })\psi ({\frac {t-a}{b^{\prime }}})\,da\,db+\int _{-\infty }^{\infty }{\frac {1}{{b_{0}}^{3/2}}}{LX}_{w}(a,b_{0})\varphi ({\frac {t-a}{b_{0}}})\,da\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8456bf9f1862e769464cf6a666df895329471451)
