查看“維格納分佈”的源代码

'''維格納分布'''(又名'''韋格納分佈'''，英文: '''Wigner Distribution Function''',縮寫為'''WDF''') 是由1963年的[[諾貝爾物理學獎]]得主[[尤金·维格纳]]，于1932年首次引用的一個新的[[方程式]]。

眾所皆知，[[傅立葉變換]]對於研究穩態(時間獨立)的訊號(波形)是一項非常有用的工具，然而，訊號(波形)一般來說在時間上並非是獨立的，這樣的訊號或是波形傅立葉變換並無法有效地完全分析其特性，因此對於一個非穩態的訊號完全分析需要測量出時間以及頻率上的表現。本頁面介紹的數學函數是[[時頻分析]]中的基礎方法，在1980年，Claasen，Mecklenbrauker對WDF做了更進一步的研究。除此之外，線性時頻分析中，STFT、Gabor transform和WDF扮演了相當重要的角色，其中WDF對於分析很多非穩態的隨機訊號都有很好的表現，例如：[[量子力學]]、[[光學]]、[[聲學]]、[[通訊]]、[[生物工程]]、[[訊號處理]]和[[影像處理]]。有時也被用在分析地震的資料，以及處理聲音的相位失真。

==定義==
維格納分布有許多不同的定義，而此處的定義是特別針對時頻分析而定的。若給定一時間序列<math>x[t]</math>，它的非平穩[[自相关函数|自相關函數]]如下公式所列

　　　<math> C_x(t_1, t_2) = \left\langle \left(x[t_1] - \mu[t_1]\right) \left(x[t_2] - \mu[t_2]\right)^* \right\rangle ,</math>

其中<math>\langle \cdots \rangle</math>代表所有可能實驗的程序的平均，<math>\mu[t]</math>代表平均，其可能是時間的函數也有可能不是。維格納函數'''<math>W_x(t,f)</math>'''起初是以包含時間平均<math>t = (t_1+t_2)/2</math>與時間差<math>\tau = t_1 - t_2</math>的自相關函數和時間差進行傅立葉轉換來表示，如下：

　　<math> W_x(t,f)=\int_{-\infty}^{\infty} C_x\left(t + \frac{\tau}{2}, t - \frac{\tau}{2}\right) \, e^{-2\pi i\tau f} \, d\tau .</math>

對於單一零平均的時間序列，維格納函數可以簡化如下：

===定義一===
:<math>W_x(t,f)= \int_{-\infty}^{\infty} x(t+\frac{\tau}{2})x^*( t-\frac{\tau}{2})e^{-j2\pi\tau f}\, d\tau </math>.....(1)

===定義二===
:<math>W_x(t,\omega)= \int_{-\infty}^{\infty} x(t+\frac{\tau}{2})x^*( t-\frac{\tau}{2})e^{-j\omega\tau}\, d\tau </math>.....(2)

定義二與定義一之間的關係 : <math>\omega=2\pi f</math>
===其他定義===
:<math>\sqrt{\frac{1}{2\pi}}W_x(t,\omega)= \int_{-\infty}^{\infty} x(t+\frac{\tau}{2})x^*( t-\frac{\tau}{2})e^{-j\omega\tau}\, d\tau </math>.....(3)

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

:<math>A_x(t,\omega)= \int_{-\infty}^{\infty} x(\tau +\frac{t}{2})x^*( \tau-\frac{t}{2})e^{-j\omega\tau}\, d\tau </math>

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

而維格納分布亦為[[科恩系列分佈|科恩系列分布]]的其中一種特例，當科恩系列分布中的<math>\Phi(\eta,\tau)=1  </math>時，科恩系列分布會是維格納分布。

==WDF和STFT的比較==
WDF、STFT和Gabor transform 都佔了時頻分析中非常重要的地位，在這邊比較一下它們之間的差別。
{| class="wikitable" 
! !!WDF!!STFT
|-
|清晰度||較好
|較差
|-
|相交項的問題||嚴重
|無
|-
|複雜度||高
|低
|-
|處理隨機程序
|可
|不可
|}
相交項其實就是處理的過程中產生的額外訊號，是不想要的，WDF的清晰度和複雜度是彼此做取捨的，可以依不同的情況或是不同的方法來決定是否要使用WDF或是另外兩種。

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

優點 : 

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

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

3.可用於分析隨機程序(見WDF與隨機程序的關係)。

缺點 : 

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

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

3.不是一對一函數，無法辨別相位部分，例如: <math>WDF[x(t)] = WDF[x(t){e^{j\phi }}]</math>

4.不適合分析瞬時頻率變化為2次式以上的型態，即<math>{e^{j{t^n}}},n \ne 0,1,2</math>。

==相交項特性==
WDF不是一個線性的轉換，由於<math>x(t) </math>的signal auto-correlation function  <math>x(t+\frac{\tau}{2})x^*( t-\frac{\tau}{2}) </math>，如果有兩種以上不同性質的訊號疊加，會產生相交項。然而相交項卻有重要且有用的物理意義，像是可以用來分析期望值，相反的，短時距傅立葉轉換就沒有此特性，詳見維格納分佈與隨機程序的關係。以下數學方程式對於WDF後會產生相交項。<br />

*<math>x(t)=\begin{cases} \cos(2\pi t) & t\le-2 \\ \cos(4\pi t) & -2 < t \le 2 \\ \cos(3\pi t) & t>2 \end{cases}</math>

<!-- Image with unknown copyright status removed: [[File:WDF4_080115.jpg|none|Gabor vs WDF|{{Deletable image-caption|1=Sunday, 13 April 2008|date=May 2012}}]] -->

*<math>x(t)=e^{it^3}</math>

<!-- Image with unknown copyright status removed: [[File:WDF5_080115.jpg|none|Gabor vs WDF|{{Deletable image-caption|1=Sunday, 13 April 2008|date=May 2012}}]] -->== WDF與隨機程序的關係 ==
對於一個隨機程序x(t)，我們無法得知其確切的值，因此會將其值表示為一個機率函數，通常E[x(t)] = 0 for any t

將x(t)的維格納分布取期望值後可得其譜密度(Power spectral density，PSD)，如下公式所列

<math>E[W_x(t,f)]=\textstyle \int\limits_{ -\infty }^{\infty} \displaystyle E[x(t+\tau/2)x^*(t- \tau/2)]\cdot e^{-j2 \pi f \tau}\cdot d \tau</math>

<math>=\textstyle \int\limits_{ -\infty }^{\infty} \displaystyle R_x(t, \tau)\cdot e^{-j2 \pi f \tau}\cdot d \tau=S_x(t,f)</math>

當x(t)的統計特性不隨時間變化時，可稱x(t)為平穩的隨機程序，其譜密度也可簡化為<math>E[W_x(t,f)]=S_x(f)</math>，也就是說維格納函數能初略的告訴我們譜密度如何隨時間進行變化。維格納函數能在平穩程序對所有時間ｔ都簡化成譜密度，然而也等同於非平穩的自相關函數，這也是維格納分布的動機。

下圖為一個平穩的隨機程序進行維格納分析後的例子，可明顯看出此信號不隨時間變化，也就是時頻分析結果為水平線。反之，亦可利用時頻分析結果是否為水平線判斷該訊號是否為一平穩的訊號。

[[File:Stationary_random_process.png|200x200像素]]

而在訊號處理中常見的[[白雜訊]]，其譜密度<math>S_x(f)=\sigma</math>，其中<math>\sigma</math>為一個常數。白雜訊的維格納分布如下圖，可看出此雜訊在所有時間及頻率都存在著。
[[File:White_noise_wdf.png|200x200像素]]

維格納分布的相交項在處理隨機程序時派上用場，相對的，沒有相交項的短時距傅立葉轉換，則無法用於隨機程序，如下公式所示，只有在零平均隨機程序時，<math>E[X(t,f)]=0</math>

<math>E[X(t,f)]=E[\textstyle \int\limits_{t-B}^{t+B} \displaystyle x(\tau)w(t- \tau)e^{-j2 \pi f \tau }d \tau]=\textstyle \int\limits_{t-B}^{t+B} \displaystyle E[x(\tau)]w(t- \tau)e^{-j2 \pi f \tau }d \tau</math>

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

===常數訊號===
輸入訊號為常數，則時頻分佈為一條線重合於時軸，如果'x''(''t'')&nbsp;=&nbsp;1，則：
:<math>W_x(t,f)=\int_{-\infty}^\infty e^{-i2\pi\tau\,f}\,d\tau=\delta(f).</math>

===弦波訊號===
輸入訊號為弦波，則時頻分佈為一條線平行於時軸，如果{{math|''x''(''t'') {{=}} e <sup>i2π''kt''</sup>}}，則：
:<math>\begin{align}
  W_x(t,f) &= \int_{-\infty}^{\infty}e^{ i2\pi k    \left(t + \frac{\tau}{2}\right)}e^{-i2\pi k\left(t - \frac{\tau}{2}\right)}e^{-i2\pi\tau\,f}\,d\tau \\
           &= \int_{-\infty}^{\infty}e^{-i2\pi \tau \left(f - k\right)}\,d\tau\\
           &= \delta(f - k).
\end{align}</math>

===啁啾聲信號===
啁啾聲訊號的瞬時頻率隨時間線性，表示時頻分佈為一條斜值線，例如
:<math>x(t) = e^{i2\pi kt^2}</math> , 

則瞬時頻率為：
:<math>\frac{1}{2\pi}\frac{d(2\pi kt^2)}{dt} = 2kt~,</math>

故WDF為：
:<math>\begin{align}
  W_x(t,f) &= \int_{-\infty}^\infty e^{i2\pi k\left(t + \frac{\tau}{2}\right)^2}e^{-i2\pi k\left(t - \frac{\tau}{2}\right)^2}e^{-i2\pi\tau\,f} \, d\tau \\
           &= \int_{-\infty}^\infty e^{i4\pi kt\tau}e^{-i2\pi\tau f}\,d\tau \\
           &= \int_{-\infty}^\infty e^{-i2\pi\tau(f - 2kt)}\,d\tau\\
           &= \delta(f - 2kt) ~.
\end{align} </math>

===單位脈衝訊號===
因為單位脈衝包含所有的頻率分佈，且在時間不等於零時沒值，故WDF為通過原點的且與時軸垂直的線
:<math>\begin{align}
  W_x(t,f) &= \int_{-\infty}^{\infty}\delta\left(t + \frac{\tau}{2}\right)\delta\left(t - \frac{\tau}{2}\right) e^{-i2\pi\tau\,f}\,d\tau \\
           &= 4\int_{-\infty}^{\infty}\delta(2t + \tau)\delta(2t - \tau)e^{-i2\pi\tau f}\,d\tau \\
           &= 4\delta(4t)e^{i4\pi tf} \\
           &= \delta(t)e^{i4\pi tf} \\
           &= \delta(t).
\end{align}
</math>


===方波訊號===
:<math>x(t) = \begin{cases} 1 & |t|<1/2 \\ 0 & \text{otherwise} \end{cases} \qquad</math> ,

:<math>W_x(t,f) = \frac{1}{\pi f}\sin (f[1 - 2|t|])</math> .

==WDF的數學性質==
{| class="wikitable" 
|-
|(1)投射特性||<math>|x(t)|^2= \int_{-\infty}^{\infty} W_x(t,f)\, df </math>, <math>|X(f)|^2= \int_{-\infty}^{\infty} W_x(t,f)\, dt </math> 
|-
|(2)能量特性||<math>|x(t)|^2= \iint_{-\infty}^{\infty} W_x(t,f)\, dt \, df =\int_{-\infty}^{\infty} |x(t)|^2 \, dt=\int_{-\infty}^{\infty} |X(f)|^2 \, df</math> 
|-
|(3)回覆特性||<math>\int_{-\infty}^{\infty}W_x(\frac{t}{2},f)e^{j2\pi ft}\, df=x(t)\bullet x^*(0) </math>, <math>\int_{-\infty}^{\infty}W_x(t,\frac{f}{2})e^{-j2\pi ft}\, dt=X(f)\bullet X^*(0) </math>
|-
|(4)Mean 條件
|If <math>x(t)=\left \vert x(t) \right \vert e^{j2\pi\phi(t)}, 
\ X(f)=\left \vert X(f) \right \vert e^{j2\pi\Psi(f)}  </math>
then <math>\phi'(t)={\left\vert x(t) \right\vert}^{-2}\textstyle \int\limits_{-\infty}^{\infty} \displaystyle f\times W_x(t,f) \ df  </math>,  <math>-\Psi'(f)={\left\vert X(f) \right\vert}^{-2}\textstyle \int\limits_{-\infty}^{\infty} \displaystyle t \times W_x(t,f) \ dt  </math>
|-
|(5)Moment特性||<math>\iint_{-\infty}^{\infty} t^nW_x(t,f)\, dt\, df=\int_{-\infty}^{\infty} t^n|x(t)|^2\, dt</math>, <math>\iint_{-\infty}^{\infty} f^nW_x(t,f)\, dt\, df=\int_{-\infty}^{\infty} f^n|X(f)|^2\, dt</math>
|-
|(6)<math>W_x(t,f)</math>是實數||<math>\overline{W_x(t,f)}=W_x(t,f)</math>
|-
|(7)區域特性||If <math>x(t)=0</math> for <math>t>t_0</math> then <math>W_x(t,f)=0</math> for <math>t>t_0</math>, If <math>x(t)=0</math> for <math>t<t_0</math> then <math>W_x(t,f)=0</math> for <math>t<t_0</math>
|-
|(8)乘法特性||If <math>y(t)=x(t)h(t)</math>,then <math>W_y(t,f)=\int_{-\infty}^{\infty} W_x(t,\rho)W_h(t-\rho,f)\, d\rho</math>
|-
|(9)摺積特性||If <math>y(t)=\int_{-\infty}^{\infty} x(t-\tau)h(\tau)\, d\tau</math>,then <math>W_y(t,f)=\int_{-\infty}^{\infty} W_x(\rho,f)W_h(t-\rho,f)\, d\rho</math>
|-
|(10)相關特性||If <math>y(t)=\int_{-\infty}^{\infty} x(t+\tau)h^*(\tau)\, d\tau</math>,then <math>W_y(t,\omega)=\int_{-\infty}^{\infty} W_x(\rho,\omega)W_h(-t+\rho,\omega)\, d\rho</math>
|-
|(11)時間平移特性||If <math>y(t)=x(t-t_0)</math>, then <math>W_y(t,f)=W_x(t-t_0,f)</math>
|-
|(12)調變特性||If <math>y(t)=e^{j2\pi f_0t}x(t)</math>, then <math>W_y(t,f)=W_x(t,f-f_0)</math>
|}

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

# 直接運算(暴力法)            複雜度：<math>TF(2Q+1)  </math>
# 使用[[離散傅立葉變換]]        複雜度：<math>TN{\log _2}N</math>
# 使用[[Chirp-Z 轉換]] 複雜度 : <math>TN{\log _2}N</math>，通常為使用[[離散傅立葉變換]]的2~3倍，但限制比使用[[離散傅立葉變換]]少

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

從定義一出發

<math>{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau</math>

令 <math>\tau'=\tau/2</math>

<math>{W_x}\left( {t,f} \right) = 2\int_{ - \infty }^\infty {x\left( {t + \tau '} \right) \cdot } {x^*}\left( {t - \tau '} \right)\,{e^{ - j4\pi \,\tau '\,f}} \cdot d\tau '</math>

再令 <math>t = n\Delta_t, f = m\Delta_f,\tau'= p\Delta_t </math> ，則上述式子則為

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2\sum\limits_{p = - \infty }^\infty {x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right)\exp \left( { - j4\pi \,mp{\Delta _t}{\Delta _f}} \right){\Delta _t}}</math>

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

注意事項 :
 
若x(t)是無限長的訊號，則p要從負無限加到正無限，這點不易實現。

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


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

我們假設<math>\ x(t)=0 \ \ for \ \ t < n_1\Delta_t \ \ and \ \ t > n_2\Delta_t </math>

<BR>
===直接運算(暴力法) ===
<BR>
限制條件 : 

只有一個 : 要滿足Nyquist criterion

<math>{\Delta _t} < \frac{1}{{2B}}</math>，其中B是<math>x\left( {t + \tau } \right){x^*}\left( {t - \tau } \right)\,</math>的頻寬，大約是x(t)的兩倍。

推導 :

<math>\ x(t)=0 \ \ for \ \ t < n_1\Delta_t \ \ and \ \ t > n_2\Delta_t </math>

所以當<math>n + p \notin [{n_1},{n_2}]{\rm{ or n - p}} \notin {\rm{[}}{n_1},{n_2}{\rm{]}}</math> 時，

<math>x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right) = 0</math>

固定中間的n值 (<math>n\Delta_t</math>) 來探討p的範圍

<math>{n_1} \le n + p \le {n_2} \to {n_1} - n \le p \le {n_2} - n</math>

即 <math>\max ({n_1} - n,n - {n_2}) \le p \le \min ({n_2} - n,n - {n_1}) </math>-– (1)

<math>{n_1} \le n - p \le {n_2} \to {\rm{ }}{{\rm{n}}_1} - n \le - p \le {n_2} - n \to n - {n_2} \le p \le n - {n_1}</math>

即 <math> - \min ({n_2} - n,n - {n_1}) \le p \le \min ({n_2} - n,n - {n_1})</math>-- (2)

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

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

<math> - \min ({n_2} - n,n - {n_1}) \le p \le \min ({n_2} - n,n - {n_1})</math>

且可以發現 <math>({n_2} - n){\Delta _t},(n - {n_1}){\Delta _t}</math> 代表 <math>n\Delta_t</math> 離兩個邊界的距離

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

最後推導出的式子如下

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2\sum\limits_{p = - Q}^Q {x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right)\exp \left( { - j4\pi \,mp{\Delta _t}{\Delta _f}} \right){\Delta _t}} </math>

其中 <math>Q = \min ({n_2} - n,n - {n_1})\,,p \in [ - Q,Q],n \in [{n_1},{n_2}]</math>

===使用離散傅立葉變換===
<BR>
限制條件 :

(1)要滿足Nyquist criterion

<math>{\Delta _t} < \frac{1}{{2B}}</math>，其中B是<math>x\left( {t + \tau } \right){x^*}\left( {t - \tau } \right)\,</math>的頻寬，大約是x(t)的兩倍。

(2)<math>{\Delta _t}{\Delta _f} = {\textstyle{1 \over {2N}}}</math>

(3) <math>N \ge 2Q+1</math>

推導 :

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

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2{\Delta _t}\sum\limits_{p = - Q}^Q {x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right){e^{ - j{\textstyle{{2\pi \,mp} \over N}}}}}</math>

令<math>q=p+Q \to p=q-Q</math>

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2{\Delta _t}{e^{j{\textstyle{{2\pi \,mQ} \over N}}}}\sum\limits_{q = 0}^{2Q} {x\left( {(n + q - Q){\Delta _t}} \right){x^ * }\left( {(n - q + Q){\Delta _t}} \right){e^{ - j{\textstyle{{2\pi \,mq} \over N}}}}} </math>

針對中間<math>x\left( {(n + q - Q){\Delta _t}} \right){x^ * }\left( {(n - q + Q){\Delta _t}} \right)</math>項

令

<math>{c_1}\left( q \right) = x\left( {(n + q - Q){\Delta _t}} \right){x^ * }\left( {(n - q + Q){\Delta _t}} \right){\rm{ }} , for{\rm{ 0}} \le q \le 2{\rm{Q}}</math>

<math>{c_1}\left( q \right) = 0{\rm{ }}\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad, for{\rm{ 2}}Q + 1 \le q \le N - 1</math>


最後得出的式子如下

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2{\Delta _t}{e^{j{\textstyle{{2\pi \,mQ} \over N}}}}\sum\limits_{q = 0}^{N - 1} {{c_1}\left( q \right){e^{ - j{\textstyle{{2\pi \,mq} \over N}}}}}</math>

其中

<math>Q = \min ({n_2} - n,n - {n_1})\,,n \in [{n_1},{n_2}]</math>

<math>{c_1}\left( q \right) = x\left( {(n + q - Q){\Delta _t}} \right){x^ * }\left( {(n - q + Q){\Delta _t}} \right){\rm{ }} , for{\rm{ 0}} \le q \le 2{\rm{Q}}</math>

<math>{c_1}\left( q \right) = 0{\rm{ }}\qquad \qquad \qquad \qquad \qquad \qquad \qquad \qquad, for{\rm{ 2}}Q + 1 \le q \le N - 1</math>

===使用Chirp-Z 轉換===
<BR>
限制條件 : 

只有一個 : 要滿足Nyquist criterion

<math>{\Delta _t} < \frac{1}{{2B}}</math>，其中B是<math>x\left( {t + \tau } \right){x^*}\left( {t - \tau } \right)\,</math>的頻寬，大約是x(t)的兩倍。

推導 :

前提討論的式子可改寫為

<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2{\Delta _t}\,{e^{ - j2\pi \,{m^2}{\Delta _t}{\Delta _f}}}\sum\limits_{p = - Q}^Q {x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right){e^{ - j2\pi \,{p^2}{\Delta _t}{\Delta _f}}}{e^{j2\pi \,{{(p - m)}^2}{\Delta _t}{\Delta _f}}}} 
</math>

計算分成3步驟

STEP 1 : <math>{x_1}\left( {n,p} \right) = x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right){e^{ - j2\pi \,{p^2}{\Delta _t}{\Delta _f}}}</math>

STEP 2 : <math>{X_2}\left[ {n,m} \right] = \sum\limits_{p = n - Q}^{n + Q} {{x_1}\left[ p \right]\,c\left[ {m - p} \right]}</math>  , 其中<math>c\left[ m \right] = {e^{j2\pi \,{m^2}{\Delta _t}{\Delta _f}}}</math>

STEP 3 : <math>X\left( {n{\Delta _t},m{\Delta _f}} \right) = 2{\Delta _t}\,{e^{ - j2\pi \;{m^2}{\Delta _t}{\Delta _f}}}{X_2}\left[ {n,m} \right]</math>

==延伸變化==

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

==== 定義 ====
<math>{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {w\left( \tau \right)x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau</math> ， 其中<math>w(\tau)</math>為實數且為有限長訊號 

原始[[韋格納分佈]]定義 <math>{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau</math>

==== 優缺點 ====

* 優點 :  

# 降低運算時間，因為<math>w(\tau)</math>為有限長函數。 
# 可以有效降低相交項(cross term)問題，但不能完全消除(詳見下方說明)。 

* 缺點 :  

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

==== 實現方法 ====
從定義出發 
<math>{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {w\left( \tau \right)x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau </math> 
令<math>\tau=\tau'/2</math> 
<math>{W_x}\left( {t,f} \right) = 2\int_{ - \infty }^\infty {w\left( {2\tau '} \right)x\left( {t + \tau '} \right) \cdot } {x^*}\left( {t - \tau '} \right)\,{e^{ - j4\pi \,\tau '\,f}} \cdot d\tau '</math> 
再令<math>t = n\Delta_t, f = m\Delta_f,\tau'= p\Delta_t </math> 
<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2\sum\limits_{p = - \infty }^\infty {w\left( {2p{\Delta _t}} \right)x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right){e^{ - j4\pi \,mp{\Delta _t}{\Delta _f}}}{\Delta _t}} </math> 
假設w(t) = 0 for |t| > B 
即 <math>w\left( {2p{\Delta _t}} \right) = 0{\rm{ }} \ for{\rm{ }} \ p < - Q{\rm{ }}\ \land \ p > Q</math> 
其中<math>Q = \frac{B}{{2{\Delta _t}}}</math> 
如此一來，p範圍便可縮小。 
<math>{W_x}\left( {n{\Delta _t},m{\Delta _f}} \right) = 2\sum\limits_{p = - Q}^Q {w\left( {2p} \right)x\left( {(n + p){\Delta _t}} \right){x^ * }\left( {(n - p){\Delta _t}} \right){e^{ - j4\pi \,mp{\Delta _t}{\Delta _f}}}{\Delta _t}} </math>

==== 避免相交項的原因 ====
從定義出發 <math>{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {w\left( \tau \right)x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau </math>，其中<math>w(\tau)</math>為實數且為有限長訊號  

假設<math>x(t){\rm{ = }}\delta (t - {t_1}) + \delta (t - {t_2})</math>的情況下，比較有無mask function所產生的不同結果 
[[File:X(t).jpg|none|x(t)示意圖|替代=|缩略图|550x550像素]]
理想情形 : <math>{W_x}(t,f){\rm{ = 0 }} \ for{\rm{ }} \ t \ne {t_1},{t_2}</math>

<br />

===== 沒有使用mask function =====
即mask function <math>w(\tau ) = 1</math> 
 
<math>\begin{array}{l}
{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty {x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau \\
 = \int_{ - \infty }^\infty {\left[ {\delta \left( {t + \frac{\tau }{2} - {t_1}} \right) + \delta \left( {t + \frac{\tau }{2} - {t_2}} \right)} \right] \cdot } \left[ {\delta \left( {t - \frac{\tau }{2} - {t_1}} \right) + \delta \left( {t - \frac{\tau }{2} - {t_2}} \right)} \right]{e^{ - j2\pi \,\tau \,f}} \cdot d\tau \\
 = 4\int_{ - \infty }^\infty \underbrace{{\left[ {\delta \left( {\tau + 2t - 2{t_1}} \right) + \delta \left( {\tau + 2t - 2{t_2}} \right)} \right]  }}_{first \ term}\cdot  \underbrace{ \left[ {\delta \left( {\tau - 2t + 2{t_1}} \right) + \delta \left( {\tau - 2t + 2{t_2}} \right)} \right]}_{second \ term}{e^{ - j2\pi \,\tau \,f}} \cdot d\tau 
\end{array} 
</math> 
[[File:Wx(t,f).jpg|none|Wx(t,f)示意圖|替代=|缩略图|550x550像素]]
 
總共有3種情況要討論，如下圖，可見cross term在沒有使用mask function時，無法被消除

[[File:Without mask.jpg|none|ideal x(t)。Auto term 為自相關項。 Cross term為相交項|替代=|缩略图|500x500像素]]<br />

===== 使用mask function =====

<math>\begin{array}{l}
{W_x}\left( {t,f} \right) = \int_{ - \infty }^\infty w(\tau ) {x\left( {t + \tau /2} \right) \cdot } {x^*}\left( {t - \tau /2} \right)\,{e^{ - j2\pi \,\tau \,f}} \cdot d\tau \\
= 4\int_{ - \infty }^\infty w(\tau ) \underbrace{{\left[ {\delta \left( {\tau + 2t - 2{t_1}} \right) + \delta \left( {\tau + 2t - 2{t_2}} \right)} \right]  }}_{first \ term}\cdot  \underbrace{ \left[ {\delta \left( {\tau - 2t + 2{t_1}} \right) + \delta \left( {\tau - 2t + 2{t_2}} \right)} \right]}_{second \ term}{e^{ - j2\pi \,\tau \,f}} \cdot d\tau 
\end{array} 
</math> 

假設<math>w(\tau ) = 0{\rm{ }} \ for|\tau | > B > 0</math>，且<math>B{\rm{ < }}{{\rm{t}}_2}{\rm{ - }}{{\rm{t}}_1}</math>

由於<math>w(\tau )</math>只在-B到B有值，故乘上<math>w(\tau )</math>就能去除相交項(Cross term)，只保留下圖中兩條紅線中間的區域，也就是Auto terms。 
[[File:With mask.jpg|none|ideal x(t)。|替代=|缩略图|500x500像素]]
但上述其實是理想的情況，x(t)為窄頻信號Delta function  
 
如果X(τ)寬度太寬或是有ripple的話，Cross term仍會有殘留，示意圖如下 
[[File:With mask nonideal.jpg|none|non ideal x(t)。|替代=|缩略图|490x490像素]]
藍色線為X(τ)的訊號，若X(τ)的寬度太寬或是有ripple產生，就可能會跑進<math>w(\tau )</math>的範圍裡面，進而導致無法完全濾除Cross term。

==== 總結 ====
cross term 只有在訊號每個成分的寬度都小於2B，且時間差<math>t_2-t_1</math>都大於B時，才能被消除

===[[改進型韋格納分佈]]===
此方法可以消除相交項(cross term)。

=== [[加伯–韋格納轉換]] ===
消除相交項(cross term)問題，在某些情況下比加伯轉換擁有更好的清晰度。

== 参见 ==
* [[闵可夫斯基空间]]
* [[柯西不等式]]
* [[三角不等式]]
* [[完备空间]]

== 參考書目、資料來源 ==

* Jian-Jiun Ding, Time frequency analysis and wavelet transform class note, the Department of Electrical Engineering, National Taiwan University (NTU), Taipei, Taiwan, 2018.

[[Category:连续分布]]