用于寬帶頻譜感知的全盲亞奈奎斯特采樣方法

蓋建新 付 平 喬家慶 孟升衛(wèi)③

①(哈爾濱工業(yè)大學(xué)自動(dòng)化測(cè)試與控制系 哈爾濱 150080)

②(哈爾濱理工大學(xué)測(cè)控技術(shù)與儀器黑龍江省高校重點(diǎn)實(shí)驗(yàn)室 哈爾濱 150080)

③(中國(guó)科學(xué)院電子學(xué)研究所 北京 100190)

1 引言

認(rèn)知無(wú)線電通過(guò)感知周圍頻譜環(huán)境自主發(fā)現(xiàn)“頻譜空穴”并對(duì)其進(jìn)行有效利用，在解決無(wú)線通信中頻譜資源緊張、頻譜利用率低等問(wèn)題上表現(xiàn)出巨大的優(yōu)勢(shì)。寬帶頻譜感知技術(shù)可以在較短的時(shí)間內(nèi)為認(rèn)知無(wú)線電提供更多的頻譜接入機(jī)會(huì)[1-3]，以壓縮感知(CS)理論[4,5]為基礎(chǔ)的寬帶頻譜感知技術(shù)以其采樣率低并可精確重構(gòu)等特點(diǎn)受到了廣泛關(guān)注[6-8]。這些方法以離散化頻譜的稀疏性為前提，實(shí)現(xiàn)了基于亞奈奎斯特采樣的寬帶頻譜感知，緩解了采樣率高的壓力。然而，當(dāng)需要實(shí)現(xiàn)較高的頻率分辨率時(shí)，由于測(cè)量矩陣維數(shù)過(guò)大，數(shù)據(jù)處理負(fù)擔(dān)繁重，導(dǎo)致感知速度較慢；另外，在頻譜重構(gòu)階段需要確切的稀疏度信息，而稀疏度在實(shí)際的頻譜環(huán)境中是無(wú)法準(zhǔn)確預(yù)知的。因此找到一種測(cè)量矩陣較小，不需要準(zhǔn)確稀疏度的亞奈奎斯特采樣方法，是寬帶頻譜感知技術(shù)當(dāng)前急待解決的問(wèn)題。

近年來(lái)出現(xiàn)的調(diào)制寬帶轉(zhuǎn)換器(MWC)采樣方法[9]，不需對(duì)頻譜進(jìn)行離散化，利用低維的測(cè)量矩陣即可實(shí)現(xiàn)對(duì)多帶信號(hào)的亞奈奎斯特采樣和精確重構(gòu)。若采用MWC作為采樣前端可以緩解CS方法的計(jì)算負(fù)擔(dān)。但現(xiàn)有的MWC理論須已知頻帶數(shù)量和最大帶寬才能構(gòu)建MWC系統(tǒng)并進(jìn)行頻譜重構(gòu)。實(shí)際的無(wú)線頻譜中，特定時(shí)間內(nèi)占用的頻帶數(shù)量是無(wú)法準(zhǔn)確預(yù)知的，而且由于相鄰信道可能同時(shí)被占用，導(dǎo)致最大頻帶寬度也是未知的。針對(duì)無(wú)線頻譜上述的不可預(yù)知性特點(diǎn)，本文將無(wú)線電發(fā)射信道模型和頻譜多帶模型相結(jié)合，對(duì)MWC適用的信號(hào)模型進(jìn)行了重新定義，在此基礎(chǔ)上提出了一個(gè)不需最大頻帶寬度和確切頻帶數(shù)量的重構(gòu)充分條件。在重構(gòu)算法方面，將稀疏度自適應(yīng)匹配追蹤算法進(jìn)行推廣并應(yīng)用到頻譜重構(gòu)中以消除對(duì)頻帶數(shù)量的依賴性。最終實(shí)現(xiàn)了既不需要各頻帶寬度，也不需要準(zhǔn)確頻帶數(shù)量的寬帶頻譜全盲亞奈奎斯特采樣，數(shù)值實(shí)驗(yàn)驗(yàn)證了該方法的有效性。

2 MWC采樣理論

2.1 采樣原理

MWC是一種多帶信號(hào)亞奈奎斯特采樣方法，其具體采樣原理如下。如圖1所示，輸入信號(hào)x(t)同時(shí)進(jìn)入m個(gè)通道，在第i個(gè)通道被周期為Tp(頻率為fp)的偽隨機(jī)符號(hào)序列pi(t)混頻，混頻后采用截止頻率為1/2Ts的理想低通濾波器h(t)進(jìn)行濾波，最后通過(guò)采樣率fs=1/Ts的ADC獲得m組低速采樣yi(n)。由經(jīng)典傅里葉分析思想可以推導(dǎo)出第i個(gè)通道輸出序列yi(n)的離散時(shí)間傅里葉變換(DTFT)與x(t)的傅里葉變換X(f)之間有如下關(guān)系：

其中cin表示序列pi(t)傅里葉級(jí)數(shù)的系數(shù)，L0=

圖1 MWC采樣系統(tǒng)框圖

其中Φ是m×L矩陣，Φin=cin且m＜L。若對(duì)式(2)兩端同時(shí)進(jìn)行DTFT的逆運(yùn)算，可得未知序列Ζ(n)與測(cè)量值Y(n)之間的線性關(guān)系：

其中Y(n)= [y1(n),y2(n),…,ym(n)]T,Z(n)=[z1(n),z2(n),…,zL(n)]T,n∈Z。由于m＜L，式(2)和式(3)均是欠定的，無(wú)法通過(guò)求逆的方法獲得唯一解?？紤]到多帶信號(hào)在頻域的稀疏性，z(f)中只有少量的非零元素，當(dāng)滿足如下定理給出的充分條件時(shí)，式(2)具有唯一的最稀疏解。

定理1[9]設(shè)多帶信號(hào)x(t)由N個(gè)頻帶組成，各頻帶中最大帶寬為B，按照?qǐng)D1所示的MWC結(jié)構(gòu)進(jìn)行采樣，如果以下條件成立：

(1)fs≥fp≥B并且fs/fp數(shù)值不是很大；

(2)一個(gè)周期內(nèi)序列pi(t)的符號(hào)(±1)間隔數(shù)M

(3)m≥2N；

(4)矩陣Φ的任意2N列線性無(wú)關(guān)。

則對(duì)于 ?f∈Fs,z(f)是式(2)的唯一的N-稀疏解。

2.2 重構(gòu)方法

當(dāng)多帶信號(hào)中各個(gè)頻帶的位置已知時(shí)，聯(lián)合支撐集Ω= s upp(Z(n))是確定的，如果矩陣ΦΩ滿足列滿秩則可以通過(guò)式(4)從采樣值Y(n)中恢復(fù)出Z(n)：

其中ΦΩ表示以Ω中的元素為索引的Φ的列子集，代表矩陣ΦΩ的偽逆矩陣且。當(dāng)Ω未知時(shí)需要將式(2)無(wú)限測(cè)量向量(IMV)問(wèn)題[10]變換成與之具有相同支撐的多測(cè)量向量(MMV)問(wèn)題[11]，然后，求解該MMV問(wèn)題的支撐集Ω，最后按式(4)完成重構(gòu)。

2.3 利用MWC理論進(jìn)行寬帶頻譜感知需解決的問(wèn)題

由定理1可知，構(gòu)建MWC時(shí)，采樣率fs及符號(hào)序列的頻率fp選擇的依據(jù)是最大頻帶寬度B。而在無(wú)線頻譜環(huán)境中各個(gè)信道是可以相鄰的，由于各信道占用的隨機(jī)性，很難預(yù)測(cè)不相鄰頻帶(每個(gè)頻帶可能由一個(gè)或多個(gè)占用的相鄰信道組成)的最大帶寬。另一方面，信道相鄰導(dǎo)致無(wú)法預(yù)知特定時(shí)間互不相鄰的頻帶數(shù)量N，而重構(gòu)階段需要以頻帶數(shù)(稀疏度)作為迭代的停止條件，才能通過(guò)貪婪算法(如同步正交匹配追蹤)找到支撐集。針對(duì)上述問(wèn)題，本文對(duì)認(rèn)知頻譜環(huán)境中信號(hào)模型進(jìn)行了重新定義，并改進(jìn)了頻譜重構(gòu)充分條件。

3 頻譜稀疏信號(hào)全盲亞奈奎斯特采樣

3.1 頻譜稀疏信號(hào)模型

假設(shè)信號(hào)x(t)為帶限于F=[-fNYQ/ 2,fNYQ/2]的實(shí)值連續(xù)時(shí)間信號(hào)，其中F由彼此相鄰但不重疊的多個(gè)無(wú)線發(fā)射信道組成。如果信號(hào)x(t)的頻譜X(f)由若干個(gè)帶寬為Bi彼此不相鄰的非零頻帶組成，每個(gè)頻帶包含一個(gè)或多個(gè)相鄰發(fā)射信道，且頻帶可以任意地分布在F內(nèi)，則稱這樣的信號(hào)集M為頻譜稀疏信號(hào)。圖2給出了M中典型的信號(hào)頻譜。

圖2 典型頻譜稀疏信號(hào)的頻譜示意圖

3.2 盲譜重構(gòu)充分條件

頻譜稀疏信號(hào)模型弱化了傳統(tǒng)多帶信號(hào)模型中最大頻帶寬度B=m ax(Bi)和頻帶數(shù)量這兩個(gè)參數(shù)的作用，將相鄰子頻帶(信道)看成了一個(gè)頻帶。為了實(shí)現(xiàn)適合上述模型的亞奈奎斯特采樣，本文將聯(lián)合稀疏度引入重構(gòu)條件中。

設(shè)參數(shù)化向量U(f)對(duì)于?f∈Fs支撐集supp(U(f))= {k|Uk(f)≠ 0}，式中Uk表示向量U的第k個(gè)元素，則U(Fs)的聯(lián)合支撐集可表示為：

改進(jìn)定理假設(shè)x(t)是M中的一個(gè)任意信號(hào)，采用MWC系統(tǒng)進(jìn)行采樣。對(duì)于式(2)，如果以下條件成立：

(1)fs≥fp，并且fs/fp數(shù)值不是很大；

(2)一個(gè)周期內(nèi)pi(t)的符號(hào)數(shù)

(3)m≥2K，其中；

(4)矩陣Φ的任意2K列線性無(wú)關(guān)。則z(Fs)是式(2)唯一的K-稀疏解。

證明由式(1)，式(2)可知，z(f)的每個(gè)分量都是X(f)的nfp移位、濾波所得到的寬度為fs的頻譜片段，其中fp是每次移位的步進(jìn)值。若fp≤fs成立，則每次移位的步進(jìn)足夠小，以至于z(Fs)中的頻譜片段可以覆蓋整個(gè)頻譜。因此如果式(2)能夠重構(gòu)，則從重構(gòu)的未知向量z(Fs)可以精確地恢復(fù)出整個(gè)頻譜X(f)及原信號(hào)x(t)。

條件(2)是為保證符號(hào)序列pi(t)在一個(gè)周期內(nèi)有足夠的偽隨機(jī)符號(hào)數(shù)，進(jìn)而使所構(gòu)成的矩陣有較大的Kruskal秩，具體證明參見文獻(xiàn)[9]。

考慮式(2)與式(3)的等價(jià)性，只要式(3)具有唯一K-稀疏解，則對(duì)應(yīng)可得到式(2)的K-稀疏解。易知向量組Y(n)所張成的空間的維數(shù)dim(span(Y(n)))≤K，因此可以找到空間span(Y(n))的一個(gè)框架V(n)，如果用V(n)代替式(3)中Y(n)則V(n)=ΦZ(n)與式(3)的唯一K-稀疏解具有相同的支撐集[10]。依據(jù)文獻(xiàn)[12]中定理 2.2，有以下等價(jià)命題成立：對(duì)于滿足V=ΦZ的聯(lián)合K-稀疏的向量Z，如果Φ的任意2K列均線性無(wú)關(guān)，則Z為此方程唯一的K-稀疏解。于是通過(guò)此方程可以求得式(3)的支撐集，進(jìn)而通過(guò)偽逆得到式(3)的唯一K-稀疏解。 證畢

本文所提出的改進(jìn)定理相當(dāng)于對(duì)定理 1的推廣。定理 1通過(guò)要求fp≥B= m ax(Bi)將稀疏度與頻帶數(shù)量聯(lián)系起來(lái)。由上述分析可知該條件不是必要的，當(dāng)fp＜ m ax(Bi)時(shí)在改進(jìn)定理?xiàng)l件下問(wèn)題式(2)仍然可以精確重構(gòu)。可見改進(jìn)定理的提出使 MWC系統(tǒng)的構(gòu)建擺脫了最大頻帶寬度的束縛。

3.3 頻譜感知中MWC通道數(shù)的確定

本文所提出的改進(jìn)定理使構(gòu)建MWC時(shí)不再需要頻帶數(shù)量和最大頻帶寬度，但需要根據(jù)實(shí)際頻譜對(duì)聯(lián)合稀疏度進(jìn)行估算，進(jìn)而估計(jì)出所需的通道數(shù)m。

考慮fs=fp的基本配置情況，設(shè)在實(shí)際頻譜中第i個(gè)頻帶Bi貢獻(xiàn)到向量z(f)的非零元素?cái)?shù)量記作I(Bi)。根據(jù)MWC采樣的頻域行為并考慮頻帶位置的隨機(jī)性可得：所以整個(gè)頻譜貢獻(xiàn)到z(f)中的非零元素?cái)?shù)即聯(lián)合稀疏度滿足如下關(guān)系：

應(yīng)該指出，MWC通道數(shù)的確定并不需要準(zhǔn)確的稀疏度K，實(shí)際中可以根據(jù)頻譜監(jiān)管部門頻譜統(tǒng)計(jì)特征確定一個(gè)保守的范圍(m≥ 2K)即可。

4 稀疏度自適應(yīng)重構(gòu)算法

在頻譜重構(gòu)的過(guò)程中如果不知道準(zhǔn)確的稀疏度則容易出現(xiàn)稀疏度過(guò)估計(jì)或欠估計(jì)的問(wèn)題，導(dǎo)致重構(gòu)結(jié)果具有嚴(yán)重誤差[13]。本文在原有MWC重構(gòu)框架基礎(chǔ)上將CS中的稀疏度自適應(yīng)匹配追蹤(SAMP)算法引入到了無(wú)限測(cè)量向量(IMV)問(wèn)題的求解中，在預(yù)先不知道確切稀疏度的前提下實(shí)現(xiàn)了支撐集的精確重構(gòu)。圖3給出了重構(gòu)的總體框圖。

圖3 盲譜重構(gòu)支撐集獲取框圖

如圖3所示，重構(gòu)時(shí)首先將IMV問(wèn)題轉(zhuǎn)化成具有相同支撐集的MMV問(wèn)題。為此需先找到Y(jié)(n)所張成的空間的一個(gè)框架V(n)，具體求解方法如下：計(jì)算，并利用奇異值分解原理找到滿足Fm×m=V?VT的矩陣V，則矩陣V即是所求框架[9]。然后將V作為測(cè)量值矩陣構(gòu)建MMV問(wèn)題。最后使用所提出的SAMP MMV算法求解該MMV問(wèn)題的支撐集。

4.1 SAMP MMV算法描述

SAMP算法充分利用了“自底向上”和“自頂向下”重構(gòu)思想的優(yōu)點(diǎn)，可以在未知稀疏度的前提下實(shí)現(xiàn)精確重構(gòu)[13]。本文把適用于CS的SAMP算法推廣到了MMV問(wèn)題的求解中，改進(jìn)了SAMP算法中的初步測(cè)試和最終測(cè)試思想。改進(jìn)后，用初步測(cè)試找出Φ中與Rk-1最相關(guān)的I個(gè)列的索引值Sk，最終測(cè)試將本次短支撐集Sk與上次最終集Fk-1進(jìn)行合并形成2I個(gè)元素的候選集Ck，然后通過(guò)回溯思想從Ck中選出I個(gè)最優(yōu)的向量作為本次迭代的最終集Fk。整個(gè)過(guò)程分為若干個(gè)階段，且每個(gè)階段均有若干次迭代，具體流程如表1所示。

4.2 SAMP MMV算法有效性分析

SAMP MMV算法采用回溯的思想使每個(gè)階段都具有如下特點(diǎn)，即當(dāng)前迭代的余項(xiàng)能量總比上一次迭代的余項(xiàng)能量小[13]。步驟 6的階段轉(zhuǎn)換條件保證了迭代中余項(xiàng)能量的單調(diào)遞減性，使算法在每個(gè)階段均向收斂的方向發(fā)展。當(dāng)算法執(zhí)行到階段j時(shí)，最終集F的元素?cái)?shù)量為js，可選的元素總數(shù)為L(zhǎng)，因此F最多有種可能形式。如果迭代次數(shù)超過(guò)，則最終集開始重復(fù)，余項(xiàng)能量不再單調(diào)遞減，算法轉(zhuǎn)入下一階段。可見該算法不會(huì)在任何階段陷入無(wú)限次迭代。當(dāng)算法執(zhí)行到第階段時(shí)便等效成估計(jì)稀疏度為的子空間追蹤算法，其收斂性在文獻(xiàn)[14]中已得到證明，故該算法在有限迭代次數(shù)內(nèi)一定收斂。算法收斂時(shí)的稀疏度與真實(shí)的稀疏度K近似相等，從而在未知稀疏度的前提下實(shí)現(xiàn)了重構(gòu)。當(dāng)然重構(gòu)中算法需經(jīng)過(guò)個(gè)階段的嘗試，使計(jì)算復(fù)雜度略有增加。

表1 SAMP MMV算法流程

5 仿真實(shí)驗(yàn)及結(jié)果分析

為了驗(yàn)證本文方法的有效性，本節(jié)設(shè)計(jì)了 3個(gè)仿真實(shí)驗(yàn)：5.1節(jié)驗(yàn)證在定理1條件下SAMP MMV算法的有效性；5.2節(jié)在5.1節(jié)的基礎(chǔ)上驗(yàn)證在改進(jìn)定理?xiàng)l件下SAMP MMV算法的有效性；5.3節(jié)結(jié)合應(yīng)用背景考察該方法的重構(gòu)效果。實(shí)驗(yàn)中頻譜稀疏信號(hào)均采用下式產(chǎn)生：

其中En,Bn,fn,τn分別代表所產(chǎn)生的第n個(gè)頻帶的能量系數(shù)、帶寬、載波頻率和延遲時(shí)間，n(t)為高斯白噪聲。為了仿真實(shí)際MWC采樣過(guò)程，用5倍奈奎斯特率的離散信號(hào)來(lái)表示連續(xù)信號(hào)，采用數(shù)字乘法和數(shù)字濾波運(yùn)算來(lái)仿真模擬乘法器和模擬濾波器的實(shí)際處理效果，利用抽取的方法實(shí)現(xiàn)濾波后的采樣過(guò)程。每個(gè)實(shí)驗(yàn)中，以下過(guò)程重復(fù)500次，將成功次數(shù)的百分率作為成功概率。

(1)采樣系統(tǒng)的符號(hào)波形pi(t)按均勻分布隨機(jī)產(chǎn)生；

(2)各個(gè)頻帶的載波頻率fn在區(qū)間 [-fNYQ/2,fNYQ/2]內(nèi)按均勻分布隨機(jī)產(chǎn)生；

(3)重構(gòu)時(shí)先將IMV問(wèn)題轉(zhuǎn)換成MMV問(wèn)題，然后采用 OMP(正交匹配追蹤)MMV 算法[12]或SAMP MMV算法來(lái)估計(jì)支撐集；

需要指出的是，與 SAMP MMV算法不同，OMP MMV算法需要利用稀疏度來(lái)控制迭代的次數(shù)。因此，涉及到OMP MMV算法時(shí)則使用確切稀疏度作為最大迭代次數(shù)，盡管該數(shù)值在實(shí)際中無(wú)法準(zhǔn)確預(yù)知。

5.1 定理1條件下SAMP MMV算法的有效性驗(yàn)證

本節(jié)通過(guò)對(duì)比SAMP MMV算法和OMP MMV算法重構(gòu)支撐集的成功率來(lái)說(shuō)明在定理 1條件下即fs≥fp≥ m ax(Bi)時(shí)SAMP MMV算法的有效性。不失一般性，本實(shí)驗(yàn)以 4個(gè)頻帶的信號(hào)為例，具體參數(shù)設(shè)置如下：En={1,2};Bn= { 40,40}MHz;τn={0.75,1.88} μs ，載波頻率fn隨機(jī)地分布在 [-fNYQ/2,fNYQ/2]內(nèi),fNYQ= 1 0GHz。采樣參數(shù)按照定理1設(shè)置：fs=fp=fNYQ/197= 5 0.76MHz;L0=98,L=2L0+1 =1 97;m=50。

圖4分別給出了在不同信噪比條件下，當(dāng)通道數(shù)m在區(qū)間[10, 50]內(nèi)以1為步進(jìn)變化時(shí)兩種算法的重構(gòu)成功率情況。從圖中可以看出，當(dāng)m＜20時(shí)，在無(wú)噪聲情況下，SAMP MMV算法略優(yōu)于OMP MMV算法，在有噪聲情況下，這種優(yōu)勢(shì)變得不明顯。當(dāng)m≥20時(shí)，隨著信噪比的降低，SAMP MMV算法的重構(gòu)率優(yōu)勢(shì)趨于明顯，如當(dāng)m=30，信噪比為30 dB, 15 dB, 5 dB時(shí)，該算法的重構(gòu)率分別比OMP MMV算法提高了1%, 2%, 4%。為了具體地說(shuō)明重構(gòu)效果，圖5給出了當(dāng)m=20,fn={3.53,4.71}GHz，信噪比為15 dB時(shí)該算法重構(gòu)前后信號(hào)的時(shí)域和頻域形式。如圖所示，時(shí)域和頻域的重構(gòu)結(jié)果均較好地再現(xiàn)了原信號(hào)的波形與頻譜。上述結(jié)果表明，在定理1條件下，SAMP MMV算法可以在不需要稀疏度的情況下實(shí)現(xiàn)精確重構(gòu)，且總體重構(gòu)率情況與OMP MMV算法相當(dāng)。

5.2 改進(jìn)定理?xiàng)l件下SAMP MMV算法的有效性驗(yàn)證

圖4 定理1條件下算法重構(gòu)成功率比較

圖5 SAMP MMV算法盲重構(gòu)時(shí)、頻域波形比較

5.1節(jié)驗(yàn)證了在定理1條件下SAMP MMV算法的有效性，那么在改進(jìn)定理?xiàng)l件下算法是否有效呢？本節(jié)針對(duì)fs≥fp＜ m ax(Bi)條件下的重構(gòu)情況進(jìn)行數(shù)值仿真。此時(shí)MWC原有重構(gòu)算法不再適用，因此實(shí)驗(yàn)只觀察SAMP MMV算法的重構(gòu)率情況。本實(shí)驗(yàn)在 5.1節(jié)的基礎(chǔ)上，提高頻帶寬度使Bn={80,80}MHz，信號(hào)其它參數(shù)及系統(tǒng)采樣參數(shù)保持不變。

圖6給出了信噪比為{-5,0,5,30} dB時(shí)重構(gòu)成功率隨采樣通道數(shù)目變化的情況。可以看出，成功重構(gòu)所需要的通道數(shù)隨著信噪比的提高而減少。在信噪比為30 dB時(shí)，正如改進(jìn)定理所述，只要滿足m≥ 2K≈ 2 4則重構(gòu)成功的概率接近于1。當(dāng)信噪比降低時(shí)，則需要增加MWC的通道數(shù)來(lái)提高重構(gòu)的成功率，如信噪比為5 dB, 0 dB, -5 dB時(shí)，分別至少需要28, 31, 41個(gè)通道。綜上所述，SAMP MMV算法在未知稀疏度的情況下仍然可以實(shí)現(xiàn)高概率重構(gòu)，表明該算法在改進(jìn)定理?xiàng)l件下是有效的。

5.3 寬帶頻譜感知應(yīng)用舉例

本節(jié)應(yīng)用全盲亞奈奎斯特采樣方法，對(duì)認(rèn)知無(wú)線電頻譜進(jìn)行感知，考察該方法在頻譜感知應(yīng)用中的盲重構(gòu)效果。假設(shè)在某一時(shí)刻有10個(gè)主用戶同時(shí)在帶寬Sn= { 6,8,10,20,60,10,15,6,6,80}MHz的10個(gè)信道內(nèi)進(jìn)行通信，由于某些用戶信道相鄰，一共形成了12個(gè)(6對(duì))頻帶，各信道的中心頻率(載波)fn={1,1.007,1.016,2.5,2.54,2.575,3,3.1,3.3,4.58} GHz，信號(hào)的能量系數(shù)En= { 4,4,4,2,2,2,4,3,4,4}，環(huán)境噪聲為加性高斯白噪聲，將整個(gè)頻帶的信噪比設(shè)置在10 dB，圖7(a)所示為其奈奎斯特采樣率下頻譜形式。

圖6 改進(jìn)定理?xiàng)l件下不同信噪比時(shí)盲重構(gòu)成功率情況

采用的MWC參數(shù)配置如下：fs=fp=fNYQ/197=50.76MHz;L0=98,L=2L0+1 =1 97;m=50。

圖7(b)給出了盲重構(gòu)的頻譜結(jié)果。從圖中可以看出，盡管原信號(hào)中的最大頻帶寬度(80 MHz)超過(guò)了fp，在改進(jìn)定理?xiàng)l件下，SAMP MMV算法在未知稀疏度的前提下仍然實(shí)現(xiàn)了較好的重構(gòu)。由于重構(gòu)算法只對(duì)找出的支撐集所對(duì)應(yīng)的分量進(jìn)行了恢復(fù)，而其它分量強(qiáng)制為零，因此在恢復(fù)的頻譜中消除了大量噪聲。此外，由于該算法具有稀疏度自適應(yīng)的特點(diǎn)能夠較好地克服過(guò)匹配問(wèn)題，因而在恢復(fù)的頻譜中只存在較少的雜波，達(dá)到了較好的恢復(fù)效果。

圖7 寬帶頻譜感知應(yīng)用實(shí)例

需要說(shuō)明的是，為了直觀地說(shuō)明本文所提出的采樣方法的重構(gòu)性能，本節(jié)才恢復(fù)了完整的頻譜。當(dāng)該采樣方法應(yīng)用于寬帶頻譜感知時(shí)，只需重構(gòu)出支撐集及其對(duì)應(yīng)的低速采樣值，然后根據(jù)其能量就可以進(jìn)一步判斷出各個(gè)頻帶是否真正存在主用戶，而無(wú)需恢復(fù)出完整的頻譜。限于篇幅，將在其他文章中繼續(xù)研究。

6 結(jié)論

本文通過(guò)改進(jìn) MWC重構(gòu)理論并引入 SAMP MMV重構(gòu)算法，提出了一種可用于寬帶頻譜感知的全盲亞奈奎斯特采樣方法。所提出的MWC重構(gòu)充分條件消除了信號(hào)最大頻帶寬度對(duì)采樣的約束。所提出的SAMP MMV算法在重構(gòu)中可以不需要確切的稀疏度信息。實(shí)驗(yàn)結(jié)果表明，在未知稀疏度的前提下，該算法的性能與已知確切稀疏度時(shí)的OMP MMV算法相當(dāng)，可以較好地實(shí)現(xiàn)盲譜重構(gòu)。本文所提出的全盲亞奈奎斯特采樣方法為寬帶內(nèi)多信道快速頻譜感知提供了一種有效途徑。

[1]Taherpour A, Gazor S, and Nasiri-Kenari M. Wideband spectrum sensing in unknown white Gaussian noise[J].IET Communications, 2008, 2(6): 763-771.

[2]趙春暉, 馬爽, 楊偉超. 基于分形盒維數(shù)的頻譜感知技術(shù)研究[J]. 電子與信息學(xué)報(bào), 2011, 33(2): 475-478.Zhao C H, Ma S, and Yang W C. Spectrum sensing in cognitive radios based on fractal box dimension[J].Journal of Electronics&Information Technology, 2011, 33(2): 475-478.

[3]Hwang C H, Lai G L, and Chen S C. Spectrum sensing in wideband OFDM cognitive radios[J].IEEE Transactions on Signal Processing, 2010, 58(2): 709-719.

[4]Donoho D L. Compressed sensing[J].IEEE Transactions on Information Theory, 2006, 52(4): 1289-1306.

[5]Baraniuk R. Compressive sensing[J].IEEE Signal ProcessingMagazine, 2007, 24(4): 118-121.

[6]Polo Y L, Wang Y, Pandharipande A,et al.. Compressive wide-band spectrum sensing[C]. IEEE Int. Conf. on Acoustics,Speech, and Signal Processing (ICASSP), Taipei, China,April 2009: 2337-2340.

[7]Yu Z Z, Chen X, Hoyos S,et al.. Mixed-signal parallel compressive spectrum sensing for cognitive radios[J].International Journal of Digital Multimedia Broadcasting,2010, 2010: 1-10.

[8]Tian Z and Giannakis G B. Compressed sensing for wideband cognitive radios[C]. IEEE Int. Conf. on Acoustics, Speech,and Signal Processing (ICASSP), Honolulu, HI, April 2007, 4:IV-1357-IV-1360.

[9]Mishali M and Eldar Y C. From theory to practice:sub-Nyquist sampling of sparse wideband analog signals[J].IEEE Journal of Selected Topics in Signal Processing, 2010,4(2): 375-391.

[10]Mishali M and Eldar Y C. Reduce and boost: recovering arbitrary sets of jointly sparse vectors[J].IEEE Transactions on Signal Processing, 2008, 56(10): 4692-4702.

[11]Cotter F S, Rao D B, Engan K,et al.. Sparse solutions to linear inverse problems with multiple measurement vectors[J].IEEE Transactions on Signal Processing, 2005, 53(7):2477-2488.

[12]Chen J and Huo X M. Theoretical results on sparse representations of multiple-measurement vectors[J].IEEE Transactions on Signal Processing, 2006, 54(12): 4634-4643.

[13]Do T T, Gan L, Nguyen N,et al.. Sparsity adaptive matching pursuit algorithm for practical compressed sensing[C].Asilomar Conference on Signals, Systems, and Computers,Pacific Grove, California, October 2008: 581-587.

[14]Dai W and Milenkovic O. Subspace pursuit for compressive sensing signal reconstruction[J].IEEE Transactions on Information Theory, 2009, 55(5): 2230-2249.