實測海雜波數據的多重分形特性分析

王 磊，金 丹，察 豪

(1.華中科技大學自動化學院 海裝信息系統局,北京 100072) (2.海軍工程大學 電子工程學院，武漢 430033)

?

實測海雜波數據的多重分形特性分析

王 磊1，金 丹2，察 豪2

(1.華中科技大學自動化學院 海裝信息系統局,北京 100072) (2.海軍工程大學 電子工程學院，武漢 430033)

為充分描述海雜波時變特性和局部奇異性，將多重分形理論應用到雷達海雜波特性分析中，理論分析了海雜波的多重分形參數，并進一步分析了實測海雜波數據的多重分形譜、質量指數和奇異指數，結果表明海雜波具有明顯的多重分形特性，該分析方法可以很好地解釋海雜波內在非線性特性，并比傳統的統計特性分析方法更加能確切地描述雜波的產生機理，為基于多重分形特性的雷達目標檢測提供了依據。

海雜波；多重分形；數據分析；質量指數；奇異指數

0 引 言

隨著雷達系統分辨率的提高，我們觀測到的海雜波的特性也不再符合傳統的基于統計概率的分布規律，對海雷達面臨的海雜波表現出更加復雜的非高斯性、非平穩性和非線性等特性。1982年Mandelbrot[1]提出了研究自相似性和標度不變性的分形理論，作為非線性科學研究領域的重要分支，并應用于海雜波研究領域。文獻[2]通過對實測海雜波數據的分析，證實了海雜波的時間序列具有分形特性，但是，文章僅用單純的分維數來描述時變的海雜波及其特性，顯然不具有全面性。多重分形測度分析可以從多角度描述海雜波的復雜特性，并可以揭示海雜波系統復雜的非均勻結構。文獻[3]對實測海雜波數據進行了多重分形測度分析，并進一步提取了海雜波復雜特性的多重分形參數，但是，對多重分形譜函數的描述及分析過于簡單。

本文基于實測雷達數據，分析了其多重分形特性，針對單純分維數的簡單分形理論分析海雜波的局限性，將多重分形理論應用到海雜波分析中，研究了不同條件下的多重分形參數及其特征的變化特點，并分析其內在原因，為多重分形理論在雷達目標檢測中的應用奠定基礎。

1 多重分形理論與建模分析

1.1 多重分形理論

在分形理論中，待研究對象是由單一分形維數來描述的，但是在大多數實際現象和生活中，單一的分形維數對研究對象的描述過于簡單化。多重分形是分形的推廣,多重分形是將分形結構上的多個標度指數進行綜合表述，主要用來描述信號局部和全局的奇異性，可以全面細致的刻畫海雜波的分形特性[4-5]。

1.2 多重分形模型

我們把研究對象劃分為N個不同的區域si(i=1,2,…,N)，設ri為第i個區域線度的大小，pi為該區域si的生成概率，則不同的區域si，pi也不相同，可以用不同的標度指數αi來表示[6]

(1)

若線度大小趨于0，則上式化為

(2)

式中：α是研究對象在某區域的分維數，即局部分維數。

多重分形分析中，α表示研究對象在小區域內的分維數，當小區域數目很大時，可以得到一個隨α變化而變化的函數，我們稱之為多重分形譜[7-8]，并用f(α)表示。

從信息論角度也可以選用q和Dq來描述多重分形。把式(1)兩邊各自乘q次方并取和得

(3)

式中：取區域線性大小ri=L。

q次信息維Dq的定義為

(4)

這兩組參量之間的聯系為Legendre變換[9]

(5)

f(α)=qα-τ(q)

(6)

其中

τ(q)=(q-1)Dq

(7)

(8)

式中：τ(q)為質量指數，與q呈非線性關系，多重分形性的可以用此非線性程度表示。α(q)為奇異指數，是質量指數τ(q)與q非線性程度的表示，是多重分形性程度強弱的表示。f(α)為多重分形譜，反映了待研究對象在分形結構上的非均勻分布的特性，它可以給出比單純分形分維更豐富的內部結構信息。多重分形譜的寬度Δα=αmax-αmin反映了在整個分形結構上的待研究對象的概率測度分布的非均勻性、非線性的程度。其中，α值較大反映的是小概率測度區域的性質，α值較小則反映的是大概率測度區域的性質。多重分形譜函數可以全面的描述待研究對象在整個分形結構上的不同層次結構、不同區域特點、不同局域條件的特性。因此多重分形譜函數f(α)是對待研究對象在整個分形結構上的內部復雜程度、非均勻程度以及非線性程度的一種度量。

1.3 多重分形譜函數的參數特征分析

多重分形譜函數f(α)反映了待研究對象在整個分形結構上的非線性、非均勻分布的性質。多重分形譜函數通常為一光滑的單峰函數，其形狀為鉤狀或鐘罩狀。為更好的分析多重分形譜函數，用最小二乘法對f(α)-α進行擬合

f(α)=A(α-α0)2+B(α-α0)+C

(9)

其中

f(α0)=fmax

(10)

式中：由式(9)、式(10)可得，αmax反映的是最小概率測度區域的性質；αmin反映的是最大概率測度區域的性質；α0越大，分析過程細節越豐富。多重分形譜寬Δα反映了待研究對象在整個分形結構上的概率測度分布的非線性、非均勻性的程度，Δα越大，譜所包含的信息量越大，譜的結構越豐富。若B<0，則譜的頂點向右側傾斜，呈右鉤狀，這表示α較大(小概率)部分占優勢；若B>0，則譜的頂點向左側傾斜，呈左鉤狀，這表示α較小(大概率)部分占優勢。

2 實測數據介紹及仿真

本文的研究基于實測海雜波數據，數據來源于McMaster University公布的IPIX雷達數據，其測量時間為1993年11月。關于測量數據和IPIX雷達的詳細情況，可參考文獻[10]和McMaster University的網站(soma.ece.mcmaster.ca)。表1給出了所選用的數據的基本條件和相關參數，由于本文研究的是海雜波特性，因此只對不含有目標的雜波距離單元進行研究分析。

表1 海雜波數據基本參數

圖1給出了表1中一組海雜波數據的歸一化幅度波形圖。對海雜波時域波形進行多重分形分析，圖2是q為-30～30，步長為5的雙對數關系圖，如圖2所示，配分函數和對數尺度整體上表現出了良好的線性關系，并滿足一定的冪率關系，這說明在特定的尺度范圍內實測海雜波數據具有良好的無標度性，也就是說通過分析表明實測海雜波具有良好的多重分形特性。

圖1 歸一化海雜波時域波形圖

圖2 配分函數lb2X(q)與對數尺度lb2L的關系

由式(12)和圖3可得，圖像頂點右偏，即B<0，小概率部分占優勢，海雜波測度場的粗糙程度較高。

圖3 多重分形譜函數

圖4 質量指數關系圖

從圖4可以看出，海雜波數據的質量指數函數τ(q)是q上的一個上凸的函數，即質量指數τ(q)與q之間存在著較強的非線性關系，通過分析可以得出與圖2一致的結論，即實測海雜波數據具有一定的多重分形特性。而圖5奇異指數α(q)是質量指數τ(q)與q呈非線性程度的度量，由圖可得，實測海染波數據具有一定強度的多重分形特性，這與我們理論分析的情形一致。

圖5 奇異指數關系圖

3 結束語

本文重點討論了海雜波的多重分形特性，特別是對其多重分形譜、質量指數和奇異指數進行了分析。得出了實測海雜波數據的多重分形特性圖，并結合理論模型和實測數據分析了參數變化特點，結果表明海雜波的多重分形特性分析比傳統的統計特性分析模型(瑞利分布、對數正態分布、K分布等)更加能夠描述海雜波的產生機理，揭示了海雜波復雜的非線性的結構，這對于雷達海雜波的建模與仿真以及海雜波背景下的小目標檢測具有現實意義。

下一步的研究重點是對海雜波的入射角、極化方式、風速、風向以及海情等不確定因素的全面深入研究，為精確建立海雜波模型奠定理論基礎，對研究海雜波特性及目標檢測提供理論支撐。

[1] MANDELBROT B B.The fractal geometry of nature[M].New York:W.H.Freeman and Company,1982.

[2] HAYKIN S,BAKKER R,CURRIE B W.Uncovering nonlinear dynamics the case study of sea clutter[J].Proceedings of the IEEE,2002,90(5):860-881.

[3] 石志廣,周建雄,趙宏鐘,等.海雜波的多重分形特性分析[J].數據采集與處理,2006,21(2):168-173.

SHI Zhiguang,ZHOU Jianxiong,ZHAO Hongzhong,et al.Multifractal analysis of radar sea clutter[J].Journal of Data Acquisition and Processing,2006,21(2):168-173.

[4] 盛 文,任 吉.高頻雷達海雜波的多重分形建模方法研究[J].電波科學學報,2011,26(5):983-989.

SHENG Wen,REN Ji.Reaearch on multifractal modeling of high frequency radar sea clutter[J].Chinese Journal of Radio Science,2011,26(5):983-989.

[5] 關 鍵,劉寧波,張 建,等.海雜波的多重分形關聯特性與微弱目標檢測[J].電子與信息學報,2010,32(1):54-61.

GUAN Jian,LIU Ningbo,ZHANG Jian,et al.Multifractal correlation characteristic of real sea clutter and low-observable targets detection[J].Journal of Electronics & Information Technology,2010,32(1):54-61.

[6] 王祖林,周蔭清.多重分形譜及其計算[J].北京航空航天大學學報,2000,26(3):256-258.

WANG Zulin,ZHOU Yinqing.Multifractal spectrum and calculation[J].Journal of Beijing University of Aeronautics and Astronautics,2000,26(3):256-258.

[7] 劉元永,羅曉曙,陳全斌,等.多重分形譜在葉片圖像處理中的應用[J].計算機工程與應用,2008,44(28):190-192.

LIU Yuanyong,LUO Xiaoshu,CHEN Quanbin,et al.Application of multifractal spectrum in leaf images processing[J].Computer Engineering and Applications,2008,44(28):190- 192.

[8] 劉小剛,趙慧潔,李 娜.基于多重分形譜的高光譜數據特征提取[J].光學學報,2009,29(3):844-848.

LIU Xiaogang,ZHAO Huijie,LI Na.Feature extraction based on multifractal spectrum for hyperspectral data[J].Acta Optica Sinica,2009,29(3):844-848.

[9] 熊 剛,張淑寧,趙慧昌.基于小波leaders的海雜波時變奇異譜分布分析[J].物理學報,2014,63(15):150503,1-7.

XIONG Gang,ZHANG Shuning,ZHAO Huichang.Time-varying singularity spectrum distribution of sea clutter based on wavelet leaders[J].Acta Physica Sinica,2014,63(15):150503,1-7.

[10] DROSOPOULUS A.Description of the OHGR database[R].DREO-TN-94-14.Canada.Defense Research Establishment Ottawa Technical Note,1994.

王 磊 男，1982年生，碩士。研究方向為雷達工程。

金 丹 女，1986年生，博士。研究方向為信號與信息處理。

察 豪 男，1966年生，博士，教授。研究方向為雷達工程。

Multifractal Analysis of Measured Sea Clutter Data

WANG Lei1，JIN Dan2，CHA Hao2

(1.Huazhong University of Science and Technology,Marine Information System,Beijing 100071,China) (2.College of Electronic Engineering,Naval University of Engineering,Wuhan 430033,China)

In order to fully describe the time-varying characteristic of the sea clutter and its local singularity,the multifractal theory is introduced to the analysis of radar sea clutter.The multifractal parameters of sea clutter are extracted theoretically,analyses the multifractal spectrum,quality index and singularity index of measured sea clutter data are analysed.The results show that the sea clutter has the characteristics of multifractal features and it can well explain its inherent nonlinear characteristics.Compared with the traditional statistical method,the multifractal analysis method is more accurate than the traditional method in describing the generation mechanism of sea clutter,which provides a basis for radar target detection based on multifractal theory.

sea clutter; multifractal theory; data analysis; quality index; singularity index

??處理·

10.16592/j.cnki.1004-7859.2016.10.008

國家自然科學基金資助項目(42415009)

金丹 Email：jindan667788@126.com

2016-07-20

2016-09-22

TN959.1

A

1004-7859(2016)10-0032-04