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

通信化雷達(dá)信息分段波形設(shè)計(jì)與處理

2023-05-25 09:12:04王漻鯤施龍飛劉甲磊關(guān)一夫
關(guān)鍵詞:信號檢測信息

王漻鯤,施龍飛,劉甲磊,關(guān)一夫

(1.國防科技大學(xué)電子科學(xué)學(xué)院,湖南長沙 410073;2.國防科技大學(xué)電子信息系統(tǒng)復(fù)雜電磁環(huán)境效應(yīng)國家重點(diǎn)實(shí)驗(yàn)室,湖南長沙 410073)

0 引 言

通過在發(fā)射信號中嵌入發(fā)射站位置、發(fā)射時(shí)刻等輔助信息,通信化雷達(dá)能夠?qū)崿F(xiàn)無通信鏈路情況下雙基地探測能力,廣泛適用于分布式探測場景。文獻(xiàn)[1]提出了通信化雷達(dá)架構(gòu),給出了基本的波形設(shè)計(jì)、信號處理和目標(biāo)定位方法。

典型的通信化雷達(dá)發(fā)射波形的嵌入信息一般包括校驗(yàn)碼、發(fā)射站位置、發(fā)射時(shí)刻、發(fā)射波束指向等。現(xiàn)有雷達(dá)通信一體化技術(shù)中的波形設(shè)計(jì)方法和思路可供借鑒。如2015年,Liu 等人[2]將最小頻移鍵控技術(shù)(MSK)與LFM 技術(shù)結(jié)合,提出MSKLFM 信號。文獻(xiàn)[3]通過在LFM 信號中嵌入一組正交FM 項(xiàng)來標(biāo)識通信符號,并通過引入加權(quán)系數(shù)來平衡雷達(dá)探測和通信性能;文獻(xiàn)[4]通過前導(dǎo)碼長度設(shè)置取代功率分配方法來平衡雷達(dá)探測和通信傳輸?shù)男阅埽晃墨I(xiàn)[5]將信息序列用連續(xù)相位調(diào)制的方法調(diào)制在多相調(diào)頻編碼-線性調(diào)頻信號(PCFM)上[6-7],通過調(diào)整指數(shù)來調(diào)整因信息調(diào)制對原始雷達(dá)波形頻率變化的影響;2017年,Liu等人[8]設(shè)計(jì)了一種自適應(yīng)OFDM 雷達(dá)通信一體化波形。總的來說,雷達(dá)通信一體化技術(shù)的出發(fā)點(diǎn)是使雷達(dá)同時(shí)具備目標(biāo)探測與通信的功能[9-10],而這兩種功能在頻譜、能量資源及性能上是不可避免的相互競爭關(guān)系;通信化雷達(dá)的設(shè)計(jì)出發(fā)點(diǎn)是雷達(dá)探測,其有限的通信功能是服務(wù)于雷達(dá)探測功能的,故雷達(dá)通信一體化的波形設(shè)計(jì)難以滿足通信化雷達(dá)的探測需求。

尤其是考慮到機(jī)載平臺(tái)中、高重頻應(yīng)用條件下,脈沖寬度受限,文獻(xiàn)[1]中所設(shè)計(jì)的通信化雷達(dá)波形無法在一個(gè)脈沖內(nèi)嵌入完整的信息量(文中為52 比特,實(shí)際應(yīng)用中可能需要傳輸更大的信息量),因此需要完善通信化雷達(dá)波形的設(shè)計(jì)方法,使其能夠在多種應(yīng)用場合具有更強(qiáng)的適用性。

本文提出了通信化雷達(dá)信息分段波形設(shè)計(jì)思路以及相應(yīng)的信號處理方法,實(shí)現(xiàn)了中、高重頻場景下通信化雷達(dá)波形的信息傳輸和目標(biāo)探測。

1 通信化雷達(dá)信息分段波形設(shè)計(jì)及其信號處理

為了兼顧信息傳遞和目標(biāo)探測需求,信息分段波形采用內(nèi)外兩層正交復(fù)合調(diào)制,將信息序列映射到內(nèi)層的正負(fù)線性調(diào)頻信號,并通過在外層增加13位巴克碼的相位調(diào)制使波形滿足目標(biāo)探測需求。

1.1 信息分段波形信號模型

信息分段波形的信號模型為

式中,N為信息分段數(shù),每段包含13 比特信息,Ts為每比特信息對應(yīng)的子脈寬,Iin(i=1,…,N;n=1,…,13)表示信息段Ii的第n位信息序列。信息分段波形的外層調(diào)制φn為13 位巴克碼bc:[+++++--++-+-+]對應(yīng)的相位序列,其中“+”對應(yīng)相位0,“-”對應(yīng)相位π;內(nèi)層調(diào)制s0(t)、s1(t)為一對正負(fù)線性調(diào)頻信號,如式(2)所示,信息“0”對應(yīng)負(fù)線性調(diào)頻信號s0(t),信息“1”對應(yīng)正線性調(diào)頻信號s1(t),正負(fù)線性調(diào)頻信號具有良好的正交性,同時(shí)線性調(diào)頻體制也有利于信號處理的信噪比提升。

式中,K為正負(fù)線性調(diào)頻信號調(diào)頻斜率,

信息分段波形的結(jié)構(gòu)示意圖如圖1所示,假設(shè)完整的信息序列需要分為4 段進(jìn)行傳輸。為了在后續(xù)信號處理中能夠通過積累提高信息提取及目標(biāo)探測的準(zhǔn)確性,發(fā)射時(shí)將每段信號依次重復(fù)發(fā)射M次,即先發(fā)射M次嵌入信息段I1的脈沖,然后發(fā)射M次嵌入信息段I2的脈沖,以此類推直至將所有信息段發(fā)射完。

圖1 信息分段波形示意圖

1.2 信號處理方法

信號處理具體的流程框圖如圖2所示,完整的信息序列分為N段進(jìn)行傳輸,r11,…,r1M為信息段I1的M次回波脈沖,rN1,…,rNM為信息段IN的M次回波脈沖。

圖2 信息分段波形信號處理流程框圖

由圖2可以得到信號處理流程主要包括以下4個(gè)步驟:

步驟1(第一層匹配) 將嵌入同一信息段Ii,i=1,2,…,N的M個(gè)回波rij(t),j=1,2,…,M分別與正負(fù)線性調(diào)頻信號s0(t)、s1(t)進(jìn)行匹配,得到單次回波信息提取結(jié)果,以信息段I1的單次信息提取為例:

式中,r1j(t)代表了嵌入信息段I1的第j次回波脈沖。r1j0代表回波r1j(t)中信息位“0”的提取,同理r1j1代表信息位“1”的提取。即對信息位“0”“1”分別進(jìn)行了M次提取。

步驟2(第一次積累) 將同一信息段Ii的M次回波單次信息提取結(jié)果進(jìn)行積累,綜合得到信息段Ii中信息位“0”和信息位“1”的提取結(jié)果rIi0,rIi1。

步驟3(第二層匹配濾波) 信息段Ii的第一層匹配濾波輸出結(jié)果rij1,rij0相加得到sij(t),將sij(t)與由13 位巴克碼組成的梳齒濾波器h(t)進(jìn)行第二層匹配濾波,得到單次目標(biāo)檢測結(jié)果sijo(t)。

sij(t)=rij1+rij0,i=1,2,…,N,j=1,2,…,M(5)以嵌入信息段I1回波的單次目標(biāo)檢測結(jié)果s1jo(t)為例:

式中,Ts為每比特信息對應(yīng)的子脈寬。

步驟4(第二次積累) 將同一信息段Ii的M次回波單次目標(biāo)檢測結(jié)果進(jìn)行積累得到對目標(biāo)的檢測結(jié)果soi(t)。

2 性能分析

為了驗(yàn)證信息分段波形的信息傳輸與目標(biāo)探測性能,本節(jié)從以下兩個(gè)方面進(jìn)行分析:

1)對信息分段波形與二進(jìn)制振幅鍵控系統(tǒng)(2ASK)、二進(jìn)制頻移鍵控系統(tǒng)(2FSK)的誤碼率進(jìn)行了對比分析,以此驗(yàn)證信息分段波形的信息傳輸性能。

2)對信息分段波形與傳統(tǒng)匹配濾波(加相參積累)的峰值信噪比進(jìn)行了對比分析,以此驗(yàn)證信息分段波形的目標(biāo)探測性能。

2.1 誤碼率

圖3為隨機(jī)生成的信息段I1=[0111101001101],經(jīng)過16 次積累后的信息提取結(jié)果rI10,rI11,在圖中分別用藍(lán)色和橙色的線表示,對應(yīng)信息位“0”和“1”。

圖3 信息段I1提取結(jié)果圖

可以得到:通過對比輸出rI10和rI11在峰值處的大小就可以提取出嵌入的信息,判決準(zhǔn)則為

這種判決準(zhǔn)則與常用的二進(jìn)制數(shù)字調(diào)制系統(tǒng)—2ASK、2FSK 的判決準(zhǔn)則(通過與判決門限/抽樣判決器兩路信號比較來解調(diào)信息)類似,故可作為誤碼率性能的比較對象。

2ASK、2FSK 系統(tǒng)常用的解調(diào)方法有非相干解調(diào)和相干解調(diào)[11],誤碼率分別如式(8)、式(9)所示:

式中,Eb為每比特信號能量,N0為噪聲功率譜密度,Eb/N0為比特信噪比。通過與2FSK 信號誤碼率的對比來分析信息分段波形的誤碼率性能。

圖4為信息分段波形與2FSK、2ASK 相干解調(diào)、非相干解調(diào)(在單次回波的情況下)的誤碼率經(jīng)過106次蒙特卡洛仿真的結(jié)果。可以得到在單次回波下,信息分段波形的誤碼率性能優(yōu)于2ASK 非相干解調(diào),在信噪比大于5 dB 時(shí)優(yōu)于2ASK 相干解調(diào);同時(shí)也可以得到信息分段波形的誤碼率曲線與2FSK 非相干解調(diào)基本吻合,誤碼率性能與2FSK 非相干解調(diào)相近,劣于2FSK 相干解調(diào)。

圖4 單次回波誤碼率隨信噪比的變化

圖5為信息分段波形與2FSK、2ASK 相干解調(diào)、非相干解調(diào)在積累的情況下,經(jīng)過106次蒙特卡洛仿真的結(jié)果,從中可以看出:隨著脈沖積累次數(shù)M的增加,信息分段波形及2FSK、2ASK 信號的誤碼率降低,誤碼率性能均提高。

圖5 積累后誤碼率隨信噪比的變化

經(jīng)過4次積累后,信息分段波形的誤碼率性能在信噪比大于0 dB 時(shí)優(yōu)于2FSK 非相干解調(diào),在信噪比大于2 dB 時(shí)優(yōu)于2ASK 相干解調(diào);16 次積累后,信息分段波形的誤碼率性能優(yōu)于2FSK、2ASK非相干解調(diào),在信噪比大于-3 dB 時(shí)優(yōu)于2ASK 相干解調(diào);32 次積累后,信息分段波形的誤碼率性能優(yōu)于2FSK 非相干解調(diào)和2ASK 相干解調(diào)、非相干解調(diào),同時(shí)趨近于2FSK相干解調(diào)。

通過前面的分析我們可以得到:在不積累的情況下,信息分段波形的誤碼率性能與2FSK 非相干解調(diào)相近,優(yōu)于2ASK 非相干解調(diào);在積累的情況下,信息分段波形的誤碼率性能隨著積累次數(shù)的增加而提高,經(jīng)過32次積累后,信息分段波形的誤碼率性能優(yōu)于2FSK 非相干解調(diào)和2ASK 相干解調(diào)、非相干解調(diào),故信息分段波形及其信號處理方法能夠滿足信息傳輸需求。

2.2 目標(biāo)探測性能

以信息段I1的單次回波為例,對信息分段波形的目標(biāo)探測性能進(jìn)行分析,信息段I1的單次回波模型如式(10)所示:

式中,I1n為信息段I1的第n位信息序列,t0為回波時(shí)延,w(t)表示功率譜密度為N0/2 的高斯白噪聲。將r11(t)分別與負(fù)線性調(diào)頻信號s0(t)、正線性調(diào)頻信號s1(t)進(jìn)行第一次匹配濾波,根據(jù)第1節(jié)中介紹的信息分段波形處理方法可以得到單次信息提取結(jié)果:

將式(10)代入式(11)可得

r110(t)和r111(t)分別代表信息序列中的“0”和“1”,信息序列提取結(jié)果的正負(fù)由外層相位ejφn決定,由式(12)可以得到,當(dāng)I1n為0時(shí),

當(dāng)I1n為1時(shí),

故單次信息提取結(jié)果的峰值信噪比與脈寬為Ts的線性調(diào)頻信號進(jìn)行匹配濾波的峰值信噪比相同,即輸出峰值信噪比為2EsN0,Es為信息分段波形子脈寬對應(yīng)的能量。

r110(t)和r111(t)相加后與由13 位巴克碼組成的梳齒濾波器h(t)進(jìn)行第二次匹配濾波,得到目標(biāo)檢測結(jié)果:

梳齒濾波器h(t)自相關(guān)函數(shù)的主峰峰值為旁瓣峰值的13 倍,且由式(13)和式(14)可知,r110(t)和r111(t)相加后,噪聲功率譜密度翻倍,對應(yīng)信息處輸出峰值信噪比為EsN0,故經(jīng)過第二次匹配濾波后目標(biāo)檢測輸出的峰值信噪比為13EsN0。

信號脈寬與子脈寬的關(guān)系為T=13Ts,與信號相同脈寬的線性調(diào)頻信號輸出峰值信噪比為2EN0,其中E=13Es,故與傳統(tǒng)的匹配濾波相比,目標(biāo)檢測的輸出峰值信噪比損失為1 2(約3 dB)。

由式(15)得到,目標(biāo)檢測輸出的距離分辨率主要由分別對應(yīng)信息位“1”“0”的正、負(fù)線性調(diào)頻信號的自相關(guān)函數(shù)R0(t)、R1(t)決定,故信息分段波形的距離分辨率為c(2B),其中B為信號的帶寬。

接下來對多脈沖情況下進(jìn)行分析。在多脈沖情況下本文方法使用了相干積累處理,理想情況下得到的峰值信噪比為單個(gè)脈沖處理時(shí)的M倍。下面在輸入信噪比為-5 dB 的條件下,對不同積累次數(shù)時(shí)信息分段波形和傳統(tǒng)匹配濾波的輸出峰值信噪比進(jìn)行了103次蒙特卡洛仿真,得到的輸出峰值信噪比隨積累次數(shù)的變化曲線如圖6所示。

圖6 峰值信噪比分析

由圖6可以得到,單個(gè)脈沖情況下,信息分段波形相比于傳統(tǒng)匹配濾波峰值信噪比損失約為3 dB;隨著積累次數(shù)的增加,信息分段波形的峰值信噪比損失略小于3 dB,這與前面的理論分析是一致的。

如果場景中存在多個(gè)目標(biāo),只要信息分段波形能夠在距離維上將各目標(biāo)準(zhǔn)確分辨出來(若距離分辨率不滿足,可以通過增加信號帶寬B,來提高距離分辨率),則本文提出的信息分段波形也能實(shí)現(xiàn)信息傳遞和目標(biāo)探測的功能。具體思路如下:通過記錄檢測得到各目標(biāo)對應(yīng)的距離單元,依次保留各目標(biāo)第一次匹配濾波(信息提取)、第二次匹配濾波(目標(biāo)檢測)后對應(yīng)距離單元的信號,分別對對應(yīng)不同目標(biāo)的多個(gè)回波信號進(jìn)行積累處理,從而得到信息提取結(jié)果和對各目標(biāo)的探測結(jié)果。對多目標(biāo)場景下信息分段波形的性能分析將作為后續(xù)研究內(nèi)容。

3 仿真實(shí)驗(yàn)

為了驗(yàn)證信息分段波形在中高重頻場景下的適用性,本節(jié)主要從以下兩個(gè)方面進(jìn)行仿真分析:

1)仿真分析中重頻場景下積累處理對分段波形信息提取及目標(biāo)檢測的性能改善情況。

2)為了驗(yàn)證信息分段波形在高重頻場景下的適用性,在與中重頻場景相等的回波積累時(shí)間這一條件的約束下,對比分析了積累處理對信息分段波形的信息提取與目標(biāo)檢測性能的改善情況。

3.1 仿真參數(shù)設(shè)置

假設(shè)目標(biāo)最遠(yuǎn)距離為150 km,最大速度為300 m/s,雷達(dá)載頻為10 GHz,對應(yīng)重頻范圍[12]如表1所示。

表1 脈沖重復(fù)頻率分類

首先對中、高重頻場景的仿真參數(shù)進(jìn)行設(shè)置,信息分段波形的基本參數(shù)設(shè)置如表2所示。

表2 基本參數(shù)設(shè)置

在表3中對信息分段波形的重頻參數(shù)進(jìn)行設(shè)置,分別設(shè)置了對應(yīng)中、高重頻場景的脈寬、占空比、目標(biāo)運(yùn)動(dòng)的速度等參數(shù)。

表3 重頻參數(shù)設(shè)置

3.2 信息分段波形在中重頻場景下的性能分析

式(16)為隨機(jī)生成的52 位“01”傳輸信息序列,13位一組(一個(gè)脈沖),分為4組進(jìn)行傳輸:

圖7為發(fā)射信號的仿真圖,為提高信息提取和目標(biāo)檢測的準(zhǔn)確度,在發(fā)射時(shí)將嵌入每個(gè)信息段的脈沖重復(fù)發(fā)射M次。

圖7 信息分段波形發(fā)射信號仿真圖

為信息提取的結(jié)果如圖8所示,可以看到4 次回波積累的信息提取結(jié)果淹沒在噪聲和旁瓣中,信息提取出錯(cuò)。經(jīng)過32 次的積累,信息對應(yīng)位置幅值增大,能準(zhǔn)確提取出各段嵌入的信息序列。

圖8 中重頻場景下信息段一提取結(jié)果

在提取出信息序列后,進(jìn)行對目標(biāo)的檢測。以信息段I1為例,將第一次匹配濾波輸出r110(t)和r111(t)相加,再與由13位巴克碼組成的梳齒濾波器h(t)進(jìn)行匹配,得到單次目標(biāo)檢測結(jié)果。將單次目標(biāo)檢測結(jié)果進(jìn)行32次積累,得到目標(biāo)檢測結(jié)果,如圖9所示。

圖9 目標(biāo)檢測結(jié)果

從圖9可以得到:經(jīng)過32 次脈沖積累之后,能夠準(zhǔn)確地檢測出目標(biāo)。同時(shí)速度檢測為99.16 m/s,與目標(biāo)的真實(shí)速度誤差在可接受范圍內(nèi)。

在積累次數(shù)為4,8,16,32,64次的情況下,對目標(biāo)檢測進(jìn)行103次蒙特卡洛仿真,結(jié)果如圖10所示。

圖10 中重頻場景下目標(biāo)檢測結(jié)果隨積累次數(shù)的變化

結(jié)合圖10的仿真結(jié)果可以將信息分段波形在中重頻場景中的性能總結(jié)如下:

1)積累次數(shù)增加后,目標(biāo)檢測峰值提高,有利于提高目標(biāo)檢測的準(zhǔn)確度。

2)隨著積累次數(shù)的增加,測速誤差減小。故通過增加積累次數(shù)M,能夠提高目標(biāo)速度檢測結(jié)果準(zhǔn)確度,但也帶來了時(shí)間成本的增加。

3.3 信息分段波形在高重頻場景下的性能分析

高重頻場景參數(shù)設(shè)置如3.1 節(jié)所示,發(fā)射的信息序列與中重頻場景下一致,在積累時(shí)間相同的條件下,高重頻場景的積累次數(shù)與中重頻積累次數(shù)的關(guān)系如表4所示。

表4 積累次數(shù)設(shè)置

經(jīng)過20 次、160 次積累后信息提取結(jié)果如圖11所示。

圖11 高重頻場景下信息段一提取結(jié)果

從圖11可以看到,經(jīng)過20 次積累后的信息提取結(jié)果出錯(cuò),不能準(zhǔn)確提取出信息序列。經(jīng)過160次的積累,信息對應(yīng)位置峰值增大,能夠準(zhǔn)確地提取出嵌入的信息序列。

圖12為高重頻場景下,積累160 次后的目標(biāo)檢測結(jié)果。

圖12 目標(biāo)檢測結(jié)果

從圖12可以看到,經(jīng)過160 次回波積累后能夠準(zhǔn)確地檢測出目標(biāo),目標(biāo)速度檢測結(jié)果為279.45 m/s,與真實(shí)速度的誤差在可接受范圍內(nèi)。

對表4中高重頻場景下的各積累次數(shù)進(jìn)行103次目標(biāo)檢測的蒙特卡洛仿真,結(jié)果如圖13所示,從圖中可以看到,在積累時(shí)間相同的條件下,高重頻場景積累帶來的增益與中重頻場景下得到的積累增益基本相同,信息分段波形在高重頻場景中的性能與中重頻場景相似,即隨著積累次數(shù)的增加,目標(biāo)速度的檢測精度提高,與真實(shí)速度的偏差越來越小,同時(shí)目標(biāo)檢測的峰值提高。

圖13 高重頻場景下目標(biāo)檢測結(jié)果隨積累次數(shù)的變化

4 結(jié)束語

為了在中、高重頻場景下確保通信化雷達(dá)的信息提取與目標(biāo)檢測性能,本文在文獻(xiàn)[1]的基礎(chǔ)上對通信化雷達(dá)的波形設(shè)計(jì)方法進(jìn)行了優(yōu)化,提出了信息分段波形的概念及其信號處理方法,該波形能夠在無通信鏈路的情況下,實(shí)現(xiàn)信息傳遞和目標(biāo)探測。通過將信息分段波形的誤碼率性能與二進(jìn)制數(shù)字調(diào)制系統(tǒng)2ASK、2FSK 系統(tǒng)進(jìn)行對比,得到了信息分段波形的誤碼率性能與2FSK 非相干解調(diào)相近,優(yōu)于2ASK 非相干解調(diào)的結(jié)論,同時(shí)經(jīng)過積累處理后,信息分段波形的誤碼率性能進(jìn)一步提升,能夠滿足信息傳輸需求;通過將信息分段波形的目標(biāo)探測性能與傳統(tǒng)線性調(diào)頻波形進(jìn)行對比,得到了信息分段波形及其處理方法輸出信噪比損失略小于3 dB的結(jié)論。最后,通過仿真實(shí)驗(yàn)分別驗(yàn)證了信息分段波形在中、高重頻場景下的適用性。信息分段波形內(nèi)層采用正負(fù)線性調(diào)頻信號,外層采用13 位巴克碼的相位調(diào)制,具有恒包絡(luò)特性,有利于發(fā)射機(jī)功率放大器工作在線性區(qū),提高能量效率。本文所做工作為通信化雷達(dá)在機(jī)載平臺(tái)上應(yīng)用奠定了重要基礎(chǔ),具有理論與工程意義。

猜你喜歡
信號檢測信息
“不等式”檢測題
“一元一次不等式”檢測題
“一元一次不等式組”檢測題
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
基于FPGA的多功能信號發(fā)生器的設(shè)計(jì)
電子制作(2018年11期)2018-08-04 03:25:42
訂閱信息
中華手工(2017年2期)2017-06-06 23:00:31
小波變換在PCB缺陷檢測中的應(yīng)用
基于LabVIEW的力加載信號采集與PID控制
展會(huì)信息
主站蜘蛛池模板: 色综合久久88| 国产小视频a在线观看| 91原创视频在线| 久久永久视频| 在线视频亚洲欧美| 99视频精品在线观看| 在线观看国产小视频| 精品视频福利| 亚洲毛片一级带毛片基地| 67194成是人免费无码| 青青青草国产| 综合网天天| 欧美成人综合视频| 国产成人精品三级| 亚洲香蕉久久| 亚洲码在线中文在线观看| 青青青国产免费线在| 91在线免费公开视频| A级毛片无码久久精品免费| 亚洲最黄视频| 亚洲娇小与黑人巨大交| 国产成人精品优优av| 好紧太爽了视频免费无码| 国产97视频在线| 久久黄色视频影| 九九热这里只有国产精品| 国产免费久久精品44| 免费看一级毛片波多结衣| 国产美女在线观看| 欧美中文一区| 草草影院国产第一页| 一区二区在线视频免费观看| 91精品啪在线观看国产| 午夜国产精品视频| 少妇精品久久久一区二区三区| 亚洲精品国产综合99| 亚洲人成在线精品| 精品视频一区在线观看| 亚洲黄色视频在线观看一区| 亚洲91精品视频| 熟女视频91| 亚洲美女视频一区| 国产综合无码一区二区色蜜蜜| 在线视频97| 亚洲精品成人福利在线电影| 四虎成人免费毛片| 亚洲成aⅴ人在线观看| 一级一级一片免费| 国产黄在线观看| 一本久道久久综合多人| 18禁黄无遮挡网站| julia中文字幕久久亚洲| 国产爽妇精品| 亚洲人成网7777777国产| 亚洲婷婷在线视频| 国产成人三级| 免费av一区二区三区在线| 亚洲最大情网站在线观看 | 成人午夜视频免费看欧美| 欧美亚洲国产一区| 成人免费视频一区| 最新无码专区超级碰碰碰| 久久免费观看视频| 国产精品永久在线| 永久免费AⅤ无码网站在线观看| 亚洲黄色激情网站| 久久人人妻人人爽人人卡片av| 亚国产欧美在线人成| 久久不卡国产精品无码| 日韩国产无码一区| 久久大香伊蕉在人线观看热2| 91福利一区二区三区| av天堂最新版在线| 91www在线观看| 精品人妻一区无码视频| 成人国产小视频| 国产亚洲精| 日韩在线播放欧美字幕| 日韩高清无码免费| 欧美综合激情| 国产福利拍拍拍| 成人午夜网址|