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

    球載雷達多雜波自適應(yīng)抑制

    2019-12-23 09:01:24宿文濤劉潤華張小涵
    中國電子科學研究院學報 2019年10期
    關(guān)鍵詞:云雨雜波矢量

    宿文濤,汪 楓,劉潤華,張小涵

    (空軍預警學院, 湖北 武漢 430019)

    0 引 言

    近年來,球載雷達以預警能力強、低空性能好、留空時間長和效費比高等優(yōu)勢,越來越受到世界各國軍隊的青睞[1-2]。面對低空、超低空突防的威脅,球載雷達是有效的抵抗手段之一。而球載雷達檢測這些目標存在兩個主要難點:一是面臨的雜波環(huán)境與常規(guī)雷達有所區(qū)別[3],波束俯視導致雜波強度大,雷達平臺的移動使固定的地雜波頻譜發(fā)生移動和展寬,難以確切知道雜波的頻譜特征[4]。二是低空低速目標受到強大的地雜波的干擾[5],由于速度較慢,其特性在多普勒域與雜波非常接近,導致傳統(tǒng)的頻域濾波方法難以處理。

    對于低空低速目標的檢測,能否精準的去除雜波,保留低速目標是關(guān)鍵。這就要求濾波器的阻帶要與雜波的頻譜展寬相匹配,通帶要盡可能的平坦以保留目標信息。傳統(tǒng)的延遲線對消器是現(xiàn)役球載雷達采用的濾波方法,雖然延遲線對消器的改善因子隨對消脈沖數(shù)的提升而提升,但是其阻帶的寬度以及通帶的平坦度是不可控的,這就限制其對低空低速目標的檢測性能。

    針對上述問題,本文采用改進的特征矢量法設(shè)計濾波器來匹配雜波頻譜特征,使用相關(guān)函數(shù)法來估計頻譜特征,針對估計誤差過大的問題運用覆蓋法盡可能消除誤差的影響。通過計算機仿真,驗證了方法的有效性。

    1 雜波模型

    球載雷達接收的回波信號中可能包含地物、海浪以及云雨等多種雜波。雷達位于陸地上空時主要是地物和云雨雜波,位于海面上空時主要是海雜波[6]和云雨雜波。球載雷達面臨的雜波一般均可以用高斯模型表示[7],對于多個雜波來說其歸一化形式為

    (1)

    式中:fdk為第k個雜波的中心頻率;σck為第k個雜波的標準差,代表頻譜的展寬程度。

    (2)

    式中:vk為第k個雜波的相對水平速度;ψ為俯視角;λ為雷達波長。

    (3)

    例如波長為0.3 m,重復頻率為300 Hz,天線轉(zhuǎn)速為6 r/min,θ3=2°時,表1給出了一般情況下幾種雜波譜的典型譜寬。

    表1 幾種典型雜波的標準差

    根據(jù)維納濾波理論,雜波的自相關(guān)函數(shù)與功率譜密度互為傅里葉變換,因而雜波的自相關(guān)函數(shù)可以表示為

    Rc(i,j)=F-1[S(f)]=

    (4)

    令τij=ti-tj為相關(guān)時間,將(1)式帶入(4)式化簡得

    Rc(i,j)=

    (5)

    由式(5)可以看出,多個高斯譜的混合雜波,其自相關(guān)函數(shù)由對應(yīng)的多雜波分量之和構(gòu)成。由隨機過程相關(guān)理論可知,平穩(wěn)隨機過程的自相關(guān)函數(shù)就是它的協(xié)方差矩陣。因此,由Rc(i,j)構(gòu)成的多雜波協(xié)方差矩陣為

    (6)

    2 基于特征矢量法的實時處理技術(shù)

    濾波器的設(shè)計目標就是設(shè)計一組合適的濾波器系數(shù),使其有效地抑制雜波,并保證目標信號能無損失地通過。濾波器的輸出[9]為

    y(t)=wTx(t)=w(i)x(t-iTr)

    (7)

    式中:w=(w0,w1,…wN)T為權(quán)矢量,x(t)為信號矢量,Tr為球載雷達的脈沖重復周期。工程上常以改善因子來衡量權(quán)系數(shù)的性能優(yōu)劣。改善因子定義為I=(So/Co)/(Si/Ci),顯然,改善因子越大,其對雜波的抑制效果越好。

    用特征適量法設(shè)計的濾波器對多個雜波組成的混合雜波有較好的濾波效果[10]。特征矢量法是以平均改善因子最大為準則的雜波抑制方法。它屬于有限脈沖響應(yīng)(FIR)濾波器的一種,其系統(tǒng)函數(shù)為

    (8)

    式中:N為濾波器階數(shù)。

    已經(jīng)證明[11-12],要使改善因子達到最大,濾波器的最佳權(quán)矢量應(yīng)為輸入雜波的協(xié)方差矩陣Rc的最小特征值λmin對應(yīng)的特征向量。

    Rcwopt=λminwopt

    (9)

    此時,平均改善因子達到最大。

    (10)

    由式(9)可知,由特征矢量法計算濾波器系數(shù)w取決于雜波協(xié)方差矩陣Rc,由式(5)、式(6)可知Rc由fdk和σck所決定,所以最佳權(quán)系數(shù)w是fdk和σck的二元多維復函數(shù)矢量,記為w(fdk,σck)。

    如果已知濾波器的權(quán)系數(shù),則特征矢量法構(gòu)造的濾波器的頻率響應(yīng)函數(shù)與其它FIR濾波器類似,可以統(tǒng)一由式(11)表示:

    (11)

    雷達脈沖重復頻率300 Hz時特征矢量法濾波器的頻率響應(yīng)如下圖。為了便于闡述,不失一般性,本文以兩個雜波(地雜波和云雨雜波)為例闡述多雜波自適應(yīng)抑制。

    表2 仿真參數(shù)設(shè)置

    經(jīng)計算其中各權(quán)矢量數(shù)值如下:

    w3=[-0.60+0.06i,0.27+0.51i,0.61]

    w4=[0.38+0.14i,0.47-0.32i,

    0.34-0.46i,-0.41]

    w5=[-0.15-0.20i,-0.53-0.05i,

    -0.24+0.50i,0.37+0.39,0.26]

    wg=[-0.27+0.24i,-0.56+0.22i,

    0.27+0.54,0.36]

    由圖1、圖2可知特征矢量法設(shè)計的濾波器的凹口位置和寬度都隨雜波的數(shù)量和譜寬的變化而自適應(yīng)變化。這樣的特性比較適合用于球載雷達面臨的多個不同雜波的情形。但是特征矢量法構(gòu)造的濾波器處理多個譜寬不同的雜波時性能會有不同程度的惡化,濾波器階數(shù)較少時,特征矢量法求解的零點較少,分辨率不高導致兩個凹口寬度近似相等,濾波器階數(shù)較多時凹口較淺且寬度較大。因此要對該方法進行改進以適應(yīng)多個不同頻譜特征的雜波抑制。

    圖1 多凹口幅頻響應(yīng) (譜寬不同時)

    圖2 改進濾波器幅頻響應(yīng)

    針對該方法在抑制多個具有不同頻譜特性雜波時性能惡化的問題,可以采用設(shè)置虛擬雜波的方式來增加譜寬。由圖1可知特征矢量法處理譜寬相同的多個雜波性能較好,所以首先選擇譜寬較窄的地雜波的譜寬為參考作為其它所有雜波的譜寬值,然后在譜寬較寬的云雨雜波的頻率附近設(shè)置合適的虛擬雜波來增加云雨雜波處的阻帶寬度。例如地雜波譜中心20 Hz、譜寬2 Hz,云雨雜波譜中心130 Hz、譜寬20 Hz,生成濾波器系數(shù)時以譜中心分別為20 Hz、130 Hz、125 Hz(虛擬設(shè)置的雜波),它們的譜寬均為2 Hz,這樣虛擬設(shè)置的雜波(fd=125 Hz,σc=2 Hz)就起到了拓寬云雨雜波譜寬的作用。效果如圖2所示,虛擬雜波設(shè)置示意如圖3所示。注意濾波器的階數(shù)應(yīng)大于總的雜波數(shù)(真實雜波加虛擬設(shè)置的雜波)

    圖3 虛擬雜波設(shè)置示意圖

    圖4 特征矢量法處理多雜波時的改善因子

    圖5 改善因子對比

    由圖4清晰地可以看出特征矢量法處理兩個頻譜特征存在差異的情況下性能會惡化,且差異越大,性能惡化越嚴重。由圖5可以看出本文對特征矢量法的改進明顯會提高改善因子,性能的改善與設(shè)置的雜波的位置有關(guān),但是改善效果不會超過特征矢量法處理兩個譜寬相同的雜波的改善因子值。

    3 雜波特征提取

    要實現(xiàn)自適應(yīng)雜波抑制需要提取雜波頻譜特征(雜波中心頻率和雜波譜寬)來計算濾波器系數(shù)。傳統(tǒng)的自適應(yīng)濾波通常采用兩級級聯(lián)的濾波器:地雜波(零多普勒頻移)濾波器和運動雜波(海雜波或云雨雜波)濾波器,這種情況下通常只能估計一種雜波的頻譜特征,因此需要先通過地雜波濾波器濾除地雜波,以減少地雜波對運動雜波估計的影響。本文第一節(jié)已提到球載雷達的地雜波并不是固定的,隨著平臺的移動會展現(xiàn)類似運動雜波的特性,因此對球載雷達雜波而言,需要同時提取包括地雜波在內(nèi)的多個雜波的頻譜特征。由式(5)可知,雜波的頻譜特性由譜寬σc和功率譜中心fd決定。對于確定的σc和fd基于特征矢量法的濾波器權(quán)矢量也就確定了。對于球載雷達檢測低空低速目標而言,由于目標與雜波頻譜特征相近,所以雜波特征估計的準確性,決定了檢測的性能。

    因此,有必要采取方法同時對多個雜波的頻譜特征進行估計,同時也要采取措施盡可能減少估計誤差。使用相關(guān)函數(shù)法[13-14]可以估計多個雜波的特征,注意進行多雜波特征同時估計時,需要滿足多個雜波在不同距離,即在距離上雜波不重合。

    多雜波可表示為如下形式:

    (12)

    式中A(t)為復包絡(luò);Tk為第k個雜波存在范圍對應(yīng)的時間寬度;fdk為第k個雜波的多普勒頻率;φ0為初相;n(t)為加性白噪聲。延遲一個周期可寫成:

    u(t-Tr)=

    (j2πfdk(t-Tr)+φ0)+n(t-Tr)

    (13)

    一般情況下,噪聲與雜波不相關(guān),注意到各雜波之間由于時間上不同也均不相關(guān),所以它們的相關(guān)函數(shù)為:

    R(Tr)=E[u(t)u*(t-Tr)]=

    (14)

    式中B(Tr)=E[A(t)A(t-Tr)]。A(t)為窄帶信號,即A(t)=A(t-Tr),那么B(Tr)=E[|A(t)|2]為一實數(shù)。由式(14)可得

    (15)

    對一般遍歷平穩(wěn)隨機過程可由時間平均代替統(tǒng)計平均,得到估計式

    (16)

    式中i表示雜波在不同脈沖上的獨立采樣序列號。結(jié)合(15)(16)式得到雜波的多普勒頻率估計值

    (17)

    假設(shè)地雜波位于5~15 km處,譜中心為65 Hz,譜寬為11 Hz,云雨雜波位于25~80 km處,譜中心為130 Hz,譜寬為26 Hz,用相關(guān)函數(shù)法進行估計,結(jié)果如圖6所示。

    圖6 相關(guān)函數(shù)法估計結(jié)果

    由圖6可以看出利用相關(guān)函數(shù)法對雜波進行估計,頻率估計在相鄰單元有一定的差距,但均在預設(shè)值上下波動,也即存在一定的譜寬。取鄰近單元的測頻值求平均作為譜中心估計值,標準差作為譜寬估計值。表3給出了相關(guān)函數(shù)法估計多雜波譜寬譜中心的仿真統(tǒng)計結(jié)果。

    表3 相關(guān)函數(shù)法頻譜估計結(jié)果

    由表3的結(jié)果可知,相鄰函數(shù)法的譜中心估計效果較為理想,但是譜寬估計有很大誤差地雜波譜寬預設(shè)值為10 Hz,估計值為8.4~12.1 Hz,誤差為-1.6~2.1 Hz。而云雨雜波譜寬預設(shè)值為26 Hz,估計值為20.7~31.2 Hz,誤差為-5.3~5.2 Hz,相對誤差較大且可能偏大也可能偏小,這也是大多數(shù)頻譜估計方法的共性。

    4 球載雷達多雜波自適應(yīng)抑制

    4.1 譜寬覆蓋法

    考慮到譜寬σck的估計誤差太大,這對低速小目標的檢測是不能容忍的。因此要盡可能地減小誤差帶來的影響。在陸地上空球載雷達主要面臨地雜波和云雨雜波的干擾,由前面分析可知地雜波的綜合譜寬較窄一般為10 Hz左右,估計誤差在±2 Hz左右,云雨雜波的譜寬則一般為26 Hz左右,估計誤差在±5 Hz左右,可以看出譜寬越大,誤差越大。

    設(shè)相對誤差

    (18)

    地雜波的相對誤差δσc=0.19,云雨雜波的相對誤差為δσc=0.25,因此可以在譜寬估計的基礎(chǔ)上,對譜寬的取值進行量化拓展,形成多組濾波器同時進行濾波。譜寬的量化拓展與第二節(jié)對特征矢量法的改進是匹配的,通過設(shè)置虛擬雜波的形式實現(xiàn)。各個雜波的譜寬差異很大,理論上,地雜波的譜寬最小,云雨雜波的譜寬較大一般是地雜波的1~3倍,箔條等其他雜波一般是地雜波的1~2倍。因此選擇地雜波的頻譜特征值作為參考值。頻寬覆蓋法的步驟:

    (2)設(shè)置虛擬零點拓寬其他雜波的譜寬,設(shè)步進頻率Δfd=σc1/4,以云雨雜波fd2為例,考慮到計算量,在fd2左右各設(shè)置一個虛擬雜波,共設(shè)置三組fd2±2Δfd、fd2±3Δfd、fd2±4Δfd,三組虛擬雜波的設(shè)置對應(yīng)形成三組譜寬不同的濾波器,這就對雜波譜寬進行了覆蓋,很有可能將譜寬的真值包含在內(nèi)。虛擬雜波的個數(shù)和頻寬覆蓋的組數(shù)可以靈活設(shè)置,在這里考慮到計算量,設(shè)置的組數(shù)較少。依據(jù)前面所述,這三組可以求出三個濾波器系數(shù)。如某一種情況下的濾波器的頻率響應(yīng)如圖7所示。

    圖7 多組濾波器頻率響應(yīng)

    4.2 系統(tǒng)組成

    針對多雜波同時抑制問題,采用特征矢量法求解權(quán)系數(shù),自適應(yīng)雜波抑制需要估計多個雜波頻譜特征,由于估計譜寬存在較大誤差,該系統(tǒng)運用頻寬覆蓋法盡可能減小估計誤差的影響。系統(tǒng)組成如圖8所示。

    圖8 系統(tǒng)組成框圖

    系統(tǒng)由雜波存儲單元、雜波特征估計、權(quán)系數(shù)求解及擴展權(quán)系數(shù)、多組濾波器和一個比較器組成。

    具體步驟如下所示:

    (1)在本系統(tǒng)中,首先對雜波進行N個脈沖采樣得到ui(t),i=1,2,…N并進行存儲。

    (3)考慮到譜寬σck的估計誤差太大,利用頻寬覆蓋法形成多組濾波器,然后將雜波數(shù)據(jù)并行輸入到濾波器中,比較每一個濾波器的輸出結(jié)果,以輸出雜波剩余最小的濾波器作為最佳輸出。

    5 仿真分析

    為驗證多雜波抑制方法的有效性,分別對特征矢量法,改進特征矢量法雜波抑制性能進行對比分析。

    球載雷達高度H=3000 m;脈沖重復頻率fr=3000 Hz;載頻f0=3 GHz;天線俯仰波束寬度θ3=10°;天線水平波束寬度φ3=1.8°;設(shè)置為球載雷達平臺速度為vr=5 m/s,云雨雜波速度vc=15 m/s,對應(yīng)地雜波多普勒頻率與云雨雜波多普勒頻率分別為100 Hz與300 Hz,其中地雜波的譜寬為設(shè)置為10 Hz,云雨雜波譜寬為30Hz。在5 km處設(shè)置一目標,目標的高度為500 m,速度為10 m/s,對應(yīng)多普勒頻率為200 Hz,平均信雜比為設(shè)置為-10 dB。

    圖9為濾波前雜波功率譜示意圖,圖9(a)可以看到地雜波與云雨雜波,由于目標功率小于地雜波,所以目標淹沒于地雜波中無法分辨,圖9(b)為目標所在距離單元的功率譜,可以看到目標基本位于雜波副瓣內(nèi),較難分辨。

    圖9 濾波前結(jié)果示意圖

    圖10 特征矢量法濾波后結(jié)果示意圖

    圖10是采用特征矢量法濾波后結(jié)果示意圖,從圖10(a),圖10(b)均無法分辨目標信息,這是由于該方法處理多個不同雜波時性能惡化導致,在云雨雜波處形成了寬的凹口對目標也進行了削弱,導致信雜比的改善并不理想。

    圖11是采用改進的特征矢量法濾波后結(jié)果示意圖,從圖11(a)可以看到一點目標的信息,從圖11(b)可以較為清晰地分辨目標信息。由前面分析已知,改進的特則矢量法可以較好的匹配雜波頻率與譜寬,盡可能的保留目標信息,對信雜比改善的效果較好。

    圖11 改進特征矢量法濾波后結(jié)果示意圖

    6 結(jié) 語

    球載雷達雜波環(huán)境的特殊性導致了傳統(tǒng)的自適應(yīng)雜波抑制方法不適應(yīng),本文分析了球載雷達雜波特性,針對多個雜波同時抑制的情況采用基于最大平均改善因子的特征矢量法計算濾波器系數(shù),特征矢量法在處理多個頻譜特性差異較大的情況下性能會惡化,提出了通過設(shè)置虛擬雜波的方法改善其性能,通過仿真驗證其性能有較好的改善;基于相關(guān)函數(shù)法推導了多雜波頻譜估計的相關(guān)理論,針對估計誤差較大的問題,結(jié)合改進的特征矢量法具有可以自由設(shè)置虛擬雜波來拓展譜寬的特點,在雜波估計的譜中心附近形成多種譜寬覆蓋。該系統(tǒng)比傳統(tǒng)雜波抑制相比可以同時處理多個未知的雜波,且具備減小誤差的優(yōu)點。本文方法含義明確,計算也不復雜,可操作性強,對球載雷達具有較強的使用價值。

    猜你喜歡
    云雨雜波矢量
    STAR2000型空管一次雷達雜波抑制淺析
    矢量三角形法的應(yīng)用
    青山一道同云雨,明月何曾是兩鄉(xiāng)
    學生天地(2020年15期)2020-08-25 09:21:50
    青山一道同云雨 明月何曾是兩鄉(xiāng)(書法)
    陽朔云雨
    基于矢量最優(yōu)估計的穩(wěn)健測向方法
    三角形法則在動態(tài)平衡問題中的應(yīng)用
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    相關(guān)廣義復合分布雷達海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    微波雷達海雜波混沌特性分析
    青春草亚洲视频在线观看| 最近最新中文字幕大全免费视频| 自线自在国产av| 国产无遮挡羞羞视频在线观看| 一本综合久久免费| 国产精品99久久99久久久不卡| 国产精品免费视频内射| 日日摸夜夜添夜夜添小说| 美女福利国产在线| 一区二区三区四区激情视频| 女人精品久久久久毛片| 91精品国产国语对白视频| 中文字幕最新亚洲高清| 亚洲成国产人片在线观看| www.熟女人妻精品国产| 大片免费播放器 马上看| 热99国产精品久久久久久7| 国产成人欧美| 两个人免费观看高清视频| videosex国产| 黄网站色视频无遮挡免费观看| 天天躁夜夜躁狠狠躁躁| 国产真人三级小视频在线观看| 美女中出高潮动态图| 午夜福利视频在线观看免费| 国产成人欧美| 交换朋友夫妻互换小说| 大陆偷拍与自拍| 交换朋友夫妻互换小说| 永久免费av网站大全| 桃红色精品国产亚洲av| 国产欧美日韩一区二区精品| 午夜精品久久久久久毛片777| 高清黄色对白视频在线免费看| 一本综合久久免费| 自拍欧美九色日韩亚洲蝌蚪91| 久久香蕉激情| 亚洲国产欧美在线一区| 国产成人啪精品午夜网站| 亚洲av成人一区二区三| 激情视频va一区二区三区| 97在线人人人人妻| 1024香蕉在线观看| 精品国产乱码久久久久久男人| 日韩大码丰满熟妇| 国产成人欧美在线观看 | 久久久国产成人免费| 一区二区日韩欧美中文字幕| 国产在视频线精品| 永久免费av网站大全| 最近中文字幕2019免费版| 在线看a的网站| 久久精品人人爽人人爽视色| 国产av精品麻豆| 午夜福利在线观看吧| 电影成人av| 丝袜美腿诱惑在线| 国产精品av久久久久免费| 国产精品一区二区在线观看99| 国产一区二区三区综合在线观看| 国产三级黄色录像| 91九色精品人成在线观看| 一级,二级,三级黄色视频| 狠狠婷婷综合久久久久久88av| 女人久久www免费人成看片| 日本a在线网址| 欧美xxⅹ黑人| 中国美女看黄片| 亚洲精品日韩在线中文字幕| 50天的宝宝边吃奶边哭怎么回事| 爱豆传媒免费全集在线观看| 欧美大码av| 91大片在线观看| 免费在线观看日本一区| 午夜福利视频精品| 久久影院123| 亚洲avbb在线观看| 亚洲avbb在线观看| 高清在线国产一区| 国产高清国产精品国产三级| 99久久精品国产亚洲精品| 黑人猛操日本美女一级片| 欧美日韩国产mv在线观看视频| 亚洲专区国产一区二区| 精品人妻1区二区| 一本一本久久a久久精品综合妖精| 欧美激情高清一区二区三区| 国内毛片毛片毛片毛片毛片| 手机成人av网站| 国产精品熟女久久久久浪| 少妇粗大呻吟视频| 男女午夜视频在线观看| 中文字幕人妻丝袜一区二区| 国产av精品麻豆| 青草久久国产| 女人精品久久久久毛片| 亚洲av片天天在线观看| 精品少妇一区二区三区视频日本电影| 丝袜在线中文字幕| 男女午夜视频在线观看| 国产淫语在线视频| 老司机影院毛片| 国产精品偷伦视频观看了| 国产精品免费大片| 午夜免费鲁丝| 黄色毛片三级朝国网站| 在线 av 中文字幕| 12—13女人毛片做爰片一| 婷婷成人精品国产| 久久国产精品男人的天堂亚洲| 19禁男女啪啪无遮挡网站| 老鸭窝网址在线观看| 在线av久久热| 亚洲国产精品一区三区| 91大片在线观看| 色精品久久人妻99蜜桃| 亚洲成人国产一区在线观看| 中文字幕高清在线视频| av免费在线观看网站| 精品一品国产午夜福利视频| 纵有疾风起免费观看全集完整版| 成年女人毛片免费观看观看9 | 天堂俺去俺来也www色官网| avwww免费| 亚洲性夜色夜夜综合| 十八禁人妻一区二区| 日韩有码中文字幕| 777米奇影视久久| 欧美大码av| 99久久综合免费| 午夜福利在线观看吧| 在线精品无人区一区二区三| 性色av乱码一区二区三区2| 波多野结衣一区麻豆| 色综合欧美亚洲国产小说| 国产成人精品久久二区二区免费| 色老头精品视频在线观看| 日韩 亚洲 欧美在线| 午夜视频精品福利| 日韩中文字幕欧美一区二区| 亚洲专区国产一区二区| 黄色 视频免费看| 精品亚洲乱码少妇综合久久| 久久av网站| av网站免费在线观看视频| 亚洲五月婷婷丁香| 精品一品国产午夜福利视频| 日韩精品免费视频一区二区三区| 国产亚洲欧美在线一区二区| 久久久国产精品麻豆| 91国产中文字幕| 黄片播放在线免费| 天堂俺去俺来也www色官网| 国产精品香港三级国产av潘金莲| 狠狠婷婷综合久久久久久88av| 欧美日韩福利视频一区二区| 黄色怎么调成土黄色| 婷婷色av中文字幕| 欧美黄色片欧美黄色片| 欧美精品一区二区大全| 国产麻豆69| 黄片大片在线免费观看| 91国产中文字幕| 在线看a的网站| 一区二区三区精品91| 免费不卡黄色视频| 久久青草综合色| 午夜精品国产一区二区电影| 国产淫语在线视频| 色婷婷久久久亚洲欧美| 男女午夜视频在线观看| 五月开心婷婷网| 亚洲美女黄色视频免费看| 色婷婷av一区二区三区视频| 一本大道久久a久久精品| 免费少妇av软件| av网站免费在线观看视频| 国产亚洲午夜精品一区二区久久| 又紧又爽又黄一区二区| 久久精品国产亚洲av香蕉五月 | 日本猛色少妇xxxxx猛交久久| 下体分泌物呈黄色| 天天躁狠狠躁夜夜躁狠狠躁| 19禁男女啪啪无遮挡网站| 18在线观看网站| 国产人伦9x9x在线观看| 动漫黄色视频在线观看| 亚洲精品美女久久久久99蜜臀| 欧美xxⅹ黑人| 国产黄频视频在线观看| 99久久精品国产亚洲精品| av在线播放精品| 午夜91福利影院| 日本91视频免费播放| 久久人妻熟女aⅴ| 亚洲av男天堂| 少妇猛男粗大的猛烈进出视频| 亚洲国产精品成人久久小说| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕人妻丝袜制服| 精品国产一区二区久久| 丁香六月欧美| 飞空精品影院首页| 久久av网站| 久久久国产成人免费| 91精品国产国语对白视频| 精品亚洲成国产av| 欧美日韩视频精品一区| 一级黄色大片毛片| 免费观看人在逋| 久久国产精品人妻蜜桃| 精品视频人人做人人爽| 黄色片一级片一级黄色片| 国产精品99久久99久久久不卡| 亚洲中文av在线| 熟女少妇亚洲综合色aaa.| 久久久久久人人人人人| 亚洲黑人精品在线| 丁香六月欧美| 国产亚洲精品第一综合不卡| 蜜桃在线观看..| 亚洲成人免费av在线播放| 国产亚洲精品第一综合不卡| 久久精品国产综合久久久| 男人舔女人的私密视频| 叶爱在线成人免费视频播放| av网站在线播放免费| 性色av一级| 亚洲av美国av| 中文欧美无线码| 亚洲av电影在线观看一区二区三区| 国产成人欧美| 国产精品自产拍在线观看55亚洲 | 国产xxxxx性猛交| 欧美激情高清一区二区三区| 两性夫妻黄色片| 亚洲精品国产一区二区精华液| 一区二区三区四区激情视频| 国产精品久久久久久精品电影小说| 国产深夜福利视频在线观看| 免费在线观看完整版高清| 亚洲美女黄色视频免费看| 国产1区2区3区精品| 中文字幕制服av| 国产精品.久久久| 免费在线观看黄色视频的| 操美女的视频在线观看| 人妻人人澡人人爽人人| 国产精品麻豆人妻色哟哟久久| 精品国产乱码久久久久久男人| 久久人人爽av亚洲精品天堂| 一区二区日韩欧美中文字幕| 亚洲第一欧美日韩一区二区三区 | 精品亚洲乱码少妇综合久久| 最黄视频免费看| 97人妻天天添夜夜摸| 黑丝袜美女国产一区| 我的亚洲天堂| 黄片小视频在线播放| 各种免费的搞黄视频| 美女视频免费永久观看网站| 一区二区三区精品91| 久久久久网色| 老汉色av国产亚洲站长工具| 日韩一卡2卡3卡4卡2021年| 成年人午夜在线观看视频| 久久狼人影院| 一级,二级,三级黄色视频| 一区二区av电影网| 国产精品九九99| 免费观看a级毛片全部| 亚洲欧美日韩高清在线视频 | 啦啦啦中文免费视频观看日本| 纯流量卡能插随身wifi吗| 日韩三级视频一区二区三区| 免费高清在线观看视频在线观看| 久久国产精品人妻蜜桃| 在线观看舔阴道视频| 嫩草影视91久久| 国产一区二区三区在线臀色熟女 | 亚洲精品乱久久久久久| 国产亚洲精品一区二区www | 午夜福利在线观看吧| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 国产色视频综合| a 毛片基地| 精品少妇黑人巨大在线播放| 免费在线观看黄色视频的| 在线观看舔阴道视频| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久毛片微露脸 | 免费久久久久久久精品成人欧美视频| 精品少妇一区二区三区视频日本电影| 日本精品一区二区三区蜜桃| 夫妻午夜视频| 午夜影院在线不卡| 亚洲激情五月婷婷啪啪| 国产99久久九九免费精品| av在线app专区| 91成年电影在线观看| 好男人电影高清在线观看| 叶爱在线成人免费视频播放| 亚洲精品乱久久久久久| 亚洲精品国产av成人精品| 在线观看人妻少妇| 亚洲色图 男人天堂 中文字幕| 99久久99久久久精品蜜桃| 久久天堂一区二区三区四区| 国产97色在线日韩免费| 亚洲中文字幕日韩| 成人亚洲精品一区在线观看| 日韩一卡2卡3卡4卡2021年| 又黄又粗又硬又大视频| 久久久久精品人妻al黑| 777米奇影视久久| 韩国高清视频一区二区三区| 五月开心婷婷网| 国产成人欧美在线观看 | 国产精品二区激情视频| 国产av又大| netflix在线观看网站| 青春草视频在线免费观看| 成年人黄色毛片网站| 人妻人人澡人人爽人人| 久久久欧美国产精品| 一区二区av电影网| 水蜜桃什么品种好| 国产在线免费精品| 大片电影免费在线观看免费| 动漫黄色视频在线观看| 美女脱内裤让男人舔精品视频| 国产精品 欧美亚洲| 成人av一区二区三区在线看 | 不卡av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 十分钟在线观看高清视频www| 97精品久久久久久久久久精品| 亚洲中文av在线| 熟女少妇亚洲综合色aaa.| 激情视频va一区二区三区| 成年人免费黄色播放视频| 久久久久久久国产电影| 国产不卡av网站在线观看| a在线观看视频网站| 亚洲国产欧美网| 亚洲国产成人一精品久久久| 国产成人精品久久二区二区91| 国产99久久九九免费精品| 国产一区有黄有色的免费视频| 男女下面插进去视频免费观看| 人人妻人人添人人爽欧美一区卜| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 国产99久久九九免费精品| 美女中出高潮动态图| 免费在线观看黄色视频的| 水蜜桃什么品种好| 久久女婷五月综合色啪小说| 丰满人妻熟妇乱又伦精品不卡| 嫩草影视91久久| 18在线观看网站| 国产亚洲精品一区二区www | 黑人猛操日本美女一级片| 狂野欧美激情性bbbbbb| 国产又色又爽无遮挡免| 久久精品成人免费网站| 欧美大码av| 少妇粗大呻吟视频| 女人久久www免费人成看片| 国产欧美日韩综合在线一区二区| 夜夜骑夜夜射夜夜干| 精品视频人人做人人爽| 亚洲一卡2卡3卡4卡5卡精品中文| 国产1区2区3区精品| 国产精品影院久久| 纵有疾风起免费观看全集完整版| 亚洲精品国产色婷婷电影| 国产精品99久久99久久久不卡| 日韩精品免费视频一区二区三区| a在线观看视频网站| 久久久久久久久免费视频了| 久久久久久免费高清国产稀缺| 91成年电影在线观看| 欧美另类一区| 亚洲精品国产区一区二| 国产免费福利视频在线观看| 一级毛片精品| kizo精华| 18禁裸乳无遮挡动漫免费视频| 国产无遮挡羞羞视频在线观看| 男女免费视频国产| 精品少妇内射三级| 深夜精品福利| 欧美另类亚洲清纯唯美| 亚洲 国产 在线| a在线观看视频网站| 免费日韩欧美在线观看| 日韩熟女老妇一区二区性免费视频| 欧美精品av麻豆av| 国产亚洲av片在线观看秒播厂| 男人操女人黄网站| 老司机在亚洲福利影院| 叶爱在线成人免费视频播放| 啦啦啦中文免费视频观看日本| 午夜福利一区二区在线看| 在线av久久热| 好男人电影高清在线观看| 大香蕉久久成人网| 国产精品久久久久久精品电影小说| 国产男女超爽视频在线观看| 我的亚洲天堂| 考比视频在线观看| 老司机午夜福利在线观看视频 | 侵犯人妻中文字幕一二三四区| 天堂中文最新版在线下载| 中文字幕av电影在线播放| 1024香蕉在线观看| 高清av免费在线| 精品人妻一区二区三区麻豆| 欧美人与性动交α欧美精品济南到| 777久久人妻少妇嫩草av网站| 丝袜美腿诱惑在线| 在线观看免费高清a一片| 深夜精品福利| 伦理电影免费视频| 男女无遮挡免费网站观看| 国产成人免费观看mmmm| 老熟妇仑乱视频hdxx| 伊人亚洲综合成人网| av天堂久久9| 欧美人与性动交α欧美精品济南到| 美国免费a级毛片| kizo精华| 一本—道久久a久久精品蜜桃钙片| 久久久久久人人人人人| 日韩人妻精品一区2区三区| 大片免费播放器 马上看| 日韩,欧美,国产一区二区三区| 国产人伦9x9x在线观看| 波多野结衣av一区二区av| 国产精品成人在线| 高清欧美精品videossex| 日韩三级视频一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 国产野战对白在线观看| 日本一区二区免费在线视频| 欧美亚洲日本最大视频资源| 亚洲精品久久久久久婷婷小说| 欧美日韩精品网址| 亚洲av日韩在线播放| 人妻人人澡人人爽人人| 久久精品人人爽人人爽视色| 国产99久久九九免费精品| 美女午夜性视频免费| 黄色视频不卡| 欧美少妇被猛烈插入视频| 精品人妻1区二区| 成年人黄色毛片网站| 9热在线视频观看99| 欧美日韩中文字幕国产精品一区二区三区 | 日本av免费视频播放| 亚洲精品第二区| 99国产精品一区二区蜜桃av | 成人手机av| 狂野欧美激情性bbbbbb| 日本av免费视频播放| 国产在线一区二区三区精| 午夜激情久久久久久久| 99精品欧美一区二区三区四区| 午夜福利视频精品| 国产成人精品无人区| 建设人人有责人人尽责人人享有的| 美女主播在线视频| 欧美精品高潮呻吟av久久| 正在播放国产对白刺激| 欧美日韩中文字幕国产精品一区二区三区 | 黄色怎么调成土黄色| 久久国产精品大桥未久av| 国产精品久久久久久精品电影小说| 久久性视频一级片| av天堂在线播放| 国产精品免费大片| 91精品三级在线观看| 天天躁日日躁夜夜躁夜夜| 欧美精品人与动牲交sv欧美| 欧美 日韩 精品 国产| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| 啦啦啦中文免费视频观看日本| 欧美激情高清一区二区三区| 亚洲熟女毛片儿| 操出白浆在线播放| 美女中出高潮动态图| 国产成人av激情在线播放| 国产一区二区三区在线臀色熟女 | 欧美午夜高清在线| 免费一级毛片在线播放高清视频 | 夜夜夜夜夜久久久久| 日韩欧美一区二区三区在线观看 | 国产成人一区二区三区免费视频网站| 美女中出高潮动态图| 高潮久久久久久久久久久不卡| 欧美日韩亚洲综合一区二区三区_| 久久国产精品影院| 热99re8久久精品国产| av欧美777| 我的亚洲天堂| 男女之事视频高清在线观看| 黄网站色视频无遮挡免费观看| 极品人妻少妇av视频| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 亚洲黑人精品在线| 999久久久国产精品视频| av网站在线播放免费| 国产欧美日韩精品亚洲av| 高清欧美精品videossex| 一区二区三区精品91| 新久久久久国产一级毛片| 久久ye,这里只有精品| 成人亚洲精品一区在线观看| 国产精品1区2区在线观看. | 日韩视频一区二区在线观看| 亚洲专区国产一区二区| 精品久久久久久电影网| 考比视频在线观看| 亚洲国产欧美网| 麻豆av在线久日| 老熟女久久久| 成人影院久久| 丝袜美腿诱惑在线| 热re99久久精品国产66热6| 日本撒尿小便嘘嘘汇集6| 欧美午夜高清在线| 欧美乱码精品一区二区三区| 又大又爽又粗| 午夜免费成人在线视频| 成人18禁高潮啪啪吃奶动态图| 99精国产麻豆久久婷婷| 久久人人爽人人片av| 国产麻豆69| 欧美午夜高清在线| 一本综合久久免费| 国产一区二区三区av在线| 丝袜喷水一区| 欧美久久黑人一区二区| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| 美女扒开内裤让男人捅视频| 成人免费观看视频高清| av免费在线观看网站| 国产麻豆69| 国产91精品成人一区二区三区 | 欧美av亚洲av综合av国产av| 精品国产乱子伦一区二区三区 | 999精品在线视频| 亚洲综合色网址| 国产精品九九99| av免费在线观看网站| 亚洲伊人久久精品综合| 亚洲国产精品999| 亚洲专区字幕在线| 久久天躁狠狠躁夜夜2o2o| 国产日韩欧美视频二区| 久久久国产成人免费| 久久人人97超碰香蕉20202| 国产一区二区在线观看av| 精品亚洲成a人片在线观看| 国产男女内射视频| 久久久欧美国产精品| 国产精品一区二区免费欧美 | 亚洲色图综合在线观看| 90打野战视频偷拍视频| 亚洲精品国产av成人精品| 满18在线观看网站| 亚洲精品一卡2卡三卡4卡5卡 | 久久99热这里只频精品6学生| 国产亚洲欧美在线一区二区| 又黄又粗又硬又大视频| 亚洲精品自拍成人| 国产无遮挡羞羞视频在线观看| 久久香蕉激情| 男女午夜视频在线观看| 久久久精品94久久精品| 免费一级毛片在线播放高清视频 | 久久 成人 亚洲| 色94色欧美一区二区| 国产1区2区3区精品| 免费观看av网站的网址| 精品少妇一区二区三区视频日本电影| 免费观看a级毛片全部| 国产在线一区二区三区精| 夜夜夜夜夜久久久久| 日韩电影二区| 9色porny在线观看| 欧美日本中文国产一区发布| 国产精品久久久久成人av| 成年美女黄网站色视频大全免费| 欧美精品亚洲一区二区| 久久性视频一级片| 男女午夜视频在线观看| 中文字幕最新亚洲高清| 免费高清在线观看日韩| 国产精品免费视频内射| 视频区图区小说| 精品视频人人做人人爽| 青青草视频在线视频观看| 视频区欧美日本亚洲| 免费日韩欧美在线观看| 国产成人影院久久av| 婷婷成人精品国产| 三级毛片av免费| 久久国产精品男人的天堂亚洲| 亚洲欧美日韩高清在线视频 |