亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于改進希爾伯特-黃變換的進動目標周期估計

        2022-07-25 03:51:54金家偉阮懷林
        火力與指揮控制 2022年4期
        關鍵詞:進動希爾伯特時頻

        金家偉,阮懷林,孫 兵

        (國防科技大學電子對抗學院,合肥 230037)

        0 引言

        物體或物體結(jié)構(gòu)的微動可能會對返回的信號產(chǎn)生頻率調(diào)制,從而產(chǎn)生多普勒頻移的邊帶,稱為微多普勒效應。彈頭目標為了保持飛行的穩(wěn)定性,不僅會繞自身對稱軸作自旋運動,還會受到?jīng)_擊力矩的作用,力矩消失后對稱軸將在平衡位置作圓錐運動,即進動。進動是目標保持穩(wěn)定的一種特殊微動,而誘餌一般不具備姿態(tài)控制裝置,因此,進動是中段彈頭特有的微動方式。通過提取具有進動錐體目標的幾何特征,可以區(qū)分它是彈頭還是誘餌。

        至今為止,由于時頻分析可將一維信號映射到二維時頻平面,從而展現(xiàn)信號的時頻信息,因此時頻分析技術被廣泛應用于分析時變的微多普勒頻率特征,因為時頻分析可以通過將一維信號映射到二維時頻平面來同時產(chǎn)生信號的時頻信息。文獻[6]結(jié)合Gabor 時頻分布和變分模態(tài)分解,來估計目標的自旋頻率和錐旋頻率;文獻[7]利用短時分數(shù)階傅立葉變換分離肢體和軀干的微多普勒信號;文獻[8]采用基于加權迭代自適應的時頻分析方法結(jié)合逆Radon 變換,分離重構(gòu)不同散射點的微多普勒分量。然而,在實際應用中,大多數(shù)現(xiàn)有的時頻方法在提取微多普勒信息方面仍然存在一些不可避免的不足,包括邊界失真和干擾項。另外,當目標含有多個強散射中心時,微動目標返回的雷達信號是多分量的。因此,所需的時頻分析方法應該同時抑制交叉項和自動項。

        近年來,一種新的分析非線性、非平穩(wěn)數(shù)據(jù)的方法——希爾伯特- 黃變換(hilbert-huang transform,HHT)越來越受到人們的重視。該方法的關鍵部分是經(jīng)驗模態(tài)分解(empirical mode decomposition,EMD),其主要思想是將需要進行處理的非線性非平穩(wěn)信號分解,得到一組頻率分量各有不同的本征模態(tài)函數(shù)(intrinsic mode functions,IMF)。通過希爾伯特變換(hilbert transform,HT),IMF 能夠產(chǎn)生隨時間變化的瞬時頻率,其最終表示形式是能量-頻率- 時間分布,稱為希爾伯特- 黃譜(hilberthuang spectrum,HHS)。希爾伯特-黃變換的優(yōu)點之一是希爾伯特- 黃譜不涉及頻率分辨率和時間分辨率的概念,而是瞬時頻率。因其良好的適應性,已被廣泛應用于微多普勒效應的研究,文獻[12]使用了希爾伯特-黃變換對振動目標進行雷達微多普勒特征分析;文獻[13]提出了一種基于經(jīng)驗模態(tài)分解的圓錐目標雷達微動特征提取方法。然而,在實際應用中,EMD 算法并不完美,同樣存在一些問題,如模態(tài)混疊、端點效應等。

        針對基于希爾伯特- 黃變換的進動目標微多普勒特征提取問題,由于雷達信號不可避免地含有較高水平的隨機噪聲,其瞬時頻率和瞬時幅度非常分散,常規(guī)希爾伯特-黃變換難以將目標信號與噪聲信號分離。本文采用EMD 的衍生算法CEEMDAN先對雷達回波進行分解,然后再進行希爾伯特譜分析,以此提取進動目標的微多普勒周期。

        1 目標模型

        空間進動錐體目標模型如圖1 所示,坐標系為雷達坐標系(U,V,W),雷達靜止于原點Q。O 為目標質(zhì)心,以O 為原點、目標對稱軸Oz 為z 軸建立目標本體坐標系O-xyz。以O 為原點建立參考坐標系O-XYZ,以初始時刻與目標對稱軸Oz、進動軸OZ共面且垂直于OZ 的方向為Y 軸,X 軸根據(jù)右手準則確定。目標在平動的同時,以角速度ω繞對稱軸z軸作自旋運動,同時以角速度ω繞OZ 軸錐旋運動(ω和ω均采用參考坐標系中的表達式),自旋軸和錐旋軸之間的夾角為進動角θ。

        圖1 目標模型

        在光學區(qū),雷達目標的整體散射特性通??梢缘刃槿舾蓚€散射中心的疊加。不失一般性,假設目標的微動是周期性進動,目標等效為K 個散射中心,各散射中心各向同性,雷達發(fā)射的電磁波為連續(xù)單頻波,載頻為f。雷達向由i 個散射點組成的目標發(fā)射電磁波,不考慮目標平動帶來的影響,返回的基帶信號可表示為

        2 改進的希爾伯特-黃變換

        2.1 希爾伯特-黃變換

        2.2 完備總體經(jīng)驗模態(tài)分解

        Torres 等提出了一種對噪聲有自適應分解特性的完備總體經(jīng)驗模態(tài)分解(CEEMDAN),在每一階段的分解中都加入一個特定的白噪聲,再通過計算殘差得到該階段的模態(tài)。CEEMDAN 是一種噪聲輔助數(shù)據(jù)分析方法,所加入的白噪聲在整個時頻空間分布相對均勻,以提供時頻空間中濾波的參考尺度,而不同尺度的信號自動投影到由白噪聲建立的合適的參考尺度上。雖然每次試驗產(chǎn)生的結(jié)果都會包含噪聲,但是當試驗的數(shù)量足夠時,其結(jié)果的總體平均當中的噪聲將消失,即總體平均當中所添加的白噪聲往往會相互抵消,從而只包含所分解得到的信號。

        3 仿真與分析

        為驗證希爾伯特- 黃變換在進動目標微多普勒特征提取方面的效果,本文分別以EMD 和CEEMDAN 來進行信號分解,然后進行希爾伯特譜分析,并對比兩種方法的提取效果。

        仿真條件:假設雷達工作在10 GHz,且雷達坐標系中本體坐標系原點O 的坐標為(400,500,100)km,本地坐標系和參考坐標系之間的初始歐拉角(x-y-z 序列)為(30°,60°,45°)。假設目標繞z 軸旋轉(zhuǎn),目標上有兩個散射點:第1 散射點A 位于錐頂,在本體坐標系中的坐標為(0,0,1)m;第2 個散射點B 位于錐底的尾翼,在本體坐標系中的坐標為(0.5,0,-0.5)m。雷達照射時間為8 s。旋轉(zhuǎn)頻率為f=1 Hz,圓錐運動頻率為f=0.5 Hz。假設環(huán)境噪聲為高斯白噪聲。

        當不存在噪聲時,基于EMD 和CEEMDAN 的進動周期提取結(jié)果分別如圖2 和圖3 所示,圖2(a)和圖3(a)分別是兩種算法獲取的回波信號希爾伯特譜,圖2(b)和圖3(b)分別是兩種算法獲取的IMF前3 個分量的希爾伯特譜。

        圖2 基于EMD 的估計結(jié)果(無噪聲)

        希爾伯特譜是時間、頻率和能量的三維譜,當沒有噪聲時,圖2(a)和圖3(a)都展現(xiàn)了回波信號時頻分布的細節(jié),其中的頻率都是取模后的頻率,也就是正頻率。微多普勒信號經(jīng)過EMD 后都被分解在第1 個IMF 分量中,如圖2(b)所示;而微多普勒信號經(jīng)過CEEMDAN 后被分解在各個IMF 分量中,如圖3(b)所示。由圖2 和圖3 可看出,兩種算法在不存在噪聲時,都能夠準確估計出進動目標的微多普勒周期為2 s。

        圖3 基于CEEMDAN 的估計結(jié)果(無噪聲)

        當雷達回波信號中混雜著噪聲時,同樣的仿真條件下,假設信噪比SNR=10 dB,圖4(a)為基于EMD 的回波信號希爾伯特譜,圖4(b)是基于EMD得到的IMF 前3 個分量的希爾伯特譜。

        圖4 基于EMD 的估計結(jié)果(SNR=10 dB)

        從圖4 可以看出,隨著信噪比的降低,信號的希爾伯特譜中的微多普勒頻率逐漸淹沒在噪聲中,無法用以提取進動目標微多普勒周期,如圖4(a)所示。由于EMD 未能將微多普勒信號和噪聲分離開,而是都混在第1 個IMF 分量中,如圖4(b)所示,隨著信噪比SNR 的降低,將同樣無法通過第1 個IMF分量的希爾伯特譜來提取目標信息。但在其他條件不變的情況下,改用CEEMDAN 對雷達回波進行處理,然后進行希爾伯特譜分析。圖5(a)為回波信號希爾伯特譜,圖5(b)為采用CEEMDAN 所獲取的IMF 前3 個分量的希爾伯特譜。

        圖5 基于CEEMDAN 的估計結(jié)果(SNR=10 dB)

        從圖5(a)中可以看出,當信噪比持續(xù)降低時,與圖4(a)相似,其中的微多普勒信號伴隨著大量的噪聲,難以直接通過回波信號地希爾伯特譜來提取微多普勒信息。但通過CEEMAN 分解后所得到的一組IMF,噪聲頻率主要集中于第1 個IMF 分量中,雖然第1 個IMF 分量中也包含了一部分的微多普勒頻率信號,但更多的微多普勒信號被分解到了第2 個和第3 個分量中,從而達到將噪聲和目標信號分離的目的,如圖5(b)所示,則從該圖中的第2 個和第3 個IMF 分量的希爾伯特譜中,可以準確提取出目標的微多普勒周期。

        進一步降低信噪比進行對比分析,圖6(a)和圖6(b)分別給出了信噪比分別為5 dB 和1 dB 時,基于EMD 得到的IMF 前3 個分量的希爾伯特譜。

        圖6 EMD 獲取的IMF 前3 個分量的希爾伯特譜

        由于EMD 算法難以從噪聲分離出微多普勒信號,在其IMF 分量的希爾伯特譜中,即使信噪比為5 dB,仍然無法估計出進動周期,如圖6(a)所示。但采用CEEMAN 時,圖7(a)和圖7(b)分別給出了信噪比分別為5 dB 和1 dB 時,基于CEEMDAN 得到的IMF 前3 個分量的希爾伯特譜。

        圖7 CEEMDAN 獲取的IMF 前3 個分量的希爾伯特譜

        將圖7 與圖6 對比可知,隨著信噪比的進一步降低,即使信噪比低至1 dB,仍能通過其IMF 分量的希爾伯特譜提取出進動周期,如圖7(b)所示。這也進一步驗證了在低信噪比環(huán)境中,以CEEMDAN改進的希爾伯特-黃變換具備更加良好的性能。仿真中假設環(huán)境噪聲為高斯白噪聲,旨在驗證CEEMDAN 具備比EMD 更好的抗噪性能,但在實際工程中的噪聲不一定是高斯白噪聲,算法具體的抗噪性能需要以實際雷達回波數(shù)據(jù)進行驗證,這是本文算法存在的局限性。

        4 結(jié)論

        本文將CEEMDAN 算法替換EMD 算法,以改進希爾伯特-黃變換,用以提取進動目標的微多普勒周期。仿真實驗表明,改進的希爾伯特-黃變換能夠獲取進動目標微多普勒信號的希爾伯特譜,從而提取出目標的進動周期;相對于EMD 算法,CEEMDAN 算法能夠?qū)⑽⒍嗥绽招盘柡驮肼曅盘栠M行分離,從而有效克服了常規(guī)希爾伯特-黃變換在噪聲環(huán)境中性能較差的不足。

        猜你喜歡
        進動希爾伯特時頻
        一個真值函項偶然邏輯的希爾伯特演算系統(tǒng)
        邏輯學研究(2021年3期)2021-09-29 06:54:34
        基于窄帶雷達網(wǎng)的彈道目標三維進動特征提取
        雷達學報(2017年1期)2017-05-17 04:48:57
        下一個程序是睡覺——數(shù)學家希爾伯特的故事
        進動錐體目標平動補償及微多普勒提取
        基于希爾伯特-黃變換和小波變換的500kV變電站諧振數(shù)據(jù)對比分析
        電測與儀表(2016年7期)2016-04-12 00:22:14
        基于希爾伯特- 黃變換的去噪法在外測數(shù)據(jù)處理中的應用
        基于時頻分析的逆合成孔徑雷達成像技術
        基于雷達距離像的錐體目標進動參數(shù)估計方法
        對采樣數(shù)據(jù)序列進行時頻分解法的改進
        雙線性時頻分布交叉項提取及損傷識別應用
        亚洲av午夜福利精品一区不卡 | 尤物AV无码色AV无码麻豆 | 精品国产av一区二区三区 | 中国人妻与老外黑人| 国产精品成人免费视频网站京东| 日韩精品欧美激情亚洲综合| 日本一区二区午夜视频| 91九色免费视频网站 | 亚洲精品无码高潮喷水在线 | 狠狠躁天天躁无码中文字幕图| 亚洲视频综合在线第一页 | 夜夜躁狠狠躁日日躁视频| 国产精品福利视频一区| 久久道精品一区二区三区| 视频一区精品中文字幕| 欧美肥妇毛多水多bbxx水蜜桃| 日韩少妇激情一区二区| 日韩久久久黄色一级av| 草青青视频手机免费观看| 久久精品国产字幕高潮| 欧美日韩精品一区二区在线观看 | 亚洲av无码片在线观看| 大屁股少妇一区二区无码| 久久久人妻精品一区bav| 亚洲欧美国产精品久久| 久久国产精品久久精品国产| 亚洲一二三四五区中文字幕 | 无码人妻丰满熟妇区五十路百度| 久久国产成人午夜av影院| 国产一区二区三区白浆肉丝| 国产日韩精品suv| 九九视频在线观看视频6| 亚洲无线码1区| 91偷自国产一区二区三区| 亚洲精品久久久久成人2007| 欧美成人免费高清视频| 日韩精品少妇专区人妻系列| 久久午夜av一区二区三区| 国产在线观看www污污污| 国产69精品一区二区三区| av手机免费在线观看高潮|