楊帆 李新霞

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

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

關(guān)鍵詞：希爾伯特變換 信號(hào)提取 復(fù)數(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ù)學(xué)原理

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

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

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

3 實(shí)驗(yàn)測(cè)試結(jié)果與分析

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

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

4.1復(fù)數(shù)道分析原理和數(shù)學(xué)基礎(chǔ)

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

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

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

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

4.2.1 瞬時(shí)振幅

瞬時(shí)振幅也稱瞬時(shí)包絡(luò)或反射強(qiáng)度，它相當(dāng)于地震道的包絡(luò)函數(shù)。實(shí)際上，它是在某一給定時(shí)刻對(duì)地震道能量作粗略平滑的一個(gè)度量[14]。實(shí)驗(yàn)測(cè)試了復(fù)雜函數(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 瞬時(shí)頻率

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

4.2.3 瞬時(shí)相位

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

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

5 結(jié)論

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

參考文獻(xiàn)

[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] 李小舟，金海彬.基于希爾伯特變換的信號(hào)解調(diào)算法及其在飛機(jī)供電特性參數(shù)測(cè)試系統(tǒng)中的應(yīng)用[J].計(jì)量學(xué)報(bào)，2020，41（03）：344-348.

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

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

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

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

[11] 陳常樂. 地震波場(chǎng)構(gòu)造導(dǎo)向?yàn)V波關(guān)鍵技術(shù)及應(yīng)用[D].吉林大學(xué)，2015.55-60

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

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

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

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

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目（11775108）

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

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