查看“︁迭代稀疏漸近最小方差算法”︁的源代码

{{translating|[[:en:SAMV (algorithm)]]|||time=2019-06-27UTC01:11:00}}
{{NoteTA|G1=Signals and Systems}}
'''迭代稀疏漸近最小方差算法'''<ref name=AbeidaZhang/>是用於信號處理中的譜估計和到達方向（DOA）估計的無參數超分辨率算法。 這個名稱是為了強調漸近最小方差（AMV）標準的創造基礎。 它是在惡劣環境下恢復多個高相關源的幅度和頻率特性的有力工具，例如有限數量的快照，低信噪比。 它可以用於[[合成孔徑雷達]]<ref name=SAR/><ref name=Yang/>。

== 概述 ==
迭代稀疏漸近最小方差算法是一種基於[[壓縮感知]]的[[超高解析度成像]]程式, 可以用於[[合成孔徑雷達]]成像, [[訊號處理]], [[核磁共振成像]]等[[醫學影像]]領域。

== 定義 ==

SAMV算法的公式在DOA估計的背景下作為[[反問題]]給出。假設 <math> M </math> -元素 [[傳感器陣列|均勻線性陣列]]（ULA）分別接收從位於<math>\mathbf{\theta} = \{\theta_a, \ldots, \theta_K \}</math> 位置發出的<math>K</math>窄帶信號。 ULA中的傳感器在特定時間累積<math>N</math>快照。 <math>M\times 1</math>維快照向量是

: <math>\mathbf{y}(n) = \mathbf{A} \mathbf{x}(n) + \mathbf{e}(n), n = 1, \ldots, N</math>

其中 <math>\mathbf{A} = [ \mathbf{a}(\theta_1), \ldots, \mathbf{a}(\theta_K) ]</math> 是[[相控陣|轉向矩陣]], <math>{\bf x}(n)=[{\bf x}_1(n), \ldots, {\bf x}_K(n)]^T </math> 包含源波形， 和 <math>{\bf e}(n)</math> 是噪音詞。假設<math>\mathbf{E}\left({\bf e}(n){\bf e}^H(\bar{n})\right)= \sigma{\bf I}_M\delta_{n,\bar{n}}</math>, <math>\delta_{n,\bar{n}}</math> 是 [[Dirac delta function|Dirac delta 函數]] 並且它僅等於1，唯一存在 <math>n=\bar{n}</math> 否則為0。並且假設 <math>{\bf e}(n)</math> and <math>{\bf x}(n)</math> 是獨立的，而 <math>\mathbf{E}\left({\bf x}(n){\bf x}^H(\bar{n})\right)={\bf P}\delta_{n,\bar{n}}</math>, where <math>{\bf P}= \operatorname{Diag}( {p_1,\ldots,p_K})</math>. Let <math>{\bf p}</math> 是包含未知信號功率和噪聲方差的向量， <math>{\bf p} =  [p_1,\ldots,p_K, \sigma]^T</math>.

<math>{\bf y}(n)</math>的協方差矩陣，其中有關<math>\boldsymbol{\bf p}</math>的是
: <math>{\bf R} = {\bf A}{\bf P}{\bf A}^H+\sigma{\bf I}. </math>
該協方差矩陣可以通過樣本協方差矩陣進行傳統估計 <math>{\bf R}_{N} = {\bf Y}{\bf Y}^H/N</math> ，其中 <math>{\bf Y} =[{\bf y}(1), \ldots,{\bf y}(N)]</math> 。將向量化運算符應用於矩陣 <math>{\bf R}</math> 後，獲取的向量 <math>{\bf r}(\boldsymbol{\bf p}) = \operatorname{vec}({\bf R})</math> 與未知參數線性相關<math>\boldsymbol{\bf p}</math> 當

<math> {\bf r}(\boldsymbol{\bf p}) = \operatorname{vec}({\bf R})={\bf S}\boldsymbol{\bf p}</math>,

其中  <math>{\bf S}=  [{\bf S}_1,\bar{\bf a}_{K+1}]</math>, <math>{\bf S}_1 =[\bar{\bf a}_1,\ldots,\bar{\bf a}_K]</math>, <math>\bar{\bf a}_k = {\bf a}^{*}_k \otimes{\bf a}_k</math>, <math>k=1,\ldots, K</math>, 和使 <math>\bar{\bf a}_{K+1} = \operatorname{vec}({\bf I})</math>.

== SAMV演算法 == 
要從統計的<math>{\bf r}_N</math>去估算<math>\boldsymbol{\bf p}</math> ，我們基於漸近最小方差準則開發了一系列迭代SAMV方法。從<ref name=AbeidaZhang/>開始，從協方差矩陣 <math>\operatorname{Cov}^\operatorname{Alg}_{\boldsymbol{p}}</math> 的任意一致的估計值 <math>\boldsymbol{p}</math> ，基於二階統計值 <math>{\bf r}_N</math> ，以實數對稱-正定矩陣為界

: <math>\operatorname{Cov}^\operatorname{Alg}_{\boldsymbol{p}}\geq[{\bf S}^H_d {\bf C}^{-1}_r{\bf S}_d]^{-1},</math>

其中 <math>{\bf S}_d =  {\rm d}{\bf r}(\boldsymbol{p})/ {\rm d}\boldsymbol{p}</math>。此外，這個下界是通過最小化得到的<math>\hat{\bf p}</math> 的漸近分佈的協方差矩陣得到的。 ，

: <math>\hat{\boldsymbol{p}} =\arg \min_{\boldsymbol{p}} f(\boldsymbol{p}),</math>

其中
<math>f(\boldsymbol{p}) = [{\bf r}_N-{\bf r}(\boldsymbol{p})]^H {\bf C}_r^{-1} [{\bf r}_N-{\bf r}(\boldsymbol{p})].</math>

因此，可以迭代地獲 <math>\boldsymbol{\bf p}</math> 的估計值。
<math>\{\hat{p}_k\}_{k=1}^K</math> 和最小化<math>f(\boldsymbol{p})</math>的<math>\hat{\sigma}</math> 可藉由以下計算獲得。

假設 <math>\hat{p}^{(i)}_k</math> 和 <math>\hat{\sigma}^{(i)}</math> 在第 <math>i</math> 迭代中已被估算到某種程度, 第 <math>(i+1)</math>迭代可以被精簡成，

: <math>
\hat{p}^{(i+1)}_k = \frac{{\bf a}^H_k{\bf R}^{-1{(i)}}{\bf R}_N {\bf R}^{-1{(i)}}{\bf a}_k}{ ({\bf a}^H_k{\bf R}^{-1{(i)}}{\bf a}_k)^2}+\hat{p}^{(i)}_k-\frac{1}{{\bf a}^H_k{\bf R}^{-1{(i)}}{\bf a}_k}, \quad k=1, \ldots,K</math>

: <math>
\hat{\sigma}^{(i+1)} = \left(\operatorname{Tr}({\bf R}^{-2^{(i)}}{\bf R}_N) + \hat{\sigma}^{(i)}\operatorname{Tr}({\bf R}^{-2^{(i)}}) -\operatorname{Tr}({\bf R}^{-1^{(i)}})\right)/{\operatorname{Tr}{({\bf R}^{-2^{(i)}})}},</math>

其中 <math>{\bf R}</math> 的估計值在第 <math>i</math> 迭代是 <math>{\bf R}^{(i)}={\bf A}{\bf P}^{(i)}{\bf A}^H+\hat{\sigma}^{(i)}{\bf I}</math> with <math>{\bf P}^{(i)}=\operatorname{Diag}(\hat{p}^{(i)}_1, \ldots, \hat{p}^{(i)}_K)</math>.

== 超越掃描網格精度 ==
基於大多數壓縮感知的源定位技術的分辨率受到覆蓋位置參數空間的方向網格的精細度的限制。<ref name="Malioutov Cetin Willsky 2005 pp. 3010–3022"/> 在稀疏信號恢復模型中，真值信號的稀疏性 <math>\mathbf{x}(n)</math> 取決於超完備字典<math>{\bf A}</math>中相鄰元素之間的距離因此, 會出現選擇最佳超完備字典的難度。計算複雜度與方向網格的精細度成正比，高密度網格在計算上不實用。為了克服網格強加的分辨率限制，提出了無網格SAMV-SML（迭代稀疏漸近最小方差 - 隨機最大似然<ref name=AbeidaZhang/>, 它藉由迭代的最小化隨機[[最大似然估計]]的消耗函數，相對於單一純數<math>\theta_k</math>，改進了位置估計 <math>\boldsymbol{\bf \theta}=(\theta_1,\ldots,\theta_K)^T</math>。

==應用於距離多普勒成像==

[[File:IEEE transaction on Signal Processing Paper Results Sample.jpg|thumb|upright=1.5|[[Single-input single-output system|SISO]]範圍多普勒成像結果與三個5 dB和六個25 dB目標進行比較。 （a）基礎事實，（b）匹配濾波器（MF），（c）IAA算法，（d）SAMV-0算法。所有功率水平均以dB為單位。相對於多普勒軸，MF和IAA方法的分辨率都受到限制。 SAMV-0在範圍和多普勒方面提供卓越的分辨率。 <ref name = AbeidaZhang />]]

在[[單輸入單輸出系統| SISO]] [[雷達]] / [[聲納]] [[脈衝多普勒雷達|距離 - 多普勒成像]]問題中使用SAMV算法的典型應用。該成像問題是單快照應用，並且包括與單快照估計兼容的算法，即[[匹配濾波器]]（MF，類似於[[週期圖]]或[[Radon變換|反投影]]，這通常被有效地實現為[[快速傅立葉變換]]（FFT）），IAA <ref name="Yardibi Li Stoica Xue 2010 pp. 425–443"/>和SAMV算法的變體（SAMV-0）。 模擬條件與之相同<ref name="Yardibi Li Stoica Xue 2010 pp. 425–443"/>: 一個 <math>30</math>-元素的多項 [[脈沖壓縮|pulse compression]]使用P3代碼相同作為發射脈衝，模擬總共九個運動目標。在所有移動目標中，三個是<math>5</math> dB功率，其餘六個是<math>25</math> dB功率。假設接收信號被<math> 0 </math> dB功率的均勻高斯白噪聲污染。

[[匹配濾波器]]檢測結果在多普勒和範圍域都受到嚴重的拖尾和[[光譜洩漏]]影響，因此無法區分<math>5</math> dB目標。相反，IAA算法提供增強的成像結果，具有可觀察的目標範圍估計和多普勒頻率。 SAMV-0方法提供高度稀疏的結果並完全消除拖尾效應，但它錯過了弱<math>5</math> dB目標。


== 參考資料 ==
{{reflist|refs=<ref name=AbeidaZhang>{{cite journal | last=Abeida | first=Habti | last2=Zhang | first2=Qilin | last3=Li | first3=Jian | last4=Merabtine | first4=Nadjim | title=Iterative Sparse Asymptotic Minimum Variance Based Approaches for Array Processing | journal=IEEE Transactions on Signal Processing | volume=61 | issue=4 | year=2013 | issn=1053-587X | doi=10.1109/tsp.2012.2231676 | pages=933–944 | url=https://qilin-zhang.github.io/_pages/pdfs/SAMVpaper.pdf | bibcode=2013ITSP...61..933A | arxiv=1802.03070 | access-date=2018-09-14 | archive-date=2021-02-24 | archive-url=https://web.archive.org/web/20210224162708/https://qilin-zhang.github.io/_pages/pdfs/SAMVpaper.pdf | dead-url=no }}</ref>

<ref name=SAR>{{cite journal | last=Glentis | first=George-Othon | last2=Zhao | first2=Kexin | last3=Jakobsson | first3=Andreas | last4=Abeida | first4=Habti | last5=Li | first5=Jian | title=SAR imaging via efficient implementations of sparse ML approaches | journal=Signal Processing | volume=95 | year=2014 | doi=10.1016/j.sigpro.2013.08.003 | pages=15–26}}</ref>

<ref name=Yang>{{cite journal | last=Yang | first=Xuemin | last2=Li | first2=Guangjun | last3=Zheng | first3=Zhi | title=DOA Estimation of Noncircular Signal Based on Sparse Representation | journal=Wireless Personal Communications | volume=82 | issue=4 | date=2015-02-03 | doi=10.1007/s11277-015-2352-z | pages=2363–2375}}</ref>

<ref name="Malioutov Cetin Willsky 2005 pp. 3010–3022">{{cite journal | last=Malioutov | first=D. | last2=Cetin | first2=M. | last3=Willsky | first3=A.S. | title=A sparse signal reconstruction perspective for source localization with sensor arrays | journal=IEEE Transactions on Signal Processing | volume=53 | issue=8 | year=2005 | doi=10.1109/tsp.2005.850882 | pages=3010–3022| bibcode=2005ITSP...53.3010M }}</ref>

<ref name="Yardibi Li Stoica Xue 2010 pp. 425–443">{{cite journal | last=Yardibi | first=Tarik | last2=Li | first2=Jian | last3=Stoica | first3=Petre | last4=Xue | first4=Ming | last5=Baggeroer | first5=Arthur B. | title=Source Localization and Sensing: A Nonparametric Iterative Adaptive Approach Based on Weighted Least Squares | journal=IEEE Transactions on Aerospace and Electronic Systems | volume=46 | issue=1 | year=2010 | doi=10.1109/taes.2010.5417172 | pages=425–443| bibcode=2010ITAES..46..425Y }}</ref>}}

[[Category:醫學成像]]
[[Category:信號處理]]
[[Category:斷層掃描]]
[[Category:估计理论]]
[[Category:電信理論]]
[[Category:信号处理]]
[[Category:滤波器理论]]