现假设有一组基元
b
0
(
t
)
,
b
1
(
t
)
,
b
2
(
t
)
,
b
3
(
t
)
.
.
.
.
.
.
b
N
(
t
)
{\displaystyle b_{0}(t),b_{1}(t),b_{2}(t),b_{3}(t)......b_{N}(t)}
源自于字典
D
{\displaystyle D}
,这组
b
n
(
t
)
{\displaystyle b_{n}(t)}
组成一个过完备性 (Over-complete)、非正交集合。
压缩感知的问题在此即为:
问题1:完全相等
能不能找到最小的
|
|
a
|
|
0
{\displaystyle ||a||_{0}}
(
a
k
{\displaystyle a_{k}}
不为0的个数,
∀
k
=
{
1
,
2
,
.
.
.
,
N
}
{\displaystyle \forall k=\{1,2,...,N\}}
),使得近似讯号
f
^
{\displaystyle {\hat {f}}}
等于
f
{\displaystyle f}
。
f
(
t
)
=
∑
m
a
m
b
m
(
t
)
{\displaystyle f(t)=\sum _{m}a_{m}b_{m}(t)}
问题2:近似问题
很明显的,自然界中讯号大多只能以近似方式逼近,因此改为求
min
|
|
a
|
|
0
∫
|
|
f
(
t
)
−
∑
m
a
m
b
m
(
t
)
|
|
2
d
t
<
t
h
r
e
s
h
o
l
d
{\displaystyle \min _{||a||_{0}}\int ||f(t)-\sum _{m}a_{m}b_{m}(t)||^{2}dt<threshold}
然而,这个最佳化问题是NP困难 ,无法在线性多项式时间内找到解答,因此使用匹配追求的近似解来求解。
问题3:匹配追求
min
∫
|
|
f
(
t
)
−
∑
m
a
m
b
m
(
t
)
|
|
2
d
t
, such that
|
|
a
|
|
0
≤
N
{\displaystyle \min \int ||f(t)-\sum _{m}a_{m}b_{m}(t)||^{2}dt{\text{ , such that }}||a||_{0}\leq N}
算法
有鉴于匹配追求的限制条件以及贪婪算法在特定情况上造成不适当的基元选择,S.S.Chen等人于1998年提出的基底追求。[ 1] 若将匹配追求想成,从一个"空"的集合中,在每次迭代中慢慢增加一个基元来构成近似讯号,基底追求则可以想成,从一个完整的基元集合,慢慢减少基元数量。
现假设有一组基元
b
0
(
t
)
,
b
1
(
t
)
,
b
2
(
t
)
,
b
3
(
t
)
.
.
.
.
.
.
b
N
(
t
)
{\displaystyle b_{0}(t),b_{1}(t),b_{2}(t),b_{3}(t)......b_{N}(t)}
源自于字典
D
{\displaystyle D}
,这组
b
n
(
t
)
{\displaystyle b_{n}(t)}
组成一个过完备性(Over-complete)、非正交集合。不同于匹配追求,将原先所求的0级范数 (
|
|
a
|
|
0
{\displaystyle ||a||_{0}}
),改为求其绝对值
|
|
a
|
|
1
=
|
a
0
|
+
|
a
1
|
+
|
a
2
|
+
|
a
3
|
+
.
.
.
.
.
.
+
|
a
N
|
{\displaystyle ||a||_{1}=|a_{0}|+|a_{1}|+|a_{2}|+|a_{3}|+......+|a_{N}|}
。
在此情况下,基底追求对于压缩感知问题的解法亦能细分为三大问题:
问题1:完全相等
能不能找到最小的
|
|
a
|
|
1
{\displaystyle ||a||_{1}}
,使得近似讯号
f
^
{\displaystyle {\hat {f}}}
等于
f
{\displaystyle f}
。
f
(
t
)
=
∑
m
a
m
b
m
(
t
)
{\displaystyle f(t)=\sum _{m}a_{m}b_{m}(t)}
问题2:近似问题
同理,自然界中讯号大多只能以近似方式逼近,因此改为求
min
|
|
a
|
|
1
∫
|
|
f
(
t
)
−
∑
m
a
m
b
m
(
t
)
|
|
2
d
t
<
t
h
r
e
s
h
o
l
d
{\displaystyle \min _{||a||_{1}}\int ||f(t)-\sum _{m}a_{m}b_{m}(t)||^{2}dt<threshold}
这个最佳化问题依旧是NP困难,因此可使用基底追求的近似解来求解。
问题3:基底追求
min
∫
|
|
f
(
t
)
−
∑
m
a
m
b
m
(
t
)
|
|
2
d
t
, such that
|
|
a
|
|
1
≤
N
{\displaystyle \min \int ||f(t)-\sum _{m}a_{m}b_{m}(t)||^{2}dt{\text{ , such that }}||a||_{1}\leq N}
相关型态
Three-Parameter Atoms
f
(
t
)
≈
f
^
(
t
)
=
∑
a
t
0
,
f
0
,
σ
φ
t
0
,
f
0
,
σ
(
t
)
{\displaystyle f(t)\approx {\hat {f}}(t)=\sum a_{t_{0},f_{0},\sigma }\varphi _{t_{0},f_{0},\sigma }(t)}
where
φ
t
0
,
f
0
,
σ
(
t
)
=
2
1
/
4
σ
1
/
2
e
j
2
π
f
0
t
−
π
(
t
−
t
0
)
2
σ
2
{\displaystyle \varphi _{t_{0},f_{0},\sigma }(t)={\frac {2^{1/4}}{\sigma ^{1/2}}}e^{j2\pi f_{0}t-{\frac {\pi (t-t_{0})^{2}}{\sigma ^{2}}}}}
近似讯号
f
^
{\displaystyle {\hat {f}}}
是
N
=
3
{\displaystyle N=3}
的特例,由三个基元组合而成,每个基元由
t
0
{\displaystyle t_{0}}
决定在时间轴上的中心位置,
f
0
{\displaystyle f_{0}}
决定在频率轴上的中心位置,以及由
σ
{\displaystyle \sigma }
选择缩放尺度的大小。
由于
φ
t
0
,
f
0
,
σ
{\displaystyle \varphi _{t_{0},f_{0},\sigma }}
是非正交集合,
a
t
0
,
f
0
,
σ
{\displaystyle a_{t_{0},f_{0},\sigma }}
需要透过匹配追求来求得。
Four-Parameter Atoms(Chirplet) [ 2]
f
(
t
)
≈
f
^
(
t
)
=
∑
a
t
0
,
f
0
,
σ
,
η
⋅
φ
t
0
,
f
0
,
σ
,
η
(
t
)
{\displaystyle f(t)\approx {\hat {f}}(t)=\sum a_{t_{0},f_{0},\sigma ,\eta }\cdot \varphi _{t_{0},f_{0},\sigma ,\eta }(t)}
where
φ
t
0
,
f
0
,
σ
,
η
(
t
)
=
2
1
/
4
σ
1
/
2
e
j
2
π
(
f
0
t
+
η
2
t
2
)
−
π
(
t
−
t
0
)
2
σ
2
{\displaystyle \varphi _{t_{0},f_{0},\sigma ,\eta }(t)={\frac {2^{1/4}}{\sigma ^{1/2}}}e^{j2\pi (f_{0}t+{\frac {\eta }{2t^{2}}})-{\frac {\pi (t-t_{0})^{2}}{\sigma ^{2}}}}}
近似讯号
f
^
{\displaystyle {\hat {f}}}
是
N
=
4
{\displaystyle N=4}
的特例,由四个基元组合而成,每个基元由
t
0
{\displaystyle t_{0}}
决定在时间轴上的中心位置,
f
0
{\displaystyle f_{0}}
决定在频率轴上的中心位置,以及
σ
{\displaystyle \sigma }
选择缩放尺度的大小和
η
{\displaystyle \eta }
控制啁啾率。
同理,由于
φ
t
0
,
f
0
,
σ
,
η
{\displaystyle \varphi _{t_{0},f_{0},\sigma ,\eta }}
是非正交集合,
a
t
0
,
f
0
,
σ
,
η
{\displaystyle a_{t_{0},f_{0},\sigma ,\eta }}
需要透过匹配追求来求得。
参考资料
1. Jian-Jiun Ding, Time frequency analysis and wavelet transform class note, Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2017.
外部链接
^ Chen, S. S.; Donoho, D. L.; Saunders, M. A. Atomic Decomposition by Basis Pursuit. SIAM Journal on Scientific Computing. 1998, 20 (1): 33–61. doi:10.1137/S003614450037906X .
^ Bultan, A. A four-parameter atomic decomposition of chirplets. IEEE Transactions on Signal Processing. 1999, 47 (3): 731–745. doi:10.1109/78.747779 .