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

    統(tǒng)計MIMO雷達目標跟蹤技術(shù)研究

    2014-10-13 07:59:40
    艦船電子對抗 2014年3期
    關(guān)鍵詞:參數(shù)估計雷達噪聲

    盧 毅

    (解放軍77108部隊,成都611233)

    0 引 言

    多輸入多輸出(MIMO)雷達是近幾年提出的一種新體制雷達[1]。與傳統(tǒng)雷達相比,MIMO雷達在目標檢測、參數(shù)估計、雜波抑制、抗干擾等方面具有顯著優(yōu)勢[1-4],因而成為學術(shù)界研究的熱點。國內(nèi)外對MIMO雷達的研究大致可以分為2條技術(shù)路線:第1類是收發(fā)天線大間距配置的統(tǒng)計MIMO雷達[2];第2類是收發(fā)單元緊密分布的相干MIMO雷達[3]。統(tǒng)計 MIMO雷達的陣元采用空間分集配置方式,發(fā)射波從不同角度照射目標,充分利用空域信息來克服目標的閃爍效應,獲得空間分集增益,進而提高閃爍目標的檢測性能,而且發(fā)射波形采用正交波形,還可同時獲得波形分集增益。統(tǒng)計MIMO雷達的信號處理方式有相參處理和非相參處理,2種處理方式都要求收發(fā)單元時間同步,其中相參處理還要求相位同步。這里主要研究非相參統(tǒng)計MIMO雷達對運動目標的跟蹤技術(shù)。

    目前,對MIMO雷達的研究主要集中在目標檢測[4-5]、參數(shù)估計[6-7]、目標跟蹤[8-10]和優(yōu)化布站[7]等方面。文獻[6]研究了MIMO雷達對目標參數(shù)的聯(lián)合估計,推導了聯(lián)合估計的最大似然表達式和相應的克拉美羅界(CRLB),理論分析和仿真驗證后得出陣元數(shù)越多估計性能越好的結(jié)論。文獻[7]對分布式非相參MIMO雷達的目標跟蹤技術(shù)進行了研究,在MIMO雷達估計目標參數(shù)的基礎上,分別基于卡爾曼算法和粒子濾波算法對運動目標進行了跟蹤,但是文章考慮的情況忽略了目標運動對回波信號造成的多普勒頻移,跟蹤時也沒有將速度的估值做為觀測數(shù)據(jù)加以考慮。文獻[8]提出了一種將信號時延和多普勒頻移的估計值做為觀測結(jié)果對目標進行跟蹤的方法,推導了相應的后驗克拉美羅界(PCRLB),分析了不同布站方式對 MIMO雷達跟蹤性能的影響,但因需要被估計的參數(shù)較多,不易實現(xiàn),文章僅僅分析了性能并沒有進行實際的跟蹤,也沒有跟蹤誤差的分析。文獻[9]基于MIMO雷達特有的物理結(jié)構(gòu)和信號處理方式提出了2種跟蹤方法:集中式跟蹤和分布式跟蹤,推導了基于這2種方法的PCRLB,分析了這2種跟蹤技術(shù)的最優(yōu)性能。

    對于固定布站的雷達,目標運動是造成雷達橫截面積(RCS)閃爍的直接原因,因此研究MIMO雷達對運動目標的探測和跟蹤是非常迫切的。目前,統(tǒng)計MIMO雷達均采用全向發(fā)射,能量利用率較低,研究MIMO跟蹤技術(shù)可以為實現(xiàn)雷達定向發(fā)射提供支持,提高雷達的探測威力;現(xiàn)階段,MIMO雷達的參數(shù)估計大都采用最大似然搜索,搜索范圍很大,效率較低,跟蹤濾波不僅可以提高參數(shù)估計精度,而且對每一步的預測值加以利用可以極大地減輕下一步參數(shù)估計的復雜度。本文深刻地分析了MIMO雷達在目標跟蹤中出現(xiàn)的問題,從解決實際問題的角度出發(fā),研究了MIMO雷達對運動目標的參數(shù)估計和跟蹤技術(shù),提出了基于MIMO雷達特有結(jié)構(gòu)和信號處理方式的跟蹤方法。

    1 信號模型

    假設1部分布式MIMO雷達有M個發(fā)射站、N個接收站。為了抑制目標RCS閃爍,收發(fā)站間滿足文獻[11]中的空間分集條件,第k個發(fā)射站的平面位置坐標為(xtk,ytk),第l號接收站的位置位于(xrl,yrl),其中k=1,2,…,M ;l=1,2,…,N 。

    設m 時刻在坐標 (xm,ym)處有一目標正以(vx,vy)的速度勻速運動,此時從k號發(fā)射機發(fā)射的信號經(jīng)目標反射到達l號接收機所經(jīng)歷的時間延遲為τkl,多普勒頻移為fkl,則:

    假設目標的散射點數(shù)目非常多,且沒有哪一散射點占主導地位,則從k號發(fā)射機發(fā)射的信號經(jīng)目標反射被l號接收機接收后可以表示為sk(t)的延遲和一復高斯隨機變量的乘積。這個隨機變量就是k-l收發(fā)通道的通道系數(shù)αkl,它是傳輸通道衰減、目標RCS閃爍和相位移動影響的總和,服從αkl~CN(0,σ2)分布,其中σ2為通道系數(shù)的方差。

    第l號接收機的接收信號是目標所有散射點向這個方向反射信號與噪聲的疊加,用rl(t)表示l號接收機的接收信號,nl(t)表示接收機的內(nèi)部噪聲,則:

    由于MIMO雷達發(fā)射正交信號,任意2部發(fā)射機的發(fā)射信號在任意的時延情況下均能保持很好的正交性,則:

    但是由于本文考慮了目標運動給回波信號帶來的多普勒頻移,因此,在相鄰發(fā)射信號的多普勒頻移小于信號帶寬的情況下(即在滿足|fkl-f(k-1)l|<Bk時),回波信號的正交性依然不會改變。如果雷達參數(shù)選擇得當,這種情況在實際中一般是滿足的,則:

    在接收端正是利用了波形的正交性,將接收機中M個通道的信號分離開,得到每一收發(fā)通道的接收信號:

    通常噪聲為帶限噪聲N0≠,這一點在仿真時需要特別注意。

    在處理中心將各通道接收信號組合成總的接收矢量:

    2 動目標參數(shù)估計和性能分析

    2.1 最大似然法估計目標參數(shù)

    獲得接收信號矢量以后,如何從中估計出目標參數(shù)將是這一節(jié)的研究重點。下面的研究中,根據(jù)上一節(jié)建立的信號模型推導了參數(shù)估計的最大似然表達式,這也是本文結(jié)合雷達測量與跟蹤的第1步。設目標位置參數(shù)為θ= [x,y,vx,vy]T,θ∈Θ,Θ 是一個包括θ所有可能取值的四維空間。

    根據(jù)上一節(jié)建立的信號模型,可以得到k-l通道接收信號關(guān)于目標參數(shù)θ和通道系數(shù)αkl的條件似然比函數(shù):

    式中:F為與目標參數(shù)無關(guān)的常數(shù)。

    通過計算式(8)關(guān)于通道系數(shù)的期望,可以得到k-l通道接收信號關(guān)于目標參數(shù)θ的似然比函數(shù):

    由式(10)可以看出,k-l通道接收信號的似然比函數(shù)表達式中已不再包含通道系數(shù)αkl。考慮到通道噪聲的不相關(guān)性和各通道系數(shù)的獨立性,在下一步可以推導出接收矢量的聯(lián)合似然函數(shù):

    因此最大似然估計的表達式為:

    2.2 CRLB推導

    參數(shù)估計另一個重要方面就是分析估計性能,而參數(shù)估計性能最好的衡量標準就是CRLB,這一節(jié)將在上一節(jié)的基礎上推導各模型下的CRLB,這里的推導,而且在本文所提的算法中,克拉美羅矩陣是測量噪聲協(xié)方差矩陣的最佳選擇。

    假定最大似然估計結(jié)果為未知參數(shù)矢量θ的無偏估計:

    式中:J(θ)為費什信息矩陣(FIM)。 根據(jù)推導[6]費什信息矩陣為:

    3 MIMO雷達目標跟蹤

    式中:I2為二維單位矩陣;?表示Kronecker積;T為觀測間隔;Vk為動態(tài)模型中的過程噪聲,在這里認為它是零均值高斯白噪聲,用Q表示其協(xié)方差矩陣:

    式中:q為過程噪聲的功率譜密度。

    因為在運動過程中目標的RCS閃爍特別嚴重,因此在跟蹤過程中將MIMO雷達的估計(測量)數(shù)據(jù)作為測量方程的觀測值,這樣不僅可以抵抗跟蹤過程中的RCS閃爍,還可以得到精度較高的觀測數(shù)據(jù)。下面結(jié)合MIMO雷達的參數(shù)估計建立測量方程:

    式中:H為觀測矩陣,與文獻[7]中不同的是測量值和測量噪聲均為四維向量,測量矩陣為四維單位矩陣;Wk為MIMO雷達的測量噪聲,用Rk表示其協(xié)方差矩陣,根據(jù)雷達工作在漸近區(qū)的假設,取k時刻的測量噪聲協(xié)方差矩陣為MIMO雷達估計目標k時刻運動狀態(tài)的克拉美羅矩陣與一個系數(shù)β的乘積,通常1<β≤2。

    要注意的是,在實際計算測量噪聲協(xié)方差矩陣時,通常取目標k時刻運動狀態(tài)的預測值帶入CRLB公式算得或者用k時刻運動目標參數(shù)的估計值計算:

    需要說明的是,在后面的仿真中MIMO雷達工作在文獻[6]定義的漸近區(qū)域,因此 MIMO雷達所估計出的參數(shù)是有效的無偏估計值,可以將觀測噪聲視為高斯白噪聲。綜上所述,目標的動態(tài)方程和雷達的觀測方程均是線性高斯的,可以直接采用經(jīng)典的卡爾曼算法實施跟蹤。

    上面的方法可以實現(xiàn)MIMO雷達對目標跟蹤的良好性能,但是實現(xiàn)MIMO雷達跟蹤的意義遠不限于此。在跟蹤時采用交互式信號處理方法,利用狀態(tài)的預測值Xk|k-1為k時刻的目標參數(shù)估計設定一個搜索范圍(稱為置信區(qū)域),讓MIMO雷達在信號處理時在這個置信區(qū)域內(nèi)搜索目標,這樣可以減輕MIMO雷達參數(shù)估計的復雜度,而且不會影響參數(shù)估計的精度,具有很高的適用價值,這就是文章里一直提到的交互式信號處理方法,即跟蹤與參數(shù)估計的交互。

    試驗中具體的做法是,以狀態(tài)預測值為中心建立一個矩陣,將矩陣離散成許多孤立的點,在下一步的MIMO雷達參數(shù)估計時,僅通過比較這些點的似然比來得到目標位置的估計值。在實際應用時置信區(qū)域的范圍設定顯得很重要,經(jīng)過多次試驗和分析建議采用變范圍搜索:起始時刻采用大范圍、多點數(shù)搜索,隨著跟蹤步數(shù)的遞增,精度逐步提高,預測值也越貼近真實位置,在不影響估計精度的前提下縮小范圍、減少估計點數(shù),直至穩(wěn)定,具體操作情況在后面的仿真中有進一步說明。

    綜上,具備了MIMO雷達實施跟蹤所需的所有條件,下面利用Kalman算法對目標實施跟蹤,具體步驟為:

    第1步:首先在目標RCS閃爍的假設下,利用MIMO雷達在監(jiān)視空域搜索到目標,估計出其運動參數(shù)矢量Y0,并根據(jù)目標參數(shù)測量值算得該值估計的CRLB矩陣J(θ0)。

    第2步:將J(θ0)帶入公式(24),得到濾波的初始協(xié)方差矩陣R0。

    第3步:令X0=Y(jié)0,R0=J(θ0)。則:

    第4步:根據(jù)目標的預測狀態(tài)X1|0,以它的值為中心建立該時刻MIMO雷達參數(shù)估計的置信區(qū),并將置信區(qū)離散成孤立的點,對雷達接收信號在這些點進行最大似然估計,得到目標該時刻參數(shù)的測量矢量Y1,并將Y1帶入公式(20)計算該時刻雷達估計(測量)誤差的協(xié)方差矩陣R1。

    第5步:在得到雷達的估計(測量)數(shù)據(jù)和誤差矩陣后,利用濾波算法進行數(shù)據(jù)處理:

    第6步:根據(jù)上面的方法繼續(xù)遞歸,直至跟蹤結(jié)束。

    4 仿真驗證

    使用Matlab軟件對前面的研究內(nèi)容進行仿真,更直觀地展現(xiàn)本文算法的優(yōu)勢。設二維平面上有一目標做勻速直線運動,速度為(80m/s,20m/s),初始位置位于坐標原點。因為實際中目標的運動受到各種環(huán)境因素的影響,運動并不是嚴格勻速的,在仿真時給勻速運動的目標加一個零均值高斯分布的加速度,服從a~CN(0,0.9)分布,得到其真實軌跡,如圖1所示。實驗中設定MIMO雷達收發(fā)站均勻分布在第一象限內(nèi)一段半徑為90km的圓弧上,近似為收發(fā)共址配置。同時,為了減少參數(shù)估計時的運算量,在后面的仿真實驗中均采用3×3的典型布站方式 (具體分布情況見圖2),并假定各收發(fā)通道的信噪比相同,在仿真時均設定為10dB。

    圖1 收發(fā)站與目標運動軌跡的空間分布情況

    圖2 目標的真實運動軌跡

    發(fā)射信號采用頻率擴展的高斯脈沖,第k個發(fā)射站發(fā)射信號的復包絡形式如下:

    式中:T為決定高斯脈沖有效持續(xù)時間的參數(shù)。

    在下面的仿真中取T=5×10-2s,這與文獻[7]中T=1.125 4×10-7s的取值差別甚大,主要原因是本文中MIMO雷達估計的參數(shù)是位置和速度,二者的估計精度與T的關(guān)系正好是相反的,如果繼續(xù)按文獻[7]中的T取值可以得到很高的位置估計精度,但是速度估計精度會變得很差,下面的仿真圖展示了這一情況。Δf是相鄰發(fā)射信號之間的頻移參數(shù),只要Δf大于發(fā)射信號的帶寬,就可以保證各發(fā)射波形之間相互正交,因為高斯脈沖的帶寬與T值緊密相關(guān),所以仿真中Δf值的設定要與T相對應。

    分析圖3、圖4可以發(fā)現(xiàn),MIMO雷達在進行速度和位置的聯(lián)合估計時,隨著T的減小,發(fā)射信號帶寬增大,位置估計的誤差降低,而速度估計的誤差迅速提高,以至于在T=10-7s、SNR=10dB時速度估計的誤差遠遠大于能接受的范圍(如圖3所示),因此就必須增加發(fā)射脈沖寬度來使位置和速度聯(lián)合估計的誤差取到合適的值。由于在試驗中假設MIMO雷達各通道的信噪比為10dB,所以圖4所仿真的雷達參數(shù)是比較合適的,這個參數(shù)也是實際雷達可以很容易實現(xiàn)的。

    圖3 T=1.125 4×10-7s時參數(shù)估計的CRLB

    上一節(jié)提出的變限搜索方法,具體操作時所取的初始搜索范圍、搜索點數(shù)和范圍變化值會根據(jù)雷達監(jiān)視區(qū)域、對精度要求的差異而有所不同。初始的搜索范圍應該是雷達的監(jiān)視區(qū)域,在第1次捕捉到目標以后便實施跟蹤,利用跟蹤的預測值逐步縮小下一步最大似然估計的搜索點數(shù),直到趨于穩(wěn)定。在本實驗的仿真中,根據(jù)相應的精度要求,設定初始搜索范圍為[-100m,100m]×[-100m,100m],初始搜索點數(shù)為10 000。隨著跟蹤的進行,依據(jù)跟蹤誤差的變化規(guī)律,搜索范圍以10m的步長減小,并相應地減少最大似然參數(shù)估計的搜索點數(shù),直至搜索范圍為[-10m,10m]×[-10m,10m]。

    圖4 T=0.05s時參數(shù)估計的CRLB

    圖5描述了MIMO雷達交互式跟蹤中的變限參數(shù)估計法,具體的做法為:利用跟蹤濾波的預測值為目標下一時刻的參數(shù)估計設定一個置信范圍(圖中的矩形框),并將矩形離散成許多點,在最大似然估計時只需比較矩形內(nèi)離散點的似然比就可以得到該時刻目標參數(shù)的雷達測量值,極大地減輕了參數(shù)估計的運算量。由圖中可以看出從第2步跟蹤開始隨著跟蹤誤差的逐步減小,相應地縮小了MIMO雷達最大似然估計的搜索范圍(置信區(qū)域大小),參數(shù)估計的運算量銳減,提高了跟蹤和定位的速度。從圖6中可以發(fā)現(xiàn),目標的真實軌跡和跟蹤軌跡非常接近,說明MIMO雷達在考慮運動目標RCS閃爍的情況下仍然能夠?qū)δ繕吮3趾芎玫母櫺阅堋?/p>

    圖5 參數(shù)估計時的變限搜索

    圖7、圖8描述了MIMO雷達的跟蹤誤差、每一時刻的參數(shù)估計誤差和參數(shù)估計的CRLB,由誤差曲線圖分析發(fā)現(xiàn):由于在本文的算法中,目標的運動狀態(tài)矢量和MIMO雷達的測量參數(shù)是相同的,所以本文的跟蹤算法對雷達測量數(shù)據(jù)起到了很好的去噪作用。隨著跟蹤步數(shù)的推進,MIMO雷達的跟蹤誤差逐步減小,到后來降低至CRLB以下,當然也低于參數(shù)估計誤差。

    圖6 跟蹤結(jié)果和目標的真實軌跡

    圖7 X坐標的跟蹤濾波誤差和參數(shù)估計誤差

    圖8 Y坐標的跟蹤濾波誤差和參數(shù)估計誤差

    表1比較了傳統(tǒng)跟蹤方法與本文算法對目標實施跟蹤時每一步參數(shù)估計時的搜索點數(shù)和計算機的運行時間,直觀地說明本文方法對算法復雜度的改善和搜索效率的提高。

    表1 CPU運行時間比較

    5 結(jié)束語

    本文根據(jù)MIMO雷達特有的物理結(jié)構(gòu)和信號處理方式,對MIMO雷達的目標跟蹤技術(shù)進行了研究,提出了目標跟蹤與參數(shù)估計相結(jié)合的交互式跟蹤方法。經(jīng)過理論分析和仿真實驗結(jié)果發(fā)現(xiàn):文中所提的方法在目標RCS閃爍嚴重的情況下仍能很好地跟蹤目標,獲得很高的跟蹤精度,且不易丟失目標;同時,在交互式數(shù)據(jù)處理時所采用變限參數(shù)估計法可以有效地減輕MIMO雷達參數(shù)估計的運算復雜度,提高搜索效率,解決了參數(shù)估計中一大難題。而且,MIMO雷達跟蹤技術(shù)的實現(xiàn)為分布式MIMO雷達定向發(fā)射能量提供了理論支持,可以較好地提高雷達探測威力,這也是分布式MIMO雷達亟需須改進的地方。下一步需要繼續(xù)努力的方向是解決低信噪比、陣元數(shù)少的情況下的MIMO跟蹤問題。

    [1]Fisher E,Haimovich A,Blum R S,et al.MIMO radar:an idea whose time has come[A].IEEE Radar Conference[C].Philadelphia,2004:71-78.

    [2]Haimovich A,Blum R S,Lenard J,et al.MIMO radar with widely separated antennas[J].IEEE Transactions on Signal Processing,2008,25(1):116-129.

    [3]Li Jian,Stoica P.MIMO radar with collocated antennas-review of some recent work[J].IEEE Transactions on Signal Processing,2007,24(5):106-114.

    [4]肖文書.MIMO雷達中的信號檢測[J].電子學報,2010,38(3):626-631.

    [5]Ali Tajer,Jajamovich Guido H,Xiaodong Wang,et al.Optimal joint target detection and parameter estimation by MIMO radar[J].IEEE Transactions on Signal Processing,2010,4(1):127-145.

    [6]Qian He,Blum R S,Alexander M,Haimovich.Non-coherent MIMO radar for target estimation:more antennas means better performance[A].43rd CISS[C],2009:108-113.

    [7]Qian He,Blum R S,Hana Godrich,Haimovich Alexander M.Target velocity estimation and antenna placement for MIMO radar with widely separated antennas[J].IEEE Transactions on Signal Process,2010,4(1):79-99.

    [8]Ruixin Niu,Blum Rick S,Varshney Pramod K,et al.Target tracking in widely separated non-coherent multiple-input multiple-output radar systems[A].43rd ACSSC[C],2009:1181-1185.

    [9]Hana Godrich,Haimovich Alexander M,Blum Rick S.A MIMO radar system approach to target tracking[A].43rd ACSSC[C],2009:1186-1190.

    [10]Hana Godrich,Chiriac Vlad M,Haimovich Alexander M,Blum Rick S.Target tracking in MIMO radar systems:techniques and performance analysis[A].IEEE Radar Conference[C].2010:111-1116.

    [11]Fisher E,Haimovich A,Blum R S,et al.Spatial diversity in radars-models and detection performance[J].IEEE Transactions on Signal Processing,2006,54(3):823-838.

    猜你喜歡
    參數(shù)估計雷達噪聲
    有雷達
    大自然探索(2023年7期)2023-08-15 00:48:21
    基于新型DFrFT的LFM信號參數(shù)估計算法
    噪聲可退化且依賴于狀態(tài)和分布的平均場博弈
    雷達
    控制噪聲有妙法
    Logistic回歸模型的幾乎無偏兩參數(shù)估計
    基于向前方程的平穩(wěn)分布參數(shù)估計
    基于競爭失效數(shù)據(jù)的Lindley分布參數(shù)估計
    基于空時二維隨機輻射場的彈載雷達前視成像
    現(xiàn)代“千里眼”——雷達
    午夜视频国产福利| 亚洲精品日韩在线中文字幕| 国产成人freesex在线| 日韩欧美 国产精品| 中国美白少妇内射xxxbb| 高清在线视频一区二区三区| 久久久欧美国产精品| 天天一区二区日本电影三级| 中文天堂在线官网| 少妇人妻久久综合中文| 欧美性感艳星| 午夜爱爱视频在线播放| av女优亚洲男人天堂| 亚洲人与动物交配视频| 国产91av在线免费观看| 人人妻人人爽人人添夜夜欢视频 | 校园人妻丝袜中文字幕| 一本久久精品| 九九爱精品视频在线观看| 久热这里只有精品99| 亚洲人成网站在线观看播放| 22中文网久久字幕| 丰满人妻一区二区三区视频av| 国产视频首页在线观看| 中文乱码字字幕精品一区二区三区| 国产精品蜜桃在线观看| 人妻一区二区av| 夜夜看夜夜爽夜夜摸| 精品人妻视频免费看| 亚洲国产欧美在线一区| 久久6这里有精品| 在线免费十八禁| 内地一区二区视频在线| 精品少妇黑人巨大在线播放| 亚洲精品aⅴ在线观看| av国产免费在线观看| 欧美精品国产亚洲| 久久久久久久精品精品| 精品久久久久久久末码| 一级二级三级毛片免费看| 久久精品国产亚洲av天美| 嫩草影院入口| 亚洲美女视频黄频| 1000部很黄的大片| 欧美日韩综合久久久久久| 精品酒店卫生间| 亚洲国产日韩一区二区| 精品少妇黑人巨大在线播放| 国产又色又爽无遮挡免| 天天一区二区日本电影三级| 国产成人freesex在线| 在线精品无人区一区二区三 | 亚洲精品国产av蜜桃| 亚洲美女搞黄在线观看| 免费黄网站久久成人精品| av黄色大香蕉| 成人高潮视频无遮挡免费网站| 欧美老熟妇乱子伦牲交| 在线亚洲精品国产二区图片欧美 | 久久久久久久久久成人| av在线app专区| 成年av动漫网址| 欧美 日韩 精品 国产| 女的被弄到高潮叫床怎么办| 欧美丝袜亚洲另类| 久久精品夜色国产| 日日摸夜夜添夜夜爱| 久久久久久久精品精品| 国产免费视频播放在线视频| 校园人妻丝袜中文字幕| 青春草国产在线视频| 欧美一级a爱片免费观看看| 免费大片黄手机在线观看| 久久精品人妻少妇| 精品久久久久久久久av| 国语对白做爰xxxⅹ性视频网站| 国产人妻一区二区三区在| 色播亚洲综合网| 在线免费十八禁| av黄色大香蕉| 网址你懂的国产日韩在线| 成人黄色视频免费在线看| 久久久久久久午夜电影| 午夜福利视频1000在线观看| 中文字幕久久专区| 午夜福利高清视频| 欧美少妇被猛烈插入视频| 日本一二三区视频观看| 亚洲精品乱码久久久v下载方式| 人妻一区二区av| 日韩av不卡免费在线播放| 亚洲精品国产成人久久av| 18禁裸乳无遮挡免费网站照片| 国产欧美日韩一区二区三区在线 | 日本午夜av视频| 舔av片在线| 国产片特级美女逼逼视频| 国产日韩欧美亚洲二区| 亚洲经典国产精华液单| 国产国拍精品亚洲av在线观看| 97热精品久久久久久| 亚洲在线观看片| 亚洲欧美日韩东京热| 国产黄频视频在线观看| 嫩草影院精品99| 下体分泌物呈黄色| 精品久久久久久久人妻蜜臀av| 边亲边吃奶的免费视频| 一级毛片电影观看| 秋霞在线观看毛片| www.av在线官网国产| 精品人妻视频免费看| 男女边吃奶边做爰视频| 精品一区二区三区视频在线| 男人添女人高潮全过程视频| 纵有疾风起免费观看全集完整版| 国产男人的电影天堂91| 欧美日韩国产mv在线观看视频 | 黑人高潮一二区| 免费观看无遮挡的男女| 国产乱人偷精品视频| 五月开心婷婷网| 免费播放大片免费观看视频在线观看| 在现免费观看毛片| 美女被艹到高潮喷水动态| 亚洲激情五月婷婷啪啪| 久久综合国产亚洲精品| 嘟嘟电影网在线观看| 亚洲欧美日韩东京热| 欧美bdsm另类| 97热精品久久久久久| 直男gayav资源| 国产精品一区二区三区四区免费观看| 国产成人91sexporn| 免费观看无遮挡的男女| 观看免费一级毛片| 全区人妻精品视频| 国产成人福利小说| 欧美成人午夜免费资源| 一级毛片黄色毛片免费观看视频| 国产伦精品一区二区三区视频9| 18禁裸乳无遮挡免费网站照片| 97在线人人人人妻| 五月开心婷婷网| 日韩三级伦理在线观看| 亚洲国产精品成人久久小说| 一级毛片aaaaaa免费看小| 在线a可以看的网站| 亚洲国产欧美在线一区| 99热6这里只有精品| a级一级毛片免费在线观看| 午夜福利在线观看免费完整高清在| 精品久久国产蜜桃| 成人鲁丝片一二三区免费| 欧美老熟妇乱子伦牲交| 日本与韩国留学比较| 狂野欧美激情性xxxx在线观看| 亚洲欧美日韩另类电影网站 | 18禁动态无遮挡网站| 伦精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 搡老乐熟女国产| 少妇 在线观看| 可以在线观看毛片的网站| 免费av毛片视频| 国产成人91sexporn| 亚洲欧洲国产日韩| 最新中文字幕久久久久| 国产伦精品一区二区三区视频9| 国产av不卡久久| 一级毛片久久久久久久久女| 国模一区二区三区四区视频| 欧美 日韩 精品 国产| 天堂网av新在线| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在| 久久国产乱子免费精品| 欧美极品一区二区三区四区| 成年版毛片免费区| 成人国产av品久久久| 亚洲激情五月婷婷啪啪| 日韩强制内射视频| 国产一级毛片在线| 欧美激情国产日韩精品一区| 深夜a级毛片| 嫩草影院新地址| 国产精品蜜桃在线观看| 国产精品久久久久久久电影| 色吧在线观看| 色婷婷久久久亚洲欧美| 久久精品熟女亚洲av麻豆精品| 日本黄大片高清| 国产成人精品婷婷| 亚洲成人久久爱视频| 秋霞伦理黄片| 91午夜精品亚洲一区二区三区| 国国产精品蜜臀av免费| 麻豆乱淫一区二区| 亚洲在久久综合| 少妇被粗大猛烈的视频| 国产成人a∨麻豆精品| 内地一区二区视频在线| 三级国产精品欧美在线观看| 夫妻午夜视频| 国产69精品久久久久777片| 国产老妇女一区| 免费电影在线观看免费观看| 精品人妻熟女av久视频| 国产精品一区二区三区四区免费观看| 精品少妇黑人巨大在线播放| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产淫语在线视频| 国产亚洲5aaaaa淫片| 久久99热这里只频精品6学生| 欧美成人精品欧美一级黄| 日本-黄色视频高清免费观看| 国产成人精品婷婷| 美女视频免费永久观看网站| 欧美最新免费一区二区三区| 亚洲国产欧美在线一区| 婷婷色av中文字幕| 国产精品.久久久| 在线观看一区二区三区| 久久久色成人| 大又大粗又爽又黄少妇毛片口| 深夜a级毛片| 久久6这里有精品| 直男gayav资源| 人妻夜夜爽99麻豆av| 日本wwww免费看| 亚洲va在线va天堂va国产| 欧美一级a爱片免费观看看| 日韩人妻高清精品专区| 久久6这里有精品| av在线app专区| 国内精品美女久久久久久| 男女边吃奶边做爰视频| 狠狠精品人妻久久久久久综合| av天堂中文字幕网| 久久久久久国产a免费观看| 美女高潮的动态| 九草在线视频观看| www.色视频.com| 五月天丁香电影| 国产成人午夜福利电影在线观看| 欧美bdsm另类| 亚洲国产色片| 尾随美女入室| 国产精品一及| 亚洲国产色片| 成人免费观看视频高清| 一级a做视频免费观看| 建设人人有责人人尽责人人享有的 | 嫩草影院精品99| 午夜精品国产一区二区电影 | 男男h啪啪无遮挡| 成年女人在线观看亚洲视频 | 美女cb高潮喷水在线观看| 老司机影院成人| 中文字幕免费在线视频6| 七月丁香在线播放| 一级毛片久久久久久久久女| 成年人午夜在线观看视频| 欧美极品一区二区三区四区| 成年av动漫网址| 久久鲁丝午夜福利片| 制服丝袜香蕉在线| 男女边摸边吃奶| 哪个播放器可以免费观看大片| 日韩三级伦理在线观看| 毛片一级片免费看久久久久| 国产精品国产三级国产专区5o| 欧美高清性xxxxhd video| 小蜜桃在线观看免费完整版高清| 国产白丝娇喘喷水9色精品| 亚洲欧美一区二区三区黑人 | 午夜免费观看性视频| 久久精品国产亚洲av涩爱| 在线免费观看不下载黄p国产| 六月丁香七月| 欧美丝袜亚洲另类| 青春草国产在线视频| 亚洲在线观看片| 日本免费在线观看一区| 18禁裸乳无遮挡动漫免费视频 | 亚洲欧美日韩东京热| 自拍偷自拍亚洲精品老妇| 国产精品精品国产色婷婷| 99热全是精品| 在线播放无遮挡| 精品人妻视频免费看| 校园人妻丝袜中文字幕| 建设人人有责人人尽责人人享有的 | 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 中文字幕亚洲精品专区| 2021天堂中文幕一二区在线观| 亚洲av免费高清在线观看| www.av在线官网国产| 亚洲精品国产色婷婷电影| 新久久久久国产一级毛片| 中文欧美无线码| 男人爽女人下面视频在线观看| 欧美最新免费一区二区三区| 久久久久久久精品精品| 亚洲国产精品999| 啦啦啦在线观看免费高清www| 国产爽快片一区二区三区| 国产精品国产三级专区第一集| 在线观看美女被高潮喷水网站| 国产人妻一区二区三区在| 高清在线视频一区二区三区| 免费少妇av软件| 久久精品国产亚洲av涩爱| 亚洲自拍偷在线| 丰满乱子伦码专区| 午夜精品一区二区三区免费看| 亚洲美女视频黄频| 中文字幕制服av| 精品人妻偷拍中文字幕| 禁无遮挡网站| 在线a可以看的网站| 视频中文字幕在线观看| 特级一级黄色大片| 久久精品综合一区二区三区| 亚洲三级黄色毛片| 亚洲av福利一区| 别揉我奶头 嗯啊视频| 91aial.com中文字幕在线观看| 下体分泌物呈黄色| 亚洲精品aⅴ在线观看| 一级毛片我不卡| 五月开心婷婷网| 在线观看美女被高潮喷水网站| 午夜免费男女啪啪视频观看| 自拍欧美九色日韩亚洲蝌蚪91 | av在线播放精品| 成人毛片a级毛片在线播放| 大片电影免费在线观看免费| 亚洲欧美一区二区三区国产| 看黄色毛片网站| 人体艺术视频欧美日本| 26uuu在线亚洲综合色| 亚洲av福利一区| 亚洲最大成人手机在线| 不卡视频在线观看欧美| 一级毛片久久久久久久久女| 欧美老熟妇乱子伦牲交| 欧美性感艳星| 日韩欧美一区视频在线观看 | 亚洲色图av天堂| 精品少妇久久久久久888优播| 欧美高清成人免费视频www| 麻豆国产97在线/欧美| 99热这里只有是精品50| 亚洲精品色激情综合| 久久久久久久久久久丰满| 成人国产麻豆网| 嘟嘟电影网在线观看| 蜜臀久久99精品久久宅男| 黄色视频在线播放观看不卡| 1000部很黄的大片| 男人爽女人下面视频在线观看| 成人午夜精彩视频在线观看| 美女被艹到高潮喷水动态| 国产精品久久久久久久久免| 精品久久久久久久人妻蜜臀av| 亚洲精品视频女| 欧美日韩国产mv在线观看视频 | 51国产日韩欧美| 国产白丝娇喘喷水9色精品| 亚洲综合色惰| 欧美成人一区二区免费高清观看| 毛片一级片免费看久久久久| 亚洲精品第二区| 又黄又爽又刺激的免费视频.| 交换朋友夫妻互换小说| 校园人妻丝袜中文字幕| 久久久久久久久久久丰满| 在线看a的网站| 国产一区亚洲一区在线观看| 一本一本综合久久| 亚洲人成网站在线播| 黄色日韩在线| 亚洲天堂av无毛| 国产综合精华液| 精品久久久噜噜| 91午夜精品亚洲一区二区三区| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看| 香蕉精品网在线| 91狼人影院| 视频区图区小说| 精品久久久精品久久久| 国产91av在线免费观看| av网站免费在线观看视频| 黑人高潮一二区| 久久精品久久精品一区二区三区| 在线免费十八禁| 欧美激情久久久久久爽电影| av在线老鸭窝| 成年女人看的毛片在线观看| 国产欧美另类精品又又久久亚洲欧美| 69av精品久久久久久| 国产欧美日韩一区二区三区在线 | 看黄色毛片网站| 一级片'在线观看视频| 干丝袜人妻中文字幕| 22中文网久久字幕| 视频中文字幕在线观看| 日本欧美国产在线视频| 久久精品国产亚洲av涩爱| 午夜福利在线在线| 观看美女的网站| 亚洲av在线观看美女高潮| 国产毛片在线视频| 成人国产av品久久久| 久久久久久久国产电影| 免费av不卡在线播放| 丰满乱子伦码专区| 国产精品国产三级国产专区5o| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 蜜桃久久精品国产亚洲av| 天天躁夜夜躁狠狠久久av| 天美传媒精品一区二区| 国产大屁股一区二区在线视频| 婷婷色麻豆天堂久久| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 久久97久久精品| 亚洲性久久影院| 激情五月婷婷亚洲| 免费人成在线观看视频色| 2018国产大陆天天弄谢| 97在线人人人人妻| 亚洲av成人精品一二三区| 亚洲美女搞黄在线观看| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 亚洲av国产av综合av卡| 久久精品夜色国产| 日韩人妻高清精品专区| 国产淫片久久久久久久久| 精品熟女少妇av免费看| 丰满少妇做爰视频| 亚洲最大成人中文| 激情 狠狠 欧美| 欧美97在线视频| 中文精品一卡2卡3卡4更新| 亚洲在线观看片| 草草在线视频免费看| 日本欧美国产在线视频| 免费观看在线日韩| 一本一本综合久久| 国内精品美女久久久久久| 99久久精品国产国产毛片| 日本-黄色视频高清免费观看| 成人无遮挡网站| 国产精品99久久久久久久久| 精品久久久久久久久av| 成人欧美大片| 有码 亚洲区| 亚洲av中文字字幕乱码综合| 大码成人一级视频| 成人美女网站在线观看视频| 97超视频在线观看视频| 久久久久久久久久久丰满| 精品久久久久久电影网| 亚洲久久久久久中文字幕| 又大又黄又爽视频免费| 国产爽快片一区二区三区| 久久99蜜桃精品久久| 秋霞在线观看毛片| 亚洲精品自拍成人| 天天躁日日操中文字幕| 三级经典国产精品| 久久97久久精品| 色5月婷婷丁香| 欧美+日韩+精品| 只有这里有精品99| 成年人午夜在线观看视频| 精品国产露脸久久av麻豆| 精品久久久久久久末码| 熟女电影av网| 69人妻影院| 91狼人影院| 大片电影免费在线观看免费| 精品一区二区三卡| 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 免费看a级黄色片| 乱系列少妇在线播放| 交换朋友夫妻互换小说| 国产 精品1| 日韩一本色道免费dvd| 亚洲欧美一区二区三区国产| av网站免费在线观看视频| 99久久精品国产国产毛片| 夜夜看夜夜爽夜夜摸| 国产在视频线精品| 色综合色国产| 91精品伊人久久大香线蕉| 欧美日韩综合久久久久久| 亚洲av免费高清在线观看| 纵有疾风起免费观看全集完整版| 久久久色成人| 永久免费av网站大全| 日韩 亚洲 欧美在线| 欧美亚洲 丝袜 人妻 在线| 三级国产精品片| 久久久午夜欧美精品| 狂野欧美白嫩少妇大欣赏| 久久久久久九九精品二区国产| 男人和女人高潮做爰伦理| 成人亚洲欧美一区二区av| 赤兔流量卡办理| 国产av不卡久久| 免费少妇av软件| 各种免费的搞黄视频| 男女那种视频在线观看| 中国三级夫妇交换| 中文在线观看免费www的网站| 视频区图区小说| 高清在线视频一区二区三区| 国产成人免费观看mmmm| 老女人水多毛片| 日韩av免费高清视频| 成人欧美大片| 亚洲精品国产成人久久av| 欧美人与善性xxx| 欧美日韩视频高清一区二区三区二| 亚洲最大成人中文| 卡戴珊不雅视频在线播放| av天堂中文字幕网| 禁无遮挡网站| 高清午夜精品一区二区三区| 亚洲国产色片| 亚洲激情五月婷婷啪啪| 国产精品久久久久久精品古装| 少妇被粗大猛烈的视频| 国产在线一区二区三区精| 精品一区二区免费观看| 国产亚洲av片在线观看秒播厂| 欧美精品国产亚洲| 精品久久久久久久人妻蜜臀av| 噜噜噜噜噜久久久久久91| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 免费av不卡在线播放| 国产精品嫩草影院av在线观看| 午夜日本视频在线| 极品教师在线视频| 久久久成人免费电影| 日日啪夜夜爽| 日韩大片免费观看网站| 中文字幕久久专区| 久久久久久久午夜电影| 日韩在线高清观看一区二区三区| 亚洲欧洲日产国产| 久久久久久久久久成人| 女人十人毛片免费观看3o分钟| 网址你懂的国产日韩在线| 亚洲综合精品二区| 天堂中文最新版在线下载 | 国产一级毛片在线| 精华霜和精华液先用哪个| 亚洲久久久久久中文字幕| 中国三级夫妇交换| 国产伦精品一区二区三区四那| 国产高清国产精品国产三级 | av在线亚洲专区| 一级二级三级毛片免费看| 特级一级黄色大片| 国产淫语在线视频| 免费少妇av软件| 丰满乱子伦码专区| 午夜免费男女啪啪视频观看| 国产成人91sexporn| 国产在线男女| 色播亚洲综合网| 日本-黄色视频高清免费观看| 久久久色成人| 免费高清在线观看视频在线观看| 久久久久久久大尺度免费视频| 国产欧美另类精品又又久久亚洲欧美| 国产成人一区二区在线| 亚洲精品国产av成人精品| 亚洲欧美日韩无卡精品| 亚洲三级黄色毛片| 国产精品久久久久久av不卡| 亚洲精品久久久久久婷婷小说| 尤物成人国产欧美一区二区三区| 在线观看国产h片| 夜夜看夜夜爽夜夜摸| 18禁动态无遮挡网站| 一级a做视频免费观看| 三级男女做爰猛烈吃奶摸视频| 热re99久久精品国产66热6| 自拍欧美九色日韩亚洲蝌蚪91 | 日本黄色片子视频| 欧美精品国产亚洲| 波野结衣二区三区在线| 欧美丝袜亚洲另类| 国产成人精品婷婷| av天堂中文字幕网| 国产精品伦人一区二区| 国国产精品蜜臀av免费| 99精国产麻豆久久婷婷| 免费看日本二区| av一本久久久久| 一级毛片黄色毛片免费观看视频| 中文字幕制服av| 91久久精品国产一区二区三区|