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

        ?

        不完整復(fù)制信號在雷達(dá)中的檢測性能

        2017-07-12 17:47:34楊立明蔚建利
        探測與控制學(xué)報(bào) 2017年3期
        關(guān)鍵詞:濾波器雷達(dá)概率

        楊立明,韋 高,蔚建利

        (1.西北工業(yè)大學(xué)電子信息學(xué)院,陜西 西安 710072;2.西北工業(yè)大學(xué)航天學(xué)院,陜西 西安 710072)

        不完整復(fù)制信號在雷達(dá)中的檢測性能

        楊立明1,韋 高1,蔚建利2

        (1.西北工業(yè)大學(xué)電子信息學(xué)院,陜西 西安 710072;2.西北工業(yè)大學(xué)航天學(xué)院,陜西 西安 710072)

        針對缺少不完整復(fù)制信號在雷達(dá)接收機(jī)上的檢測性能的研究,提出了對不完整復(fù)制信號在雷達(dá)中的檢測性能進(jìn)行分析的方法。該方法對不完整復(fù)制的常規(guī)信號和LFM信號進(jìn)行了建模,分析了上述信號在雷達(dá)中的檢測性能,為客觀評估該類信號在雷達(dá)中產(chǎn)生的干擾效果提供了理論依據(jù)。仿真驗(yàn)證表明,兩種不完整復(fù)制的信號都可以在雷達(dá)中被檢測,在適當(dāng)?shù)膹?fù)制長度和信噪比條件下,不完整復(fù)制信號可以獲得和原信號十分接近的發(fā)現(xiàn)概率。

        不完整復(fù)制信號;LFM信號;檢測性能;干擾效果

        0 引言

        現(xiàn)代電子戰(zhàn)技術(shù)中,通過復(fù)制信號的一部分(為方便,本文簡稱為“不完整復(fù)制信號”)進(jìn)行轉(zhuǎn)發(fā),轉(zhuǎn)發(fā)信號被雷達(dá)接收時(shí)能夠獲得部分處理增益,進(jìn)而產(chǎn)生影響雷達(dá)正常接收的假目標(biāo)信號[1]。假目標(biāo)信號用于干擾時(shí)使用靈活,效果明顯,對雷達(dá)的正常接收產(chǎn)生了明顯的影響。

        國外針對該類信號主要開展了形成干擾的機(jī)理研究[1],但未見針對該類信號的檢測性能研究的文獻(xiàn)。在國內(nèi),對此類信號形成的假目標(biāo)干擾的效果比如功率損失、信噪比損失進(jìn)行了研究[2],對與此相關(guān)的脈內(nèi)疊加信號的性能進(jìn)行了分析[3],還提出了部分拓展的干擾方法[4-5],另外,在抗干擾方面也進(jìn)行了一定的研究[6]。

        但是,目前的研究尚未涉及到不完整復(fù)制信號在雷達(dá)接收機(jī)上的檢測性能,本文主要針對不完整復(fù)制信號的特點(diǎn),分析了該信號通過雷達(dá)匹配濾波后的變化特征,研究了該信號在雷達(dá)中的檢測性能,從而為客觀評估該類信號產(chǎn)生的干擾效果提供了理論依據(jù)。

        1 不完整復(fù)制雷達(dá)信號模型

        1.1 常規(guī)信號模型

        為了分析不完整復(fù)制常規(guī)信號在雷達(dá)中的檢測性能,需要對不完整復(fù)制常規(guī)信號建模,并分析其在通過原常規(guī)信號(這里指完整的常規(guī)信號)匹配濾波后的檢測性能。

        (1)

        (2)

        式中,ΔT為復(fù)制時(shí)間長度,通常ΔT≤T。

        1.2LFM信號模型

        雷達(dá)的發(fā)射信號一般可以表示為:

        v(t)=u(t)ej2πf0t

        (3)

        現(xiàn)代雷達(dá)為了實(shí)現(xiàn)探測距離和分辨率的同時(shí)實(shí)現(xiàn),大多采用了脈沖壓縮技術(shù)。脈沖壓縮技術(shù)的實(shí)現(xiàn)可以分為線性調(diào)頻和相位編碼,下面的分析以線性調(diào)頻脈沖壓縮技術(shù)為例。采用線性調(diào)頻脈沖壓縮技術(shù)的雷達(dá)發(fā)射信號可以寫為[7]:

        (4)

        (5)

        式中,f1=f0+μ(τ/2+τs)為不完整復(fù)制LFM信號的中心頻率。

        2 不完整復(fù)制雷達(dá)信號檢測性能

        2.1 不完整復(fù)制常規(guī)信號檢測性能分析

        不完整復(fù)制常規(guī)信號被接收后,首先要經(jīng)過原常規(guī)信號的匹配濾波器處理,原常規(guī)信號的匹配濾波器脈沖響應(yīng)滿足[7]

        (6)

        t0是濾波器的延遲時(shí)間。則不完整復(fù)制常規(guī)信號通過原信號的匹配濾波器后的輸出為:

        (7)

        令t-t0=t′,則

        (8)

        當(dāng)0?t′?T時(shí),

        (9)

        當(dāng)t′=0,即t=t0時(shí),式(9)輸出值最大,

        (10)

        式(9)積分部分是原常規(guī)信號與其共軛信號的乘積,代表原信號的模平方。由于原常規(guī)信號已知,所以其模平方可以認(rèn)為是已知的常量,設(shè)為M。

        則式(9)可以簡化為:

        (11)

        經(jīng)過同樣的計(jì)算過程可知,-T?t′?0,式(8)的輸出與式(10)相同。

        則在t=t0時(shí)不完整復(fù)制常規(guī)信號經(jīng)過原常規(guī)信號的匹配濾波器后輸出的最大信噪比可以表示為:

        (12)

        在ΔT=T時(shí),相當(dāng)于原常規(guī)信號經(jīng)過匹配濾波器,則經(jīng)過匹配濾波后的輸出信噪比可以表示為:

        (13)

        回波信號可以看成是幅度和相位均隨機(jī)變化的,在此情況下,發(fā)現(xiàn)概率和虛警概率存在如下關(guān)系[8]:

        (14)

        將式(12)、式(13)分別代入式(14)即可得到雷達(dá)對不完整復(fù)制常規(guī)信號和原常規(guī)信號的發(fā)現(xiàn)概率分別為:

        (15)

        (16)

        2.2 不完整復(fù)制LFM信號檢測性能分析

        2.2.1 不完整復(fù)制LFM信號的匹配濾波特性

        為分析方便,把復(fù)制LFM信號的矩形脈沖的中心作為時(shí)間原點(diǎn),則影響不完整復(fù)制LFM信號匹配濾波特性的因素僅有復(fù)制長度。此時(shí),f1=f0,進(jìn)而可以推算出τs=-τ/2,式(5)可以改寫為:

        (17)

        不完整復(fù)制LFM信號之所以能有效干擾雷達(dá),是因?yàn)樵撔盘柦?jīng)過原發(fā)射信號的匹配濾波器后,可以獲得與原發(fā)射信號相近的處理增益。根據(jù)匹配濾波理論,信號S(t)的匹配濾波器頻率特性可以寫為[7]:

        (18)

        f表示頻率。雷達(dá)發(fā)射信號S(t)通過該濾波器后的輸出為[4]:

        (19)

        而式(17)所示的不完整復(fù)制LFM信號經(jīng)過同樣的匹配濾波器后的輸出為:

        (20)

        式中,B′=μτ為不完整復(fù)制LFM信號的帶寬。由式(20)可以看出,不完整復(fù)制LFM信號經(jīng)過雷達(dá)發(fā)射信號的匹配濾波器后,輸出的信號包絡(luò)仍近似為sinc函數(shù),但是幅度和主瓣寬度都發(fā)生了變化。

        幅度變化:

        (21)

        主瓣寬度變化:

        (22)

        可以看出,不完整復(fù)制LFM的匹配濾波輸出幅度降為雷達(dá)發(fā)射信號的匹配濾波器輸出幅度的1/K0,而主瓣寬度則展寬到原來的K0倍。

        取雷達(dá)發(fā)射LFM信號的帶寬為30MHz,時(shí)寬為10μs,不完整復(fù)制LFM長度為2.5μs。這兩個(gè)信號經(jīng)過同一個(gè)匹配濾波器后的輸出波形如圖1所示。

        由圖1中可以看出,不完整復(fù)制LFM信號在原匹配濾波器中會(huì)有失配損失,如果考慮加窗、信噪比等損失[2],不完整復(fù)制LFM信號的損失可達(dá)33dB。但是,干擾信號比目標(biāo)回波信號具有傳播路徑的優(yōu)勢,而且干擾機(jī)功率一般能滿足在雷達(dá)接收機(jī)處的功率比回波大12~13dB。因此,在不完整復(fù)制LFM干擾信號作用下,大部分干擾信號可以通過匹配濾波器而進(jìn)入雷達(dá)接收機(jī),完全可以用于對雷達(dá)實(shí)施干擾,造成對雷達(dá)檢測的影響。

        2.2.2 不完整復(fù)制LFM信號對雷達(dá)檢測性能的影響

        不完整復(fù)制LFM信號被雷達(dá)接收時(shí),該信號和回波信號同時(shí)被檢測器從噪聲中檢測。不完整復(fù)制LFM信號的存在會(huì)影響到回波信號的檢測,這是因?yàn)椴煌暾麖?fù)制LFM信號能夠獲得雷達(dá)接收機(jī)的處理增益,在雷達(dá)顯示器上形成密集假目標(biāo),使得操作員對目標(biāo)的識(shí)別判斷困難。首先分析只存在一個(gè)不完整復(fù)制LFM信號的情況,此時(shí)雷達(dá)接收信號模型可以寫為:

        0≤t

        (23)

        (24)

        (25)

        a1為干擾信號的復(fù)幅度,也是隨機(jī)變量,設(shè)其相位分布與a0相同。由式(21)可知,雷達(dá)信號和不完整復(fù)制LFM信號通過匹配濾波器后的幅度比K0=T/τ,則不完整復(fù)制LFM信號的幅度分布密度函數(shù)可以寫為:

        (26)

        Neyman-Pearson準(zhǔn)則下,對目標(biāo)回波的最優(yōu)檢測器為似然比檢測器,其檢測統(tǒng)計(jì)過程為[6]

        (27)

        在該檢測過程中,經(jīng)過包線檢波后的采樣時(shí)刻t0是指匹配濾波器輸出信號噪聲比達(dá)到最大值的時(shí)刻,對一個(gè)物理可實(shí)現(xiàn)的系統(tǒng)而言,該時(shí)刻是指在輸入信號全部結(jié)束后[7]。

        設(shè)LFM信號調(diào)頻帶寬為100 MHz,信號時(shí)長100 μs。干擾信號是原信號的20 μs重復(fù)復(fù)制形成的,干擾信號時(shí)延為1 μs。干信比為20 dB時(shí),LFM信號和不完整復(fù)制信號經(jīng)過匹配濾波和包絡(luò)檢波后的輸出波形如圖2、圖3所示。

        該檢測過程的門限由下式確定[8]

        (28)

        式中Pf為虛警概率,Pf設(shè)為固定值。此時(shí)的發(fā)現(xiàn)概率可以根據(jù)下式得到[8]

        (29)

        (30)

        當(dāng)信號幅度服從式(25)所示的分布密度函數(shù)的瑞利分布時(shí),回波信號的統(tǒng)計(jì)平均能量為:

        (31)

        同理,可以計(jì)算出由式(5)給出的不完整復(fù)制LFM信號的能量為:

        (32)

        不完整復(fù)制LFM信號被雷達(dá)接收時(shí),其幅度經(jīng)歷了和回波信號相同的傳播過程,可以認(rèn)為其幅度也服從瑞利分布,則可以計(jì)算出不完整復(fù)制LFM信號的統(tǒng)計(jì)平均能量為:

        (33)

        (34)

        將式(21)、式(32)、式(33)代入式(34)可得

        (35)

        為了計(jì)算方便,假設(shè)噪聲功率譜密度N0=1,回波信號脈沖寬度100 μs。計(jì)算不同復(fù)制長度時(shí)的發(fā)現(xiàn)概率隨比值k的變化情況如圖4所示。

        由圖4中可以看出,復(fù)制長度越短,雷達(dá)檢測器對干擾信號的發(fā)現(xiàn)概率就越小,反應(yīng)在圖中曲線上就是k值的增大。反之,復(fù)制長度越長,k值越小,干擾信號被檢測到的概率就越大。在相同復(fù)制長度時(shí),虛警概率越低,k值越大。這也符合實(shí)際情況,在雷達(dá)檢測中,Neyman-Pearson準(zhǔn)則以虛警概率來設(shè)置檢測門限,虛警概率低,相當(dāng)于檢測門限高。在同樣復(fù)制長度時(shí),干擾信號獲取匹配濾波的處理增益相同,門限抬高后,對應(yīng)的檢測概率自然降低。

        雷達(dá)發(fā)射的LFM信號通過上述過程后,可以得到匹配濾波的全部處理增益,在經(jīng)過包線檢波,采樣和門限判定后,可以被有效檢測。

        由于回波信號受目標(biāo)反射特性和電波傳播的影響,回波信號可以看成是幅度和相位均隨機(jī)變化的,在此情況下,發(fā)現(xiàn)概率和虛警概率的關(guān)系如式(29)所示。

        (36)

        (37)

        3 仿真驗(yàn)證

        3.1 部分復(fù)制常規(guī)信號的檢測性能仿真

        設(shè)原常規(guī)信號時(shí)長為20μs,虛警概率為10-6,噪聲功率譜密度變化范圍為0.02~2μW/Hz。在信號時(shí)長和復(fù)制長度一定時(shí),對應(yīng)的信號功率和復(fù)制信號的功率也就確定,由式(12)和式(13)可知,噪聲功率譜密度的變化對應(yīng)的即為信噪比的變化。不完整復(fù)制信號的長度以占原信號時(shí)長的比例k=ΔT/T表示,取值范圍為0.1~1。

        計(jì)算雷達(dá)對不同復(fù)制長度的不完整復(fù)制常規(guī)信號的發(fā)現(xiàn)概率在不同噪聲功率譜密度下的變化趨勢如圖5所示。

        從圖5可以看出,不完整復(fù)制常規(guī)信號在雷達(dá)中的發(fā)現(xiàn)概率與復(fù)制長度密切相關(guān)。復(fù)制長度為原信號時(shí)長0.1倍的信號要達(dá)到與原常規(guī)信號相同的發(fā)現(xiàn)概率(以Pd=0.3為例),噪聲功率譜密度相差可達(dá)2倍。換言之,在達(dá)到0.3的發(fā)現(xiàn)概率時(shí),原常規(guī)信號可以在1.8 μW/Hz噪聲功率譜密度水平下被發(fā)現(xiàn),而復(fù)制長度為原常規(guī)信號時(shí)長0.1倍的信號則只能在0.9 μW/Hz的噪聲功率譜密度水平下被發(fā)現(xiàn)。隨著噪聲功率譜密度的減小,雷達(dá)對不完整復(fù)制常規(guī)信號的發(fā)現(xiàn)概率逐漸與對原常規(guī)信號的發(fā)現(xiàn)概率接近。這與實(shí)際情況也吻合,因?yàn)樵肼暪β首V密度減小,無論對復(fù)制信號還是原信號,匹配濾波器的輸出信噪比都會(huì)增大,信噪比增大到一定程度時(shí),由式(14)可知,雷達(dá)對復(fù)制常規(guī)信號和原常規(guī)信號的發(fā)現(xiàn)概率將趨于相當(dāng)。

        3.2 部分復(fù)制LFM信號的檢測性能仿真

        設(shè)線性調(diào)頻信號中心頻率1 GHz,帶寬30 MHz,脈寬20 μs;噪聲功率譜0.02~2 μW/Hz?;夭ㄐ盘柵c不完整復(fù)制LFM信號的幅度均服從均值為1,方差為1的瑞利分布。計(jì)算雷達(dá)對不同復(fù)制長度的不完整復(fù)制LFM信號的發(fā)現(xiàn)概率如圖6所示。其中,復(fù)制長度以占信號時(shí)長的比例表示,比例系數(shù)記為k,取值范圍為0.1~1。

        從圖6可以看出,與常規(guī)信號類似,不完整復(fù)制LFM信號在雷達(dá)中的發(fā)現(xiàn)概率與截取長度密切相關(guān)。復(fù)制長度為原信號時(shí)長0.1倍的信號要達(dá)到與原信號相同的發(fā)現(xiàn)概率(設(shè)Pd=0.55),噪聲功率譜密度相差可達(dá)2倍。換言之,在達(dá)到0.55的發(fā)現(xiàn)概率時(shí),原信號可以在0.8 μW/Hz噪聲功率譜密度水平下被發(fā)現(xiàn),而截取長度為原信號時(shí)長0.1倍的信號則只能在0.4 μW/Hz的噪聲功率譜密度水平下被發(fā)現(xiàn)。隨著噪聲功率譜密度的減小,雷達(dá)對不完整復(fù)制信號的發(fā)現(xiàn)概率逐漸與對原信號的發(fā)現(xiàn)概率接近。總體上,不完整復(fù)制的LFM信號在雷達(dá)中的檢測性能變化趨勢與第一節(jié)分析的情況基本相同,證明了分析結(jié)果的合理性。

        (38)

        (39)

        式(39)說明,不完整復(fù)制LFM信號經(jīng)過不斷轉(zhuǎn)發(fā)后,輸出為相互間隔為τ/2的近似目標(biāo)回波信號,該信號的幅度和主瓣寬度相對于原發(fā)射信號通過匹配濾波器器的比值仍可以用式(21)、式(22)表示。由于復(fù)制信號的長度遠(yuǎn)小于發(fā)射信號,通過重復(fù)發(fā)射截取信號,可以在雷達(dá)顯示器上形成比直接轉(zhuǎn)發(fā)整個(gè)LFM信號密集的多的近似回波信號,造成雷達(dá)檢測目標(biāo)困難。

        由式(32)可知,在存在部分復(fù)制LFM信號時(shí),其被雷達(dá)檢測到的概率與原信號的發(fā)現(xiàn)概率呈一定的比例關(guān)系。這主要是因?yàn)椴煌暾麖?fù)制LFM信號的截取長度限制了匹配濾波器后的輸出,在實(shí)際使用中,有意發(fā)射的不完整復(fù)制信號相對于目標(biāo)回波的路徑少一倍,帶來的功率損失以1/R2衰減,R為不完整復(fù)制信號發(fā)射機(jī)與離雷達(dá)的距離。而回波信號功率按原信號的1/R4衰減,回波信號比不完整復(fù)制信號的衰減要大得多。因此,在雷達(dá)接收的信號中,不完整復(fù)制信號要比回波信號強(qiáng)很多,相應(yīng)的被發(fā)現(xiàn)的概率也大大增加。設(shè)雷達(dá)檢測器對不完整復(fù)制信號和回波信號的發(fā)現(xiàn)概率相同,回波信號的發(fā)現(xiàn)概率降為原來的1/N,則雷達(dá)檢測器對目標(biāo)的檢測性能將急劇惡化。設(shè)存在5個(gè)不完整復(fù)制信號,虛警概率為10-11,雷達(dá)對某個(gè)目標(biāo)的發(fā)現(xiàn)概率曲線和存在干擾時(shí)發(fā)現(xiàn)概率曲線的對比如圖7所示,可見雷達(dá)對目標(biāo)的發(fā)現(xiàn)概率下降十分嚴(yán)重,幾乎不能有效發(fā)現(xiàn)目標(biāo),實(shí)際使用中,不完整復(fù)制信號密集程度更高,雷達(dá)的對目標(biāo)的檢測性能下降嚴(yán)重,必須采取有效措施。

        4 結(jié)論

        本文提出了不完整復(fù)制信號在雷達(dá)中的檢測性能的分析方法。該方法建立了不完整復(fù)制常規(guī)信號和LFM信號的模型,分別分析了兩種信號在雷達(dá)中的檢測性能。仿真驗(yàn)證表明,兩種不完整復(fù)制的信號都可以在雷達(dá)中被檢測,在適當(dāng)?shù)膹?fù)制長度和信噪比條件下,不完整復(fù)制信號可以獲得和原信號十分接近的發(fā)現(xiàn)概率。在多次進(jìn)行復(fù)制后,可以形成可觀的干擾效果。不完整復(fù)制信號極大地影響了雷達(dá)對目標(biāo)的檢測,必須采取措施,及時(shí)消除不完整復(fù)制對雷達(dá)的影響,保證對目標(biāo)的檢測性能。

        [1]施萊赫D C.信息時(shí)代的電子戰(zhàn)[M].成都:信息產(chǎn)業(yè)部電子第二十九研究所,2000.

        [2]李宏,鄭光勇,楊英科,等.部分截取雷達(dá)信號多假目標(biāo)干擾性能分析[J].電子信息對抗技術(shù),2010(3):39-44.

        [3]舒玉貴,林選鋒,郭燕.脈內(nèi)疊加信號性能分析[J].電子信息對抗技術(shù),2010(2):64-68.

        [4]勞國超,葉偉,阮航.基于壓縮感知的LFM雷達(dá)轉(zhuǎn)發(fā)干擾方法[J].電子信息對抗技術(shù),2013(6):40-44.

        [5]潘小義,王偉,馮德軍,等.基于全脈沖分段轉(zhuǎn)發(fā)的LFM雷達(dá)干擾方法[J].國防科技大學(xué)學(xué)報(bào),2013(3):119-125.

        [6]粘朋雷,李國林,路翠華,等.一種線性調(diào)頻引信抗卷積調(diào)制干擾處理方法[J].電訊技術(shù),2013(8):1028-1032.

        [7]林茂庸,柯有安.雷達(dá)信號理論[M].北京:國防工業(yè)出版社,1984.

        [8]齊國清.信號檢測與估計(jì)-原理及應(yīng)用[M].北京:電子工業(yè)出版社,2010.

        Radar Detection Performance of Partly Copied Signal

        YANG Liming1, WEI Gao1, WEI Jianli2

        (1.School of Electronics and Information, Northwestern Polytechnic University, Xi’an 710072,China; 2.School of Astronautics,Northwestern Polytechnieal University, Xi’an 710072,China)

        The radar detection performance of the partly copied signal has not be studied at present. The partly copied signal detection performance analysis method was presented in this paper. The routine partly copied signal and LFM signal were modeled and the detection performance of radar were studied. The results of the simulation indicated that the two kind signals could be detected by radar. The finding probability of the partly copied signal was very close to the original signal at a likely copied length and SNR.

        partly copied signal; LFM signal; detection performance; jamming effective

        2016-12-17

        國家青年科學(xué)基金項(xiàng)目資助(61401361)

        楊立明(1976—),男,河北趙縣人,博士,研究方向:雷達(dá)信號處理。E-mail:ylmw@163.com。

        TN957.51

        A

        1008-1194(2017)03-0041-07

        猜你喜歡
        濾波器雷達(dá)概率
        基于無擾濾波器和AED-ADT的無擾切換控制
        有雷達(dá)
        大自然探索(2023年7期)2023-08-15 00:48:21
        第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
        第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
        概率與統(tǒng)計(jì)(一)
        概率與統(tǒng)計(jì)(二)
        從濾波器理解卷積
        電子制作(2019年11期)2019-07-04 00:34:38
        雷達(dá)
        開關(guān)電源EMI濾波器的應(yīng)用方法探討
        電子制作(2018年16期)2018-09-26 03:26:50
        基于空時(shí)二維隨機(jī)輻射場的彈載雷達(dá)前視成像
        乱码1乱码2美美哒| 成 人 色综合 综合网站| 日日摸夜夜添狠狠添欧美| 青青草综合在线观看视频| 毛片av在线尤物一区二区| 亚洲av色香蕉一区二区三区软件 | 色老头一区二区三区| 国产精品久久婷婷婷婷| 日本av一区二区三区四区| 日本一区二区三区亚洲| 国产精品爽爽v在线观看无码| 日韩乱码人妻无码中文字幕视频| 国产欧美va欧美va香蕉在线观| 精品一区二区三区在线观看l| 免费人妻精品区一区二区三| 日韩精品一区二区三区乱码| 午夜精品久久久久久久| 99国产精品无码| 456亚洲老头视频| 国产区高清在线一区二区三区| 青青草在线免费观看视频| 亚洲av高清一区二区三| 777国产偷窥盗摄精品品在线| 国产成人综合在线视频| 欧美激情国产亚州一区二区| 91麻豆精品久久久影院| 国产乱码一区二区三区精品| 久久久久无码国产精品一区| 国外精品视频在线观看免费| 国产精品久久久久尤物| 久久99热精品免费观看麻豆| 亚洲精品国产精品乱码视色| 2020年国产精品| 漂亮人妻被黑人久久精品| 喷潮出白浆视频在线观看| 亚洲国产一区二区网站| 亚洲av中文无码乱人伦在线视色| 亚洲永久精品ww47| 又色又爽又黄的视频网站| 亚洲国产一区一区毛片a| 亚洲av午夜精品无码专区|