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

    一種基于徑向加速度的Singer-EKF機(jī)動(dòng)目標(biāo)跟蹤算法

    2014-07-12 12:00:23賈舒宜張赟

    賈舒宜,張赟

    (海軍航空工程學(xué)院a.電子信息工程系;b.飛行器工程系,山東煙臺(tái)264001)

    劉衛(wèi)華

    (海軍航空工程學(xué)院科研部,山東煙臺(tái)264001)

    一種基于徑向加速度的Singer-EKF機(jī)動(dòng)目標(biāo)跟蹤算法

    賈舒宜a,張赟b

    (海軍航空工程學(xué)院a.電子信息工程系;b.飛行器工程系,山東煙臺(tái)264001)

    針對(duì)雷達(dá)均不能提供目標(biāo)加速度信息,在目標(biāo)機(jī)動(dòng)時(shí)會(huì)出現(xiàn)跟蹤精度差甚至跟蹤發(fā)散的問(wèn)題,提出一種基于徑向加速度的Singer-EKF算法。該算法在信號(hào)處理階段利用Radon-Ambiguity變換(RAT)估計(jì)出目標(biāo)的徑向加速度,并通過(guò)坐標(biāo)轉(zhuǎn)換將其引入量測(cè)向量中,然后采用基于Singer模型的擴(kuò)展卡爾曼濾波(EKF)算法實(shí)現(xiàn)機(jī)動(dòng)目標(biāo)的跟蹤。仿真驗(yàn)證了該方法的有效性,并與傳統(tǒng)的不帶徑向加速度的擴(kuò)展卡爾曼濾波(EKF)方法進(jìn)行了比較,結(jié)果表明該方法在徑向距離、位置、加速度和速度估計(jì)精度方面都有所提高。

    機(jī)動(dòng)目標(biāo);徑向加速度;擴(kuò)展卡爾曼濾波;Singer模型;Radon-Ambiguity變換

    機(jī)動(dòng)目標(biāo)跟蹤一直是雷達(dá)、聲納和紅外應(yīng)用領(lǐng)域的難點(diǎn)。目前許多機(jī)動(dòng)目標(biāo)跟蹤方法,主要是利用機(jī)動(dòng)目標(biāo)模型和濾波算法實(shí)現(xiàn)機(jī)動(dòng)目標(biāo)跟蹤[1-8]。但現(xiàn)階段雷達(dá)均不能提供目標(biāo)加速度信息,因而在目標(biāo)機(jī)動(dòng)時(shí)會(huì)出現(xiàn)跟蹤精度差甚至跟蹤發(fā)散的問(wèn)題。如果雷達(dá)在提供目標(biāo)位置信息同時(shí),也能夠提供加速度信息,就可以提高機(jī)動(dòng)目標(biāo)的跟蹤性能。

    目前,關(guān)于信號(hào)處理階段提取徑向加速度的方法已有大量研究。當(dāng)雷達(dá)發(fā)射恒定載頻信號(hào)時(shí),機(jī)動(dòng)目標(biāo)回波信號(hào)可近似看做線性調(diào)頻(LFM)信號(hào),徑向加速度信息體現(xiàn)在雷達(dá)回波信號(hào)的調(diào)頻率上,因而對(duì)雷達(dá)目標(biāo)徑向加速度估計(jì)就是對(duì)回波信號(hào)調(diào)頻率的提取。目前,基于時(shí)頻分析方法的LFM信號(hào)檢測(cè)與估計(jì)技術(shù)不斷出現(xiàn)。例如Wigner-Ville分布(WVD)[9]是一種常用的LFM信號(hào)檢測(cè)估計(jì)方法,但在多目標(biāo)情況下容易受交叉項(xiàng)的影響使得估計(jì)性能降低。為了解決這一問(wèn)題,許多雙線性變換的時(shí)頻分析方法提出來(lái),例如自適應(yīng)chirplet變換(ACT)[10]、Radon-Ambiguity變換(RAT)[11]、分?jǐn)?shù)階傅里葉變換(FRFT)[12]等。其中,RAT的線性積分過(guò)程可以有效抑制交叉項(xiàng)干擾,并且計(jì)算量較小。因此,本文采用RAT方法在信號(hào)處理階段估計(jì)目標(biāo)徑向加速度,再將徑向加速度估計(jì)值進(jìn)行坐標(biāo)轉(zhuǎn)換引入濾波算法的量測(cè)向量中,有效地提高了機(jī)動(dòng)目標(biāo)的跟蹤精度。

    本文采用Singer模型作為機(jī)動(dòng)模型,原因在于Singer模型是一種典型的全局統(tǒng)計(jì)模型,兼顧了機(jī)動(dòng)模型和非機(jī)動(dòng)模型的精度。在目標(biāo)跟蹤系統(tǒng)中,本文采用了擴(kuò)展卡爾曼濾波算法(EKF)實(shí)現(xiàn)了徑向加速度量測(cè)到狀態(tài)向量的轉(zhuǎn)換。仿真驗(yàn)證了該方法的有效性,與傳統(tǒng)的不帶徑向加速度量測(cè)的Singer-EKF方法相比,該方法的跟蹤精度更高,并且在徑向距離、方位角、加速度和速度估計(jì)精度方面也有所改善。

    1 基于徑向加速度的Singer-EKF算法

    1.1 機(jī)動(dòng)目標(biāo)回波模型分析

    假設(shè)脈沖多普勒雷達(dá)發(fā)射恒定載頻信號(hào),對(duì)于遠(yuǎn)距離單機(jī)動(dòng)點(diǎn)目標(biāo),R(t)=R0+v0t+at2/2為其徑向運(yùn)動(dòng)規(guī)律。這只考慮加速運(yùn)動(dòng)情況,忽略短時(shí)間內(nèi)回波幅度的變化和距離徙動(dòng)影響,其多普勒變化率為式(1)中:R0為目標(biāo)到雷達(dá)的初始距離;λ為雷達(dá)發(fā)射信號(hào)波長(zhǎng);為初始速度,fd為多普勒頻率;為目標(biāo)徑向加速度,k為目標(biāo)加速度引起的多普勒變化率,可以看出其多普勒呈線性變化。

    不考慮傳播過(guò)程的衰減、低頻調(diào)制及混進(jìn)的噪聲,回波可以看成一個(gè)線性調(diào)頻信號(hào)。為了得到只與目標(biāo)有關(guān)的頻率信息,將信號(hào)與發(fā)射信號(hào)混頻,并省略掉相位中的常量,回波信號(hào)可表示為

    式(2)中:0≤t≤T,T為脈沖持續(xù)時(shí)間長(zhǎng)度;A為個(gè)信號(hào)分量的幅度,這里假設(shè)為1。

    由此可見(jiàn),回波信號(hào)的調(diào)頻率k包含了雷達(dá)的徑向加速度a的信息,a的計(jì)算公式為

    由式(3)可以看出,通過(guò)估計(jì)信號(hào)的調(diào)頻率k就可得到目標(biāo)的徑向加速度。

    1.2基于RAT的徑向加速度估計(jì)

    將持續(xù)時(shí)間有限的LFM信號(hào)(2)代入模糊函數(shù)定義中可得:

    對(duì)于任意時(shí)刻τ∈(-T,T),AFs(τ,f)在f=kτ處取最大值,其模糊函數(shù)呈背鰭狀,模糊函數(shù)及在時(shí)延—頻移平面上的投影分別如圖1和圖2所示。由此可見(jiàn),對(duì)于有限長(zhǎng)的線性調(diào)頻信號(hào)在時(shí)延—頻移平面上投影為過(guò)原點(diǎn)的一條直線f=kτ。由于有限長(zhǎng)的線性調(diào)頻信號(hào)的模糊函數(shù)為過(guò)原點(diǎn)的一條直線,而Radon變換具有沿直線積分累加的作用,因此我們將兩者結(jié)合,采用|AFs(τ,f)|在模糊平面(τ,f)上沿通過(guò)模糊原點(diǎn)的θ角方向的直線進(jìn)行積分,得到的結(jié)果即為Radon-Ambiguity變換[11],即:

    當(dāng)沿直線f=kτ(圖2)作直線積分時(shí),只需考慮搜索角度θ。因而對(duì)信號(hào)的模糊圖進(jìn)行Radon變換時(shí),可以把RAT(θ)看成關(guān)于搜索角度θ的線性累加器,只須進(jìn)行關(guān)于θ的一維搜索:

    經(jīng)過(guò)一系列搜索角度θ的計(jì)算和積累,最后在二維平面上出現(xiàn)一個(gè)峰值,這個(gè)峰值對(duì)應(yīng)的角度θ,就是估計(jì)調(diào)頻率所需的角度θ。徑向加速度估計(jì)具體步驟如下:

    1)回波信號(hào)序列為

    N為序列長(zhǎng)度,A是信號(hào)的幅度,采樣周期為Ts,w(n)為均值為零方差為σ2的高斯白噪聲。

    2)計(jì)算s(n)的模糊函數(shù)AF(τ,f),并求出模|AF(τ,f)。

    3)對(duì)||AF(τ,f)進(jìn)行Radon變換,并且進(jìn)行關(guān)于θ的一維搜索,并求得峰值對(duì)應(yīng)的θ角。

    圖1模糊函數(shù)圖Fig.1 Ambiguity function

    圖2模糊函數(shù)投影圖Fig.2 Projection of ambiguity function

    1.3 狀態(tài)方程

    狀態(tài)方程是由目標(biāo)的運(yùn)動(dòng)模型決定的,基于Singer模型的狀態(tài)方程可以表示為

    式(7)中,狀態(tài)向量

    xk、x˙k、x¨k和yk、y˙k、y¨k分別表示k時(shí)刻目標(biāo)x方向和y方向的位置、速度、加速度。

    狀態(tài)轉(zhuǎn)移矩陣

    式(9)中:α為機(jī)動(dòng)時(shí)間常數(shù);T為采樣周期。其離散時(shí)間過(guò)程噪聲w具有協(xié)方差:

    式(10)中,Q精確表達(dá)式為:

    1.4 量測(cè)方程

    設(shè)雷達(dá)在極坐標(biāo)系下位置量測(cè)(距離和方位角)、徑向速度量測(cè)和徑向加速度量測(cè)分別為zρ、zφ和za,則

    式(11)中:ρ、φ和a分別表示目標(biāo)距離、方位角和徑向加速度的真值;和分別表示方差為σ2ρ、σ2φ和σ2a的零均值高斯白噪聲。

    非線性量測(cè)方程為

    式中:Z(k)表示量測(cè)向量;H表示非線性量測(cè)轉(zhuǎn)移矩陣;W(k)表示方差為R(k)的高斯白噪聲。

    量測(cè)矩陣

    式中:

    量測(cè)噪聲協(xié)方差

    2 仿真驗(yàn)證

    在仿真過(guò)程中,假設(shè)雷達(dá)波長(zhǎng)為λ=8mm,雷達(dá)在信號(hào)處理階段采樣間隔為Ts=5×10-5s,觀測(cè)時(shí)間為T=50 ms,信噪比為SNR=0 dB,RAT算法中,搜索步長(zhǎng)為Δθ=0.125°,搜索半徑為Δρ=1.1Ts。雷達(dá)在數(shù)據(jù)處理階段采樣間隔為T′s=1 s,距離量測(cè)誤差為ρr=100 m,角度量測(cè)誤差為ρφ=0.5°,加速度測(cè)量誤差為ρa(bǔ)=1 m/s2,目標(biāo)的起始狀態(tài)為目標(biāo)運(yùn)動(dòng)過(guò)程歷時(shí)90 s,目標(biāo)發(fā)生機(jī)動(dòng)時(shí)刻及加速度見(jiàn)表1,目標(biāo)直角坐標(biāo)系下運(yùn)動(dòng)軌跡如圖3所示。

    表1目標(biāo)機(jī)動(dòng)運(yùn)動(dòng)情況表Tab.1 Movement of maneuvering target

    圖3目標(biāo)運(yùn)動(dòng)軌跡Fig.3 Trajectory of the target

    由上述條件下得到信號(hào)模型為

    式中,w(t)為均值為1的高斯白噪聲。在信號(hào)處理階段利用RAT方法得到目標(biāo)徑向加速度估計(jì),估計(jì)值結(jié)果如圖3和圖4所示,并將其作為量測(cè)向量中徑向加速度量測(cè)值。在數(shù)據(jù)處理階段,利用Singer模型和EKF算法對(duì)目標(biāo)進(jìn)行跟蹤,結(jié)果如圖5~12所示。從仿真結(jié)果中可以看出,目標(biāo)跟蹤精度明顯提高,原因在于徑向加速度的引入增加了量測(cè)向量的維數(shù),因而提高了跟蹤性能,同時(shí)目標(biāo)位置、加速度和速度的估計(jì)精度也有明顯改善。由圖5~12還可以看出,在x軸方向的跟蹤精度要明顯優(yōu)于y方向的跟蹤精度。這是由于當(dāng)目標(biāo)處于不同位置時(shí)引入徑向加速度對(duì)跟蹤精度影響程度是不一樣的,目標(biāo)跟蹤精度的改善效果與其所處位置有關(guān),目標(biāo)越靠近的哪個(gè)坐標(biāo)軸方向,哪個(gè)軸的跟蹤精度越高。

    圖4徑向加速度真實(shí)和估計(jì)值Fig.4 True values and estimations of radial acceleration

    圖5徑向距離誤差曲線Fig.5 Error curve of radial range

    圖6方位角誤差曲線Fig.6 Error curve of bearing

    圖7x軸位置誤差曲線Fig.7 Error curve of position inxaxis

    圖8y軸位置誤差曲線Fig.8 Error curve of position inyaxis

    圖9x軸速度誤差曲線Fig.9 Error curve of velocity inxaxis

    圖10y軸速度誤差曲線Fig.10 Error curve of velocity inyaxis

    圖11x軸加速度誤差曲線Fig.11 Error curve of acceleration inxaxis

    圖12y軸加速度誤差曲線Fig.12 Error curve of acceleration inyaxis

    3 結(jié)論

    本文提出一種帶徑向加速度的Singer-EKF機(jī)動(dòng)目標(biāo)跟蹤算法。由于徑向加速度是在信號(hào)處理階段利用RAT方法能精確地提取出來(lái)的,因而將其引入測(cè)量向量可提高機(jī)動(dòng)目標(biāo)跟蹤精度。仿真結(jié)果驗(yàn)證了該方法在機(jī)動(dòng)條件下的跟蹤效果要優(yōu)于傳統(tǒng)的Singer-EKF算法。利用該方法的研究成果,可以在雷達(dá)提供目標(biāo)位置和速度信息的同時(shí),提供目標(biāo)的徑向加速度信息,從而實(shí)現(xiàn)對(duì)機(jī)動(dòng)目標(biāo)的精確跟蹤。

    [1]TIMOTHY L MOLLOY,JASON J FORD.HMM triangle relative entropy concepts in sequential change detection applied to vision-based dim target maneuver detection[C]//15th International Conference on Information Fusion.Singapore:Interational Society of Information Fusion,2012:1120-1125.

    [2]GAO L,XING J P,MA Z L.Improved IMM algorithm for nonlinear maneuvering target tracking[J].Proedia Engineering,2012,29(6):4117-4123.

    [3]FANTACCI C,BATTISTELLI G,CHISCI L,et al.Multiple-model algorithms for distributed tracking of a maneuvering target[C]//15th International Conference on Information Fusion.Singapore:Interational Society of Information Fusion,2012:1231-1235.

    [4]JIAO LIANMENG,PAN QUAN,LIANG YAN,et al.A nonlinear tracking algorithm with range-rate measurements based on unbiased measurement conversion[C]// 15th International Conference on Information Fusion.Singapore:Interational Society of Information Fusion,2012:1400-1405.

    [5]PENG D L,GU Y.IMM algorithm for a 3D high maneuvering target tracking[J].Lecture Notes in Computer Science,2011,6729(1):529-536.

    [6]GAO CAICAI,CHEN WEI.Ground moving target tracking with vs-imm using mean shift unscented particle filter [J].Chinese Journal of Aeronautics,2011,24(5):622-630.

    [7]XIA X G.Discrete chirp-fourier transform and its application to chirp rate estimation[J].IEEE Transactions on signal processing,2000,47(11):3122-3133.

    [8]YANG X J,XING K Y,F(xiàn)ENG X L.Maneuvering target tracking in dense clutter based on particle filtering[J]. Chinese Journal ofAeronautics,2011,24(2):171-180.

    [9]KUMAR K ANOOP,ARVIND MANOJ,DIVAKAR K,et al.A novel time-frequency approach for acceleration estimation from a single PR[C]//Fifth International Symposium Signal Processing and it’s applications.Brisbane:Signal Processing Research Centre,1999:531-534.

    [10]XIA X G.Discrete chirp-fourier transform and its application to chirp rate estimation[J].IEEE Transactions on signal processing,2000,47(11):3122-3133.

    [11]孫軍海,李亞超,刑蒙道,等.基于Radon-Ambiguity變換的多目標(biāo)檢測(cè)[J].西安交通大學(xué)學(xué)報(bào),2009,43(4):81-89. SUN JUNHAI,LI YACHAO,XING MENGDAO,et al. Detection of multiple targets based on radon-ambiguity transformation[J].Journal of Xi’an Jiaotong University,2009,43(4):81-89.(in Chinese)

    [12]杜文超,王國(guó)宏,高學(xué)強(qiáng).低信噪比條件下在短時(shí)信號(hào)中提取目標(biāo)徑向加速度方法研究[J].中國(guó)科學(xué):技術(shù)科學(xué),2007,37(7):923-943. DU WENCHAO,WANG GUOHONG,GAO XUEQIANG.Acceleration estimation from short time duration in low SNR[J].Science China:Technology Scinces,2007,37(7):923-943.(in Chinese)

    JIA Shu-yia,ZHANG Yunb

    (Naval Aeronautical and Astronautical university a.Department of Electronic and Information Engineering; b.Department of Airborne Vehicle Engineering,Yantai Shandong 264001,China)

    海雜波中廣義符號(hào)恒虛警檢測(cè)算法性能分析

    劉衛(wèi)華

    (海軍航空工程學(xué)院科研部,山東煙臺(tái)264001)

    摘 要:分析了廣義符號(hào)檢測(cè)算法在仿真的高斯雜波背景和實(shí)測(cè)海雜波背景下,對(duì)2種目標(biāo)(Swerling 0型和Swerling II型)的檢測(cè)性能,以及對(duì)實(shí)際漁船目標(biāo)的檢測(cè)性能。研究表明,隨著脈沖數(shù)、參考單元數(shù)和信雜比的提高,該檢測(cè)算法的檢測(cè)性能有所提高;在低信雜比條件下,GS檢測(cè)算法對(duì)Swerling II型目標(biāo)的檢測(cè)性能優(yōu)于對(duì)Swerling 0型目標(biāo)的檢測(cè)性能,在高信雜比的條件下,對(duì)Swerling 0型目標(biāo)的檢測(cè)性能優(yōu)于對(duì)Swerling II型目標(biāo)的檢測(cè)性能。

    關(guān)鍵詞:非參量;廣義符號(hào);海雜波;信雜比

    在復(fù)雜的噪聲背景下保持恒定的虛警概率是設(shè)計(jì)目標(biāo)檢測(cè)算法必須考慮的問(wèn)題[1]。當(dāng)雜波包絡(luò)分布類型已知時(shí),通過(guò)采集數(shù)據(jù)來(lái)估計(jì)某些參數(shù),從而對(duì)信號(hào)進(jìn)行有效的檢測(cè),這種檢測(cè)方法被稱為參量檢測(cè)方法。然而,如果實(shí)際雜波分布和假設(shè)差距較大,檢測(cè)性能就會(huì)大幅下降。因此,設(shè)計(jì)虛警概率與背景分布無(wú)關(guān)的檢測(cè)算法很有必要,這種檢測(cè)算法被稱為非參量檢測(cè)算法。在背景分布已知的情況下,非參量檢測(cè)算法的性能往往不如參量檢測(cè)算法,但是在背景分布未知的情況下,非參量檢測(cè)算法的性能往往好于參量檢測(cè)算法[2]。在雷達(dá)實(shí)際工作環(huán)境當(dāng)中,雜波分布類型往往無(wú)法準(zhǔn)確預(yù)知,因而非參量檢測(cè)方法受到越來(lái)越廣泛的關(guān)注[3-5]。常見(jiàn)的非參量檢測(cè)算法有符號(hào)檢測(cè)算法[6],Wilcoxon檢測(cè)算法[7],廣義符號(hào)(GS)檢測(cè)算法[2]和Mann-Whitney檢測(cè)算法[8]。GS檢測(cè)算法處理的是經(jīng)包絡(luò)檢波器輸出的回波信號(hào),將檢測(cè)單元與周圍參考單元進(jìn)行比較,并得到相應(yīng)的非參量化檢驗(yàn)統(tǒng)計(jì)量[7],是易實(shí)現(xiàn)的兩樣本非參量檢測(cè)算法之一。

    為了研究GS檢測(cè)算法對(duì)不同類型目標(biāo)的檢測(cè)性能,本文在仿真雜波數(shù)據(jù)、實(shí)測(cè)海雜波環(huán)境中,對(duì)GSCFAR檢測(cè)算法的檢測(cè)性能和海雜波抑制能力進(jìn)行了仿真分析。由于GS檢測(cè)算法的性能與脈沖數(shù)和參考單元數(shù)有關(guān)[2],因而本文還研究了脈沖數(shù)和參考單元數(shù)對(duì)GS檢測(cè)算法性能的影響。

    1 算法模型

    1.1 GS檢測(cè)算法原理

    GS檢測(cè)算法是對(duì)接收信號(hào)進(jìn)行多脈沖檢測(cè)的兩樣本符號(hào)檢測(cè)算法。該算法將檢測(cè)單元的幅度值與周圍參考單元的幅度值進(jìn)行比較,得到檢測(cè)單元的檢驗(yàn)統(tǒng)計(jì)量,再將檢驗(yàn)統(tǒng)計(jì)量與對(duì)應(yīng)的檢測(cè)閾值比較。圖1給出了GS檢測(cè)算法的原理框圖[7],檢測(cè)單元采樣用yj(j=1,2,…,N),參考單元采樣xji(i=1,2,…,M)。檢測(cè)單元兩側(cè)各有一個(gè)保護(hù)單元,防止檢測(cè)單元中可能的目標(biāo)信號(hào)泄漏到鄰近的參考單元中而影響檢測(cè)。

    圖1 GS檢測(cè)算法原理框圖Fig.1 Schematic diagram of GS detection algorithm

    檢驗(yàn)統(tǒng)計(jì)量

    式中:U(·)是單位階躍函數(shù);N為脈沖數(shù);M為參考單元數(shù)。

    將檢驗(yàn)統(tǒng)計(jì)量TGS和檢測(cè)閾值K進(jìn)行比較,如果TGS≥K,則判決存在目標(biāo);如果TGS<K,則判決不存在目標(biāo)。在無(wú)目標(biāo)條件下,假設(shè)yj滿足獨(dú)立同分布;在有目標(biāo)條件下,假設(shè)yj和xji滿足獨(dú)立同分布[2]。

    1.2 檢測(cè)閾值K的確定

    虛警概率Pfa可根據(jù)脈沖數(shù)N、參考單元數(shù)M和檢測(cè)閾值K的大小來(lái)設(shè)定[2],因而Pfa的值可以N、M和K進(jìn)行函數(shù)表示,具體關(guān)系如式(2)所示:

    式中,

    從式(2)可以看出,Pfa不依賴于雜波和噪聲的分布形式,只與K、N和M有關(guān),因而GS檢測(cè)算法可使虛警概率保持恒定[7]。

    當(dāng)Pfa、N和M已知時(shí),K可以由Pfa反解出來(lái),那么K也不依賴于雜波和噪聲的分布形式。當(dāng)Pfa=10-4、N=10、M=16時(shí),根據(jù)上述公式得出K=127。

    2 GS檢測(cè)算法性能分析

    2.1 仿真雜波+目標(biāo)回波條件下的性能分析

    本節(jié)采用高斯分布來(lái)模擬產(chǎn)生海雜波。分別產(chǎn)生I、Q2路標(biāo)準(zhǔn)高斯分布(均值為0,方差為1)序列,求包絡(luò)后作為背景雜波數(shù)據(jù)。

    產(chǎn)生不起伏Swerling 0型目標(biāo)回波和快起伏Swerling II型目標(biāo)回波,將其與雜波數(shù)據(jù)合成,并調(diào)整信雜比(SCR)。

    信雜比的計(jì)算公式為

    當(dāng)SCR已知時(shí),可以根據(jù)上式反推出目標(biāo)信號(hào)的幅度,從而計(jì)算出在不同信雜比情況下的檢測(cè)概率[9]。

    圖2給出了Pfa=10-3,脈沖數(shù)和參考單元數(shù)分別改變時(shí),2種目標(biāo)模型條件下,SCR與檢測(cè)概率關(guān)系曲線。

    圖2基于仿真數(shù)據(jù)的檢測(cè)概率與信雜比關(guān)系曲線Fig.2 Relationship curve of detection probability and SNR based on simulated data

    由圖2可以看出:

    1)當(dāng)脈沖數(shù)和參考單元數(shù)一定,信雜比在一定范圍內(nèi)變化時(shí),檢測(cè)概率隨著信雜比的增大而不斷增大;

    2)在信雜比和參考單元數(shù)一定的條件下,隨著脈沖數(shù)的增大,檢測(cè)概率有所增大;

    3)在信雜比和脈沖數(shù)一定的條件下,隨著參考單元數(shù)的增加,檢測(cè)概率有所增大;

    4)脈沖數(shù)變化對(duì)檢測(cè)Swerling 0型目標(biāo)和Swerling II型目標(biāo)的影響要強(qiáng)于參考單元數(shù)對(duì)檢測(cè)2種目標(biāo)的影響;

    5)在低信雜比條件下,GS檢測(cè)算法對(duì)Swerling II型目標(biāo)的檢測(cè)性能優(yōu)于對(duì)Swerling 0型目標(biāo)的檢測(cè)性能;在高信雜比條件下,對(duì)Swerling 0型目標(biāo)的檢測(cè)性能優(yōu)于對(duì)Swerling II型目標(biāo)的檢測(cè)性能。

    2.2 真實(shí)雜波+仿真目標(biāo)回波條件下的性能分析

    為進(jìn)一步驗(yàn)證不同參數(shù)對(duì)GS檢測(cè)算法性能的影響情況,本節(jié)基于某實(shí)測(cè)數(shù)據(jù)對(duì)GS檢測(cè)算法性能進(jìn)行分析,實(shí)際海雜波背景較為復(fù)雜,且分布類型一般是未知的。設(shè)Pfa=10-3,雜波數(shù)據(jù)是通過(guò)雷達(dá)采集得到的,目標(biāo)類型仍采用Swerling 0型和Swerling II型。圖3給出N和M分別變化時(shí),2種目標(biāo)模型的SCR與檢測(cè)概率的對(duì)應(yīng)關(guān)系曲線。

    圖3基于實(shí)測(cè)數(shù)據(jù)的檢測(cè)概率與信雜比關(guān)系曲線Fig.3 Relationship curve of detection probability and SNR based on real data

    由圖3可以看出:

    1)在其他條件相同的情況下,檢測(cè)概率隨脈沖數(shù)的增大而有所增大。對(duì)于實(shí)測(cè)數(shù)據(jù),在高信雜比條件下,對(duì)Swerling 0型目標(biāo)的檢測(cè)能力優(yōu)于對(duì)Swerling II型目標(biāo)的檢測(cè)能力;在低信雜比條件下,對(duì)Swerling II型目標(biāo)的檢測(cè)能力優(yōu)于對(duì)Swerling 0型目標(biāo)的檢測(cè)能力;

    2)在其他條件相同的情況下,檢測(cè)概率隨著參考單元數(shù)的增大而有所增大,但變化不明顯;

    3)脈沖數(shù)變化對(duì)檢測(cè)Swerling 0型目標(biāo)和Swerling II型目標(biāo)的影響要強(qiáng)于參考單元數(shù)對(duì)檢測(cè)Swerling 0型目標(biāo)和Swerling II型目標(biāo)的影響。

    上述分析進(jìn)一步說(shuō)明,GS檢測(cè)算法是一種與雜波分布類型無(wú)關(guān)的檢測(cè)方法,即在不同的雜波分布類型情況下,幾種參數(shù)對(duì)非參量檢測(cè)算法性能的影響趨勢(shì)保持一致。

    分析原因:

    1)Swerling II型目標(biāo)為快起伏目標(biāo),脈沖和脈沖間的起伏是統(tǒng)計(jì)獨(dú)立的,若前一個(gè)脈沖沒(méi)有超過(guò)門限,則相繼脈沖有可能超過(guò)門限;而不起伏目標(biāo)的相關(guān)性較強(qiáng),如果第一個(gè)脈沖振幅小于檢測(cè)門限,則相繼脈沖也不會(huì)超過(guò)門限。在低信雜比條件下,信號(hào)幅度較低,快起伏目標(biāo)經(jīng)過(guò)檢測(cè)算法處理之后得到的檢驗(yàn)統(tǒng)計(jì)量超過(guò)檢測(cè)閾值的概率相對(duì)較大而更容易被檢測(cè)到;在高信雜比條件下,信號(hào)幅度相對(duì)較大,相對(duì)穩(wěn)定的不起伏目標(biāo)經(jīng)過(guò)檢測(cè)算法處理之后得到的檢驗(yàn)統(tǒng)計(jì)量要高于快起伏目標(biāo),從而更容易被檢測(cè)到;

    2)GS檢測(cè)算法將檢測(cè)單元和周圍的參考單元進(jìn)行比較以得到檢驗(yàn)統(tǒng)計(jì)量,其本質(zhì)是對(duì)目標(biāo)優(yōu)勢(shì)的積累。若增大脈沖數(shù),則在脈沖間相對(duì)穩(wěn)定的有用信號(hào)的優(yōu)勢(shì)可得到進(jìn)一步的積累,檢測(cè)效果會(huì)明顯改善。

    2.3 真實(shí)雜波+真實(shí)目標(biāo)回波條件下的性能分析

    本節(jié)基于雷達(dá)對(duì)海上目標(biāo)檢測(cè)的實(shí)測(cè)數(shù)據(jù)進(jìn)行分析,以驗(yàn)證GS檢測(cè)算法對(duì)漁船目標(biāo)的檢測(cè)性能。圖4給出了虛警概率Pfa=10-3、脈沖數(shù)N和參考單元數(shù)M分別變化時(shí)的二維平面圖。圖4中縱坐標(biāo)表示距離單元,橫坐標(biāo)表示重復(fù)掃描的次數(shù),每一個(gè)像素點(diǎn)表示一次掃描的回波強(qiáng)度[10],圖4 a)為原始數(shù)據(jù)平面圖,圖4 b)~d)為GS檢測(cè)算法處理后的顯示結(jié)果。

    圖4實(shí)測(cè)數(shù)據(jù)顯示圖和檢測(cè)算法處理結(jié)果Fig.4 Radar original video data and the result of detection algorithm

    由圖4 a)可知,在60 km處有1個(gè)漁船目標(biāo)湮沒(méi)在海雜波當(dāng)中,較難分辨;經(jīng)過(guò)GS檢測(cè)檢測(cè)算法處理后的結(jié)果如圖4 b)~d)所示,絕大部分雜波被抑制掉,畫(huà)面非常清晰。從圖4可以看出:

    1)相比于原始數(shù)據(jù)平面圖,GS檢測(cè)算法成功檢測(cè)出了目標(biāo),該檢測(cè)算法對(duì)雜波有很好的抑制作用。當(dāng)N增大時(shí),顏色不一的間斷條狀回波逐漸變?yōu)檫B續(xù)的條狀回波;

    2)對(duì)于1 000次觀測(cè),當(dāng)M為20,N為10時(shí),檢測(cè)到目標(biāo)259次;當(dāng)M為20,N為20時(shí),檢測(cè)到目標(biāo)422次;當(dāng)M為30,N為10時(shí),檢測(cè)到目標(biāo)272次。由此可知,隨著脈沖數(shù)N的增多檢測(cè)效果明顯變好,但參考單元數(shù)M對(duì)檢測(cè)效果的影響并不明顯。

    3 總結(jié)

    本文利用Swerling 0和Swerling II型目標(biāo)回波疊加仿真雜波和實(shí)測(cè)海雜波數(shù)據(jù),研究了不同參數(shù),即信雜比、參考單元數(shù)和脈沖數(shù)對(duì)GS檢測(cè)算法性能的影響,并采用真實(shí)海雜波+目標(biāo)回波數(shù)據(jù)進(jìn)行驗(yàn)證。結(jié)果表明:

    1)在高斯雜波背景和實(shí)測(cè)海雜波背景下,隨著脈沖數(shù)、參考單元數(shù)和信雜比的增大檢測(cè)概率不斷增大,但脈沖數(shù)對(duì)檢測(cè)概率的影響強(qiáng)于參考單元數(shù)對(duì)檢測(cè)概率的影響;在高信雜比條件下,對(duì)Swerling 0型目標(biāo)的檢測(cè)能力優(yōu)于對(duì)Swerling II型目標(biāo)的檢測(cè)能力;在低信雜比條件下,對(duì)Swerling II型目標(biāo)的檢測(cè)能力優(yōu)于對(duì)Swerling 0型目標(biāo)的檢測(cè)能力。

    2)通過(guò)對(duì)實(shí)測(cè)數(shù)據(jù)驗(yàn)證表明,GS檢測(cè)算法可以檢測(cè)到海雜波中的漁船,在虛警概率保持不變的條件下,隨著脈沖數(shù)和參考單元數(shù)的增大,算法的檢測(cè)性能得到改善,但是參考單元數(shù)對(duì)檢測(cè)性能的影響并不明顯。

    參考文獻(xiàn):

    [1]馬克·理查滋.雷達(dá)信號(hào)處理基礎(chǔ)[M].3版.北京:電子工業(yè)出版社,2008:260-286. MARK A RICHARDS.Fundamentals of radar signal processing[M].3rd ed.Beijing:Publishing House of Electronics Industry,2008:260-286.(in Chinese)

    [2]VILHELM GREGERS HENSEN.Nonparametic radar extraction using a generalized sign test[J].IEEE Transactions on Aerospace&Electronic Systems,1971,7(5):941-950.

    [3]SEYEF B,SHARAFAT A R.Signed-rank nonparametric multiuser detection in non-gaussian channels[J].IEEE Transactions on Information Theory,2005,51(4):1478-1486.

    [4]HAO CHEN.Noise enhanced nonpatametric detection[J]. IEEE Transactions on Information Theory,2009,55(2):499-506.

    [5]ORTON M.Nonparametric rank detectors under k-distributed clutter in radar application[J].IEEE Transactions on Aerospace&Electronic Systems,2005,41(2):702-710.

    [6]CARLYLE J W.On nonparametric siganl detectors[J]. IEEE Transactions on Information Theory,1964,10(2):146-152.

    [7]何友,關(guān)鍵.雷達(dá)目標(biāo)檢測(cè)與恒虛警處理[M].2版.北京:清華大學(xué)出版社,2011:230-250. HE YOU,GUAN JIAN.Radar target detection and cfar processing[M].2rd ed.Beijing:Tsinghua University Press,2011:230-250.(in Chinese)

    [8]ZEORI G W.Performance of a two-sample mann-whitney nonparametric detector in a radar application[J]. IEEE Transactions on Aerospace&Electronic Systems,1971,7(5):951-959.

    [9]趙志堅(jiān).海雜波中非參量恒虛警檢測(cè)算法性能分析[J].雷達(dá)科學(xué)與技術(shù),2010,8(1):65-73. ZHAO ZHIJIAN.Performance analysis of nonparametric CFAR detectors in sea clutter[J].Radar Science and Technology,2010,8(1):65-73.(in Chinese)

    [10]張林.三種兩樣本非參量檢測(cè)器CFAR檢測(cè)算法分析[J].海軍航空工程學(xué)院學(xué)報(bào),2011,26(3):296-300. ZHANG LIN.Performance analysis of three kinds of twosample nonparametric CFAR detectors[J].Journal of Navel Aeronautical and Astronautical University,2011,26(3):296-300.(in Chinese)

    Abstract:The performance of GS detector for two kinds of target(Swerling 0 and Swerling II)was studied under Gaussian clutter and real sea clutter.The analysis results showed that,with the numbers of pulse,reference cells and signal-toclutter ratio(SCR)increasing,the performance of GS detector rose.In the condition of low SCR,the performance of GS detector detecting Swerling II targets was better than that of GS detector detecting Swerling 0 targets;and in the condition of high SCR,the performance of GS detector detecting Swerling 0 targets was better than that of GS detector detecting Swerling II targets.

    Key words:nonparametric;generalized sign;sea clutter;SCR

    Tracking Algorithm of Singer-EKF Maneuvering Target Based on Radial Acceleration

    The performance of maneuvering tracking can be improved if acceleration is derived.Therefore,a Singer-EKF method based on Radon-Ambiguity transformation(RAT)was proposed for enhancing the tracking of a maneuvering target.In the proposed method,the radial acceleration was derived based on Radon-Ambiguity transformation in signal processing and then brought into the measurement vector coordinates transform model.In the filtering approach,a method of Extended Kalman Filter(EKF)based on Singer model was adopted to resolve the problem of the non-linearity of the measurement equation.The tracking performance of the proposed method is evaluated and compared with the traditional Extended Kalman Filter(EKF)algorithm based on the Singer model without acceleration measurement.It is showned that the proposed algorithm outperformed better in maneuvering scenario,and the estimate accuracies of radial distance,position, acceleration and velocity were improved at the same time.

    maneuvering target;radial acceleration;EKF;Singer model;Radon-Ambiguity transformation

    Performance Analysis of Generalized Sign CFAR Detection Algorithms Within Sea Clutterr

    LIU Wei-Hua
    (Department of Scientific Research,NAAU,Yantai Shandong 264001,China)

    TN911.23

    A

    1673-1522(2014)02-0141-05

    10.7682/j.issn.1673-1522.2014.02.009

    TN957文獻(xiàn)標(biāo)志碼:A

    ??編號(hào):1673-1522(2014)02-0146-05DOI:10.7682/j.issn.1673-1522.2014.02.010

    2013-10-09;

    2014-02-12

    國(guó)家自然科學(xué)基金資助項(xiàng)目(61102165;61102167)

    賈舒宜(1984-),女,博士。

    收稿日期:2013-12-16;修回日期:2014-02-10

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(61002045,61179017,61201445,61302008)

    作者簡(jiǎn)介:劉衛(wèi)華(1973-),男,講師,碩士。

    亚洲 欧美 日韩 在线 免费| 亚洲成av人片免费观看| 天堂动漫精品| 色综合亚洲欧美另类图片| 久久精品国产清高在天天线| 欧美午夜高清在线| 丰满的人妻完整版| 香蕉av资源在线| 久久99热这里只有精品18| 久久婷婷人人爽人人干人人爱| avwww免费| 国产精品永久免费网站| 精品熟女少妇八av免费久了| 久久午夜亚洲精品久久| 欧美色视频一区免费| 毛片女人毛片| 国产真实乱freesex| 午夜精品久久久久久毛片777| 麻豆成人午夜福利视频| 国产真实乱freesex| 亚洲人与动物交配视频| 亚洲人成伊人成综合网2020| 国产视频内射| 99久久精品国产亚洲精品| av片东京热男人的天堂| 精品无人区乱码1区二区| 久久精品国产99精品国产亚洲性色| 欧美另类亚洲清纯唯美| 免费看日本二区| 婷婷丁香在线五月| 美女扒开内裤让男人捅视频| 99热这里只有精品一区 | 少妇粗大呻吟视频| 一级黄色大片毛片| 国产精品99久久99久久久不卡| 国产精品1区2区在线观看.| 在线观看66精品国产| 精品国产超薄肉色丝袜足j| 999精品在线视频| 欧美又色又爽又黄视频| av中文乱码字幕在线| 精品人妻1区二区| 美女免费视频网站| 国产三级在线视频| 美女午夜性视频免费| 大型黄色视频在线免费观看| 国内少妇人妻偷人精品xxx网站 | 亚洲九九香蕉| 久久久久久国产a免费观看| 淫妇啪啪啪对白视频| 国产av一区在线观看免费| 欧美中文日本在线观看视频| 欧美丝袜亚洲另类 | 久久这里只有精品中国| 久久99热这里只有精品18| 亚洲成人国产一区在线观看| 免费av毛片视频| 日本 av在线| 男女做爰动态图高潮gif福利片| 亚洲成人精品中文字幕电影| 国产精品久久久久久亚洲av鲁大| 男人舔女人的私密视频| 国产三级在线视频| 精品人妻1区二区| 欧美在线黄色| 一边摸一边做爽爽视频免费| 女人爽到高潮嗷嗷叫在线视频| 成人精品一区二区免费| 男女那种视频在线观看| 欧美一级a爱片免费观看看 | 久久99热这里只有精品18| 亚洲18禁久久av| 国产精品日韩av在线免费观看| 久久欧美精品欧美久久欧美| 我的老师免费观看完整版| www.自偷自拍.com| 97人妻精品一区二区三区麻豆| 亚洲国产看品久久| 亚洲国产精品sss在线观看| 美女午夜性视频免费| 国产av一区二区精品久久| 亚洲国产欧洲综合997久久,| 视频区欧美日本亚洲| 欧美日韩中文字幕国产精品一区二区三区| 我要搜黄色片| 久久精品亚洲精品国产色婷小说| 成人国语在线视频| 免费在线观看成人毛片| 国产亚洲欧美98| 色av中文字幕| 脱女人内裤的视频| 麻豆一二三区av精品| 欧美乱码精品一区二区三区| 国产精品一区二区三区四区免费观看 | 亚洲一码二码三码区别大吗| 1024视频免费在线观看| 精品一区二区三区四区五区乱码| 特大巨黑吊av在线直播| 蜜桃久久精品国产亚洲av| 精品国产乱子伦一区二区三区| 亚洲最大成人中文| 又黄又爽又免费观看的视频| 国产亚洲av高清不卡| 嫩草影视91久久| 国产一区在线观看成人免费| 一二三四社区在线视频社区8| 一个人免费在线观看的高清视频| 黄色视频,在线免费观看| 欧美成狂野欧美在线观看| 美女午夜性视频免费| 亚洲熟女毛片儿| 亚洲av片天天在线观看| 国内精品久久久久精免费| 国产亚洲欧美98| 国产精品久久视频播放| 久久精品国产清高在天天线| 欧美一区二区精品小视频在线| 女人爽到高潮嗷嗷叫在线视频| 国产人伦9x9x在线观看| 可以在线观看的亚洲视频| 国产黄a三级三级三级人| 日韩 欧美 亚洲 中文字幕| 一进一出抽搐gif免费好疼| 又黄又粗又硬又大视频| 国产高清videossex| e午夜精品久久久久久久| 老鸭窝网址在线观看| 90打野战视频偷拍视频| 国产精品,欧美在线| 一二三四社区在线视频社区8| 91成年电影在线观看| 丝袜人妻中文字幕| 97超级碰碰碰精品色视频在线观看| 在线播放国产精品三级| 人人妻人人看人人澡| 国产黄片美女视频| 欧美一区二区精品小视频在线| 久久精品aⅴ一区二区三区四区| av有码第一页| av视频在线观看入口| 亚洲人与动物交配视频| 美女扒开内裤让男人捅视频| 亚洲午夜理论影院| 激情在线观看视频在线高清| 日本免费一区二区三区高清不卡| 久久精品夜夜夜夜夜久久蜜豆 | 十八禁人妻一区二区| 精品一区二区三区四区五区乱码| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕最新亚洲高清| 男人舔女人下体高潮全视频| 久久久水蜜桃国产精品网| 99久久无色码亚洲精品果冻| 狂野欧美白嫩少妇大欣赏| 成年免费大片在线观看| 日韩高清综合在线| 久久久久性生活片| 国产片内射在线| 亚洲中文日韩欧美视频| 在线观看免费视频日本深夜| 欧美日本亚洲视频在线播放| 深夜精品福利| 制服丝袜大香蕉在线| 欧美成人一区二区免费高清观看 | 别揉我奶头~嗯~啊~动态视频| 人妻久久中文字幕网| 欧美一级a爱片免费观看看 | 叶爱在线成人免费视频播放| 91国产中文字幕| 欧美极品一区二区三区四区| 99精品在免费线老司机午夜| 手机成人av网站| 巨乳人妻的诱惑在线观看| 不卡一级毛片| 中文字幕人成人乱码亚洲影| 亚洲国产欧洲综合997久久,| 亚洲精华国产精华精| 免费在线观看亚洲国产| 精品第一国产精品| 18禁观看日本| 亚洲中文字幕一区二区三区有码在线看 | 久久99热这里只有精品18| 日韩欧美 国产精品| 免费无遮挡裸体视频| 日韩国内少妇激情av| 精品高清国产在线一区| 国产精品一区二区免费欧美| 国产成人欧美在线观看| 99精品久久久久人妻精品| 亚洲成av人片免费观看| 久久久久九九精品影院| 一级片免费观看大全| 亚洲av熟女| 国产高清有码在线观看视频 | 久久九九热精品免费| 国产视频内射| 亚洲午夜理论影院| 国产爱豆传媒在线观看 | 成熟少妇高潮喷水视频| 精品高清国产在线一区| 日韩中文字幕欧美一区二区| 巨乳人妻的诱惑在线观看| 久久久久久大精品| 日韩有码中文字幕| 黄色丝袜av网址大全| 无限看片的www在线观看| 日韩欧美 国产精品| 国产精品,欧美在线| av超薄肉色丝袜交足视频| 97超级碰碰碰精品色视频在线观看| 亚洲自拍偷在线| 美女高潮喷水抽搐中文字幕| av视频在线观看入口| 国产亚洲欧美98| 一级毛片高清免费大全| 亚洲精品中文字幕在线视频| 欧美午夜高清在线| 日本 欧美在线| 亚洲男人的天堂狠狠| 亚洲成av人片在线播放无| 日韩免费av在线播放| 天堂√8在线中文| 日韩高清综合在线| www.精华液| 99国产精品一区二区三区| 国产野战对白在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日韩大码丰满熟妇| 亚洲国产精品久久男人天堂| 欧美 亚洲 国产 日韩一| 午夜视频精品福利| 日韩欧美三级三区| 母亲3免费完整高清在线观看| 亚洲av电影不卡..在线观看| 国产免费男女视频| 99国产极品粉嫩在线观看| 亚洲人成77777在线视频| 精品熟女少妇八av免费久了| 九色国产91popny在线| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品在线美女| 91麻豆精品激情在线观看国产| 99久久无色码亚洲精品果冻| 久久中文字幕人妻熟女| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 97碰自拍视频| 日韩 欧美 亚洲 中文字幕| 99国产综合亚洲精品| 午夜老司机福利片| 国产激情偷乱视频一区二区| 日韩欧美 国产精品| 久久香蕉国产精品| 国产真实乱freesex| 久久热在线av| 日韩中文字幕欧美一区二区| 久久精品aⅴ一区二区三区四区| 成人午夜高清在线视频| 亚洲一区高清亚洲精品| 熟妇人妻久久中文字幕3abv| 午夜免费成人在线视频| 国产成年人精品一区二区| 男女视频在线观看网站免费 | 999久久久国产精品视频| 日本精品一区二区三区蜜桃| 欧美成狂野欧美在线观看| 精品久久久久久久久久久久久| 岛国在线观看网站| 久久国产乱子伦精品免费另类| 久久中文字幕一级| 成人国语在线视频| 亚洲自偷自拍图片 自拍| 巨乳人妻的诱惑在线观看| 波多野结衣高清作品| 国模一区二区三区四区视频 | 久久伊人香网站| 777久久人妻少妇嫩草av网站| 亚洲片人在线观看| а√天堂www在线а√下载| 特大巨黑吊av在线直播| 亚洲av五月六月丁香网| 男人舔女人下体高潮全视频| 啪啪无遮挡十八禁网站| 欧美成人免费av一区二区三区| www日本在线高清视频| 亚洲片人在线观看| 久久精品成人免费网站| 国内久久婷婷六月综合欲色啪| 黄色成人免费大全| 黄片小视频在线播放| 国产欧美日韩一区二区精品| 怎么达到女性高潮| 日韩欧美 国产精品| or卡值多少钱| 久久精品国产99精品国产亚洲性色| 在线十欧美十亚洲十日本专区| 久久亚洲精品不卡| 国产精品久久视频播放| xxx96com| 日韩精品青青久久久久久| 99精品欧美一区二区三区四区| 日本一二三区视频观看| 欧美一级a爱片免费观看看 | 精品欧美国产一区二区三| 国产探花在线观看一区二区| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| 久久久久久亚洲精品国产蜜桃av| 在线十欧美十亚洲十日本专区| 婷婷六月久久综合丁香| 日韩高清综合在线| 黄色a级毛片大全视频| 亚洲一区二区三区色噜噜| 老汉色av国产亚洲站长工具| 在线看三级毛片| 岛国在线观看网站| 久久99热这里只有精品18| 午夜福利视频1000在线观看| 欧美久久黑人一区二区| 999精品在线视频| 成在线人永久免费视频| 老司机福利观看| 熟女电影av网| 亚洲精品美女久久av网站| 午夜老司机福利片| 99国产精品一区二区三区| 久久精品国产亚洲av香蕉五月| 丝袜人妻中文字幕| 一级毛片女人18水好多| 国产av在哪里看| 男女之事视频高清在线观看| 高清在线国产一区| 我的老师免费观看完整版| 黄色视频,在线免费观看| 久久精品国产99精品国产亚洲性色| 国产精品一区二区三区四区免费观看 | 国产99白浆流出| 久久久久久国产a免费观看| 欧美中文日本在线观看视频| 国产精品99久久99久久久不卡| 黄色女人牲交| 亚洲精品久久国产高清桃花| 国内揄拍国产精品人妻在线| 亚洲熟妇熟女久久| 成人av在线播放网站| 成人永久免费在线观看视频| 黄色女人牲交| 青草久久国产| 日日爽夜夜爽网站| 亚洲国产精品999在线| 99国产精品一区二区蜜桃av| 亚洲第一电影网av| 国产爱豆传媒在线观看 | 国产精品久久久久久亚洲av鲁大| 国产黄a三级三级三级人| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区免费观看 | 此物有八面人人有两片| 日日干狠狠操夜夜爽| 黄色a级毛片大全视频| 无遮挡黄片免费观看| 18禁裸乳无遮挡免费网站照片| 一级毛片高清免费大全| 国产精品香港三级国产av潘金莲| 成年人黄色毛片网站| 搡老熟女国产l中国老女人| 久久这里只有精品中国| 成年人黄色毛片网站| 亚洲专区字幕在线| 国产免费av片在线观看野外av| 18禁黄网站禁片午夜丰满| 岛国在线观看网站| 波多野结衣高清作品| 性欧美人与动物交配| 亚洲人成网站在线播放欧美日韩| 国产黄a三级三级三级人| 狂野欧美激情性xxxx| 麻豆成人av在线观看| 舔av片在线| 国产精品 欧美亚洲| 国产爱豆传媒在线观看 | 国产熟女xx| 婷婷六月久久综合丁香| 午夜激情福利司机影院| 精品福利观看| 90打野战视频偷拍视频| www.熟女人妻精品国产| 我的老师免费观看完整版| 亚洲九九香蕉| 国产成人av激情在线播放| 男女午夜视频在线观看| 日本一二三区视频观看| 欧美乱码精品一区二区三区| 国产精品久久久久久精品电影| 久久久国产成人免费| 亚洲av五月六月丁香网| 亚洲国产欧美人成| 老司机午夜福利在线观看视频| 亚洲国产高清在线一区二区三| 九色成人免费人妻av| 91麻豆av在线| 少妇粗大呻吟视频| 色播亚洲综合网| 天堂av国产一区二区熟女人妻 | 嫩草影院精品99| 久久伊人香网站| 神马国产精品三级电影在线观看 | 手机成人av网站| 欧美一区二区国产精品久久精品 | 黄色成人免费大全| 免费av毛片视频| 在线观看午夜福利视频| 免费一级毛片在线播放高清视频| 美女午夜性视频免费| 日韩有码中文字幕| 波多野结衣巨乳人妻| 久久久久久人人人人人| 男女视频在线观看网站免费 | 12—13女人毛片做爰片一| 欧美在线一区亚洲| 精品高清国产在线一区| 1024香蕉在线观看| 亚洲一区高清亚洲精品| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 一级毛片精品| 国内揄拍国产精品人妻在线| 非洲黑人性xxxx精品又粗又长| 搡老熟女国产l中国老女人| 99精品欧美一区二区三区四区| 欧美一级毛片孕妇| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久久毛片微露脸| 首页视频小说图片口味搜索| 亚洲欧美激情综合另类| 又紧又爽又黄一区二区| 成在线人永久免费视频| 欧美一级a爱片免费观看看 | 久久久久久亚洲精品国产蜜桃av| 国产一区二区三区视频了| 69av精品久久久久久| netflix在线观看网站| 国产精品九九99| 九色国产91popny在线| 国产野战对白在线观看| 欧美日韩一级在线毛片| 久久香蕉精品热| 欧美成人性av电影在线观看| 在线国产一区二区在线| 久久中文看片网| 久久伊人香网站| 少妇人妻一区二区三区视频| 三级毛片av免费| 精品国产超薄肉色丝袜足j| 人成视频在线观看免费观看| 在线观看舔阴道视频| АⅤ资源中文在线天堂| 日本三级黄在线观看| 亚洲美女黄片视频| 久久午夜亚洲精品久久| 国产精品98久久久久久宅男小说| 丝袜美腿诱惑在线| 国产高清videossex| 亚洲美女黄片视频| 黄色片一级片一级黄色片| 小说图片视频综合网站| 欧美性猛交╳xxx乱大交人| 成人欧美大片| 亚洲激情在线av| 精品日产1卡2卡| 天天躁狠狠躁夜夜躁狠狠躁| bbb黄色大片| 精品一区二区三区视频在线观看免费| 色播亚洲综合网| 亚洲最大成人中文| 一二三四社区在线视频社区8| 热99re8久久精品国产| 亚洲中文日韩欧美视频| 亚洲一区二区三区色噜噜| 国产真人三级小视频在线观看| 中出人妻视频一区二区| 日本在线视频免费播放| 亚洲成av人片在线播放无| 久久久久久人人人人人| 99久久精品热视频| 久久这里只有精品19| 欧美又色又爽又黄视频| 国产av又大| 又紧又爽又黄一区二区| 欧美在线一区亚洲| 香蕉久久夜色| videosex国产| 可以在线观看毛片的网站| 中文字幕人妻丝袜一区二区| 观看免费一级毛片| 18禁国产床啪视频网站| 亚洲第一欧美日韩一区二区三区| 老熟妇仑乱视频hdxx| 国产亚洲av嫩草精品影院| 老熟妇乱子伦视频在线观看| 18禁美女被吸乳视频| 韩国av一区二区三区四区| 亚洲国产欧美人成| 99精品在免费线老司机午夜| 88av欧美| 久久久久精品国产欧美久久久| 国产一区二区在线观看日韩 | 一本大道久久a久久精品| 国产一区二区三区在线臀色熟女| 国产99久久九九免费精品| 免费在线观看完整版高清| 国产av一区在线观看免费| av欧美777| 精品久久蜜臀av无| 国产野战对白在线观看| 午夜老司机福利片| 一区二区三区高清视频在线| 国产熟女午夜一区二区三区| 欧美色视频一区免费| 亚洲色图av天堂| 亚洲七黄色美女视频| 国产成+人综合+亚洲专区| 最近视频中文字幕2019在线8| 久久 成人 亚洲| 亚洲国产欧美人成| 欧美最黄视频在线播放免费| 国产成年人精品一区二区| 99热这里只有是精品50| 美女黄网站色视频| 黄色丝袜av网址大全| 日韩欧美一区二区三区在线观看| 国产久久久一区二区三区| 色噜噜av男人的天堂激情| 国产成人aa在线观看| 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看 | 在线观看免费视频日本深夜| 99久久无色码亚洲精品果冻| 亚洲欧美日韩东京热| 亚洲欧美一区二区三区黑人| 色综合站精品国产| 国产激情偷乱视频一区二区| 色综合站精品国产| 日韩欧美在线乱码| 天天添夜夜摸| 欧美日韩国产亚洲二区| 91av网站免费观看| 欧美性长视频在线观看| 国产不卡一卡二| 舔av片在线| 国产亚洲精品av在线| 久久人妻福利社区极品人妻图片| 91国产中文字幕| 日本成人三级电影网站| 免费搜索国产男女视频| 一级黄色大片毛片| 在线观看午夜福利视频| 国产亚洲精品久久久久久毛片| 这个男人来自地球电影免费观看| 国产精品免费视频内射| 又大又爽又粗| 亚洲精品在线美女| 成人av在线播放网站| 亚洲av成人不卡在线观看播放网| 久久这里只有精品19| 最好的美女福利视频网| 国产97色在线日韩免费| 啦啦啦免费观看视频1| 蜜桃久久精品国产亚洲av| 美女免费视频网站| 精品久久久久久成人av| 男女下面进入的视频免费午夜| 国产久久久一区二区三区| 日韩欧美一区二区三区在线观看| 观看免费一级毛片| 欧美午夜高清在线| 日韩大尺度精品在线看网址| 亚洲av日韩精品久久久久久密| 亚洲av片天天在线观看| 极品教师在线免费播放| 亚洲va日本ⅴa欧美va伊人久久| а√天堂www在线а√下载| 老司机在亚洲福利影院| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看 | 国产亚洲精品久久久久久毛片| 男插女下体视频免费在线播放| 亚洲一区二区三区色噜噜| 欧美人与性动交α欧美精品济南到| 少妇粗大呻吟视频| 亚洲人成伊人成综合网2020| 日韩欧美精品v在线| 成人三级做爰电影| 精品国产美女av久久久久小说| 九色国产91popny在线| 床上黄色一级片| 91九色精品人成在线观看| 中文字幕熟女人妻在线| 欧美绝顶高潮抽搐喷水| 亚洲精华国产精华精| 日本黄大片高清| 热99re8久久精品国产| 国产精品久久久av美女十八| 亚洲国产中文字幕在线视频| 色综合亚洲欧美另类图片| 九色国产91popny在线| 亚洲中文av在线| 免费在线观看视频国产中文字幕亚洲| av视频在线观看入口| 又紧又爽又黄一区二区| 久久久国产成人免费| 看免费av毛片| 午夜福利免费观看在线| 亚洲欧美日韩东京热| 日韩高清综合在线|