楊帆 李新霞

摘要：希爾伯特變換廣泛應用于在數(shù)學與信號處理領域，通過希爾伯特變換，能夠從復雜的信號中提取到瞬時參數(shù)——瞬時振幅、瞬時頻率和瞬時相位，使得我們對短信號和復雜信號的瞬時參數(shù)的定義及計算成為可能，能夠?qū)崿F(xiàn)真正意義上的瞬時信號的提取。通過在MATLAB平臺上實現(xiàn)簡單的希爾伯特變換，觀察對比原信號與變換后的信號之間關系，結(jié)合數(shù)學推導，能對希爾伯特變換有清楚認識。

在通信理論中最早引入了復數(shù)道分析方法，利用它可以比較方便地從地震記錄中分離出波的瞬時振幅、瞬時頻率、瞬時相位等參數(shù)。本文在熟悉希爾伯特變換基礎上，基于MATLAB平臺，繪制三瞬曲線圖，有利于掌握該變換及基本理論。

關鍵詞：希爾伯特變換 信號提取 復數(shù)道分析 三瞬曲線

Abstract： The Hilbert transform is widely used in the field of mathematics and signal processing. The Hilbert transform enables the extraction of instantaneous parameters - instantaneous amplitude， instantaneous frequency and instantaneous phase - from complex signals， making it possible to define and calculate the instantaneous parameters of short signals and complex signals， enabling the real sense of instantaneous signal extraction. By implementing a simple Hilbert transform on MATLAB platform， observing and comparing the relationship between the original signal and the transformed signal， combined with mathematical derivation， we can have a clear understanding of Hilbert transform.

The complex channel analysis method was first introduced in the communication theory， using which the instantaneous amplitude， instantaneous frequency， instantaneous phase and other parameters of the wave can be separated from the seismic record relatively easily. In this paper， based on the familiarity with Hilbert transform， we plot the three-instantaneous curve based on MATLAB platform， which is beneficial to master the transform and the basic theory.

Keywords： Hilbert transform， signal extraction， complex channel analysis， three-transient curve

1 希爾伯特變換的數(shù)學原理

把某一實函數(shù)， 的希爾伯特變換定義[1]為

2 在MATLAB平臺上希爾伯特變換的實現(xiàn)

為了掌握希氏變換，我們利用MATLAB平臺，借助環(huán)境自帶的變換函數(shù)，實現(xiàn)了基本的希爾伯特變換[7]。給定一個正弦信號y=sin（2πx+π/6），下面繪制出其經(jīng)過Hilbert變換后的信號曲線。

3 實驗測試結(jié)果與分析

實驗中，我們繪制了原始sin函數(shù)信號，通過簡單的Hilbert變換后，得到了變換后的圖像，很明顯可以看出，經(jīng)過Hilbert變換后圖像翻轉(zhuǎn)，其函數(shù)的相位變化了π/2，實驗圖像十分符合預期分析，具體情況如圖1所示。最后我們將兩個圖像橫坐標對齊，比較在同一個相位下，y值的情況，為了直觀地看到sin圖像變成了cos函數(shù)圖像。

4 復數(shù)道分析及三瞬曲線圖

4.1復數(shù)道分析原理和數(shù)學基礎

由傅里葉分析得知，任何一種振動都可以把它視為無數(shù)個頻率不同的簡諧振動的疊加結(jié)果[9]。質(zhì)內(nèi)傳播時，由于實際介質(zhì)的影響，它的振幅、頻率和相位都將發(fā)生變化。過來說根據(jù)這些參數(shù)的變化也就可能了解實際介質(zhì)性質(zhì)的變化。因此，引入復道分析法在一個振動中提取參數(shù)[10]。

在坐標系中設橫軸 f（t）為實軸，代表復數(shù)道實部，縱軸為虛軸，p（t）代表復數(shù)道虛部， u（t）為復數(shù)地震道。今矢量 u（t）與實軸夾角為θ，則復數(shù)道的實部和虛部p（t）可以分別表示為

4.2 三瞬曲線物理意義及圖像

原始地震剖面包含豐富的地質(zhì)信息。通過不同的數(shù)字變換手段將原始地震剖面中的多種可利用的信息賦予明確概念，并且單獨地提取出來[12，13]。地震信息剖面，揭示出原始地震剖面中不易被發(fā)現(xiàn)的地質(zhì)異?，F(xiàn)象、油氣存在情況，它們的具體的物理意義如下。

4.2.1 瞬時振幅

瞬時振幅也稱瞬時包絡或反射強度，它相當于地震道的包絡函數(shù)。實際上，它是在某一給定時刻對地震道能量作粗略平滑的一個度量[14]。實驗測試了復雜函數(shù)x=（1+0.5* cos （2* pi *5*t））.* cos （2* pi *50*t+A* sin （2* pi *10*t）），繪制了三瞬曲線圖，具體情況如圖3所示。

4.2.2 瞬時頻率

瞬時頻率則是根據(jù)地震道中心頻率所作的逐點抽樣所得，它實際上相當于瞬時相位的時間導數(shù)[15]。實驗上測試結(jié)果如圖4所示。

4.2.3 瞬時相位

瞬時相位是某一瞬間的相位，即所反映的是在某一時刻地震道與其虛地震道的夾角，它是同一時刻子波真實相位的度量。

基于以上分析，前面實驗都是簡單的三角函數(shù)信號模型，我們進一步進行測試，給定任意函數(shù)包絡，以指數(shù)函數(shù)為例，并利用Hilbert求解包絡以及瞬時頻率，并給出對應的Hilbert譜。具體情況如圖8所示。

5 結(jié)論

本文從希爾伯特變換的數(shù)學過程推導和圖像分析解釋，并簡單介紹了當前希爾伯特變換廣泛應用的場景，從理論上說明了其受廣泛應用的原因，鑒于其自身強大的對信號提取能力，希爾伯特變換將會在更多的領域結(jié)合前沿的理論技術，發(fā)揮出更大的潛力。為了掌握基本的希爾伯特變換原理，本文在MATLAB環(huán)境中進行了希氏變換，分析了實驗結(jié)果。最后又著重介紹了基于復數(shù)道分析法的地震波信號過程，結(jié)合三瞬曲線的物理意義，繪制出了三瞬曲線圖，增加了理論和實驗的認識。

參考文獻

[1] Kukker Amit，Sharma Rajneesh. A Genetic Algorithm Assisted Fuzzy Q-Learning epileptic seizure classifier[J]. Computers and Electrical Engineering，2021，92.117-119

[2] Cai Ning，Li Chenxu，Shi Chao. Pricing Discretely Monitored Barrier Options： When Malliavin Calculus Expansions Meet Hilbert Transforms[J]. Journal of Economic Dynamics and Control，2021.62.213-215

[3] Shamlou Alireza，Reza Feyzi Mohammad，Behjat Vahid. Winding deformation classification in a power transformer based on the time-frequency image of frequency response analysis using Hilbert-Huang transform and evidence theory[J]. International Journal of Electrical Power and Energy Systems，2021，129.41-45

[4] Zang L.G.，Ohshima S.，Qu Y.F.，Shi P.W.，Zhong W.L.，Hou Y.M.，Yan L.W.，Ji X.Q.，Li J.Q.，Yu D.L.，Shi Z.B.，Liu Yi，Yang Q.W.，Xu M. Phase tracking with Hilbert transform and nonlinear wave-wave coupling analysis on the HL-2A tokamak[J]. Nuclear Fusion，2021，61（2）.31-40

[5] Yasser H. El-Sharkawy，Sherif Elbasuney. Laser induced fluorescence with 2-D Hilbert transform edge detection algorithm and 3D fluorescence images for white spot early recognition[J]. Spectrochimica Acta Part A： Molecular and Biomolecular Spectroscopy，2020，240（5）.57-59

[6] 李小舟，金海彬.基于希爾伯特變換的信號解調(diào)算法及其在飛機供電特性參數(shù)測試系統(tǒng)中的應用[J].計量學報，2020，41（03）：344-348.

[7] 呂丙南，陳學華，李泂，瞿雷，羅鑫. 基于二維希爾伯特變換的三維地震資料體邊緣檢測[A]. 中國石油學會（CPS）、國際勘探地球物理學家學會（SEG）.CPS/SEG北京2018國際地球物理會議暨展覽電子論文集[C].中國石油學會（CPS）、國際勘探地球物理學家學會（SEG）：石油地球物理勘探編輯部，2018：4.

[8] 梁波. 基于小波包和希爾伯特包絡分析的盾構機主軸承故障診斷方法研究[D].蘭州理工大學，2018.34-41

[9] 宮寧寧. 非線性轉(zhuǎn)子裂紋檢測及延緩裂紋擴展機理的研究[D].天津理工大學，2017.22-25

[10]劉譯膠. 一種時域希爾伯特變換方法[D].重慶大學，2016.45-50

[11] 陳常樂. 地震波場構造導向濾波關鍵技術及應用[D].吉林大學，2015.55-60

[12] 趙航，陶重陽，李高磊.水下目標特征提取中希爾伯特黃變換方法的應用[J].中國新通信，2014，16（07）：88.

[13] 張靜怡.基于希爾伯特變換原理與應用的研究[J].科技與創(chuàng)新，2014（05）：133-134.

[14] 王宏巖.希爾伯特變換在齒輪故障診斷中的應用[J].設備管理與維修，2014（01）：29-30.

[15] 曲鑫，劉財，馮晅，楊慶節(jié)，盧勇旭.利用加權積分能譜方法和復數(shù)道分析技術對薄層的識別與分析[J].世界地質(zhì)，2013，32（01）：114-122.

基金項目：國家自然科學基金項目（11775108）

第一作者：楊帆（1996-），男，漢族，江蘇人，碩士研究生，主要從事核科學技術研究。

* 通訊作者：李新霞（1978-），男，漢族，博士生導師，教授。研究方向：核聚變