維格納分佈

维基百科,自由的百科全书
跳转至: 导航搜索

維格納分布(又名韋格納分佈,英文: Wigner Distribution Function,縮寫為WDF) 是由1963年的諾貝爾物理學獎得主尤金·维格纳,于1932年首次引用的一個新的方程式。 眾所皆知,傅立葉變換對於研究穩態(時間獨立)的訊號(波形)是一項非常有用的工具,然而,訊號(波形)一般來說在時間上並非是獨立的,這樣的訊號或是波形傅立葉變換並無法有效地完全分析其特性,因此對於一個非穩態的訊號完全分析需要測量出時間以及頻率上的表現。本頁面介紹的數學函數是時頻分析中的基礎方法,在1980年,Claasen,Mecklenbrauker對WDF做了更進一步的研究。除此之外,線性時頻分析中,STFT、Gabor transform和WDF扮演了相當重要的角色,其中WDF對於分析很多非穩態的隨機訊號都有很好的表現,例如:量子力學光學聲學通訊生物工程訊號處理影像處理。有時也被用在分析地震的資料,以及處理聲音的相位失真。

定義[编辑]

定義一[编辑]

.....(1)

定義二[编辑]

.....(2)

定義二與定義一之間的關係 :

其他定義[编辑]

.....(3)

聲納雷達系統中,傳送出去的聲波的反射波可以用來偵測目標物的位置跟速度,在很多情形下,收到的訊號因為都普勒位移,所以跟原本的訊號並不一樣。Woodward(1953) 改寫了原本的公式

這個公式被稱為Woodward ambiguity function,這個式子在雷達系統的訊號處理和設計上扮演重要的角色。

WDF和STFT的比較[编辑]

WDF、STFT和Gabor transform 都佔了時頻分析中非常重要的地位,在這邊比較一下它們之間的差別。

WDF STFT
清晰度 較好 較差
相交項的問題 嚴重
複雜度

相交項其實就是處理的過程中產生的額外訊號,是不想要的,WDF的清晰度和複雜度是彼此做取捨的,可以依不同的情況或是不同的方法來決定是否要使用WDF或是另外兩種。

WDF的優缺點[编辑]

在這裡列出WDF主要的優缺點

優點 :

1.有良好的解析度,尤其是對單一成分,且瞬時頻率變化不為2次式以上。

2.有好的數學運算性質(見WDF的數學性質)。

3.可用於分析隨機程序。

缺點 :

1.有相交項(cross term)的問題,改進方法請見 改進型韋格納分佈

2.需要更多的時間去計算,若訊號時間越長,則需要更久的時間。

3.不是一對一函數,無法辨別相位部分,例如:

4.不適合分析瞬時頻率變化為2次式以上的型態,即

相交項特性[编辑]

WDF不是一個線性的轉換,如果有兩種以上不同性質的訊號疊加,會產生相交項。 然而對於物理,相交項是重要且有用的,可以用來分析期望值,詳見維格納分佈。以下數學方程式對於WDF後會產生相交項。



常見的時頻分析例子[编辑]

以下的例子說明如何用WDF來做時頻分析

常數訊號[编辑]

輸入訊號為常數,則時頻分佈為一條線重合於時軸,如果'x(t) = 1,則:

弦波訊號[编辑]

輸入訊號為弦波,則時頻分佈為一條線平行於時軸,如果x(t) = e i2πkt,則:

啁啾聲信號[编辑]

啁啾聲訊號的瞬時頻率隨時間線性,表示時頻分佈為一條斜值線,例如

,

則瞬時頻率為:

故WDF為:

單位脈衝訊號[编辑]

因為單位脈衝包含所有的頻率分佈,且在時間不等於零時沒值,故WDF為通過原點的且與時軸垂直的線


方波訊號[编辑]

,
.


WDF的數學性質[编辑]

(1)投射特性 ,
(2)能量特性
(3)回覆特性 ,
(4)Moment特性 ,
(5)是實數
(6)區域特性 If for then for , If for then for
(7)乘法特性 If ,then
(8)摺積特性 If ,then
(9)相關特性 If ,then
(10)時間平移特性 If , then
(11)調變特性 If , then

WDF實現方法[编辑]

以下為電腦計算WDF的實現方式

(i) 直接運算(暴力法)

複雜度 : TF(2Q+1)

(ii) 使用離散傅立葉變換

複雜度 :

(iii) 使用Chirp-Z 轉換

複雜度 : ,通常為使用離散傅立葉變換的2~3倍,但限制比使用離散傅立葉變換

在使用這三個方法前,先來做個前提討論

從定義一出發

再令 ,則上述式子則為

下面介紹的三種方法都是從這條式子開始推導

注意事項 :

若x(t)是無限長的訊號,則p要從負無限加到正無限,這點不易實現。

若x(t)為有限長的訊號,則p範圍可以縮小,就可能實現。


故下面三種方法都是在第2種情況下討論,即x(t)為有限長訊號,p範圍可以縮小

我們假設


直接運算(暴力法)[编辑]


限制條件 :

只有一個 : 要滿足Nyquist criterion

,其中B是的頻寬,大約是x(t)的兩倍。

推導 :

所以當 時,

固定中間的n值 () 來探討p的範圍

-– (1)

-- (2)

其中 (1) & (2) 的下限是同義的

故(1) & (2)皆可改寫為

且可以發現 代表 離兩個邊界的距離

注意事項: 當 n > n2 或 n < n1 時,將沒有 p 能滿足上面的不等式

最後推導出的式子如下

其中

使用離散傅立葉變換[编辑]


限制條件 :

(1)要滿足Nyquist criterion

,其中B是的頻寬,大約是x(t)的兩倍。

(2)

(3)

推導 :

前提討論的式子可以改寫為

針對中間


最後得出的式子如下

其中

使用Chirp-Z 轉換[编辑]


限制條件 :

只有一個 : 要滿足Nyquist criterion

,其中B是的頻寬,大約是x(t)的兩倍。

推導 :

前提討論的式子可改寫為

計算分成3步驟

STEP 1 :

STEP 2 :

其中

STEP 3 :

延伸變化[编辑]

視窗型韋格納分佈[编辑]

視窗型韋格納分佈Windowed Wigner Distribution Function),在韋格納分佈中,當x(t)為無限長訊號時,WDF很難去實現它。所以在積分中加入一個新的函數 ,目的是擷取x(t)中的片段來計算,不需從負無限積分到正無限。

定義

其中為實數且為有限長訊號

原始韋格納分佈定義

優缺點

優點 : (1) 降低運算時間,因為為有限長函數。 (2) 可以有效降低相交項(cross term)問題,但不能完全消除(詳見下方說明)。

缺點 : (1) 一些相交項(cross term)問題仍被保留。 (2) 可能不符合譜密度(Power spectral density)的定義。 (3) 一些好用的數學運算性質會消失。

實現方法

從定義出發 再令 假設w(t) = 0 for |t| > B 即 其中 如此一來,p範圍便可縮小。

避免相交項的原因

從定義出發 其中為實數且為有限長訊號 假設的情況下 比較有無mask function的不同

x(t)示意圖

理想情形 :

沒有使用mask function

即mask function

Wx(t,f)示意圖

總共有3種情況要討論

ideal x(t)。Auto term 為自相關項。 Cross term為相交項
使用mask function

假設 如果

ideal x(t)。

由於只在-B到B有值,故乘上就能去除相交項(Cross term),而保留上圖中兩條紅線中間的區域,也就是Auto term。

但上述其實是理想的情況,x(t)為Delta function是窄頻信號。 如果X(τ)寬度太寬或是有ripple的話,Cross term仍會有殘留,示意圖如下

non ideal x(t)。

藍色線為X(τ)的訊號,若X(τ)的寬度太寬或是有ripple產生,就可能會跑進的範圍裡面,進而導致無法完全濾除Cross term。

改進型韋格納分佈[编辑]

此方法可以消除相交項(cross term)。