• <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)代“千里眼”——雷達
    一区二区三区四区激情视频| 国产精品国产三级国产专区5o| 婷婷色av中文字幕| 国产91精品成人一区二区三区 | 欧美精品高潮呻吟av久久| 成年人免费黄色播放视频| 美女大奶头黄色视频| 午夜老司机福利片| 黑人欧美特级aaaaaa片| 另类精品久久| 国产主播在线观看一区二区 | 久久国产精品影院| 人妻 亚洲 视频| 99国产综合亚洲精品| 久久精品久久久久久噜噜老黄| 国产一区二区 视频在线| 午夜激情av网站| 免费少妇av软件| 日韩熟女老妇一区二区性免费视频| 只有这里有精品99| 欧美中文综合在线视频| 亚洲欧美一区二区三区国产| 两个人免费观看高清视频| 你懂的网址亚洲精品在线观看| 中文字幕人妻熟女乱码| 啦啦啦中文免费视频观看日本| 又紧又爽又黄一区二区| 考比视频在线观看| 免费看十八禁软件| 中文字幕最新亚洲高清| 亚洲精品成人av观看孕妇| 国产精品香港三级国产av潘金莲 | 日韩 亚洲 欧美在线| 亚洲第一av免费看| 在线看a的网站| 最新在线观看一区二区三区 | 精品少妇黑人巨大在线播放| 欧美日韩视频精品一区| 啦啦啦中文免费视频观看日本| 久久精品久久久久久久性| 18禁国产床啪视频网站| 欧美成狂野欧美在线观看| 中文字幕最新亚洲高清| 美女视频免费永久观看网站| 日本黄色日本黄色录像| 黄色片一级片一级黄色片| 中文字幕av电影在线播放| 国产成人欧美| 亚洲成色77777| 最黄视频免费看| 国产成人av激情在线播放| 精品国产超薄肉色丝袜足j| 亚洲精品一二三| 中国美女看黄片| 欧美少妇被猛烈插入视频| 久久久久网色| 欧美变态另类bdsm刘玥| 亚洲五月婷婷丁香| 黄色怎么调成土黄色| 久久久久久免费高清国产稀缺| 亚洲第一青青草原| 最黄视频免费看| 国产成人免费无遮挡视频| 国产一级毛片在线| 黑丝袜美女国产一区| 精品久久久久久电影网| 美国免费a级毛片| 国产亚洲av片在线观看秒播厂| 夫妻性生交免费视频一级片| 久久久久国产一级毛片高清牌| 热99国产精品久久久久久7| 巨乳人妻的诱惑在线观看| 人体艺术视频欧美日本| 一本—道久久a久久精品蜜桃钙片| 午夜91福利影院| 免费黄频网站在线观看国产| 人成视频在线观看免费观看| 久久国产精品人妻蜜桃| 亚洲 国产 在线| 人人妻人人澡人人爽人人夜夜| 日韩制服骚丝袜av| 美女国产高潮福利片在线看| 1024香蕉在线观看| 久久精品久久久久久噜噜老黄| 热re99久久国产66热| 黑人巨大精品欧美一区二区蜜桃| 51午夜福利影视在线观看| xxx大片免费视频| 国产老妇伦熟女老妇高清| 亚洲精品久久成人aⅴ小说| 亚洲av欧美aⅴ国产| 欧美另类一区| 国产亚洲av高清不卡| 欧美黑人欧美精品刺激| 免费观看a级毛片全部| 黄色毛片三级朝国网站| 亚洲国产毛片av蜜桃av| 国产黄色免费在线视频| 国产成人免费无遮挡视频| 乱人伦中国视频| 叶爱在线成人免费视频播放| 777米奇影视久久| 亚洲午夜精品一区,二区,三区| 亚洲精品久久成人aⅴ小说| 婷婷色麻豆天堂久久| 欧美久久黑人一区二区| 纯流量卡能插随身wifi吗| 欧美黑人欧美精品刺激| 日本a在线网址| 国产黄色免费在线视频| 久9热在线精品视频| 国产激情久久老熟女| 女性被躁到高潮视频| 精品人妻在线不人妻| 日韩av不卡免费在线播放| 成人18禁高潮啪啪吃奶动态图| 超碰97精品在线观看| 欧美激情极品国产一区二区三区| 99久久精品国产亚洲精品| 老司机深夜福利视频在线观看 | av网站免费在线观看视频| 搡老乐熟女国产| 激情视频va一区二区三区| 欧美日本中文国产一区发布| 爱豆传媒免费全集在线观看| 97精品久久久久久久久久精品| 成年人黄色毛片网站| 黑人欧美特级aaaaaa片| 男人舔女人的私密视频| 亚洲av欧美aⅴ国产| 两个人看的免费小视频| 少妇粗大呻吟视频| 国产高清不卡午夜福利| 成人国产av品久久久| 男女高潮啪啪啪动态图| 丝瓜视频免费看黄片| 免费观看人在逋| 夫妻午夜视频| 丰满饥渴人妻一区二区三| 国产欧美日韩一区二区三 | 十分钟在线观看高清视频www| 欧美性长视频在线观看| 亚洲成人国产一区在线观看 | 国产深夜福利视频在线观看| 纵有疾风起免费观看全集完整版| 精品一区在线观看国产| 精品亚洲乱码少妇综合久久| 最新的欧美精品一区二区| 亚洲色图综合在线观看| 精品一区二区三卡| 美女主播在线视频| 男人爽女人下面视频在线观看| 亚洲,欧美精品.| 如日韩欧美国产精品一区二区三区| 丁香六月欧美| 免费看不卡的av| 如日韩欧美国产精品一区二区三区| 精品国产国语对白av| 精品一区二区三卡| 亚洲欧美成人综合另类久久久| 欧美日韩国产mv在线观看视频| 天天影视国产精品| 成年人黄色毛片网站| 亚洲av片天天在线观看| 国产精品久久久人人做人人爽| 欧美日韩视频精品一区| 在线 av 中文字幕| videosex国产| 亚洲综合色网址| 侵犯人妻中文字幕一二三四区| 捣出白浆h1v1| 一边亲一边摸免费视频| 妹子高潮喷水视频| 男女之事视频高清在线观看 | 日本欧美视频一区| 观看av在线不卡| 免费在线观看视频国产中文字幕亚洲 | 自线自在国产av| 国产免费福利视频在线观看| 欧美日韩视频高清一区二区三区二| 乱人伦中国视频| 一级片免费观看大全| 青青草视频在线视频观看| 十八禁高潮呻吟视频| av网站免费在线观看视频| 久久久久久久精品精品| 国产真人三级小视频在线观看| 久久精品成人免费网站| 久久精品久久精品一区二区三区| 人妻 亚洲 视频| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 80岁老熟妇乱子伦牲交| 亚洲av成人精品一二三区| 久久精品亚洲熟妇少妇任你| 国产免费现黄频在线看| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 精品福利永久在线观看| 亚洲欧美一区二区三区久久| 在线精品无人区一区二区三| 一级毛片 在线播放| 亚洲精品第二区| 岛国毛片在线播放| 国产精品一区二区在线不卡| 极品少妇高潮喷水抽搐| 两性夫妻黄色片| 9191精品国产免费久久| 精品少妇黑人巨大在线播放| 亚洲欧洲日产国产| 波多野结衣av一区二区av| 两性夫妻黄色片| 亚洲欧美清纯卡通| 久久鲁丝午夜福利片| 在线av久久热| 黑人巨大精品欧美一区二区蜜桃| 久久免费观看电影| 免费在线观看视频国产中文字幕亚洲 | 国产精品九九99| 91麻豆av在线| 久久久久久久国产电影| 性色av乱码一区二区三区2| 黑人欧美特级aaaaaa片| 午夜影院在线不卡| 精品国产乱码久久久久久小说| 男人爽女人下面视频在线观看| 日韩,欧美,国产一区二区三区| 国产男女内射视频| 中文字幕制服av| 欧美在线一区亚洲| 99久久99久久久精品蜜桃| 91麻豆精品激情在线观看国产 | 人妻人人澡人人爽人人| 一二三四社区在线视频社区8| 欧美黄色淫秽网站| 亚洲专区中文字幕在线| 1024视频免费在线观看| 老司机亚洲免费影院| 亚洲国产日韩一区二区| 中国国产av一级| 国产精品久久久av美女十八| 亚洲国产av新网站| 三上悠亚av全集在线观看| 久久久亚洲精品成人影院| 午夜日韩欧美国产| 国产一区二区激情短视频 | 一级a爱视频在线免费观看| 国产三级黄色录像| 欧美精品av麻豆av| av在线老鸭窝| 日韩一卡2卡3卡4卡2021年| 99热全是精品| 成年美女黄网站色视频大全免费| 真人做人爱边吃奶动态| 尾随美女入室| 国产视频首页在线观看| 日本欧美视频一区| 国产成人欧美在线观看 | 大片电影免费在线观看免费| 涩涩av久久男人的天堂| 女警被强在线播放| 欧美亚洲日本最大视频资源| 亚洲欧美一区二区三区黑人| 久久国产亚洲av麻豆专区| 欧美黑人精品巨大| 婷婷丁香在线五月| 精品亚洲乱码少妇综合久久| bbb黄色大片| 水蜜桃什么品种好| 天堂中文最新版在线下载| 国产1区2区3区精品| 欧美老熟妇乱子伦牲交| 欧美另类一区| 成人亚洲精品一区在线观看| 国产一区二区 视频在线| 国产精品香港三级国产av潘金莲 | 热99久久久久精品小说推荐| 大片电影免费在线观看免费| 国产xxxxx性猛交| bbb黄色大片| 一边摸一边做爽爽视频免费| 免费看十八禁软件| 精品福利观看| 免费观看人在逋| 真人做人爱边吃奶动态| 亚洲成av片中文字幕在线观看| 午夜激情久久久久久久| 日本五十路高清| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 久久99一区二区三区| 国产成人a∨麻豆精品| 国产一卡二卡三卡精品| 中文字幕最新亚洲高清| 考比视频在线观看| 美女扒开内裤让男人捅视频| 欧美在线黄色| 成人国产av品久久久| 丰满少妇做爰视频| videosex国产| 高潮久久久久久久久久久不卡| 午夜福利乱码中文字幕| 久久久国产一区二区| 国产亚洲精品第一综合不卡| 性少妇av在线| 中文字幕制服av| 日本av免费视频播放| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av高清一级| 少妇猛男粗大的猛烈进出视频| 高清视频免费观看一区二区| 亚洲精品国产色婷婷电影| 日韩视频在线欧美| 婷婷丁香在线五月| 日本欧美视频一区| 大片电影免费在线观看免费| 久久国产精品男人的天堂亚洲| 日本wwww免费看| 黄色 视频免费看| 男女边摸边吃奶| www日本在线高清视频| 90打野战视频偷拍视频| 免费一级毛片在线播放高清视频 | 免费高清在线观看视频在线观看| 一边摸一边抽搐一进一出视频| 色94色欧美一区二区| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 国产片内射在线| 免费黄频网站在线观看国产| 可以免费在线观看a视频的电影网站| 天天躁夜夜躁狠狠躁躁| 女警被强在线播放| av线在线观看网站| 精品国产一区二区久久| 精品人妻一区二区三区麻豆| 99久久精品国产亚洲精品| 999久久久国产精品视频| 日本av免费视频播放| 一级黄色大片毛片| 亚洲精品日韩在线中文字幕| 精品少妇一区二区三区视频日本电影| 美女大奶头黄色视频| 美女国产高潮福利片在线看| 天天操日日干夜夜撸| 亚洲精品日本国产第一区| 久久精品国产亚洲av高清一级| 天天躁日日躁夜夜躁夜夜| 夫妻性生交免费视频一级片| 精品国产国语对白av| 亚洲av电影在线进入| 丁香六月欧美| videosex国产| 美国免费a级毛片| 亚洲av电影在线进入| 精品高清国产在线一区| 最近中文字幕2019免费版| 国产免费视频播放在线视频| 视频区欧美日本亚洲| 啦啦啦在线免费观看视频4| 99国产精品免费福利视频| 搡老乐熟女国产| 国产一卡二卡三卡精品| 丰满迷人的少妇在线观看| 国产一卡二卡三卡精品| 欧美久久黑人一区二区| 日韩,欧美,国产一区二区三区| 人人澡人人妻人| 亚洲成av片中文字幕在线观看| 久久国产精品人妻蜜桃| 日韩,欧美,国产一区二区三区| 国产成人av激情在线播放| 中文字幕人妻丝袜一区二区| 50天的宝宝边吃奶边哭怎么回事| 日韩,欧美,国产一区二区三区| 精品一区二区三区av网在线观看 | 每晚都被弄得嗷嗷叫到高潮| 观看av在线不卡| 三上悠亚av全集在线观看| 国产极品粉嫩免费观看在线| 好男人电影高清在线观看| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 国产97色在线日韩免费| 日韩电影二区| 多毛熟女@视频| 免费看不卡的av| 男女高潮啪啪啪动态图| 欧美黑人欧美精品刺激| 国产免费视频播放在线视频| 精品久久久久久久毛片微露脸 | 久久国产精品影院| 一区二区三区精品91| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 久久久精品区二区三区| 国产熟女午夜一区二区三区| www.自偷自拍.com| 日韩大片免费观看网站| 精品国产一区二区三区四区第35| 国产精品欧美亚洲77777| 欧美变态另类bdsm刘玥| 色视频在线一区二区三区| 十八禁网站网址无遮挡| 少妇猛男粗大的猛烈进出视频| 捣出白浆h1v1| 午夜91福利影院| 51午夜福利影视在线观看| 国产av国产精品国产| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看 | 久久久国产欧美日韩av| 首页视频小说图片口味搜索 | 亚洲欧美日韩高清在线视频 | 国产伦理片在线播放av一区| 人人妻人人添人人爽欧美一区卜| 捣出白浆h1v1| 亚洲国产欧美网| 国产有黄有色有爽视频| 国产片特级美女逼逼视频| 1024视频免费在线观看| 晚上一个人看的免费电影| 国产欧美日韩综合在线一区二区| 国产野战对白在线观看| 欧美日韩亚洲国产一区二区在线观看 | 天天躁夜夜躁狠狠久久av| 九色亚洲精品在线播放| 成人手机av| 天天躁夜夜躁狠狠躁躁| 亚洲免费av在线视频| 精品视频人人做人人爽| 中国国产av一级| 看免费成人av毛片| 又黄又粗又硬又大视频| 国产女主播在线喷水免费视频网站| 在线观看免费午夜福利视频| 日韩熟女老妇一区二区性免费视频| 国产av国产精品国产| 黑丝袜美女国产一区| 色婷婷久久久亚洲欧美| 视频在线观看一区二区三区| 免费在线观看影片大全网站 | 黄色 视频免费看| 国产极品粉嫩免费观看在线| 免费黄频网站在线观看国产| 精品一区在线观看国产| 国产在线免费精品| 国语对白做爰xxxⅹ性视频网站| 精品免费久久久久久久清纯 | 少妇人妻久久综合中文| 欧美av亚洲av综合av国产av| av一本久久久久| 国产伦理片在线播放av一区| 色播在线永久视频| 高清黄色对白视频在线免费看| 亚洲精品一卡2卡三卡4卡5卡 | 各种免费的搞黄视频| 欧美在线黄色| 亚洲国产精品一区三区| 日韩熟女老妇一区二区性免费视频| 一本大道久久a久久精品| 亚洲av欧美aⅴ国产| 日韩中文字幕欧美一区二区 | 亚洲第一av免费看| 国产人伦9x9x在线观看| 九草在线视频观看| 久久精品亚洲熟妇少妇任你| 丝袜人妻中文字幕| 嫩草影视91久久| 久久久久视频综合| 赤兔流量卡办理| 一区二区三区四区激情视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人午夜精品| 人人妻人人澡人人看| 欧美在线黄色| 国产精品一区二区免费欧美 | 久久这里只有精品19| 午夜影院在线不卡| 欧美大码av| 大香蕉久久网| 久久精品亚洲熟妇少妇任你| 国产色视频综合| 亚洲精品国产av蜜桃| 国产精品亚洲av一区麻豆| 一区福利在线观看| 青青草视频在线视频观看| 国产高清不卡午夜福利| 午夜激情av网站| 99香蕉大伊视频| 女人高潮潮喷娇喘18禁视频| 国产片内射在线| 19禁男女啪啪无遮挡网站| 亚洲精品国产av成人精品| 麻豆av在线久日| 亚洲国产成人一精品久久久| 中文字幕色久视频| 国产一卡二卡三卡精品| 热re99久久国产66热| 男的添女的下面高潮视频| 性色av一级| 亚洲伊人久久精品综合| 一区二区日韩欧美中文字幕| 久久国产亚洲av麻豆专区| 欧美xxⅹ黑人| 蜜桃国产av成人99| 国产精品亚洲av一区麻豆| 宅男免费午夜| 91精品伊人久久大香线蕉| 人成视频在线观看免费观看| 蜜桃国产av成人99| 国产男人的电影天堂91| 中文字幕人妻丝袜一区二区| 人妻人人澡人人爽人人| 国产精品偷伦视频观看了| 亚洲av成人不卡在线观看播放网 | 人妻人人澡人人爽人人| 久久人妻熟女aⅴ| 亚洲自偷自拍图片 自拍| 黄色视频不卡| 色精品久久人妻99蜜桃| 欧美精品啪啪一区二区三区 | 久久久久久久大尺度免费视频| 欧美成狂野欧美在线观看| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 亚洲国产成人一精品久久久| 亚洲五月婷婷丁香| 男女国产视频网站| 亚洲精品久久成人aⅴ小说| 婷婷丁香在线五月| 青春草视频在线免费观看| 免费av中文字幕在线| 日韩欧美一区视频在线观看| 91麻豆av在线| 国产亚洲av片在线观看秒播厂| 观看av在线不卡| 国产亚洲av高清不卡| 十八禁人妻一区二区| 五月天丁香电影| 天天添夜夜摸| 亚洲色图综合在线观看| 国产免费现黄频在线看| 夜夜骑夜夜射夜夜干| 美女大奶头黄色视频| 国产日韩欧美视频二区| 亚洲精品日韩在线中文字幕| 久久精品久久久久久久性| 男人操女人黄网站| av天堂在线播放| 一级片'在线观看视频| 国产91精品成人一区二区三区 | 国产女主播在线喷水免费视频网站| 老司机在亚洲福利影院| a级毛片在线看网站| 亚洲精品在线美女| 国产成人91sexporn| 亚洲成人免费av在线播放| 欧美日韩亚洲国产一区二区在线观看 | 欧美亚洲日本最大视频资源| 一区在线观看完整版| 一区二区三区激情视频| 精品福利永久在线观看| 婷婷色av中文字幕| 免费不卡黄色视频| 精品视频人人做人人爽| 国产高清国产精品国产三级| 欧美变态另类bdsm刘玥| 色视频在线一区二区三区| 黄色一级大片看看| 女人爽到高潮嗷嗷叫在线视频| 日本a在线网址| 久久人人爽人人片av| 人人妻人人澡人人看| 亚洲国产毛片av蜜桃av| 黑人欧美特级aaaaaa片| 国产免费现黄频在线看| 超碰成人久久| 1024视频免费在线观看| 欧美日韩亚洲高清精品| 欧美黄色淫秽网站| 亚洲精品一二三| 在线观看免费高清a一片| 国产有黄有色有爽视频| 国产老妇伦熟女老妇高清| 亚洲精品自拍成人| 国产成人欧美| 欧美黄色淫秽网站| kizo精华| 国产1区2区3区精品| 黄频高清免费视频| 国产亚洲一区二区精品| 只有这里有精品99| av线在线观看网站| 多毛熟女@视频| 91精品三级在线观看| 啦啦啦视频在线资源免费观看| 亚洲人成77777在线视频| 亚洲自偷自拍图片 自拍| 亚洲精品日韩在线中文字幕| 亚洲av综合色区一区| 久久久精品国产亚洲av高清涩受| av天堂在线播放| 2018国产大陆天天弄谢| 人人妻人人爽人人添夜夜欢视频| 一区二区三区四区激情视频| 免费观看人在逋| 国产精品99久久99久久久不卡| 日本午夜av视频| 少妇精品久久久久久久| 免费av中文字幕在线| 在线精品无人区一区二区三| 久久精品久久精品一区二区三区| 黑丝袜美女国产一区| 不卡av一区二区三区|