• <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)代“千里眼”——雷達
    一级毛片精品| av免费在线观看网站| 国产亚洲欧美98| 亚洲成av片中文字幕在线观看| 亚洲欧美精品综合一区二区三区| 亚洲综合色网址| 免费人成视频x8x8入口观看| 国产一区在线观看成人免费| 悠悠久久av| 成熟少妇高潮喷水视频| 国产伦人伦偷精品视频| 日韩欧美一区视频在线观看| 少妇猛男粗大的猛烈进出视频| 国产精品美女特级片免费视频播放器 | 女人高潮潮喷娇喘18禁视频| 一级毛片精品| 中文字幕人妻熟女乱码| 老司机影院毛片| 亚洲国产精品sss在线观看 | 亚洲专区国产一区二区| 日本vs欧美在线观看视频| 久久人妻福利社区极品人妻图片| 亚洲av成人不卡在线观看播放网| 色综合婷婷激情| 国产单亲对白刺激| 国产亚洲精品久久久久5区| 亚洲第一av免费看| 免费在线观看视频国产中文字幕亚洲| 亚洲精品av麻豆狂野| 欧美成人免费av一区二区三区 | 亚洲精品久久成人aⅴ小说| 国产精品综合久久久久久久免费 | 亚洲全国av大片| 成年人午夜在线观看视频| 日本撒尿小便嘘嘘汇集6| 91成年电影在线观看| 美女午夜性视频免费| xxx96com| 视频在线观看一区二区三区| 亚洲成人国产一区在线观看| 国产精品综合久久久久久久免费 | 精品人妻1区二区| 成人特级黄色片久久久久久久| 免费一级毛片在线播放高清视频 | 老汉色∧v一级毛片| 香蕉久久夜色| 国产精品免费大片| tocl精华| 国产伦人伦偷精品视频| 高清视频免费观看一区二区| 99国产极品粉嫩在线观看| 亚洲自偷自拍图片 自拍| 亚洲精品国产色婷婷电影| 国产极品粉嫩免费观看在线| 51午夜福利影视在线观看| 99re6热这里在线精品视频| av线在线观看网站| 美女国产高潮福利片在线看| 人妻丰满熟妇av一区二区三区 | 免费一级毛片在线播放高清视频 | 伊人久久大香线蕉亚洲五| 欧美 亚洲 国产 日韩一| 我的亚洲天堂| 91麻豆av在线| 色综合欧美亚洲国产小说| 午夜老司机福利片| 国产黄色免费在线视频| 香蕉丝袜av| 天堂动漫精品| 午夜福利视频在线观看免费| 欧美成人午夜精品| 丝袜美足系列| 男男h啪啪无遮挡| 久久精品亚洲熟妇少妇任你| 国产精品永久免费网站| 无人区码免费观看不卡| 亚洲五月天丁香| 99在线人妻在线中文字幕 | 无人区码免费观看不卡| 亚洲国产欧美一区二区综合| 成人免费观看视频高清| 国产免费av片在线观看野外av| 亚洲成人免费电影在线观看| 好看av亚洲va欧美ⅴa在| xxx96com| av网站在线播放免费| bbb黄色大片| 精品免费久久久久久久清纯 | 国产激情久久老熟女| 制服人妻中文乱码| 人成视频在线观看免费观看| 中文字幕人妻丝袜一区二区| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲高清精品| 国产精品免费视频内射| 丁香六月欧美| 熟女少妇亚洲综合色aaa.| 欧美乱妇无乱码| 亚洲精品中文字幕一二三四区| 国产午夜精品久久久久久| 一级a爱片免费观看的视频| 国产精品亚洲一级av第二区| 亚洲av片天天在线观看| 久久久久久久久免费视频了| 久久精品国产99精品国产亚洲性色 | 精品亚洲成a人片在线观看| 一边摸一边抽搐一进一出视频| 我的亚洲天堂| 欧美日韩亚洲国产一区二区在线观看 | 90打野战视频偷拍视频| 国产主播在线观看一区二区| 中文字幕制服av| 成年人黄色毛片网站| 欧美日韩福利视频一区二区| 免费高清在线观看日韩| 大片电影免费在线观看免费| 99国产精品一区二区蜜桃av | 亚洲欧洲精品一区二区精品久久久| 三级毛片av免费| 99香蕉大伊视频| 久久午夜亚洲精品久久| 18在线观看网站| 精品福利永久在线观看| 视频区欧美日本亚洲| 精品人妻熟女毛片av久久网站| 久久国产乱子伦精品免费另类| 成人av一区二区三区在线看| 午夜两性在线视频| 国产精品98久久久久久宅男小说| 窝窝影院91人妻| av一本久久久久| 国产成人欧美在线观看 | 国产欧美日韩一区二区三区在线| 欧美成人午夜精品| 老熟妇乱子伦视频在线观看| av不卡在线播放| 日日摸夜夜添夜夜添小说| av国产精品久久久久影院| 欧美久久黑人一区二区| av国产精品久久久久影院| 淫妇啪啪啪对白视频| 丰满迷人的少妇在线观看| 黄片播放在线免费| 最近最新中文字幕大全电影3 | 少妇裸体淫交视频免费看高清 | 久久九九热精品免费| av超薄肉色丝袜交足视频| 啪啪无遮挡十八禁网站| 久久亚洲真实| 18禁裸乳无遮挡免费网站照片 | 亚洲精品久久成人aⅴ小说| 一边摸一边做爽爽视频免费| 欧美人与性动交α欧美软件| 天天躁日日躁夜夜躁夜夜| 99久久99久久久精品蜜桃| 王馨瑶露胸无遮挡在线观看| 丁香六月欧美| 一级毛片高清免费大全| 在线国产一区二区在线| 在线观看舔阴道视频| 精品国产国语对白av| 黑人欧美特级aaaaaa片| 欧美一级毛片孕妇| 国产极品粉嫩免费观看在线| 大香蕉久久成人网| 国产精品久久久久久人妻精品电影| 欧美日韩中文字幕国产精品一区二区三区 | 久久久精品免费免费高清| 精品亚洲成国产av| 欧美 日韩 精品 国产| 成人特级黄色片久久久久久久| 亚洲精品成人av观看孕妇| 精品久久久精品久久久| 精品高清国产在线一区| 国产日韩一区二区三区精品不卡| 国产精品 欧美亚洲| 精品久久久久久电影网| 丝袜人妻中文字幕| 一级片免费观看大全| 久久人妻福利社区极品人妻图片| 一二三四社区在线视频社区8| 日韩欧美一区视频在线观看| 两个人免费观看高清视频| 亚洲精品在线观看二区| 中文字幕制服av| 99久久人妻综合| 精品电影一区二区在线| 免费在线观看影片大全网站| 老汉色∧v一级毛片| 91字幕亚洲| 国产深夜福利视频在线观看| 男女午夜视频在线观看| 欧美精品高潮呻吟av久久| 国产高清激情床上av| 露出奶头的视频| 久久久精品国产亚洲av高清涩受| aaaaa片日本免费| 久久草成人影院| 免费久久久久久久精品成人欧美视频| 国产精品永久免费网站| 国内久久婷婷六月综合欲色啪| 久久国产精品人妻蜜桃| 最新美女视频免费是黄的| 两个人看的免费小视频| 看免费av毛片| 午夜福利影视在线免费观看| 又黄又爽又免费观看的视频| 变态另类成人亚洲欧美熟女 | 夜夜爽天天搞| 麻豆国产av国片精品| 99热网站在线观看| www.精华液| 亚洲精品久久午夜乱码| 看片在线看免费视频| 国产成+人综合+亚洲专区| 亚洲专区国产一区二区| 男人舔女人的私密视频| 国产一区在线观看成人免费| 亚洲成av片中文字幕在线观看| 久久久国产精品麻豆| 精品福利永久在线观看| 国产精品偷伦视频观看了| 极品人妻少妇av视频| 欧洲精品卡2卡3卡4卡5卡区| 啦啦啦视频在线资源免费观看| 老司机在亚洲福利影院| 12—13女人毛片做爰片一| 女人高潮潮喷娇喘18禁视频| 两个人免费观看高清视频| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区久久久樱花| 一级片免费观看大全| 欧美精品一区二区免费开放| 中文欧美无线码| 国产精品亚洲av一区麻豆| 国产亚洲欧美精品永久| 亚洲欧美日韩高清在线视频| 色综合婷婷激情| 一区二区三区激情视频| 丝瓜视频免费看黄片| 国产又爽黄色视频| 欧美激情极品国产一区二区三区| 亚洲av日韩在线播放| 亚洲人成电影免费在线| 91麻豆av在线| 精品少妇久久久久久888优播| 91在线观看av| 变态另类成人亚洲欧美熟女 | 91国产中文字幕| av超薄肉色丝袜交足视频| 欧美成狂野欧美在线观看| 视频在线观看一区二区三区| 欧美日韩福利视频一区二区| 亚洲第一欧美日韩一区二区三区| www.熟女人妻精品国产| 免费一级毛片在线播放高清视频 | 美女高潮到喷水免费观看| 狠狠婷婷综合久久久久久88av| 日韩人妻精品一区2区三区| 一级作爱视频免费观看| 老司机福利观看| 成人影院久久| 久久精品亚洲精品国产色婷小说| 亚洲精品一卡2卡三卡4卡5卡| 久久人人97超碰香蕉20202| 在线观看免费视频网站a站| 黄色视频,在线免费观看| av线在线观看网站| 女人高潮潮喷娇喘18禁视频| 啦啦啦免费观看视频1| 91麻豆精品激情在线观看国产 | 久久中文看片网| 成人免费观看视频高清| 午夜福利在线观看吧| 国产97色在线日韩免费| 成人18禁高潮啪啪吃奶动态图| 亚洲av片天天在线观看| 97人妻天天添夜夜摸| 午夜福利视频在线观看免费| 中文字幕人妻丝袜一区二区| 99精国产麻豆久久婷婷| 欧美日韩视频精品一区| 亚洲免费av在线视频| 欧美大码av| 欧美日韩国产mv在线观看视频| 国产乱人伦免费视频| www.自偷自拍.com| 国产精品久久久人人做人人爽| 国产午夜精品久久久久久| 亚洲欧美一区二区三区黑人| 青草久久国产| 精品午夜福利视频在线观看一区| 国产激情欧美一区二区| 色婷婷久久久亚洲欧美| 日本黄色视频三级网站网址 | 日本vs欧美在线观看视频| 在线播放国产精品三级| 精品久久久精品久久久| 水蜜桃什么品种好| 精品国产亚洲在线| 日韩免费av在线播放| 91精品三级在线观看| 视频在线观看一区二区三区| 大码成人一级视频| 人人妻人人添人人爽欧美一区卜| 老熟妇乱子伦视频在线观看| 日本五十路高清| 99热只有精品国产| 色精品久久人妻99蜜桃| 热99re8久久精品国产| 亚洲欧美日韩高清在线视频| 亚洲五月色婷婷综合| 亚洲精品在线观看二区| 欧美日本中文国产一区发布| cao死你这个sao货| 国产亚洲欧美98| 一级毛片精品| 亚洲精品在线美女| 国产免费av片在线观看野外av| 黑丝袜美女国产一区| 国产精品偷伦视频观看了| 丰满人妻熟妇乱又伦精品不卡| 国产在视频线精品| 午夜影院日韩av| 成年人免费黄色播放视频| 欧美亚洲 丝袜 人妻 在线| 首页视频小说图片口味搜索| 一级毛片精品| 一边摸一边抽搐一进一出视频| 午夜日韩欧美国产| e午夜精品久久久久久久| 女性生殖器流出的白浆| 天天躁狠狠躁夜夜躁狠狠躁| 97人妻天天添夜夜摸| av一本久久久久| 亚洲av成人一区二区三| 国产亚洲一区二区精品| 黄色女人牲交| 夜夜爽天天搞| www.精华液| 一二三四在线观看免费中文在| 精品一区二区三卡| 少妇裸体淫交视频免费看高清 | 午夜日韩欧美国产| 成人av一区二区三区在线看| 国产一区二区三区综合在线观看| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 久久香蕉精品热| 欧美黑人精品巨大| 国产精品 国内视频| 国产成人啪精品午夜网站| 成人永久免费在线观看视频| 村上凉子中文字幕在线| 天天添夜夜摸| 三级毛片av免费| 757午夜福利合集在线观看| 国产精品久久久av美女十八| 国产男女超爽视频在线观看| 很黄的视频免费| 亚洲成人手机| 国产精品影院久久| 国内毛片毛片毛片毛片毛片| 免费不卡黄色视频| 亚洲成人手机| 波多野结衣av一区二区av| 久久精品人人爽人人爽视色| 精品一区二区三区四区五区乱码| 一级毛片精品| 精品电影一区二区在线| 国产三级黄色录像| 极品人妻少妇av视频| 91九色精品人成在线观看| 天天影视国产精品| 久久人妻熟女aⅴ| 亚洲一区高清亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 国产又爽黄色视频| av超薄肉色丝袜交足视频| av一本久久久久| 午夜视频精品福利| 日韩视频一区二区在线观看| 首页视频小说图片口味搜索| 中文字幕色久视频| 亚洲成人手机| 精品人妻1区二区| 亚洲av第一区精品v没综合| 夫妻午夜视频| 天堂俺去俺来也www色官网| 天堂√8在线中文| 亚洲人成电影免费在线| 欧美乱码精品一区二区三区| 国产一区有黄有色的免费视频| 亚洲欧美一区二区三区黑人| 大片电影免费在线观看免费| 国产乱人伦免费视频| 高清视频免费观看一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 王馨瑶露胸无遮挡在线观看| 精品久久久久久,| 视频在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 午夜免费成人在线视频| 俄罗斯特黄特色一大片| 美女扒开内裤让男人捅视频| 国产亚洲精品久久久久久毛片 | 久久久国产成人精品二区 | 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 中亚洲国语对白在线视频| 国产在线精品亚洲第一网站| 欧美+亚洲+日韩+国产| 少妇裸体淫交视频免费看高清 | 精品一区二区三区视频在线观看免费 | 亚洲国产中文字幕在线视频| 日韩中文字幕欧美一区二区| 露出奶头的视频| 精品久久蜜臀av无| 亚洲欧美精品综合一区二区三区| 美女高潮喷水抽搐中文字幕| 日本wwww免费看| 国产精品99久久99久久久不卡| 99精品欧美一区二区三区四区| 很黄的视频免费| 九色亚洲精品在线播放| 捣出白浆h1v1| 国产区一区二久久| 一级黄色大片毛片| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 黄色 视频免费看| 丝袜美腿诱惑在线| 国产人伦9x9x在线观看| 亚洲伊人色综图| 国产精品久久久久久人妻精品电影| 久久影院123| 18禁观看日本| 交换朋友夫妻互换小说| 精品第一国产精品| av天堂在线播放| 最近最新免费中文字幕在线| 精品第一国产精品| 亚洲国产精品sss在线观看 | 嫩草影视91久久| 黄色视频,在线免费观看| 国产成人欧美| 热99久久久久精品小说推荐| 18禁黄网站禁片午夜丰满| 高清欧美精品videossex| 国产蜜桃级精品一区二区三区 | 亚洲欧美激情在线| 露出奶头的视频| 法律面前人人平等表现在哪些方面| 老熟妇仑乱视频hdxx| 男女免费视频国产| 脱女人内裤的视频| 久久国产精品男人的天堂亚洲| 国产精品久久视频播放| 看黄色毛片网站| 热99久久久久精品小说推荐| 很黄的视频免费| 少妇 在线观看| 俄罗斯特黄特色一大片| 9色porny在线观看| 精品久久久久久,| tocl精华| 在线十欧美十亚洲十日本专区| 在线免费观看的www视频| 大码成人一级视频| 老熟妇仑乱视频hdxx| 亚洲五月天丁香| 亚洲一区二区三区不卡视频| 天堂动漫精品| 国产亚洲欧美精品永久| 成人国产一区最新在线观看| 亚洲精品自拍成人| 好看av亚洲va欧美ⅴa在| 日本黄色日本黄色录像| 搡老乐熟女国产| 国产在线观看jvid| 操出白浆在线播放| 色老头精品视频在线观看| 女性生殖器流出的白浆| 欧美中文综合在线视频| 精品人妻1区二区| 亚洲人成电影观看| 国产精品乱码一区二三区的特点 | 国产精品一区二区免费欧美| 国产真人三级小视频在线观看| 国产极品粉嫩免费观看在线| 又紧又爽又黄一区二区| 中文字幕另类日韩欧美亚洲嫩草| 国产精品电影一区二区三区 | 老司机午夜福利在线观看视频| 老熟妇仑乱视频hdxx| 欧美日韩成人在线一区二区| 色婷婷久久久亚洲欧美| 午夜福利免费观看在线| 我的亚洲天堂| 国产欧美日韩一区二区三| 狂野欧美激情性xxxx| 国产精品久久久久久精品古装| 久久久久久久久久久久大奶| 欧美日韩国产mv在线观看视频| 精品人妻熟女毛片av久久网站| 免费av中文字幕在线| 在线观看一区二区三区激情| 操出白浆在线播放| 高清av免费在线| 丁香欧美五月| 成年女人毛片免费观看观看9 | 狠狠狠狠99中文字幕| 热99国产精品久久久久久7| 免费观看人在逋| 国产在视频线精品| 叶爱在线成人免费视频播放| 99国产精品99久久久久| 老汉色av国产亚洲站长工具| 亚洲情色 制服丝袜| 国产精品免费一区二区三区在线 | 妹子高潮喷水视频| 桃红色精品国产亚洲av| 中文亚洲av片在线观看爽 | 中文亚洲av片在线观看爽 | 亚洲成人免费电影在线观看| 亚洲人成电影观看| 热re99久久精品国产66热6| 欧美色视频一区免费| 一进一出抽搐gif免费好疼 | 午夜福利影视在线免费观看| 亚洲五月婷婷丁香| 亚洲国产毛片av蜜桃av| 精品国产乱子伦一区二区三区| 精品免费久久久久久久清纯 | 亚洲精品乱久久久久久| 日本精品一区二区三区蜜桃| 久久精品亚洲熟妇少妇任你| 免费观看a级毛片全部| 巨乳人妻的诱惑在线观看| 久久青草综合色| 男女下面插进去视频免费观看| 999久久久精品免费观看国产| 久久久久视频综合| 日本黄色日本黄色录像| 激情在线观看视频在线高清 | 精品视频人人做人人爽| 看片在线看免费视频| 亚洲精华国产精华精| 乱人伦中国视频| 正在播放国产对白刺激| 叶爱在线成人免费视频播放| 精品一区二区三区视频在线观看免费 | 一级片免费观看大全| 日韩制服丝袜自拍偷拍| 国产高清videossex| 亚洲熟妇中文字幕五十中出 | 老汉色∧v一级毛片| 成人国语在线视频| 久热这里只有精品99| 久久久久久久久免费视频了| 丝袜美腿诱惑在线| 999精品在线视频| 亚洲avbb在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲五月色婷婷综合| 十分钟在线观看高清视频www| 男女床上黄色一级片免费看| 欧美激情久久久久久爽电影 | 亚洲熟女精品中文字幕| 一a级毛片在线观看| 美女福利国产在线| 18禁裸乳无遮挡免费网站照片 | 久久这里只有精品19| 亚洲精品一二三| 高清欧美精品videossex| 一区在线观看完整版| 国产高清激情床上av| 亚洲欧美一区二区三区久久| 久久精品亚洲av国产电影网| 国产一区在线观看成人免费| 成人精品一区二区免费| 欧美日韩成人在线一区二区| 久久国产精品男人的天堂亚洲| 免费一级毛片在线播放高清视频 | 久久99一区二区三区| 亚洲 国产 在线| 制服人妻中文乱码| 中文字幕人妻丝袜制服| 欧美亚洲日本最大视频资源| 欧美中文综合在线视频| 黄色视频,在线免费观看| 妹子高潮喷水视频| 国产男女内射视频| 极品教师在线免费播放| 免费在线观看影片大全网站| 国产精品 欧美亚洲| 国产精品av久久久久免费| 久久久久久免费高清国产稀缺| 国产亚洲欧美精品永久| 亚洲av电影在线进入| 国产99白浆流出| 欧美中文综合在线视频| 叶爱在线成人免费视频播放| 每晚都被弄得嗷嗷叫到高潮| 黄色丝袜av网址大全| 男女之事视频高清在线观看| 91字幕亚洲| 国产在线精品亚洲第一网站| 亚洲黑人精品在线| 国产亚洲精品第一综合不卡| 亚洲 欧美一区二区三区| 这个男人来自地球电影免费观看| 天天影视国产精品| 在线国产一区二区在线| 亚洲色图 男人天堂 中文字幕|