吳君欽+李寧+王加莉

摘 要： 對(duì)于MIMO?OFDM系統(tǒng)信道估計(jì)的研究在國(guó)內(nèi)外已經(jīng)取得了一定的成果，其中采用CoSaMP算法進(jìn)行信道估計(jì)時(shí)支撐候選集的選擇尤為重要。采用CoSaMP算法對(duì)MIMO?OFDM信道進(jìn)行估計(jì)，并在此基礎(chǔ)上引入原子弱選擇標(biāo)準(zhǔn)對(duì)候選集進(jìn)行更進(jìn)一步的優(yōu)化，使改進(jìn)后的CoSaMP算法在達(dá)到自適應(yīng)的同時(shí)提高運(yùn)算速率，達(dá)到信道估計(jì)算法自適應(yīng)的目的。實(shí)驗(yàn)結(jié)果表明，對(duì)CoSaMP算法的候選集選擇添加閾值后進(jìn)行重構(gòu)信號(hào)時(shí)復(fù)雜度降低，提高了運(yùn)算速率，同時(shí)均方根誤差較CoSaMP算法和OMP算法都有一定的提高。

關(guān)鍵詞： 信道估計(jì)； 壓縮感知； MIMO?OFDM； 原子弱選擇； CoSaMP算法； OMP算法

中圖分類號(hào)： TN92?34 文獻(xiàn)標(biāo)識(shí)碼： A 文章編號(hào)： 1004?373X（2018）01?0009?04

Abstract： The study on channel estimation of MIMO?OFDM system has gotten a certain achievement at home and abroad， in which it is particularly important for channel estimation to use the CoSaMP algorithm to support candidate set selection. The CoSaMP algorithm is adopted for estimation of the MIMO?OFDM channel. And on this basis， the atomic weak selection criteria is introduced to optimize the candidate set further， which can improve the computing speed while making the improved CoSaMP algorithm realize self?adaptation. The experimental results show that the improved CoSaMP algorithm can reduce the complexity of signal reconstruction while the threshold is selected and added into the candidate set of the algorithm and improve the computing speed， and its root?mean?square error is better than that of the CoSaMP algorithm and OMP algorithm.

Keywords： channel estimation； compressed sensing； MIMO?OFDM； atomic weak selection； CoSaMP algorithm； OMP algorithm

0 引 言

正交頻分復(fù)用（OFDM）技術(shù)在無(wú)線通信傳輸中的運(yùn)用提高了無(wú)線傳輸效率。與多輸入多輸出（MIMO）技術(shù)相結(jié)合，能提高系統(tǒng)的通信質(zhì)量和系統(tǒng)容量[1]。在現(xiàn)有的無(wú)線通信環(huán)境中，無(wú)線頻譜資源相當(dāng)匱乏，MIMO?OFDM系統(tǒng)能有效提高資源利用率和對(duì)抗信號(hào)頻率的選擇性衰落。信號(hào)在傳輸過(guò)程中會(huì)產(chǎn)生符號(hào)間干擾（Inter Symbol Interference，ISI），因此，要對(duì)信道的狀態(tài)信息進(jìn)行估計(jì)，從而減少符號(hào)間的干擾，提高傳輸質(zhì)量，發(fā)揮其優(yōu)越性[2]。傳統(tǒng)的信道估計(jì)算法沒(méi)有充分考慮信號(hào)的特性，文獻(xiàn)[3]對(duì)信道頻域自相關(guān)矩陣采用矩陣低秩逼近得到信道響應(yīng)，該算法對(duì)精度有一定的提高，但是運(yùn)算量較大，存在一定的復(fù)雜度問(wèn)題，現(xiàn)實(shí)中運(yùn)用存在一定的難度。近年來(lái)實(shí)驗(yàn)研究表明，信號(hào)在傳輸過(guò)程中大多數(shù)能量主要集中在少數(shù)幾個(gè)信道路徑上，體現(xiàn)信號(hào)的傳輸具有稀疏性[4]。根據(jù)信號(hào)的這一先驗(yàn)知識(shí)，DL Donoho等人在2006年提出了壓縮感知（Compressed Sensing，CS）[5]技術(shù)，極大地提高了信道估計(jì)的速率和精確度。壓縮感知技術(shù)的提出和證明，隨即在生物傳感、圖像處理、信號(hào)處理等領(lǐng)域取得了一定的研究成果。文獻(xiàn)[6]采用匹配追蹤算法（MP），相比于傳統(tǒng)利用最小二乘估計(jì)算法（LS）信道估計(jì)復(fù)雜度降低，以較少的導(dǎo)頻達(dá)到較高的信道恢復(fù)效果。文獻(xiàn)[7]運(yùn)用正交匹配追蹤算法（OMP）對(duì)OFDM系統(tǒng)進(jìn)行信道估計(jì)，仿真結(jié)果顯示具有很好的效果，但是存在復(fù)雜度的問(wèn)題，并且需要稀疏度已知。

本文在匹配追蹤算法的基礎(chǔ)上進(jìn)行優(yōu)化，并通過(guò)系統(tǒng)仿真證明，利用CoSaMP算法對(duì)MIMO?OFDM系統(tǒng)進(jìn)行信道估計(jì)，在保證精確的基礎(chǔ)上，算法的運(yùn)算速率還有一定的提升空間。

1 MIMO?OFDM系統(tǒng)原理

假設(shè)MIMO?OFDM系統(tǒng)在發(fā)射端有[Nt]根天線，在接收端有[Nr]根天線，子載波個(gè)數(shù)為[K。]將發(fā)送的數(shù)據(jù)[X]分成[Nt]個(gè)數(shù)據(jù)塊[Xi，]其中[i=1，2，…，Nt，]原理圖如圖1所示。信道參數(shù)為選擇性衰落信道，并且該信道滿足壓縮感知的稀疏性，且稀疏信道的相干時(shí)間遠(yuǎn)大于OFDM的符號(hào)周期。則第[m]根發(fā)送天線與第[n]根接收天線之間的響應(yīng)為：

[hnm=i=0N-1hnmiδτ-τi] （1）

式中：[N]為信道的多徑數(shù)；[hnm]為第[i]徑的復(fù)增益；[τi]為第[i]徑的時(shí)延。其在子載波[k]處的頻域信道響應(yīng)為：

[Hnmk=t=0N-1hnmtWktK] （2）

式中[WktK=exp-2πjktK。]

對(duì)發(fā)射端的數(shù)據(jù)塊[Xi]分別進(jìn)行IFFT變換，然后插入循環(huán)前綴（Cyclic Prefix，CP），其中假設(shè)循環(huán)長(zhǎng)度要大于最大的路徑時(shí)延，這樣就可以忽略系統(tǒng)的載波間干擾和符號(hào)間干擾。最后進(jìn)行并/串轉(zhuǎn)換后通過(guò)天線進(jìn)行無(wú)線傳輸。在接收端，接收天線對(duì)數(shù)據(jù)進(jìn)行發(fā)送前的逆處理得到輸出信號(hào)。假設(shè)OFDM系統(tǒng)采取[N]點(diǎn)傅里葉變換，系統(tǒng)中有[P]個(gè)導(dǎo)頻信號(hào)分別位于子載波[k1，k2，…，kp]上，用于信道估計(jì)。則在第[n]根接收天線接收到的信號(hào)表示為：endprint

[ynk=m=1NtHnmkXmk+Znk] （3）

式中：[ynk]和[Xmk]為第[n]根接收天線和第[m]根發(fā)送天線在第[k]個(gè)子載波上接收和發(fā)送的信號(hào)；[Hnmk]為第[m]根天線與第[n]根天線在[k]上的頻域沖擊響應(yīng)。

由式（3）可知，第[n]根天線上接收的導(dǎo)頻信號(hào)表示為：

[yn=m=1NtdiagXmHnm+Zn] （4）

式中：[yn=ynk1，ynk2，…，ynkpT]為第[n]根天線接收到的導(dǎo)頻信號(hào)；[Xm=Xmk1，Xmk2，…，Xmkp]為第[m]根天線發(fā)送的導(dǎo)頻信號(hào)；[Zn]是均值為零，方差為[σ2n]的高斯白噪聲；[Hnm=Hnmk1，Hnmk2，…，Hnmkp]代表第[m]根天線與第[n]根天線間的信道頻域沖擊響應(yīng)。令[F]為傅里葉變換矩陣：

[F=1KW00KW01K…W0K-1KW10KW11K…W1K-1K????WK-10KWK-11K…WK-1K-1K] （5）

則式（4）可以表示為：

[yn=m=1NtdiagXmFLhnm+Zn] （6）

令[hn=hTn1，hTn2，…，hTNNtT]為[NtN×1]的列向量；[X=diagx1FN，diagx2FN，…，diagxNtFN]為[P×NtN]維矩陣，則式（6）為：

[yn=Xhn+Zn] （7）

式（7）為MIMO?OFDM系統(tǒng)模型，可以表示為：

[y=Xh+Z] （8）

2 壓縮感知

2.1 壓縮感知綜述

壓縮感知的提出突破了信號(hào)處理對(duì)于奈奎斯特抽樣定理的限定，其原理是利用非線性優(yōu)化從線性觀測(cè)中恢復(fù)逼近出原始信號(hào)。該方法對(duì)稀疏信號(hào)能以較低的采樣頻率采樣，并在接收端轉(zhuǎn)化為求解最優(yōu)化的問(wèn)題，能達(dá)到壓縮的目的，從而提高資源利用率。由式（8）知壓縮感知觀測(cè)恢復(fù)模型為：[y=Xh+Z。]式中，[X]為[M×N]維測(cè)量矩陣（[M?N]），[Z]為隨機(jī)高斯白噪聲，其中假設(shè)一個(gè)有限長(zhǎng)信號(hào)[h，][h∈RN×1，][h]在某個(gè)變換域[ψ]上是稀疏的，信號(hào)[h]可表示為：

[h=i=1Nθiψi 或 h=ψθ] （9）

并且這里要求[X]滿足“限制等距特性（RIP）”[8]，測(cè)量矩陣[X]與[K]稀疏信號(hào)[h]和常數(shù)[δK∈（0，1）]滿足[K]階約束等距性，即：

[1-δKh22≤Xh22≤1+δKh22] （10）

式中[h22=i=1Nhi2]表示矢量信號(hào)[h]的2階范數(shù)。

本文采取隨機(jī)高斯矩陣作為觀測(cè)矩陣[Φ，]理論證明當(dāng)[M≥cKlogNK]時(shí)高斯矩陣在很大概率下具有RIP性質(zhì)[9]，且隨機(jī)高斯矩陣與大多數(shù)矩陣具有很好的不相關(guān)性。對(duì)于接收端接收到的觀測(cè)值[y，]通過(guò)轉(zhuǎn)化為求解一個(gè)最優(yōu)[l0]范數(shù)問(wèn)題來(lái)重構(gòu)原始信號(hào)。

[minψTX0s.t. ACSX≈?ψTX=Y] （11）

本文采用貪婪算法對(duì)觀測(cè)信號(hào)進(jìn)行估計(jì)，通過(guò)每次迭代過(guò)程選出一個(gè)或多個(gè)最優(yōu)解來(lái)逐步逼近，達(dá)到最大限度恢復(fù)出原始信號(hào)的目的。這類算法包括MP算法、OMP算法、CoSaMP算法等，相比于其他類恢復(fù)算法，具有實(shí)時(shí)性和重建信號(hào)的維數(shù)較小等優(yōu)點(diǎn)。

2.2 稀疏信道估計(jì)算法

CoSaMP算法在進(jìn)行信道估計(jì)時(shí)需要已知信號(hào)的稀疏度，并且復(fù)雜度較高。本文基于壓縮感知匹配追蹤算法來(lái)降低復(fù)雜度和稀疏度自適應(yīng)改進(jìn)，并在此基礎(chǔ)上使用原子弱選擇標(biāo)準(zhǔn)[10]進(jìn)行篩選，使候選集的選擇更加精確，使算法的復(fù)雜度有所降低。

CoSaMP算法如下：

輸入：觀測(cè)值[y，]恢復(fù)矩陣[A，]稀疏度[K]

輸出：信號(hào)的稀疏逼近

① 初始化[r0=y，][i=1，][Ω0=?；]

② 求殘差向量與恢復(fù)矩陣的內(nèi)積，[θ=ATri-1；]

③ 內(nèi)積的前[2K]為索引，[P=supθ2K]；

④ 更新索引集，[P=P?Ωi-1；]

⑤ 最小二乘法，[θP=A+：，Py；]

⑥ 選取前[K]個(gè)值，[θ=θK；]

⑦ 更新殘差向量，[ri=y-Aθ。]

從傳統(tǒng)的CoSaMP算法中可以看出，該算法每次在對(duì)索引進(jìn)行選擇時(shí)都為[2K，]而在進(jìn)一步的選擇時(shí)為[K，]因而這些下標(biāo)存在冗余，并且隨著殘差能量的減少，冗余效果更加明顯，因此，該算法的預(yù)選方案會(huì)使算法的復(fù)雜度提高。另外，傳統(tǒng)CoSaMP算法在稀疏度已知情況下進(jìn)行信道估計(jì)，限制了該算法的實(shí)際運(yùn)用。本文在不影響算法效果的同時(shí)，通過(guò)對(duì)該算法預(yù)選階段的改進(jìn)對(duì)該算法進(jìn)行優(yōu)化和提升，并實(shí)現(xiàn)稀疏度自適應(yīng)，以達(dá)到提高算法運(yùn)算速率的目的。

改進(jìn)后的算法在CoSaMP算法基礎(chǔ)上采用自適應(yīng)恢復(fù)算法，并在預(yù)選階段采用原子弱選擇標(biāo)準(zhǔn)進(jìn)行篩選，降低算法的復(fù)雜度，原子弱選擇標(biāo)準(zhǔn)如下：

[i：gj≥α?maxjgj] （12）

式中：[α∈（0，1]]；[gjgj=j=1，2，…，N]為觀測(cè)矩陣[A]與殘差矩陣的內(nèi)積。

在求得殘差與每一列內(nèi)積后的索引值時(shí)，通過(guò)式（12）進(jìn)行進(jìn)一步篩選，減少索引的冗余，降低算法的復(fù)雜度。本算法的優(yōu)點(diǎn)實(shí)現(xiàn)了算法自適應(yīng)，并在此基礎(chǔ)上提高了算法的運(yùn)算速率。

CoSaMP算法閾值改進(jìn)后的流程如下：

輸入：觀測(cè)值[y，]恢復(fù)矩陣[A，]初始化支撐集尺寸[step_size]

輸出：信號(hào)估計(jì)值

1） 算法初始化。

初始化殘差[r0=y；]初始化支撐集[S=step_size；]索引集[Ω0=?；]迭代次數(shù)[i=1；]階段次數(shù)[n=1。]

2） 循環(huán)迭代過(guò)程。

① 計(jì)算殘差[rn-1]與矩陣[A]的每一列的內(nèi)積，[u=ATrn-1，]并將[u]的絕對(duì)值中[2S]個(gè)最大值對(duì)應(yīng)的索引存入[P]中。

② 原子弱選擇：從[2S]個(gè)原子中選取所有滿足[i：gj≥α?maxjgj]的腳標(biāo)存入候選集[Jn-1]中。

③ 更新候選集[Ωn-1=Ωn?Jn-1]。

④ 獲取索引內(nèi)新的估計(jì)值[hn=X+Ωny]。

⑤ 選取最大的[S]個(gè)主要抽頭系數(shù)，其他非主要抽頭系數(shù)置0，[ΩnD=suphnS]。

⑥ 更新殘差向量。若新的殘差能量大于上次殘差能量，[n=n+1]，[S=S×n]。

⑦ 循環(huán)①～⑥，直到滿足迭代停止條件。其中，本文的停止條件為[hn-h22≤10-3S≥64]。

改進(jìn)后的CoSaMP算法在得到原子候選集后，通過(guò)原子弱選擇的閾值選擇優(yōu)化準(zhǔn)則后，降低了原算法在預(yù)選階段的冗余，在實(shí)際應(yīng)用中原子的預(yù)選階段對(duì)算法精確度的影響不是很大，而且提高了算法的復(fù)雜度，因此，在預(yù)選階段對(duì)原子的選擇引用原子弱選擇標(biāo)準(zhǔn)對(duì)CoSaMP算法進(jìn)行閾值改進(jìn)，在不影響精確度的同時(shí)能夠降低算法的復(fù)雜度。

3 實(shí)驗(yàn)分析

MIMO?OFDM系統(tǒng)信道為選擇性衰落信道，采用Simulink繪制，其參數(shù)為：[Nt=2，][Nr=2，]調(diào)制方式為QAM調(diào)制，子載波個(gè)數(shù)[N=1 024，]其中導(dǎo)頻子載波個(gè)數(shù)為56，信道長(zhǎng)度[L=32，]其中原子弱選擇標(biāo)準(zhǔn)中的[α]取0.5，假設(shè)系統(tǒng)的參數(shù)在一個(gè)OFDM符號(hào)周期內(nèi)不變化。本文采用均方根誤差來(lái)反應(yīng)信道估計(jì)誤差，[αRMSE]越小則說(shuō)明估計(jì)誤差越小。均方根誤差表示為：

[αRMSE=1Mm=1Mh-hm22] （13）

式中[M]為蒙特卡洛仿真次數(shù)。

仿真結(jié)果驗(yàn)證了CoSaMP算法在引入原子弱選擇標(biāo)準(zhǔn)進(jìn)行閾值改進(jìn)，在不影響復(fù)雜度的前提下降低了算法的復(fù)雜度，提高了算法的運(yùn)算速率。圖2為不同信噪比下的LS算法、OMP算法、CoSaMP算法及改進(jìn)后的CoSaMP算法的效果仿真圖。

在MIMO?OFDM系統(tǒng)中，采用OMP算法、CoSaMP算法及改進(jìn)后的CoSaMP算法較傳統(tǒng)的LS算法的均方根誤差值有很大的提高，并且改進(jìn)后的CoSaMP算法與傳統(tǒng)的CoSaMP算法的均方根誤差相似，說(shuō)明CoSaMP算法在通過(guò)原子弱選擇標(biāo)準(zhǔn)進(jìn)行閾值改進(jìn)后對(duì)算法的精確度沒(méi)有影響。

由于算法的復(fù)雜度決定了運(yùn)算速率，對(duì)算法的實(shí)際應(yīng)用具有很大的影響，因此算法的復(fù)雜度在進(jìn)行信道估計(jì)應(yīng)用中至關(guān)重要。圖3是各算法復(fù)雜度與LS算法比值的仿真圖，本文仿真硬件參數(shù)為Intel i5?3210M，CPU為2.50 GHz，RAM為4.00 GB，Microsoft Windows 7操作系統(tǒng)。本文中用運(yùn)算時(shí)間對(duì)算法的復(fù)雜度進(jìn)行近似估計(jì)，即表示為：

[R1=O（OMP）O（LS）R2=O（CoSaMP）O（LS）R3=O（改進(jìn)CoSaMP）O（LS）]

式中：[O（?）]表示復(fù)雜度符號(hào)；[R1，][R2，][R3]為比值。

從圖3中可以看出，改進(jìn)后的CoSaMP算法在復(fù)雜度上較CoSaMP算法有很大的提高。其中改進(jìn)后的算法在運(yùn)算時(shí)間上大概是LS算法的2倍，而CoSaMP算法是LS算法的10倍左右。結(jié)合圖2、圖3可以證明，改進(jìn)后的算法在不影響精確度的情況下提高了算法的運(yùn)算速率，減少了運(yùn)算時(shí)間。

4 結(jié) 語(yǔ)

本文在MIMO?OFDM信道估計(jì)研究成果的基礎(chǔ)上對(duì)CoSaMP算法進(jìn)行閾值改進(jìn)，并對(duì)LS，OMP，CoSaMP算法進(jìn)行仿真比較，結(jié)果表明改進(jìn)后的CoSaMP算法在運(yùn)算復(fù)雜度上有較大降低，提高了算法的運(yùn)算速率。進(jìn)一步研究表明，改進(jìn)后的CoSaMP算法在降低復(fù)雜度的前提下，算法的均方誤差的性能還有提升的空間，這也是以后研究的目標(biāo)。

參考文獻(xiàn)

[1] STUBER G L， BARRY J R， MCLAUGHLIN S W， et al. Broadband MIMO?OFDM wireless communications [J]. Proceedings of the IEEE， 2004， 92（2）： 271?294.

[2] 朱立君，周圍.MIMO?OFDM系統(tǒng)中的信道估計(jì)算法綜述[J].鄭州輕工業(yè)學(xué)院學(xué)報(bào)（自然科學(xué)版），2009，24（2）：87?90.

ZHU Lijun， ZHOU Wei. Overview of channel estimation algorithms in MIMO?OFDM systems [J]. Journal of Zhengzhou University of Light Industry （natural science）， 2009， 24（2）： 87?90.

[3] WANG Z J， HAN Z， LIU K J R. A MIMO?OFDM channel estimation approach using time of arrivals [J]. IEEE transactions on wireless communications， 2005， 4（3）： 1207?1213.

[4] LI W， PREISIG J C. Estimation of rapidly time?varying sparse channels [J]. IEEE journal of oceanic engineering， 2007， 32（4）： 927?939.

[5] DONOHO D L. Compressed sensing [J]. IEEE transactions on information theory， 2006， 52（4）： 1289?1306.

[6] COTTER S F， RAO B D. Sparse channel estimation via mat?ching pursuit with application to equalization [J]. IEEE transactions on communications， 2002， 50（3）： 374?377.

[7] 何雪云，宋榮方，周克琴.基于壓縮感知的OFDM系統(tǒng)稀疏信道估計(jì)新方法研究[J].南京郵電大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2010，30（2）：60?65.

HE Xueyun， SONG Rongfang， ZHOU Keqin. A novel method for sparse channel estimation in OFDM systems based on compressed sensing [J]. Journal of Nanjing University of Posts and Telecommunications （natural science）， 2010， 30（2）： 60?65.

[8] BAJWA W U， SAYEED A， NOWAK R. Sparse multipath channels： modeling and estimation [C]// Proceedings of the 13th Digital Signal Processing Workshop. Marco Island： IEEE， 2009： 320?325.

[9] BARANIUK R G. A lecture on compressive sensing [J]. IEEE signal processing magazine， 2007（7）： 1?9.

[10] BLUMENSATH T， DAVIES M E. Stagewise weak gradient pursuits [J]. IEEE transactions on signal processing， 2009， 57（11）： 4333?4346.endprint