999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于深度學(xué)習的調(diào)制識別綜述

2021-06-20 03:54:52孫姝君彭盛亮姚育東楊喜
電信科學(xué) 2021年5期
關(guān)鍵詞:特征信號

孫姝君,彭盛亮,姚育東,楊喜

(1. 華僑大學(xué)信息科學(xué)與工程學(xué)院,福建 廈門 361021;2. 史蒂文斯理工學(xué)院電子與計算機工程系,美國 新澤西州 霍博肯 07030;3. 吉首大學(xué)信息科學(xué)與工程學(xué)院,湖南 吉首 416000)

1 引言

調(diào)制識別是要識別出信號采用的調(diào)制方案,有著廣泛的應(yīng)用。在認知無線電中,調(diào)制識別有助于次級用戶了解主用戶信號的細節(jié)信息,是無線環(huán)境全面感知的重要組成部分;在智能通信中,調(diào)制識別能夠減少調(diào)制信息的傳輸開銷,實現(xiàn)信號的盲解調(diào);在無線電監(jiān)管中,調(diào)制識別有利于更精準地判斷信號是否滿足無線電管理規(guī)定;在電子對抗中,識別敵方信號的調(diào)制方案是有效實施無線電干擾和欺騙的前提[1]。

近年來,隨著深度學(xué)習(deep learning, DL)的蓬勃發(fā)展,基于DL的調(diào)制識別技術(shù)越來越受到人們的關(guān)注。與傳統(tǒng)調(diào)制識別方法不同,基于DL的方法采用深度神經(jīng)網(wǎng)絡(luò)(deep neural network,DNN)作為分類器,具有3方面的優(yōu)勢:第一,能夠充分利用通信系統(tǒng)中存在的海量數(shù)據(jù),有效改善調(diào)制識別精度;第二,能夠?qū)崿F(xiàn)自動提取特征,克服手動提取特征過程依賴專業(yè)背景和人工經(jīng)驗的難題;第三,能夠持續(xù)受益于深度學(xué)習工具的快速迭代和演進,具備在復(fù)雜場景下解決復(fù)雜調(diào)制識別問題的潛能[2]。

2 相關(guān)基礎(chǔ)

2.1 傳統(tǒng)算法

傳統(tǒng)的調(diào)制識別算法可分為兩大類,即基于似然(likelihood based,LB)的算法和基于特征(feature based,F(xiàn)B)的算法[1]。LB算法將調(diào)制識別問題視為一個多元假設(shè)檢驗問題,它可以得到最優(yōu)解,但其計算復(fù)雜度相對較高。參考文獻[3]研究了平均似然比和廣義似然比在調(diào)制識別中的應(yīng)用。參考文獻[4]設(shè)計了兩種分別基于Gauss-Legendre和Gauss-Hermite求積規(guī)則的LB算法,用于識別線性調(diào)制。FB算法從信號中提取特征用以區(qū)別不同的調(diào)制方案,如果特征和分類器選擇恰當,該算法可以在保證較低復(fù)雜度的同時獲得近似最優(yōu)的性能,因而逐漸占據(jù)了主導(dǎo)地位[5]。參考文獻[6]提出了一種以高階累積量(high-order cumulant,HOC)為特征的FB算法,在多徑衰落條件下具有較好的魯棒性。參考文獻[7]以近似熵為特征,實現(xiàn)了連續(xù)相位調(diào)制的識別。在參考文獻[8]中,四階譜特征和壓縮感知被用來區(qū)分星座圖類似的調(diào)制信號。此外,在分類器方面,參考文獻[9]的研究表明,決策樹比人工神經(jīng)網(wǎng)絡(luò)更適用于多輸入多輸出系統(tǒng)。參考文獻[10-12]中也提出了基于其他分類器的FB算法,如支持向量機和K最近鄰等。

2.2 深度學(xué)習

DL是當前人工智能浪潮中強有力的工具之一,其本質(zhì)是一種基于多層神經(jīng)網(wǎng)絡(luò)、以海量數(shù)據(jù)為輸入、自動組合低層特征形成高層特征的規(guī)則學(xué)習方法。最基本的神經(jīng)網(wǎng)絡(luò)包括3部分,分別為輸入層、隱藏層和輸出層[13]。其中,輸入層的神經(jīng)元用于接收輸入數(shù)據(jù),并將數(shù)據(jù)傳輸?shù)诫[藏層的各神經(jīng)元;隱藏層中的神經(jīng)元利用激活函數(shù)處理數(shù)據(jù),并將處理后的數(shù)據(jù)傳到輸出層;最后由輸出層的神經(jīng)元輸出結(jié)果。在DNN中,隱藏層的層數(shù)越多,結(jié)構(gòu)越復(fù)雜。此時,整個網(wǎng)絡(luò)的參數(shù)也變得更多,每一層相對上一層的抽象表示也更加深入,因而DNN相較于基本的神經(jīng)網(wǎng)絡(luò)具有更強大的性能。

近年來,隨著DL的蓬勃發(fā)展,DNN也被引入通信領(lǐng)域,用于解決調(diào)制識別等傳統(tǒng)通信問題。基于DL的調(diào)制識別技術(shù),能夠充分發(fā)揮DNN的性能優(yōu)勢,彌補原有LB和FB算法的不足。

3 基于深度學(xué)習的調(diào)制識別

3.1 系統(tǒng)架構(gòu)

在通信系統(tǒng)中,接收信號通常可表示為如下形式:

其中,n(k)代表加性高斯白噪聲(additive white Gaussian noise,AWGN),x(l)代表符號序列,A代表幅度因子,f0代表載波頻率偏移,kθ代表相位抖動,T代表符號間隔,h(·)代表殘留信道影響,T∈為定時誤差。

考慮由N個接收信號采樣點組成的接收向量和M種調(diào)制方案組成的候選集合,調(diào)制識別的任務(wù)是根據(jù)接收向量確定調(diào)制方案Si。

基于DL的調(diào)制識別的系統(tǒng)結(jié)構(gòu)如圖1所示。接收到的調(diào)制信號首先經(jīng)過數(shù)據(jù)預(yù)處理,轉(zhuǎn)換為適當?shù)男盘柋碚鳎恍盘柋碚髟偎腿隓NN中處理,并完成調(diào)制方案的識別。在實際應(yīng)用中,根據(jù)調(diào)制候選集和無線環(huán)境的不同,有必要采用適當?shù)臄?shù)據(jù)預(yù)處理方法,設(shè)計合理的DNN結(jié)構(gòu),以滿足不同的應(yīng)用需求。

圖1 基于DL的調(diào)制識別系統(tǒng)結(jié)構(gòu)

3.2 數(shù)據(jù)預(yù)處理

數(shù)據(jù)預(yù)處理作為基于DL的調(diào)制識別中的關(guān)鍵步驟,目的在于對原始數(shù)據(jù)進行變換,將其表征成恰當?shù)男问剑怪軌蜻m于DNN處理。如果數(shù)據(jù)預(yù)處理方式得當,DNN就能夠更好地從數(shù)據(jù)中提取典型特征,達到更佳的識別性能。根據(jù)預(yù)處理后數(shù)據(jù)表征形式的不同,數(shù)據(jù)預(yù)處理方式可以分為三大類:特征表征、序列表征和圖像表征。

特征表征是將接收信號預(yù)處理成一個或者多個特征值,例如HOC特征、高階矩特征、頻譜特征及其他統(tǒng)計特征等[1]。參考文獻[6]將接收信號轉(zhuǎn)換成HOC特征,以對抗多徑衰落的影響;參考文獻[10]使用四階譜特征表征接收信號,實現(xiàn)了4種調(diào)相信號的識別;參考文獻[14]研究了多達28種特征表征,并討論了不同特征值組合對調(diào)制識別性能的影響。

序列表征是將接收信號預(yù)處理成一個一維向量,例如正交同相(in-phase and quadrature,IQ)序列[15-17]、幅度相位(amplitude and phase,AP)序列[18-19]、快速傅里葉變換(fast Fourier transform,F(xiàn)FT)序列[20]和幅值直方圖序列[21]等。參考文獻[22]最早將IQ表征用于調(diào)制識別,性能較傳統(tǒng)方法具有較大優(yōu)勢;參考文獻[16]利用IQ序列表征信號,討論了不確定噪聲條件下的調(diào)制識別問題;參考文獻[18]分析了IQ、AP和FFT 3種序列表征,同時指出在不同的信噪比(signal-to-noise ratio,SNR)條件下,不同的表征形式具有不同的調(diào)制識別精度。

圖像表征是將接收信號預(yù)處理成一個二維矩陣,例如星座圖[2,5,23]、眼圖[24]、特征點圖[25]以及經(jīng)各種變換后產(chǎn)生的譜圖[26-28]等。參考文獻[2]利用接收信號的星座圖表征,將調(diào)制識別問題轉(zhuǎn)換成圖像識別問題,并利用DNN解決;參考文獻[27]和參考文獻[28]分別將接收信號轉(zhuǎn)換成雙譜圖和循環(huán)熵譜圖,也取得了較高的識別精度;參考文獻文獻[29]采用等勢星座圖替代原始的星座圖,進一步提升了調(diào)制識別的性能。

3.3 深度神經(jīng)網(wǎng)絡(luò)

DNN是基于DL的調(diào)制識別的核心部分,它負責處理信號表征結(jié)果,并完成調(diào)制方案的推斷和輸出。DNN一般分為三大類,即深度前饋網(wǎng)絡(luò)、卷積神經(jīng)網(wǎng)絡(luò)(convolutional neural network,CNN)和循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural network,RNN)[13]。

3.3.1 深度前饋網(wǎng)絡(luò)

深度前饋網(wǎng)絡(luò),也叫前饋神經(jīng)網(wǎng)絡(luò)或多層感知機。深度前饋網(wǎng)絡(luò)的特點是整個網(wǎng)絡(luò)中無反饋,信號從輸入層到輸出層單向傳播。其每層神經(jīng)元與下一層神經(jīng)元之間完全互連,各神經(jīng)元之間不存在同層連接和跨層連接。深度置信網(wǎng)絡(luò)(deep belief network,DBN)是深度前饋網(wǎng)絡(luò)的典型代表之一,由多層受限玻爾茲曼機(restricted Boltzmann machine,RBM)組成[30]。其中每個RBM有兩層,即上層隱藏層和下層可見層。

最近幾年,國內(nèi)外的研究者已經(jīng)提出了不少基于DBN的調(diào)制識別算法。參考文獻[31]設(shè)計了一種包括3個常規(guī)RBM和1個高斯-伯努利RBM的DBN,用于處理接收信號的譜圖表征,在SNR大于或等于2 dB時,對5種調(diào)制方案的識別精度達到了90%以上。參考文獻[32]提出了一種基于DBN和星座圖投影的調(diào)制識別算法,利用DBN學(xué)習星座圖投影的潛在特征,其識別精度接近LB算法,但復(fù)雜度更低。

3.3.2 卷積神經(jīng)網(wǎng)絡(luò)

CNN通常指包含卷積運算的神經(jīng)網(wǎng)絡(luò)。在CNN中,相鄰層局部連通,多層之間利用空間局部相關(guān)性實現(xiàn)權(quán)重共享,因此同一個參數(shù)可以被不同的神經(jīng)元多次使用。此外,隱藏層內(nèi)的卷積核使得CNN的輸入層與輸出層之間具有稀疏連通性,從而大大提高了計算效率,降低了網(wǎng)絡(luò)存儲需求[33]。代表性的CNN結(jié)構(gòu)有AlexNet、VGG-Net、GoogLeNet、ResNet等。

AlexNet包含5個卷積層和3個全連接層,其中卷積層之后還連接了最大池化層,一共有6 000萬個參數(shù)和65萬個神經(jīng)元。參考文獻[2]提出了一種基于AlexNet和星座圖表征的調(diào)制識別算法,利用AlexNet學(xué)習星座圖的特征,完成了8種調(diào)制方案的識別。參考文獻[28]提出了一種基于AlexNet和譜圖的調(diào)制識別算法,將譜圖送入AlexNet進行訓(xùn)練,并用訓(xùn)練好的AlexNet模型推斷接收信號的調(diào)制方案。為了解決數(shù)據(jù)不足導(dǎo)致的調(diào)制識別性能惡化的問題,參考文獻[34]在AlexNet的基礎(chǔ)上,引入生成對抗網(wǎng)絡(luò)對訓(xùn)練數(shù)據(jù)進行擴充,使得整體識別精度提高了0.1%~6%。為了便于在邊緣設(shè)備和實時場景中部署基于DL的調(diào)制識別技術(shù),參考文獻[29]提出采用網(wǎng)絡(luò)剪枝技術(shù)對AlexNet進行裁剪。與原始的AlexNet相比,剪枝后網(wǎng)絡(luò)僅使用了1.5%~5%的參數(shù),僅耗費了33%~35%的識別時間。

VGG-Net由牛津大學(xué)計算機視覺組和Google DeepMind公司聯(lián)合提出,通過反復(fù)堆疊大小為3×3的小型卷積核和2×2的最大池化層,構(gòu)筑了一種16~19層深的CNN結(jié)構(gòu),具有更強的泛化能力和較少的迭代次數(shù)。參考文獻[15]提出使用VGG-Net對接收信號的IQ序列進行學(xué)習和識別,在AWGN信道下24種調(diào)制方案的平均識別精度約為33.5%,其中11種調(diào)制方案的平均識別精度約為54%。

自2012年AlexNet取得歷史性突破直到GoogLeNet橫空出世以前,主流的CNN結(jié)構(gòu)都是朝著更深、更寬的方向發(fā)展,但增大網(wǎng)絡(luò)也會使得網(wǎng)絡(luò)計算復(fù)雜度急劇增加。GoogLeNet由若干Inception模塊堆疊而成,采用全局平均池化層代替了最后的全連接層。這種結(jié)構(gòu)在增加CNN的深度和寬度的同時,保證了網(wǎng)絡(luò)計算量基本恒定。參考文獻[2]研究了基于GoogLeNet的調(diào)制識別,對8種調(diào)制方案的平均識別精度約為89%。

傳統(tǒng)的CNN在信息傳遞時或多或少地存在信息丟失和損耗的問題,同時還可能存在梯度消失或梯度爆炸的風險,導(dǎo)致網(wǎng)絡(luò)無法訓(xùn)練。ResNet在一定程度上解決了這個問題。它通過在網(wǎng)絡(luò)中增加直連通道,將輸入信息同時繞道傳送給輸出,既保護了信息的完整性,又降低了網(wǎng)絡(luò)訓(xùn)練的難度,使得構(gòu)建超深層網(wǎng)絡(luò)成為現(xiàn)實。參考文獻[27]就提出了一種基于ResNet和譜圖的調(diào)制識別算法,用于應(yīng)對非高斯噪聲影響。該算法在不同脈沖噪聲條件下均能保持較好的識別性能。參考文獻[17]在ResNet的基礎(chǔ)上,引入了3種融合方法用于解決實際接收信號長度不固定的問題。

3.3.3 循環(huán)神經(jīng)網(wǎng)絡(luò)

RNN是一種擅于處理序列數(shù)據(jù)的神經(jīng)網(wǎng)絡(luò)。在RNN中,網(wǎng)絡(luò)會存儲之前的信息并用于計算當前的輸出,即隱藏層的輸入不僅包括上一層的輸出,還包括上一時刻自身的輸出。當RNN的相關(guān)輸入與期望輸出在時間軸上間隔很長時,網(wǎng)絡(luò)訓(xùn)練容易產(chǎn)生梯度消失或者梯度爆炸,使得RNN的長時記憶失效。為了解決這個問題,長短期記憶(long short-term memory,LSTM)網(wǎng)絡(luò)應(yīng)運而生。LSTM網(wǎng)絡(luò)將數(shù)據(jù)維護在獨特的cell結(jié)構(gòu)中,通過輸入門、遺忘門和輸出門這3個門結(jié)構(gòu)實現(xiàn)信息的添加和刪除,以克服傳統(tǒng)RNN存在的長期依賴問題[35]。

參考文獻[16]設(shè)計了一種由3個LSTM層堆疊的RNN,用于對抗調(diào)制識別中的瑞利衰落和不確定噪聲影響,并增強頻率偏移場景下的魯棒性。參考文獻[36]提出了一種基于LSTM網(wǎng)絡(luò)的調(diào)制識別算法,用來學(xué)習接收信號的IQ序列和HOC表征,相較于CNN,它具有更好的識別性能。參考文獻[18]利用LSTM網(wǎng)絡(luò)處理接收信號的AP序列和FFT序列表征,在衰落信道下11種調(diào)制方案的平均識別精度接近90%。

3.4 常用數(shù)據(jù)集

為便于評估和對比不同調(diào)制識別算法的性能,研究者們創(chuàng)建了一些公開數(shù)據(jù)集,其中最為流行的是RadioML數(shù)據(jù)集[15,37],包含RadioML2016.04c、RadioML2016.10a、RadioML 2016.10b和RadioML2018.01a等多個版本。

RadioML2016.04c是一個利用GNU Radio生成的數(shù)據(jù)集,共涉及11種常用信號調(diào)制方案,包含8種數(shù)字調(diào)制和3種模擬調(diào)制,分別是二進制相移鍵控(binary phase shift keying,BPSK)、正交相移鍵控(quadrature phase shift keying,QPSK)、8PSK、十六進制正交幅度調(diào)制(16quadrature amplitude modulation,16QAM)、64QAM、二進制頻移鍵控(binary frequency shift keying,BFSK)、連續(xù)相位頻移鍵控(continuous phase frequency shift keying , CPFSK)、四進制脈沖幅度調(diào)制(pulse amplitude modulation 4,PAM4)、寬帶頻率調(diào)制(wide band frequency modulation,WBFM)、單邊帶調(diào)幅(single-sideband amplitude modulation,AM-SSB)、雙邊帶調(diào)幅(double-sideband amplitude modulation,AM-DSB)。其采樣速率為每符號8個采樣點,包括220 000個調(diào)制信號樣本,信道類型為多徑衰落信道。此外,它還是一個具有中等時鐘漂移、輕度衰落的可變SNR數(shù)據(jù)集,可用于不同信號和噪聲場景下的調(diào)制識別性能測試。

RadioML2016.10a是RadioML2016.04c的一個更干凈規(guī)范的版本,同樣包含了上述11種調(diào)制方案,信噪比變化范圍為-20~18 dB。該數(shù)據(jù)集考慮了AWGN和多徑衰落影響,并引入了采樣率偏移和中心頻率偏移等惡劣條件,以與實際無線通信環(huán)境更為相似。

RadioML2016.10b是RadioML2016.10a的擴展版本,一共包含1 200 000個信號樣本,但只涉及10種調(diào)制方案,分別為BPSK、QPSK、8PSK、16QAM、64QAM、BFSK、CPFSK、PAM4、WB-FM和 AM-DSB。

RadioML2018.01a對RadioML2016.10a中的調(diào)制類型進行了增廣,共涉及24種調(diào)制方案,包括通斷鍵控(on-off keying,OOK)、四進制幅移鍵控(4 amplitude shift keying,4ASK)、8ASK、BPSK、QPSK、偏移QPSK 、8PSK、16PSK、32PSK、十六進制幅度相移鍵控(16 amplitude and phase shift keying,16APSK)、32APSK、64APSK、128APSK、16QAM、32QAM、64QAM、128QAM、256QAM、單邊帶載波調(diào)幅(single-sideband amplitude modulation with carrier,AM-SSB-WC)、單邊帶抑制載波調(diào)幅(single-sideband amplitude modulation with suppressed-carrier,AM-SSB-SC)、AM-DSB-WC、AM-DSB-SC、調(diào)頻(frequency modulation,F(xiàn)M)、高斯最小頻移鍵控(Gaussian filtered minimum shift keying,GMSK)。數(shù)據(jù)以浮點復(fù)數(shù)的形式存儲,有200萬個信號樣本,每個樣本的長度為1 024個采樣值。

3.5 評價指標

基于DL的調(diào)制識別性能主要從識別精度和模型復(fù)雜度兩個方面評價。

3.5.1 識別精度

識別精度的第一個衡量指標為混淆矩陣(confusion matrix)。混淆矩陣也稱誤差矩陣,能夠直觀反映識別結(jié)果的分布狀況,通常為M×M的方陣,其第p行、第q列的值M(p,q)表示有M(p,q)個調(diào)制方案為Sp的信號樣本被識別為調(diào)制方案Sp。假設(shè)候選集中有2種調(diào)制方案,即M=2,那么,其混淆矩陣如圖1所示。

圖1 混淆矩陣

其中,TP(truepositive,TP)和TN(true negative,TN)表示預(yù)測結(jié)果正確,F(xiàn)P(false positive,F(xiàn)P)和FN(false negative,F(xiàn)N)表示預(yù)測結(jié)果錯誤。

識別精度的另一個衡量指標為識別準確率(accuracy,Acc),表示識別正確的信號樣本數(shù)量在所有信號樣本中所占的百分比。以M=2為例,Acc可具體計算如下:

3.5.2 模型復(fù)雜度

DNN的模型復(fù)雜度主要包括時間復(fù)雜度和空間復(fù)雜度。時間復(fù)雜度可以通過訓(xùn)練時間或識別時間衡量;空間復(fù)雜度取決于模型的參數(shù)數(shù)量和模型所需的存儲空間大小。

表1 基于DL的調(diào)制識別算法對比

為了進一步展示基于DL的調(diào)制識別算法的性能,表1列舉了本文所提及的部分算法識別精度和模型復(fù)雜度,并給出了它們的數(shù)據(jù)預(yù)處理方式和DNN類型。從表1可以看出,大部分算法均具有較好的識別性能。其中,候選調(diào)制方案越簡單,SNR越高,識別性能越好。

4 未來研究方向

雖然基于DL的調(diào)制識別研究已經(jīng)取得了一系列的成果,但仍有進一步研究的空間。

首先,從數(shù)據(jù)預(yù)處理的角度看,在同一種信號表征形式下,如果候選調(diào)制方案不同,或者無線環(huán)境不同,算法的識別精度會有明顯差異,因此有必要研究信號表征與調(diào)制方案及信道條件之間的關(guān)聯(lián)性,以便根據(jù)實際場景選擇恰當?shù)臄?shù)據(jù)預(yù)處理方式。同時,綜合利用信號的特征表征、序列表征和圖像表征,已被證實有利于提高調(diào)制識別精度,但如何融合這些表征形式還沒有得到深入的研究。

其次,從DNN的角度看,不同的DNN結(jié)構(gòu)擅長處理不同形式的數(shù)據(jù),因此未來需要進一步根據(jù)數(shù)據(jù)預(yù)處理結(jié)果,設(shè)計匹配的網(wǎng)絡(luò)結(jié)構(gòu),以充分學(xué)習信號的典型特征,實現(xiàn)高精度的調(diào)制識別。

最后,從模型復(fù)雜度的角度看,基于DL的調(diào)制識別的精度雖然在不斷提升,但其模型復(fù)雜度也在不斷增加,導(dǎo)致其難以部署到資源有限的移動平臺和邊緣節(jié)點上。因此,如何在保證識別精度不明顯下降的同時減少網(wǎng)絡(luò)參數(shù)、改善模型效率,值得未來繼續(xù)研究。

5 結(jié)束語

本文首先回顧了傳統(tǒng)的調(diào)制識別方法及DL的相關(guān)基礎(chǔ);然后重點討論了基于DL的調(diào)制識別技術(shù),給出了系統(tǒng)架構(gòu),闡述了特征表征、序列表征和圖像表征3種數(shù)據(jù)預(yù)處理方式,著重分析了3種常用的DNN,包括深度前饋網(wǎng)絡(luò)、CNN和RNN,介紹了RadioML數(shù)據(jù)集,列舉了識別精度和模型復(fù)雜度兩方面的性能指標,并對部分算法進行了總結(jié);接著從數(shù)據(jù)預(yù)處理、DNN和模型復(fù)雜度3個角度展望了下一步的研究工作。未來隨著DL的不斷演進,基于DL的調(diào)制識別技術(shù)也將不斷完善,取得更大的發(fā)展。

猜你喜歡
特征信號
抓住特征巧觀察
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
新型冠狀病毒及其流行病學(xué)特征認識
如何表達“特征”
不忠誠的四個特征
當代陜西(2019年10期)2019-06-03 10:12:04
孩子停止長個的信號
抓住特征巧觀察
基于LabVIEW的力加載信號采集與PID控制
一種基于極大似然估計的信號盲抽取算法
主站蜘蛛池模板: 国产高清在线精品一区二区三区 | 亚洲国产日韩一区| 国产成人亚洲欧美激情| 精品久久蜜桃| 92精品国产自产在线观看 | 天堂成人av| 国产91小视频在线观看| 一本大道香蕉中文日本不卡高清二区 | 一边摸一边做爽的视频17国产| 女同久久精品国产99国| 亚洲成人黄色在线观看| 国产真实乱子伦视频播放| 国产成人精品一区二区不卡| 国产成人综合亚洲欧美在| 精品国产美女福到在线不卡f| 亚洲天堂精品在线观看| 亚洲午夜福利在线| 沈阳少妇高潮在线| 青草91视频免费观看| 99色亚洲国产精品11p| 日本a∨在线观看| 精品伊人久久久大香线蕉欧美| 国产精品久久久久久久伊一| 蜜桃视频一区二区三区| 国产精品嫩草影院av| 成人日韩视频| 日日碰狠狠添天天爽| 91最新精品视频发布页| 日韩无码精品人妻| 日韩无码真实干出血视频| 亚洲日韩精品无码专区| 香蕉视频在线观看www| 91精品久久久无码中文字幕vr| 男女男免费视频网站国产| 素人激情视频福利| 亚洲男女在线| 国产精品视频第一专区| 亚洲伊人天堂| 98精品全国免费观看视频| 极品私人尤物在线精品首页| 91网红精品在线观看| 国产国语一级毛片在线视频| 不卡午夜视频| 伊人久久久大香线蕉综合直播| 国产网站免费看| 久久永久精品免费视频| 欧美国产精品拍自| 毛片免费观看视频| 欧美一区日韩一区中文字幕页| 亚洲精品波多野结衣| 婷婷色中文网| 人妻无码中文字幕第一区| 色噜噜狠狠色综合网图区| 福利国产微拍广场一区视频在线 | 午夜精品久久久久久久无码软件 | 国产区精品高清在线观看| 国产成人亚洲综合A∨在线播放| 亚洲欧美日韩成人在线| a毛片基地免费大全| 超清无码一区二区三区| 亚洲无码91视频| 久久综合成人| 久久中文电影| 日韩人妻无码制服丝袜视频| 成年人福利视频| 色国产视频| 国产精品无码制服丝袜| 国产精欧美一区二区三区| 91www在线观看| 五月综合色婷婷| 亚洲精品va| 全免费a级毛片免费看不卡| 精品国产中文一级毛片在线看 | 美女国内精品自产拍在线播放 | 天天综合网在线| 热99re99首页精品亚洲五月天| 成人福利在线视频免费观看| 国产精品亚洲一区二区三区z| 国产69囗曝护士吞精在线视频| 中文字幕调教一区二区视频| 国产欧美在线观看一区| 国产本道久久一区二区三区|