• <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-),男,講師,碩士。

    精品福利观看| 十八禁网站免费在线| 叶爱在线成人免费视频播放| 黄色片一级片一级黄色片| 国产精品永久免费网站| 侵犯人妻中文字幕一二三四区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲七黄色美女视频| 日本撒尿小便嘘嘘汇集6| 最新美女视频免费是黄的| 亚洲色图av天堂| 99国产极品粉嫩在线观看| 很黄的视频免费| √禁漫天堂资源中文www| 91九色精品人成在线观看| 日韩免费av在线播放| 亚洲,欧美精品.| 电影成人av| 婷婷丁香在线五月| 国产伦一二天堂av在线观看| 神马国产精品三级电影在线观看 | 国产精品秋霞免费鲁丝片| 伦理电影免费视频| 欧美精品亚洲一区二区| 丰满的人妻完整版| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 欧美日韩精品网址| 久久天堂一区二区三区四区| 青草久久国产| 免费一级毛片在线播放高清视频 | 亚洲男人天堂网一区| 大型av网站在线播放| 三级毛片av免费| 人妻久久中文字幕网| 日本撒尿小便嘘嘘汇集6| 人妻久久中文字幕网| 久久人妻av系列| 夜夜夜夜夜久久久久| 亚洲欧美激情综合另类| 成年女人毛片免费观看观看9| 精品久久蜜臀av无| 久久国产精品人妻蜜桃| 91在线观看av| 亚洲成人免费电影在线观看| 亚洲一码二码三码区别大吗| 交换朋友夫妻互换小说| 丝袜人妻中文字幕| 中文字幕人妻熟女乱码| aaaaa片日本免费| 久久这里只有精品19| 一级片'在线观看视频| av天堂久久9| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区不卡视频| 水蜜桃什么品种好| 女人被躁到高潮嗷嗷叫费观| 国产xxxxx性猛交| 99精品欧美一区二区三区四区| 久久久精品欧美日韩精品| 啪啪无遮挡十八禁网站| 男人操女人黄网站| 一级片'在线观看视频| 午夜免费鲁丝| 岛国在线观看网站| 国产蜜桃级精品一区二区三区| 99在线人妻在线中文字幕| 1024视频免费在线观看| 一级毛片女人18水好多| 亚洲美女黄片视频| 久久人妻熟女aⅴ| 久久久久国产一级毛片高清牌| 免费在线观看视频国产中文字幕亚洲| av视频免费观看在线观看| 亚洲五月天丁香| 日韩欧美在线二视频| 在线观看66精品国产| 亚洲精品一二三| 亚洲五月婷婷丁香| 成人精品一区二区免费| 国产三级在线视频| 免费观看人在逋| 免费不卡黄色视频| 日本wwww免费看| 巨乳人妻的诱惑在线观看| 亚洲少妇的诱惑av| 亚洲精品在线美女| 国产免费av片在线观看野外av| 男女之事视频高清在线观看| 黄色片一级片一级黄色片| 成人18禁在线播放| 欧美激情久久久久久爽电影 | 中文字幕精品免费在线观看视频| 美女国产高潮福利片在线看| 999精品在线视频| 亚洲成人久久性| 亚洲自拍偷在线| 18禁国产床啪视频网站| 午夜福利,免费看| 在线观看www视频免费| 岛国在线观看网站| 日本欧美视频一区| 俄罗斯特黄特色一大片| 久久九九热精品免费| 国产成人啪精品午夜网站| 香蕉久久夜色| 三上悠亚av全集在线观看| 亚洲第一青青草原| 欧美成人性av电影在线观看| 国产成人av教育| 久久久久久亚洲精品国产蜜桃av| 天堂影院成人在线观看| 中出人妻视频一区二区| xxxhd国产人妻xxx| 欧美 亚洲 国产 日韩一| 亚洲国产精品合色在线| 嫩草影院精品99| 亚洲国产看品久久| 精品国产一区二区三区四区第35| 久久亚洲精品不卡| 久久久国产精品麻豆| 亚洲av日韩精品久久久久久密| 国产免费av片在线观看野外av| 国产一卡二卡三卡精品| 欧美黑人精品巨大| 日本黄色日本黄色录像| 久久精品国产亚洲av高清一级| 亚洲狠狠婷婷综合久久图片| 母亲3免费完整高清在线观看| 97超级碰碰碰精品色视频在线观看| 国产高清国产精品国产三级| 好看av亚洲va欧美ⅴa在| 午夜视频精品福利| 精品国产美女av久久久久小说| 国产精品自产拍在线观看55亚洲| 女性生殖器流出的白浆| 色播在线永久视频| ponron亚洲| 99国产精品免费福利视频| 久久午夜亚洲精品久久| 黄片播放在线免费| 久久香蕉精品热| 少妇被粗大的猛进出69影院| 亚洲一码二码三码区别大吗| 成熟少妇高潮喷水视频| 一本综合久久免费| 国产野战对白在线观看| 欧美日韩乱码在线| 欧美中文日本在线观看视频| 欧美乱色亚洲激情| 99久久人妻综合| 女性生殖器流出的白浆| 一区在线观看完整版| svipshipincom国产片| 欧美性长视频在线观看| 成人国产一区最新在线观看| 黄色女人牲交| 极品教师在线免费播放| 免费在线观看黄色视频的| 国产亚洲欧美精品永久| 黄色毛片三级朝国网站| 成人特级黄色片久久久久久久| avwww免费| 欧美日韩亚洲综合一区二区三区_| 99精品久久久久人妻精品| 国产aⅴ精品一区二区三区波| 97人妻天天添夜夜摸| 久久国产乱子伦精品免费另类| 国产一区二区在线av高清观看| 免费看十八禁软件| 日本三级黄在线观看| 在线十欧美十亚洲十日本专区| 欧美日本亚洲视频在线播放| 最近最新中文字幕大全免费视频| 性少妇av在线| 一级毛片精品| 精品久久久久久久久久免费视频 | 久久狼人影院| 日韩成人在线观看一区二区三区| 在线永久观看黄色视频| 99精国产麻豆久久婷婷| 中文字幕人妻熟女乱码| 久久国产精品人妻蜜桃| 久久国产精品影院| 在线观看一区二区三区| 老汉色av国产亚洲站长工具| a级毛片在线看网站| 熟女少妇亚洲综合色aaa.| 国产成人精品久久二区二区免费| 国产无遮挡羞羞视频在线观看| 两个人免费观看高清视频| а√天堂www在线а√下载| 久久这里只有精品19| 久久 成人 亚洲| 欧美亚洲日本最大视频资源| 午夜福利一区二区在线看| 性色av乱码一区二区三区2| 亚洲成国产人片在线观看| 国产精品乱码一区二三区的特点 | 成人特级黄色片久久久久久久| 国产单亲对白刺激| 久久香蕉激情| svipshipincom国产片| 丝袜美足系列| 丝袜在线中文字幕| 国产精品二区激情视频| 欧美乱妇无乱码| 女人被狂操c到高潮| 一区二区三区国产精品乱码| 丝袜在线中文字幕| 99热只有精品国产| 久久精品亚洲av国产电影网| 老司机午夜福利在线观看视频| 99久久99久久久精品蜜桃| 一进一出抽搐gif免费好疼 | 色婷婷久久久亚洲欧美| 日韩三级视频一区二区三区| av天堂久久9| 国产亚洲精品第一综合不卡| 一区二区三区精品91| 亚洲专区字幕在线| 国产一区二区激情短视频| √禁漫天堂资源中文www| 97碰自拍视频| 男人操女人黄网站| 国产视频一区二区在线看| 极品教师在线免费播放| 日本三级黄在线观看| 国产免费男女视频| 麻豆久久精品国产亚洲av | 好男人电影高清在线观看| 国产成年人精品一区二区 | 热99国产精品久久久久久7| 三上悠亚av全集在线观看| 人人妻人人爽人人添夜夜欢视频| 久久亚洲真实| www.精华液| 亚洲国产精品合色在线| 欧美人与性动交α欧美精品济南到| 精品欧美一区二区三区在线| av电影中文网址| 91字幕亚洲| 黄色视频不卡| tocl精华| 成人精品一区二区免费| 久久欧美精品欧美久久欧美| 久久精品人人爽人人爽视色| 99久久精品国产亚洲精品| 国产精品电影一区二区三区| 在线观看一区二区三区| 国产精品亚洲av一区麻豆| 妹子高潮喷水视频| 18美女黄网站色大片免费观看| 亚洲午夜理论影院| 亚洲中文av在线| 国产片内射在线| 大型黄色视频在线免费观看| 久久九九热精品免费| 久久久久国产精品人妻aⅴ院| 国产欧美日韩精品亚洲av| 一级片'在线观看视频| 97碰自拍视频| 别揉我奶头~嗯~啊~动态视频| 九色亚洲精品在线播放| 日韩欧美免费精品| 亚洲欧美一区二区三区久久| av有码第一页| 欧美丝袜亚洲另类 | 国产激情欧美一区二区| 国产深夜福利视频在线观看| 亚洲精品美女久久av网站| av片东京热男人的天堂| 水蜜桃什么品种好| 欧美日本中文国产一区发布| 桃色一区二区三区在线观看| 欧美人与性动交α欧美软件| 久久人人97超碰香蕉20202| 亚洲人成77777在线视频| bbb黄色大片| 91老司机精品| 欧美成人午夜精品| 啦啦啦 在线观看视频| 国产区一区二久久| 美女国产高潮福利片在线看| 黄频高清免费视频| 久久精品亚洲精品国产色婷小说| 久久久精品欧美日韩精品| 午夜免费鲁丝| 亚洲成av片中文字幕在线观看| 日韩欧美在线二视频| 亚洲自拍偷在线| 国产欧美日韩综合在线一区二区| 一进一出抽搐gif免费好疼 | 国产精品影院久久| 一进一出抽搐动态| 悠悠久久av| 国产精品久久视频播放| 欧美日韩精品网址| 精品免费久久久久久久清纯| 亚洲专区字幕在线| 少妇裸体淫交视频免费看高清 | 十八禁网站免费在线| 在线播放国产精品三级| 午夜福利一区二区在线看| 中文字幕色久视频| 美女大奶头视频| 久久这里只有精品19| 成人av一区二区三区在线看| 大香蕉久久成人网| 91在线观看av| 免费少妇av软件| a在线观看视频网站| 国产精品 欧美亚洲| 国产xxxxx性猛交| 午夜精品在线福利| 正在播放国产对白刺激| av免费在线观看网站| 久久久水蜜桃国产精品网| 91麻豆av在线| 欧美日韩黄片免| 国产欧美日韩精品亚洲av| 性欧美人与动物交配| 久久香蕉精品热| 色综合站精品国产| 亚洲欧美精品综合久久99| 99久久久亚洲精品蜜臀av| 欧美在线一区亚洲| 国产激情欧美一区二区| 亚洲伊人色综图| 母亲3免费完整高清在线观看| 麻豆一二三区av精品| 国产xxxxx性猛交| 在线观看日韩欧美| 免费日韩欧美在线观看| 久久精品影院6| 欧美日韩亚洲国产一区二区在线观看| 国产高清视频在线播放一区| 精品福利永久在线观看| 黑人欧美特级aaaaaa片| 久久婷婷成人综合色麻豆| 日韩欧美三级三区| 亚洲免费av在线视频| 国产无遮挡羞羞视频在线观看| 久久中文字幕人妻熟女| 人人妻人人添人人爽欧美一区卜| 国产高清激情床上av| 亚洲欧美激情综合另类| 亚洲一区二区三区色噜噜 | 后天国语完整版免费观看| 超碰成人久久| 日本黄色日本黄色录像| 99久久综合精品五月天人人| 亚洲精品一卡2卡三卡4卡5卡| 成在线人永久免费视频| 欧美久久黑人一区二区| 日本五十路高清| 国产精品久久久久久人妻精品电影| 可以在线观看毛片的网站| 在线视频色国产色| 真人一进一出gif抽搐免费| 精品卡一卡二卡四卡免费| 亚洲精品久久成人aⅴ小说| 日韩有码中文字幕| 高清欧美精品videossex| 免费在线观看完整版高清| 精品国产一区二区久久| 久热爱精品视频在线9| 女生性感内裤真人,穿戴方法视频| 成人三级做爰电影| www.自偷自拍.com| 国产97色在线日韩免费| 啦啦啦免费观看视频1| 中出人妻视频一区二区| 亚洲片人在线观看| 在线播放国产精品三级| 欧美黄色片欧美黄色片| av国产精品久久久久影院| 日韩三级视频一区二区三区| 午夜日韩欧美国产| 国产精品av久久久久免费| 欧美乱码精品一区二区三区| а√天堂www在线а√下载| 欧美精品啪啪一区二区三区| 最近最新免费中文字幕在线| 亚洲 欧美 日韩 在线 免费| 国产熟女午夜一区二区三区| 99久久人妻综合| 两人在一起打扑克的视频| 精品日产1卡2卡| 精品久久久久久,| 两人在一起打扑克的视频| 91成人精品电影| 成人国语在线视频| 男女下面进入的视频免费午夜 | 最近最新中文字幕大全免费视频| 国产精品一区二区精品视频观看| 日韩人妻精品一区2区三区| 又紧又爽又黄一区二区| 男人的好看免费观看在线视频 | 久久精品国产99精品国产亚洲性色 | √禁漫天堂资源中文www| 日本 av在线| 国产精品 欧美亚洲| 欧美午夜高清在线| 国产成人av激情在线播放| 久久久水蜜桃国产精品网| 在线观看免费高清a一片| 午夜福利免费观看在线| 欧美性长视频在线观看| 视频在线观看一区二区三区| 不卡一级毛片| 久久这里只有精品19| 亚洲av成人不卡在线观看播放网| 欧美精品一区二区免费开放| 99国产精品99久久久久| 久久影院123| www.自偷自拍.com| 中文欧美无线码| 老司机深夜福利视频在线观看| 久久久久久久午夜电影 | 国产高清视频在线播放一区| 99精国产麻豆久久婷婷| 黄色怎么调成土黄色| 成人国产一区最新在线观看| 亚洲人成电影观看| 欧美国产精品va在线观看不卡| 美女高潮到喷水免费观看| 久久伊人香网站| 欧美日韩亚洲高清精品| 国产av精品麻豆| 99精品在免费线老司机午夜| 一进一出好大好爽视频| 丰满迷人的少妇在线观看| 在线天堂中文资源库| 日日摸夜夜添夜夜添小说| 男女床上黄色一级片免费看| 脱女人内裤的视频| 欧美日韩亚洲综合一区二区三区_| 老司机亚洲免费影院| 中亚洲国语对白在线视频| 亚洲人成网站在线播放欧美日韩| 久99久视频精品免费| 亚洲全国av大片| 国产精品亚洲av一区麻豆| 一区福利在线观看| 黑人欧美特级aaaaaa片| 女生性感内裤真人,穿戴方法视频| 久久精品国产综合久久久| 亚洲国产精品一区二区三区在线| 久久国产精品影院| 久99久视频精品免费| 乱人伦中国视频| 国产av精品麻豆| 国产99久久九九免费精品| 久久久国产精品麻豆| 免费在线观看视频国产中文字幕亚洲| 免费看十八禁软件| 91字幕亚洲| 精品久久蜜臀av无| a级片在线免费高清观看视频| 欧美最黄视频在线播放免费 | 国产aⅴ精品一区二区三区波| 成年人免费黄色播放视频| 欧美日韩国产mv在线观看视频| 日本a在线网址| 欧美乱色亚洲激情| 国产一区二区三区在线臀色熟女 | 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲精品一区二区精品久久久| 久久久国产欧美日韩av| 日韩有码中文字幕| 悠悠久久av| 欧美日韩一级在线毛片| 999久久久精品免费观看国产| av网站免费在线观看视频| 久久精品国产综合久久久| 色综合站精品国产| 国产aⅴ精品一区二区三区波| 国产一区在线观看成人免费| 久久久久久亚洲精品国产蜜桃av| 久久人人97超碰香蕉20202| 午夜免费激情av| 一级a爱视频在线免费观看| 亚洲九九香蕉| 国产在线精品亚洲第一网站| 80岁老熟妇乱子伦牲交| 大码成人一级视频| 99国产极品粉嫩在线观看| 午夜精品久久久久久毛片777| 超碰97精品在线观看| 欧美精品一区二区免费开放| 97人妻天天添夜夜摸| 亚洲国产精品一区二区三区在线| 亚洲一区二区三区不卡视频| 亚洲av电影在线进入| 热re99久久国产66热| 黑丝袜美女国产一区| 国产精品一区二区精品视频观看| 亚洲人成电影免费在线| 亚洲欧美日韩无卡精品| xxx96com| 亚洲性夜色夜夜综合| 777久久人妻少妇嫩草av网站| 久久香蕉激情| 99re在线观看精品视频| 免费观看精品视频网站| 免费女性裸体啪啪无遮挡网站| 黄色成人免费大全| 免费在线观看黄色视频的| 99re在线观看精品视频| 看免费av毛片| 亚洲一码二码三码区别大吗| 99久久久亚洲精品蜜臀av| 精品久久蜜臀av无| 十八禁网站免费在线| 中文字幕精品免费在线观看视频| a在线观看视频网站| 久久久久亚洲av毛片大全| 国产激情久久老熟女| 日韩三级视频一区二区三区| 亚洲一区中文字幕在线| 日韩有码中文字幕| 久久精品aⅴ一区二区三区四区| 日日夜夜操网爽| 精品免费久久久久久久清纯| 亚洲精品国产区一区二| 性欧美人与动物交配| 亚洲av美国av| 老司机深夜福利视频在线观看| 91国产中文字幕| 精品久久久久久,| 久久人人爽av亚洲精品天堂| 久久人妻福利社区极品人妻图片| 亚洲 欧美 日韩 在线 免费| 十分钟在线观看高清视频www| a级片在线免费高清观看视频| 欧美乱码精品一区二区三区| 一区二区三区精品91| 91成人精品电影| 国产aⅴ精品一区二区三区波| 午夜福利,免费看| 国产99久久九九免费精品| 国产精品av久久久久免费| 亚洲熟女毛片儿| 国产国语露脸激情在线看| 亚洲 国产 在线| av片东京热男人的天堂| 午夜福利,免费看| 免费在线观看视频国产中文字幕亚洲| 后天国语完整版免费观看| 免费少妇av软件| av国产精品久久久久影院| 亚洲黑人精品在线| 91字幕亚洲| 久久久久精品国产欧美久久久| 嫩草影院精品99| 亚洲中文字幕日韩| 老司机午夜福利在线观看视频| 日本黄色日本黄色录像| 国产乱人伦免费视频| 久久这里只有精品19| 18禁国产床啪视频网站| 免费在线观看视频国产中文字幕亚洲| 99在线视频只有这里精品首页| 99re在线观看精品视频| 真人做人爱边吃奶动态| www.www免费av| 波多野结衣av一区二区av| 精品国产乱子伦一区二区三区| 国产成人精品久久二区二区91| 国产亚洲精品久久久久久毛片| 欧美激情极品国产一区二区三区| 九色亚洲精品在线播放| 亚洲成a人片在线一区二区| 日韩欧美三级三区| www日本在线高清视频| 一进一出抽搐动态| 91麻豆精品激情在线观看国产 | 久久天躁狠狠躁夜夜2o2o| 久久欧美精品欧美久久欧美| 丰满迷人的少妇在线观看| 日本a在线网址| 999久久久国产精品视频| 亚洲精品国产色婷婷电影| 在线天堂中文资源库| 欧美日韩av久久| 操出白浆在线播放| 天天躁狠狠躁夜夜躁狠狠躁| √禁漫天堂资源中文www| 日韩 欧美 亚洲 中文字幕| 美女高潮到喷水免费观看| 麻豆一二三区av精品| 叶爱在线成人免费视频播放| 99国产精品一区二区三区| 又大又爽又粗| 国产精品影院久久| 狂野欧美激情性xxxx| 男人舔女人下体高潮全视频| 久久午夜综合久久蜜桃| 88av欧美| 看黄色毛片网站| 一夜夜www| 午夜两性在线视频| 99久久国产精品久久久| 亚洲第一av免费看| 精品第一国产精品| 国产亚洲欧美在线一区二区| 免费日韩欧美在线观看| 高清毛片免费观看视频网站 | 免费av毛片视频| 欧美久久黑人一区二区| 国产黄a三级三级三级人| 亚洲av片天天在线观看| 叶爱在线成人免费视频播放| 欧美激情久久久久久爽电影 |