維格納分佈

維格納分佈(又名韋格納分佈，英文: 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_9...(1)formula_10...(2)定義二.定義二與定義一之間的關係 : formula_11formula_12...(3)其他定義.在聲納和雷達系統中，傳送出去的聲波的反射波可以用來偵測目標物的位置跟速度，在很多情形下，收到的訊號因為都普勒位移，所以跟原本的訊號並不一樣。Woodward(1953) 改寫了原本的公式formula_13這個公式被稱為Woodward ambiguity function，這個式子在雷達系統的訊號處理和設計上扮演重要的角色。而維格納分佈亦為科恩系列分佈的其中一種特例，當科恩系列分佈中的formula_14時，科恩系列分佈會是維格納分佈。WDF和STFT的比較.WDF、STFT和Gabor transform 都佔了時頻分析中非常重要的地位，在這邊比較一下它們之間的差別。相交項其實就是處理的過程中產生的額外訊號，是不想要的，WDF的清晰度和複雜度是彼此做取捨的，可以依不同的情況或是不同的方法來決定是否要使用WDF或是另外兩種。WDF的優缺點.在這裡列出WDF主要的優缺點優點 :1.有良好的解析度，尤其是對單一成分，且瞬時頻率變化不為2次式以上。2.有好的數學運算性質(見WDF的數學性質)。3.可用於分析隨機程式(見WDF與隨機程式的關係)。缺點 :1.有相交項(cross term)的問題，改進方法請見 改進型韋格納分佈。2.需要更多的時間去計算，若訊號時間越長，則需要更久的時間。3.不是一對一函式，無法辨別相位部分，例如: formula_154.不適合分析瞬時頻率變化為2次式以上的型態，即formula_16。相交項特性.WDF不是一個線性的轉換，由於formula_17的signal auto-correlation function formula_18，如果有兩種以上不同性質的訊號疊加，會產生相交項。然而相交項卻有重要且有用的物理意義，像是可以用來分析期望值，相反的，短時距傅立葉轉換就沒有此特性，詳見維格納分佈與隨機程式的關係。以下數學方程式對於WDF後會產生相交項。WDF與隨機程式的關係.對於一個隨機程式x(t)，我們無法得知其確切的值，因此會將其值表示為一個機率函式，通常E[x(t)] = 0 for any t將x(t)的維格納分佈取期望值後可得其譜密度(Power spectral density，PSD)，如下公式所列formula_21formula_22當x(t)的統計特性不隨時間變化時，可稱x(t)為平穩的隨機程式，其譜密度也可簡化為formula_23，也就是說維格納函式能初略的告訴我們譜密度如何隨時間進行變化。維格納函式能在平穩程式對所有時間ｔ都簡化成譜密度，然而也等同於非平穩的自相關函式，這也是維格納分佈的動機。下圖為一個平穩的隨機程式進行維格納分析後的例子，可明顯看出此訊號不隨時間變化，也就是時頻分析結果為水平線。反之，亦可利用時頻分析結果是否為水平線判斷該訊號是否為一平穩的訊號。而在訊號處理中常見的白雜訊，其譜密度formula_24，其中formula_25為一個常數。白雜訊的維格納分佈如下圖，可看出此雜訊在所有時間及頻率都存在著。維格納分佈的相交項在處理隨機程式時派上用場，相對的，沒有相交項的短時距傅立葉轉換，則無法用於隨機程式，如下公式所示，只有在零平均隨機程式時，formula_26formula_27常見的時頻分析例子.以下的例子說明如何用WDF來做時頻分析常數訊號.輸入訊號為常數，則時頻分佈為一條線重合於時軸，如果'x"("t") = 1，則：formula_28弦波訊號.輸入訊號為弦波，則時頻分佈為一條線平行於時軸，如果"x"("t")e i2π"kt"，則：formula_29啁啾聲訊號.啁啾聲訊號的瞬時頻率隨時間線性，表示時頻分佈為一條斜值線，例如formula_30 ,則瞬時頻率為：formula_31故WDF為：formula_32餘弦訊號.x(t) = cos(440formula_33t), 當 t 小於 0.5， 頻率 f = 220Hzx(t) = cos(660formula_33t), 當 0.5 小於等於 t 小於 1， 頻率 f = 330Hzx(t) = cos(524formula_33t), 當 t 大於等於 1， 頻率 f = 262Hz單位脈衝訊號.因為單位脈衝包含所有的頻率分佈，且在時間不等於零時沒值，故WDF為透過原點的且與時軸垂直的線formula_36formula_37 ,formula_38 .WDF實現方法.以下為電腦計算WDF的實現方式在使用這三個方法前，先來做個前提討論從定義一出發formula_42令 formula_43formula_44再令 formula_45 ，則上述式子則為formula_46下面介紹的三種方法都是從這條式子開始推導注意事項 :若x(t)是無限長的訊號，則p要從負無限加到正無限，這點不易實現。若x(t)為有限長的訊號，則p範圍可以縮小，就可能實現。故下面三種方法都是在第2種情況下討論，即x(t)為有限長訊號，p範圍可以縮小我們假設formula_47直接運算(暴力法).限制條件 :只有一個 : 要滿足Nyquist criterionformula_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) &amp; (2) 的下限是同義的故(1) &amp; (2)皆可改寫為formula_57且可以發現 formula_59 代表 formula_53 離兩個邊界的距離注意事項： 當 n &gt; n2 或 n  B即 formula_92其中formula_93如此一來，p範圍便可縮小。formula_94避免相交項的原因.從定義出發 formula_87，其中formula_84為實數且為有限長訊號假設formula_97的情況下，比較有無mask function所產生的不同結果理想情形 : formula_98沒有使用mask function.即mask function formula_99formula_100總共有3種情況要討論，如下圖，可見cross term在沒有使用mask function時，無法被消除使用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)問題，在某些情況下比加伯轉換擁有更好的清晰度。