瞬時頻率
瞬時頻率
在信號處理中,觀察信號的瞬時頻率是很重要的課題。假設一信號可寫成指數信號的N項相加(有無線多種表示法,以小的為宜),即,其中為常數。則瞬時頻率(以頻率表示), k=1,...,N瞬時頻率 (Hz)以頻率來表示(單位為赫茲):, k=1,...,N。
以角速度來表示(單位為弧度每秒): , k=1,...,N
直觀上,瞬時頻率為相位的微分。對於自然界中的實數訊號,如何定義訊號的相位。Gabor提出解析訊號法(Analytic Signal Method),將實數訊號表示為對應的複數訊號,即可定義複數訊號的大小與相位,將實數訊號的瞬時頻率求出。實數訊號的解析訊號(Analytic Signal)定義為
解析函數的極坐標表示
瞬時相位
瞬時頻率
其中虛數項為實數訊號的希爾伯特轉換(Hilbert Transform),將它定義為。稱作解析函數的理由是,此型式的複數函數滿足柯西-里曼(Cauchy-Riemann)的可微分條件,稱之為解析函數(Analytic Function)。因此,解析訊號可以表示為
其中
; 如果是沒有任何限制條件的時間訊號,計算出來的瞬時頻率可能不是正確的結果。對於平均值為零的局部對稱訊號而言,前述定義的瞬時頻率才具有物理意義。在1998年,黃鍔(Norden E. Huang)博士提出一個有效的演演演算法,將訊號先行分解成具有局部對稱之分量,以正確地求得資料的瞬時頻率。這個方法稱為希爾伯特-黃轉換(Hilbert Huang Transform, HHT)。
以下簡單的例子來說明,對於平均值為零的訊號,此瞬時頻率的定義才具有物理意義。對於一個弦波訊號,
考慮三種情況: (1) β = 0 (2) 0 < β < 1 (3) β > 1
(1) β = 0: 當弦波訊號平均值為零時,在複數平面上的描述是以座標原點為中心的單位圓,它的相位角θ(t)則是以座標原點為中心,逆時鐘方向呈線性遞增,其圖形為斜率1的直線,而瞬時頻率是一個常數值。
(2) 0 < β < 1: 在複數平面上仍然是一個單位圓,但其圓心從原點偏移了β個單位,其相角θ(t)不再呈現線性遞增,瞬時頻率出現震蕩的現象,而不是應有的常數值。
(3) β > 1: 因為β值超過單位圓的半徑,因此的圓心在單位圓之外。如此相位角θ(t)在[ − π/2, π/2]震蕩,瞬時頻率出現負值,與原訊號的特性有極大的差別。
因為在目前許多數位信號處理的應用上都與信號的頻譜或信號的頻寬有很大的關係。
若能確實地偵測信號的瞬時頻率,則通道頻寬可以被可適性(adaptive)的決定,如此一來能更有效地利用系統資源,提高系統效能。
相關應用調變(modulation) 多工方式(multiplexing) 濾波器的設計(filter design) 信號壓縮(data compression) 信號分析(signal analysis) 信號辨識(signal identification) 語音信號處理(acoustical signal processing) 製作系統的模型(system modling) 雷達系統的分析(rader system analysis) 取樣(sampling)
瞬時頻率為常數-使用傅立葉轉換
當瞬時頻率為常數即 為一階時間函數,使用傅立葉轉換做信號分析。
由於從傅立葉轉換中是無法觀察出信號頻譜隨著時間改變的變化。
故只有當瞬時頻率為常數,不是時間的函數時,便可使用傅立葉轉換做信號分析。
瞬時頻率不為常數-使用時頻分析
當瞬時頻率不為常數即 為高階時間函數,使用時頻分析做信號分析。
從時頻分析可觀察出信號頻率隨著時間變化的改,這是傅立葉轉換無法做到的。
因此當瞬時頻率為時間的函數,使用時頻分析做信號分析,如下圖,可以確切地觀察到信號瞬時頻率的變化。