維格納分佈
維格納分布(又名韋格納分佈，英文: Wigner Distribution Function,縮寫為WDF) 是由1963年的諾貝爾物理學獎得主尤金·维格纳，于1932年首次引用的一個新的方程式。
眾所皆知，傅立葉變換對於研究穩態(時間獨立)的訊號(波形)是一項非常有用的工具，然而，訊號(波形)一般來說在時間上並非是獨立的，這樣的訊號或是波形傅立葉變換並無法有效地完全分析其特性，因此對於一個非穩態的訊號完全分析需要測量出時間以及頻率上的表現。本頁面介紹的數學函數是時頻分析中的基礎方法，在1980年，Claasen，Mecklenbrauker對WDF做了更進一步的研究。除此之外，線性時頻分析中，STFT、Gabor transform和WDF扮演了相當重要的角色，其中WDF對於分析很多非穩態的隨機訊號都有很好的表現，例如：量子力學、光學、聲學、通訊、生物工程、訊號處理和影像處理。有時也被用在分析地震的資料，以及處理聲音的相位失真。
維格納分布有許多不同的定義，而此處的定義是特別針對時頻分析而定的。若給定一時間序列formula_1，它的非平穩自相關函數如下公式所列
formula_2
其中formula_3代表所有可能實驗的程序的平均，formula_4代表平均，其可能是時間的函數也有可能不是。維格納函數formula_5起初是以包含時間平均formula_6與時間差formula_7的自相關函數和時間差進行傅立葉轉換來表示，如下：
formula_8
對於單一零平均的時間序列，維格納函數可以簡化如下：
定義二與定義一之間的關係 : formula_11
在聲納和雷達系統中，傳送出去的聲波的反射波可以用來偵測目標物的位置跟速度，在很多情形下，收到的訊號因為都普勒位移，所以跟原本的訊號並不一樣。Woodward(1953) 改寫了原本的公式
這個公式被稱為Woodward ambiguity function，這個式子在雷達系統的訊號處理和設計上扮演重要的角色。
而維格納分布亦為科恩系列分布的其中一種特例，當科恩系列分布中的formula_14時，科恩系列分布會是維格納分布。
WDF、STFT和Gabor transform 都佔了時頻分析中非常重要的地位，在這邊比較一下它們之間的差別。
相交項其實就是處理的過程中產生的額外訊號，是不想要的，WDF的清晰度和複雜度是彼此做取捨的，可以依不同的情況或是不同的方法來決定是否要使用WDF或是另外兩種。
在這裡列出WDF主要的優缺點
優點 :
1.有良好的解析度，尤其是對單一成分，且瞬時頻率變化不為2次式以上。
2.有好的數學運算性質(見WDF的數學性質)。
3.可用於分析隨機程序(見WDF與隨機程序的關係)。
缺點 :
1.有相交項(cross term)的問題，改進方法請見 改進型韋格納分佈。
2.需要更多的時間去計算，若訊號時間越長，則需要更久的時間。
3.不是一對一函數，無法辨別相位部分，例如: formula_15
4.不適合分析瞬時頻率變化為2次式以上的型態，即formula_16。
WDF不是一個線性的轉換，由於formula_17的signal auto-correlation function formula_18，如果有兩種以上不同性質的訊號疊加，會產生相交項。然而相交項卻有重要且有用的物理意義，像是可以用來分析期望值，相反的，短時距傅立葉轉換就沒有此特性，詳見維格納分佈與隨機程序的關係。以下數學方程式對於WDF後會產生相交項。 對於一個隨機程序x(t)，我們無法得知其確切的值，因此會將其值表示為一個機率函數，通常E[x(t)] = 0 for any t
將x(t)的維格納分布取期望值後可得其譜密度(Power spectral density，PSD)，如下公式所列
formula_21
formula_22
當x(t)的統計特性不隨時間變化時，可稱x(t)為平穩的隨機程序，其譜密度也可簡化為formula_23，也就是說維格納函數能初略的告訴我們譜密度如何隨時間進行變化。維格納函數能在平穩程序對所有時間ｔ都簡化成譜密度，然而也等同於非平穩的自相關函數，這也是維格納分布的動機。
下圖為一個平穩的隨機程序進行維格納分析後的例子，可明顯看出此信號不隨時間變化，也就是時頻分析結果為水平線。反之，亦可利用時頻分析結果是否為水平線判斷該訊號是否為一平穩的訊號。
而在訊號處理中常見的白雜訊，其譜密度formula_24，其中formula_25為一個常數。白雜訊的維格納分布如下圖，可看出此雜訊在所有時間及頻率都存在著。 維格納分布的相交項在處理隨機程序時派上用場，相對的，沒有相交項的短時距傅立葉轉換，則無法用於隨機程序，如下公式所示，只有在零平均隨機程序時，formula_26
formula_27
以下的例子說明如何用WDF來做時頻分析
輸入訊號為常數，則時頻分佈為一條線重合於時軸，如果'x"("t") = 1，則：
輸入訊號為弦波，則時頻分佈為一條線平行於時軸，如果，則：
啁啾聲訊號的瞬時頻率隨時間線性，表示時頻分佈為一條斜值線，例如
則瞬時頻率為：
故WDF為： x(t) = cos(440formula_33t), 当 t 小於 0.5， 頻率 f = 220Hz x(t) = cos(660formula_33t), 当 0.5 小於等於 t 小於 1， 頻率 f = 330Hz x(t) = cos(524formula_33t), 当 t 大於等於 1， 頻率 f = 262Hz
因為單位脈衝包含所有的頻率分佈，且在時間不等於零時沒值，故WDF為通過原點的且與時軸垂直的線
以下為電腦計算WDF的實現方式在使用這三個方法前，先來做個前提討論
從定義一出發
formula_42
令 formula_43
formula_44
再令 formula_45 ，則上述式子則為
formula_46
下面介紹的三種方法都是從這條式子開始推導
注意事項 : 若x(t)是無限長的訊號，則p要從負無限加到正無限，這點不易實現。
若x(t)為有限長的訊號，則p範圍可以縮小，就可能實現。 故下面三種方法都是在第2種情況下討論，即x(t)為有限長訊號，p範圍可以縮小
我們假設formula_47 限制條件 :
只有一個 : 要滿足Nyquist criterion
formula_48，其中B是formula_49的頻寬，大約是x(t)的兩倍。
推導 :
formula_47
所以當formula_51 時，
formula_52
固定中間的n值 (formula_53) 來探討p的範圍
formula_54
即 formula_55-– (1)
formula_56
即 formula_57-- (2)
其中 (1) & (2) 的下限是同義的
故(1) & (2)皆可改寫為
formula_57
且可以發現 formula_59 代表 formula_53 離兩個邊界的距離
注意事項： 當 n > n2 或 n < n1 時，將沒有 p 能滿足上面的不等式
最後推導出的式子如下
formula_61
其中 formula_62限制條件 :
(1)要滿足Nyquist criterion
formula_48，其中B是formula_49的頻寬，大約是x(t)的兩倍。
(2)formula_65
(3) formula_66
推導 :
前提討論的式子可以改寫為
formula_67
令formula_68
formula_69
針對中間formula_70項
令
formula_71
formula_72 最後得出的式子如下
formula_73
其中
formula_74
formula_71
formula_72限制條件 :
只有一個 : 要滿足Nyquist criterion
formula_48，其中B是formula_49的頻寬，大約是x(t)的兩倍。
推導 :
前提討論的式子可改寫為
formula_79
計算分成3步驟
STEP 1 : formula_80
STEP 2 : formula_81
STEP 3 : formula_82
視窗型韋格納分佈（），在韋格納分佈中，當x(t)為無限長訊號時，WDF很難去實現它。所以在積分中加入一個新的函數 ，目的是擷取x(t)中的片段來計算，不需從負無限積分到正無限。
formula_83 ， 其中formula_84為實數且為有限長訊號
原始韋格納分佈定義 formula_42從定義出發  formula_87  令formula_88  formula_89  再令formula_45  formula_91  假設w(t) = 0 for |t| > B  即 formula_92  其中formula_93  如此一來，p範圍便可縮小。  formula_94
從定義出發 formula_87，其中formula_84為實數且為有限長訊號
假設formula_97的情況下，比較有無mask function所產生的不同結果
理想情形 : formula_98
即mask function formula_99  formula_100  總共有3種情況要討論，如下圖，可見cross term在沒有使用mask function時，無法被消除
formula_101
假設formula_102，且formula_103
由於formula_104只在-B到B有值，故乘上formula_104就能去除相交項(Cross term)，只保留下圖中兩條紅線中間的區域，也就是Auto terms。
但上述其實是理想的情況，x(t)為窄頻信號Delta function  如果X(τ)寬度太寬或是有ripple的話，Cross term仍會有殘留，示意圖如下
藍色線為X(τ)的訊號，若X(τ)的寬度太寬或是有ripple產生，就可能會跑進formula_104的範圍裡面，進而導致無法完全濾除Cross term。
cross term 只有在訊號每個成分的寬度都小於2B，且時間差formula_107都大於B時，才能被消除
此方法可以消除相交項(cross term)。
消除相交項(cross term)問題，在某些情況下比加伯轉換擁有更好的清晰度。