頻譜檢測(cè)門(mén)限與n-out-of-K融合規(guī)則的聯(lián)合優(yōu)化?

張曉格，張士兵，邱恭安（南通大學(xué)電子信息學(xué)院，江蘇南通226019）

頻譜檢測(cè)門(mén)限與n-out-of-K融合規(guī)則的聯(lián)合優(yōu)化?

張曉格，張士兵，邱恭安
（南通大學(xué)電子信息學(xué)院，江蘇南通226019）

為了減少計(jì)算復(fù)雜度，提出聯(lián)合優(yōu)化檢測(cè)門(mén)限λ和n-out-of-K融合規(guī)則的算法。以λ和n為參數(shù)建立目標(biāo)函數(shù)，并將參數(shù)以二進(jìn)制形式表示，從而把算法轉(zhuǎn)化為組合優(yōu)化問(wèn)題。接著，采用基于樣值修改的互熵優(yōu)化方法漸次逼近最優(yōu)的參數(shù)。仿真表明，該算法在獲得與已有算法幾乎相當(dāng)?shù)目傚e(cuò)誤率情況下，可有效降低平均搜索次數(shù)，且隨著K的增加搜索次數(shù)增加更平緩。

認(rèn)知無(wú)線電；協(xié)作感知；頻譜檢測(cè)；融合規(guī)則；虛警率；漏檢率；互熵

1 引言

當(dāng)前，認(rèn)知無(wú)線電［1］已經(jīng)成為緩解主用戶的低頻譜占用率與無(wú)線頻譜資源短缺之間矛盾的重要技術(shù)。次用戶通過(guò)各種方法尋找頻譜空洞并接入，以充分利用閑置頻譜資源。尋找頻譜空洞的前提是高效的頻譜感知［2］。有關(guān)頻譜感知的研究可以分為兩類：本地感知和協(xié)作感知。本地感知方面以能量檢測(cè)為代表，因其實(shí)現(xiàn)簡(jiǎn)單且無(wú)需主用戶的先驗(yàn)知識(shí)而得到廣泛使用。但能量檢測(cè)也存在明顯不足，即單個(gè)節(jié)點(diǎn)的感知能力有限（尤其是在低信噪比和深衰落環(huán)境中），認(rèn)知用戶可能無(wú)法獲得足夠的主用戶信息。常用的解決辦法是協(xié)作感知［3］，利用認(rèn)知用戶間的相互協(xié)作，使得單個(gè)認(rèn)知用戶可以從其他認(rèn)知用戶獲得協(xié)助，從而提高感知能力。

協(xié)作感知中有一類重要的方式即基于決策融合［4］，其分為本地感知、向融合中心傳遞感知結(jié)果以及融合中心根據(jù)某種規(guī)則作出全局判決3個(gè)步驟。在認(rèn)知無(wú)線電中，虛警概率和漏檢概率是兩個(gè)重要的指標(biāo)，兩者越低越好，但這兩個(gè)指標(biāo)通常相互排斥，為了盡可能同時(shí)降低兩者，文獻(xiàn)［5］對(duì)決策融合下的感知系統(tǒng)進(jìn)行了優(yōu)化，提出了總錯(cuò)誤率（漏檢概率+虛警概率之和）最小的優(yōu)化目標(biāo)，并分別給出了n-out-of-K融合規(guī)則下最小化總錯(cuò)誤率的最優(yōu)n值，以及最佳檢測(cè)門(mén)限值，有效提高了感知效果。針對(duì)文獻(xiàn)［5］分別優(yōu)化檢測(cè)門(mén)限和n值的不足，文獻(xiàn)［6］在其基礎(chǔ)上，在考慮感知用戶到融合中心存在傳輸錯(cuò)誤情形同時(shí)，進(jìn)一步提出聯(lián)合優(yōu)化檢測(cè)門(mén)限與n值的算法，從而進(jìn)一步降低總錯(cuò)誤率。

但是，聯(lián)合優(yōu)化算法不易采用類似文獻(xiàn)［5］中的解析方法，因而在文獻(xiàn)［6］所提改進(jìn)算法中n值的優(yōu)化采用全搜索方法，且針對(duì)每一個(gè)n值采用二分法搜索方法優(yōu)化檢測(cè)門(mén)限。雖然二分法搜索方法大大降低了檢測(cè)門(mén)限的優(yōu)化復(fù)雜度，但由于融合規(guī)則中n值的優(yōu)化還是采用了全搜索方法，因此當(dāng)K值較大時(shí)，搜索次數(shù)仍然較高。為此，需要研究更好的搜索方法。將目標(biāo)函數(shù)中的參數(shù)λ與n以二進(jìn)制形式表示，則參數(shù)λ與n的取值分別由一連串二進(jìn)制數(shù)0或1表示。不同的0、1組合，得到不同的參數(shù)λ與n，從而把參數(shù)的優(yōu)化轉(zhuǎn)化為對(duì)應(yīng)二進(jìn)制數(shù)串中0或1的組合問(wèn)題。采用互熵優(yōu)化方法求解該組合問(wèn)題，可得到近似最優(yōu)的優(yōu)化參數(shù)。由于不同0、1組合下參數(shù)的二進(jìn)制表達(dá)與原參數(shù)并不完全對(duì)應(yīng)，在計(jì)算過(guò)程中還對(duì)樣值同步修改。通過(guò)上述處理，參數(shù)的優(yōu)化避免了全搜索方法，故可降低計(jì)算復(fù)雜度。

2 系統(tǒng)模型

如圖1所示，考慮含有一個(gè)主用戶、一個(gè)融合中心和K個(gè)次用戶的認(rèn)知無(wú)線網(wǎng)絡(luò)。采用決策融合，每個(gè)次用戶分別獨(dú)立地進(jìn)行基于能量檢測(cè)的本地感知，然后把本地判決結(jié)果通過(guò)存在傳輸差錯(cuò)的正交控制信道發(fā)送到融合中心，融合中心則根據(jù)n-out

of-K融合規(guī)則得出全局的判決結(jié)果。假定主用戶到這K個(gè)用戶的瞬時(shí)信噪比服從獨(dú)立指數(shù)同分布，均值為ρ，這K個(gè)用戶到融合中心的瞬時(shí)信噪比也服從獨(dú)立指數(shù)同分布，均值為γ。同時(shí)，假定各個(gè)次用戶采用相同的判決門(mén)限λ。

在瑞利衰落信道下，單個(gè)用戶平均虛警概率和檢測(cè)概率分別為［7］

假定主用戶存在的概率為Pr（H1），不存在的概率為Pr（H0）=1-Pr（H1），則總的錯(cuò)誤率為

顯然，總錯(cuò)誤率與檢測(cè)門(mén)限λ以及n的值有關(guān)，在其他參數(shù)固定的情況下，調(diào)節(jié)這兩個(gè)參數(shù)，將得到不同的頻譜感知能力，因此存在如下的聯(lián)合優(yōu)化問(wèn)題：

3 聯(lián)合優(yōu)化算法

3.1 問(wèn)題轉(zhuǎn)化

3.2 算法分析

顯然，式（7）是一個(gè)組合優(yōu)化問(wèn)題，即通過(guò)選擇合適的組合系數(shù)ωi使得總錯(cuò)誤率最低。最優(yōu)的方法是采用全搜索方法，窮盡所有可能組合，但這樣一種方法需要大量的計(jì)算，一般不可取。注意到最優(yōu)的組合系數(shù)集ω*是稀少事件，因此可以用最近受到越來(lái)越多關(guān)注的互熵優(yōu)化算法［8］獲得優(yōu)化的ω設(shè)置。作為一種近年來(lái)興起的具有全局優(yōu)化特點(diǎn)的數(shù)學(xué)工具，已經(jīng)在通信領(lǐng)域得到越來(lái)越多的關(guān)注［9-12］。

3.2.1 互熵方法

互熵方法是一個(gè)不斷逼近最優(yōu)的過(guò)程，每個(gè)巡回過(guò)程主要包含兩個(gè)步驟，考慮第t個(gè)過(guò)程。

（1）ω［t］中的每一項(xiàng)可以認(rèn)為是Bernoulli隨機(jī)變量，即Pr（1（ωi［t］）=1）=pi［t］，Pr（1（ωi［t］）=0）=1-pi［t］。指示函數(shù)1（ωi［t］）∈{0，1}，指示ω［t］中的第i項(xiàng)是否為1，若為1則1（ωi［t］）=1，否則1（ωi［t］）=0。ω［t］的概率密度函數(shù)為

根據(jù)該概率密度函數(shù)產(chǎn)生J個(gè)樣值ω1［t］，…，

（2）計(jì)算Qeω［t］），…，Qe［t］），并按照從小到大排列為Qe［t］）≤…≤Qe（］），并認(rèn)為位于前「θ（J）位的樣值，即小于或等于優(yōu)選門(mén)限

的樣值為好樣值（低總錯(cuò)誤率），參數(shù)θ用來(lái)界定前多少位為好樣值。根據(jù)確定的好樣值計(jì)算pi［t］以在下一巡回過(guò)程產(chǎn)生更佳的樣值與優(yōu)選門(mén)限：

隨著巡回過(guò)程的進(jìn)行，最后pi［t］的值必然為0或者為1。pi［t］等于0對(duì)應(yīng)ωi=0，pi［t］等于1對(duì)應(yīng)ωi=1。在如何終止巡回過(guò)程方面，天然的方法是當(dāng)pi［t］∈0，{} 1，i=u+v，u+v-1，…，1時(shí)，算法終止。不過(guò)在后面的仿真過(guò)程中，發(fā)現(xiàn)在巡回過(guò)程中當(dāng)pi［t］大于一個(gè)接近1的值后，其最終值為1的概率很大，而當(dāng)pi［t］小于一個(gè)接近0的值后，則其最終值為0的概率很大。為減少搜索次數(shù)而幾乎不降低性能，算法采用如下的終止條件，即當(dāng)所有pi［t］≥ε或pi［t］≤1-ε時(shí)算法終止，其中ε是非常接近或者等于1的數(shù)。

3.2.2 樣值修改

如前所述，在根據(jù)概率密度函數(shù)f（ω［t］，p［t -1］）產(chǎn)生的J個(gè)樣值ω1［t］，…，ωJ［t］中，其中的部分樣值可能并不符合要求，不妨設(shè)該樣值為ωj［t］所謂不符合要求，是指其使得M，在算法中i=u+v應(yīng)該通過(guò)修改不符合要求的樣值排除這種情況的出現(xiàn)。修改的目的是降低的值使得該樣值符合要求。因?yàn)楦怕蕄i［t］低意味著最終概率為0的概率高，所以修改算法就是分別逐次把與｛pu［t-1］，…，p1［t -1］｝和｛pu+v［t-1］，…，pu+1［t-1］｝中的概率最低項(xiàng)對(duì)應(yīng)的組合系數(shù)置0，直到使得樣值符合要求。

0.5 ；

（2）循環(huán)

①按式（8）產(chǎn)生J個(gè)樣值ω1［t］，…，ωJ［t］；

②如果u≠lb（K+1），則修改樣值；

③如果v≠lb（M+1），則修改樣值；

④計(jì)算各樣值的總錯(cuò)誤率，從小到大排列，根據(jù)式（9）得到優(yōu)選門(mén)限r(nóng)［t］；

⑤按式（10）計(jì)算p［t］；

⑥按式（11）平滑p［t］；

⑦如果不滿足終止條件，則t=t+1，返回①；

（3）結(jié)束

①對(duì)所有i，若pi［t］≥ε則ωi［t］=1，否則ωi［t］=0；

②輸出優(yōu)化后的組合系數(shù)ωopt=｛ωu+v［t］，ωu+v-1［t］，…，ω1［t］｝，進(jìn)而得到優(yōu)化的檢測(cè)門(mén)限λopt與nopt值。

3.2.3 算法步驟

綜合上面的分析，總結(jié)出基于互熵方法的聯(lián)合優(yōu)化算法如下。

（1）初始化

t=1，p［0］=0.5，0.5，…，[ 0.5，0.5，…，0.5]；

3.3 復(fù)雜度

精確地確定算法的復(fù)雜度比較困難，故以算法的搜索次數(shù)近似衡量算法的復(fù)雜性，則所提算法的復(fù)雜度為t次搜索，復(fù)雜度取決于收斂速度的快慢。已有算法的復(fù)雜度也以算法的搜索次數(shù)衡量，設(shè)尋找最優(yōu)門(mén)限的二分法搜索方法的搜索次數(shù)為?t，則總的搜索次數(shù)為K?t，復(fù)雜度取決于最優(yōu)門(mén)限的搜索速度及K值的大小。盡管兩種算法每次搜索的實(shí)際計(jì)算量并不相同，但通過(guò)搜索次數(shù)的仿真比較，可以大致給出兩種算法在復(fù)雜度上的差別。

4 仿真結(jié)果

在以下所有仿真中，無(wú)論是所提算法還是已有算法，均設(shè)μ=4，Pr（H1）=0.4，ρ=γ，檢測(cè)門(mén)限的范圍為0～35，檢測(cè)門(mén)限的精度為0.35。在上述參數(shù)確定的情況下，已有算法結(jié)果具有唯一性，但所提算法由于是一種隨機(jī)優(yōu)化算法，其優(yōu)化結(jié)果存在一定的小范圍波動(dòng)，在此取其平均值，其中對(duì)搜索次數(shù)平均結(jié)果還要進(jìn)一步四舍五入處理。另外，假定在所提算法中，如果沒(méi)有額外說(shuō)明，則仿真中J=15，δ =0.7，θ=0.1，ε=1。

圖2給出了不同信噪比ρ下兩種算法獲得的總錯(cuò)誤率比較，由圖可見(jiàn)，兩者獲得的總錯(cuò)誤率幾乎一樣，所提算法的總錯(cuò)誤率略高。

圖3 給出了K值不同情形下兩種算法需要的搜索次數(shù)比較。結(jié)果顯示原算法隨著K的增加而線性增加。當(dāng)K較小時(shí)，原算法需要的搜索次數(shù)并不是很高，如K=10時(shí)搜索次數(shù)為70；K=20時(shí)搜索次數(shù)為140；但當(dāng)K=40時(shí)搜索次數(shù)為280；當(dāng)K =50時(shí)搜索次數(shù)為350。反觀所提算法，其需要的搜索次數(shù)明顯低于原算法，如K=10時(shí)搜索次數(shù)為7，K=30時(shí)搜索次數(shù)為10。更為突出的是隨著K的增加，所提算法需要的搜索次數(shù)增加極其緩慢，K從10增加到50，搜索次數(shù)僅從7增加到12。因此，K越大，所提算法的搜索次數(shù)優(yōu)勢(shì)越大。這種現(xiàn)象出現(xiàn)的原因是由于K從10到50時(shí)，根據(jù)u=lb（K +1）」，u僅從4增加到6，導(dǎo)致組合優(yōu)化問(wèn)題的規(guī)模并未明顯增加。

圖4 和圖5分別給出了ε變化情況下所提算法獲得的總錯(cuò)誤率及所需搜索次數(shù)比較。ε=1時(shí)的算法總錯(cuò)誤率最低，但搜索次數(shù)也最高；ε=0.7時(shí)的算法總錯(cuò)誤率最高低，但搜索次數(shù)也最低；ε等于0.9、0.8時(shí)的算法總錯(cuò)誤率增加很不明顯，幾乎可以忽略，搜索次數(shù)則有一定的下降。正如在確定終止條件時(shí)指出的那樣，出現(xiàn)這種現(xiàn)象的原因是因?yàn)樵谘不剡^(guò)程中當(dāng)組合系數(shù)ωi的值為1或0的趨勢(shì)確定以后，最終變化的可能性很低。因此，通過(guò)設(shè)置合適的ε，可保證低總錯(cuò)誤率的前提下明顯降低搜索次數(shù)。

圖6 和圖7分別給出了不同數(shù)量樣值下總錯(cuò)誤率、搜索次數(shù)的比較。

根據(jù)仿真結(jié)果，發(fā)現(xiàn)J的大小與總錯(cuò)誤率和搜索次數(shù)均有關(guān)。J越大，獲得的總錯(cuò)誤率越小，所需搜索次數(shù)也越少。這是因?yàn)楫?dāng)J增加時(shí)，每次搜索的準(zhǔn)確性增加了，所以，不僅總錯(cuò)誤率與優(yōu)化結(jié)果的波動(dòng)性降低，而且需要的搜索次數(shù)也因?yàn)樗阉鳒?zhǔn)確度的提高而降低了。之所以準(zhǔn)確性提高，一是因?yàn)闃又翟黾?，有利于下一巡回過(guò)程形成更精確的概率密度函數(shù)，從而也有利于下一巡回過(guò)程形成更好的樣值；二是由于可能要修改樣值，使得最終的樣值并不完全符合上一巡回過(guò)程確定的概率密度函數(shù)，增加樣值使得需要修改的樣值所占全部樣值比例降低，在一定程度上減少了這種不利影響。

5 結(jié)束語(yǔ)

針對(duì)已有的聯(lián)合優(yōu)化檢測(cè)門(mén)限λ和n-out-of-K融合規(guī)則算法復(fù)雜度仍然較高的不足，提出了基于互熵方法的漸次逼近最優(yōu)參數(shù)的算法，避免了采用全搜索方法優(yōu)化n值。與已有算法相比，總錯(cuò)誤率相當(dāng)，但搜索次數(shù)明顯降低，且隨著K的增加，搜索次數(shù)增加較已有算法更平緩，因此該算法在協(xié)作認(rèn)知用戶數(shù)量較多時(shí)更能體現(xiàn)出其低復(fù)雜度特性。同時(shí)，進(jìn)一步研究在不同樣值數(shù)量J以及不同ε下所提算法的總錯(cuò)誤率與搜索次數(shù)，結(jié)果發(fā)現(xiàn)合適的ε或J值可以更好地權(quán)衡總錯(cuò)誤率與搜索次數(shù)，因此該算法還可以在實(shí)際使用中根據(jù)性能與復(fù)雜度的不同要求而設(shè)置不同的ε或（和）J值，具有比已有算法更佳的適應(yīng)能力。

［1］Mitola J，Maguire GQ.Cognitive radio：making software radios more personal［J］.IEEE Personal Communications，1999，6（4）：13-18.

［2］Akyildiz IF，LeeWY，Mehmet C V，et al.Next generation／dynamic spectruMaccess／cognitive radio wireless networks：A survey［J］.Computer Networks，2006，50（13）：2127-2159.

［3］Ganesan G，Li Y G.Cooperative spectruMsensing in cognitive radio-part I：two user networks［J］.IEEE Transactions on Wireless Communications，2007，6（6）：2204-2213.

［4］Myung K B，Jin Y K.Effective signal detection using cooperative spectruMsensing in cognitive radio systems［C］／／Proceedings of2009 International Conference on Advanced Communication Technology.Phoenix Park，Ireland：IEEE，2009：1746-1750.

［5］Wei Z，Ranjan KM，Khaled B L.Optimization of cooperative spectruMsensing with energy detection in cognitive radio networks［J］.IEEE Transactions on Wireless Communications，2009，8（12）：5761-5766.

［6］Ha N V，Tran TD，Kong H Y.An optimal cooperative spectruMsensingmethod in cognitive radio network over Rayleigh fading channel［C］／／Proceedings of 2011 International Federation for Information Processing Wireless Days.Niagara Falls，USA：IEEE，2011：1-5.

［7］DighaMF F，AlouiniMS，Simon MK.On the energy detection of unknown signals over fading channels［C］／／Proceedings of2003 International Conference on Communication.Anchorage，AK，USA：IEEE，2003：3575-3579.

［8］Rubinstein R Y，Kroese D P.The cross-entropymethod：A unified approach to combinatorial optimization，Monte-Carlo simulation，and machine learning［M］.Berlin，Germany：Springer，2004.

［9］Zhang Y Y，JiC L，WasiMQM，etal.Receive antenna selection for MIMO systems over correlated fading channels［J］.IEEE Transactions onWireless Communications，2009，8（9）：4393-4399.

［10］Zhang Y Y，Zheng G，Wong K K，etal.Near-optimal joint antenna selection for amplify-and-forward relay networks［J］.IEEETransactions onWireless Communications，2010，9（8）：2401-2407.

［11］Chen JC，Wen CK.Near-optimal relay subsetselection for two-way amplify-and-forward MIMO relaying systems［J］. IEEE Transactions on Wireless Communications，2011，10（1）：37-42.

［12］Wang Y J，ChenW，Chintha T.PAPR reductionmethod based on parametric minimuMcross entropy for OFDMsignals［J］. IEEECommunications Letters，2010，14（6）：563-565.

ZHANG Xiao-ge was born in Nantong，Jiangsu Province，in 1975.He received the Ph.D.degree in Signaland Information Processing froMNanjing University of Posts and Telecommunications in 2010.He is noWa lecturer.His research concerns signaland information processing inmodern communications.

Email：zxg＠ntu.edu.cn

張士兵（1962—），男，江蘇南通人，2006年于南京郵電大學(xué)獲信號(hào)與信息處理專業(yè)博士學(xué)位，現(xiàn)為教授，主要研究方向?yàn)閷拵?shù)字通信；

ZHANG Shi-bing was born in Nantong，Jiangsu Province，in 1962.He received the Ph.D.degree in Signaland Information Processing froMNanjing University of Posts and Telecommunications in 2006.He is noWa professor.His research direction is wideband digital communications.

Email：zhangshb＠ntu.edu.cn

邱恭安（1973—），男，湖北浠水人，2008年于南京郵電大學(xué)獲信號(hào)與信息處理專業(yè)博士學(xué)位，現(xiàn)為副教授，主要研究方向?yàn)闊o(wú)線通信網(wǎng)絡(luò)。

QIUGong-an was born in Xishui，Hubei Province，in 1973. He received the Ph.D.degree in Communications and Information Systems froMNanjing University of Posts and Telecommunications in 2008.He is noWan associate professor.His research direction is wireless communication networks.

Email：qiugongan＠ntu.edu.cn

Joint OptiMization of SpectruMDetection Threshold and the n-out-of-K Fusion Rule

ZHANGXiao-ge，ZHANGShi-bing，QIU Gong-an

（College of Electronic and Information，Nantong University，Nantong 226019，China）

To reduce the computational complexity，an algorithMoptimizing the spectruMdetection thresholdλ and the n-out-of-K fusion rule jointly is proposed.A target function is builtwith the parametersλand n.By binary encoding the parameters，this algorithMis converted to a combinatorial optimization problem.Then，the algorithMgains the near optimal results step by step by use of the modifying samples-based cross entropy method.Simulations shoWthe proposed algorithMdecreases the average search times effectively with almost the same total error rates as that of the existing one and the search times of the proposed algorithMincreasemore evenly than that of the existing one.

cognitive radio；cooperative sensing；spectruMdetection；fusion rule；missed-detection probability；false-alarMprobability；cross entropy

s：The National Natural Science Foundation of China（No.61071086）；Collegiate Natural Science Fund of Jiangsu Province（11KJB510021）；Nantong Applied Research Project（BK2011064）

TN929.53

A

10.3969／j.issn.1001-893x.2012.11.008

張曉格（1975—），男，江蘇南通人，2010年于南京郵電大學(xué)獲信號(hào)與信息處理專業(yè)博士學(xué)位，現(xiàn)為講師，主要研究方向?yàn)楝F(xiàn)代通信中的信號(hào)與信息處理；

1001-893X（2012）11-1746-06

2012-06-11；

2012-08-17

國(guó)家自然科學(xué)基金資助項(xiàng)目（61071086）；江蘇省高校自然科學(xué)基金資助項(xiàng)目（11KJB510021）；南通市應(yīng)用研究計(jì)劃項(xiàng)目（BK2011064）