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

    微多普勒理論在海面目標檢測中的應用及展望

    2013-07-25 06:26:10陳小龍
    雷達學報 2013年1期
    關鍵詞:微動時頻雜波

    陳小龍 關 鍵* 何 友

    ?

    微多普勒理論在海面目標檢測中的應用及展望

    陳小龍 關 鍵何 友

    (海軍航空工程學院電子信息工程系 煙臺 264001)

    微動特性是目標物理特性之一,微動目標的雷達特征包含了對目標形狀、結構和運動的精細刻畫,同時微多普勒反映了信號的非平穩(wěn)特性,因此,在高海況條件下采用微多普勒理論分析海雜波及檢測海面目標具有很大的優(yōu)越性。該文首先從微多普勒機理和特點出發(fā),對動態(tài)海面散射雜波建模和海雜波多普勒特性分析方法等相關研究進行歸納與分析,表明微多普勒理論應用的必要性;然后,從海面微動目標回波建模和微動特征分析與檢測方法等方面重點介紹了微多普勒理論在海面目標檢測領域的應用和主要技術途徑;最后,針對現(xiàn)有研究中存在的問題,闡述了有待于進一步研究的方向。

    微動;微多普勒;海雜波;目標檢測

    1 引言

    海雜波中微弱目標尤其是“低(低掠射角)、慢(靜止或慢速運動)、小(目標尺寸小)”目標的檢測技術始終是雷達信號處理領域的難題,不僅具有重要的理論意義,而且在軍用和民用領域均占有非常重要的地位,如海面目標的檢測、船舶的安全航行、浮冰規(guī)避和海洋環(huán)境的監(jiān)測等。傳統(tǒng)的基于統(tǒng)計理論的目標檢測方法將海雜波視為隨機過程,但目標模型和雜波模型均呈多樣化發(fā)展趨勢,尤其是雜波分布模型,在復雜環(huán)境中往往不成立或不完全成立,這就使經(jīng)典目標檢測方法由于模型失配而不能取得預期的檢測結果。Haykin等人從混沌和分形角度出發(fā),采用關聯(lián)維數(shù)和盒維數(shù)等特征量從海雜波中檢測微弱目標,但當信雜比持續(xù)降低時,檢測結果卻不盡人意,而且也很難實現(xiàn)對運動目標的檢測。

    研究表明,強海雜波背景下,雜波幅度概率密度函數(shù)表現(xiàn)出較長的“拖尾”現(xiàn)象,并且此時的多普勒譜展寬,體現(xiàn)出明顯的時變和非平穩(wěn)特性,它的統(tǒng)計特性和多普勒頻率隨時間變化,因此海雜波的頻譜是時間和頻率兩個變量的函數(shù),應從時間和頻率兩方面對海雜波數(shù)據(jù)進行分析和處理。同時,海上目標隨海面顛簸導致姿態(tài)變化,引起雷達回波功率調制效應,不僅存在平動,艦船還繞參考點作三軸轉動,導致散射點的多普勒頻率隨時間非線性變化。因此,新的檢測算法應能較好地處理時變、非平穩(wěn)和非線性信號,同時反映信號的精細特征。

    在這種情況下,人們開始從其它方向尋找解決該問題的有效途徑。近年來,微多普勒理論的研究已成為信號處理領域的一個新的技術熱點。微多普勒反映了多普勒的變化特性,目標的幾何結構和精細運動對雷達回波的綜合調制特征則體現(xiàn)在目標的微動特征上,從而為雷達目標檢測和特征提取提供了新的途徑。目標距離的周期性將引起目標多普勒頻率的周期性,多普勒譜出現(xiàn)展寬,目標姿態(tài)變化將對回波產(chǎn)生調制特性,反映出目標的瞬時速度變化特性,而海面目標,如艦船目標,其回波多普勒譜實質是平動和微動共同作用的結果。因此,微多普勒非常適于分析海雜波的非平穩(wěn)特性以及海面目標回波信號。目前,雖然微多普勒的理論研究及應用處在起步階段,仍有很多問題亟待解決,但已成為國內外的一大熱點研究內容。本文從微多普勒效應的概念出發(fā),首先介紹了海雜波多普勒建模及特性分析方法的國內外研究現(xiàn)狀,在此基礎上結合海面目標電磁散射特性和運動模型,對目標微動特征提取與檢測方法進行總結和闡述,尤其對目前研究過程中出現(xiàn)的新方法進行了詳細介紹,最后,針對現(xiàn)有研究中存在的問題,展望了微多普勒理論在對海雷達回波信號分析中的發(fā)展趨勢,指出了進一步的研究方向。

    2 微多普勒研究概述

    微動最早出現(xiàn)在相干激光雷達中,美國海軍研究實驗室(Naval Research Laboratory)的Victor C. Chen最早提出了微動(Micro-motion)和微多普勒(Micro-Doppler)的概念并將其引入到雷達觀測中,將目標或目標部件除質心平動以外的振動、轉動和加速運動等微小運動統(tǒng)稱為微動,而由目標微動所引起的多普勒頻移稱為微多普勒。國防科技大學的陳行勇等人又在此定義的基礎上,將微動推廣為“目標或目標組成部分在徑向相對雷達的小幅非勻速運動或運動分量(相對于目標與雷達的徑向距離)”。從多普勒效應來看,若點目標在雷達徑向存在非勻速運動,則其多普勒隨時間變化,表現(xiàn)為頻譜展寬;同時,目標往往存在復雜的非剛體運動形式,各質元的相對運動也會使其回波頻譜出現(xiàn)旁瓣或展寬,稱為微多普勒效應(Micro-Doppler effect)。微多普勒反映了頻率的瞬時特性,描述了微動引起的瞬時多普勒變化特性,因此微多普勒具有時變性,即頻率隨時間變化。自然界中普遍存在著微動現(xiàn)象, 如行人手臂和四肢的擺動,人體的心跳和呼吸時胸腔的運動,電動機,履帶車、裝甲車和汽車車輪,直升機旋翼,彈道導彈彈頭的顫動(進動和章動),船舶的顛簸和擺動以及天線的轉動等。微動目標的雷達回波特征反映了目標的精細特征,如結構特性、電磁散射特性和運動特征等,因此,對于空間、空中、地面目標、海面艦船的探測和識別具有重要的意義。

    目前,國外以美國海軍研究實驗室、加拿大防御技術研究中心、美國Texas大學、新加坡南洋理工大學為主的研究機構和院校對雷達目標微多普勒開展了較為系統(tǒng)的研究,逐漸形成了從目標調制模型、特征分析與提取到方法驗證、實驗測量的研究體系。在工程應用上,微多普勒測量雷達系統(tǒng)由單一的多普勒測量功能發(fā)展到同時具有高多普勒分辨率和高距離分辨率的測量能力,美國在此方面的研究水平一直處于國際領先地位。2000年,Philip Gatt等人成功研制雙脈沖連續(xù)波激光雷達,成功探測地面和空中目標的微多普勒信號,將其用于特征提取和目標識別以及戰(zhàn)場態(tài)勢評估。2001年,美國Georgia技術研究院采用連續(xù)波雷達觀測車輛目標的微動特征并進行特征提取,已應用于汽車安全預警系統(tǒng)的研發(fā)中;2002年,美國研制的THAAD GBR X波段雷達能夠精確測量彈道導彈彈頭的微動特征,并對真假彈頭進行識別;2007年,美國Texas大學研制了X波段3個發(fā)射頻率的多頻連續(xù)波體制的微動測量雷達系統(tǒng),能夠在測量微多普勒的同時測距,采集并分析了行人、車輛和動物在不同機動情況下,多角度的微多普勒特征,其結果可用于目標自動檢測和識別。鑒于有限的微波波段資源和昂貴的實驗設備,2011年,英國倫敦大學的Alessio Balleri等人采用頻率捷變超聲雷達采集微多普勒特征,降低了試驗成本,但受聲信號在大氣中極易衰減的限制,僅能用于近程探測,如室內和機場旅客的監(jiān)視等。

    在國內,國防科技大學、空軍工程大學、中國航天二院二十三所和西安電子科技大學等多家科研機構均在微多普勒方面開展了一定的研究,在此方面有了一定的研究基礎,但有關微多普勒理論和技術的研究尚未成熟,研究內容主要集中在微多普勒效應的應用研究,包括微多普勒特征建模、微動特征提取、微動目標雷達SAR和ISAR成像等方面;研究對象包括高速運動目標,如地面車輛如坦克、裝甲車、導彈發(fā)射架等振動目標、海面角反射器、直升飛機為代表的旋翼目標、彈道導彈為代表的進動和章動目標以及行人等非剛體微動目標等。由此可見微多普勒已成為國內外的一大熱點研究內容。

    3 海雜波多普勒研究現(xiàn)狀

    由于海雜波對來自海面或接近海面的目標(包括低空掠海飛行的飛機、小型軍艦、航海浮標以及漂浮在海上的冰塊)的雷達回波的可檢測性形成嚴重制約,因此對海雜波的研究不僅具有理論上的重要性,而且具有實踐上的重要性。海雜波譜是海表面單個距離門內連續(xù)相干時間序列信號自相關函數(shù)的傅里葉變換,通常稱為功率譜。由于海表面是運動的,海雜波譜將產(chǎn)生多普勒頻移,所以也稱為多普勒譜。海雜波譜是雷達相參處理、多普勒域雜波抑制和動目標檢測的基礎。下面分別從動態(tài)海面散射雜波建模、多普勒譜特性以及在海雜波中微弱目標檢測中的應用等方面作回顧和總結。

    3.1 動態(tài)海面散射雜波建模

    由于海雜波譜與海表面的運動和海面的擾動狀態(tài)密切相關,海面復雜的運動引起散射的電磁波產(chǎn)生多普勒頻移、展寬、強度變化等特性,因此海雜波多普勒譜研究與海面的散射機理密不可分。近年來,人們對動態(tài)粗糙海面散射模型以及多普勒理論進行了廣泛而深入的研究,在諸如雷達目標成像、海洋環(huán)境監(jiān)測、電波傳播、遙感、海面目標檢測與跟蹤等領域中都有極其重要的意義。目前隨機粗糙面散射理論大致上可分為兩大類:第1類是數(shù)值方法,該方法雖然計算精度高,但計算復雜、耗時;第2類是近似的可數(shù)值求解的方法,較為經(jīng)典的有適用于大尺度隨機起伏的Kirchhoff近似法、小尺度起伏的微擾法以及獨立疊加這兩種起伏的雙尺度方法。第2類方法可以定量地解釋粗糙面產(chǎn)生散射的物理現(xiàn)象,但是其適用范圍、理論基礎以及實際應用中的一些問題,如低掠入射情況和高海情等情況,仍然有待于進一步研究。另外,文獻[24]將分形幾何用于自然粗糙面的模擬,利用隨時間變化的帶限分形模型來模擬海面,研究了具有自仿射分形特征的1維、2維隨機粗糙面的散射,但分形模型未考慮實際粗糙面的譜分布,需要選擇合適的尺度描述海面起伏的劇烈程度,也不能反映海雜波多普勒的變化情況。

    海浪由大尺度重力波和小尺度張力波組成,可將海面簡化為僅含有兩種尺度粗糙度的表面,即大尺度粗糙面和小尺度粗糙面,根據(jù)這一特性,提出了粗糙面電磁散射的雙尺度模型。然而,雙尺度法依賴粗糙面的劃分方式,且該模型基于海面散射回波信號是非時變的,即頻率不隨時間發(fā)生變化。 近年來,更多的研究表明,在高分辨率雷達對海觀測中,當以低掠射角照射粗糙海面或高海況時,海面回波強度會明顯增強,相應的雜波幅度概率密度函數(shù)表現(xiàn)出較長的“拖尾”現(xiàn)象,海面出現(xiàn)白浪等破碎波,雷達回波表現(xiàn)為海尖峰;文獻[26]指出此時的多普勒譜包括由“慢變信號”引起的Bragg散射,以及“快變信號”產(chǎn)生的非Bragg散射,非Bragg散射導致Bragg譜展寬,海面散射回波信號由非時變向時變轉變。海尖峰是海雜波信號非平穩(wěn)特性的一個重要體現(xiàn),它表現(xiàn)為隨機分布在不同距離、不同角度上的零星的運動或者靜止目標,但與一般目標不同,它們具有較為強烈的回波起伏特性。由于海尖峰出現(xiàn)時間較短,使得海雜波由穩(wěn)態(tài)向非穩(wěn)態(tài),非時變向時變轉換,體現(xiàn)出微多普勒特性,雷達有可能將海尖峰判斷為一個具有一定速度的機動目標,進而導致虛警概率的增加。因此,采用微多普勒理論對它進行深入研究,分析其對海雜波的影響是非常必要的?,F(xiàn)有文獻多是對海尖峰的定性描述和研究,對其進行的系統(tǒng)分析還很欠缺。M. Greco等人對海尖峰的統(tǒng)計特性進行了分析,認為幅度門限、最小間隔時間和最小尖峰寬度是描述海尖峰的3個重要因素,給出了海尖峰的判定方法,但是3個因素的確定僅通過統(tǒng)計分析得到,具有不確定性。圖1(a)給出了IPIX雷達19931107_135603數(shù)據(jù)在HH極化方式下的海尖峰判定情況,圖中紅色實線表示判別為海尖峰的海雜波數(shù)據(jù),其它的藍色虛線背景表示非海尖峰的海雜波數(shù)據(jù)。為了方便分析海尖峰特性,將海尖峰和背景海雜波分離,分析兩者的時頻分布,如圖1(b)和圖1(c)所示,可以明顯看出海尖峰幅度高于背景海雜波幅度,且頻率分布不同,具有時變特性。

    3.2海雜波多普勒特性分析

    海雜波的多普勒譜特征反映了海面本身的動態(tài)特性,一方面依賴于雷達工作參數(shù),另一方面和引起散射現(xiàn)象的海洋環(huán)境參數(shù)密切相關,主要包括極化、頻率、雷達相對風向的觀測角度、擦地角、海況、波浪速度及波浪類型等。眾多學者從雷達實測數(shù)據(jù)分析和隨機粗糙面散射理論兩方面對海雜波多普勒譜特性進行了大量的研究。然而,前者僅通過大量的數(shù)據(jù)統(tǒng)計得出經(jīng)驗模型,未從理論上對多普勒譜頻移及展寬機理進行分析;后者針對的是特定背景或環(huán)境,粗略的假設與實際海面不符,不能很好地反映真實海面的時變性和非平穩(wěn)特征,而且所提出的越來越復雜的建模方法帶來的是實時性的急劇降低或者缺乏可實現(xiàn)性。對于“慢變信號”引起的Bragg散射,通過傳統(tǒng)的Fourier變換(Fourier Transform, FT)得到海面散射回波的時頻分析,能夠給出特定時間和特定頻率范圍的能量分布;對于“快變信號”產(chǎn)生的非Bragg散射,回波信號頻率隨時間發(fā)生變化,若仍采用Fourier變換不能很好地反映和提取頻率的變化信息。

    圖1 IPIX雷達數(shù)據(jù)海尖峰判別及時頻分析(19931107_135603)

    時頻分析方法是研究時變、非平穩(wěn)信號的有力工具,作為時間和頻率的2維函數(shù),時頻分布給出了特定時間和特定頻率范圍的能量分布,也描述了非平穩(wěn)信號的頻率隨時間變化的過程,適于分析海雜波的多普勒特性并應用到弱目標檢測中。主要方法包括基于FT的譜相減法、基于短時Fourier變換(Short Time Fourier Transform, STFT)和Wigner- Vill變換(Wigner-Vill Distribution, WVD)的峰值估計法、平滑偽WVD法(Smoothed Pseudo WVD Distribution, SPWVD)、小波變換(Wavelet Transform, WT)和自適應chirplet分解法等,各種方法的基本原理和性能比較如表1所示。然而,由于海雜波信號往往表現(xiàn)為非高斯、非平穩(wěn)、時變、多分量信號,經(jīng)典的時頻變換方法和濾波方法有很大的局限性。尤其是高海況時,海面起伏變得劇烈、粗糙,海雜波的幅度和多普勒均隨時間變化,“快變信號”產(chǎn)生非Bragg譜使得多普勒的中心頻率偏移或展寬,降低了雷達目標的檢測性能。近年來,人們開始研究新的變換方法研究海雜波特性,并用于弱目標檢測。

    表1基于經(jīng)典時頻分析的海雜波的多普勒分析方法比較

    Tab. 1 Comparison of Doppler analysis methods of sea clutter based on classical time-frequency representation

    分數(shù)階Fourier變換(FRactional Fourier Transform, FRFT)作為一種統(tǒng)一的時頻變換,將信號分解在一組正交的Chirp基上,更適于分析或處理時變的非平穩(wěn)信號,能夠反映多普勒的變化規(guī)律。與常用二次型時頻分布不同的是FRFT采用單一變量表示時頻信息,降低了運算量,又是一種線性變換,沒有交叉項干擾。文獻[4,23]采用FRFT積累動目標回波信號能量,用于解決海雜波中背景下的目標檢測問題,結果表明,在FRFT域能夠很好地反映海雜波的時變特性。圖2給出了IPIX雷達19931107_135603純海雜波數(shù)據(jù)的FRFT譜分布,可見,海雜波在FRFT域幅值起伏變化劇烈,頻率較高,在變換階數(shù)=1(頻域)周圍能量分布相對集中,但海雜波也具有微弱變化的加速度,能夠很好地體現(xiàn)和描述較高海況情況下海雜波中的非Bragg散射回波特性。FRFT算法從不同程度上對非平穩(wěn)信號的時變性給予了恰當?shù)拿枋?,彌補了Fourier分解的不足,但仍屬于全局分析的范疇,究其原因在于基函數(shù)決定了對信號的分析能力。1998年,Huang等人提出了一種基于經(jīng)驗的模式分解算法(Empirical Mode Decomposition, EMD)和基于Hilbert變換的時頻譜圖,即Hilbert-Huang變換(Hilbert-Huang Transformation, HHT),具有優(yōu)秀的時頻分辨性能。EMD是基于數(shù)據(jù)時域局部特征的,它能將復雜信號分解成若干個具有不同頻率的固有模態(tài)函數(shù)(Intrinsic Mode Function, IMF),得到由信號瞬時頻率和瞬時幅值構成的完整時頻譜,從而使瞬時頻率具有實際的物理意義。由于分解是基于信號時域局部特征的,因此分解是自適應的,特別適合分析非平穩(wěn)、非線性隨機過程。文獻[36,37]采用HHT分析雷達目標和海雜波信號,進行了初步探討。采用EMD方法,對海雜波幅值信號進行篩選,可獲得其IMF, IPIX雷達的19931118_023604和19931118_162658純海雜波幅值經(jīng)過EMD分解后得到的IMF如圖3所示,()表示原始海雜波的幅值信號,(),,()表示海雜波的各個IMF, IMF成份的頻率從()到()依次降低,故利用EMD可在時域內將信號按頻率高低進行分解,由于IMF每一時刻僅包含單一頻率成份,因此,通過Hilbert變換可以得出各IMF分量的時頻圖,定量地描述時間與頻率的關系。而STFT和WT等方法只能定性地描述時間和頻率的關系,而不能定量、精確地描述時間與頻率的關系。

    4 海面微動目標回波信號建模及檢測研究現(xiàn)狀

    海面目標特性分析與建模是目標檢測的一項重要內容,不過與海雜波特性分析相比,基于雷達實測回波數(shù)據(jù)的目標特性研究相對較少,原因在于,實際雷達觀測條件下難以獲得只有目標信號的回波數(shù)據(jù)。目前應用較多的目標模型是起伏模型,包括Swerling起伏模型 (也稱為Rayleigh起伏)、Rice起伏模型、起伏模型、對數(shù)正態(tài)起伏模型以及Nakagami起伏模型等。然而,起伏模型僅從統(tǒng)計角度對目標的幅度起伏特征進行建模,不能反映海面目標隨海面顛簸導致的姿態(tài)變化,引起的雷達回波功率的調制效應。因此,眾多國內外學者從實際出發(fā),對運動目標及其運動特性進行建模:一是海面運動目標的雷達回波幅度建模,包括運動目標電磁散射特性、回波的幅度調制、本身形狀、大小等靜態(tài)幾何參量的確定;二是對海面目標運動特征的建模,即描述目標在海面上的運動方程。

    圖2 IPIX雷達海雜波數(shù)據(jù)的FRFT譜(19931107_135603)

    圖3 IPIX雷達海雜波數(shù)據(jù)幅值的IMF

    4.1海面微動目標回波信號建模

    4.1.1海面目標電磁散射特性建模 對海面目標的電磁散射特性的建模是一項非常復雜的工作。首先,從艦船目標本身考慮,由于其具有很大的電尺寸及復雜的結構,電磁散射機理十分復雜;其次,海表面是時變的,由于其波浪起伏使得海面成為巨大的粗糙表面,時變海面對電磁波的反射也隨時間變化;再次,海面艦船目標同海面之間的電磁耦合十分復雜,加之海面與艦船之間的相互作用使得艦船在海面的姿態(tài)不斷變化,因而海面艦船目標的電磁散射計算問題變得更為復雜。Victor C. Chen研究了微動目標的雷達后向散射特性,指出目標微動對電磁波具有回波功率調制特征,可根據(jù)微動目標的RCS時間序列反演微動規(guī)律,提取目標微動信息,從而為研究微動目標的回波功率調制特性奠定了理論基礎。2010年,ángel F等人首次基于ISAR圖像和微動特征分析剛體和非剛體運動目標的雷達后向散射特性,并通過仿真的艦船、直升機和人體驗證算法,能夠提取目標的微多普勒信息。由此可見,通過計算微動目標的電磁散射特性,研究目標復雜運動對雷達回波的RCS調制作用,能夠進一步揭示目標姿態(tài)變化引起的散射強度變化特性,從而更好地分析和提取微動特征。

    4.1.2海面目標運動建模 在海面目標運動建模方面,國內外許多學者從多種角度進行了研究,運動模型也經(jīng)歷了由簡單到復雜,由單散射目標到多散射目標,由剛體到非剛體,由平動到轉動,由周期運動到非周期運動,由時域到頻域再到變換域等幾個主要研究階段。傳統(tǒng)的對海面運動目標的建模僅認為目標在3維空間中平動,即艦船在運動過程中,坐標系的各坐標軸永遠相互平行,表現(xiàn)為勻速運動、勻加速運動和變加速運動等,根據(jù)Weierstrass近似原理,其回波信號可由足夠階次的多項式相位信號近似表示,而運動狀態(tài)比較復雜的目標在一段短的時間范圍內,??捎镁€性調頻(Linear Frequency Modulation, LFM)信號作為其一階近似?;诖四P?,文獻[45]研究了兩種海雜波背景下的微弱勻加速運動目標的檢測方法,在最大程度保留目標信號能量的同時盡量地抑制海雜波,達到低信雜比下有效地檢測出海面運動目標的效果,但算法的通用性不強。由于艦船隨海面波動而產(chǎn)生顛簸,一方面,在高海況條件下,海雜波的存在降低了回波信號的信雜比;另一方面,艦船姿態(tài)變化復雜,不僅存在平動,艦船還繞參考點作3軸轉動以及天線的轉動,導致散射點的多普勒頻率隨時間非線性變化。引起的姿態(tài)變化包括艦船以其前進方向為軸作滾動,以在水平面內與前進方向垂直的方向為軸作俯仰運動,在由船身確定的平面內以塔臺的方向為軸作偏航運動等3種運動形式,如圖4所示。根據(jù)微多普勒的定義,海面艦船目標的微動表現(xiàn)為平動中的非勻速運動、3軸轉動以及天線的轉動。由此可見,傳統(tǒng)的海面目標平動模型在復雜的海洋環(huán)境和日益提高的目標檢測要求下,越來越顯得捉襟見肘,主要體現(xiàn)為目標模型難以描述轉動運動形式以及對目標運動狀態(tài)的時變性、周期性考慮不足等。

    目前,對艦船目標的3軸轉動的研究主要集中于ISAR艦船成像方面,對于岸基高分辨雷達觀測海面艦船卻少有報道。ISAR對艦船成像依靠艦船與雷達視線之間的相對轉動,它由艦船航行和海浪作用下艦船自身搖擺產(chǎn)生的轉動組成。只有首先了解艦船運動的特點,才能更有效地設計出合適的ISAR成像技術和算法。邢孟道等人對外場艦船目標進行了ISAR成像試驗,分析表明由于海水波動等原因,使艦船作3維運動,散射點子回波為調幅-調頻信號(Amplitude Modulation-Frequency Modulation, AM-FM),在一定程度上,可近似為調幅-線性調頻(AM-LFM)信號,或者近似分段的AM-LFM信號,則對艦船目標的瞬時成像問題轉化為噪聲和雜波背景下的多分量AM-LFM信號瞬時參數(shù)估計問題。針對多分量的AM-FM信號給信號的參數(shù)估計帶來的極大困難,人們相繼提出了基于自適應Chirplet分解、匹配傅里葉變換(Matching Fourier Transform, MFT和Chirp傅里葉變換(Chirp Fourier Transform, CFT)的艦船目標ISAR成像方法。然而,已有的文獻并未詳細推導艦船的3維轉動與AM-FM信號參數(shù)的對應關系,也沒有考慮復雜運動情況下的微動特征,難以對時變的多普勒分量進行分析和提取。

    圖4 雷達觀測海面目標幾何關系圖

    4.2目標微動特征分析與檢測方法

    微動本質上是一種非勻速運動或周期運動,微動目標雷達回波受到調制,其多普勒隨時間非線性變化,稱為一類頻率調制的時變信號。微多普勒信號具有非線性和非平穩(wěn)的特點,同時微動目標回波一般為多分量信號,因此,微動目標雷達特征分析和提取以非平穩(wěn)、時變信號、多分量信號處理技術為主要工具。

    4.2.1基于Fourier變換的頻譜分析方法 基于FT的頻譜分析方法是最早用于微多普勒信號分析的工具,人們相繼提出了利用多普勒譜寬和多普勒譜峰值比的振動振幅和相位估計方法、包絡檢測法、多諧波微多普勒信號分析方法以及循環(huán)平穩(wěn)特征調制相位信號參數(shù)估計方法等,并應用于醫(yī)學成像、高速公路上車輛檢測、進動彈頭回波分析以及噪聲抑制等方面。然而,頻譜分析方法為一種全域變換,缺乏時間局域性,不能有效處理時變的非平穩(wěn)微多普勒信號。

    4.2.2經(jīng)典時頻分析方法 時頻分析方法被廣泛用于微多普勒的特性分析以及微動特征提取與參數(shù)估計。第一類為線性時頻表示,如STFT, WT和Gabor變換等,線性時頻分布無交叉項,但受不確定原理的制約,時頻分辨率較低;另一類為非線性時頻表示,主要包括Cohen類時頻分布,如WVD等。目前,基于Cohen類方法分析微多普勒信號存在以下不足:(1)非線性頻率調制信號或多元信號的時頻分布中存在交叉項影響;(2)典型的微動目標雷達回波模型為多分量正弦調頻信號(Sinusoidal Frequency-Modulated, SFM),其相位是無限階可導的,因此采用WVD處理SFM信號不能得到理想的時頻聚集性;(3)通常微動信號是微弱的調頻信號,而Cohen類方法對信雜比要求較高,從而影響參數(shù)估計的準確性。

    微多普勒信號的特點要求信號分析方法應具有低交叉項、高時頻分辨率和大動態(tài)范圍的特點,以便更好地揭示目標微多普勒特征。從本質上說,微多普勒是由于目標運動引起的多普勒頻移,因此,微多普勒信號處理也即對調頻信號的處理。在這種情況下,國內外學者開始從其它方向尋找解決該問題的有效途徑。

    4.2.3相位匹配分析方法 相位匹配處理方法是根據(jù)微多普勒信號的形式,設計與其相對應的指數(shù)基函數(shù),并將微多普勒信號在此基函數(shù)上進行分解。根據(jù)不同基函數(shù)形式可設計出不同的基函數(shù)分解方法,例如MFT、自適應Chirplet分解法(線性調頻小波基函數(shù)分解法)和正弦調頻基函數(shù)分解法等。MFT使用的變換因子是一組正交函數(shù)集,故不存在多分量信號交叉項的影響。Chirplet分解與時頻變換、FT等方法相比,能夠方便地實現(xiàn)信號在二階相位上的分解與重構。正弦調頻基函數(shù)分解法采用調頻基函數(shù)提取振動或旋轉頻率以及微多普勒成分的相位信息,不需要信號成分的個數(shù),但需已知微多普勒信號的最大瞬時頻移。Victor C. Chen指出獨立成分分析(Independent Component Analysis, ICA)技術更適合分析具有空間局域性的微多普勒譜圖特征,他提出了時空聯(lián)合獨立成分分析(Spatial-temporal ICA),將2維空間ICA推廣至3維時間-空間域,可將微多普勒譜分解為對應信號不同微動特征的基函數(shù)的線性組合。由于微動信號在一定程度上可近似為AM-LFM信號,而FRFT對LFM有良好的能量聚集性和檢測性能,同時根據(jù)分段信號長度自適應地選擇FRFT的長度,可近似得到階段的微動信號,達到分析和提取微多普勒信號的目的。目前,少有文獻將FRFT作為微多普勒信號分析的工具,主要受限于FRFT缺少時域定位的功能。為此,陶然等人在FRFT中加入滑動的短時窗函數(shù),得到短時FRFT(Short-Time FRFT, STFRFT),通過窗函數(shù)的滑動完成整個時間上的信號局部性質分析,可得到在任意時刻的該段信號的頻率變化,極大地擴展了FRFT的應用范圍。隨著國內外學者對FRFT理論的研究不斷深入,這種新的時頻分析工具必定成為分析微多普勒信號的理想手段。

    4.2.4參數(shù)模型分析方法 信號時頻處理方法可分為非參數(shù)模型法和參數(shù)模型法兩大類。上述方法均為非參數(shù)模型法,其作為微多普勒特征分析工具具有不可比擬的優(yōu)勢,但估計性能受時頻分辨率的限制。而對目標微動特征進行參數(shù)化估計,在建立目標微動參數(shù)估計模型的基礎上,僅需要求解最優(yōu)模型參數(shù),便能較好地分析和提取微動特征。目前,基于參數(shù)模型的微多普勒信號分析的研究尚處于起步階段,研究的方法主要有時變參數(shù)模型法和基于稀疏表示的參數(shù)估計模型法。時變參數(shù)模型法是近年來應用于非平穩(wěn)隨機信號分析與處理的一種新方法,其優(yōu)點在于將非平穩(wěn)時變問題轉化為平穩(wěn)的線性時不變問題,進一步提高了參數(shù)估計的精確度,其中具有代表性的是時變自回歸(Time- Varying AutoRegression, TVAR)模型法。研究表明,單分量AM-FM信號可由一階TVAR模型表示,多分量AM-FM信號可由高階TVAR建模,因此,TVAR模型適用于對微動目標回波進行參數(shù)化建模和估計。基于時頻分析的參數(shù)估計方法可以看作將信號在時頻基函數(shù)上的分解,如果信號的特性與分解的基函數(shù)相匹配,就可以采用某幾個基函數(shù)的組合來表示原始信號,即信號可以被稀疏表示。基于時頻分析的方法的基本思想是將信號在一組完備的時頻基上展開,而如果將能夠很好地刻畫信號局部時頻結構的時頻原子構成的過完備字典替代完備基函數(shù),使得信號的自適應表示稱為可能,則參數(shù)估計問題轉化為信號的稀疏表示問題,并且對頻率具有超分辨能力,從而更有利于獲得目標精細特征。目標雷達回波可視為少數(shù)強散射中心回波的疊加,回波具有稀疏特性,因此,采用稀疏表示的方法分析微動信號,并進行參數(shù)估計是非常適合的。

    5 研究展望

    由國內外研究現(xiàn)狀可知,微多普勒理論在精細雷達特征信號處理中具有良好的表現(xiàn),可以很好地刻畫信號本身的頻率變化。微動目標的RCS調制從功率上描述了微動目標雷達特征,表征了目標姿態(tài)的變化。海面目標的非勻速平動和3軸轉動導致散射點的多普勒頻率隨時間非線性變化,因此,海面目標的微動特征能夠反映目標的雷達回波調制特性、幾何結構和運動特征,同時展示了時間、空間、頻率3維特征空間的變化特性。這些均為微多普勒理論在對海雷達回波信號分析中的應用提供了很好的基礎與借鑒。從現(xiàn)有研究狀況分析,未來研究中需要進一步解決的問題包括以下幾個方面。

    5.1海雜波建模及特性分析方面

    隨著高分辨力雷達的發(fā)展,對海面微弱目標的檢測與識別成為可能。但在低掠射角及高海況的條件下,“快變信號”(海尖峰)導致海面起伏變劇烈、粗糙,其產(chǎn)生的非Bragg譜導致多普勒的中心頻率頻移或展寬,海面回波信號表現(xiàn)出時變特性。此時,海面回波信號已不能采用傳統(tǒng)的雙尺度模型建模。經(jīng)典時頻變換方法和濾波方法有很大的局限性,很少考慮各個因素之間的相互作用(包括雷達的工作狀態(tài)(入射角、發(fā)射頻率、極化、分辨率等),以及背景狀況(如海況,風速、風向等)),也沒有揭示出海雜波產(chǎn)生的物理機制。目前尚未有學者從微多普勒的角度對海雜波進行建模和特性分析,而海雜波的“三非”特性(非線性、非平穩(wěn)、非高斯),與微動信號特征存在相似性,但也有不同。采用微動信號分析方法對海雜波建模,能夠更好地揭示海雜波的起伏特性及頻率特征,充分掌握海雜波的特征信息,一方面,從海表面運動特性的角度出發(fā),研究白浪、破碎波的姿態(tài)變化對雷達回波信號的調制作用(RCS調制和相位調制);另一方面,根據(jù)海尖峰的微動特性,建立相應的微多普勒譜模型,并研究微多普勒譜的偏移和展寬特性,從而為海雜波抑制和微弱目標檢測提供必要的先驗信息。

    5.2海面微動目標回波建模及特性分析方面

    艦船等微動目標在海面上將產(chǎn)生由轉動引起的姿態(tài)變化包括偏航、滾動和俯仰3種周期運動,這3種轉動之間可能同時存在,而且,由于受非線性策動力和非線性阻尼力的作用,3維轉動都呈現(xiàn)出隨機性和多倍周期的特點,海面的波動起伏以及海面反射引起的多徑效應,使得艦船姿態(tài)變化較為復雜。此時,需要分別研究海面剛體目標和非剛體目標的運動特征,由簡單到復雜,由一維到多維,由僅考慮目標到同時考慮目標和海面起伏,建立海面微動目標回波模型。對微動特征的提取和識別是建立在微動的散射中心和微動物體本身一致的假設基礎上的,然而,隨著觀測角度的變化,微動散射中心也隨之變化,頻率變化不滿足正弦調頻的形式,以上這些因素都會對海面微動目標的建模造成困難,有待于進一步研究。同時,在微動目標特性分析方面,應分別從理論模型和實測數(shù)據(jù)兩方面分析海雜波和微動目標的微多普勒特性,尋找回波信號的微多普勒區(qū)域,從時域、極化域、變換域和稀疏域多個角度進行分析,擴展數(shù)據(jù)的利用維度,判定信號是否具有微動特性,目的在于充分掌握海雜波及海面目標的微動特征信息,為后續(xù)的微動目標檢測、特征提取與參數(shù)估計打下基礎。

    5.3微多普勒信號特征分析和檢測方面

    目前,國內外對雷達目標微動特征的分析和研究還主要集中在點目標微多普勒方面,在對微動引起的雷達信號調制效應分析方面,如電磁波、RCS、時域、變換域和極化域調制等,還未見系統(tǒng)的闡述。通過建立完善的雷達目標微動特征體系,將會為雷達目標的檢測和識別提供新的途徑。微多普勒是估計目標微動參數(shù)的前提,關鍵在于瞬時頻率的高精度估計和提取。微動特征的提取集中于微動參數(shù)估計,并未考慮微動目標姿態(tài)變化對回波產(chǎn)生的影響,聯(lián)合微動參數(shù)估計和散射強度變化引起的周期性調制效應等,可揭示更加豐富的目標微動特征信息。目前,國內外對微多普勒信號的分析方法主要應用了非平穩(wěn)信號處理方法如時頻分析等,而在時頻分辨率、多分量信號分析、對信噪比的要求以及魯棒性等方面,經(jīng)典的時頻分析方法有待改進。FRFT, HHT及其改進方法和基于參數(shù)模型的微多普勒信號分析方法(時變參數(shù)模型法和稀疏表示法)為微動目標的精細描述提供了新的途徑,但其理論體系仍還不夠完善,在參數(shù)估計精度及算法計算量方面還存在不足。另外,由于艦船目標的轉動速度比較慢,要觀測到精細的微多普勒信息需要較長的積累時間,而長時間里會產(chǎn)生較大的高階相位項。并且在長相參積累時間內也常會出現(xiàn)跨距離單元的問題,導致目標信號能量無法有效積累。因此,多目標、長時間積累、距離徙動條件下微動特征增強與提取方法仍需進一步研究。

    6 結束語

    本文著眼于微多普勒理論在非平穩(wěn)信號處理領域中的應用,回顧了微多普勒的研究發(fā)展歷程,總結了動態(tài)海面散射雜波建模方法和海雜波多普勒特性分析方法,重點介紹了微多普勒在海面微動目標檢測中的應用現(xiàn)狀、主要技術途徑、發(fā)展前景以及當前存在的問題。通過對其發(fā)展歷程和應用場景的回顧與總結,可以發(fā)現(xiàn)微多普勒理論已在海雜波建模和目標檢測領域中發(fā)揮著重要作用,是繼多普勒理論后又一重要的信號分析方法。微多普勒拓展了信號利用的維度,可提取更多有用的信息,因此還可以用于目標識別和成像等領域??梢灶A期,微多普勒技術的成熟和發(fā)展將彌補已有方法的不足,能夠進一步提升雷達目標探測和識別能力。

    [1] Ward K D and Watts S. Use of sea clutter models in radar design and development[J].,&, 2010, 4(2): 146-157.

    [2] Haykin S, Bakker R, and Currie B W. Uncovering nonlinear dynamics: the case study of sea clutter[J]., 2002, 90(5): 860-881.

    [3] Greco M, Stinco P, and Gini F. Identi?cation and analysis of sea radar clutter spikes[J].,&, 2010, 4(2): 239-250.

    [4] Guan Jian, Chen Xiao-long, and He You. Adaptive fractional Fourier transform-based detection algorithm for moving target in heavy sea clutter[J].,&, 2012, 6(5): 389-401.

    [5] Chen V C, Li Fa-yin, Ho Shen-shyang,.. Micro-Doppler effect in radar: phenomenon, model, and simulation study[J]., 2006, 42(1): 2-21.

    [6] Chen V C, Li F, Ho S S,.. Analysis of micro-Doppler signatures[J].-,, 2003, 150(4): 271-276.

    [7] Renhorn I, Karlsson C, and Letalick D. Coherent laser radar for vibrometry: robust design and adaptive signal processing[C]. Proceedings of SPIE, 1995, 2472: 23-30.

    [8] 陳行勇. 微動目標雷達特征提取技術研究[D]. [博士論文], 國防科學技術大學, 2006: 20-22.

    Chen Xing-yong. Research on radar signature extraction from target with micro-motions[D]. [Ph.D. dissertation], National University of Defense Technology, 2006: 20-22.

    [9] Thayaparan T, Abrol S, Riseborough E,.. Analysis of radar micro-Doppler signatures from experimental helicopter and human data[J].,&, 2007, 1(4): 289-299.

    [10] Cai Cheng-jie, Liu Wei-xian, Jeffrey Shiang Fu,.. Radar micro-Doppler signature analysis with HHT[J]., 2010, 46(2): 929-938.

    [11] Philip Gatt, Henderson S W, Thomson J A,.. Micro- Doppler lidar signals and noise mechanisms: theory and experiment[C]. Proceedings of SPIE, 2000, 4035: 422-435.

    [12] Gene Greneker, Jon Geisheimer, and David Asbell. Extraction of micro-Doppler data from vehicle targets at X-band frequencies[C]. Proceedings of SPIE, 2001, 4374: 1-9.

    [13] Lemnios W Z and Grometstein A A. Overview of the Lincoln Laboratory ballistic missile defense program[J]., 2002, 3(1): 9-32.

    [14] Anderson M G and Rogers R L. Micro-Doppler analysis of multiple frequency continuous wave radar signatures[C]. Proceedings of SPIE, 2007, 6547: 65470A.

    [15] Alessio Balleri, Karl Woodbridge, and Kevin Chetty. Frequency-agile non-coherent ultrasound radar for collection of micro-Doppler signatures[C]. IEEE Radar Conference, Washington D.C., USA, May 2011: 23-27.

    [16] 莊釗文, 劉永祥, 黎湘. 目標微動特性研究進展[J]. 電子學報, 2007, 35(3): 520-525.

    Zhuang Zhao-wen, Liu Yong-xiang, and Li Xiang. The achievements of target characteristic with micro-motion[J]., 2007, 35(3): 520-525.

    [17] 黃孟俊, 陳建軍, 趙宏鐘, 等. 海面角反射器干擾微多普勒建模與仿真[J]. 系統(tǒng)工程與電子技術, 2012, 34(9): 1781-1787.

    Huang Meng-jun, Chen Jian-jun, Zhao Hong-zhong,.. Micro-Doppler modeling and simulating of corner reflector in sea surface[J]., 2012, 34(9): 1781-1787.

    [18] Li Kai-ming, Liang Xian-jiao, Zhang Qun,.. Micro- Doppler signature extraction and ISAR imaging for target with micromotion dynamics[J]., 2011, 8(3): 411-415.

    [19] 高紅衛(wèi), 謝良貴, 文樹梁, 等. 基于微多普勒分析的彈道導彈目標進動特性研究[J]. 系統(tǒng)工程與電子技術, 2008, 30(1): 50-52.

    Gao Hong-wei, Xie Liang-gui, Wen Shu-liang,.. Research on precession of ballistic missile warhead based on micro- Doppler analysis[J]., 2008, 30(1): 50-52.

    [20] Bai Xue-ru, Zhou Feng, Xing Meng-dao,.. High resolution ISAR imaging of targets with rotating parts[J]., 2011, 47(4): 2530-2543.

    [21] Jeong-Won Park and Joong-Sun Won. An efficient method of Doppler parameter estimation in the time–frequency domain for a moving object from TerraSAR-X data[J]., 2011, 49(12): 4771-4787.

    [22] Zhang Min, Chen Hui, and Yin Hong-cheng. Facet-based investigation on EM scattering from electrically large sea surface with two-Scale pro?les: theoretical model[J]., 2011, 49(6): 1967-1975.

    [23] 陳小龍, 劉寧波, 宋杰, 等. 海雜波FRFT域分形特征判別及動目標檢測方法[J]. 電子與信息學報, 2011, 33(4): 823-830.

    Chen Xiao-long, Liu Ning-bo, Song Jie,.. Fractal feature discriminant of sea clutter in FRFT domain and moving target detection algorithm[J].&, 2011, 33(4): 823-830.

    [24] Marco Martorella, Fabrizio Berizzi, and Enzo Dalle Mese. On the fractal dimension of sea surface backscattered signal at low grazing angle[J]., 2004, 52(5): 1193-1204.

    [25] Walker D. Experimentally motivated model for low grazing angle radar Doppler spectra of the sea surface[J].-, 2000, 147(3): 114-120.

    [26] Chen Xiao-long, Wang Guo-qing, Huang Yong,.. Study on the fractional power spectrum of the scattered echoes from 1-D time-varying sea surface[C]. Asia-Pacific Conference on Antennas and Propagation, Singapore, August 2012: 321-322.

    [27] Guan Jian, Liu Ning-bo, Huang Yong,.. Fractal Poisson model for target detection within spiky sea clutter[J]., 2013, 10(2): 411-415.

    [28] Farina A, Gini F, Greco M V,.. High resolution sea clutter data statistical analysis of recorded live data[J].,, 1997, 144(3): 2808-2811.

    [29] Wang Yun-hua, Guo Li-xin, and Wu Zhen-sen. Modified two-scale model for electromagnetic scattering from the non-Gaussian oceanic surface[J]., 2005, 22(11): 2808-2811.

    [30] Allan N, Trizna D B, and McLaughlin D J. Numerical comparison of techniques for estimating Doppler velocity time series from coherent sea surface scattering measurements[J].-,, 1998, 145(6): 367-373.

    [31] Yasotharan A and Thayaparan T. Time-frequency method for detecting an accelerating target in sea Clutter[J]., 2006, 42(4): 1289-1310.

    [32] Hao Yan-ling, Tang Yan-hong, and Zhu Yi. A combined denoising algorithm approach to sea clutter in wave monitoring[C]. 2nd International Symposium on Intelligent Information Technology Application, Shanghai, China, December 2008: 99-103.

    [33] Zhang Y, Qian S, and Thayaparan T. Detection of a maneuvering air target in strong sea clutter via joint time-frequency representation[J].,&, 2008, 2(3): 216-222.

    [34] 陶然, 鄧兵, 王越. 分數(shù)階傅里葉變換及其應用[M]. 北京: 清華大學出版社, 2009, 第2章.

    Tao Ran, Deng Bing, and Wang Yue. Fractional Fourier Transform and Its Applications[M]. Beijing: Tsinghua

    University Press, 2009, Ch. 2.

    [35] Huang N E, Zheng S, Steven R,.. The empirical mode decomposition and the Hilbert spectrum for nonlinear non-stationary time series analysis[J].,&, 1998, 454(1971): 903-995.

    [36] 關鍵, 張建. 基于固有模態(tài)能量熵的微弱目標檢測算法[J]. 電子與信息學報, 2011, 33(10): 2494-2499.

    Guan Jian and Zhang Jian. Weak target detection based on intrinsic mode energy entropy[J].&, 2011, 33(10): 2494-2499.

    [37] 張建, 關鍵, 董云龍, 等. 基于局部Hilbert譜平均帶寬的微弱目標檢測算法[J]. 電子與信息學報, 2012, 34(1): 121-127.

    Zhang Jian, Guan Jian, Dong Yun-long,.. Weak target detection based on the average bandwidth of the partial Hilbert spectrum[J].&, 2012, 34(1): 121-127.

    [38] Derakhtian M, Khaliliazad Z, Masnadi Shirazi M A,.. Generalized likelihood ratio test-based multiple-target detection for fluctuating targets with unknown parameters[J].,&, 2011, 5(6): 613-625.

    [39] Chen V C. Doppler signatures of radar backscattering from objects with micro-motions[J]., 2008, 2(3): 291-300.

    [40] ángel F García-fernández, Omar A Yeste-Ojeda, and Jesús Grajal. Facet model of moving targets for ISAR imaging and radar back-scattering simulation[J]., 2010, 46(3): 1455-1467.

    [41] Guan Jian and Zhang Xiao-li. Subspace detection for range and Doppler distributed targets with Rao and Wald tests[J]., 2011, 91(1): 51-56.

    [42] Lei Peng, Sun Jin-ping, Wang Jun,.. Micromotion parameter estimation of free rigid targets based on radar micro-Doppler[J]., 2012, 50(10): 3776-3786.

    [43] Chen Xiao-long, Guan Jian, and He You. Sea clutter suppression and moving target detection method based on clutter map cancellation in FRFT domain[C]. CIE International Conference on Radar, Chengdu, China, October 2011: 438-441.

    [44] 郭海燕, 董云龍, 關鍵. 基于分數(shù)階模糊函數(shù)的海面運動弱目標檢測[J]. 系統(tǒng)工程與電子技術, 2011, 33(6): 1212-1216.

    Guo Hai-yan, Dong Yun-long, and Guan Jian. Weak target detection in sea clutter based on fractional ambiguity function[J]., 2011, 33(6): 1212-1216.

    [45] 關鍵, 李寶, 劉加能, 等. 兩種海雜波背景下的微弱勻加速運動目標檢測方法[J]. 電子與信息學報, 2009, 31(8): 1898-1902.

    Guan Jian, Li Bao, Liu Jia-neng,.. Two approaches of detecting weak moving target with constant acceleration in sea clutter[J].&, 2009, 31(8): 1898-1902.

    [46] Li Y, Wu R, Xing M D,.. Inverse synthetic aperture radar imaging of ship target with complex motion[J].,&, 2008, 2(6): 395-403.

    [47] 邢孟道, 保錚. 外場實測數(shù)據(jù)的艦船目標ISAR成像[J]. 電子與信息學報, 2001, 23(12): 1271-1277.

    Xing Meng-dao and Bao Zheng. ISAR ship imaging of real data[J].&, 2001, 23(12): 1271-1277.

    [48] 王勇, 姜義成. 基于自適應Chirplet分解的艦船目標ISAR成像[J]. 電子與信息學報, 2006, 28(6): 982-984.

    Wang Yong and Jiang Yi-cheng. The ISAR imaging of ship based on adaptive Chirplet decomposition[J].&, 2006, 28(6): 982-984.

    [49] 黃雅靜, 曹敏, 付耀文, 等. 基于匹配傅里葉變換的勻加速旋轉目標成像[J]. 信號處理, 2009, 25(6): 864-867.

    Huang Ya-jing, Cao Min, Fu Yao-wen,.. ISAR imaging of equably accelerative rotating targets based on matching Fourier transform[J]., 2009, 25(6): 864-867.

    [50] Wu Liang, Wei Xi-zhang, Yang De-gui,.. ISAR imaging of targets with complex motion based on discrete chirp Fourier transform for cubic chirps[J]., 2012, 50(10): 4201-4212.

    [51] Huang S, Lerner R M, and Parker K J. On estimating the amplitude of harmonic vibration from the Doppler spectrum of reflected signals[J]., 1990, 88(6): 2702-2712.

    [52] Greneker E F and Sylvester V B. Use of the envelope detection method to detect micro-Doppler[C]. Proceedings of SPIE, 2003, 5077: 167-174.

    [53] 蘇婷婷, 孔令講, 楊建宇. 基于多諧波微多普勒信號分析的目標攝動參數(shù)提取方法[J]. 電子與信息學報, 2008, 30(11): 2646-2649.

    Su Ting-ting, Kong Ling-jiang, and Yang Jian-yu. Target parameters extracted from multi-harmonic micro-Doppler signal analysis[J].&, 2008, 30(11): 2646-2649.

    [54] 霍凱, 李康樂, 姜衛(wèi)東, 等. 基于循環(huán)平穩(wěn)特征的正弦調制相位信號參數(shù)估計[J]. 電子與信息學報, 2010, 32(2): 355-359.

    Huo Kai, Li Kang-le, Jiang Wei-dong,.. Parameters estimation of signals with sinusoid modulated phase based on cyclostationary character[J].&, 2010, 32(2): 355-359.

    [55] 李開明, 李長棟, 李松, 等. 基于Gabor變換的微動目標微多普勒分析與仿真[J]. 空軍工程大學學報(自然科學版), 2010, 11(1): 40-43, 94.

    Li Kai-ming, Li Chang-dong, Li Song,.. Analysis and simulation for micro-Doppler information of micro-motion target based on Gabor transformation[J].(), 2010, 11(1): 40-43, 94.

    [56] 張云, 姜義成. 基于匹配傅里葉變換的艦船目標成像算法[J]. 大連海事大學學報, 2009, 35(1): 47-52.

    Zhang Yun and Jiang Yi-cheng. A method of ISAR ship imaging based on match Fourier transform[J]., 2009, 35(1): 47-52.

    [57] Pawan Setlur, Moeness Amin, and Fauzia Ahmad. Analysis of micro-Doppler signals using linear FM basis decomposition[C]. Proceedings of SPIE, 2006, 6210: 62100M.

    [58] Chen V C. Spatial and temporal independent component analysis of micro-Doppler features[C]. IEEE International Radar Conference, Virginia, USA, May 2005: 348-353.

    [59] Tao Ran, Li Yan-lei, and Wang Yue. Short-time fractional Fourier transform and its applications[J]., 2010, 58(5): 2568-2580.

    [60] Abramovich Y I, Spencer N K, and Turley M D E. Order estimation and discrimination between stationary and time- varying (TVAR) autoregressive models[J]., 2007, 55(6): 2861-2876.

    [61] 李松, 朱豐, 劉昌云, 等. 基于壓縮感知的彈道導彈微多普勒提取方法[J]. 電波科學學報, 2011, 26(5): 990-995.

    Li Song, Zhu Feng, Liu Chang-yun,.. Extraction of micro-Doppler of ballistic missile based on compressive sensing[J]., 2011, 26(5): 990-995.

    [62] Guo Jin-ku, Zou Hong-xing, Yang Xiao-jun,.. Parameter estimation of multicomponent chirp signals via sparse representation[J]., 2011, 47(3): 2261-2268.

    [63] Ma Liang, Liu Jin, Wang Tao,.. Micro-Doppler characteristics of sliding-type scattering center on rotationally symmetric target[J]., 2011, 54(9): 1957-1967.

    [64] Xu Jia, Yu Ji, Peng Ying-ning,.. Radon-Fourier transform for radar target detection, I: generalized Doppler filter bank[J]., 2011, 47(2): 1186-1202.

    Applications and Prospect of Micro-motion Theory in the Detection of Sea Surface Target

    Chen Xiao-long Guan Jian He You

    (Department of Electronic and Information Engineering, Naval Aeronautical and Astronautical University, Yantai 264001, China)

    Micro-Doppler signature is one of the physical characteristics of the target. The radar signature of a target with micro-motion can make fine characterizations of the shape, structure, and moving state of target, which reflects the nonstationary property of a radar signal. Hence, it has great superiority in the analysis of sea clutter and target detection in the case of high sea states based on the micro-Doppler theory. In this paper, to show the need for micro-Doppler, the modeling of scattering clutter from time-varying sea surface and analysis methods of sea clutter Doppler are first reviewed based on the principles and characteristics of micro-Doppler. Then, applications and technological approaches of micro-Doppler in sea surface target detection are introduced from the perspective of micro-motion target modeling and detection methods of micro-motion signatures. Finally, future research interests are highlighted based on problems experienced in present studies.

    Micro-motion; Micro-Doppler; Sea clutter; Target detection

    TN957.51

    A

    2095-283X(2013)01-0123-12

    10.3724/SP.J.1300.2013.20102

    陳小龍(1985-),男,山東省煙臺市,2010年獲海軍航空工程學院信息與通信工程學科碩士學位,現(xiàn)就讀于海軍航空工程學院,博士研究生,海軍航空工程學院研究生三隊,主要研究方向為雷達信號處理,海雜波中微弱目標檢測。獲全軍優(yōu)秀碩士論文獎。E-mail: cxlcxl1209@163.com

    關 鍵(1968-),男,遼寧省錦州市,2000年獲清華大學信息與通信工程學科博士學位,教授,博士生導師,海軍航空工程學院電子信息工程系主任,主要研究方向為雷達目標檢測與跟蹤、偵察圖像處理和信息融合。獲全國優(yōu)秀博士學位論文獎,新世紀百千萬人才工程國家級人選。E-mail: guanjian96@tsinghua.org.cn

    2012-12-17收到,2013-01-04改回;2013-01-09網(wǎng)絡優(yōu)先出版

    國家自然科學基金資助項目(61002045, 61179017, 61201445)和“泰山學者”建設工程專項經(jīng)費資助課題

    關鍵 guanjian96@tsinghua.org.cn

    猜你喜歡
    微動時頻雜波
    STAR2000型空管一次雷達雜波抑制淺析
    基于RID序列的微動目標高分辨三維成像方法
    雷達學報(2018年5期)2018-12-05 03:13:08
    基于稀疏時頻分解的空中目標微動特征分析
    密集雜波環(huán)境下確定性退火DA-HPMHT跟蹤算法
    微動橋橋足距離對微動裂紋萌生特性的影響
    相關廣義復合分布雷達海雜波仿真
    遙測遙控(2015年2期)2015-04-23 08:15:22
    基于時頻分析的逆合成孔徑雷達成像技術
    對采樣數(shù)據(jù)序列進行時頻分解法的改進
    雙線性時頻分布交叉項提取及損傷識別應用
    微波雷達海雜波混沌特性分析
    欧美少妇被猛烈插入视频| videosex国产| 亚洲,欧美,日韩| 老熟女久久久| 人妻 亚洲 视频| 久久国产精品男人的天堂亚洲| 一区二区日韩欧美中文字幕| 尾随美女入室| 国产片特级美女逼逼视频| 亚洲av日韩在线播放| av在线app专区| 一个人免费看片子| 久久久久久人人人人人| 亚洲av电影在线进入| 纯流量卡能插随身wifi吗| 99久久精品国产国产毛片| 国产黄色视频一区二区在线观看| 日韩一区二区三区影片| 国产精品蜜桃在线观看| a 毛片基地| 久久久久久久大尺度免费视频| 丝瓜视频免费看黄片| 国产精品香港三级国产av潘金莲 | 日韩欧美一区视频在线观看| 亚洲精品一区蜜桃| 国产成人一区二区在线| 精品国产一区二区三区久久久樱花| 日本-黄色视频高清免费观看| 男女啪啪激烈高潮av片| 欧美精品一区二区免费开放| 日韩中文字幕视频在线看片| 丰满乱子伦码专区| 一二三四中文在线观看免费高清| 男男h啪啪无遮挡| 在线 av 中文字幕| 各种免费的搞黄视频| 母亲3免费完整高清在线观看 | 亚洲av男天堂| 国产成人精品久久久久久| 精品福利永久在线观看| 一级片免费观看大全| 久久99精品国语久久久| 天堂中文最新版在线下载| 亚洲第一区二区三区不卡| 国产成人精品久久久久久| 久久99精品国语久久久| 又大又黄又爽视频免费| 人人妻人人澡人人看| 一级毛片 在线播放| 欧美日韩精品成人综合77777| 五月伊人婷婷丁香| 色视频在线一区二区三区| 黄网站色视频无遮挡免费观看| 成人午夜精彩视频在线观看| 18在线观看网站| 成人影院久久| 亚洲精品成人av观看孕妇| 国产精品国产av在线观看| a级毛片在线看网站| 免费av中文字幕在线| 人妻少妇偷人精品九色| 秋霞伦理黄片| 大香蕉久久成人网| 国产野战对白在线观看| 精品视频人人做人人爽| 2022亚洲国产成人精品| 国产av国产精品国产| 久热久热在线精品观看| 中文字幕人妻丝袜一区二区 | 亚洲欧洲日产国产| 黄色视频在线播放观看不卡| 在线观看一区二区三区激情| 制服丝袜香蕉在线| 亚洲欧洲精品一区二区精品久久久 | 啦啦啦在线观看免费高清www| 成人漫画全彩无遮挡| 欧美日韩亚洲国产一区二区在线观看 | 欧美日韩视频精品一区| 这个男人来自地球电影免费观看 | 97在线视频观看| 如日韩欧美国产精品一区二区三区| 婷婷色麻豆天堂久久| 极品少妇高潮喷水抽搐| 亚洲三级黄色毛片| 免费少妇av软件| 日韩欧美一区视频在线观看| 亚洲av日韩在线播放| 2018国产大陆天天弄谢| 久久久精品区二区三区| 黑丝袜美女国产一区| 亚洲av男天堂| 亚洲精品第二区| 精品国产一区二区三区四区第35| 亚洲精品日韩在线中文字幕| 中文字幕av电影在线播放| 免费黄网站久久成人精品| 最近中文字幕2019免费版| 可以免费在线观看a视频的电影网站 | 在线观看www视频免费| 伦精品一区二区三区| 亚洲熟女精品中文字幕| 天天躁日日躁夜夜躁夜夜| 亚洲国产毛片av蜜桃av| freevideosex欧美| 久久国产亚洲av麻豆专区| 如何舔出高潮| 国产有黄有色有爽视频| 国产片内射在线| 日韩中文字幕欧美一区二区 | 欧美变态另类bdsm刘玥| 如日韩欧美国产精品一区二区三区| 亚洲精品av麻豆狂野| 九九爱精品视频在线观看| 午夜福利在线免费观看网站| 2018国产大陆天天弄谢| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花| 大码成人一级视频| 国产精品免费视频内射| 久久女婷五月综合色啪小说| 免费av中文字幕在线| 亚洲精品在线美女| 免费观看av网站的网址| 男女高潮啪啪啪动态图| 成人二区视频| 人人妻人人澡人人爽人人夜夜| 亚洲av.av天堂| av不卡在线播放| 欧美在线黄色| 这个男人来自地球电影免费观看 | 亚洲精品国产av蜜桃| 色94色欧美一区二区| 成人国语在线视频| 水蜜桃什么品种好| 看免费成人av毛片| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| 久久毛片免费看一区二区三区| 国产精品成人在线| 午夜日韩欧美国产| 久久毛片免费看一区二区三区| 免费观看a级毛片全部| av卡一久久| 欧美精品高潮呻吟av久久| 亚洲精品自拍成人| 啦啦啦在线免费观看视频4| 久久精品国产鲁丝片午夜精品| 欧美xxⅹ黑人| 看免费av毛片| 在线观看人妻少妇| 午夜福利网站1000一区二区三区| 观看美女的网站| 丝瓜视频免费看黄片| 成人毛片60女人毛片免费| 精品一区在线观看国产| 亚洲av在线观看美女高潮| 成人毛片60女人毛片免费| 亚洲av电影在线观看一区二区三区| 日韩一区二区三区影片| 国产精品.久久久| 七月丁香在线播放| 韩国高清视频一区二区三区| 国产成人精品久久久久久| 制服人妻中文乱码| av线在线观看网站| 1024视频免费在线观看| 精品国产国语对白av| 成年女人在线观看亚洲视频| av免费观看日本| 又黄又粗又硬又大视频| 精品少妇黑人巨大在线播放| 成人国语在线视频| 国产无遮挡羞羞视频在线观看| 国产精品一区二区在线观看99| 在线观看三级黄色| 另类精品久久| 国产成人精品久久二区二区91 | 国产在视频线精品| 日韩,欧美,国产一区二区三区| 中文字幕人妻熟女乱码| 国产成人aa在线观看| 十八禁网站网址无遮挡| 欧美 亚洲 国产 日韩一| 黄色 视频免费看| 亚洲第一区二区三区不卡| 久久久久人妻精品一区果冻| 波多野结衣av一区二区av| 丝袜人妻中文字幕| 国产精品熟女久久久久浪| 亚洲,欧美,日韩| 日韩大片免费观看网站| 女的被弄到高潮叫床怎么办| 国产成人精品无人区| 久久人人爽av亚洲精品天堂| 久久国产精品男人的天堂亚洲| 制服诱惑二区| av电影中文网址| 午夜免费男女啪啪视频观看| 看非洲黑人一级黄片| 亚洲一区中文字幕在线| 色94色欧美一区二区| 欧美中文综合在线视频| av线在线观看网站| 亚洲天堂av无毛| 欧美日韩国产mv在线观看视频| 女人高潮潮喷娇喘18禁视频| 人人妻人人澡人人爽人人夜夜| 少妇被粗大的猛进出69影院| 麻豆精品久久久久久蜜桃| 久久精品国产鲁丝片午夜精品| 久久毛片免费看一区二区三区| 满18在线观看网站| 日韩一本色道免费dvd| 国产一级毛片在线| 亚洲av综合色区一区| 性高湖久久久久久久久免费观看| 日本-黄色视频高清免费观看| 亚洲av.av天堂| 寂寞人妻少妇视频99o| 欧美日韩亚洲国产一区二区在线观看 | 久久国产亚洲av麻豆专区| 一本久久精品| 亚洲欧美精品自产自拍| 久久久a久久爽久久v久久| 在线看a的网站| 色哟哟·www| 国产黄色视频一区二区在线观看| 欧美在线黄色| 大香蕉久久成人网| av一本久久久久| 国产精品熟女久久久久浪| 少妇的逼水好多| 免费高清在线观看视频在线观看| 午夜免费鲁丝| 欧美少妇被猛烈插入视频| 麻豆乱淫一区二区| 巨乳人妻的诱惑在线观看| 久久国产精品大桥未久av| 午夜激情久久久久久久| 麻豆精品久久久久久蜜桃| 久久久国产精品麻豆| 啦啦啦在线观看免费高清www| 亚洲三区欧美一区| 久久久国产精品麻豆| 国产成人精品婷婷| 久久久久久免费高清国产稀缺| 美女高潮到喷水免费观看| 亚洲熟女精品中文字幕| 亚洲,一卡二卡三卡| 亚洲国产欧美日韩在线播放| 在线天堂最新版资源| 午夜久久久在线观看| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 成人毛片a级毛片在线播放| 亚洲美女黄色视频免费看| 97在线视频观看| 婷婷成人精品国产| 亚洲第一青青草原| 精品国产一区二区久久| 春色校园在线视频观看| 国产精品国产av在线观看| 最黄视频免费看| 热re99久久国产66热| 制服人妻中文乱码| 久热久热在线精品观看| 国产又色又爽无遮挡免| 亚洲精品,欧美精品| 9色porny在线观看| 亚洲国产精品国产精品| 亚洲精品久久久久久婷婷小说| 99九九在线精品视频| 麻豆乱淫一区二区| 午夜福利一区二区在线看| 69精品国产乱码久久久| 婷婷色综合www| 人妻 亚洲 视频| 叶爱在线成人免费视频播放| 国产国语露脸激情在线看| 一级毛片电影观看| 免费在线观看视频国产中文字幕亚洲 | 免费在线观看视频国产中文字幕亚洲 | 人妻人人澡人人爽人人| 一个人免费看片子| 欧美最新免费一区二区三区| 亚洲精品自拍成人| 日本av手机在线免费观看| 日韩制服丝袜自拍偷拍| 五月开心婷婷网| 美女大奶头黄色视频| 亚洲欧洲日产国产| 黄色视频在线播放观看不卡| 国产日韩欧美在线精品| 免费观看av网站的网址| 一级片'在线观看视频| 一二三四中文在线观看免费高清| 午夜精品国产一区二区电影| 亚洲经典国产精华液单| 日本免费在线观看一区| 国产一区二区三区av在线| 天美传媒精品一区二区| 亚洲精品成人av观看孕妇| 国产xxxxx性猛交| 最近最新中文字幕大全免费视频 | 美女国产视频在线观看| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 国产深夜福利视频在线观看| 久久久久精品久久久久真实原创| 久久久久久久国产电影| 一区在线观看完整版| 青春草国产在线视频| 老鸭窝网址在线观看| 免费观看无遮挡的男女| 国产淫语在线视频| 亚洲人成电影观看| 亚洲av日韩在线播放| 老司机影院毛片| 欧美人与性动交α欧美精品济南到 | 制服人妻中文乱码| 飞空精品影院首页| 国产高清不卡午夜福利| 成人影院久久| 精品亚洲成a人片在线观看| 只有这里有精品99| 午夜老司机福利剧场| 亚洲国产精品999| 亚洲国产av新网站| 亚洲欧美日韩另类电影网站| 蜜桃在线观看..| 精品国产乱码久久久久久小说| 一二三四在线观看免费中文在| 亚洲中文av在线| 亚洲伊人久久精品综合| 中国三级夫妇交换| 成人影院久久| 亚洲美女搞黄在线观看| 午夜福利在线免费观看网站| 波多野结衣一区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 精品卡一卡二卡四卡免费| 亚洲精品在线美女| 久久精品人人爽人人爽视色| 国产亚洲欧美精品永久| 最近2019中文字幕mv第一页| 国产伦理片在线播放av一区| 少妇人妻 视频| 人妻人人澡人人爽人人| 另类亚洲欧美激情| 大香蕉久久成人网| 久久99热这里只频精品6学生| 九色亚洲精品在线播放| 欧美老熟妇乱子伦牲交| 国产成人精品无人区| 国产亚洲欧美精品永久| a 毛片基地| www.av在线官网国产| 成人国产av品久久久| 欧美人与性动交α欧美精品济南到 | 女性被躁到高潮视频| 超碰97精品在线观看| 夫妻性生交免费视频一级片| 欧美日韩精品成人综合77777| 综合色丁香网| 国产欧美日韩一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 国产精品免费视频内射| 精品国产露脸久久av麻豆| 日本爱情动作片www.在线观看| 国产日韩欧美在线精品| 亚洲欧美成人综合另类久久久| 9191精品国产免费久久| 亚洲精品久久成人aⅴ小说| 日本vs欧美在线观看视频| 国产亚洲一区二区精品| 大香蕉久久网| 十分钟在线观看高清视频www| 国精品久久久久久国模美| 人体艺术视频欧美日本| 亚洲av.av天堂| 亚洲精品第二区| 国产精品一区二区在线观看99| 看非洲黑人一级黄片| 2021少妇久久久久久久久久久| 美女视频免费永久观看网站| 色94色欧美一区二区| 精品久久蜜臀av无| 黄网站色视频无遮挡免费观看| 欧美激情 高清一区二区三区| 国产精品偷伦视频观看了| 美女福利国产在线| 精品人妻熟女毛片av久久网站| xxx大片免费视频| 国产在线一区二区三区精| 久久久久久伊人网av| 免费黄频网站在线观看国产| 一级片'在线观看视频| 一区福利在线观看| 久久热在线av| 五月开心婷婷网| 久久午夜福利片| 欧美日韩综合久久久久久| 亚洲婷婷狠狠爱综合网| 成年美女黄网站色视频大全免费| 国产精品99久久99久久久不卡 | √禁漫天堂资源中文www| 日韩欧美一区视频在线观看| 黄片播放在线免费| 性高湖久久久久久久久免费观看| 国产在视频线精品| 精品人妻偷拍中文字幕| 国产精品一二三区在线看| 亚洲天堂av无毛| 999久久久国产精品视频| 亚洲欧美成人综合另类久久久| 女人高潮潮喷娇喘18禁视频| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| 久久久久久久精品精品| 亚洲av电影在线进入| 午夜91福利影院| 永久网站在线| 亚洲欧美日韩另类电影网站| 观看av在线不卡| 女的被弄到高潮叫床怎么办| 女人精品久久久久毛片| 男女无遮挡免费网站观看| 男女高潮啪啪啪动态图| 国产在线视频一区二区| 午夜91福利影院| 亚洲欧美一区二区三区黑人 | 日韩中字成人| 亚洲少妇的诱惑av| 午夜激情久久久久久久| 亚洲第一青青草原| 美女福利国产在线| 满18在线观看网站| 国产精品无大码| 久久人人97超碰香蕉20202| 99久久中文字幕三级久久日本| 国产成人精品一,二区| 日本猛色少妇xxxxx猛交久久| 国产又色又爽无遮挡免| 一二三四中文在线观看免费高清| 免费女性裸体啪啪无遮挡网站| 777久久人妻少妇嫩草av网站| 国产极品天堂在线| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕免费大全7| 成年动漫av网址| 国产一区二区 视频在线| 久久久国产一区二区| 中国国产av一级| 超碰97精品在线观看| 成人黄色视频免费在线看| 婷婷成人精品国产| 日韩av在线免费看完整版不卡| 水蜜桃什么品种好| 女人久久www免费人成看片| 观看美女的网站| 久久亚洲国产成人精品v| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 日本vs欧美在线观看视频| 97在线人人人人妻| 久久这里只有精品19| 国产深夜福利视频在线观看| 日韩精品有码人妻一区| 亚洲国产成人一精品久久久| 如日韩欧美国产精品一区二区三区| 亚洲精品在线美女| 免费高清在线观看视频在线观看| 黄网站色视频无遮挡免费观看| 免费人妻精品一区二区三区视频| 建设人人有责人人尽责人人享有的| 男女啪啪激烈高潮av片| 五月天丁香电影| 久久人妻熟女aⅴ| 不卡av一区二区三区| 欧美日韩综合久久久久久| 咕卡用的链子| 亚洲人成77777在线视频| 久久久久久久久免费视频了| 亚洲人成网站在线观看播放| 黑人猛操日本美女一级片| 考比视频在线观看| 日本vs欧美在线观看视频| av在线老鸭窝| 99热网站在线观看| 久久久国产欧美日韩av| 国产成人精品婷婷| 久热这里只有精品99| 精品一区二区免费观看| 999精品在线视频| 五月伊人婷婷丁香| www.自偷自拍.com| 视频区图区小说| 久久精品国产鲁丝片午夜精品| 菩萨蛮人人尽说江南好唐韦庄| av线在线观看网站| 国产成人91sexporn| 精品亚洲成a人片在线观看| 色吧在线观看| 欧美在线黄色| 日韩av免费高清视频| 精品人妻在线不人妻| 免费久久久久久久精品成人欧美视频| 80岁老熟妇乱子伦牲交| 免费av中文字幕在线| 深夜精品福利| 久久99热这里只频精品6学生| 黄色一级大片看看| 亚洲国产欧美网| 亚洲成色77777| 日日摸夜夜添夜夜爱| 国产精品一区二区在线观看99| 国产成人a∨麻豆精品| 国产在线免费精品| 一二三四中文在线观看免费高清| www.av在线官网国产| 精品国产国语对白av| 日本免费在线观看一区| 精品国产国语对白av| 九草在线视频观看| 国产av国产精品国产| www.av在线官网国产| 午夜91福利影院| 欧美日本中文国产一区发布| 大片免费播放器 马上看| 国产高清不卡午夜福利| 亚洲国产欧美网| 日本免费在线观看一区| 午夜激情av网站| 2022亚洲国产成人精品| 飞空精品影院首页| 视频区图区小说| 亚洲伊人色综图| 性少妇av在线| tube8黄色片| 激情五月婷婷亚洲| a级片在线免费高清观看视频| 黄色怎么调成土黄色| 久久精品人人爽人人爽视色| 亚洲综合精品二区| 天天躁夜夜躁狠狠久久av| 精品国产超薄肉色丝袜足j| 精品国产乱码久久久久久小说| 精品亚洲成国产av| 又粗又硬又长又爽又黄的视频| 少妇的逼水好多| 看免费av毛片| 国产欧美亚洲国产| videossex国产| 免费观看av网站的网址| 午夜激情av网站| 男人操女人黄网站| 99久久人妻综合| 亚洲一级一片aⅴ在线观看| 青春草视频在线免费观看| 青春草亚洲视频在线观看| 国产伦理片在线播放av一区| 亚洲精品av麻豆狂野| 欧美最新免费一区二区三区| 美女国产高潮福利片在线看| 一区二区三区精品91| 午夜精品国产一区二区电影| 日日撸夜夜添| 80岁老熟妇乱子伦牲交| 久久99热这里只频精品6学生| 美女脱内裤让男人舔精品视频| 午夜激情av网站| 制服诱惑二区| 天天躁夜夜躁狠狠躁躁| 国产一区亚洲一区在线观看| 亚洲欧美色中文字幕在线| 丰满迷人的少妇在线观看| videosex国产| 美国免费a级毛片| 亚洲av电影在线进入| 丝袜脚勾引网站| 久久精品亚洲av国产电影网| 久久精品国产亚洲av天美| 亚洲av男天堂| 国产又色又爽无遮挡免| 热re99久久国产66热| 美女大奶头黄色视频| 狠狠精品人妻久久久久久综合| 国产一区二区在线观看av| 成年女人在线观看亚洲视频| 青青草视频在线视频观看| 99国产综合亚洲精品| 国产日韩欧美视频二区| 26uuu在线亚洲综合色| 久久狼人影院| 久久热在线av| 国产有黄有色有爽视频| 亚洲精品成人av观看孕妇| 一区在线观看完整版| 日韩制服丝袜自拍偷拍| 色哟哟·www| 亚洲精品美女久久av网站| 热99国产精品久久久久久7| 国产精品一国产av| 免费黄频网站在线观看国产| 丰满迷人的少妇在线观看| 久久精品人人爽人人爽视色| 亚洲欧美一区二区三区国产| 亚洲欧洲国产日韩| 亚洲第一av免费看| 在线观看美女被高潮喷水网站| 亚洲中文av在线| 曰老女人黄片| 26uuu在线亚洲综合色| 久久婷婷青草| 欧美日韩亚洲国产一区二区在线观看 |