亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        降低RS碼算法復(fù)雜度的改進(jìn)KV算法

        2014-08-10 08:09:54江南
        宜賓學(xué)院學(xué)報(bào) 2014年6期
        關(guān)鍵詞:重?cái)?shù)碼元碼字

        江南

        (福建信息職業(yè)技術(shù)學(xué)院軟件工程系,福建福州350003)

        降低RS碼算法復(fù)雜度的改進(jìn)KV算法

        江南

        (福建信息職業(yè)技術(shù)學(xué)院軟件工程系,福建福州350003)

        提出對(duì)Koetter-Vardy(KV)算法進(jìn)行改進(jìn)后的重編碼算法,利用Reed-Solomon(RS)碼的線性性質(zhì)對(duì)重?cái)?shù)矩陣進(jìn)行預(yù)處理,改進(jìn)了插值算法的初始多項(xiàng)式條件,降低插值算法的復(fù)雜度,從而降低了KV算法的總體復(fù)雜度,帶來的復(fù)雜度的節(jié)省因子是n2(n-k)2.對(duì)該算法的軟件實(shí)現(xiàn)以及仿真結(jié)果顯示,對(duì)高碼率的RS碼,重編碼算法幾乎不犧牲譯碼性能.

        KV算法;重編碼算法;RS碼;復(fù)雜度

        KV算法是Koetter和Vardy于2003年在Guruswami-Sudan算法的基礎(chǔ)上提出來的[1-2]、實(shí)現(xiàn)Reed-Solomon(RS)碼代數(shù)軟譯碼的算法簡(jiǎn)稱.它將信道軟信息轉(zhuǎn)化成代數(shù)約束條件,然后再利用GS算法的部分步驟,相對(duì)于傳統(tǒng)的Berlekamp-Massey硬判決算法有著0.28~1.23 dB的增益.假設(shè)RS碼碼長(zhǎng)為n,信息長(zhǎng)度為k,KV算法復(fù)雜度與n成多項(xiàng)式關(guān)系,是目前較為有效的RS碼軟譯碼算法.經(jīng)過改進(jìn)后的Koetter-Vardy算法[3-4]可以達(dá)到更好的譯碼效果,但其復(fù)雜度也不能忽視.重編碼算法[5]是降低KV算法復(fù)雜度的重要算法,能夠在幾乎不損失性能的情況下將復(fù)雜度較高的插值算法的復(fù)雜度降低n2(n-k)2因子,有著廣泛的潛在應(yīng)用[6].

        1 Koetter-Vardy算法實(shí)現(xiàn)步驟

        1.1 重?cái)?shù)分配

        發(fā)送符號(hào)從有限域GF(Q)取出,在無記憶信道傳輸.信道輸入輸出是隨機(jī)變量X和Y.觀察信道輸出的軟信息是后驗(yàn)概率:

        其中1≤i≤Q,1≤j≤n.這個(gè)可靠性矩陣Π就是KV算法的輸入.

        從可靠性矩陣計(jì)算出一個(gè)(Q×n)的整數(shù)矩陣M,此矩陣稱為重?cái)?shù)矩陣,矩陣中的元素用mi,j表示.其Q行代表一個(gè)符號(hào)的Q種不同取值,即αi,1≤i≤Q.n列代表一個(gè)碼字中n個(gè)符號(hào)位置,即支持集中的xj,1≤j≤n.文本只考慮Gross提出的簡(jiǎn)單重?cái)?shù)分配方案[7],即這是目前KV算法研究中較為常用的方案,更好的方案見文獻(xiàn)[8].

        1.2 插值

        尋找一個(gè)二元多項(xiàng)式Q(x,y),使其滿足下面兩個(gè)條件:①對(duì)重?cái)?shù)矩陣中每個(gè)非零的元素mi,j,Q(x,y)以重?cái)?shù)mi,j通過GF(Q)×GF(Q)平面上的點(diǎn)(xj,αi),或者說(xj,αi)是Q(x,y)的mi,j重零點(diǎn);②Q(x,y)的(1,k-1)-加權(quán)階數(shù)盡可能?。蠹訖?quán)階數(shù)盡可能小的原因是這樣算法的搜索半徑最大,糾錯(cuò)能力最強(qiáng).該步驟是KV算法復(fù)雜度的主要來源,本文所要介紹的重編碼算法也就是要降低這一步算法的復(fù)雜度.不難看出,上述條件①說明二元多項(xiàng)式Q(x,y)需要滿足C(M)=個(gè)線性約束條件.Koetter提出一種插值算法[9],其基本步驟如下,假定C=C(M):

        (1)將二元多項(xiàng)式的集合FL[x,y]按照首單項(xiàng)式的y-階數(shù)進(jìn)行劃分,劃分成L個(gè)子集,初始化L個(gè)多項(xiàng)式.這里FL[x,y]指的是所有y階數(shù)低于L的二元多項(xiàng)式.

        (2)每次迭代,更新分別來自上面劃分出來的L個(gè)集合的L個(gè)多項(xiàng)式,i次迭代后,L個(gè)多項(xiàng)式都能滿足第1到i個(gè)線性約束,而且在各自的所在的劃分集合中,它是最小的.

        (3)C次迭代后,L個(gè)多項(xiàng)式都能滿足C個(gè)條件,并且在各自集合中加權(quán)階數(shù)最低,在這L個(gè)里面再選出最小的那個(gè),作為結(jié)果輸出.

        1.3 因式分解

        找出二元多項(xiàng)式Q(x,y)的所有形式為y-f(x)的因式,degf(x)<k,這些因式是預(yù)選的消息多項(xiàng)式,將他們重新編碼,就得到預(yù)選碼字列表.

        1.4 ML判決

        根據(jù)極大似然估計(jì)Maximum Likelihood準(zhǔn)則,在列表上尋找一個(gè)與接收矢量y→最接近的碼字,作為譯碼結(jié)果輸出.

        從以上步驟中可以看出,KV算法的復(fù)雜度在很大程度上取決于插值的復(fù)雜度,而插值的復(fù)雜度則取決于所需要滿足的線性條件的總個(gè)數(shù).因此,降低插值算法的復(fù)雜度就可以降低KV算法的復(fù)雜度,從而提高算法效率.

        2 重編碼算法

        重編碼算法是降低KV算法復(fù)雜度的重要算法,它能夠在幾乎不損失性能的情況下將插值算法的復(fù)雜度降低n2(n-k)2因子.下面先介紹其步驟和原理,再通過算例和仿真結(jié)果做進(jìn)一步了解.

        重編碼算法利用了RS碼的線性性質(zhì)和廣義系統(tǒng)碼編碼,線性分組碼具有線性性質(zhì),即如果是碼字集合,則這個(gè)性質(zhì)在譯碼中可以得到應(yīng)用,比如在硬判決譯碼中,可以將接收碼字加上一個(gè)碼字,或者可以叫做偏移量,則得到的結(jié)果仍然是一個(gè)碼字,對(duì)其進(jìn)行譯碼,譯碼輸出再加上(GF(2m)上加減法是一樣的),就可以得到原接收碼字的譯碼輸出.這是因?yàn)?,和的錯(cuò)誤圖樣是一樣的.設(shè)發(fā)送碼字為,錯(cuò)誤圖樣為,則:

        下面討論RS碼的系統(tǒng)碼編碼.所謂系統(tǒng)碼,就是消息碼元“顯式”地在編碼后碼字中出現(xiàn).常提到的系統(tǒng)碼,消息碼元都是連續(xù)地在碼字的前k位或者后k位出現(xiàn),這叫做嚴(yán)格的系統(tǒng)碼.這里要用到的是,k個(gè)編碼前碼元任意地出現(xiàn)在碼字中任何k個(gè)位置的系統(tǒng)碼,可以稱之為廣義系統(tǒng)碼.這種廣義系統(tǒng)碼編碼也可以看作已知一個(gè)碼字任意k個(gè)位置的碼元,要求出這個(gè)碼字.由于RS碼是Maximum Distance Separable碼,總是可以用一個(gè)糾刪譯碼器來求出這個(gè)碼字.

        重編碼算法的步驟是:

        (1)當(dāng)譯碼器接收到軟信息時(shí),首先做硬判決,得到硬判決碼字r→.

        (2)由可靠性矩陣Π求出重?cái)?shù)矩陣M,這里采用簡(jiǎn)單的分配方案,即

        (3)根據(jù)的M元素值,選擇出含有最大重?cái)?shù)mmax(即可靠性較高)的k個(gè)碼元位置,設(shè)為集合R,其余碼元位置設(shè)為集合U(假定可以找到k個(gè)這樣的位置,找不滿k個(gè)的情況在后面討論).

        (4)將r→中k個(gè)對(duì)應(yīng)于R中元素位置的碼元看作消息碼元,對(duì)其進(jìn)行重編碼,使編碼后碼字在這k個(gè)位置上與相同,將此碼字記為(用了廣義系統(tǒng)碼編碼).

        (6)對(duì)重?cái)?shù)矩陣插值,這樣處理的重?cái)?shù)矩陣可以減少插值算法的迭代次數(shù).

        (7)對(duì)插值多項(xiàng)式進(jìn)行因式分解,后面的步驟和原KV算法相同.

        分析重編碼算法的原理:第(5)步中對(duì)重?cái)?shù)矩陣的偏移類似于上面提到的對(duì)接收硬判碼字進(jìn)行偏移.只不過,KV是軟判決譯碼算法,這里要偏移的不僅是一個(gè)碼字,而是整個(gè)重?cái)?shù)矩陣.重?cái)?shù)矩陣的Q行代表著GF(Q)上的符號(hào),n列代表一個(gè)碼字中n個(gè)碼元的位置.矩陣中某列各行代表的符號(hào)加上ψ→中對(duì)應(yīng)位置的符號(hào),就得到新的符號(hào),該符號(hào)又對(duì)應(yīng)了新的矩陣行號(hào).矩陣中每個(gè)元素都要按照這種關(guān)系進(jìn)行變換.例如,對(duì)GF(8),(7,5)RS碼,M中某個(gè)元素:m7,3=3,說明第三個(gè)符號(hào)位置,對(duì)應(yīng)于α6的重?cái)?shù)為3,而ψ→的第三個(gè)符號(hào)為α1,α1+α6=α5,則m7,3的值移到m6,3,原先m6,3的值也根據(jù)此規(guī)則進(jìn)行移動(dòng),如果把行號(hào)視為正整數(shù)集合:NQ:{0,1,…,Q-1},則可以把這種變換關(guān)系視為一個(gè)映射NQ?NQ,這個(gè)映射是一一映射,也就是說不會(huì)出現(xiàn)某個(gè)值被其它值覆蓋的情況.

        再看這種偏移帶來的效果.重?cái)?shù)矩陣中最大的重?cái)?shù)為mmax,用m來簡(jiǎn)化表示.重?cái)?shù)為m的位置說明了該位置可靠性較高.在重?cái)?shù)矩陣對(duì)應(yīng)的k列中,m所在的行對(duì)應(yīng)的符號(hào)就是硬判決的符號(hào).對(duì)矩陣做偏移時(shí),偏移碼字ψ→在這k個(gè)可靠符號(hào)位置上和硬判決碼字是一樣的(因?yàn)樽龅氖菑V義系統(tǒng)碼編碼).這樣,偏移之后,這k個(gè)m就被移到了第0行,即對(duì)應(yīng)GF(Q)零元所在的行上.注意到,由于m是最大的重?cái)?shù),把C(M)'/C(M)個(gè)m移到零元位置就相當(dāng)于把相當(dāng)多的重?cái)?shù)移到了y分量為零的插值點(diǎn)上,這個(gè)插值點(diǎn)集合可以記為:

        對(duì)這些插值點(diǎn)的插值,可以直接求得滿足這些點(diǎn)的多項(xiàng)式.

        對(duì)PR中的k個(gè)重?cái)?shù)均為m點(diǎn),一個(gè)滿足這些插值重?cái)?shù)要求且次數(shù)最低的多項(xiàng)式是p(x)m,其中

        同時(shí),不難看出,p(x)m-jyj也是滿足各點(diǎn)的重?cái)?shù)要求,且次數(shù)最低,因此可以將本文第一部分插值算法中的初始條件設(shè)為:

        這樣,這些初始條件已經(jīng)滿足(1-k n) C(M)個(gè)重?cái)?shù)為m的插值點(diǎn)的要求,后面的插值就不需要再在這些點(diǎn)上進(jìn)行迭代,這就是重編碼算法降低復(fù)雜度的原理.

        上面假定了重?cái)?shù)矩陣中存在k個(gè)重?cái)?shù)為m的點(diǎn),如果沒有這么多個(gè)點(diǎn)該怎么辦?在此做個(gè)近似,將重?cái)?shù)m-1的點(diǎn)(如果還不夠的話就采用重?cái)?shù)m-2的點(diǎn))的重?cái)?shù)升為m,這樣就一樣能使用重編碼算法了.當(dāng)然,這樣的修改重?cái)?shù)可能會(huì)帶來一些性能上的降低,但是從仿真結(jié)果上看,幾乎沒有什么影響,特別是在高信噪比的情況下,這種修正重?cái)?shù)的情況是很少的.

        重編碼算法能降低插值算法的復(fù)雜度,而插值算法的復(fù)雜度使用線性約束條件的個(gè)數(shù)C(M)來表示的.采用Gross的重?cái)?shù)分配方案,C(M)的上界是:

        這實(shí)際上是將每個(gè)硬判決符號(hào)分配重?cái)?shù)C(M)'的GS算法的C(M).而現(xiàn)在,重?cái)?shù)矩陣?yán)锶サ袅薻個(gè)m,C(M)'的上界是:

        這樣可以估計(jì)出迭代次數(shù)的節(jié)省,在高碼率情況下是很可觀的.后面通過仿真能夠看到更直觀的數(shù)據(jù).

        另外,這里還討論一下“可靠性較高”的含義.選取重編碼位置的集合R時(shí),選取k個(gè)可靠性較高的位置作為重編碼位置.這里的可靠性較高,并不是要求這些位置上的碼元實(shí)現(xiàn)了無誤傳輸,而只是說在接收機(jī)看來,可靠性較高而已.從前文提到的線性分組碼性質(zhì)可以知道,其實(shí)隨便設(shè)定一個(gè)偏移碼字ψ→,對(duì)重?cái)?shù)矩陣按照這個(gè)偏移碼字進(jìn)行偏移,再進(jìn)入KV算法譯碼,譯碼后再偏移回來,如果原先碼字是KV算法可譯的,這樣偏移一下也一樣還是可以譯出來的.選擇可靠性較高是因?yàn)檫@些位置對(duì)應(yīng)有最大的重?cái)?shù)m,能夠最大程度地把最大重?cái)?shù)移動(dòng)到y(tǒng)坐標(biāo)為零的點(diǎn)上,降低插值復(fù)雜度.換句話說,如果這些“可靠性高”的碼字實(shí)際上是有誤的,但通過KV算法,仍然能夠把它糾正.

        3 算例及性能仿真

        下面的算例突出了和重編碼有關(guān)的步驟.考慮一個(gè)(7,5)碼,消息碼字是:

        編碼后發(fā)送碼字為:

        經(jīng)過BPSK調(diào)制,AWGN信道,SNR=4.0.

        從可靠性矩陣能得到硬判決碼字:

        可以看到,這里有2個(gè)錯(cuò)誤.

        設(shè)C(M),用Gross重?cái)?shù)分配方法得到重?cái)?shù)矩陣:

        可以選擇前5個(gè)重?cái)?shù)均為4的位置作為重編碼位置集合R,即

        以此為偏移量,對(duì)重?cái)?shù)矩陣偏移,得到:

        同時(shí),可以計(jì)算出:

        將M'首行5個(gè)4全部置為0,同時(shí)將(2)作為插值初始條件,進(jìn)行插值,可以得到多項(xiàng)式QM(x,y),該多項(xiàng)式較長(zhǎng),這里不再寫出.進(jìn)行因式分解,可以得到2個(gè)預(yù)選消息碼字:

        將它們編碼后:

        這時(shí)候,要將它們恢復(fù)到偏移前的碼字,然后再進(jìn)行ML判決:

        進(jìn)行ML判決選出c2作為譯碼輸出,譯碼成功.

        圖1和圖2分別仿真了RS(15,11)和RS(63,55)兩種碼型,ReEn表示重編碼算法,KV表示原KV算法,復(fù)雜度因子m一律取4.可以看到,重編碼和原KV算法的性能確實(shí)是十分接近的,在長(zhǎng)碼字表現(xiàn)得甚至略好些,這說明,使用重編碼算法在性能上是不會(huì)受到損失的.

        圖1 BPSK調(diào)制,AWGN信道下,RS(15,11)碼重編碼算法性能

        圖2 BPSK調(diào)制,AWGN信道下,RS(63,55)碼重編碼算法性能

        仿真表明了復(fù)雜度降低的情況.設(shè)RS碼為(63,55)的RS碼,碼率為0.87,通過BPSK調(diào)制和AWGN信道后用KV算法進(jìn)行譯碼,在每個(gè)信噪比上仿真100個(gè)碼字,取插值迭代次數(shù)的平均值來估計(jì)實(shí)際迭代次數(shù),得到表1(C(M)為原KV算法迭代次數(shù),C(M)'為重編碼算法迭代次數(shù)).

        表1 重編碼算法與原KV算法迭代次數(shù)對(duì)比

        從表1看出,實(shí)際上,迭代次數(shù)的節(jié)省比上面用上界得出的結(jié)論還要多,基本達(dá)到了90%的節(jié)省,這是很可觀的,當(dāng)然,要在碼率較高的時(shí)候才會(huì)有如此可觀的節(jié)省.插值的復(fù)雜度和迭代次數(shù)C呈二次方關(guān)系[10],注意到重編碼以后C已經(jīng)減少了1-k n因子,因此,重編碼帶來的復(fù)雜度的節(jié)省因子是n2(n-k)2.

        4 結(jié)語

        本文討論了KV算法的一種簡(jiǎn)單而且有效的改進(jìn)——重編碼算法.對(duì)高碼率的RS碼來說,使用重編碼算法可以大大降低復(fù)雜度,同時(shí)幾乎不犧牲性能.KV算法的復(fù)雜度雖然是多項(xiàng)式時(shí)間的,但和傳統(tǒng)的硬判決算法相比還是復(fù)雜許多,這成為了它實(shí)際應(yīng)用的一個(gè)瓶頸.重編碼算法大大降低了插值算法的復(fù)雜度,可以想象,在大部分應(yīng)用KV算法的實(shí)際系統(tǒng)中,該算

        法將有廣泛的應(yīng)用.最新研究成果表明[11-12],KV算法除了在信道編碼方面的應(yīng)用以外,還可以用于分布式信源編碼,上述論文并未提及重編碼算法.未來研究方向可以朝包括研究如何將重編碼算法應(yīng)用到分布式信源編碼上的方向發(fā)展.

        [1]Koetter R,Vardy A.Algebraic soft-decision decoding of Reed-Solomon codes[J].Information Theory,IEEE Transactions on,2003,49(11): 2809-2825.doi:10.1109/TIT.2003.819332.

        [2]丁溯泉,楊知行,潘長(zhǎng)勇.RS碼軟判決譯碼算法研究的最新進(jìn)展[J].電子科學(xué)技術(shù)評(píng)論,2005(2):37-41.

        [3]Jiang J,Narayanan K R.Algebraic soft-decision decoding of Reed-Solomon codes using bit-level soft information[C].Proc Allerton Conference on Communications,Control and Computing,2006.doi:10.1109/ TIT.2008.928238.

        [4]El-Khamy M,McEliece R J.Iterative Algebraic soft-decision List decoding of Reed-Solomon codes[J].IEEE Journal on Selected Areas in Communications,2006,24(3):481-490.doi:10.1109/JSAC.2005.862399.

        [5]Koetter R,Vardy A.A complexity reducing transformation in algebraic list-decoding of Reed-Solomon codes[C].Proc.IEEE Inform.Theory Workshop,Paris,France,April 2003.doi:10.1109/ITW.2003.1216682.

        [6]Koetter R,Ma J,Vardy A.The re-encoding transformation in Algebraic list-decoding of Reed-Solomon codes[J].Information Theory,IEEE Transactions on,2011,57(2):633-647.

        [7]Gross W J,Kschischang F R,Koetter R,et al.Applications of algebraic soft-decision decoding of Reed-Solomon codes[J].IEEE Trans Communication,2006,54(7):1224-1234.doi:10.1109/TCOMM.2006.877972.

        [8]Huang Q,Wu J,Zhao C M,et al.Waterfilling-like multiplicity assignment algorithm for algebraic soft-decision decoding of Reed-Solomon codes[C].IEEE International Conference on Communications(ICC), 2007.doi:10.1109/ICC.2007.1028.

        [9]Koetter R.Fast generalized minimum-distance decoding of algebraicgeometry and Reed-Solomon codes[J].IEEE Trans Inform Theory, 1996,42(3):721-737.doi:10.1109/18.490540.

        [10]Cassuto Y,Bruck J;McEliece R J.On the average complexity of Reed-Solomon list decoders[J].Information Theory,IEEE Transactions on, 2013,59(4):2336-2351.doi:10.1109/TIT.2012.2235522.

        [11]Li S Z,Ramamoorthy A.Algebraic code design for Slepian-Wolf codes [C].IEEE International Symposium on Information Theory(ISIT),Saint Petersburg,Russia,2011:1861-1865.

        [12]Ali M,Kuijper M.Source coding with side information using list decoding[J].Information Theory Proceedings(ISIT),IEEE International Symposium on,2010:91-95.doi:10.1109/ISIT.2010.5513280.

        【編校:王露】

        Improved KV Algorithm with RS Code Algorithm Complexity Reduction

        JIANG Nan
        (Fujian Polytechnic of Information Technology,Fuzhou,Fujian 350003,China)

        The reencoding algorithm,an improvement on the Koettery-Vardy(KV)algorithm for Reed-Solomon(RS)soft decoding was introduced.This algorithm takes advantage of the linearity of the Reed-Solomon codes.It preprocesses the multiplicity matrix,improve the initial polynomials in the interpolation algorithm and reduces the complexity of the interpolation algorithm. The overall complexity of the KV algorithm is reduced by a factor ofn2(n-k)2.Observed from the software implementation and simulations,for high rate RS codes,the algorithm reduces the complexity of the KV algorithm without sacrificing much performance.

        Koetter-Vardy(KV)algorithm;reencoding algorithm;Reed-Solomon(RS)codes;complexity

        TP301.6

        A

        1671-5365(2014)06-0094-05

        2013-09-10修回:2013-10-20

        2012年度福建省A類科技項(xiàng)目計(jì)劃立項(xiàng)(JA12385)

        江南(1972-),女,副教授,工學(xué)學(xué)士,研究方向?yàn)橛?jì)算機(jī)軟件及項(xiàng)目管理

        時(shí)間:2013-10-30 15:06

        http://www.cnki.net/kcms/detail/51.1630.Z.20131030.1506.007.html

        猜你喜歡
        重?cái)?shù)碼元碼字
        C3型李代數(shù)的張量積分解
        微分在代數(shù)證明中的兩個(gè)應(yīng)用
        A3型李代數(shù)的張量積分解
        LFM-BPSK復(fù)合調(diào)制參數(shù)快速估計(jì)及碼元恢復(fù)
        以較低截?cái)嘀財(cái)?shù)分擔(dān)超平面的亞純映射的唯一性問題
        放 下
        數(shù)據(jù)鏈系統(tǒng)中軟擴(kuò)頻碼的優(yōu)選及應(yīng)用
        放下
        基于極大似然準(zhǔn)則的短猝發(fā)信號(hào)盲解調(diào)
        一種碼元同步時(shí)鐘信號(hào)的提取方法及單片機(jī)實(shí)現(xiàn)
        日本一区二区精品88| 亚洲成人av在线蜜桃| 西西人体大胆视频无码| 中文字幕久久国产精品| 色呦呦九九七七国产精品| 久久综合九色综合97欧美| 国产成人精品日本亚洲18| 亚洲色欲色欲欲www在线| 亚洲日本精品一区久久精品| 亚洲丰满熟女乱一区二区三区| 人妻激情偷乱视频一区二区三区| 人妻无码中文字幕| 亚洲一区爱区精品无码| 一区二区三区蜜桃在线视频| 日本视频在线观看一区二区| 国产特级毛片aaaaaa高潮流水| 亚洲美腿丝袜 欧美另类| 狠狠色综合播放一区二区| 粉嫩的18在线观看极品精品| 亚洲av网一区二区三区| 欧美日韩一区二区综合| 亚洲国产精品悠悠久久琪琪| 少妇我被躁爽到高潮在线影片| 人妻体内射精一区二区三区| 国产精品国产成人国产三级| 男女视频在线一区二区| 我揉搓少妇好久没做高潮| 亚洲av综合av成人小说| 91spa国产无码| 久久精品国产亚洲av蜜桃av| 免费一级淫片日本高清| 午夜无码国产理论在线| 精品国偷自产在线不卡短视频| 久久少妇高潮免费观看| 综合色就爱涩涩涩综合婷婷| 99久久人人爽亚洲精品美女| 久久青草伊人精品| 久久这黄色精品免费久| 亚洲αv在线精品糸列| 欧美肥胖老妇做爰videos| 正在播放淫亚洲|