劉開華，陳偉凱，馬永濤

(天津大學(xué)電子信息工程學(xué)院，天津 300072)

正交頻分復(fù)用(orthogonal frequency division multiplexing，OFDM)是下一代移動(dòng)通信的核心技術(shù)，所采用的多載波高速傳輸技術(shù)允許子信道的頻譜相互重疊，從而可以實(shí)現(xiàn)極高的頻譜利用率．然而，OFDM 系統(tǒng)對(duì)相位噪聲和載波頻偏非常敏感，這對(duì)子載波之間的正交性提出了嚴(yán)格的要求．因此，在OFDM 系統(tǒng)接收端，往往需要知道精確的信道狀態(tài)信息(channel state information，CSI)才能有效對(duì)抗信道衰落，實(shí)現(xiàn)信號(hào)的準(zhǔn)確接收．

在實(shí)際的無線信號(hào)傳輸中，雙選擇性的多徑信道通常由少數(shù)的主要路徑簇所主導(dǎo)，因此所呈現(xiàn)的物理信道常具有稀疏特性．而當(dāng)信號(hào)的傳輸帶寬較大或天線個(gè)數(shù)較多時(shí)，信道的稀疏特性尤為明顯[1]． 由于稀疏信道只有少數(shù)非零抽頭，傳統(tǒng)的基于導(dǎo)頻序列的方法(如基于導(dǎo)頻的線性插值和 FFT插值法)極有可能采樣到信道的零抽頭，而無法準(zhǔn)確地插值出信道響應(yīng)[2]．為了突破傳統(tǒng)算法的瓶頸，近幾年來，國內(nèi)外相繼有一些學(xué)者將壓縮感知技術(shù)[3-6]應(yīng)用到對(duì)稀疏信道的估計(jì)[7-12]中，其中包括 OFDM 系統(tǒng)[9-10]、超寬帶(ultra wide band，UWB)系統(tǒng)[11]和水聲通信系統(tǒng)[12]等．作為應(yīng)用數(shù)學(xué)和信號(hào)處理領(lǐng)域的一大突破，壓縮感知(compressive sensing，CS)理論表示當(dāng)信號(hào)是可壓縮的或是在某個(gè)變換域具有稀疏性時(shí)，通過采樣少量信號(hào)就能實(shí)現(xiàn)信號(hào)的準(zhǔn)確或近似重構(gòu)[3]．在該理論框架下，壓縮感知技術(shù)可以充分挖掘信道的稀疏特性，能利用非常有限的導(dǎo)頻有效地恢復(fù)稀疏的信道脈沖響應(yīng)，在實(shí)現(xiàn)對(duì)稀疏信道準(zhǔn)確估計(jì)的同時(shí)，還能極大減少導(dǎo)頻數(shù)目，提高系統(tǒng)的頻譜利用率．

在已有的研究成果中，文獻(xiàn)[7]最先提出將壓縮感知的方法應(yīng)用到對(duì)稀疏信道的估計(jì)中去，并具體采用匹配追蹤(matching pursuit，MP)算法實(shí)現(xiàn)了信道估計(jì)．文獻(xiàn)[8]研究了在頻率選擇性條件下對(duì)信道進(jìn)行壓縮感知，但是并未將其應(yīng)用到多載波系統(tǒng)，也并未考慮到信道的時(shí)間選擇性．文獻(xiàn)[9]和文獻(xiàn)[10]分別應(yīng)用群匹配追蹤(group matching pursuit，GMP)和子空間追蹤(subspace pursuit，SP)的方法對(duì) OFDM 系統(tǒng)進(jìn)行信道估計(jì)，但是二者都沒有考慮到時(shí)間和頻率的雙選擇性對(duì)信道估計(jì)的影響．為此，本文重點(diǎn)研究了 OFDM 系統(tǒng)在雙選擇性衰落下的信道模型，并利用基追蹤(basis pursuit，BP)算法[13]實(shí)現(xiàn)對(duì)稀疏信道的準(zhǔn)確重構(gòu)．在模型的建立過程中，將 OFDM 系統(tǒng)的雙選擇性信道模型轉(zhuǎn)化為 CS可解的 BPIC(basis pursuit inequality constraint)模型，并利用導(dǎo)頻處的信道響應(yīng)作為 CS重構(gòu)所需的測(cè)量值．仿真結(jié)果顯示，新方法可以大大減少 OFDM 系統(tǒng)的導(dǎo)頻數(shù)，且在FFT-Linear和 FFT-FFT二維插值算法無法準(zhǔn)確估計(jì)出信道響應(yīng)的同時(shí)，BP算法仍能實(shí)現(xiàn)較高的估計(jì)精度．

1 壓縮感知理論

壓縮感知技術(shù)是針對(duì)稀疏信號(hào)的新型數(shù)據(jù)采樣和壓縮技術(shù)，其理論核心可以簡(jiǎn)述為：對(duì)于一個(gè)長度為N的一維離散信號(hào) x，若 x是K(K?N)稀疏的(即只有K個(gè)非零元素)，則通過選取合適的測(cè)量矩陣以獲得x在該矩陣下的線性測(cè)量值為

那么僅利用y中的M個(gè)測(cè)量值就能以很大概率重構(gòu)原始信號(hào)x．由于y的維數(shù)要遠(yuǎn)低于x的維數(shù)，因此式(1)有無窮多個(gè)解，即該問題是不適定的．然而，Candès等[6]證明，若測(cè)量次數(shù)M能滿足M＝O(Klog(N))，且測(cè)量矩陣Φ能滿足約束等距條件，則可以通過求解一個(gè)最小l0范數(shù)問題來重構(gòu)信號(hào)x，即

但 Donoho[14]指出，最小l0范數(shù)問題仍然是個(gè)NP-Hard問題，因而難以求解．為此，研究人員提出了一系列可求得近似最優(yōu)解的算法，包括MP系列算法[15-17]、BP算法[13]、迭代閾值法[18]以及專門用于圖像處理的最小全變分法[3]等．其中，以 MP系列算法和 BP算法最為常用．MP系列算法以貪婪迭代為核心思想，其運(yùn)算量相對(duì)較小，但其理論保證比基追蹤算法要弱，并非對(duì)所有的信號(hào)都能重構(gòu)，且對(duì)測(cè)量矩陣的要求比約束等距性更加嚴(yán)格．再者，MP系列算法對(duì)維數(shù)較低的小尺度信號(hào)問題運(yùn)算速度很快，但對(duì)于噪聲存在的大尺度問題，重構(gòu)結(jié)果不是很準(zhǔn)確，也不具有魯棒性．考慮實(shí)際信道中的噪聲條件，為保證OFDM 系統(tǒng)中所有的雙選擇性衰落信道都能采用壓縮感知的方法進(jìn)行準(zhǔn)確估計(jì)，本文采用精度更高的BP算法．

BP算法用l1范數(shù)替代了l0范數(shù)，即得到如下問題：

可以證明式(3)和式(2)在理論上是等價(jià)的[13].由于問題式(3)實(shí)際上是凸優(yōu)化問題，可以轉(zhuǎn)化為線性規(guī)劃問題加以求解．如果考慮重構(gòu)誤差及噪聲，則可以將式(3)轉(zhuǎn)化為BPIC模型，即

該問題可以利用二階圓錐規(guī)劃來求解．

2 雙選擇性衰落條件下 OFDM系統(tǒng)的信道模型

考慮單天線 OFDM 通信系統(tǒng)，通常情況下該系統(tǒng)可以建模為線性時(shí)變系統(tǒng)．在不考慮噪聲的情況下，其基帶的發(fā)射和接收信號(hào)的關(guān)系為

式中：x(t)和y(t)分別表示系統(tǒng)的發(fā)射和其對(duì)應(yīng)的接收信號(hào)；h(t，τ)代表信道的時(shí)變沖激響應(yīng)；S(υ，τ)是時(shí)延-多普勒分布函數(shù)；X(f)和H(t，f)分別是對(duì)應(yīng)x(t)和h(t，τ)的傅里葉變換；τmax和υmax是刻畫信道特性的重要參數(shù)，τmax代表信道所引入的最大時(shí)延擴(kuò)展，而υmax/2則代表信道的最大單邊多普勒頻移．假設(shè)碼元的符號(hào)寬度為T，基帶信號(hào)的帶寬為W，若信道為時(shí)間和頻率雙選擇性信道，則系統(tǒng)的符號(hào)寬度和信號(hào)帶寬都將大于其對(duì)應(yīng)的相干時(shí)間 1/υmax和相干帶寬 1/τmax，因此存在Wτmax≥1 和Tυmax≥1．

離散的信道時(shí)延-多普勒分布函數(shù)可用來表征雙選擇性信道的多徑效應(yīng)和多普勒頻移特性，即

式中：Npath表示多徑信道中存在的路徑總數(shù)；αi∈C代表第i個(gè)路徑的路徑增益；υi∈[－υmax/2，υmax/2]和τi∈[0，τmax]分別代表第i個(gè)路徑中的時(shí)延和多普勒頻移．因此，雙選擇性信道的信道響應(yīng)可以表示為

雖然式(6)很好地刻畫了物理多徑信道的傳輸特性，但由于它與信道參數(shù) [αi，υi，τi]之間的非線性關(guān)系，當(dāng)存在大量的信道參數(shù)時(shí)，很難用式(6)進(jìn)行分析．一般情況下，由于傳輸系統(tǒng)的帶寬及符號(hào)時(shí)間有限，因此，運(yùn)用采樣定理，信道響應(yīng)H(t，f)可以用離散參數(shù)來近似積分，即

由式(8)可以看出，時(shí)延域和多普勒域的最小采樣間隔分別為1/T和1/W，而L和K分別代表信道中存在的最大時(shí)延數(shù)和單邊頻移數(shù)．一個(gè)最小采樣間隔內(nèi)的時(shí)延分布或頻移分布可視為一個(gè)時(shí)延簇或多普勒頻移簇．式(9)中Lτ,l∈{i：τi∈[l/W－1/2,W，l/W＋1/2,W)}和Lυ,k∈{i：υi∈[k/T－1/2,T，k/T＋1/2,T)}分別代表第l個(gè)時(shí)延簇和第k個(gè)多普勒頻移簇的路徑集合．根據(jù)文獻(xiàn)[19]所提出的方法，可利用短時(shí)傅里葉(short time Fourier，STF)變換基來表征雙選擇性信道．由于 STF脈沖的持續(xù)時(shí)間和帶寬可以完全擬合雙選擇性信道的時(shí)延和多普勒頻移，因此，可將信道響應(yīng)完全轉(zhuǎn)化為離散形式．設(shè) STF基的最小時(shí)域和頻域間隔分別為T0和F0，則在時(shí)頻域內(nèi)分別有Nt＝T/T0個(gè)時(shí)域間隔和Nf＝W/F0個(gè)頻域間隔．設(shè)t＝nT0，f＝mF0，則有

式中：?代表矩陣的克羅內(nèi)克積運(yùn)算；H是L×(2,K+1)的信道系數(shù)矩陣；h＝vec(H)∈CN是信道系數(shù)向量，其元素是由H的每行依次銜接而成，即h＝[h0,-K，…，h0,K，h1,-K，…，h1,K，…，hL-1,-K，…，hL-1,K]，由于在雙選擇衰落條件下，信道呈現(xiàn)稀疏特性，因此，h為稀疏向量，即只有少數(shù)非零元素．

基于式(10)，由 Hn,m組成的一維信道響應(yīng)向量H′可以表示為

式中：H′含有 NtNf個(gè)元素；U為 NtNf×L(2,K＋1)維矩陣，其每一行由組成．對(duì)于基于導(dǎo)頻的信道估計(jì)，假設(shè)導(dǎo)頻插入位置(n，m)∈P，則在接收端可直接估計(jì)出在P處的信道響應(yīng)

式中：Rn,m和 Sn,m分別代表接收符號(hào)和導(dǎo)頻符號(hào)；Zn,m為高斯噪聲．設(shè)Hp是由導(dǎo)頻處信道響應(yīng)組成的H′的子向量，Up為所對(duì)應(yīng)的維子矩陣，Z為由組成的|P|維向量，則結(jié)合式(15)，式(14)可轉(zhuǎn)化為

通過式(16)恢復(fù)出信道系數(shù)向量h，則最終可通過式(14)得到信道響應(yīng)向量H′．

3 基于基追蹤的信道壓縮感知

3.1 傳統(tǒng)的二維插值算法

傳統(tǒng)的二維插值算法主要分為兩步：①先估計(jì)插入導(dǎo)頻符號(hào)處的信道響應(yīng)(見式(15))；②在①的基礎(chǔ)上，根據(jù)導(dǎo)頻位置的信道估計(jì)值進(jìn)行整個(gè)信道的二維內(nèi)插重構(gòu)．目前最常用的方法包括線性插值和FFT插值．線性插值算法每次估計(jì)只用兩個(gè)相鄰的導(dǎo)頻子信道的信道估計(jì)值，內(nèi)插得到兩個(gè)導(dǎo)頻之間的數(shù)據(jù)子載波的信道響應(yīng)．線性算法原理簡(jiǎn)單且易于實(shí)現(xiàn)，在實(shí)際應(yīng)用中非常有效，但估計(jì)精度比較粗略.

基于 FFT的信道估計(jì)算法是目前計(jì)算復(fù)雜度和性能最為均衡的選擇．該算法實(shí)現(xiàn)的主要思想是利用在時(shí)域補(bǔ)零等效于在頻域進(jìn)行內(nèi)插的原理，來恢復(fù)信道的頻率響應(yīng)．設(shè)導(dǎo)頻數(shù)量為Np，接收端估計(jì)出的導(dǎo)頻子信道處的信道響應(yīng)為 Hp(k)，k＝0，1，…，Np－1，將導(dǎo)頻子信道處的頻率響應(yīng)變換到時(shí)域的沖激響應(yīng) hp(n)，n＝0，1，…，Np－1，并在序列 hp(n)的尾部補(bǔ)零得到

將h(n)變換到頻域，則得到信道的頻域響應(yīng)

式中N0為符號(hào)總數(shù)，m＝0，1，…，N-1．

3.2 基于BP算法的稀疏信道估計(jì)

本文第 2節(jié)已逐步將雙選擇性衰落條件下的OFDM信道模型轉(zhuǎn)化為基于CS的信號(hào)重構(gòu)模型．令Hp＝y(tǒng)為線性測(cè)量值，Up＝Φ 為測(cè)量矩陣，并將 Up中的每一列歸一化，使其具有單位 l2范數(shù)．令 h＝x為待重構(gòu)的稀疏向量，則式(16)可轉(zhuǎn)化為式(1)中的CS重構(gòu)模型

式中Z是未知的噪聲向量．用BP算法求解考慮了重構(gòu)誤差及噪聲的線性規(guī)劃問題，即將式(19)轉(zhuǎn)化為BPIC模型

文獻(xiàn)[14]已證明式(20)與式(2)的解是等效的，因此應(yīng)用 BP算法即可求解出稀疏向量 x，即雙選擇性衰落條件下的信道系數(shù)向量 h．由于 BP算法采用內(nèi)點(diǎn)法求解 l1范數(shù)最小化問題，其重構(gòu)精度較高，且只需少量測(cè)量值(即導(dǎo)頻處的信道響應(yīng)值)即可重構(gòu)信道響應(yīng)，因而可以有效地減少導(dǎo)頻數(shù)量，實(shí)現(xiàn)頻譜利用率的大幅提升．然而，BP算法的算法復(fù)雜度也相對(duì)較高，設(shè)待重構(gòu)信號(hào)長度為 N，測(cè)量次數(shù)為 M，則算法復(fù)雜度為 O(23/2M N )．在本模型中測(cè)量矩陣U為NtNf×L(2,K＋1)維矩陣，因此，基于BP算法的信道估計(jì)復(fù)雜度為

4 仿真結(jié)果與分析

在仿真實(shí)驗(yàn)中，將基于壓縮感知的BP算法與傳統(tǒng)的FFT-Linear和FFT-FFT二維聯(lián)合插值算法進(jìn)行比較．采用英屬哥倫比亞大學(xué)(University of British Columbia，UBC)發(fā)布的 Matlab SGPL1[19]工具箱中的spg_bpdn(·)函數(shù)模擬 BP算法，并將重構(gòu)誤差ε設(shè)為1×10-3．FFT-Linear插值算法由兩個(gè)級(jí)聯(lián)的一維插值算法組成，即先在頻率軸上采用 FFT算法插值出信道響應(yīng)，再利用線性算法在時(shí)域軸上進(jìn)行插值，從而實(shí)現(xiàn)整個(gè)時(shí)頻域的信道估計(jì)．類似地，F(xiàn)FT-FFT插值算法則是在時(shí)頻域上皆采用 FFT插值算法，其性能要略好于 FFT-Linear算法，但同時(shí)也具有較高的計(jì)算復(fù)雜度．

假設(shè)OFDM系統(tǒng)子載波數(shù)為512，循環(huán)前綴長度為 128，系統(tǒng)的載波頻率為 5,GHz，且每個(gè)符號(hào)采用64QAM 調(diào)制．利用基于空間幾何結(jié)構(gòu)的信道仿真工具包IlmProp(IlmProp工具包是由德國伊爾默瑙工業(yè)大學(xué)電子系開發(fā)的功能強(qiáng)大的信道模擬工具包，能根據(jù)現(xiàn)實(shí)場(chǎng)景參數(shù)，如移動(dòng)臺(tái)速度、用戶場(chǎng)景(農(nóng)村或城市)等，產(chǎn)生與真實(shí)信道逼近的大規(guī)模系統(tǒng)級(jí)仿真模型)[20]產(chǎn)生模擬現(xiàn)實(shí)環(huán)境的雙選擇性衰落信道，所模擬的OFDM傳輸幀每幀有24個(gè)符號(hào)．在仿真環(huán)境中，設(shè)發(fā)送端和接收端之間距離 1,500,m，有 10個(gè)多徑路線并形成3個(gè)路徑簇(信道稀疏度約為3)，每個(gè)路徑的時(shí)延在(0，τmax)上隨機(jī)分布，每個(gè)路徑復(fù)增益服從復(fù)高斯分布，而且隨著路徑時(shí)延的增大，路徑復(fù)增益的功率以指數(shù)衰減．接收端的移動(dòng)速度為20,m/s，且與路徑簇具有獨(dú)立同分布的路徑導(dǎo)向．

本次仿真實(shí)驗(yàn)中，在不同信噪比(SNR)和導(dǎo)頻數(shù)下，對(duì)以上 3種算法的歸一化均方誤差(MSE)進(jìn)行了比較，其中導(dǎo)頻在傳輸符號(hào)中的插入方式服從均勻分布．圖1～圖3分別是導(dǎo)頻開銷為6.25%、12.5%和25%時(shí)，各算法的SNR-MSE曲線．

從圖中可以看出，無論導(dǎo)頻開銷是多少，BP算法的估計(jì)性能都要遠(yuǎn)遠(yuǎn)高于 FFT-Linear和 FFT-FFT插值算法．當(dāng)導(dǎo)頻開銷較大且信噪比較小時(shí)(見圖2(導(dǎo)頻開銷12.5%)和圖3(導(dǎo)頻開銷25%))，BP算法的優(yōu)越性并非十分明顯，但隨著信噪比的不斷增加，BP算法的MSE曲線不斷下降，而傳統(tǒng)的兩種算法的估計(jì)誤差則趨于一條水平的直線，因此 BP算法的性能優(yōu)勢(shì)越來越明顯．例如，當(dāng)導(dǎo)頻開銷為 12.5%、信噪比為40,dB時(shí)，BP算法的 MSE可達(dá)-42,dB，而傳統(tǒng)的插值算法仍只能維持在-16,dB左右．當(dāng)導(dǎo)頻開銷僅為6.25%(見圖1)時(shí)，BP算法在低信噪比時(shí)就表現(xiàn)出了算法性能的優(yōu)越性．當(dāng)信噪比僅為 5,dB時(shí)，傳統(tǒng)算法的MSE僅約為-13.5,dB，而BP算法則能達(dá)到-16.8,dB．隨著信噪比的增大，F(xiàn)FT-Linear和 FFTFFT算法的估計(jì)性能仍然沒有太大提升，而BP算法在信噪比為22,dB時(shí)，其MSE就能達(dá)到-21,dB．

圖1 BP與FFT-Linear、FFT-FFT算法MSE性能比較(導(dǎo)頻開銷6.25%)Fig.1 MSE performance comparison between BP and FFTLinear、FFT-FFT algorithms(Pilot cost is 6.25%)

圖2 BP與FFT-Linear、FFT-FFT算法MSE性能比較(導(dǎo)頻開銷12.5%)Fig.2 MSE performance comparison between BP and FFTLinear，F(xiàn)FT-FFT algorithms(Pilot cost is 12.5%)

圖3 BP與FFT-Linear、FFT-FFT算法MSE性能比較(導(dǎo)頻開銷25%)Fig.3 MSE performance comparison between BP and FFT-Linear，F(xiàn)FT-FFT algorithms(Pilot cost is 25%)

上述實(shí)驗(yàn)表明，當(dāng)物理信道存在強(qiáng)烈的雙選擇性衰落的情況下(即信道呈現(xiàn)稀疏特性)，傳統(tǒng)的 FFTLinear和 FFT-FFT二維插值算法已經(jīng)不能準(zhǔn)確地估計(jì)出信道響應(yīng)，而本文所提出的 BP算法則能利用極少的導(dǎo)頻(導(dǎo)頻開銷僅為 6.25%)實(shí)現(xiàn)對(duì)雙選擇性信道的準(zhǔn)確估計(jì)．因此，可見 BP算法可以極大地減少OFDM 系統(tǒng)中的導(dǎo)頻開銷，從而實(shí)現(xiàn)頻譜利用率的大幅提高．

5 結(jié) 語

本文提出了一種基于壓縮感知技術(shù)的OFDM系統(tǒng)信道估計(jì)方法．該方法針對(duì)雙選擇性衰落信道多為稀疏信道的特性，將OFDM的信道模型轉(zhuǎn)化為CS可解的重構(gòu)模型，并利用 BP算法實(shí)現(xiàn)了對(duì)稀疏信道的準(zhǔn)確估計(jì)．該法克服了傳統(tǒng)二維插值算法無法準(zhǔn)確估計(jì)稀疏信道的瓶頸，并可大大減少 OFDM 系統(tǒng)的導(dǎo)頻數(shù)量，提高系統(tǒng)的頻譜利用率．仿真實(shí)驗(yàn)表明，BP算法的信道估計(jì)性能要遠(yuǎn)高于 FFT-Linear和FFT-FFT二維插值算法，并能在導(dǎo)頻數(shù)量極其有限的情況下實(shí)現(xiàn)對(duì)雙選擇性稀疏信道的精確估計(jì)．

［1］ Raghavan V，Hariharan G，Sayeed A M. Capacity of sparse multipath channels in the ultra-wideband regime[J].IEEE Journal of Selected Topics in Signal Processing，2007，1(3)：357-371.

［2］ Yang B，Cao Z，Letaief K B. Analysis of lowcomplexity windowed DFT-based MMSE channel estimator for OFDM systems[J].IEEE Transactions on Communications，2001，49(11)：1977-1987.

［3］ Candès E J，Romberg J，Tao T. Robust uncertainty principles：Exact signal reconstruction from highly incomplete frequency information[J].IEEE Transactions on Information Theory，2006，52(2)：489-509.

［4］ Donoho D L. Compressed sensing[J].IEEE Transactions on Information Theory，2006，52(4)：1289-1306.

［5］ Baraniuk R G. Compressive sensing [J].IEEE Signal Processing Magazine，2007，24(4)：118-120.

［6］ Candès E J，Tao T. Decoding by linear programming[J].IEEE Transactions on Information Theory，2005，51(12)：4203-4215.

［7］ Cotter S F，Rao B D. Sparse channel estimation via matching pursuit with application to equalization[J].IEEE Transactions on Communications，2002，50(3)：374-377.

［8］ Bajwa W U，Haupt J，Raz G，et al. Compressed channel sensing[C] //Proceedings of Annual Conference on Information Sciences and Systems. Princeton，NJ，2008：5-10.

［9］ Wu C J，Lin D W. A group matching pursuit algorithm for sparse channel estimation for OFDM transmission[C]//International Conference on Acoustic，Speech and Signal Processing. Toulouse，F(xiàn)rance，2006：14-19.

［10］ Qi C H，Wu L N. Digital broadcast channel estimation with compressive sensing[J].Journal of Southeast University：English Edition，2010，26(3)：389-393.

［11］ Paredes J L，Arce G R，Wang Z. Ultra-wideband compressed sensing：Channel estimation[J].IEEE Journal of Selected Topics in Signal Processing，2007，1(3)：383-395.

［12］ Berger C R，Zhou S，Preisig J C，et al. Sparse channel estimation for multicarrier underwater acoustic communication：From subspace methods to compressive sensing[J].IEEE Transactions on Signal Processing，2010，58(3)：1708-1721.

［13］ Chen S B，Donoho D L，Saunders M A. Atomic decomposition by basis pursuit[J].SIAM Journal on Scientific Computing，1998，20(1)：33-61.

［14］ Donoho D L. For most large underdetermined systems of linear equations，the minimall1-norm solution is also the sparsest solution[J].Communications on Pure and Applied Mathematics，2006，59(6)：797-829.

［15］ Tropp J，Gilbert A. Signal recovery from random measurements via orthogonal matching pursuit[J].IEEE Transactions on Information Theory，2007，53(12)：4655-4666.

［16］ Needell D，Vershynin R. Uniform uncertainty principle and signal recovery via regularized orthogonal matching pursuit [J].Foundations of Computational Mathematics，2009，9：317-334.

［17］ Donoho D L，Tsaig Y，Drori I，et al. Sparse Solution of Underdetermined Linear Equations by Stagewise Orthogonal Matching Pursuit [EB/OL]. http：//www. dsp.ece. rice. edu/cs/，2011-11-20.

［18］ Herrity K K，Gilbert A C，Tropp J A. Sparse approximation via iterative thresholding[C] //International Conference on Acoustics，Speech and Signal Processing(ICASSP). Washington，USA，2006：624-627.

［19］ Friedlander M，van den Berg E. Toolbox SPGL1[EB/OL]. http：//www. cs. ubc. ca/labs/scl/spgl1，2011-11-20.

［20］ Ilmenau University of Technology. IlmProp Project[EB/OL]. http：//tu-ilmenau. de/ilmprop，2011-11-20.