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

OFDM系統的信道與脈沖噪聲的聯合估計方法

2018-04-19 05:41:47呂新榮李有明余明宸
通信學報 2018年3期
關鍵詞:符號方法

呂新榮,李有明,余明宸

?

OFDM系統的信道與脈沖噪聲的聯合估計方法

呂新榮1,2,李有明1,余明宸1

(1. 寧波大學信息科學與工程學院,浙江 寧波 315211;2. 浙江工商職業技術學院智能電子學院,浙江 寧波 315012)

針對OFDM系統中的脈沖噪聲問題,提出一種基于壓縮感知技術的脈沖噪聲抑制方法。該方法將信道脈沖響應和脈沖噪聲聯合視作一個稀疏向量,將發射數據符號視作未知參數,利用稀疏貝葉斯學習理論聯合估計信道、脈沖噪聲和數據符號。與現有脈沖噪聲抑制方法相比,該方法不僅能夠充分利用全部子載波信息,而且不需要信道和脈沖噪聲的先驗統計信息。仿真結果表明,所提方法在信道估計及誤比特率性能上有明顯改善。

正交頻分復用;信道估計;脈沖噪聲;稀疏貝葉斯學習;壓縮感知

1 引言

在無線通信技術的許多應用領域,如車載通信、智能電網、淺海水下通信等,傳輸的數據信號會受到脈沖噪聲的干擾。脈沖噪聲的來源多種多樣,包括汽車的點火裝置[1]、變電所電氣設備的開關[2]、各種海上作業[3]等。正交頻分復用(OFDM, orthogonal frequency division multiplexing)技術已經廣泛應用于現代無線通信標準中。OFDM的每個子信道之間獨立進行數據檢測,由于脈沖噪聲在時域上具有突發性、持續時間短、幅值高且非高斯的特征,脈沖噪聲能量通過OFDM的傅里葉變換分布到所有子信道上,導致傳統基于最大似然估計方法的數據檢測性能急劇下降[4]。

有效的脈沖噪聲抑制方法成為提升OFDM通信系統性能的重要保證。由于時域上脈沖噪聲的幅值通常遠高于背景噪聲,因此,可以通過設置一個閾值來判斷接收到的信號是否含有脈沖噪聲,從而通過Clipping或Blanking方法消除脈沖噪聲的影響[5]。但時域上的脈沖噪聲的統計特性參數難以準確獲得,導致最優的閾值計算非常困難,因此,在實際中往往依靠經驗設置,這使該類方法性能不僅受到限制,而且還會破壞OFDM子載波之間的正交性[6]。近年來,基于壓縮感知技術的脈沖噪聲抑制方法受到重視[1,3,7],這些方法利用接收端OFDM符號中空子載波信息重構出時域上稀疏的脈沖噪聲采樣值,感知矩陣由傅里葉變換矩陣對應空子載波集合的行向量構成。基于結構壓縮感知技術[8]的原理,一些利用脈沖噪聲塊稀疏性或感知矩陣結構性質的脈沖噪聲估計方法被提出[9,10],這些方法在一定程度上提升了估計性能。由于這些方法沒有考慮信道和發射信號,無法利用數據子載波中的信息,導致只能利用空子載波中的信息作為觀測信息。根據壓縮感知技術原理,更多的觀測信息能顯著提升信號重構性能[11]。頻譜效率等原因使OFDM的空子載波數量受到嚴格限制,導致這些僅利用空子載波的脈沖噪聲估計方法難以獲得更多的觀測信息。

為了獲得脈沖噪聲的更多觀測信息,一些利用所有子載波信息的脈沖估計方法被提出。在接收端,數據子載波上的接收信號可以表示為經過多徑信道衰減后的有用信號與脈沖噪聲的疊加。文獻[7]在假設信道信息完全已知情況下,將有用信號視作噪聲均值的估計值,從而利用所有子載波估計脈沖噪聲。然而在脈沖噪聲環境下,獲得準確的信道狀態信息不是一件容易的事[12]。OFDM信道估計也是設計接收機的關鍵工作之一[13],一些聯合考慮信道估計和脈沖噪聲抑制的方法陸續被提出,但是這些方法設置了一些比較嚴格的條件,例如,信道的脈沖響應長度已知且信道在幾個OFDM符號持續時間內保持不變[14]或假設信道脈沖響應的非零元素和脈沖噪聲的非零元素沒有重疊[15]。文獻[16]結合因子圖(factor graph)和近似消息傳遞(AMP, approximate message passing)方法循環迭代估計信道、脈沖噪聲和解碼,但是該方法需要獲得信道和脈沖噪聲的先驗信息及信道編碼的支持,而且收斂性沒有保證。

本文利用信道脈沖響應和脈沖噪聲的聯合稀疏性,提出了一種基于信道估計和數據檢測的脈沖噪聲抑制方法。該方法首先利用接收符號中的所有子載波信息構建壓縮感知方程,然后引入稀疏貝葉斯學習理論[17]估計未知參數。與現有的方法相比,本文的貢獻主要有以下幾點。1) 提出了基于貝葉斯推理的信道估計與脈沖噪聲聯合估計算法,本文算法除了假設信道和脈沖噪聲的稀疏性以外,不需要其他先驗統計信息;2) 針對觀測矩陣中部分元素(即發射符號)未知的問題,通過將發射的數據符號視為未知參數,利用期望值最大化(EM, expectation maximization)算法進行數據符號估計,實現了信道、脈沖噪聲和數據符號聯合估計;3) 本文算法利用數據符號估計構建了基于全部子載波的數據觀測矩陣,突破了傳統基于壓縮感知的脈沖噪聲估計算法僅能利用空子載波的缺點,在提高了脈沖噪聲估計性能的同時也改善了脈沖噪聲環境下的信道估計性能。

仿真結果表明,與沒有聯合考慮信道、脈沖噪聲和數據檢測的方法相比,本文算法能使OFDM系統獲得更好的信道估計性能與誤比特率性能。

2 系統模型

假設OFDM系統總的子載波數為,其中,導頻子載波數為,空子載波數為,其余用于數據子載波。信號源產生的二進制數據流經過信道編碼與調制后,與導頻符號一起映射為OFDM頻域符號(0,1,…,x?1)T。頻域符號經過IDFT運算轉換為OFDM時域信號,時域信號插入循環前綴(CP, cyclic prefix)后送入信道。OFDM的時域接收信號可以表示為

其中,是一個循環矩陣,矩陣第一列是經過對信道脈沖響應向量(0,1,…,h?1)T進行尾部零元素填充后構成,是信道時延擴展長度;表示歸一化后的離散傅里葉變化矩陣(酉矩陣),是的共軛轉置;表示脈沖噪聲的時域采樣值;表示加性高斯白噪聲(AWGN, additive white Gaussian noise)。

去掉循環前綴后,時域信號經過FFT運算轉變為頻域符號,該過程可以表示為

根據OFDM的原理,式(2)可以寫成另外一種形式,即

其中,符號矩陣diag()是對角矩陣,對角元素由頻域符號的元素組成。

在無線通信中,離散信道脈沖響應向量(0,1,…,h?1)T可以用式(4)建模。

脈沖噪聲的出現次數和瞬時幅值都是一個隨機過程,通常使用Bernoulli-Gaussian(BG)、Gaussian-mixture(GM)、Middleton class A(MCA)這3種模型描述脈沖噪聲。由于BG和MCA這2種模型可以由GM模型表示,因此,本文采用GM模型仿真脈沖噪聲。在GM模型中,脈沖噪聲采樣值被認為是若干個均值為0、方差不同的隨機變量的總和。如用復數隨機變量表示,相應的概率密度函數()可以由式(5)表示。

3 同步符號檢測的聯合估計算法

因為向量是由信道脈沖響應向量和脈沖噪聲向量聯合構成,而這兩者均是稀疏向量,因此,可以認為向量也是稀疏向量。注意矩陣的行數少于列數,是一個欠定矩陣。這樣通過式(6)求解未知稀疏向量就轉化成了一個壓縮感知問題,為了解決這個問題,本文引入稀疏貝葉斯學習(SBL, sparse Bayesian learning)理論[17]。

3.1 SBL理論

考慮一個典型的壓縮感知模型[19]

在SBL算法中,假設未知向量的先驗概率分布服從獨立高斯分布,可表示為

其中,均值和協方差矩陣分別為

其中,均值是的最大后驗估計值。算法收斂后,大部分的γ會趨于0,從而μ也會趨于0,得到的稀疏解。在式(10)和式(11)中存在未知參數和,因此,均值和協方差矩陣無法直接計算。為了估計這些超參數,通常使用Type-Ⅱ最大似然估計方法[20],該方法通過最大化接收信號的邊緣概率求解超參數。由于無法直接獲得超參數的閉式解,需要利用EM算法迭代逼近最優解。

在壓縮感知問題中,感知矩陣的列相關性(coherence)對信號的重構性能有決定性影響。與凸優化方法和貪婪迭代方法相比,當矩陣的列相關性很強時,SBL算法仍具有良好的性能[20],因此,SBL算法在大量場合得到了廣泛應用[21]。

3.2 JCIS算法

設未知向量的先驗概率分布服從獨立高斯分布,即

根據貝葉斯理論,給定后的后驗概率分布仍舊為高斯分布,即

其中,均值和協方差矩陣分別為

其中,E{?}代表求期望,(k)表示第次迭代求得的參數值。聯合概率分布(?)求對數后的單調性保持不變,但是運算更加方便。

(k+1)和(k+1)的解分別為

其中,表示所有與無關的量,(?)表示求矩陣的跡。結合,,的定義可得

在實際中,信道與脈沖噪聲可以認為是互相獨立且不相關的,因此可將式(21)中的ΣΣ中所有元素強制置為0,將式(21)代入式(20),利用是對角陣且OFDM每個子信道獨立進行數據檢測的特性以及矩陣跡的性質,可得

獲得脈沖噪聲的估計值后,先去除接收信號中的脈沖噪聲,該過程可以表示為

本文所提JCIS算法需要提供數據符號的初始值,如果隨機選擇星座點會導致算法不收斂。可以先使用傳統的SBL算法利用導頻子載波和空子載波聯合估計出信道和脈沖噪聲,然后經過均衡和判決后得到的估計值(0),將這個估計值以及信道和脈沖噪聲的估計值一同作為JCIS的初始值輸入。

圖1是本文提出的聯合估計信道、脈沖噪聲和數據符號的接收機結構。在初始階段,數據符號是未知的,因此,接收信號利用SBL模塊進行初始的信道與脈沖噪聲估計,然后通過去除脈沖噪聲和信道均衡后得到一個初始的數據符號估計值。利用這個初始的數據符號估計值構建一個初始的數據符號矩陣,利用JCIS模塊進行聯合的信道、脈沖噪聲和數據符號迭代估計,算法收斂后給出最終的數據符號估計值。

3.3 復雜度分析

4 仿真分析

圖1 本文所提接收機結構

表1OFDM系統參數

下面先分析本文提出的JCIS算法在各種不同的導頻子載波數和空子載波數的情況下性能隨信噪比的變化情況,然后對比分析JCIS與其他算法的性能。

4.1 JCIS的性能分析

首先,分析分別采用導頻子載波數目為64、44、24時JCIS的信道估計和誤比特率性能隨信噪比的變化情況,結果如圖2和圖3所示,其中,空子載波數目固定為50。

圖2 不同導頻子載波數目下信道估計均方誤差性能變化

圖3 不同導頻子載波數目下誤比特率變化

從圖2可以看出,隨著導頻數目的減少,JCIS的性能只下降了3 dB左右。導頻數目為64和44這2種情形下,JCIS的性能幾乎一致。在圖3中,隨著導頻數目的減少,JCIS的誤比特率性能均隨之下降。當導頻數減少為44時,JCIS下降了2.5 dB左右;當導頻數減少為24時,JCIS的性能有明顯降低。

接下來,分析空子載波數目分別為20、50、80時JCIS的信道估計和誤比特率性能隨信噪比變化情況,結果如圖4和圖5所示,其中,導頻子載波數目固定為44。

從圖4可以看出,空子載波數目的變化對JCIS的信道估計性能幾乎沒有影響。在圖5中,JCIS的誤比特率性能均隨著空子載波數目的減少而下降。從圖5中可以看出,當空子載波數為50時,JCIS比空子載波為80時下降了1.5 dB左右。當空子載波數為20時,JCIS又下降了2~3 dB。

圖4 不同空子載波數目下信道估計均方誤差性能變化

圖5 不同空子載波數目下誤比特率變化

從圖2~圖5可以看出,導頻子載波數目和空子載波數目的變化對JCIS的信道估計性能影響較小,但對誤比特率性能影響較明顯,原因在于信道估計誤差和脈沖噪聲估計誤差均會影響符號檢測的性能。

4.2 JCIS與其他算法的性能對比

下面,將本文提出算法與其他算法進行對比。圖6和圖7分別對比了導頻子載波數為44時,改變空子載波數后信道估計與誤比特率性能隨信噪比的變化情況。“LS”方法[13]是假設完全去除脈沖噪聲且信道抽頭系數位置已知,利用最小二乘法估計信道,這是一種最大似然準則下最優的信道估計算法;“SBL-LS”方法[7]先基于空子載波利用SBL方法去除脈沖噪聲后,假設殘留噪聲為高斯白噪聲,然后利用“LS”方法估計信道參數,這種方法信道估計和脈沖噪聲估計獨立進行;“JCSwIN”方法先基于空子載波利用SBL方法去除脈沖噪聲后,利用文獻[11]提出的方法獨立進行信道估計,該方法利用所有子載波提高信道估計性能。

圖6 信道估計性能對比(導頻子載波數為44)

由于“LS”方法在導頻子載波數量固定的情況下是最優估計,因此,將其作為對比的基準。從圖6可以看出,當空子載波數為100時,“JCIS”和“SBL-LS”性能幾乎相同,與“LS”相差2.5 dB;當空子載波數減少為50后,“JCIS”性能幾乎不變,而“SBL-LS”和“JCSwIN”均顯著下降,與“LS”相差約10 dB。在圖7中,當空子載波數為100時,在誤比特率性能10?3處,“JCIS”“SBL-LS”“JCSwIN”與“LS”分別相差2 dB、4 dB和5 dB;當空子載波數減少為50后,3種方法性能都有下降,但“JCIS”下降明顯更少,在誤比特率性能10?3處,“JCIS”“SBL-LS”“JCSwIN”與“LS”分別相差5 dB、9 dB和9 dB。

圖7 誤比特率性能對比(導頻子載波數為44)

從對比結果可以看出,當空子載波數為50時,“JCIS”性能顯著好于“SBL-LS”和“JCSwIN”方法,而當空子載波數為100時,由于空子載波數增加引起脈沖噪聲估計性能提升,但“JCIS”性能仍優于“SBL-LS”和“JCSwIN”方法。這說明本文提出的信道和脈沖噪聲聯合估計的方法要顯著優于傳統信道和脈沖噪聲獨立估計的方法,而傳統獨立估計方法的性能要顯著依賴于空子載波的數量。當空載波數為50時,“JCSwIN”方法通過結合符號檢測可以彌補空載波數量不足帶來的影響,使在高信噪比場合性能略好于“SBL-LS”。

5 結束語

本文提出了一種基于稀疏貝葉斯學習理論的信道、脈沖噪聲和數據符號的聯合估計方法。首先,利用信道和脈沖噪聲的聯合稀疏性,構建了一個基于全部子載波信息的壓縮感知方程。為了解決觀測矩陣部分元素未知的困難,將發射符號視作未知參數,利用期望值最大化算法迭代估計信道、脈沖噪聲和數據符號。仿真結果表明,本文提出的方法通過利用所有子載波上的信號作為觀測信息,能較好地提升脈沖噪聲環境下的OFDM系統性能。

[1] LIU S, YANG F, DING W, et al. Double kill: compressive-sensing-based narrow-band interference and impulsive noise mitigation for vehicular communications[J]. IEEE Transactions on Vehicular Technology, 2016, 65(7):5099-5109.

[2] SACUTO F, LABEAU F, AGBA B L. Wide band time-correlated model for wireless communications under impulsive noise within power substation[J]. IEEE Transactions on Wireless Communications, 2014, 13(3):1449-1461.

[3] KUAI X, SUN H, ZHOU S, et al. Impulsive noise mitigation in underwater acoustic OFDM systems[J]. IEEE Transactions on Vehicular Technology, 2016, 65(10): 8190-8202.

[4] MAHMOOD A, CHITRE M, ARMAND M A. Detecting OFDM signals in alpha-stable noise[J]. IEEE Transactions on Communications, 2014, 62(10):3571-3583.

[5] JUWONO F H, GUO Q, HUANG D, et al. Deep clipping for impulsive noise mitigation in OFDM-based power-line communications[J]. IEEE Transactions on Power Delivery, 2014, 29(3):1335-1343.

[6] DARSENA D, GELLI G, MELITO F, et al. ICI-free equalization in OFDM systems with blanking preprocessing at the receiver for impulsive noise mitigation[J]. IEEE Signal Processing Letters, 2015, 22(9): 1321-1325.

[7] LIN J, NASSAR M, EVANS B L. Impulsive noise mitigation in powerline communications using sparse Bayesian learning[J]. IEEE Journal on Selected Areas in Communications, 2013, 31(7):1172-1183.

[8] DUARTE M F, ELDAR Y C. Structured compressed sensing: from theory to applications[J]. IEEE Transactions on Signal Processing, 2011, 59(9):4053-4085.

[9] AL-NAFFOURI T Y, QUADEER A A, CAIRE G. Impulse noise estimation and removal for OFDM systems[J]. IEEE Transactions on Communications, 2014, 62(3):976-989.

[10] KORKI M, ZHANG J, ZHANG C, et al. Block-sparse impulsive noise reduction in OFDM systems-a novel iterative bayesian approach[J]. IEEE Transactions on Communications, 2015, 64(1): 271-284.

[11] PRASAD R, MURTHY C R, RAO B D. Joint approximately sparse channel estimation and data detection in OFDM systems using sparse Bayesian learning[J]. IEEE Transactions on Signal Processing, 2014, 62(14): 3591-3603.

[12] ZHENG Z, HAO C, YANG X. Least squares channel estimation with noise suppression for OFDM systems[J]. Electronics Letters, 2016, 52(1): 37-39.

[13] LIU Y, TAN Z, HU H, et al. Channel estimation for OFDM[J]. IEEE Communications Surveys & Tutorials, 2014, 16(4):1891-1908.

[14] CHIEN Y R. Iterative channel estimation and impulsive noise mitigation algorithm for OFDM-based receivers with application to power-line communications[J]. IEEE Transactions on Power Delivery, 2015, 30(6): 2435-2442.

[15] MEHBOOB A, LI Z, KHANGOSSTAR J, et al. Joint channel and impulsive noise estimation for OFDM based power line communication systems using compressed sensing[C]//IEEE International Symposium on Power Line Communications and ITS Applications. 2013: 203-208.

[16] NASSAR M, SCHNITER P, EVANS B L. A factor graph approach to joint OFDM channel estimation and decoding in impulsive noise environments[J]. IEEE Transactions on Signal Processing, 2014, 62(6): 1576-1589.

[17] WIPF D P, RAO B D. Sparse Bayesian learning for basis selection[J]. IEEE Transactions on Signal Processing, 2004, 52(8): 2153-2164.

[18] DING W, YANG F, PAN C, et al. Compressive sensing based channel estimation for OFDM systems under long delay channels[J]. IEEE Transactions on Broadcasting, 2014, 60(2): 313-321.

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

[20] WIPF D D. Sparse estimation with structured dictionaries[C]//International Conference on Neural Information Processing Systems. 2011: 2016-2024.

[21] DING W, LU Y, YANG F, et al. Spectrally efficient CSI acquisition for power line communications: a Bayesian compressive sensing perspective[J]. IEEE Journal on Selected Areas in Communications, 2016,34 (7): 2022-2032.

[22] ITU-R M.1225. Test environment and deployment models-ANNEX2[S]. ITU-R, 2000: 24-30.

Joint channel and impulsive noise estimation method for OFDM systems

LYU Xinrong1,2, LI Youming1, YU Mingchen1

1. Faculty of Information Science and Engineering, Ningbo University, Ningbo 315211, China2. School of Intelligent Electronics, Zhejiang Business & Technology Institute, Ningbo 315012, China

Aiming at the impulsive noise occurring in OFDM systems, an impulsive noise mitigation algorithm based on compressed sensing theory was proposed. The proposed algorithm firstly treated the channel impulse response and the impulsive noise as a joint sparse vector by exploiting the sparsity of both them. Then the sparse Bayesian learning framework was adopted to jointly estimate the channel impulse response, the impulsive noise and the data symbols, in which the data symbols were regarded as unknown parameters. Compared with the existing impulsive noise mitigation methods, the proposed algorithm not only utilized all subcarriers but also did not use any a priori information of the channel and impulsive noise. The simulation results show that the proposed algorithm achieves significant improvement on the channel estimation and bit error rate performance.

orthogonal frequency division multiplexing, channel estimation, impulsive noise, sparse Bayesian learning, compressed sensing

TN911.4

A

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

2017-04-20;

2018-02-06

李有明,liyouming@nbu.edu.cn

國家自然科學基金資助項目(No.61571250);寧波市自然科學基金資助項目(No.2015A610121)

The National Natural Science Foundation of China (No.61571250), The Natural Science Foundation of Ningbo (No.2015A610121)

呂新榮(1976-),男,浙江永康人,寧波大學博士生,主要研究方向為無線通信技術、電力線通信、稀疏信號處理。

李有明(1963-),男,陜西扶風人,博士,寧波大學教授,主要研究方向為無線寬帶通信、電力線通信、協作中繼、認知無線電等。

余明宸(1991-),男,河南洛陽人,寧波大學碩士生,主要研究方向為電力線通信技術。

猜你喜歡
符號方法
學符號,比多少
幼兒園(2021年6期)2021-07-28 07:42:14
“+”“-”符號的由來
學習方法
變符號
用對方法才能瘦
Coco薇(2016年2期)2016-03-22 02:42:52
四大方法 教你不再“坐以待病”!
Coco薇(2015年1期)2015-08-13 02:47:34
賺錢方法
捕魚
倍圖的全符號點控制數
圖的有效符號邊控制數
主站蜘蛛池模板: 色精品视频| 日韩免费成人| 亚洲av无码成人专区| 色屁屁一区二区三区视频国产| 凹凸精品免费精品视频| 99成人在线观看| 91亚洲精品第一| 欧美亚洲国产一区| 亚洲成肉网| 精品视频一区二区三区在线播| 日韩中文无码av超清| 2022精品国偷自产免费观看| 亚洲中文字幕精品| 亚洲综合18p| 亚洲欧美日韩成人高清在线一区| 女高中生自慰污污网站| 国产午夜精品一区二区三区软件| 日a本亚洲中文在线观看| 2021国产v亚洲v天堂无码| 亚洲 日韩 激情 无码 中出| 小13箩利洗澡无码视频免费网站| 中文字幕在线观| 国产高清免费午夜在线视频| 五月天久久综合| 国产区免费精品视频| 亚洲一区二区视频在线观看| 国产成年无码AⅤ片在线| 中文字幕伦视频| 成人午夜在线播放| 国产高清无码第一十页在线观看| 国产第一页免费浮力影院| 少妇露出福利视频| 亚洲成人www| 中文字幕 欧美日韩| 国产成人资源| a天堂视频在线| 99re在线免费视频| 日本一本在线视频| 久久 午夜福利 张柏芝| 欧美激情视频一区| 91无码人妻精品一区二区蜜桃| 97在线国产视频| 99精品视频播放| 国产免费一级精品视频| 精品少妇人妻一区二区| 欧美自慰一级看片免费| 三区在线视频| 亚洲欧美另类久久久精品播放的| 国产精品林美惠子在线播放| 91久草视频| 久久精品无码国产一区二区三区| 亚洲AⅤ波多系列中文字幕 | 亚洲精品午夜天堂网页| 久久久久久尹人网香蕉| 秋霞一区二区三区| 日韩国产精品无码一区二区三区| 深夜福利视频一区二区| 国禁国产you女视频网站| 午夜视频免费一区二区在线看| 中文字幕2区| 最新国产成人剧情在线播放| 呦女精品网站| julia中文字幕久久亚洲| 91麻豆久久久| 亚洲欧美不卡中文字幕| 国产69囗曝护士吞精在线视频 | 欧美日韩国产在线播放| 国产色网站| 老司机精品99在线播放| 青青极品在线| 九九热视频精品在线| 91av成人日本不卡三区| 国产免费久久精品99re不卡| 日韩精品专区免费无码aⅴ| 欧美精品一区二区三区中文字幕| 青青草原国产一区二区| 玖玖免费视频在线观看 | 亚洲精品你懂的| 久久久精品国产亚洲AV日韩| 亚洲伊人天堂| 亚洲欧洲天堂色AV| 国产成人毛片|