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

    基于離散頻率編碼的PD雷達(dá)二維解模糊方法

    2021-07-05 00:55:30孫希平盧景月鐘衛(wèi)軍
    關(guān)鍵詞:信噪比多普勒波形

    孫希平, 盧景月, 張 磊,*, 鐘衛(wèi)軍

    (1. 中山大學(xué)電子與通信工程學(xué)院, 廣東 廣州 510006;2. 西安電子科技大學(xué)電子工程學(xué)院, 陜西 西安 710071;3. 西北工業(yè)大學(xué)自動(dòng)化學(xué)院, 陜西 西安 710072)

    0 引 言

    雷達(dá)作為一種利用電磁波輻射的主動(dòng)探測(cè)傳感器,可以對(duì)超遠(yuǎn)距目標(biāo)進(jìn)行探測(cè),這一特性決定了雷達(dá)成為戰(zhàn)機(jī)上的關(guān)鍵設(shè)備,對(duì)于戰(zhàn)機(jī)在戰(zhàn)場(chǎng)上提前發(fā)現(xiàn)目標(biāo)的能力至關(guān)重要。機(jī)載脈沖多普勒(pulse Doppler, PD)雷達(dá)的中脈沖重復(fù)頻率(medium pulse repetition frequency, MPRF)波形具有全方位的探測(cè)能力[1],因此在機(jī)載火控中備受重視。但是MPRF對(duì)目標(biāo)的探測(cè)可能同時(shí)存在距離和多普勒的二維模糊,需要對(duì)二維模糊進(jìn)行正確處理以后,才能進(jìn)行準(zhǔn)確的目標(biāo)參數(shù)提取以供戰(zhàn)場(chǎng)態(tài)勢(shì)評(píng)估。

    就雷達(dá)工作原理來(lái)說(shuō),距離模糊的本質(zhì)在于脈沖信號(hào)的周期性發(fā)射,多普勒模糊在于多普勒采樣率的不足,因此兩者是一對(duì)相互矛盾的量。對(duì)于距離模糊,可以使用距離搜索方法[2]、混合濾波[3]解模糊,或者使用脈沖間相位編碼[4-5]、脈沖頻率編碼[6]等,對(duì)發(fā)射波形序列進(jìn)行標(biāo)記以解決距離模糊的問(wèn)題。但是這種方法只能解決距離模糊。由于需要多個(gè)脈沖作為一組信號(hào)進(jìn)行距離解模糊,等效脈沖重復(fù)頻率降低,因此增加了多普勒的模糊。文獻(xiàn)[7]中,提出一種使用正交頻分復(fù)用(orthogonal frequency division multiplexing, OFDM)正交波形進(jìn)行距離解模糊的方案,為了保證一定范圍的多普勒測(cè)量,需要通過(guò)增加脈沖重復(fù)頻率獲得距離不模糊區(qū)域的擴(kuò)展。多普勒模糊通??梢酝ㄟ^(guò)變脈沖重復(fù)頻率[8]、結(jié)合測(cè)距信息[9]、利用信號(hào)多普勒敏感性[10]等方法解決。通過(guò)改變脈沖重復(fù)頻率或者信號(hào)波長(zhǎng),可以改變目標(biāo)速度對(duì)應(yīng)的多普勒頻率,通過(guò)多組信息的匹配可以對(duì)速度進(jìn)行解模糊[8]。文獻(xiàn)[10]提出了一種利用信號(hào)多普勒敏感性解多普勒模糊的方法,利用不同模糊數(shù)搜索時(shí)目標(biāo)脈壓結(jié)果的幅度,根據(jù)最大匹配濾波峰值確定多普勒模糊數(shù),以此來(lái)解決多普勒模糊問(wèn)題。經(jīng)典的多脈沖重復(fù)頻率(multiple pulse repetition frequency, Multi-PRF)方法是一種有效的距離-多普勒二維解模糊方法[11-13]。Multi-PRF在被提出后,已經(jīng)廣泛應(yīng)用在MPRF模式的信號(hào)處理中,并且演變出了多個(gè)形式。在Multi-PRF模式下,利用中國(guó)余數(shù)定理[11-12]或者距離-多普勒擴(kuò)展聚類[13]的方法,可以對(duì)目標(biāo)距離和多普勒進(jìn)行同時(shí)解模糊。但是在各脈沖重復(fù)頻率(pulse repetition frequency, PRF)探測(cè)結(jié)果匹配過(guò)程中,容易產(chǎn)生虛假目標(biāo)。且為了準(zhǔn)確匹配,需要目標(biāo)在駐留時(shí)間內(nèi)不發(fā)生越距離單元運(yùn)動(dòng),對(duì)高速目標(biāo)進(jìn)行探測(cè)時(shí)需要限制每個(gè)PRF脈沖簇的脈沖個(gè)數(shù),造成積累信噪比增益有限,在低信噪比情況下,可能無(wú)法有效檢測(cè)目標(biāo),造成漏警。

    本文針對(duì)機(jī)載火控雷達(dá)MPRF下的二維(距離-多普勒)模糊問(wèn)題以及傳統(tǒng)方法無(wú)法長(zhǎng)時(shí)間相參積累等缺點(diǎn),提出了一種基于循環(huán)發(fā)射正交離散頻率編碼波形進(jìn)行二維解模糊的方法,通過(guò)發(fā)射正交波形,不同距離模糊區(qū)域的回波對(duì)應(yīng)不同的頻率編碼波形,使用由發(fā)射波形設(shè)計(jì)的匹配濾波器組可以進(jìn)行不同距離模糊區(qū)域的分離,并且高速目標(biāo)在較長(zhǎng)時(shí)間的回波中,不可避免地產(chǎn)生距離走動(dòng),為了進(jìn)行相參積累,利用Keystone進(jìn)行距離走動(dòng)校正。但是當(dāng)回波存在多普勒模糊時(shí),Keystone無(wú)法拉直回波包絡(luò),導(dǎo)致相參處理積累增益降低,因此可以搜索最佳相參積累增益下的補(bǔ)償多普勒模糊數(shù)[14-15],確定回波的多普勒模糊次數(shù),實(shí)現(xiàn)固定重頻下的多普勒解模糊。對(duì)距離走動(dòng)校正后,目標(biāo)處于相同的距離單元,可以進(jìn)行駐留時(shí)間內(nèi)所有回波的相參積累,獲得較高的相參處理增益,增加了在低信噪比回波情況下的檢測(cè)及參數(shù)估計(jì)性能。

    DFC波形與相位編碼波形類似,具有圖釘型的模糊函數(shù)[16-17],存在多普勒敏感性,在對(duì)信號(hào)的匹配濾波中,如果回波信號(hào)的多普勒頻移較大,則會(huì)嚴(yán)重影響信號(hào)的脈壓結(jié)果。本文針對(duì)回波信號(hào)存在的多普勒敏感性,將匹配濾波過(guò)程放在多普勒域內(nèi)進(jìn)行,對(duì)多普勒域內(nèi)的每條多普勒譜線設(shè)計(jì)針對(duì)性的匹配濾波函數(shù),不同多普勒的目標(biāo)回波使用對(duì)應(yīng)脈內(nèi)的多普勒匹配濾波器,實(shí)現(xiàn)精確地匹配濾波,避免了由多普勒敏感性引起的匹配濾波增益降低或虛假目標(biāo)產(chǎn)生,保證對(duì)不同速度目標(biāo)的檢測(cè)能力。

    1 機(jī)載PD雷達(dá)回波模型

    1.1 機(jī)載PD雷達(dá)的MPRF回波模型

    機(jī)載火控PD雷達(dá)的MPRF情況下,多個(gè)位于不同距離模糊區(qū)域目標(biāo)的回波模型如圖1所示。假設(shè)雷達(dá)發(fā)射脈沖重復(fù)頻率為T(mén)r,其最大不模糊距離為Rum=cTr/2,其中c為光速。當(dāng)目標(biāo)距離大于Rum時(shí),脈沖信號(hào)的重復(fù)周期小于雷達(dá)回波到達(dá)時(shí)間,產(chǎn)生距離模糊。如圖1所示,不同模糊區(qū)域的目標(biāo)回波疊加在一起到達(dá)雷達(dá)接收機(jī)。

    圖1 MPRF模式下距離模糊模型Fig.1 Range ambiguity model of MPRF

    假設(shè)雷達(dá)使用一組脈沖進(jìn)行循環(huán)發(fā)射,脈組中第n個(gè)發(fā)射脈沖可以表示為

    sn(t)=un(t)exp(j2πfct)

    (1)

    式中:t表示時(shí)間;fc為發(fā)射信號(hào)載頻;un(t)為第n個(gè)發(fā)射脈沖的脈內(nèi)調(diào)制項(xiàng)。

    設(shè)雷達(dá)視場(chǎng)中存在一個(gè)距離模糊的目標(biāo),考慮在空空模式下,迎飛目標(biāo)具有高相對(duì)速度,雷達(dá)發(fā)射第n個(gè)脈沖時(shí)接收到的目標(biāo)回波可以表示為

    (2)

    (3)

    式中:k為目標(biāo)的距離模糊次數(shù),k滿足

    Rumk

    (4)

    式中:R=R0+vtm表示目標(biāo)的距離,R0為觀測(cè)起始時(shí)間的目標(biāo)距離,v為目標(biāo)速度。

    式(2)中,τ(tm)為目標(biāo)在tm時(shí)刻的時(shí)延,考慮光速遠(yuǎn)大于目標(biāo)速度,該時(shí)延可以表示為

    (5)

    目標(biāo)回波信號(hào)被接收后,使用載頻fc下變頻至基頻信號(hào),基頻信號(hào)可以表示為

    (6)

    (7)

    脈內(nèi)多普勒會(huì)影響脈沖脈內(nèi)調(diào)制項(xiàng)的快時(shí)間頻域,在匹配濾波過(guò)程產(chǎn)生失配。

    如上所述,在機(jī)載火控雷達(dá)中,距離-多普勒二維模糊必須得到解決,且在空空迎飛情況下,還需要考慮脈內(nèi)多普勒調(diào)制問(wèn)題。

    1.2 離散頻率編碼信號(hào)

    雷達(dá)使用正交離散頻率編碼(discrete frequency coding, DFC)波形在空域探測(cè)目標(biāo),循環(huán)發(fā)射N個(gè)脈沖信號(hào),每個(gè)發(fā)射脈沖具有不同的頻率編碼,N個(gè)脈沖信號(hào)構(gòu)成一個(gè)脈組,在圖1中展示了使用5個(gè)不同頻率編碼脈沖組成的一個(gè)發(fā)射脈組。圖2(a)為一個(gè)典型DFC波形的示意圖,發(fā)射脈沖由一系列點(diǎn)頻窄脈沖信號(hào)隨機(jī)排列拼接組成,每一個(gè)點(diǎn)頻信號(hào)稱為一個(gè)碼片。式(1)中信號(hào)的脈內(nèi)調(diào)制項(xiàng)可以表示為

    圖2 離散頻率編碼波形Fig.2 Discrete frequency coding waveform

    (8)

    式中:rect(·)表示矩形函數(shù);M為一個(gè)脈沖內(nèi)的頻率碼片個(gè)數(shù);Tsp為每個(gè)頻率碼片的時(shí)間寬度,整個(gè)脈沖時(shí)間寬度為MTsp;發(fā)射脈沖的頻率編碼序列為fn=anΔf,n表示脈沖在脈組中的索引,Δf為頻率碼片的頻率間隔,且Δf=1/Tsp,an={a1,a2,…,aM}n為第n個(gè)脈沖頻率編碼系數(shù),且是由整數(shù){0,1,…,M-1}的亂序重排;un(t)表示第n個(gè)發(fā)射信號(hào)的脈內(nèi)調(diào)制項(xiàng)。

    文獻(xiàn)[16]通過(guò)改變信號(hào)頻率編碼獲得相互正交的DFC信號(hào)組。本文通過(guò)設(shè)定N個(gè)不同的頻率編碼,生成一組正交的DFC波形,并使用該信號(hào)組進(jìn)行循環(huán)發(fā)射。在雷達(dá)的每次回波中,包含了與目標(biāo)模糊區(qū)域?qū)?yīng)的頻率編碼信號(hào)回波,對(duì)接收到的回波信號(hào)使用各個(gè)正交波形依次進(jìn)行匹配濾波,可以準(zhǔn)確地對(duì)不同距離模糊區(qū)域的回波進(jìn)行分離。

    頻率編碼信號(hào)具有與相位編碼信號(hào)類似的圖釘型模糊函數(shù)[17],這意味著高速目標(biāo)回波的脈內(nèi)多普勒調(diào)制項(xiàng)不可忽略,由于脈內(nèi)多普勒的存在,使用發(fā)射信號(hào)設(shè)計(jì)的匹配濾波函數(shù)會(huì)引起失配,使得匹配濾波輸出峰值降低,脈沖壓縮增益損失。脈壓增益的損失可以使用模糊函數(shù)0時(shí)延(脈壓峰值)切片的多普勒響應(yīng)進(jìn)行評(píng)估,脈壓峰值的多普勒響應(yīng)可以表示為

    (9)

    式中:fd為信號(hào)的脈內(nèi)多普勒頻率;u*(t)表示共扼。將un(t)代入式(9)可得

    R(fd)=MTspsinc(fdMTsp)

    (10)

    圖2(b)為一個(gè)碼片數(shù)為20,碼片時(shí)寬為1 μs的頻率編碼信號(hào)0時(shí)間延時(shí)切片相關(guān)結(jié)果的多普勒響應(yīng)。由式(10)及圖2(b)可知,隨著脈內(nèi)多普勒頻率的增大,在回波信號(hào)使用發(fā)射信號(hào)作為匹配濾波參考函數(shù)時(shí),脈沖壓縮的峰值明顯降低,其3 dB的增益損失對(duì)應(yīng)的多普勒頻率為22.3 kHz,當(dāng)雷達(dá)波形載頻為10 GHz時(shí),對(duì)應(yīng)的速度為334.5 m/s,在機(jī)載PD火控雷達(dá)中,當(dāng)迎飛目標(biāo)相對(duì)速度大于該速度時(shí),脈沖壓縮需要考慮脈內(nèi)多普勒帶來(lái)的影響。

    2 基于DFC的兩維解模糊方法

    2.1 回波處理流程框架

    本文根據(jù)DFC波形特性,提出了一種融合了精確匹配濾波的二維解模糊方法,本文不僅可以通過(guò)波形正交性以及Keystone多普勒模糊補(bǔ)償準(zhǔn)確解決距離-多普勒二維模糊,而且將回波信號(hào)的匹配濾波放在多普勒域進(jìn)行,根據(jù)每條多普勒譜線進(jìn)行針對(duì)性的匹配濾波器設(shè)計(jì),消除空空模式可能存在的高速目標(biāo)回波的PD對(duì)脈壓的影響,從而達(dá)到最大化匹配濾波輸出。本文算法流程圖如圖3所示,以下各小節(jié)將結(jié)合算法流程圖對(duì)回波處理流程進(jìn)行詳細(xì)介紹。

    圖3 算法流程圖Fig.3 Flow diagram of algorithm

    2.2 回波信號(hào)Keystone校正及多普勒模糊補(bǔ)償

    考慮雷達(dá)對(duì)高速目標(biāo)觀測(cè)過(guò)程中,目標(biāo)回波包絡(luò)有越距離單元走動(dòng)的現(xiàn)象,在進(jìn)行相參積累時(shí),會(huì)引起距離向以及多普勒向的主瓣展寬,無(wú)法獲得最大的相參積累增益。本文使用Keystone方法,通過(guò)解決多普勒頻域與距離頻率的耦合關(guān)系,拉直信號(hào)包絡(luò),達(dá)到距離走動(dòng)校正的目的,以獲得較高的相參積累增益。且本文在多普勒域進(jìn)行匹配濾波,因此Keystone是對(duì)脈壓前的回波信號(hào)進(jìn)行處理。

    由于使用了一組正交DFC波形循環(huán)發(fā)射,因此在回波中同一目標(biāo)相鄰的回波在頻域的調(diào)制項(xiàng)上并不一致,為了保證Keystone過(guò)程中信號(hào)慢時(shí)間重采樣的正確性,需要對(duì)回波信號(hào)進(jìn)行N次抽取,以保證Keystone在相鄰脈沖之間有正確的插值信號(hào),如圖4(a)所示。

    圖4 信號(hào)抽取及多普勒模糊示意圖Fig.4 Schematic diagram of signal extract and Doppler ambiguity

    將脈內(nèi)調(diào)制項(xiàng)un(t)進(jìn)行傅里葉變換,第n個(gè)脈內(nèi)調(diào)制信號(hào)對(duì)應(yīng)的頻域形式可以表示為

    (11)

    式中:

    (12)

    根據(jù)式(6),目標(biāo)基頻回波的快時(shí)間頻域?yàn)?/p>

    (13)

    對(duì)回波基頻信號(hào)進(jìn)行抽取,信號(hào)共分為N組,各個(gè)模糊區(qū)域的回波在相鄰的脈沖間具有相同的頻率編碼,且不隨慢時(shí)間變化。Keystone可以通過(guò)插值得到正確的慢時(shí)間重采樣。由于信號(hào)抽取,各個(gè)脈組的起始距離會(huì)發(fā)生變化。設(shè)抽取后慢時(shí)間為tr=Ntm,第n個(gè)抽取分組表示為

    (14)

    對(duì)抽取后的脈沖分組在慢時(shí)間維度tr進(jìn)行重采樣,令τr=(f+fc)tr/fc,重采樣后第n個(gè)脈沖分組的回波信號(hào)為

    (15)

    對(duì)于多普勒模糊的回波信號(hào),由于式(15)中第二個(gè)指數(shù)項(xiàng)的存在,信號(hào)包絡(luò)不完全處于同一個(gè)距離單元,需要進(jìn)行補(bǔ)償。進(jìn)行準(zhǔn)確的多普勒補(bǔ)償后,信號(hào)包絡(luò)才能處于同一距離單元,在相參積累后,獲得最大的積累增益。利用這一特性,通過(guò)模糊數(shù)F1和F2的搜索遍歷,根據(jù)相參積累最大積累增益對(duì)應(yīng)的模糊數(shù),達(dá)到多普勒解模糊的目的。多普勒模糊數(shù)補(bǔ)償通過(guò)下式完成:

    (16)

    式中:F=F1N+F2為對(duì)應(yīng)的多普勒模糊數(shù)。通過(guò)遍歷一定范圍內(nèi)的多普勒模糊數(shù),并根據(jù)后續(xù)相參積累結(jié)果的峰值大小,可以找到對(duì)應(yīng)的多普勒模糊數(shù),完成多普勒解模糊。

    2.3 多普勒敏感性匹配濾波與解距離模糊

    由于DFC波形存在多普勒敏感性,本文將匹配濾波放在多普勒域進(jìn)行,針對(duì)每一個(gè)多普勒譜線進(jìn)行匹配濾波。在上一小節(jié)中,已經(jīng)針對(duì)每組波形進(jìn)行了速度模糊數(shù)估計(jì),在不同的模糊數(shù)下,每條多普勒譜線對(duì)應(yīng)不同的多普勒值。

    對(duì)式(16)進(jìn)行慢時(shí)間傅里葉變換,變換到快時(shí)間頻率-多普勒域?yàn)?/p>

    (17)

    式中:TR表示回波的相干處理時(shí)間。

    考慮多普勒模糊,上式可以重寫(xiě)為

    (18)

    式中:F為目標(biāo)的多普勒模糊數(shù),

    (19)

    (20)

    此時(shí),對(duì)回波信號(hào)快時(shí)間頻域-多普勒域信號(hào)進(jìn)行匹配濾波,同時(shí)得到距離解模糊結(jié)果。針對(duì)每一條多普勒譜線,設(shè)計(jì)匹配濾波器組:

    (21)

    式中:r∈[1,2,…,N]表示不同的脈內(nèi)調(diào)制項(xiàng),對(duì)應(yīng)于不同的距離模糊區(qū)域;(·)*表示取共軛,該匹配濾波器同時(shí)補(bǔ)償了由抽取帶來(lái)的起始距離的變化。

    (22)

    使用考慮了脈內(nèi)多普勒的匹配濾波器組,將不同距離模糊區(qū)域的信號(hào)進(jìn)行分離,并且該匹配濾波器與帶有脈內(nèi)多普勒調(diào)制的回波匹配,可以進(jìn)行精確的匹配濾波。

    2.4 多普勒譜濾波融合

    進(jìn)行精確匹配濾波后,補(bǔ)償了脈內(nèi)多普勒的影響,獲得了每個(gè)抽取后脈組的距離-多普勒結(jié)果。為了獲得更好的相參積累效果以及簡(jiǎn)化多普勒模糊數(shù)對(duì)比估計(jì),對(duì)各個(gè)抽取后的脈沖分組結(jié)果進(jìn)行融合,并在融合過(guò)程補(bǔ)償?shù)粲沙槿?lái)的多普勒模糊問(wèn)題,此時(shí)只需要估計(jì)目標(biāo)多普勒相對(duì)于fp的多普勒模糊數(shù)。

    由于發(fā)射脈沖是在固定重復(fù)頻率下發(fā)射的,因此各個(gè)發(fā)射脈沖之間的時(shí)間間隔相同,在一個(gè)發(fā)射脈組內(nèi),多普勒采樣是等間隔連續(xù)的,利用這一特性,可以進(jìn)行多普勒濾波重構(gòu),如圖5所示。利用多個(gè)抽取脈組的距離-多普勒結(jié)果,根據(jù)各個(gè)抽取脈組之間該多普勒頻點(diǎn)之間的關(guān)系(式(22)最后一個(gè)相位項(xiàng)),重構(gòu)出目標(biāo)真實(shí)多普勒,解決由抽取帶來(lái)的目標(biāo)多普勒模糊。

    圖5 多普勒模糊及融合示意圖Fig.5 Doppler ambiguity and fusion diagram

    A=[a-I,a-I+1,…,ai,…,aI-1,aI]

    (23)

    (24)

    (25)

    令q為不同分組信號(hào)所具有的相同的距離像包絡(luò),p=[p-I,p-I+1,…,pF,…,pI-1,p-I],其中多普勒模糊數(shù)對(duì)應(yīng)的pF為1,其他為0。

    (26)

    式(25)可以表示為

    SB=(ApTq)T

    (27)

    (28)

    式中:Wi為第i次多普勒模糊的權(quán)重系數(shù),且Wi=A-1hT;h=[h-I,h-I+1,…,hi,…h(huán)I]為各個(gè)模糊分量的選擇向量,令hi=1,其他為0。

    將各個(gè)多普勒頻點(diǎn)在不同多普勒模糊次數(shù)的結(jié)果,進(jìn)行結(jié)果的拼接即可得到最終的融合后距離-多普勒結(jié)果:

    (29)

    此步融合了整個(gè)抽取分組的能量,即不同距離模糊區(qū)域的多個(gè)分組也完成了相參積累,在此結(jié)果之上可以進(jìn)行目標(biāo)檢測(cè)以及參數(shù)估計(jì)。

    2.5 目標(biāo)檢測(cè)及多普勒模糊數(shù)確定

    在經(jīng)過(guò)以上步驟后,得到了回波在不同距離模糊區(qū)域以及不同多普勒模糊單元的距離-多普勒結(jié)果,對(duì)該結(jié)果進(jìn)行二維恒虛警率(constant false-alarm rate, CFAR)檢測(cè)[18]及參數(shù)估計(jì)。由于距離模糊根據(jù)信號(hào)正交性已經(jīng)分離,此處檢測(cè)結(jié)果的距離信息加上距離模糊區(qū)域?qū)?yīng)的起始距離即為真實(shí)距離信息。對(duì)于多普勒模糊,當(dāng)信噪比較高時(shí),多普勒模糊數(shù)錯(cuò)誤補(bǔ)償?shù)姆e累記過(guò)也會(huì)被檢測(cè)出來(lái),其多普勒中心與最優(yōu)相參積累結(jié)果中心一致,利用不同多普勒模糊數(shù)補(bǔ)償后的檢測(cè)結(jié)果,對(duì)比不同模糊次數(shù)下多普勒頻點(diǎn)的相參積累輸出,找出最大值對(duì)應(yīng)的多普勒模糊次數(shù)以及多普勒值,即可得到目標(biāo)對(duì)應(yīng)多普勒以及目標(biāo)真實(shí)速度。對(duì)目標(biāo)檢測(cè)結(jié)果進(jìn)行進(jìn)一步處理,得到目標(biāo)的參數(shù)估計(jì)結(jié)果進(jìn)行輸出。

    2.6 算法復(fù)雜度分析

    為了便于衡量算法的工程實(shí)現(xiàn)可行性,對(duì)算法復(fù)雜度進(jìn)行分析。假設(shè)一次處理內(nèi),共發(fā)射Na個(gè)脈組,共計(jì)NaN個(gè)脈沖,每個(gè)脈沖內(nèi)的信號(hào)采樣點(diǎn)為Nr,進(jìn)行L次多普勒模糊遍歷。在算法流程中,信號(hào)分組Keystone變換,算法運(yùn)算量為NaNNr(1+log2Nr),多普勒模糊數(shù)補(bǔ)償中,算法運(yùn)算量為NrNaNL,脈沖壓縮過(guò)程算法運(yùn)算量為NrNaN2L(1+log2Nr),多普勒譜融合的運(yùn)算量為NrNaNL,實(shí)際算法中,解速度模糊次數(shù)L與脈組脈沖個(gè)數(shù)相同,即L=N,算法流程的總運(yùn)算量為

    OP=NaNNr((1+N2)(1+log2Nr)+2N)≈NaNrN3log2Nr

    (30)

    由式(30)可見(jiàn),當(dāng)脈組中脈沖個(gè)數(shù)N=1時(shí),運(yùn)算量最小為O(NaNrlog2Nr),隨著N增大,運(yùn)算量增大為O(NaNrN3log2Nr)。N增大不僅代表解距離模糊能力的提升,同時(shí)也增加了相參處理脈沖個(gè)數(shù),提高多普勒分辨與測(cè)量能力。

    3 仿真實(shí)驗(yàn)及性能分析

    為了證明利用循環(huán)發(fā)射DFC波形進(jìn)行脈沖多普勒雷達(dá)在MPRF模式下對(duì)距離-多普勒二維解模糊以及低信噪比回波的目標(biāo)檢測(cè)能力,進(jìn)行了以下兩部分實(shí)驗(yàn)進(jìn)行驗(yàn)證。① 使用正交性匹配濾波以及Keystone多普勒模糊校正,驗(yàn)證其可以準(zhǔn)確解決距離-多普勒二維模糊;② 驗(yàn)證考慮多普勒敏感性的精確匹配濾波,應(yīng)對(duì)不同速度目標(biāo)的脈壓結(jié)果輸出、固定重頻帶來(lái)的整個(gè)駐留時(shí)間及長(zhǎng)時(shí)間相參積累帶來(lái)的信噪比增益,以及獲得信噪比增益后目標(biāo)檢測(cè)和參數(shù)估計(jì)的優(yōu)勢(shì)。

    在所進(jìn)行的實(shí)驗(yàn)中,設(shè)定機(jī)載火控PD雷達(dá)循環(huán)發(fā)射5個(gè)正交DFC波形脈沖,可以進(jìn)行最大4次的距離模糊分離。具體的雷達(dá)仿真參數(shù)如表1所示。按照信號(hào)正交性為優(yōu)化目標(biāo)函數(shù),使用遺傳算法進(jìn)行頻率編碼的優(yōu)化,得到發(fā)射的5個(gè)頻率編碼脈沖對(duì)應(yīng)頻率編碼如表2所示。

    表1 仿真參數(shù)

    表2 頻率編碼

    3.1 DFC解距離-多普勒兩維模糊

    首先通過(guò)實(shí)驗(yàn)驗(yàn)證利用正交波形進(jìn)行距離模糊區(qū)域區(qū)分的有效性,根據(jù)仿真參數(shù),雷達(dá)的最大不模糊距離為15 km。在本實(shí)驗(yàn)中,設(shè)置3個(gè)分別位于不同距離模糊區(qū)域的目標(biāo),目標(biāo)距離分別為7.8 km、27 km、33 km。由于對(duì)信號(hào)的匹配濾波在多普勒域中完成,因此需要對(duì)比相同多普勒目標(biāo)的正交波形距離模糊區(qū)分結(jié)果,此處設(shè)定3個(gè)目標(biāo)速度均為0。在一個(gè)波位駐留對(duì)應(yīng)的相干處理時(shí)間內(nèi),共循環(huán)發(fā)射了256組信號(hào),共1 280個(gè)脈沖。

    圖6為目標(biāo)所在多普勒譜線的匹配濾波結(jié)果,由圖中可以發(fā)現(xiàn),對(duì)于在同一多普勒譜線的不同距離模糊區(qū)域目標(biāo),由于不同模糊區(qū)域具有對(duì)應(yīng)的正交信號(hào),只在相對(duì)應(yīng)的距離模糊區(qū)域被正確的脈沖壓縮,在其他區(qū)域由于匹配濾波器參考信號(hào)的正交,不產(chǎn)生峰值。圖6(a)為距離不模糊區(qū)域的匹配濾波結(jié)果,信號(hào)輸出主瓣附近為-32 dB左右的信號(hào)自相關(guān)旁瓣。位于目標(biāo)距離模糊對(duì)應(yīng)的區(qū)域,存在距離模糊區(qū)域回波信號(hào)與此區(qū)域編碼信號(hào)的相關(guān)結(jié)果輸出,由于信號(hào)的正交性,匹配濾波的峰值的主瓣旁瓣比大于22.87 dB,驗(yàn)證了發(fā)射波形的正交性,以及由正交性帶來(lái)的距離模糊區(qū)域的區(qū)分。第3距離模糊區(qū)域無(wú)目標(biāo),對(duì)應(yīng)的匹配濾波結(jié)果低于-22.55 dB,為信號(hào)旁瓣量級(jí)。

    以上證明使用正交DFC波形可以實(shí)現(xiàn)在一定模糊旁瓣抑制的條件下,進(jìn)行不同距離模糊區(qū)域的分離,接下來(lái)通過(guò)實(shí)驗(yàn)驗(yàn)證多普勒解模糊的有效性。

    設(shè)定目標(biāo)距離為7.8 km的目標(biāo)速度為120 m/s,相對(duì)于最大不模糊速度來(lái)說(shuō),有了1次多普勒模糊,在回波信號(hào)進(jìn)行Keystone模糊數(shù)搜索中,對(duì)不同模糊數(shù)下距離-多普勒結(jié)果中的值進(jìn)行對(duì)比,結(jié)果如圖7所示,進(jìn)行多普勒模糊數(shù)為1的補(bǔ)償后,相參積累結(jié)果具有高的峰值能量且主瓣較窄,其他模糊數(shù)補(bǔ)償?shù)南鄥⒎e累結(jié)果由于目標(biāo)距離像分布在多個(gè)距離單元,相參處理后能量未完全積累在同一個(gè)距離單元,峰值較低,且主瓣寬,影響距離分辨能力。表3為不同多普勒模糊數(shù)補(bǔ)償下,積累峰值的輸出。通過(guò)不同模糊數(shù)補(bǔ)償?shù)姆e累結(jié)果對(duì)比,可以確定回波信號(hào)的多普勒模糊數(shù)為1。目標(biāo)在距離多普勒結(jié)果中,多普勒譜線對(duì)應(yīng)的速度為-30.082 1 m/s,結(jié)合峰值比較輸出得到的目標(biāo)多普勒模糊數(shù),以及1次多普勒模糊對(duì)應(yīng)150 m/s,可以求得目標(biāo)的估計(jì)速度為119.917 9 m/s,與設(shè)定目標(biāo)基本一致。

    圖7 多普勒解模糊結(jié)果Fig.7 Result of doppler ambiguity resolution

    表3 不同多普勒模糊數(shù)補(bǔ)償?shù)姆逯?/p>

    3.2 信號(hào)能量及目標(biāo)檢測(cè)性能

    從能量以及DFC信號(hào)多普勒敏感性考慮出發(fā),使用精確匹配濾波可以避免由匹配濾波函數(shù)不匹配引入的失配脈壓增益損失。本節(jié)首先通過(guò)實(shí)驗(yàn)對(duì)比考慮脈內(nèi)多普勒的匹配濾波與不考慮脈內(nèi)多普勒的匹配濾波在不同速度目標(biāo)下的匹配濾波輸出。

    為了驗(yàn)證不同脈內(nèi)多普勒的影響,實(shí)驗(yàn)設(shè)置目標(biāo)距離為7 800 m,目標(biāo)速度范圍為0~1 800 m/s,分別進(jìn)行不同目標(biāo)速度的回波仿真,在多普勒域根據(jù)每條多普勒譜線以及多普勒模糊次數(shù),設(shè)計(jì)相應(yīng)的匹配濾波器。為了進(jìn)行對(duì)比,使用不考慮多普勒敏感性的匹配濾波對(duì)各多普勒單元進(jìn)行匹配濾波,結(jié)果如圖8(a)所示,從該結(jié)果可以明顯看出,隨著目標(biāo)速度的增加,匹配濾波結(jié)果輸出的主瓣峰值降低,旁瓣明顯升高,同時(shí)速度越大,匹配濾波結(jié)果出現(xiàn)明顯失真,無(wú)法獲得匹配濾波增益,對(duì)目標(biāo)的檢測(cè)造成嚴(yán)重影響。圖8(b)為不同速度下,使用本文方法的目標(biāo)匹配濾波輸出,由結(jié)果可見(jiàn),匹配濾波結(jié)果的峰值沒(méi)有明顯展寬,匹配濾波輸出的峰值基本保持不變,并且脈壓輸出的副瓣未產(chǎn)生明顯升高。為了更加清晰地表示考慮多普勒敏感性匹配濾波的性能,如圖8(c)和圖8(d)所示,計(jì)算了兩種匹配濾波方法在不同速度下的歸一化峰值輸出以及主瓣旁瓣比。如圖8(c)所示,考慮了多普勒敏感性的匹配濾波輸出主瓣峰值在不同速度下基本不變,但是不考慮多普勒敏感性的匹配濾波輸出峰值隨著速度升高持續(xù)下降,在1 500 m/s速度下、已經(jīng)降低了13 dB,嚴(yán)重影響了目標(biāo)的檢測(cè)。如圖8(d)所示,考慮了多普勒敏感性的匹配濾波輸出的主副比,在1 500 m/s相對(duì)速度下依舊可以保持25 dB以上的主副比,而傳統(tǒng)直接利用發(fā)射信號(hào)構(gòu)造匹配濾波的結(jié)果,隨著速度地增加,匹配濾波輸出的主副比明顯降低,在1 200 m/s的相對(duì)速度下主副比為8.76 dB,已經(jīng)嚴(yán)重影響了匹配濾波的輸出,極易造成目標(biāo)的漏檢或虛警。

    圖8 脈內(nèi)多普勒對(duì)匹配濾波的影響Fig.8 Influence of Doppler sensitivity on matched filtering

    為了驗(yàn)證本文所提算法對(duì)目標(biāo)檢測(cè)以及參數(shù)估計(jì)性能,與經(jīng)典的多脈沖重頻法[13]、脈沖間相位編碼[5]以及DFC信號(hào)不考慮多普勒敏感性匹配濾波的方法進(jìn)行對(duì)比。多脈沖重頻法使用3/5原則,在一次駐留時(shí)間內(nèi),發(fā)射5組不同PRF的脈沖,具體PRF分別為15 kHz、8.7 kHz、7.2 kHz、11 kHz、13 kHz,每組PRF分別包含256個(gè)脈沖,信號(hào)帶寬設(shè)置為與本文方法仿真一致。使用文獻(xiàn)[5]文中13脈沖的方法進(jìn)行對(duì)比,信號(hào)重頻同樣為10 kHz,信號(hào)帶寬為20 MHz,每次駐留時(shí)間內(nèi)共有100組脈沖簇。相對(duì)應(yīng)的,在實(shí)驗(yàn)對(duì)比中,本文所提算法在每次駐留時(shí)間內(nèi)進(jìn)行256個(gè)脈組信號(hào)發(fā)射。用于對(duì)比的DFC信號(hào)不考慮多普勒敏感性匹配濾波方法,除了在匹配濾波不考慮脈內(nèi)多普勒外,其他處理流程均與所提方法一致。

    實(shí)驗(yàn)中,設(shè)置起始距離為7 800 m的目標(biāo)速度為330 m/s,對(duì)仿真目標(biāo)脈沖回波加入不同信噪比的高斯白噪聲,通過(guò)150次蒙特卡羅實(shí)驗(yàn)進(jìn)行目標(biāo)檢測(cè)以及測(cè)量精度實(shí)驗(yàn)。實(shí)驗(yàn)結(jié)果如圖9所示,未多普勒補(bǔ)償DFC法表示未考慮多普勒敏感性的DFC波形傳統(tǒng)匹配濾波方法的結(jié)果。從圖9(a)可以看到在每種信噪比下進(jìn)行150次蒙特卡羅實(shí)驗(yàn)后不同算法的檢測(cè)概率。多脈沖重頻法和脈間相位編碼法的結(jié)果分別在信噪比為-38 dB和-36 dB時(shí),目標(biāo)的檢測(cè)概率已經(jīng)趨近于0。對(duì)于未多普勒補(bǔ)償DFC法,從圖8(c)中可以看出,在目標(biāo)速度為330 m/s時(shí),由于匹配濾波的失配,脈壓結(jié)果的峰值輸出相對(duì)于考慮多普勒敏感性匹配濾波的結(jié)果下降了6 dB左右,峰值下降使得目標(biāo)在回波信噪比為-40 dB時(shí),只有24%左右的檢測(cè)概率,在信噪比為-38 dB的情況下,檢測(cè)概率為64%左右。由于本文方法在整個(gè)駐留時(shí)間的脈沖回波進(jìn)行了相參積累,且考慮了脈內(nèi)多普勒對(duì)脈壓的影響,因此在回波信噪比為-44 dB的情況下,依舊有60%左右的檢測(cè)概率。從結(jié)果可以看出,本文所提算法在低信噪比目標(biāo)回波中的檢測(cè)性能明顯高于其他3種方法。低信噪比回波下的目標(biāo)檢測(cè)能力可以使雷達(dá)對(duì)相同距離目標(biāo)探測(cè)時(shí),在保持相同檢測(cè)概率的情況下降低雷達(dá)發(fā)射功率,降低雷達(dá)信號(hào)被截獲概率,增加雷達(dá)的射頻隱身性能。由文獻(xiàn)[13]可知,多脈沖重復(fù)頻率方法在進(jìn)行目標(biāo)匹配時(shí),容易產(chǎn)生虛假目標(biāo)點(diǎn),增大虛警率。圖9(b)為不同信噪比下目標(biāo)虛警率的對(duì)比結(jié)果,多脈沖重頻法具有約5%的匹配虛假點(diǎn),其他算法的目標(biāo)檢測(cè)虛警率接近0。本文所提處理方法有利于獲得較高信噪比增益,使目標(biāo)的測(cè)量結(jié)果更加精確。圖9(c)為不同算法的目標(biāo)距離測(cè)量精度,使用-30~0 dB的回波信噪比,以保證在各個(gè)算法對(duì)目標(biāo)穩(wěn)定檢測(cè)下的測(cè)量精度的對(duì)比。從結(jié)果中可以看出,本文算法在不同信噪比下,目標(biāo)距離測(cè)量精度可以保持在0.05 m以內(nèi),多脈沖重頻法和脈間相位編碼法方法的測(cè)量精度隨著信噪比的提高也有提高,但是仍然大于本文所提方法。不考慮多普勒敏感性的DFC波形方法由于脈內(nèi)多普勒引起目標(biāo)主瓣展寬,使得目標(biāo)距離分辨變差,且在CFAR固定檢測(cè)門(mén)限下,在回波高信噪比情況下更容易檢測(cè)到展寬的主瓣,使得目標(biāo)的距離測(cè)量結(jié)果變差。對(duì)于目標(biāo)速度的測(cè)量,由于目標(biāo)存在多普勒模糊,脈間相位編碼法的方法沒(méi)有解多普勒模糊的能力,因此對(duì)于目標(biāo)速度的提取錯(cuò)誤,在速度精度結(jié)果中不再體現(xiàn)錯(cuò)誤的結(jié)果。如圖9(d)所示,由于本文所提算法使用了整個(gè)駐留時(shí)間的脈沖進(jìn)行多普勒處理,多普勒分辨率明顯高于多脈沖重頻法的方法,而且在信噪比處理增益上也有提高。因此對(duì)于速度的測(cè)量精度,本文方法明顯高于多脈沖重頻法方法。對(duì)于不考慮多普勒敏感性的DFC波形方法,速度測(cè)量精度與本文所提方法相差不大,但是由于匹配濾波失配引起的信噪比損失影響了測(cè)量精度,使得其在不同信噪比下的起伏高于本文所提方法。

    圖9 目標(biāo)檢測(cè)及參數(shù)估計(jì)性能Fig.9 Target detection and parameter estimation performance

    表4為不同方法對(duì)于距離和速度均方根誤差(root mean sguare error, RMSE)的方差。從結(jié)果可以發(fā)現(xiàn),本文方法對(duì)于目標(biāo)距離和速度測(cè)量結(jié)果的方法均優(yōu)于其他方法。

    表4 目標(biāo)測(cè)量結(jié)果RMSE方差

    從以上仿真實(shí)驗(yàn)結(jié)果可以看出,本文所提方法可以準(zhǔn)確解決PD雷達(dá)在MPRF下的距離-多普勒二維模糊問(wèn)題,同時(shí)通過(guò)長(zhǎng)時(shí)間相參積累以及精確匹配濾波,可以顯著提高目標(biāo)的檢測(cè)性能以及參數(shù)測(cè)量精度。

    4 結(jié) 論

    機(jī)載脈沖多普勒雷達(dá)在MPRF工作模式下,目標(biāo)回波可能存在距離-多普勒二維模糊,本文提出了循環(huán)發(fā)射正交DFC波形脈沖的方法,利用信號(hào)的正交性,在匹配濾波階段進(jìn)行距離模糊區(qū)域的區(qū)分,通過(guò)Keystone處理中搜索最佳相參積累增益下的補(bǔ)償多普勒模糊數(shù),實(shí)現(xiàn)多普勒解模糊,同時(shí)Keystone帶來(lái)整個(gè)駐留時(shí)間的相參積累,增加了在低信噪比回波情況下的檢測(cè)性能。為了解決相對(duì)速度較快帶來(lái)的多普勒敏感性問(wèn)題,采用了多普勒域匹配濾波的方法,通過(guò)濾波器的設(shè)計(jì),可以準(zhǔn)確有效解決多普勒敏感性問(wèn)題。仿真實(shí)驗(yàn)驗(yàn)證了算法的有效性。

    猜你喜歡
    信噪比多普勒波形
    對(duì)《壓力容器波形膨脹節(jié)》2018版新標(biāo)準(zhǔn)的理解及分析
    基于深度學(xué)習(xí)的無(wú)人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    基于LFM波形的靈巧干擾效能分析
    低信噪比下LFMCW信號(hào)調(diào)頻參數(shù)估計(jì)
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    基于ARM的任意波形電源設(shè)計(jì)
    基于多普勒效應(yīng)的車隨人動(dòng)系統(tǒng)
    電子器件(2015年5期)2015-12-29 08:43:38
    大連臺(tái)使用CTS-1記錄波形特點(diǎn)
    保持信噪比的相位分解反褶積方法研究
    基于多普勒的車輛測(cè)速儀
    av有码第一页| 欧美bdsm另类| av一本久久久久| 精品久久久精品久久久| 久久久久久人妻| 国产成人一区二区在线| 国产精品一国产av| 国产高清国产精品国产三级| 欧美精品亚洲一区二区| 少妇的逼水好多| videosex国产| 91精品三级在线观看| 最近2019中文字幕mv第一页| 亚洲色图 男人天堂 中文字幕| 91国产中文字幕| 久久免费观看电影| 青草久久国产| 爱豆传媒免费全集在线观看| 免费av中文字幕在线| 男女边吃奶边做爰视频| 中国三级夫妇交换| 亚洲欧洲精品一区二区精品久久久 | 免费av中文字幕在线| 久久这里只有精品19| 精品一区在线观看国产| 夜夜骑夜夜射夜夜干| 黄色视频在线播放观看不卡| 婷婷成人精品国产| 国产精品麻豆人妻色哟哟久久| 男女无遮挡免费网站观看| 日本欧美视频一区| 成年女人在线观看亚洲视频| 国产精品一区二区在线观看99| 国产精品国产三级专区第一集| 国产 一区精品| 午夜激情久久久久久久| 王馨瑶露胸无遮挡在线观看| 免费黄色在线免费观看| 久久久久久免费高清国产稀缺| 欧美97在线视频| 亚洲第一青青草原| 1024香蕉在线观看| 精品一品国产午夜福利视频| 国产又色又爽无遮挡免| 亚洲精品美女久久久久99蜜臀 | 国产成人av激情在线播放| 精品国产国语对白av| 我要看黄色一级片免费的| 国产黄色视频一区二区在线观看| 超碰97精品在线观看| 18禁国产床啪视频网站| 亚洲精品日韩在线中文字幕| 亚洲国产精品一区三区| 九色亚洲精品在线播放| 高清视频免费观看一区二区| 成人手机av| 国产日韩一区二区三区精品不卡| 女人被躁到高潮嗷嗷叫费观| 深夜精品福利| 久久精品国产自在天天线| 在线观看三级黄色| 久久久久网色| videos熟女内射| 黄色一级大片看看| 18禁裸乳无遮挡动漫免费视频| 亚洲综合精品二区| 国产日韩欧美在线精品| 亚洲av中文av极速乱| 久久ye,这里只有精品| 最近的中文字幕免费完整| 日韩精品免费视频一区二区三区| 久久99精品国语久久久| 啦啦啦中文免费视频观看日本| 欧美成人精品欧美一级黄| 精品人妻熟女毛片av久久网站| 香蕉国产在线看| 久久久久国产一级毛片高清牌| 又大又黄又爽视频免费| 天天躁狠狠躁夜夜躁狠狠躁| 欧美在线黄色| 亚洲第一区二区三区不卡| 亚洲成国产人片在线观看| 9191精品国产免费久久| 国产亚洲av片在线观看秒播厂| 国产av精品麻豆| 成人国产麻豆网| 黄色毛片三级朝国网站| 久久精品国产鲁丝片午夜精品| 国产高清不卡午夜福利| 亚洲精品视频女| 国产成人91sexporn| 中文字幕色久视频| 我的亚洲天堂| 卡戴珊不雅视频在线播放| 一本久久精品| 国产日韩欧美视频二区| 丁香六月天网| 亚洲久久久国产精品| 制服丝袜香蕉在线| 精品国产一区二区久久| 国语对白做爰xxxⅹ性视频网站| 极品少妇高潮喷水抽搐| 最近最新中文字幕大全免费视频 | 欧美激情高清一区二区三区 | 欧美国产精品va在线观看不卡| 在线观看一区二区三区激情| 亚洲精品自拍成人| av国产精品久久久久影院| 久久久久久人妻| 精品国产一区二区三区久久久樱花| 少妇的丰满在线观看| 久久久久网色| 精品人妻熟女毛片av久久网站| 最新中文字幕久久久久| 亚洲精品日本国产第一区| 国产午夜精品一二区理论片| 国产精品成人在线| 亚洲 欧美一区二区三区| 亚洲精品一区蜜桃| 国产福利在线免费观看视频| 观看美女的网站| 国产深夜福利视频在线观看| 人成视频在线观看免费观看| 久久久久久久精品精品| 久久久国产欧美日韩av| 午夜日本视频在线| 久久精品久久精品一区二区三区| 波多野结衣av一区二区av| 精品福利永久在线观看| 黄色怎么调成土黄色| 久久人人97超碰香蕉20202| 一级毛片黄色毛片免费观看视频| 丰满乱子伦码专区| 中文字幕亚洲精品专区| 国产免费福利视频在线观看| 国产精品女同一区二区软件| 少妇精品久久久久久久| 街头女战士在线观看网站| 国产1区2区3区精品| 十八禁高潮呻吟视频| 国产成人精品久久二区二区91 | 毛片一级片免费看久久久久| 久久久久久久亚洲中文字幕| 免费人妻精品一区二区三区视频| 欧美精品人与动牲交sv欧美| 国产综合精华液| 999久久久国产精品视频| 欧美最新免费一区二区三区| 国产亚洲午夜精品一区二区久久| 伊人久久大香线蕉亚洲五| 在线天堂最新版资源| av网站在线播放免费| 秋霞在线观看毛片| 久久久久久久久久人人人人人人| 男女国产视频网站| 黄网站色视频无遮挡免费观看| 久久鲁丝午夜福利片| 日韩一卡2卡3卡4卡2021年| 国产精品嫩草影院av在线观看| a级片在线免费高清观看视频| av福利片在线| 亚洲国产色片| 久久久久精品人妻al黑| 国产成人精品婷婷| 啦啦啦啦在线视频资源| 久久人妻熟女aⅴ| 亚洲欧洲精品一区二区精品久久久 | 又黄又粗又硬又大视频| 欧美黄色片欧美黄色片| 精品久久久久久电影网| 最近2019中文字幕mv第一页| www.av在线官网国产| 一本—道久久a久久精品蜜桃钙片| 97人妻天天添夜夜摸| 亚洲视频免费观看视频| 桃花免费在线播放| 亚洲精品成人av观看孕妇| 美女高潮到喷水免费观看| 亚洲精品乱久久久久久| 国产一级毛片在线| 久久免费观看电影| 亚洲欧美一区二区三区久久| 欧美人与善性xxx| 久久久精品国产亚洲av高清涩受| 最近的中文字幕免费完整| 伦理电影大哥的女人| av又黄又爽大尺度在线免费看| 欧美成人午夜精品| 亚洲av日韩在线播放| 美女xxoo啪啪120秒动态图| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 日韩成人av中文字幕在线观看| 国产乱人偷精品视频| 久久免费观看电影| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看| 多毛熟女@视频| 18禁动态无遮挡网站| 最近最新中文字幕大全免费视频 | 国产亚洲最大av| 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人 | 免费高清在线观看日韩| 日本猛色少妇xxxxx猛交久久| 丝袜在线中文字幕| 久久久久久久久久久久大奶| 王馨瑶露胸无遮挡在线观看| 在线精品无人区一区二区三| 国产男女内射视频| 乱人伦中国视频| 美女大奶头黄色视频| 亚洲av电影在线观看一区二区三区| 99热全是精品| 秋霞伦理黄片| 熟女电影av网| 少妇人妻 视频| 亚洲国产精品成人久久小说| 大陆偷拍与自拍| 91久久精品国产一区二区三区| 久久久久久人妻| 亚洲精品视频女| 精品国产国语对白av| 国产成人午夜福利电影在线观看| 亚洲成国产人片在线观看| 天堂俺去俺来也www色官网| 国产有黄有色有爽视频| 亚洲色图综合在线观看| 国产精品国产av在线观看| 久久这里只有精品19| 在线免费观看不下载黄p国产| 90打野战视频偷拍视频| 国产综合精华液| 天天影视国产精品| 少妇人妻久久综合中文| 国产成人欧美| 国产极品天堂在线| 性色avwww在线观看| 尾随美女入室| av网站在线播放免费| 伊人久久大香线蕉亚洲五| 麻豆av在线久日| 久久精品国产综合久久久| 国产成人一区二区在线| 激情五月婷婷亚洲| 欧美97在线视频| 久久久久精品久久久久真实原创| 蜜桃国产av成人99| 亚洲精品日韩在线中文字幕| 丝袜脚勾引网站| 女人久久www免费人成看片| 欧美中文综合在线视频| 天天影视国产精品| 七月丁香在线播放| 日韩制服丝袜自拍偷拍| 欧美人与性动交α欧美精品济南到 | 在线观看人妻少妇| 成人黄色视频免费在线看| 国产 精品1| 亚洲成人av在线免费| 亚洲情色 制服丝袜| 欧美精品高潮呻吟av久久| 男人操女人黄网站| 国产精品国产三级国产专区5o| 久久人妻熟女aⅴ| 国产1区2区3区精品| 国产不卡av网站在线观看| 中国三级夫妇交换| 亚洲国产精品国产精品| 日韩伦理黄色片| 少妇熟女欧美另类| 欧美精品人与动牲交sv欧美| 老司机影院成人| 岛国毛片在线播放| 日本av手机在线免费观看| 国产淫语在线视频| 大码成人一级视频| 精品人妻一区二区三区麻豆| 亚洲欧洲国产日韩| 国产精品蜜桃在线观看| 国产精品偷伦视频观看了| 在线 av 中文字幕| 欧美激情高清一区二区三区 | 国产在线一区二区三区精| 免费人妻精品一区二区三区视频| 侵犯人妻中文字幕一二三四区| 美女脱内裤让男人舔精品视频| 精品少妇久久久久久888优播| 两个人免费观看高清视频| av免费在线看不卡| 国产在线一区二区三区精| 一区在线观看完整版| 少妇人妻精品综合一区二区| 母亲3免费完整高清在线观看 | 美女高潮到喷水免费观看| 制服人妻中文乱码| 蜜桃在线观看..| 国产成人精品婷婷| 国产xxxxx性猛交| 国产一区二区在线观看av| 欧美人与性动交α欧美软件| 精品国产一区二区三区四区第35| 一级毛片 在线播放| 晚上一个人看的免费电影| 亚洲av国产av综合av卡| 在线 av 中文字幕| 欧美激情高清一区二区三区 | 亚洲伊人色综图| 亚洲色图 男人天堂 中文字幕| 久久精品夜色国产| av.在线天堂| 大片免费播放器 马上看| 国产一区二区三区综合在线观看| 人人妻人人添人人爽欧美一区卜| 人人妻人人爽人人添夜夜欢视频| 美女xxoo啪啪120秒动态图| h视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 精品99又大又爽又粗少妇毛片| 国产精品 国内视频| 制服人妻中文乱码| av线在线观看网站| 电影成人av| 天堂8中文在线网| 18禁观看日本| 肉色欧美久久久久久久蜜桃| 亚洲精品乱久久久久久| 国产黄色视频一区二区在线观看| 国产精品二区激情视频| 五月伊人婷婷丁香| 亚洲精品国产av成人精品| 91成人精品电影| 天天躁夜夜躁狠狠久久av| 国产成人av激情在线播放| 亚洲伊人色综图| 人成视频在线观看免费观看| 最近中文字幕高清免费大全6| 国产深夜福利视频在线观看| 久久久久久久久久久久大奶| 一本久久精品| 叶爱在线成人免费视频播放| 夜夜骑夜夜射夜夜干| 欧美日韩综合久久久久久| 亚洲精品av麻豆狂野| 不卡视频在线观看欧美| 看免费av毛片| 女性生殖器流出的白浆| av卡一久久| 美女高潮到喷水免费观看| 国精品久久久久久国模美| 亚洲在久久综合| 国产成人av激情在线播放| 18禁国产床啪视频网站| 国产综合精华液| 亚洲国产精品999| 搡老乐熟女国产| 80岁老熟妇乱子伦牲交| 亚洲av电影在线观看一区二区三区| 少妇被粗大的猛进出69影院| 丰满乱子伦码专区| 欧美黄色片欧美黄色片| 精品少妇久久久久久888优播| 亚洲成av片中文字幕在线观看 | 男女无遮挡免费网站观看| 日韩中文字幕视频在线看片| 国产色婷婷99| 精品午夜福利在线看| 国产一区有黄有色的免费视频| 成人国语在线视频| 国产日韩欧美在线精品| 高清av免费在线| 精品国产一区二区三区四区第35| 欧美人与性动交α欧美精品济南到 | 国精品久久久久久国模美| 在线亚洲精品国产二区图片欧美| 在线看a的网站| 人妻人人澡人人爽人人| 在线看a的网站| 色婷婷久久久亚洲欧美| 精品人妻偷拍中文字幕| 久久精品夜色国产| freevideosex欧美| 日韩中文字幕欧美一区二区 | 成人毛片a级毛片在线播放| 亚洲四区av| 久久99一区二区三区| 亚洲国产精品成人久久小说| 亚洲激情五月婷婷啪啪| 久久久精品区二区三区| 观看av在线不卡| 免费黄色在线免费观看| 久久久久久久久久久免费av| 中文字幕最新亚洲高清| 午夜福利视频在线观看免费| 国产精品秋霞免费鲁丝片| 熟妇人妻不卡中文字幕| 亚洲av男天堂| 巨乳人妻的诱惑在线观看| 亚洲内射少妇av| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲日产国产| 赤兔流量卡办理| 国产精品一区二区在线不卡| 欧美+日韩+精品| 日韩成人av中文字幕在线观看| 日本黄色日本黄色录像| 99久久综合免费| 欧美最新免费一区二区三区| 热99国产精品久久久久久7| 最近最新中文字幕大全免费视频 | 欧美人与性动交α欧美软件| 中文天堂在线官网| 黄网站色视频无遮挡免费观看| 两性夫妻黄色片| 亚洲国产成人一精品久久久| 亚洲伊人久久精品综合| 午夜影院在线不卡| av一本久久久久| 亚洲,欧美,日韩| 黑人欧美特级aaaaaa片| 黄色配什么色好看| 丰满少妇做爰视频| 中文天堂在线官网| 久久久精品区二区三区| 亚洲欧美成人综合另类久久久| 欧美成人午夜免费资源| 成年女人毛片免费观看观看9 | 亚洲五月色婷婷综合| 亚洲精品,欧美精品| 国产一级毛片在线| 美女国产视频在线观看| 亚洲第一av免费看| 久久ye,这里只有精品| 国产探花极品一区二区| 精品国产国语对白av| 黄频高清免费视频| 人妻少妇偷人精品九色| 久久精品亚洲av国产电影网| 18禁观看日本| 日韩三级伦理在线观看| 18+在线观看网站| 美女xxoo啪啪120秒动态图| 欧美激情极品国产一区二区三区| 丰满少妇做爰视频| av网站免费在线观看视频| 在线免费观看不下载黄p国产| 制服诱惑二区| 国产精品一二三区在线看| 日韩人妻精品一区2区三区| 99热全是精品| 欧美bdsm另类| 一级爰片在线观看| 欧美精品av麻豆av| 精品酒店卫生间| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看 | 老鸭窝网址在线观看| 99热全是精品| 人人妻人人爽人人添夜夜欢视频| 99久久综合免费| 亚洲国产精品一区二区三区在线| 国产日韩欧美视频二区| 免费在线观看黄色视频的| 国产在线一区二区三区精| 少妇的丰满在线观看| 人妻 亚洲 视频| 国产精品女同一区二区软件| 国产精品一区二区在线观看99| 老司机亚洲免费影院| 日韩一本色道免费dvd| 国产成人欧美| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| 在线观看美女被高潮喷水网站| 亚洲经典国产精华液单| 美女国产高潮福利片在线看| 欧美日韩视频精品一区| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 午夜免费男女啪啪视频观看| 国产成人免费无遮挡视频| 丝袜美腿诱惑在线| 韩国av在线不卡| 久久精品久久久久久噜噜老黄| 午夜影院在线不卡| 国产一区二区 视频在线| 人人妻人人澡人人爽人人夜夜| av不卡在线播放| 一区二区三区精品91| 日韩制服丝袜自拍偷拍| 国产精品三级大全| 在线天堂中文资源库| 日韩免费高清中文字幕av| 亚洲精品视频女| 亚洲精品国产av蜜桃| av又黄又爽大尺度在线免费看| 亚洲欧美成人精品一区二区| 国产免费一区二区三区四区乱码| 母亲3免费完整高清在线观看 | 欧美日韩成人在线一区二区| 久久精品国产a三级三级三级| 亚洲国产欧美在线一区| 亚洲精品一二三| av卡一久久| 母亲3免费完整高清在线观看 | 国产极品天堂在线| 国产深夜福利视频在线观看| 夫妻午夜视频| 美女高潮到喷水免费观看| 久久久久网色| 欧美人与性动交α欧美精品济南到 | 黄色配什么色好看| 热99国产精品久久久久久7| 国产精品成人在线| 国产亚洲av片在线观看秒播厂| 搡老乐熟女国产| 国产欧美亚洲国产| 国产1区2区3区精品| 午夜福利在线免费观看网站| 久久午夜福利片| 久久久国产欧美日韩av| 精品酒店卫生间| 国产精品女同一区二区软件| 秋霞伦理黄片| √禁漫天堂资源中文www| 热99久久久久精品小说推荐| 国产日韩一区二区三区精品不卡| 宅男免费午夜| 波多野结衣一区麻豆| 丰满迷人的少妇在线观看| 秋霞伦理黄片| av一本久久久久| 久久久久精品人妻al黑| 国产1区2区3区精品| 国产白丝娇喘喷水9色精品| 亚洲国产欧美网| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区久久| 久久人妻熟女aⅴ| 狂野欧美激情性bbbbbb| 天天躁夜夜躁狠狠久久av| 青春草国产在线视频| 国产黄频视频在线观看| 亚洲av综合色区一区| 成人国产麻豆网| 午夜影院在线不卡| 欧美另类一区| 深夜精品福利| 国产免费一区二区三区四区乱码| 亚洲第一区二区三区不卡| 少妇 在线观看| 久久99精品国语久久久| 美女视频免费永久观看网站| 国产极品粉嫩免费观看在线| 亚洲国产精品国产精品| 中文精品一卡2卡3卡4更新| h视频一区二区三区| 亚洲av福利一区| 国产午夜精品一二区理论片| 人人澡人人妻人| 女人久久www免费人成看片| av国产久精品久网站免费入址| 黄频高清免费视频| 下体分泌物呈黄色| 日韩av不卡免费在线播放| 久久久久久免费高清国产稀缺| 99香蕉大伊视频| 高清在线视频一区二区三区| 可以免费在线观看a视频的电影网站 | 国精品久久久久久国模美| 最近中文字幕2019免费版| 下体分泌物呈黄色| av在线观看视频网站免费| 妹子高潮喷水视频| 26uuu在线亚洲综合色| 免费黄频网站在线观看国产| 亚洲国产欧美在线一区| 欧美激情 高清一区二区三区| 91在线精品国自产拍蜜月| 黑丝袜美女国产一区| 国产精品久久久久久精品电影小说| 视频在线观看一区二区三区| 亚洲av综合色区一区| 亚洲美女视频黄频| 97在线人人人人妻| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美成人精品一区二区| 亚洲精品久久午夜乱码| 欧美最新免费一区二区三区| 国产片特级美女逼逼视频| 国产精品 国内视频| 超碰97精品在线观看| 亚洲国产av影院在线观看| 天天影视国产精品| 亚洲欧洲精品一区二区精品久久久 | av国产精品久久久久影院| 国产在线免费精品| 美国免费a级毛片| 中文天堂在线官网| 一区二区三区激情视频| 婷婷色综合大香蕉| 国产精品久久久av美女十八| 国产人伦9x9x在线观看 | 久久影院123| 深夜精品福利| 91精品伊人久久大香线蕉| 美女脱内裤让男人舔精品视频| 人妻少妇偷人精品九色| 最新的欧美精品一区二区| 国产无遮挡羞羞视频在线观看| 久久这里只有精品19| 欧美97在线视频| 亚洲av电影在线观看一区二区三区| av有码第一页| 久久久久网色|