系列文章目录
提示:这里可以添加系列文章的所有文章的目录,目录需要自己手动添加
例如:第一章 Python 机器学习入门之pandas的使用
提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档
前言
由于傅里叶变换是一种全局的函数变换,并且损失了空间信息和时间信息。因此在后来引入了小波变换。小波变换的目的是将空域或时域的函数转换为由尺度函数和小波函数作为基的函数。换句话说,我们采用一系列的尺度函数和小波函数对原函数进行贴合近似。在相邻的尺度函数之间,我们用多个小波函数去近似原函数。
一、什么是小波函数?
顾名思义,小波函数就是在狭小区间里迅速上升又迅速下降的函数。在该区间 ( a , b ) (a,b) (a,b)里小波函数有 ∫ a b f ( x ) d x = 1 \int _a^b f(x)dx=1 ∫abf(x)dx=1的特点,也就是 X X X轴上下的面积相等。哈儿函数是最早的小波函数。
哈尔函数
u
=
2
p
+
q
u=2^p+q
u=2p+q
h
u
(
x
)
=
{
1
u
=
0
a
n
d
0
≤
x
<
1
2
p
/
2
u
>
0
a
n
d
q
/
2
p
≤
x
<
(
q
+
0.5
)
/
2
p
−
2
p
/
2
u
<
0
a
n
d
(
q
+
0.5
)
/
2
p
≤
x
<
(
q
+
1
)
/
2
p
0
e
l
s
e
h_u(x)= \begin{cases} 1& {u=0 \; and\; 0 \leq x<1 }\\ 2^{p/2}& {u>0\; and \; q/2^p \leq x<(q+0.5)/2^p}\\ -2^{p/2}& {u<0\; and \; (q+0.5)/2^p \leq x< (q+1)/2^p}\\ 0 & else \end{cases}
hu(x)=⎩⎪⎪⎪⎨⎪⎪⎪⎧12p/2−2p/20u=0and0≤x<1u>0andq/2p≤x<(q+0.5)/2pu<0and(q+0.5)/2p≤x<(q+1)/2pelse
可以看到哈儿基函数由两个变量控制,一个是
u
一
个
是
x
u一个是x
u一个是x。当
u
u
u不为
0
0
0则
h
u
(
x
)
h_u(x)
hu(x)在
[
q
/
2
p
,
x
<
(
q
+
0.5
)
/
2
p
)
[q/2^p,x<(q+0.5)/2^p)
[q/2p,x<(q+0.5)/2p)和
(
(
q
+
0.5
)
/
2
p
,
x
<
(
q
+
1
)
/
2
p
]
((q+0.5)/2^p,x< (q+1)/2^p]
((q+0.5)/2p,x<(q+1)/2p]分别是幅值为
2
p
/
2
和
−
2
p
/
2
2^{p/2}和-2^{p/2}
2p/2和−2p/2矩形波。
u
u
u可以控制小波的宽度,
u
u
u增大,矩形波就变窄。
哈儿变换的变换矩阵
A
H
A_H
AH可以由
s
(
x
,
u
)
=
1
N
h
u
(
x
/
N
)
s(x,u)=\frac{1}{N}h_u(x/N)
s(x,u)=N1hu(x/N)得到。其中,
u
=
0
,
1
,
2...
N
−
1
;
N
=
2
n
u=0,1,2...N-1;N=2^n
u=0,1,2...N−1;N=2n。
t
H
=
A
H
∗
f
t_H=A_H*f
tH=AH∗f
二、尺度函数
尺度函数通过将小波函数进行平移和拉伸得到。以
ϕ
(
x
)
\phi(x)
ϕ(x)作为父函数通过平移变换得到:
ϕ
j
,
k
(
x
)
=
2
j
/
2
ϕ
(
2
j
x
−
k
)
\phi_{j,k}(x)=2^{j/2}\phi(2^jx-k)
ϕj,k(x)=2j/2ϕ(2jx−k)
很显然这是一个函数集合,j控制幅度和宽度,k控制平移位移。
若尺度函数满足四个多分辨率分析的条件
则父函数可以表示为自身二倍分辨率副本的线性组合,即
j
=
1
j=1
j=1的情况,其尺度也为1
ϕ
(
x
)
=
∑
k
∈
Z
h
ϕ
(
k
)
2
ϕ
(
2
x
−
k
)
\phi(x)=\sum_{k\in Z}h_\phi(k)\sqrt{2}\phi(2x-k)
ϕ(x)=k∈Z∑hϕ(k)2ϕ(2x−k)
我们将 h ϕ ( k ) h_\phi(k) hϕ(k)称为尺度函数系数
三、小波变换
母小波函数
ψ
j
,
k
(
x
)
=
2
j
/
2
ψ
(
2
j
x
−
k
)
\psi_{j,k}(x)=2^{j/2}\psi(2^jx-k)
ψj,k(x)=2j/2ψ(2jx−k)
显然这是一个小波函数的平移缩放的函数集合
小波函数
ψ
(
x
)
\psi(x)
ψ(x)可以写为平移且缩放尺度为1的尺度函数的加权和。
ψ
(
x
)
=
∑
k
h
ψ
(
k
)
2
ϕ
(
2
x
−
k
)
\psi(x)=\sum_kh_\psi(k)\sqrt2\phi(2x-k)
ψ(x)=k∑hψ(k)2ϕ(2x−k)
其中
h
ϕ
(
k
)
h_\phi(k)
hϕ(k)为小波系数。
一维离散小波变换
f ( x ) = 1 ( N ) [ T ϕ ( 0 , 0 ) ϕ ( x ) + ∑ j = 0 J − 1 ∑ k = 0 2 j − 1 T ψ ( j , k ) ψ j , k ( x ) ] f(x)=\frac{1}{\sqrt(N)}[T_\phi(0,0)\phi(x)+\sum_{j=0}^{J-1}\sum_{k=0}^{2^j-1}T_\psi(j,k)\psi_{j,k}(x)] f(x)=(N)1[Tϕ(0,0)ϕ(x)+j=0∑J−1k=0∑2j−1Tψ(j,k)ψj,k(x)]
当尺度与小波为实值时:
T
ϕ
(
0.0
)
=
1
(
N
)
∑
x
=
0
N
−
1
f
(
x
)
ϕ
(
x
)
T_\phi(0.0)=\frac{1}{\sqrt(N)}\sum_{x=0}^{N-1}f(x)\phi(x)
Tϕ(0.0)=(N)1x=0∑N−1f(x)ϕ(x)
T
ψ
(
j
,
k
)
=
1
(
N
)
∑
x
=
0
N
−
1
f
(
x
)
ψ
j
,
k
(
x
)
T_\psi(j,k)=\frac{1}{\sqrt(N)}\sum_{x=0}{N-1}f(x)\psi_{j,k}(x)
Tψ(j,k)=(N)1x=0∑N−1f(x)ψj,k(x)
总结
摸鱼了,回头补上手工计算的过程。