• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于變分貝葉斯推斷的新型全局頻譜協(xié)作感知算法

    2016-10-14 11:38:07吳名宋鐵成胡靜沈連豐
    通信學(xué)報(bào) 2016年2期
    關(guān)鍵詞:后驗(yàn)全局頻段

    吳名,宋鐵成,胡靜,沈連豐

    ?

    基于變分貝葉斯推斷的新型全局頻譜協(xié)作感知算法

    吳名,宋鐵成,胡靜,沈連豐

    (東南大學(xué)移動(dòng)通信國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210096)

    為了實(shí)現(xiàn)多維動(dòng)態(tài)頻譜接入,首先給出了主用戶的全局功率譜近似模型,并構(gòu)建了新型全局頻譜協(xié)作感知算法的總體流程,以獲得主用戶網(wǎng)絡(luò)中占用頻段、功率及位置等全局信息。接著利用變分貝葉斯推斷技術(shù),設(shè)計(jì)了相應(yīng)的模型系數(shù)向量估計(jì)器。仿真結(jié)果表明,該方法采用的近似模型具有較好的準(zhǔn)確性,相應(yīng)的系數(shù)向量估計(jì)算法具有較高的有效性和收斂穩(wěn)定性,同時(shí)指明了信噪比和泄漏總虛假功率的關(guān)系以及兩者對(duì)均方誤差性能的影響。此外,還證明了該方法通過(guò)利用系數(shù)向量的稀疏性,而在均方誤差性能上具有較大優(yōu)勢(shì)。

    認(rèn)知無(wú)線電;全局頻譜協(xié)作感知;變分貝葉斯推斷;稀疏性

    1 引言

    目前,無(wú)線通信領(lǐng)域中存在頻譜資源日益匱乏而現(xiàn)有頻譜利用效率低下這一困境,認(rèn)知無(wú)線電技術(shù)正是為解決這一問(wèn)題而提出的。認(rèn)知無(wú)線電技術(shù)的基礎(chǔ)和關(guān)鍵是頻譜感知技術(shù),該技術(shù)主要用于判斷授權(quán)頻段是否被主用戶占用。目前,其主要采取協(xié)作的方式進(jìn)行頻譜感知,以利用不同從用戶的采樣點(diǎn)在時(shí)間、空間上的獨(dú)立性或不相關(guān)性,實(shí)現(xiàn)分集、增強(qiáng)感知性能,從而達(dá)到快速、可靠感知的目的[1~3]。

    但是因?yàn)閺挠脩艟W(wǎng)絡(luò)覆蓋范圍較大,主用戶通信范圍通常只占據(jù)其中一部分。而在其他地方,由于距離主用戶較遠(yuǎn)、主用戶信號(hào)功率較弱、建筑物遮蔽等原因,從用戶對(duì)授權(quán)頻段的使用往往既不會(huì)對(duì)主用戶通信產(chǎn)生有害的干擾,也不會(huì)受到主用戶的有害影響。同時(shí)由于主/從用戶一般存在移動(dòng)性,導(dǎo)致主用戶通信影響范圍和頻譜空洞所處位置也隨時(shí)間不斷變化。由此可知,為實(shí)現(xiàn)頻譜在多維度上的動(dòng)態(tài)復(fù)用,從用戶網(wǎng)絡(luò)如何感知整個(gè)覆蓋范圍中主用戶發(fā)射機(jī)所使用的頻段、功率及其所處位置等全局頻譜信息成為必須首先解決的問(wèn)題。文獻(xiàn)[4]詳細(xì)描述了這一問(wèn)題,并進(jìn)一步指出由于主/從用戶天線所處位置不同導(dǎo)致的路徑損耗差異會(huì)產(chǎn)生很大的動(dòng)態(tài)頻譜復(fù)用機(jī)會(huì)。而文獻(xiàn)[5]則介紹了如何利用現(xiàn)代空間統(tǒng)計(jì)學(xué)來(lái)定量描述主/從用戶頻譜在大范圍內(nèi)分布的統(tǒng)計(jì)模型,并指出設(shè)計(jì)全局頻譜分布信息感知算法是實(shí)現(xiàn)動(dòng)態(tài)全局頻譜共享的首要步驟。而傳統(tǒng)的協(xié)作感知技術(shù)難以完成這種在具有移動(dòng)性的認(rèn)知網(wǎng)絡(luò)中,有效感知其全局頻譜信息的任務(wù)。因此最近出現(xiàn)多種考慮主/從用戶移動(dòng)性的協(xié)作感知算法以實(shí)現(xiàn)這一目標(biāo)。為反映主/從用戶移動(dòng)性導(dǎo)致的頻譜接入機(jī)會(huì),文獻(xiàn)[6]利用兩態(tài)連續(xù)時(shí)間馬爾可夫鏈構(gòu)建信道可用性模型,并根據(jù)該模型聯(lián)合優(yōu)化頻譜感知間隔等參數(shù),以最大化時(shí)域/空域頻譜接入機(jī)會(huì)。同時(shí)設(shè)計(jì)相應(yīng)的信道接入策略,以利用主/從用戶移動(dòng)性導(dǎo)致的各種頻譜接入機(jī)會(huì)。文獻(xiàn)[7]利用主/從用戶間信道特性、感知數(shù)據(jù)間相關(guān)度和主從用戶運(yùn)動(dòng)行為構(gòu)建接近真實(shí)的移動(dòng)性模型,并基于該模型,設(shè)計(jì)出一種新型的感知節(jié)點(diǎn)選擇方法,以選擇出低相關(guān)性、高感知性能的感知節(jié)點(diǎn)子集。文獻(xiàn)[8]考慮主用戶移動(dòng)性導(dǎo)致的頻譜接入機(jī)會(huì),由此擴(kuò)展出新的空時(shí)頻譜感知模型,并定義一種新的性能評(píng)價(jià)指標(biāo),即受感知參數(shù)和移動(dòng)性影響的從用戶可達(dá)傳輸容量。該性能評(píng)價(jià)指標(biāo)同時(shí)指明了感知機(jī)制不完美時(shí),感知機(jī)制的精度對(duì)信道接入概率的影響,并由此得出最優(yōu)感知時(shí)長(zhǎng),以最大化從用戶可達(dá)傳輸容量。文獻(xiàn)[9]在不考慮陰影衰落影響的情況下,利用貝葉斯分層先驗(yàn)建模技術(shù)設(shè)計(jì)了新型的協(xié)作頻譜感知方法,以得到主用戶信號(hào)的功率譜分布估計(jì)值。

    本文利用新出現(xiàn)的變分貝葉斯推斷技術(shù),設(shè)計(jì)了一種新型的全局頻譜信息協(xié)作感知算法。通過(guò)最小化全局功率譜模型系數(shù)后驗(yàn)函數(shù)的K-L(Kullback-Leibler)散度,得到該后驗(yàn)的近似值,然后根據(jù)最大后驗(yàn)概率準(zhǔn)則得到各主用戶發(fā)射機(jī)占用頻段、功率和所處位置3種信息的最優(yōu)值,由此得到指定區(qū)域內(nèi)全局頻譜分布,獲得頻譜在時(shí)間、空間、頻譜上的多維信息,使從用戶網(wǎng)絡(luò)可以高效復(fù)用空閑頻譜。與上述文獻(xiàn)相比,這種機(jī)制擁有不易受較小局部最優(yōu)值影響、得到的表達(dá)式易于處理等優(yōu)點(diǎn)。

    2 系統(tǒng)建模

    假設(shè)認(rèn)知無(wú)線網(wǎng)絡(luò)中有N個(gè)主用戶發(fā)射機(jī)正在通信,其位置坐標(biāo)集合為。在時(shí)刻,網(wǎng)絡(luò)中任意位置x處的從用戶接收到的信號(hào)y()可以表示為

    其中,h(;)和u()分別表示第個(gè)傳輸路徑上的信道沖激響應(yīng)和主用戶的發(fā)射信號(hào),n() 是方差為的加性高斯白噪聲,L是主用戶和從用戶間的傳輸路徑數(shù)。

    在不考慮噪聲的情況下,假設(shè)信道h(;)和信號(hào)u()都是平穩(wěn)的,則從用戶接收到的單個(gè)主用戶信號(hào)的自相關(guān)函數(shù)可以表示為。由于信道h()在相干時(shí)間內(nèi)保持不變,所以,則自相關(guān)函數(shù)可表示為,而接收功率譜可表示為。雖然上式中的信道增益能夠通過(guò)訓(xùn)練方式獲得,但這需要主用戶的配合,并擁有區(qū)分各主用戶信號(hào)的能力。因此本文采用一種替代的方法,即路徑損耗模型,其中,0、是預(yù)先選定的常數(shù)。值得強(qiáng)調(diào)的是,路徑損耗僅影響信道統(tǒng)計(jì)量,所以本文考慮的信道沖激響應(yīng)h(;)都是頻率選擇性的。

    而當(dāng)存在N個(gè)主用戶和噪聲時(shí),假設(shè)信道h(;)和信號(hào)u()都是不相關(guān)的,從用戶的接收功率譜可以表示為

    由于感知機(jī)制不需要知道功率譜的準(zhǔn)確值,所以可利用基擴(kuò)展模型[10]近似主用戶信號(hào)的功率譜。如圖1所示,帶寬為的主用戶信號(hào)功率譜的基擴(kuò)展模型由N個(gè)不重疊的單位矩形函數(shù)γ()組成,其表達(dá)式可寫為

    (3)

    其中,θ是模型系數(shù)。將式(3)代入式(2)可得

    圖1 信號(hào)功率譜的基擴(kuò)展近似模型

    其中,矩陣由構(gòu)成,表示全1向量,而N表示從用戶數(shù)量。當(dāng)每個(gè)從用戶將自己感知到的功率譜信息和自身的位置發(fā)送給融合中心后,融合中心可利用其構(gòu)建從用戶系統(tǒng)接收信號(hào)的全局功率譜模型,具體表達(dá)式為

    (6)

    其中,向量、和矩陣是將各從用戶的、按順序堆疊而成的NN×1維向量和NN×NN維矩陣;而式中,。假設(shè)每個(gè)從用戶可估計(jì)出本地噪聲功率σ2,由此可改為發(fā)送修正后的功率譜信息′=2,所以全局功率譜模型可進(jìn)一步簡(jiǎn)化為

    圖2 基于全局參考系統(tǒng)的虛擬網(wǎng)絡(luò)模型

    綜上所述,全局頻譜協(xié)作感知算法的核心任務(wù)是估計(jì)出全局功率譜模型的系數(shù),以求得主用戶發(fā)射機(jī)所使用的頻段、功率及其所處位置,其整體流程如圖3所示。需要注意的是,由于從用戶對(duì)主用戶的情況一無(wú)所知,所以每個(gè)從用戶必須同時(shí)感知所有授權(quán)頻段。而與傳統(tǒng)數(shù)據(jù)融合方案不同的是,本文的數(shù)據(jù)融合方案中各從用戶發(fā)送到融合中心的是其感知到的修正功率譜信息,而不是其接收信號(hào)能量信息;融合中心不是對(duì)這些信息進(jìn)行簡(jiǎn)單加權(quán)求和以構(gòu)成判決統(tǒng)計(jì)量,而是直接利用這些信息,構(gòu)建主用戶信號(hào)的全局功率譜模型,然后再利用下文所提迭代算法進(jìn)行求解。而下文所提迭代算法本質(zhì)上是對(duì)各從用戶感知信息的一種聯(lián)合優(yōu)化,所以其求解出的主用戶發(fā)射機(jī)所使用的頻段、功率及其所處位置對(duì)于整個(gè)認(rèn)知網(wǎng)絡(luò)來(lái)說(shuō)是統(tǒng)計(jì)意義上最優(yōu)的。

    3 算法描述

    變分貝葉斯推斷技術(shù)是新近出現(xiàn)的一類近似方法,其主要通過(guò)觀測(cè)數(shù)據(jù),利用泛函方法近似出復(fù)雜的統(tǒng)計(jì)模型中不可觀測(cè)變量的后驗(yàn)概率。因?yàn)槠洳捎玫囊彩亲畲蠛篁?yàn)估計(jì),即用單個(gè)最有可能的參數(shù)值來(lái)代替完全貝葉斯估計(jì),所以也可以看作是EM算法的擴(kuò)展。所以其提供的是一種次優(yōu),但具有確定解的近似后驗(yàn)估計(jì)方法。目前該方法廣泛地應(yīng)用于信號(hào)重構(gòu)、神經(jīng)元定位等領(lǐng)域[12,13]。本文擬利用該技術(shù)設(shè)計(jì)一種基于變分貝葉斯推斷的模型系數(shù)求解算法,以給出模型系數(shù)的近似后驗(yàn)概率,并利用最大后驗(yàn)概率準(zhǔn)則得到模型系數(shù)的最優(yōu)值。

    為了表示方便,本文將全局功率譜模型中未知的模型系數(shù)向量和模型誤差精度參數(shù)統(tǒng)一用向量來(lái)表示。然后,將全局功率譜模型的對(duì)數(shù)邊緣概率密度函數(shù)做如下分解

    (8)

    其中,()表示聯(lián)合后驗(yàn)概率密度函數(shù)(|′)的某種近似,()表示()下的下限,而(||)表示()下的K-L散度。

    由于各未知變量和參數(shù)互相獨(dú)立,所以()可以實(shí)現(xiàn)因子化,分解為各未知變量和參數(shù)的近似后驗(yàn)概率密度函數(shù)的乘積,即

    其中,Z表示模型中各未知變量和參數(shù)。將式(9)代入(),可得

    (10)

    由式(11)可知,求解最優(yōu)()必須知道模型的聯(lián)合概率密度函數(shù)(′,),所以如何確定其具體形式成為本文必須首先解決的問(wèn)題。

    由式(7)可知,全局功率譜模型的聯(lián)合概率密度函數(shù)(′,)的具體形式可表示為

    由于模型的誤差向量是高斯型噪聲,相應(yīng)協(xié)方差矩陣為?1,所以模型的似然函數(shù)也是高斯的,具體表達(dá)式為

    (13)

    由此可知,模型系數(shù)的共軛先驗(yàn)()也是高斯的,而模型誤差精度參數(shù)的共軛先驗(yàn)()是伽馬分布,因此有

    (15)

    其中,向量是模型系數(shù)的誤差精度參數(shù)向量;0、0是2個(gè)超參數(shù)。需注意的是,因?yàn)檎J(rèn)知網(wǎng)絡(luò)中存在2種稀疏性:一種是主用戶信號(hào)占用的頻段相對(duì)于可利用頻譜很窄;另一種是在網(wǎng)絡(luò)覆蓋范圍內(nèi),存在的主用戶數(shù)量有限、密度很小。所以模型系數(shù)向量也是稀疏的,因此可將式(14)中的均值假設(shè)為零向量。

    此外,由于也是未知的,所以將向量進(jìn)一步擴(kuò)展,加入系數(shù)誤差精度參數(shù)向量,用向量′來(lái)表示。則全局功率譜模型的聯(lián)合概率密度函數(shù)(′,)可擴(kuò)展為(′,′),其具體形式可表示為

    (17)

    其中,所有的cd也都是超參數(shù)。由于和以及都是相互獨(dú)立的,所以式(9)也具體表示為(,)=()()()。由上述可知,()、()和()分別是高斯分布、伽馬分布和伽馬分布,所以其近似值() 、()和()也應(yīng)分別是高斯分布、伽馬分布和伽馬分布,因此有,,。

    根據(jù)式(11)可得到,()最優(yōu)值*()的計(jì)算式為

    (19)

    由于該式是高斯分布的對(duì)數(shù)形式。所以()中的,,其中,的均值,而向量的各元素均值。

    由上述公式可知,求解最優(yōu)()必須依賴其他近似后驗(yàn)()和(),所以在實(shí)際使用中,()的最優(yōu)值難以直接得到,需利用迭代算法進(jìn)行求解,即每次指定一個(gè)近似后驗(yàn)q(Z),固定其他近似后驗(yàn)為上次計(jì)算得到的最優(yōu)值,利用式(11)求得該q(Z)的本次最優(yōu)值,如此反復(fù)直到得到()的全局最優(yōu)值,其具體流程如圖4所示。

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

    本文利用Matlab軟件構(gòu)建了所提方案的仿真平臺(tái),并根據(jù)基于估計(jì)值的真實(shí)全局功率譜模型和由式(7)所示的高斯型近似全局功率譜模型分為2種仿真環(huán)境,以驗(yàn)證近似模型的準(zhǔn)確性和所提算法的有效性。同時(shí)本文還比較了2種環(huán)境下不同信噪比(SNR, signal noise ratio)時(shí)的算法性能,分析了其收斂性,討論了不同信噪比時(shí)估計(jì)的準(zhǔn)確性和虛假功率泄漏的情況。最后,本文選擇了基于非負(fù)最小二乘(NNLS, non-negativity least square)準(zhǔn)則的全局頻譜協(xié)作感知算法與所提算法進(jìn)行了對(duì)比,驗(yàn)證了所提算法在均方誤差(MSE, mean square error)性能上的優(yōu)勢(shì)。

    本文仿真條件如下:假設(shè)認(rèn)知網(wǎng)絡(luò)中存在N=8個(gè)連續(xù)可用頻段,每個(gè)可用頻段寬度為1 MHz,每個(gè)主用戶通信占據(jù)1個(gè)頻段,發(fā)射功率為1 W。而從用戶在感知時(shí)掃描的頻點(diǎn)數(shù)為=64,即從用戶在每個(gè)可用頻段上采樣8個(gè)點(diǎn)。本仿真中認(rèn)知網(wǎng)絡(luò)的參考系統(tǒng)是一個(gè)由5×5個(gè)點(diǎn)構(gòu)成的正方形,也就是說(shuō)認(rèn)知網(wǎng)絡(luò)中有N=25個(gè)候選位置。而其中假設(shè)存在2個(gè)主用戶,其位置坐標(biāo)分別為(2,2),(3,5),所占用頻段分別為6號(hào)頻段和7號(hào)頻段。同時(shí)假設(shè)存在N=4個(gè)從用戶,其位置坐標(biāo)分別為(1,1)、(2,5)、(4,2)、(5,4)。仿真所用的路徑損耗相關(guān)參數(shù)為0=10 m,=4,而參考系統(tǒng)單位坐標(biāo)代表的實(shí)際距離為20 m。由該仿真條件可知,模型系數(shù)向量中有N×N=200個(gè)元素,但是非零的只有2個(gè),因此系數(shù)向量具有較強(qiáng)的稀疏性。而這2個(gè)非零元素在向量中的具體位置序號(hào)可由主用戶的位置坐標(biāo)和所占用頻段換算出來(lái),其換算公式為元素序號(hào)=[(橫坐標(biāo)?1)×5+(縱坐標(biāo)?1)]×8+占用頻段號(hào)。也就是說(shuō),2個(gè)非零元素分別在54號(hào)和119號(hào)位置上,大小為1,所以本文仿真通過(guò)考察系數(shù)向量的估計(jì)精確度來(lái)驗(yàn)證算法對(duì)主用戶發(fā)射機(jī)所使用頻段、功率及其位置的估計(jì)能力。

    本文先仿真并對(duì)比了真實(shí)模型和高斯模型2種不同仿真環(huán)境下所提算法的性能差異,以驗(yàn)證近似模型的準(zhǔn)確性及所提算法的有效性。如圖5所示,當(dāng)系統(tǒng)信噪比=?6 dB時(shí),真實(shí)模型下算法的迭代次數(shù)為119次,系數(shù)向量和其估計(jì)值的歐式距離為0.047 1。而在高斯模型下,算法的迭代次數(shù)為93次,系數(shù)向量和其估計(jì)值的歐式距離為0.039,兩者差異不大。由此可知,在2種模型下,算法的估計(jì)準(zhǔn)確度比較接近。這是因?yàn)殡S著增大,估計(jì)誤差越來(lái)越小,也就是誤差向量的協(xié)方差矩陣中的元素越來(lái)越小。這導(dǎo)致2種模型下誤差向量的值都主要分布在其均值附近,散布很小,所以此時(shí)的2種模型總體上比較接近,因此2種模型下所提算法的性能也比較接近。這一點(diǎn)可由圖6進(jìn)一步證明,當(dāng)從?10 dB變化到?6 dB時(shí),隨著增大,2種模型的吻合程度越來(lái)越高,算法性能也越來(lái)越接近。尤其是當(dāng)接收=?6 dB時(shí),兩者的均方誤差分別是?26.547 8 dB和?28.184 2 dB,僅相差1.6 dB。但是當(dāng)繼續(xù)增加時(shí),2種模型下MSE曲線間的差異卻仍然保持在1.5 dB左右。這是因?yàn)殡m然算法泄漏出的虛假功率會(huì)隨著的下降而減少,但算法在真實(shí)環(huán)境下總會(huì)泄漏出更多的虛假功率。這點(diǎn)可以由圖7證明,即當(dāng)接收=?6 dB時(shí),算法泄漏出的歸一化總虛假功率分別為0.028 4和0.011 9,兩者相差3.7 dB。而當(dāng)接收=0 dB時(shí),雖然算法泄漏出的歸一化總虛假功率分別下降為0.008 8和0.004 1,但兩者仍相差3.3 dB。

    本文還仿真并對(duì)比了2種仿真環(huán)境下不同信噪比時(shí)的所提算法性能,并對(duì)其收斂性進(jìn)行了相應(yīng)分析,同時(shí)討論了不同信噪比時(shí)估計(jì)的準(zhǔn)確性和虛假功率泄漏的情況。由圖8可知,當(dāng)系統(tǒng)信噪比分別等于0 dB、?6 dB和?10 dB時(shí),真實(shí)模型下算法的迭代次數(shù)分別為91次、119次和129次,系數(shù)向量和其估計(jì)值的歐式距離分別為0.019 4、0.047 1、0.165 2。而由圖9可知,對(duì)應(yīng)的高斯近似模型下算法的迭代次數(shù)分別為89次、93次和105次,系數(shù)向量和其估計(jì)值的歐式距離分別為0.016 1、0.039 0、0.099 7。將上述兩圖中各迭代次數(shù)進(jìn)行對(duì)比可以發(fā)現(xiàn),當(dāng)由?10 dB變化到0 dB時(shí),2種模型下算法的迭代次數(shù)差異都不是很大,也就是說(shuō)算法的收斂性在2種模型的各種下的穩(wěn)定性都比較好。但當(dāng)小于?6 dB時(shí),其估計(jì)準(zhǔn)確性都開(kāi)始明顯變差。這是因?yàn)楫?dāng)變小時(shí),接收的主用戶信號(hào)功率不變,因此主要是本地噪聲功率變大,導(dǎo)致算法的估計(jì)準(zhǔn)確性明顯下降。這點(diǎn)從系數(shù)向量的119號(hào)元素估計(jì)值的準(zhǔn)確度也可以看出。例如,真實(shí)模型下119號(hào)元素估計(jì)值在=?10 dB時(shí)是0.893 6,在=0 dB時(shí)是0.994 7。而高斯近似模型下的估計(jì)值在= ?10 dB時(shí)是0.927 0,在=0 dB時(shí)是0.991 9。由此可以明顯地看出估計(jì)準(zhǔn)確性隨變小而快速下降。同時(shí)將圖8和圖9進(jìn)一步對(duì)比可知,真實(shí)模型下的算法收斂速度在不同信噪比時(shí)都要略慢于高斯近似模型下的算法收斂速度,而對(duì)應(yīng)的估計(jì)準(zhǔn)確性也是高斯近似模型下的準(zhǔn)確度都略高于真實(shí)模型下的準(zhǔn)確度。但兩者相差不大,這也再次證明了本文基于高斯近似模型推導(dǎo)出的算法總體上可以較好的適用于真實(shí)模型。此外,由圖10可知,當(dāng)系統(tǒng)信噪比分別等于0 dB、?6 dB和?10 dB時(shí),真實(shí)模型下算法泄漏出的總虛假功率分別為0.008 8、0.028 4、0.086 7。而由圖11可知,對(duì)應(yīng)的高斯近似模型下算法泄漏出的總虛假功率分別為0.004 1、0.011 9、0.033 4。由上述兩圖可知,的減小不僅影響估計(jì)準(zhǔn)確度,還導(dǎo)致算法泄漏出的總虛假功率也隨之增加。而且由圖10可以看到,當(dāng)=?10 dB時(shí),算法收斂時(shí)泄漏出的總虛假功率已經(jīng)不是其最小值,而是出現(xiàn)了偏離的現(xiàn)象。這是因?yàn)檩^大噪聲功率影響了估計(jì)的準(zhǔn)確性,其收斂值偏離真實(shí)值較遠(yuǎn),導(dǎo)致泄漏出了更多的總虛假功率。同時(shí)將圖10和圖11進(jìn)一步對(duì)比可知,當(dāng)接收=?10 dB時(shí),所提算法泄漏出的總虛假功率相差4.1 dB,這再次證明了MSE曲線間存在的差異受到在真實(shí)環(huán)境下算法泄漏出更多的虛假功率的影響,而這種影響基本不隨而變化。

    最后,由于傳統(tǒng)的協(xié)作頻譜感知算法僅檢測(cè)了主用戶在授權(quán)頻段中是否存在,而沒(méi)有指明主用戶所處位置及其發(fā)射功率等全局頻譜信息,所以難以直接與本文所提算法進(jìn)行性能對(duì)比,為此本文選擇了基于NNLS準(zhǔn)則的全局頻譜協(xié)作感知算法與所提算法進(jìn)行對(duì)比。該算法采用的準(zhǔn)則為

    根據(jù)利用該準(zhǔn)則,本文仿真并對(duì)比了在真實(shí)模型下NNLS算法和所提算法的MSE性能。如圖12所示,當(dāng)從?10 dB增大到0 dB時(shí),所提算法和NNLS算法的均方誤差分別從?15.639 4 dB和9.456 6 dB逐漸減小到?34.232 8 dB和?18.401 0 dB,也就是說(shuō)兩者M(jìn)SE性能間隔從25.096 dB縮小到15.831 9 dB。由此可見(jiàn),雖然隨著的增加,所提算法的MSE性能領(lǐng)先NNLS算法的程度有所減小,但是仍保持著較大的優(yōu)勢(shì)。這是因?yàn)樗崴惴ɡ昧讼禂?shù)向量存在稀疏性這一隱含條件,所以其估計(jì)向量中除了54號(hào)和119號(hào)元素比較接近真實(shí)值1外,其他元素估計(jì)值都非常接近真實(shí)值0。而NNLS算法沒(méi)有利用這一稀疏性條件,所以其估計(jì)向量中除54號(hào)和119號(hào)元素外,還存在數(shù)量較多、數(shù)值較大的非零元素,這嚴(yán)重影響了NNLS的估計(jì)精度。因此通過(guò)該仿真對(duì)比,本文證明了利用系數(shù)向量稀疏性的所提算法在MSE性能上具有較大的優(yōu)勢(shì),其估計(jì)精確度較高,可以通過(guò)忽略功率譜模型系數(shù)中較小的元素,直接利用較大非零元素來(lái)得到主用戶發(fā)射機(jī)所使用的頻段、功率及其所處位置的真實(shí)值。

    圖12 本文所提算法和NNLS算法性能對(duì)比

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

    首先,本文利用接收到的主用戶信號(hào)構(gòu)建了時(shí)間域和頻率域上的全局功率譜近似模型,并在此基礎(chǔ)上設(shè)計(jì)了一種新型的全局頻譜協(xié)作感知算法的整體流程。接著,利用認(rèn)知無(wú)線網(wǎng)絡(luò)中存在的稀疏性,設(shè)計(jì)了基于變分貝葉斯推斷技術(shù)的估計(jì)器,以求得該模型系數(shù)向量的近似后驗(yàn)概率密度函數(shù),從而利用最大后驗(yàn)概率準(zhǔn)則得到系數(shù)向量的最優(yōu)值,由此揭示出主用戶發(fā)射機(jī)所使用的頻段、功率及其所處位置。最后,本文建立了2種分別基于真實(shí)模型和高斯近似模型的仿真環(huán)境,以驗(yàn)證近似模型的準(zhǔn)確性及所提算法的有效性。同時(shí)還比較了2種環(huán)境下不同信噪比時(shí)的算法性能,分析了其收斂性,討論了不同信噪比時(shí)估計(jì)的準(zhǔn)確性和虛假功率泄漏的情況。另外,在真實(shí)模型下,本文通過(guò)將所提算法和NNLS算法的MSE性能進(jìn)行仿真對(duì)比,證明了利用系數(shù)向量稀疏性的所提算法在MSE性能上具有較大的優(yōu)勢(shì)。

    [1] GOLDSMITH A, JAFAR S, MARIC I, et al. Breaking spectrum gridlock with cognitive radios an information theoretic perspective[J]. Proceedings of the IEEE, 2009, 97(5): 894-914.

    [2] LU L, ZHOU X W, ONUNKWO U, et al. Ten years of research in spectrum sensing and sharing in cognitive radio[J]. Eurasip Journal on Wireless Communications, 2012, 28: 1-16.

    [3] ZENG Y H, LIANG Y C, HOANG A T, et al. A review on spectrum sensing for cognitive radio: challenges and solutions[J]. Eurasip Journal on Advances in Signal Processing, 2010, ID 381465.

    [4] NISHIMORI K, TARANTO R D, YOMO H, et al. Spatial opportunity for cognitive radio systems with heterogeneous path loss conditions[C]//IEEE 65th Vehicular Technology Conference VTC.c2007: 2631-2635.

    [5] RIIHIJARVI J, MAHONEN P. Exploiting spatial statistics of primary and secondary users towards improved cognitive radio networks[C]// IEEE 3rd International Conference on Cognitive Radio Oriented Wireless Networks and Communications CrownCom.c2008:1-7.

    [6] MIN A W, KIM K H, SINGH J P, et al. Opportunistic spectrum access for mobile cognitive radios[C]//IEEE INFOCOM Conference.c2011: 2993-3001.

    [7] CASO G, NARDIS L D, HOLLAND O, et al. Impact of spatio-temporal correlation in cooperative spectrum sensing for mobile cognitive radio networks[C]//The 10th International Symposium on Wireless Communication Systems ISWCS. c2013: 1-5.

    [8] PAURA L, SAVOIA R. Mobility-aware sensing enabled capacity in cognitive radio networks[C]//2013 IEEE International Workshop on Measurements and Networking Proceedings M&N.c2013:179-183.

    [9] LI F, XU Z B. Sparse bayesian hierarchical prior modeling based cooperative spectrum sensing in wideband cognitive radio networks[J]. IEEE Signal Process Lettes, 2014, 21(5): 586-590.

    [10] GIANNAKIS G B, TEPEDELENLIOGLU C. Basis expansion models and diversity techniques for blind identification and equalization of time-varying channels[J]. Proceedings of the IEEE, 1998, 86(10): 1969-1986.

    [11] CEVHER V, DUARTE M F, BARANIUK R G. Distributed target localization via spatial sparsity[C]//European Signal Processing Conference EUSIPCO.c2008.:1-5.

    [12] GIROLAMI M. A variational method for learning sparse and overcomplete representations[J]. Neural Computation, 2011, 13(11): 2517-2532.

    [13] WIPF D, OWEN J, ATTIAS H, et al. Robust Bayesian estimation of the location, orientation, and time course of multiple correlated neural sources using MEG[J]. NeuroImage, 2010, 49(1): 641-655.

    Novel cooperative globalspectrum sensing algorithm based on variational Bayesian inference

    WU Ming, SONG Tie-cheng, HU Jing, SHEN Lian-feng

    (National Mobile Communications Research Laboratory, Southeast University, Nanjing 210096, China)

    To realize multi-dimensional dynamic spectrum access, an approximate model was proposed for the global power spectral density (PSD) of primary users (PU). Based on the proposed model, a novel cooperative spectrum sensing algorithm was proposed, and its overall flow was also built to obtain global information in the network of PU. The global information included locations, occupied frequency bands and transmitting powers of the PU. Then, an estimator of model coefficient vector was designed by utilizing the theory of variational Bayesian inference (VBI). Simulation results show that the proposed approximate model has good accuracy, and the corresponding estimation algorithm of model coefficient vector has good convergence and stability. Meanwhile, the relationship between SNR and the leakage of aggregate spurious power (LASP) was pointed out, and the influence of SNR and LASP on MSE performance was also discussed. Furthermore, it is proved that the proposed algorithm has better MSE performance than another algorithm since the sparsity of model coefficient vector is utilized.

    cognitive radio, cooperative global spectrum sensing, variational Bayesian inference, sparsity

    TN914

    A

    10.11959/j.issn.1000-436x.2016037

    2015-04-08;

    2015-08-08

    國(guó)家自然科學(xué)基金資助項(xiàng)目(No.61271207, No.61372104, No.61201248)

    The National Natural Science Foundation of China (No.61271207, No.61372104, No.61201248)

    吳名(1981-),男,江蘇南京人,東南大學(xué)博士生,主要研究方向?yàn)檎J(rèn)知無(wú)線電系統(tǒng)、協(xié)作通信、傳感器網(wǎng)絡(luò)等。

    宋鐵成(1967-),男,江蘇張家港人,博士,東南大學(xué)教授、博士生導(dǎo)師,主要研究方向?yàn)橐苿?dòng)通信理論與技術(shù)、認(rèn)知無(wú)線電、物聯(lián)網(wǎng)等。

    胡靜(1975-),女,江蘇揚(yáng)州人,博士,東南大學(xué)副研究員,主要研究方向?yàn)槎叹嚯x無(wú)線通信、泛在網(wǎng)絡(luò)等。

    沈連豐(1952-),男,江蘇邳州人,東南大學(xué)教授、博士生導(dǎo)師,主要研究方向?yàn)閷拵б苿?dòng)通信、短距離無(wú)線通信與泛在網(wǎng)絡(luò)等。

    猜你喜歡
    后驗(yàn)全局頻段
    Cahn-Hilliard-Brinkman系統(tǒng)的全局吸引子
    量子Navier-Stokes方程弱解的全局存在性
    gPhone重力儀的面波頻段響應(yīng)實(shí)測(cè)研究
    地震研究(2021年1期)2021-04-13 01:04:56
    基于對(duì)偶理論的橢圓變分不等式的后驗(yàn)誤差分析(英)
    貝葉斯統(tǒng)計(jì)中單參數(shù)后驗(yàn)分布的精確計(jì)算方法
    落子山東,意在全局
    金橋(2018年4期)2018-09-26 02:24:54
    一種基于最大后驗(yàn)框架的聚類分析多基線干涉SAR高度重建算法
    推擠的5GHz頻段
    CHIP新電腦(2016年3期)2016-03-10 14:07:52
    TD—LTE在D頻段和F頻段的覆蓋能力差異
    新思路:牽一發(fā)動(dòng)全局
    村上凉子中文字幕在线| 少妇的逼好多水| 日韩av在线大香蕉| 亚洲一区高清亚洲精品| 天堂网av新在线| 色综合色国产| 99国产极品粉嫩在线观看| 国产精品一区二区三区四区久久| 免费人成视频x8x8入口观看| 亚洲一区高清亚洲精品| 尾随美女入室| 亚洲 国产 在线| 免费电影在线观看免费观看| 毛片一级片免费看久久久久 | 国产探花在线观看一区二区| 国产老妇女一区| 精品久久久久久,| 亚洲无线在线观看| 欧美色视频一区免费| 国产在视频线在精品| 熟妇人妻久久中文字幕3abv| 日日啪夜夜撸| 丰满乱子伦码专区| 简卡轻食公司| 中文字幕久久专区| 亚洲第一区二区三区不卡| 99久久精品热视频| 欧美在线一区亚洲| 99久国产av精品| 欧美日韩国产亚洲二区| 国产女主播在线喷水免费视频网站 | 精品久久久久久久久亚洲 | 国产高清不卡午夜福利| 精品国产三级普通话版| 最近中文字幕高清免费大全6 | 一区二区三区高清视频在线| 深夜精品福利| 男人的好看免费观看在线视频| 国内久久婷婷六月综合欲色啪| 韩国av一区二区三区四区| 成人国产一区最新在线观看| 国产一区二区亚洲精品在线观看| 亚洲国产色片| 最近最新中文字幕大全电影3| 精品福利观看| 男女下面进入的视频免费午夜| 男女之事视频高清在线观看| 欧美3d第一页| 国产精品亚洲美女久久久| 国产男靠女视频免费网站| 丰满乱子伦码专区| 韩国av在线不卡| 波多野结衣巨乳人妻| 22中文网久久字幕| 麻豆国产97在线/欧美| 中亚洲国语对白在线视频| 亚洲人成网站在线播| 人妻制服诱惑在线中文字幕| 亚洲欧美日韩无卡精品| 一夜夜www| 久久精品国产99精品国产亚洲性色| 在线观看舔阴道视频| 神马国产精品三级电影在线观看| 亚洲精品粉嫩美女一区| 能在线免费观看的黄片| 久久久久久久久中文| 男女啪啪激烈高潮av片| 成人特级av手机在线观看| 亚洲人与动物交配视频| 少妇的逼水好多| 波多野结衣高清无吗| 嫩草影院精品99| aaaaa片日本免费| 免费av观看视频| 国产精品久久久久久亚洲av鲁大| 九九热线精品视视频播放| 18禁黄网站禁片午夜丰满| 成人欧美大片| 日本黄大片高清| 欧美黑人欧美精品刺激| 日本欧美国产在线视频| 大又大粗又爽又黄少妇毛片口| ponron亚洲| 国产精品日韩av在线免费观看| 日韩 亚洲 欧美在线| 午夜视频国产福利| 久久久午夜欧美精品| 国产一区二区三区在线臀色熟女| 桃色一区二区三区在线观看| 内射极品少妇av片p| 97人妻精品一区二区三区麻豆| 一夜夜www| 亚洲人成网站在线播| 亚洲专区国产一区二区| 我的女老师完整版在线观看| 99久久精品国产国产毛片| 久久久久久久亚洲中文字幕| 亚洲成人久久性| 床上黄色一级片| 一进一出抽搐动态| 日日摸夜夜添夜夜添av毛片 | 蜜桃久久精品国产亚洲av| 在线a可以看的网站| 91麻豆精品激情在线观看国产| 男人的好看免费观看在线视频| 欧美一级a爱片免费观看看| av在线蜜桃| 精品久久国产蜜桃| 国产中年淑女户外野战色| 在线观看舔阴道视频| 搡老岳熟女国产| 亚洲aⅴ乱码一区二区在线播放| 亚洲人成网站在线播| 国产视频一区二区在线看| 日韩欧美三级三区| 国产精品一区www在线观看 | 欧美xxxx黑人xx丫x性爽| 在线免费观看的www视频| 桃红色精品国产亚洲av| 免费在线观看日本一区| 日韩精品青青久久久久久| 亚洲性夜色夜夜综合| 99热这里只有是精品在线观看| 亚洲欧美日韩东京热| 日韩精品有码人妻一区| 99riav亚洲国产免费| 老女人水多毛片| av在线蜜桃| a级毛片a级免费在线| 麻豆av噜噜一区二区三区| 少妇熟女aⅴ在线视频| 少妇丰满av| 成人无遮挡网站| 天堂网av新在线| 亚洲精品亚洲一区二区| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| av在线蜜桃| 乱人视频在线观看| 亚洲av电影不卡..在线观看| 国产老妇女一区| 精品久久久久久久久av| 麻豆av噜噜一区二区三区| 免费无遮挡裸体视频| 色吧在线观看| 国产成人aa在线观看| 赤兔流量卡办理| 国产午夜精品久久久久久一区二区三区 | 国产成人福利小说| 欧美丝袜亚洲另类 | 成年人黄色毛片网站| 美女被艹到高潮喷水动态| 精品人妻1区二区| 免费大片18禁| 日本在线视频免费播放| 麻豆久久精品国产亚洲av| 国产精品自产拍在线观看55亚洲| 特大巨黑吊av在线直播| 99视频精品全部免费 在线| 床上黄色一级片| 国产单亲对白刺激| 小说图片视频综合网站| 久久久久久久午夜电影| 婷婷丁香在线五月| 免费观看在线日韩| 99九九线精品视频在线观看视频| 亚洲avbb在线观看| 精品一区二区三区人妻视频| 18禁在线播放成人免费| 高清在线国产一区| 久久热精品热| 日韩av在线大香蕉| 国产欧美日韩精品一区二区| 亚洲自拍偷在线| 色精品久久人妻99蜜桃| 亚洲av中文字字幕乱码综合| 日韩精品有码人妻一区| 99精品久久久久人妻精品| 日日摸夜夜添夜夜添av毛片 | av专区在线播放| 亚洲美女黄片视频| 日韩欧美精品免费久久| 精品久久久久久久人妻蜜臀av| 少妇人妻一区二区三区视频| 亚洲va在线va天堂va国产| 看片在线看免费视频| 美女大奶头视频| a在线观看视频网站| 天堂av国产一区二区熟女人妻| 无遮挡黄片免费观看| 亚洲va在线va天堂va国产| 午夜日韩欧美国产| 精品国产三级普通话版| 国产成人av教育| 国产精品一区二区三区四区免费观看 | 色综合站精品国产| 亚洲美女搞黄在线观看 | 亚洲精品粉嫩美女一区| 精品久久久久久成人av| 日韩一本色道免费dvd| 91在线精品国自产拍蜜月| 成人美女网站在线观看视频| 少妇人妻一区二区三区视频| 真人一进一出gif抽搐免费| 91狼人影院| 香蕉av资源在线| 无人区码免费观看不卡| 波野结衣二区三区在线| 亚洲av日韩精品久久久久久密| 免费黄网站久久成人精品| 国产爱豆传媒在线观看| 中亚洲国语对白在线视频| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 男人和女人高潮做爰伦理| 色5月婷婷丁香| 99热网站在线观看| 桃色一区二区三区在线观看| 亚洲第一区二区三区不卡| 欧美日韩乱码在线| 国产大屁股一区二区在线视频| 日韩欧美国产一区二区入口| 一边摸一边抽搐一进一小说| 伦精品一区二区三区| 深夜精品福利| 色av中文字幕| 22中文网久久字幕| 嫩草影院新地址| 亚洲欧美日韩高清在线视频| eeuss影院久久| 国产淫片久久久久久久久| 天天一区二区日本电影三级| 女生性感内裤真人,穿戴方法视频| 免费看日本二区| 亚洲第一区二区三区不卡| 国产黄色小视频在线观看| 真实男女啪啪啪动态图| 日韩中字成人| 51国产日韩欧美| 网址你懂的国产日韩在线| 国产一区二区亚洲精品在线观看| 国产在视频线在精品| 亚洲乱码一区二区免费版| 午夜免费激情av| 尾随美女入室| 哪里可以看免费的av片| 蜜桃久久精品国产亚洲av| 波多野结衣高清无吗| a在线观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 免费高清视频大片| 免费黄网站久久成人精品| 久久久午夜欧美精品| 日本一二三区视频观看| 国产探花极品一区二区| 国产国拍精品亚洲av在线观看| 又爽又黄无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 最近中文字幕高清免费大全6 | 久久精品国产清高在天天线| 婷婷丁香在线五月| 午夜福利在线在线| 亚洲四区av| 久久久久国产精品人妻aⅴ院| 久久久久久久午夜电影| 国产在线精品亚洲第一网站| 最新中文字幕久久久久| 精品国产三级普通话版| 麻豆av噜噜一区二区三区| 亚洲成人久久性| 欧美成人一区二区免费高清观看| 搡老熟女国产l中国老女人| 国产亚洲91精品色在线| 尤物成人国产欧美一区二区三区| 欧美又色又爽又黄视频| 久久久国产成人免费| 国内精品一区二区在线观看| 淫秽高清视频在线观看| 天堂网av新在线| 久久久成人免费电影| 麻豆av噜噜一区二区三区| 色哟哟·www| 日本三级黄在线观看| av在线观看视频网站免费| 亚洲欧美日韩高清专用| 国产高清不卡午夜福利| 欧美激情在线99| 尾随美女入室| 色综合站精品国产| 亚洲三级黄色毛片| 久久婷婷人人爽人人干人人爱| 国产成人aa在线观看| 动漫黄色视频在线观看| 天天躁日日操中文字幕| 看免费成人av毛片| 免费观看人在逋| 免费一级毛片在线播放高清视频| 国产精品国产三级国产av玫瑰| 精品99又大又爽又粗少妇毛片 | 别揉我奶头 嗯啊视频| 精品人妻视频免费看| 能在线免费观看的黄片| 国产一区二区三区视频了| 男女视频在线观看网站免费| 老司机深夜福利视频在线观看| 久久久久久久午夜电影| 日日夜夜操网爽| 国产精品,欧美在线| 蜜桃亚洲精品一区二区三区| 在线观看av片永久免费下载| 国产亚洲av嫩草精品影院| 免费av不卡在线播放| 日本免费a在线| 真人做人爱边吃奶动态| 丰满人妻一区二区三区视频av| 亚洲av不卡在线观看| 久久欧美精品欧美久久欧美| 在线免费十八禁| 成年女人看的毛片在线观看| 一本久久中文字幕| 变态另类成人亚洲欧美熟女| 嫩草影院入口| 成人国产综合亚洲| 99国产精品一区二区蜜桃av| 内射极品少妇av片p| 最近视频中文字幕2019在线8| 国内精品久久久久精免费| 精品国产三级普通话版| 欧美另类亚洲清纯唯美| www日本黄色视频网| 欧美3d第一页| 午夜福利18| 在线看三级毛片| 国内精品美女久久久久久| 亚洲精品久久国产高清桃花| 又紧又爽又黄一区二区| 一边摸一边抽搐一进一小说| 又紧又爽又黄一区二区| 欧美成人一区二区免费高清观看| 免费av不卡在线播放| 99久久无色码亚洲精品果冻| 欧美激情久久久久久爽电影| 日日干狠狠操夜夜爽| 国产免费男女视频| 日本a在线网址| 国产免费男女视频| 精品免费久久久久久久清纯| 伦精品一区二区三区| 欧美性猛交╳xxx乱大交人| 99久久精品一区二区三区| 国产高清视频在线观看网站| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 黄色丝袜av网址大全| 亚洲国产精品sss在线观看| 精品欧美国产一区二区三| 免费看光身美女| 欧美日韩中文字幕国产精品一区二区三区| 变态另类丝袜制服| 久久久久久久久大av| 欧美成人免费av一区二区三区| av天堂在线播放| 国产精品嫩草影院av在线观看 | 2021天堂中文幕一二区在线观| 中亚洲国语对白在线视频| 国产午夜精品论理片| 国产一区二区激情短视频| 国产色爽女视频免费观看| 精品久久久久久久久亚洲 | 22中文网久久字幕| 男女视频在线观看网站免费| 波多野结衣巨乳人妻| 国产高清激情床上av| 国产淫片久久久久久久久| 欧美激情久久久久久爽电影| 有码 亚洲区| 国产精品永久免费网站| 欧美又色又爽又黄视频| 午夜免费激情av| 免费看光身美女| 久久精品综合一区二区三区| 黄色丝袜av网址大全| 色在线成人网| 男人和女人高潮做爰伦理| 国产精品一区二区三区四区免费观看 | 国产av在哪里看| 国内精品宾馆在线| 黄色配什么色好看| 国产精品自产拍在线观看55亚洲| 亚洲av成人精品一区久久| 成人三级黄色视频| 婷婷丁香在线五月| 九九久久精品国产亚洲av麻豆| 高清在线国产一区| 亚洲精品色激情综合| 日本一本二区三区精品| 有码 亚洲区| 色综合婷婷激情| 狂野欧美激情性xxxx在线观看| 人人妻人人看人人澡| 噜噜噜噜噜久久久久久91| 国产亚洲精品综合一区在线观看| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 我要搜黄色片| 国语自产精品视频在线第100页| 性插视频无遮挡在线免费观看| 日本成人三级电影网站| 免费观看的影片在线观看| 国产极品精品免费视频能看的| 精品一区二区三区视频在线| 成年女人看的毛片在线观看| 亚洲人与动物交配视频| 深夜a级毛片| 极品教师在线视频| 日韩精品有码人妻一区| 久久人妻av系列| 精品午夜福利在线看| 啦啦啦韩国在线观看视频| 男女边吃奶边做爰视频| av在线老鸭窝| 精品久久久久久久久亚洲 | 免费观看的影片在线观看| 精品国内亚洲2022精品成人| 搡女人真爽免费视频火全软件 | 超碰av人人做人人爽久久| 国产单亲对白刺激| 久久久久久九九精品二区国产| 毛片一级片免费看久久久久 | 成人av一区二区三区在线看| 国模一区二区三区四区视频| netflix在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 欧美色欧美亚洲另类二区| 身体一侧抽搐| 国产中年淑女户外野战色| 男人舔女人下体高潮全视频| 级片在线观看| 久久久久久国产a免费观看| 亚洲无线观看免费| 久久久久久久精品吃奶| 亚洲第一区二区三区不卡| 日韩 亚洲 欧美在线| 在线观看av片永久免费下载| 蜜桃亚洲精品一区二区三区| 国产v大片淫在线免费观看| www日本黄色视频网| 尤物成人国产欧美一区二区三区| 久久人人爽人人爽人人片va| 亚洲无线观看免费| 日本 欧美在线| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 亚洲精品一区av在线观看| 亚洲av电影不卡..在线观看| 美女高潮喷水抽搐中文字幕| 亚洲午夜理论影院| 啪啪无遮挡十八禁网站| 欧美三级亚洲精品| 久久精品国产亚洲网站| avwww免费| 最近最新中文字幕大全电影3| 精品久久国产蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 夜夜爽天天搞| 欧美日韩黄片免| 狂野欧美白嫩少妇大欣赏| 日本色播在线视频| 真实男女啪啪啪动态图| 亚洲最大成人av| 日本熟妇午夜| 亚洲四区av| 99国产极品粉嫩在线观看| 欧美日韩亚洲国产一区二区在线观看| 深夜精品福利| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| 深夜精品福利| 级片在线观看| avwww免费| 亚洲图色成人| 欧美日韩精品成人综合77777| 99精品久久久久人妻精品| 人妻夜夜爽99麻豆av| 美女cb高潮喷水在线观看| 亚洲av中文av极速乱 | 国产精品一及| 中文字幕高清在线视频| 少妇猛男粗大的猛烈进出视频 | 在线观看免费视频日本深夜| 能在线免费观看的黄片| eeuss影院久久| 淫妇啪啪啪对白视频| 国产一区二区三区在线臀色熟女| 亚洲av二区三区四区| 欧美色视频一区免费| 深夜a级毛片| 最近最新中文字幕大全电影3| 国产精品爽爽va在线观看网站| 伦理电影大哥的女人| 午夜精品在线福利| 成人三级黄色视频| 成人美女网站在线观看视频| 欧美日韩乱码在线| 国产午夜福利久久久久久| 两人在一起打扑克的视频| 精品国产三级普通话版| 天堂av国产一区二区熟女人妻| 一区二区三区四区激情视频 | 变态另类成人亚洲欧美熟女| 三级毛片av免费| 久久精品国产自在天天线| 成人综合一区亚洲| 一区二区三区四区激情视频 | 小说图片视频综合网站| 一级黄片播放器| 久久99热这里只有精品18| 成人精品一区二区免费| 免费搜索国产男女视频| 国产精品综合久久久久久久免费| 亚洲avbb在线观看| 久久精品综合一区二区三区| 一个人看的www免费观看视频| 麻豆成人午夜福利视频| 男女下面进入的视频免费午夜| 啦啦啦韩国在线观看视频| 欧美成人性av电影在线观看| 日本黄大片高清| 蜜桃亚洲精品一区二区三区| 十八禁网站免费在线| 国产伦在线观看视频一区| 最后的刺客免费高清国语| 日本五十路高清| 高清在线国产一区| 听说在线观看完整版免费高清| 女同久久另类99精品国产91| 久久精品国产亚洲网站| 免费电影在线观看免费观看| 欧美+日韩+精品| 午夜免费激情av| 亚洲成人精品中文字幕电影| 亚洲在线观看片| 日本免费a在线| 国产一区二区三区视频了| 极品教师在线视频| 有码 亚洲区| h日本视频在线播放| 美女 人体艺术 gogo| 亚洲美女视频黄频| 国产精品不卡视频一区二区| 99久久精品国产国产毛片| 美女大奶头视频| 99热精品在线国产| 天天躁日日操中文字幕| 99国产极品粉嫩在线观看| 国产 一区 欧美 日韩| 日本免费a在线| 国产精品一区二区性色av| 极品教师在线视频| 成年女人永久免费观看视频| 免费av不卡在线播放| xxxwww97欧美| 国产精品一区二区三区四区久久| 日韩精品中文字幕看吧| 别揉我奶头~嗯~啊~动态视频| 久久精品国产自在天天线| 看片在线看免费视频| 国模一区二区三区四区视频| 亚洲美女黄片视频| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品成人综合色| 亚洲av中文av极速乱 | 听说在线观看完整版免费高清| 日韩中文字幕欧美一区二区| 极品教师在线视频| 色综合站精品国产| 在线天堂最新版资源| 中出人妻视频一区二区| 国产一区二区亚洲精品在线观看| 又紧又爽又黄一区二区| 欧美中文日本在线观看视频| 亚洲内射少妇av| 三级男女做爰猛烈吃奶摸视频| 啦啦啦啦在线视频资源| 哪里可以看免费的av片| 欧美最新免费一区二区三区| 欧美丝袜亚洲另类 | 成人永久免费在线观看视频| 欧美色视频一区免费| 男女做爰动态图高潮gif福利片| 中文资源天堂在线| 国产欧美日韩一区二区精品| 久久精品国产亚洲av香蕉五月| 精品久久久久久久人妻蜜臀av| 制服丝袜大香蕉在线| 欧美日韩综合久久久久久 | 在线观看66精品国产| 国产淫片久久久久久久久| 久久人人精品亚洲av| av黄色大香蕉| 99九九线精品视频在线观看视频| 我的老师免费观看完整版| 国内久久婷婷六月综合欲色啪| 在线观看美女被高潮喷水网站| 国产精品综合久久久久久久免费| 精品久久久久久久久久久久久| 日韩欧美免费精品| 99久久久亚洲精品蜜臀av| 亚洲精品乱码久久久v下载方式| 最近视频中文字幕2019在线8| 精品一区二区三区人妻视频| 亚洲第一电影网av| 我要搜黄色片| 91午夜精品亚洲一区二区三区 |