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

    Polar碼編譯碼技術(shù)研究

    2021-11-24 07:39:40魯信金舒冰心
    無線電通信技術(shù) 2021年6期

    魯信金,舒冰心,雷 菁

    (國防科技大學(xué) 電子科學(xué)學(xué)院,湖南 長沙 410000)

    0 引言

    通信系統(tǒng)中信道干擾和噪聲總是影響系統(tǒng)的可靠性,信道編碼技術(shù)可用最小代價來提升通信系統(tǒng)可靠性。1948年,香農(nóng)(Claude Shannon)在《通信的數(shù)學(xué)理論》[1]中指出對于去認(rèn)定的信道容量,當(dāng)碼率小于其信道容量時,總是可以找到一種編碼方案使得碼長足夠大時,使用最大似然譯碼可使誤碼率任意小,這奠定了糾錯編碼理論的基礎(chǔ)。此后,很多學(xué)者致力于研究設(shè)計出能夠逼近信道容量(也即香農(nóng)界)的實用編碼方案,如線性分組碼[2]、卷積碼[3]和代數(shù)幾何碼[4]。進(jìn)入20世紀(jì)90年代,信道編碼進(jìn)入了第一次快速發(fā)展的階段,研究的重點集中為Turbo碼[5]和LDPC碼[6],隨后二者在實際的通信系統(tǒng)得到了廣泛的應(yīng)用和很深入的研究,雖然這兩種編碼沒有達(dá)到理論證明的香農(nóng)限,還具有很高的復(fù)雜度,但依然備受青睞。

    近年來,無線通信物理層中的編碼領(lǐng)域發(fā)展較慢。在2009年,E.Arikan等人研究了信道極化問題,根據(jù)信道極化提出了極化碼[7]。不僅證明了極化碼在二進(jìn)制離散無記憶信道(Binary Discrete Memoryless Channel,BDMC)下,可以達(dá)到信道容量,還提出了具有低復(fù)雜度的編碼和譯碼方案。R.Moil和T.Tanaka[8]針對任意對稱二進(jìn)制無記憶信道提出了一個新的被證明可實現(xiàn)的編碼構(gòu)造方法。編碼構(gòu)造方面對于不同的信道,Bhattacharyya參數(shù)的計算和信息位選擇方法是不同的,而很多學(xué)者研究了很多適用于不同信道的信息位選擇的方法,如蒙特卡洛迫近法、密度演化構(gòu)造方法等。

    對于信道編碼而言,譯碼算法的復(fù)雜度和性能至關(guān)重要。自極化碼提出后就受到大量學(xué)者對其譯碼算法進(jìn)行研究。目前polar碼譯碼算法主要有SC算法[9]、SCL算法[10]、CA-SCL算法[11]、BP算法[12]、SCAN算法[13],以及各種算法的簡化版本[14-16]。其中SC算法最初由Arikan提出,但其在碼長有限長情形下,性能一般。SCL是SC算法性能提升的改進(jìn)版本,原理是在SC的基礎(chǔ)上提供了多條路徑;而CA-SCL是在SCL的基礎(chǔ)上對信息比特進(jìn)行了循環(huán)冗余校驗,它通過簡單的校驗就可以帶來性能的極大提升。目前基于CA-SCL譯碼算法,polar碼性能已經(jīng)優(yōu)于LDPC碼。SC、SCL以及CA-SCL算法均是硬輸出算法,即最終輸出的是0、1比特序列,而不是對應(yīng)的LLR值。為了在一些聯(lián)合設(shè)計中使用polar碼,就必須提供polar碼的軟輸出算法,即輸出對應(yīng)比特的LLR值。BP和SCAN算法就是軟輸出算法,由于BP和SCAN算法在消息傳遞和遞歸規(guī)則上不同,導(dǎo)致了兩種算法的譯碼延時和收斂速度不同。譯碼延時上,BP消息傳遞采用“洪水”規(guī)則[17],SCAN采用類SC算法的串行消除規(guī)則,因此BP算法的譯碼延時好于SCAN算法。SCAN算法1次迭代就能達(dá)到稍好于SC算法的性能,一般而言,BP算法需要40~50次的迭代過程,而SCAN算法1次迭代就能達(dá)到稍好于SC算法的性能。性能上幾種算法排序為:CA-SCL>LDPC>SCL>BP=SCAN>SC。

    從極化碼被廣泛關(guān)注研究之后,極化碼逐漸應(yīng)用在現(xiàn)代通信系統(tǒng)中以及不同信道場景中,如二進(jìn)制對稱信道、混合信道、非對稱信道、非二進(jìn)制信道、多址接入信道、并行信道、高斯信道以及瑞利信道等。

    1 信道極化原理

    polar碼的基本原理是信道極化。信道極化是這樣的一種現(xiàn)象:對二進(jìn)制對稱信道進(jìn)行特定的“組合”和“拆分”,則拆分后的“比特信道”將呈現(xiàn)極化現(xiàn)象,即一部分“比特信道”的對稱信道容量趨近于1,而其余部分“比特信道”的對稱信道信道容量趨近于0。如果發(fā)送端將源信息比特放置在“好的比特信道上”,而在“壞的比特信道上”放置固定比特(如0),同時在接收端采用連續(xù)消除譯碼算法,在碼長N→∞ 時將可以達(dá)到信道容量。雖然基于二進(jìn)制離散無記憶信道提出了信道極化,但是信道極化的現(xiàn)象是普遍存在的,除B-DMC信道外其他信道也存在此現(xiàn)象,而信道極化包括信道組合和信道分解兩個重要過程,接下來對其進(jìn)行介紹。

    1.1 信道組合

    在物理信道W(y|x)上連續(xù)發(fā)送N個符號,若發(fā)送端對這N個符號進(jìn)行特定方式的組合,則等效于信道W的組合。

    如圖1所示,對發(fā)送的信息源符號以一定方式進(jìn)行組合,并記組合后的信道為W2(y1,y2|u1,u2),則:

    W2(y1,y2|u1,u2)=W(y1|u1⊕u2)W(y2|u2)。

    (1)

    進(jìn)一步將兩個W2信道組合為W4信道,依次往復(fù),直至遞歸到兩個WN/2信道組合為WN信道,如圖2所示。其中,RN是置換矩陣,其滿足RN:(v1,v2,…vN)→(v1,v3,…vN-1,v2,v4,…vN)。

    圖1 W2信道Fig.1 Channel W2

    圖2 WN信道Fig.2 Channel WN

    1.2 信道拆分

    (2)

    可以將上式迭代的形式,保證極化碼的編碼和譯碼具有線性復(fù)雜度。圖3展示了N=8時,信道W8分解的過程:獨立的8個信道W通過信道組合之后可以得到一個信道W8,信道分解時,信道W8先分解為兩個相互獨立的信道W4,之后信道W4分解為兩個相互獨立的信道W2,最后信道W2分解為兩個相互獨立的信道W,因此信道分解是一個遞歸的過程。

    圖3 N=8信道分解Fig.3 Channel decomposition(N=8)

    1.3 信道極化定理

    polar碼選擇在“比特信道”量為1的比特信道上發(fā)送信源信息(比特信道容量為0的比特信道上發(fā)送固定信息,如0)。當(dāng)信道W為二進(jìn)制刪除(Binary Erasure Channel,BEC)信道時,圖4是碼長N=210刪余概率ε=0.5時仿真結(jié)果圖,圖4(a)是信道極化現(xiàn)象的仿真,圖4(b) 是各信道容量占比統(tǒng)計柱狀圖。圖5是N=220,刪余概率ε=0.5時的仿真結(jié)果圖,圖6是N=220,刪余概率ε=0.3時的仿真結(jié)果圖。

    (a) 信道極化

    (b) 各信道容量占比統(tǒng)計柱狀圖圖4 N=210,ε=0.5時信道極化現(xiàn)象及各信道容量占比Fig.4 Channel polarization phenomenon and the proportion of each channel capacity when N=210,ε=0.5

    (a) 極化現(xiàn)象

    (b) 各信道容量占比圖5 N=220,ε=0.5時信道極化現(xiàn)象及各信道容量占比Fig.5 Channel polarization phenomenon and the proportion of each channel capacity when N=220,ε=0.5

    (a) 極化現(xiàn)象

    (b) 各信道容量占比圖6 N=220,ε=0.3時信道極化現(xiàn)象及各信道容量占比Fig.6 Channel polarization phenomenon and the proportion of each channel capacity when N=220,ε=0.3

    圖4和圖5可以看出,當(dāng)ε=0.5時,隨著N的增大,信道容量為0和1的信道所占比例增大,信道極化的現(xiàn)象越來越明顯。由圖5和圖6對比可以看出,當(dāng)N=220,信道容量為0時的信道所占比例趨向于信道刪除概率ε的值,信道容量為1的信道所占比例趨于信道刪除概率1-ε,與上述定理所述一致。

    基于信道極化的現(xiàn)象,選擇其中較好的信道,即信道容量趨向于I(W)的信道作為信息位傳輸數(shù)據(jù),而較差的信道,即信道容量趨向于1-I(W)的信道作為凍結(jié)位。運(yùn)用該思想進(jìn)行極化碼的編碼和譯碼。

    2 polar碼編碼構(gòu)造

    polar碼碼長N→∞時,polar碼可達(dá)信道容量。而實際編碼過程中碼長N不可能趨于無窮,因此polar碼也不可能做到信道完全極化??傆胁糠帧氨忍匦诺馈钡男诺廊萘拷橛?和1中間。因此確定對稱信道容量最大的K個“比特信道”是解決polar碼碼字構(gòu)造的主要問題;其次,生成矩陣決定了信道重組的形式,信道重組的形式作用于待發(fā)送信息上的函數(shù)運(yùn)算法則,經(jīng)典的極化編碼方案以二維核矩陣F=[1 0;1 1]為基底,迭代構(gòu)造生成矩陣。

    2.1 生成矩陣

    (3)

    類似地,對于重組信道W4,其生成矩陣G4可以簡化為G2的矩陣形式:

    (4)

    其中,?為Kronecker內(nèi)積。繼而推算出任意生成矩陣的遞歸表達(dá)式,即生成矩陣GN:

    GN=F?log2(N),

    (5)

    其中,F=G2??紤]到連續(xù)消除算法的譯碼順序,Arikan在生成矩陣中引入比特翻轉(zhuǎn)置換操作:對于任意正整數(shù)i,定義(b1,b2...,bn)為其二進(jìn)制表示,比特翻轉(zhuǎn)操作rvsl(i)=(bn...,b2,b1),對于任意矢量(v0,v1...,vn-1),其比特翻轉(zhuǎn)置換后的矢量為(vrvsl(0),vrvsl(1)...,vrvsl(n-1)),該操作可以用N×N維的置換矩陣BN表示。設(shè)BN中第i行第j列的元素為bi,j,則有:

    (6)

    則生成矩陣的一般形式為:

    GN=BNF?log2(N)。

    (7)

    2.2 信息集

    信息集是信息信道的索引集合,而信息信道的質(zhì)量通常由對稱容量及巴特查里亞因子所度量。對稱容量與巴特查里亞因子是信道性質(zhì)的相反描述度量,其中對稱容量的計算復(fù)雜度更高,因此巴特查里亞因子更多地應(yīng)用在實際系統(tǒng)中。在此,僅討論基于巴特查里亞因子的信息集選取準(zhǔn)則。

    (8)

    (9)

    (10)

    當(dāng)且僅當(dāng)物理信道為BEC時,式(8)及式(9)的等號成立,此時比特信道的巴特查里亞因子可迭代計算。但對于其他二進(jìn)制信道而言,比特信道的轉(zhuǎn)移概率計算復(fù)雜度較大,其巴特查里亞因子難以確定,當(dāng)前較為常用的處理方法為近似迭代法。近似迭代法針對不同物理信道重定義了巴特查里亞因子,并使用近似迭代準(zhǔn)則迭代計算比特信道的巴特查里亞因子。在此給出BSC、BAWGNC、BFC的巴特查理亞因子。

    對于交叉概率為Pε的BSC,其巴特查理亞因子按原定義計算:

    (11)

    對于連續(xù)信道,其巴特查里亞因子的定義為積分函數(shù):

    (12)

    對于二進(jìn)制加性高斯白噪聲(Binary AWGN,BAWGN)信道,根據(jù)定義計算其轉(zhuǎn)移概率,其中σ2為噪聲方差,調(diào)制方式為BPSK(Binary Phase Shift Keying)。

    (13)

    (14)

    將式(13)與式(14)帶入式(12),可得:

    (15)

    對于二進(jìn)制衰落信道(Binary Fading Channel),可將其近似為AWGN信道處理,則其轉(zhuǎn)移概率為:

    (16)

    (17)

    (18)

    2.3 polar碼編碼

    polar碼實際上就是線性分組碼,因此其編碼過程和一般的線性分組碼一樣。不同的是一般的線性分組碼根據(jù)碼字漢明重量等挑選生成矩陣某些行(RM碼根據(jù)的是最大漢明重量),而polar碼由信息比特集A確定生成矩陣的某些行。假設(shè)C為N比特長的碼字,UI為K比特信源信息,UF=0為N-K比特固定信息,則:

    C=UGN×N=UIGK×N+UFG(N-K)×N,

    (19)

    其中,GN×N是polar碼生成矩陣,GK×N是根據(jù)信息比特集A中元素從GN×N中挑選的行組成矩陣。

    3 polar碼譯碼

    3.1 SC譯碼原理

    (20)

    判定準(zhǔn)則為:

    (21)

    (22)

    如圖7所示,首先譯碼器計算標(biāo)號為1的點的左消息LLR值,而此時必須知道標(biāo)號為2和3點的LLR值,同理要計算2的LLR值應(yīng)計算4和5的LLR值,而3需要計算6和7的LLR值,而4、5、6、7的LLR值可以直接根據(jù)信道傳輸過來的LLR值直接進(jìn)行計算。通過以上遞歸過程則得到了1點的LLR值。此時先判定該點是否是固定比特,如果是固定比特則直接譯為0;如果是信息比特,則根據(jù)剛計算得到的LLR值進(jìn)行判決:LLR值大于0判為0、LLR值小于0判為1。同理直至SC譯碼器將u0~u7的LLR值全部計算出來,并判決得到比特值。

    圖7 遞歸順序Fig.7 Recursive order

    3.2 SCL譯碼算法

    通過這種多路徑譯碼策略,極化碼在中短碼長下的誤碼性能得到顯著提高。當(dāng)L足夠大時,SCL算法的誤碼性能能夠接近于最大似然(ML)算法。但是隨著L的增加,SCL算法的復(fù)雜度也呈指數(shù)倍增長。其時間復(fù)雜度和空間復(fù)雜度均為O(LNlog2N)。所以,SCL算法僅適用于中短碼長的情況,當(dāng)碼長N和譯碼路徑數(shù)L足夠大時,高復(fù)雜度的SCL算法難以實現(xiàn)。

    圖8 SCL譯碼Fig.8 SCL decoding

    3.3 CA-SCL譯碼算法

    CA-SCL算法相較于SCL算法而言,能夠大幅降低譯碼復(fù)雜度。其需要對所有得到的譯碼結(jié)果進(jìn)行校驗,如果滿足校驗,則從中選取轉(zhuǎn)移概率最大的譯碼路徑作為最終譯碼結(jié)果;CA-SCL算法本具有較低的時延,但由于引入了再編碼步驟,導(dǎo)致仍具有較大時延。然而,在時間復(fù)雜度以及空間復(fù)雜度上,CA-SCL算法相比于SCL算法都有較大改進(jìn)。

    4 仿真與分析

    此處采用經(jīng)典SC譯碼,選擇BEC、BSC、AWGN三種信道進(jìn)行仿真,仿真通過計算誤碼率來分析編譯碼方案的優(yōu)劣。

    4.1 BEC信道下極化碼的仿真

    4.1.1 不同刪除概率及不同碼率下

    BEC信道刪除概率p表示信道中傳輸符號在接收端無法識別的概率大小,在極化碼中改變碼率R通過改變極化信道的選擇比例。圖9為碼長N=64時,誤碼率在刪除概率和碼率這兩個參數(shù)變化時的變化曲線。橫坐標(biāo)為刪除概率,不同顏色曲線對應(yīng)不同碼率R。

    圖9 當(dāng)N=64時BER變化曲線Fig.9 BER change curve when N=64

    由圖9可知,每一條曲線變化趨勢一致,隨著BEC信道刪除概率增大,誤碼率變大。這是因為在碼長N和碼率R固定時,隨著BEC信道刪除概率的增大,巴氏參數(shù)Z增大,信道極化現(xiàn)象變差,選擇用來傳輸信息的子信道參數(shù)Z也增大,所以誤碼率增大。

    比較不同顏色的4條曲線,碼率R越大,誤碼率越大。這是因為當(dāng)參數(shù)N及p一定時,極化信道巴氏參數(shù)Z的值確定。當(dāng)碼率R增大時,所選的子信道個數(shù)增大,此時會將參數(shù)Z較大的信道選擇為傳輸信息的信道,將不好的信道選做信息位,造成誤碼的增加。

    4.1.2 不同碼長條件下

    BEC信道中,碼長N值越大,極化現(xiàn)象越明顯,此時碼率R=0.5、刪除概率p=0.2,如圖10所示,隨著碼長N的增大,誤碼率明顯變小。這是因為當(dāng)N值增大時,信道容量為0和1的信道總數(shù)增大,信道容量為其他值所占的總數(shù)減小,即極化的效果更加明顯,此時巴氏參數(shù)Z的值向著0和1兩個方向集中。所以當(dāng)碼長N增大時,選擇用來傳輸信息信道的巴氏參數(shù)值更趨向于0,所以誤碼率減小。

    圖10 當(dāng)p=0.2,R=0.5時BER變化曲線Fig.10 BER change curve when p=0.2,R=0.5

    4.2 BSC信道下極化碼的仿真

    4.2.1 不同交叉概率下

    BSC信道中信元符號按交叉概率的大小比例進(jìn)行翻轉(zhuǎn)。圖11為誤碼率在BSC信道交叉概率p變化時的變化曲線,此時碼長N=256和碼率R=0.1。橫坐標(biāo)為BSC信道的交叉概率,縱坐標(biāo)為誤碼率。

    由圖11可知,隨著BSC信道交叉概率p的增大,誤碼率變大。很明顯,此次仿真的誤碼率與BEC信道具有相似的特性,故造成這種結(jié)果的原因也與其相似,可以總結(jié)為:在碼長N和碼率R固定時,隨著BSC信道交叉概率在0~0.5的范圍內(nèi)增大,巴氏參數(shù)Z增大,信道極化現(xiàn)象變差,選擇用來傳輸信息的子信道參數(shù)Z也增大,所以誤碼率增大。

    圖11 當(dāng)N=256,R=0.1時BER變化曲線Fig.11 BER change curve when N=256,R=0.1

    4.2.2 不同碼率下

    取N=256,p=0.2,圖12為極化碼在不同的碼率條件下的誤碼率曲線。隨著碼率R的增大,誤碼率明顯變大。曲線的變化趨勢與BEC信道的結(jié)果相似,但從仿真結(jié)果可以看出,當(dāng)誤碼率小于0.15時,誤碼率才小于1/100。所以BSC信道中可用的碼率范圍很小,一個碼字可以傳輸?shù)男畔⒂邢?。BSC信道中不同碼長N的情況下,誤碼率與BEC信道具有相似性,隨著碼長N的增大,誤碼率明顯變小。

    圖12 當(dāng)N=256,p=0.2時BER變化曲線Fig.12 BER change curve when N=256,p=0.2

    4.3 AWGN信道下極化碼的仿真

    雖然極化碼是基于B-DMC信道提出的,但信道極化現(xiàn)象卻是普遍存在的,基于前面提到的理論方法對高斯信道中的極化碼進(jìn)行仿真分析,采用BPSK調(diào)制,調(diào)制之后變?yōu)殡p極性的信號然后再傳送至信道,選用加性高斯白噪聲。在高斯信道中,需考慮參數(shù)信噪比SNR,其為高斯信道重要的衡量指標(biāo)。

    4.3.1 不同信噪比下

    圖13為誤碼率在信噪比變化時的變化曲線,此時碼長N=1 024,由圖可知,隨著信噪比的增大,誤碼率變小。

    在高斯信道中,信噪比越大則信道條件越好,出錯的概率越小,類似于BEC中的刪除概率。信噪比越大說明原始信道的條件越好,而高斯信道的極化依賴于信噪比,所以信噪比越大,極化效果件也越好,即巴氏參數(shù)Z的值越小。所以當(dāng)碼長N和碼率R一定時,選擇用來傳輸信息的子信道個數(shù)也就確定了,此時Z的大小就決定了選擇信道的質(zhì)量好壞。因此信噪比越大,信道的方差越小,極化信道參數(shù)Z越小,極化碼性能越好,誤碼率越小。

    4.3.2 不同碼率條件下

    極化碼是一種依賴于信道的信道編碼,所以信道確定時,信道極化的參數(shù)便隨之確定,此時碼率便是一個影響誤碼性能的變量。圖14為誤碼率在碼率變化時的變化曲線,此時碼長N=1 024。

    圖14 當(dāng)N=1 024時BER變化曲線Fig.14 BER change curve when N=1 024

    由圖14可知,隨著碼率R的減小,誤碼率也隨之變小。由于信道極化的參數(shù)已經(jīng)確定,選擇其中的信道都是好的信道。但當(dāng)R增大時,選擇傳輸信息的信道也會增多,這時不太好的信道會加入進(jìn)來,造成誤碼率的增大。

    4.3.3 不同碼長條件下

    碼長的值對信道極化現(xiàn)象有很重要的影響,從而影響誤碼性能。圖15為誤碼率在碼長變化時的變化曲線,此時碼率R=0.3,由圖可知,隨著碼長N的增大,誤碼率明顯變小。這與BEC信道中的情況類似,當(dāng)N值增大時,極化的效果更加明顯,選擇用來傳輸信息的信道巴氏參數(shù)Z的值,肯定比N值小時選出的信道Z值小,所以誤碼率減小。此外,不同的碼長N下,誤碼率隨信噪比的增大差距逐漸變大。

    圖15 當(dāng)R=0.3時BER變化曲線Fig.15 BER change curve when R=0.3

    5 結(jié)論

    極化碼是通信領(lǐng)域?qū)碚撎剿鞯囊淮钨|(zhì)的飛躍,是編碼領(lǐng)域的里程碑,具有很高的研究價值。本文重點介紹了信道極化的原理和信道極化的現(xiàn)象,對信道極化中信道組合和信道分解兩部分進(jìn)行了詳細(xì)分析,并且進(jìn)一步研究了極化碼的編譯碼過程。雖然信道極化是基于二進(jìn)制對稱信道提出的,但是信道極化現(xiàn)象卻是普遍存在的,為此討論了不同信道中極化的計算方法以及編碼參數(shù)對誤碼率的影響。作為目前已知的唯一一種能夠被嚴(yán)格證明達(dá)到香農(nóng)限的信道編碼方法,極化碼可以應(yīng)用到許多領(lǐng)域,后續(xù)的研究還可以從很多方面進(jìn)行改進(jìn)。

    午夜福利,免费看| 国产激情久久老熟女| 黄色怎么调成土黄色| 9热在线视频观看99| av天堂在线播放| 黄色视频不卡| 午夜老司机福利片| 成人精品一区二区免费| 亚洲 国产 在线| www.999成人在线观看| 母亲3免费完整高清在线观看| 色婷婷av一区二区三区视频| 亚洲欧美一区二区三区久久| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人爽人人夜夜| 午夜精品久久久久久毛片777| 久久国产精品人妻蜜桃| 侵犯人妻中文字幕一二三四区| 啦啦啦免费观看视频1| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久亚洲精品国产蜜桃av| 欧美日韩国产mv在线观看视频| 狠狠精品人妻久久久久久综合| 深夜精品福利| 满18在线观看网站| 亚洲欧洲日产国产| 亚洲色图av天堂| 美女午夜性视频免费| 免费观看人在逋| av片东京热男人的天堂| 精品久久久精品久久久| 精品一区二区三卡| 一级片免费观看大全| 老司机福利观看| 超碰成人久久| 欧美黑人精品巨大| 免费人妻精品一区二区三区视频| 黄色a级毛片大全视频| 中文亚洲av片在线观看爽 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲第一欧美日韩一区二区三区 | 99re6热这里在线精品视频| 国产成人影院久久av| 国产片内射在线| 亚洲色图av天堂| 国产xxxxx性猛交| 香蕉丝袜av| 午夜老司机福利片| 欧美国产精品va在线观看不卡| 国产亚洲精品第一综合不卡| 久久精品国产亚洲av高清一级| 老司机午夜福利在线观看视频 | 女性生殖器流出的白浆| 在线av久久热| 久久久久网色| 国产又色又爽无遮挡免费看| 久久天躁狠狠躁夜夜2o2o| 法律面前人人平等表现在哪些方面| a级毛片黄视频| 国产不卡一卡二| 69av精品久久久久久 | 国产一区二区三区综合在线观看| 大香蕉久久成人网| 美女福利国产在线| 日本五十路高清| 国产男女超爽视频在线观看| 老司机福利观看| 精品福利观看| 十八禁网站免费在线| 亚洲天堂av无毛| 怎么达到女性高潮| 色在线成人网| 色精品久久人妻99蜜桃| 一级毛片女人18水好多| 韩国精品一区二区三区| 亚洲欧美一区二区三区黑人| 人妻久久中文字幕网| 岛国在线观看网站| 一本一本久久a久久精品综合妖精| 亚洲五月婷婷丁香| svipshipincom国产片| 青青草视频在线视频观看| 亚洲精品久久午夜乱码| 天天躁日日躁夜夜躁夜夜| 亚洲专区国产一区二区| 欧美日韩亚洲国产一区二区在线观看 | 免费少妇av软件| 欧美激情极品国产一区二区三区| av超薄肉色丝袜交足视频| 日本一区二区免费在线视频| 国产精品 欧美亚洲| 免费一级毛片在线播放高清视频 | 中文字幕人妻丝袜一区二区| 亚洲人成伊人成综合网2020| 在线十欧美十亚洲十日本专区| 久久国产精品男人的天堂亚洲| 久久99一区二区三区| www.精华液| 亚洲五月色婷婷综合| 亚洲第一青青草原| 久久久国产成人免费| 男男h啪啪无遮挡| 午夜免费成人在线视频| 久久国产精品大桥未久av| 一本—道久久a久久精品蜜桃钙片| 两性午夜刺激爽爽歪歪视频在线观看 | 91九色精品人成在线观看| 精品亚洲成国产av| 麻豆av在线久日| 久久人人爽av亚洲精品天堂| 亚洲av成人不卡在线观看播放网| 国产成人精品无人区| 成人18禁在线播放| 久久久久久久大尺度免费视频| 女性生殖器流出的白浆| 国产激情久久老熟女| 国产精品亚洲一级av第二区| 99国产精品免费福利视频| 中文字幕人妻熟女乱码| 精品国产国语对白av| 国产成人系列免费观看| 欧美亚洲 丝袜 人妻 在线| 色综合欧美亚洲国产小说| 午夜久久久在线观看| 满18在线观看网站| 99精品久久久久人妻精品| 日韩中文字幕视频在线看片| 丝袜美足系列| 天天躁夜夜躁狠狠躁躁| 精品久久久久久久毛片微露脸| 亚洲精品自拍成人| 黄色丝袜av网址大全| 麻豆成人av在线观看| 热99久久久久精品小说推荐| 51午夜福利影视在线观看| 精品一区二区三区四区五区乱码| 黄色怎么调成土黄色| 国产亚洲一区二区精品| 18禁国产床啪视频网站| 中文欧美无线码| 50天的宝宝边吃奶边哭怎么回事| 在线观看66精品国产| 蜜桃在线观看..| av超薄肉色丝袜交足视频| av一本久久久久| 亚洲午夜精品一区,二区,三区| 一个人免费看片子| 欧美激情高清一区二区三区| 别揉我奶头~嗯~啊~动态视频| 午夜激情久久久久久久| 91精品国产国语对白视频| 免费日韩欧美在线观看| 黑人巨大精品欧美一区二区蜜桃| 深夜精品福利| 亚洲av国产av综合av卡| 人妻久久中文字幕网| 又黄又粗又硬又大视频| 亚洲国产成人一精品久久久| 国产精品免费一区二区三区在线 | 无遮挡黄片免费观看| 91成年电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 9色porny在线观看| 老熟妇乱子伦视频在线观看| 国产在线观看jvid| 51午夜福利影视在线观看| 狂野欧美激情性xxxx| 一边摸一边抽搐一进一出视频| 亚洲国产看品久久| 男女免费视频国产| 久久99热这里只频精品6学生| 午夜福利视频在线观看免费| 免费人妻精品一区二区三区视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产老妇伦熟女老妇高清| kizo精华| av又黄又爽大尺度在线免费看| 中文亚洲av片在线观看爽 | bbb黄色大片| 女人爽到高潮嗷嗷叫在线视频| 精品人妻在线不人妻| 免费观看a级毛片全部| 一本久久精品| 国产一区二区三区综合在线观看| 午夜成年电影在线免费观看| 午夜精品国产一区二区电影| 日韩大码丰满熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产精品影院| 久久中文字幕一级| 日本av手机在线免费观看| 欧美av亚洲av综合av国产av| 久9热在线精品视频| 高清欧美精品videossex| 亚洲av美国av| 国产精品免费大片| 99国产极品粉嫩在线观看| a级毛片在线看网站| 久久热在线av| 在线观看免费视频日本深夜| 好男人电影高清在线观看| 美女扒开内裤让男人捅视频| 欧美乱码精品一区二区三区| 操出白浆在线播放| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产区一区二| 久久狼人影院| videosex国产| 视频在线观看一区二区三区| 女警被强在线播放| 精品高清国产在线一区| 老熟女久久久| 黄色视频,在线免费观看| av片东京热男人的天堂| 一区二区三区乱码不卡18| 久久久国产欧美日韩av| 精品国产乱码久久久久久男人| av片东京热男人的天堂| 国产老妇伦熟女老妇高清| 人妻 亚洲 视频| 色精品久久人妻99蜜桃| 我的亚洲天堂| 男人舔女人的私密视频| 久久久国产精品麻豆| 一本—道久久a久久精品蜜桃钙片| 岛国在线观看网站| 老熟妇仑乱视频hdxx| 精品人妻在线不人妻| 十八禁人妻一区二区| 满18在线观看网站| 丁香六月天网| av网站免费在线观看视频| 高清av免费在线| 黄色成人免费大全| 啦啦啦 在线观看视频| 狂野欧美白嫩少妇大欣赏| 老司机深夜福利视频在线观看| 国产男靠女视频免费网站| 香蕉av资源在线| 女警被强在线播放| 精品福利观看| 国产精品九九99| 国产成人啪精品午夜网站| 国产午夜精品论理片| 中文字幕人成人乱码亚洲影| 毛片女人毛片| 99精品在免费线老司机午夜| 亚洲人成网站高清观看| 我的老师免费观看完整版| 亚洲精品美女久久久久99蜜臀| 淫秽高清视频在线观看| 国产日本99.免费观看| 一本综合久久免费| 97人妻精品一区二区三区麻豆| 国产一区二区三区在线臀色熟女| 男人的好看免费观看在线视频| 99精品在免费线老司机午夜| 我的老师免费观看完整版| 噜噜噜噜噜久久久久久91| 成年人黄色毛片网站| 午夜福利成人在线免费观看| 日韩国内少妇激情av| 国产男靠女视频免费网站| 久久中文字幕人妻熟女| 禁无遮挡网站| 国产又黄又爽又无遮挡在线| 精品一区二区三区视频在线观看免费| 久久精品aⅴ一区二区三区四区| 欧美zozozo另类| 精品99又大又爽又粗少妇毛片 | 999精品在线视频| 综合色av麻豆| 91麻豆av在线| 又黄又爽又免费观看的视频| 美女午夜性视频免费| 日韩av在线大香蕉| 国产69精品久久久久777片 | 国产成人精品久久二区二区91| 国产真实乱freesex| 中文字幕av在线有码专区| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 欧美国产日韩亚洲一区| 国产黄片美女视频| 国产欧美日韩一区二区三| 亚洲国产欧美一区二区综合| e午夜精品久久久久久久| 波多野结衣巨乳人妻| 亚洲真实伦在线观看| 男女之事视频高清在线观看| 国产单亲对白刺激| 网址你懂的国产日韩在线| 亚洲av免费在线观看| 国产亚洲精品久久久com| 国产成人av激情在线播放| 久久久久久人人人人人| 九九热线精品视视频播放| 桃色一区二区三区在线观看| 好看av亚洲va欧美ⅴa在| 国产又色又爽无遮挡免费看| 亚洲欧美日韩卡通动漫| 久久久久精品国产欧美久久久| 三级男女做爰猛烈吃奶摸视频| 欧美黄色片欧美黄色片| 久久久久性生活片| 精品久久久久久久人妻蜜臀av| 国产精品久久久av美女十八| 亚洲av熟女| 国产高清三级在线| 成人特级黄色片久久久久久久| 女生性感内裤真人,穿戴方法视频| 亚洲国产色片| 亚洲欧美日韩东京热| 日韩欧美国产一区二区入口| 国产三级黄色录像| 首页视频小说图片口味搜索| 在线十欧美十亚洲十日本专区| 久久天堂一区二区三区四区| 成人av在线播放网站| 中文字幕av在线有码专区| 97人妻精品一区二区三区麻豆| 日韩高清综合在线| 亚洲av中文字字幕乱码综合| 精品人妻1区二区| 美女cb高潮喷水在线观看 | 国产精品久久久久久人妻精品电影| 可以在线观看毛片的网站| 热99在线观看视频| 久久中文字幕一级| 久久天躁狠狠躁夜夜2o2o| 99热6这里只有精品| 9191精品国产免费久久| 亚洲乱码一区二区免费版| 免费搜索国产男女视频| 亚洲自偷自拍图片 自拍| 午夜a级毛片| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 久久久成人免费电影| 毛片女人毛片| 99热这里只有是精品50| 国产精品98久久久久久宅男小说| 欧美av亚洲av综合av国产av| 成在线人永久免费视频| 国产麻豆成人av免费视频| 免费在线观看日本一区| 每晚都被弄得嗷嗷叫到高潮| 88av欧美| 成年女人毛片免费观看观看9| 成年女人看的毛片在线观看| АⅤ资源中文在线天堂| 国产高清有码在线观看视频| 热99re8久久精品国产| 免费看光身美女| 国产91精品成人一区二区三区| 国产成人av教育| 国产精品乱码一区二三区的特点| 91字幕亚洲| 精品免费久久久久久久清纯| 成人一区二区视频在线观看| 亚洲一区二区三区色噜噜| 叶爱在线成人免费视频播放| 国产精品永久免费网站| 日韩欧美国产一区二区入口| 制服丝袜大香蕉在线| 高潮久久久久久久久久久不卡| 国产三级在线视频| 国产黄色小视频在线观看| 狂野欧美激情性xxxx| 国产精品亚洲美女久久久| 国产高清视频在线观看网站| 精品久久久久久久末码| 99精品在免费线老司机午夜| 午夜免费成人在线视频| 日日干狠狠操夜夜爽| 精品久久久久久久久久久久久| 老汉色∧v一级毛片| 国产精品av视频在线免费观看| av女优亚洲男人天堂 | 天堂√8在线中文| 欧美中文日本在线观看视频| 三级男女做爰猛烈吃奶摸视频| 岛国视频午夜一区免费看| 成人性生交大片免费视频hd| 欧美日韩瑟瑟在线播放| 88av欧美| 久久精品国产99精品国产亚洲性色| 国内少妇人妻偷人精品xxx网站 | 欧美乱码精品一区二区三区| 国产精品一及| 精品久久久久久久末码| 成人一区二区视频在线观看| av国产免费在线观看| 欧美激情在线99| 欧美大码av| 人妻夜夜爽99麻豆av| 欧美成人性av电影在线观看| 中文字幕久久专区| 欧美+亚洲+日韩+国产| 午夜亚洲福利在线播放| 亚洲人成网站高清观看| av天堂在线播放| 99riav亚洲国产免费| 99热精品在线国产| 91av网一区二区| 国产真实乱freesex| h日本视频在线播放| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久亚洲av鲁大| 中国美女看黄片| 国产精品乱码一区二三区的特点| 久久香蕉国产精品| 国产av一区在线观看免费| 无遮挡黄片免费观看| 女警被强在线播放| 国产亚洲av高清不卡| 国产精品女同一区二区软件 | 国内精品一区二区在线观看| 日韩三级视频一区二区三区| 国产亚洲欧美98| 国产精品综合久久久久久久免费| 午夜免费激情av| 国产精品99久久99久久久不卡| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品久久男人天堂| 精品国产超薄肉色丝袜足j| 欧美性猛交黑人性爽| 国产 一区 欧美 日韩| 亚洲激情在线av| 亚洲一区高清亚洲精品| 全区人妻精品视频| 国产一区二区三区在线臀色熟女| 国内精品久久久久精免费| 日韩成人在线观看一区二区三区| 亚洲国产欧美一区二区综合| 成人精品一区二区免费| 国产欧美日韩精品亚洲av| 精品久久久久久成人av| 久久久久久人人人人人| 久久久国产成人免费| 三级男女做爰猛烈吃奶摸视频| 亚洲七黄色美女视频| 日本一本二区三区精品| 午夜福利在线在线| 欧美zozozo另类| 在线观看日韩欧美| 一本综合久久免费| 久久久久久国产a免费观看| 级片在线观看| 男女下面进入的视频免费午夜| 亚洲一区高清亚洲精品| 97碰自拍视频| 亚洲成av人片免费观看| 国产精品1区2区在线观看.| 亚洲av第一区精品v没综合| 亚洲av成人一区二区三| 亚洲中文av在线| 看免费av毛片| 岛国视频午夜一区免费看| 99久久成人亚洲精品观看| 亚洲国产看品久久| 亚洲第一欧美日韩一区二区三区| 亚洲精品一区av在线观看| 在线免费观看不下载黄p国产 | 香蕉av资源在线| 国产久久久一区二区三区| 午夜福利欧美成人| 亚洲第一欧美日韩一区二区三区| 男女午夜视频在线观看| 非洲黑人性xxxx精品又粗又长| 999精品在线视频| 美女高潮喷水抽搐中文字幕| 黄色丝袜av网址大全| 少妇丰满av| 色av中文字幕| 啦啦啦免费观看视频1| 制服人妻中文乱码| 一个人免费在线观看的高清视频| 手机成人av网站| 久久久色成人| 欧美av亚洲av综合av国产av| 久久人妻av系列| 久久精品综合一区二区三区| 亚洲午夜精品一区,二区,三区| 青草久久国产| 久久人人精品亚洲av| 精品久久久久久久久久久久久| 亚洲精品国产精品久久久不卡| 又黄又粗又硬又大视频| 日本一二三区视频观看| 人人妻人人澡欧美一区二区| 国产人伦9x9x在线观看| or卡值多少钱| 一二三四在线观看免费中文在| 成人av在线播放网站| 天天躁狠狠躁夜夜躁狠狠躁| 黑人欧美特级aaaaaa片| 舔av片在线| 在线观看66精品国产| 亚洲人与动物交配视频| 久久亚洲精品不卡| 国产精品野战在线观看| 91久久精品国产一区二区成人 | 亚洲乱码一区二区免费版| 伦理电影免费视频| 国产一区二区三区在线臀色熟女| 国产精品98久久久久久宅男小说| 国产精品美女特级片免费视频播放器 | 久久午夜亚洲精品久久| 欧美日本亚洲视频在线播放| 国产又黄又爽又无遮挡在线| av在线天堂中文字幕| 黄色日韩在线| 中文字幕精品亚洲无线码一区| av欧美777| 国产欧美日韩精品一区二区| 91久久精品国产一区二区成人 | 2021天堂中文幕一二区在线观| 国内久久婷婷六月综合欲色啪| 国产黄a三级三级三级人| 少妇熟女aⅴ在线视频| 免费观看的影片在线观看| 免费看a级黄色片| 天堂影院成人在线观看| 99国产综合亚洲精品| 变态另类成人亚洲欧美熟女| 很黄的视频免费| 久久这里只有精品中国| 国产伦精品一区二区三区四那| 国产成人精品久久二区二区91| 免费在线观看影片大全网站| 精品午夜福利视频在线观看一区| 亚洲精品美女久久av网站| 757午夜福利合集在线观看| 1024手机看黄色片| 国产一区二区在线观看日韩 | 国产成+人综合+亚洲专区| a在线观看视频网站| 国产三级在线视频| 黑人欧美特级aaaaaa片| 18禁美女被吸乳视频| 午夜福利在线观看免费完整高清在 | 天天添夜夜摸| 日韩中文字幕欧美一区二区| av片东京热男人的天堂| 美女高潮喷水抽搐中文字幕| 日本熟妇午夜| 1024手机看黄色片| 在线观看一区二区三区| 国产激情偷乱视频一区二区| 精品久久久久久久末码| 国产男靠女视频免费网站| 最近最新中文字幕大全电影3| 日韩欧美免费精品| 亚洲精品在线美女| 特大巨黑吊av在线直播| 久久久水蜜桃国产精品网| 成年版毛片免费区| 精品一区二区三区四区五区乱码| 色综合亚洲欧美另类图片| 欧美乱码精品一区二区三区| 男插女下体视频免费在线播放| 在线视频色国产色| 怎么达到女性高潮| 99久久无色码亚洲精品果冻| 亚洲精品在线观看二区| 国产在线精品亚洲第一网站| 性色av乱码一区二区三区2| 色哟哟哟哟哟哟| 校园春色视频在线观看| 国产亚洲欧美在线一区二区| 欧美日韩综合久久久久久 | 男人舔女人的私密视频| av黄色大香蕉| 香蕉丝袜av| av在线蜜桃| 午夜成年电影在线免费观看| 国产精品精品国产色婷婷| 亚洲中文字幕一区二区三区有码在线看 | 每晚都被弄得嗷嗷叫到高潮| 国产成人啪精品午夜网站| 亚洲狠狠婷婷综合久久图片| 国产午夜福利久久久久久| 色av中文字幕| 露出奶头的视频| 精品国产超薄肉色丝袜足j| 色av中文字幕| 国产精品爽爽va在线观看网站| 国产亚洲精品久久久久久毛片| 老司机午夜福利在线观看视频| 国产精品九九99| 在线十欧美十亚洲十日本专区| 国模一区二区三区四区视频 | 国产精品日韩av在线免费观看| 欧美乱妇无乱码| 国内精品美女久久久久久| 女警被强在线播放| 村上凉子中文字幕在线| 国产高清videossex| 露出奶头的视频| 黄色日韩在线| 女警被强在线播放| 最新中文字幕久久久久 | 国产毛片a区久久久久| 免费人成视频x8x8入口观看| 在线a可以看的网站| 亚洲色图 男人天堂 中文字幕| 男人舔女人的私密视频| 免费看日本二区| 国产精品一及| 日韩欧美在线乱码| 亚洲国产精品成人综合色| 久久精品91蜜桃| 久久这里只有精品19|