查看“︁匹配追求過程”︁的源代码

{{cleanup-jargon|time=2018-01-08T08:35:17+00:00}}
{{howto|time=2018-01-08T08:35:17+00:00}}
'''匹配追求'''，是一個將高複雜度的信號以較簡單的訊號來近似的過程，也就是在多維空間中，將高維度的資訊投影到較低維度的[[生成空間]]，藉此降低一個信號所需的展開項目。

簡單來說，就是用盡可能少的基元，進行線性組合來逼近原訊號。

一個有N個基元的近似訊號為
: <math> f(t) \approx \hat f_N(t) = \sum_{n=1}^{N} a_n b_n(t)</math>

其<math>a_n</math>為各基元的權重。而與原訊號<math>f</math>的近似誤差為<math>R</math>，表示為
: <math> R  = f(t) - \hat f(t)</math>

匹配追求並不會賦予子空間<math>D</math>的每個基元一個權重<math>a_n</math>，意即有些基底會捨棄不用，藉此降低所需空間，相反的，匹配追求會基於[[貪婪演算法]]，盡可能的用較少的基元，來降低誤差<math>R</math>。核心思想為自子空間<math>D</math>中選擇一個和訊號<math>f</math>內積值最大的基元<math>b_i</math>，將基元乘上適當的權重<math>a_i</math>並從原訊號中減去，接著再次進行以上二步驟。以此往復實行運算，直到原訊號被解析，也就是近似誤差<math>R</math>趨近於0。

==[[壓縮感知]]==
現假設有一組基元<math>b_0(t), b_1(t), b_2(t), b_3(t)......b_N(t)</math>源自於字典<math>D</math>，這組<math>b_n(t)</math>組成一個過[[完備性]](Over-complete)、非正交集合。

壓縮感知的問題在此即為：

===問題1：完全相等===

能不能找到最小的<math>||a||_0</math>(<math>a_k</math>不為0的個數，<math>\forall k=\{1, 2, ..., N\}</math>)，使得近似訊號<math>\hat f</math>等於<math>f</math>。

: <math> f(t) = \sum_m a_m b_m(t) </math>

===問題2：近似問題===

很明顯的，自然界中訊號大多只能以近似方式逼近，因此改為求

: <math> \min_{||a||_0}  \int  || f(t) - \sum_m a_m b_m(t) || ^2 dt < threshold  </math>

然而，這個最佳化問題是[[NP困難]]，無法在線性多項式時間內找到解答，因此使用匹配追求的近似解來求解。

===問題3：匹配追求===

: <math> \min  \int  || f(t) - \sum_m a_m b_m(t) || ^2 dt \text{ , such that }  ||a||_0 \leq N </math>

==演算法==

===匹配追求(貪婪演算法)===

  輸入： 原訊號 <math>f(t)</math>
  輸出： 所需權重 <math> a_n </math> 與其對應的基元 <math> b_n </math>
  初始化：
    <math>n \leftarrow 0</math>
    <math>\hat f(t) \leftarrow f(t)</math>
  反覆：
    尋找 <math> m </math> 使得內積 <math> \int \hat f(t)b^*_m(t)dt </math> 最大
    令 <math>\phi _n(t) \leftarrow b_m(t)</math>
    令 <math>\mu _n \leftarrow \int \hat f(t)b^*_m(t)dt</math>
    <math>\hat f(t) \leftarrow \hat f(t) - \mu _n \phi _n(t)</math>
    <math>n \leftarrow n + 1</math>
  直到中止條件發生：
    問題一：  <math>\hat f(t) = 0</math>
    問題二：  <math>\int \hat f(t)^2 (t) dt < threshold </math>
    問題三：  <math>n = N</math>
  結束

==[[基底追求]](Basis Pursuit)==
有鑑於匹配追求的限制條件以及貪婪演算法在特定情況上造成不適當的基元選擇，S.S.Chen等人於1998年提出的基底追求。<ref>{{cite journal | last1 = Chen | first1 = S. S. | last2 = Donoho | first2 = D. L. | last3 = Saunders | first3 = M. A.| year = 1998 | title = Atomic Decomposition by Basis Pursuit | url = | journal = SIAM Journal on Scientific Computing | volume = 20 | issue = 1 | pages = 33-61 | doi = 10.1137/S003614450037906X }}</ref>若將匹配追求想成，從一個"空"的集合中，在每次迭代中慢慢增加一個基元來構成近似訊號，基底追求則可以想成，從一個完整的基元集合，慢慢減少基元數量。

現假設有一組基元<math>b_0(t), b_1(t), b_2(t), b_3(t)......b_N(t)</math>源自於字典<math>D</math>，這組<math>b_n(t)</math>組成一個過完備性(Over-complete)、非正交集合。不同於匹配追求，將原先所求的0級[[範數]](<math>||a||_0</math>)，改為求其絕對值<math>||a||_1 = |a_0| + |a_1| + |a_2| + |a_3| + ...... + |a_N|</math>。

在此情況下，基底追求對於壓縮感知問題的解法亦能細分為三大問題：

===問題1：完全相等===

能不能找到最小的<math>||a||_1</math>，使得近似訊號<math>\hat f</math>等於<math>f</math>。

: <math> f(t) = \sum_m a_m b_m(t) </math>

===問題2：近似問題===

同理，自然界中訊號大多只能以近似方式逼近，因此改為求

: <math> \min_{||a||_1}  \int  || f(t) - \sum_m a_m b_m(t) || ^2 dt < threshold  </math>

這個最佳化問題依舊是NP困難，因此可使用基底追求的近似解來求解。

===問題3：基底追求===

: <math> \min  \int  || f(t) - \sum_m a_m b_m(t) || ^2 dt \text{ , such that }  ||a||_1 \leq N </math>

==相關型態==

===Three-Parameter Atoms===
:<math>f(t) \approx \hat f(t) = \sum a_{t_0,f_0,\sigma}  \varphi_{t_0,f_0,\sigma}(t)</math>

: where <math>\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} }</math>

近似訊號<math>\hat f</math>是<math>N = 3</math>的特例，由三個基元組合而成，每個基元由<math>t_0 </math>決定在時間軸上的中心位置，<math>f_0</math>決定在頻率軸上的中心位置，以及由<math>\sigma</math>選擇縮放尺度的大小。

由於<math>\varphi_{t_0,f_0,\sigma}</math> 是非正交集合，<math>a_{t_0, f_0, \sigma}</math> 需要透過匹配追求來求得。

===Four-Parameter Atoms(Chirplet) <ref>{{cite journal | last1 = Bultan | first1 = A. | year = 1999 | title = A four-parameter atomic decomposition of chirplets | url = | journal = IEEE Transactions on Signal Processing | volume = 47 | issue = 3 | pages = 731-745 | doi = 10.1109/78.747779 }}</ref> ===
:<math>f(t) \approx \hat f(t) = \sum a_{t_0,f_0, \sigma, \eta}\cdot \varphi_{t_0, f_0, \sigma, \eta}(t)</math>

: where <math>\varphi_{t_0, f_0, \sigma, \eta}(t) = \frac{2^{1/4}}{\sigma^{1/2}}e^{j2\pi(f_0t + \frac{\eta}{2t^2})- \frac{\pi (t-t_0)^2}{\sigma^2} }</math>

近似訊號<math>\hat f</math>是<math>N = 4</math>的特例，由四個基元組合而成，每個基元由<math>t_0 </math>決定在時間軸上的中心位置，<math>f_0</math>決定在頻率軸上的中心位置，以及<math>\sigma</math>選擇縮放尺度的大小和<math>\eta</math>控制啁啾率。

同理，由於<math>\varphi_{t_0,f_0,\sigma, \eta}</math> 是非正交集合，<math>a_{t_0, f_0, \sigma, \eta}</math> 需要透過匹配追求來求得。




==參考資料==
1. Jian-Jiun Ding, Time frequency analysis and wavelet transform class note, Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2017.

==外部連結==