Modified Wigner distribution function

• Wigner distribution

${\displaystyle W_{x}(t,f)=\int _{-\infty }^{\infty }x(t+\tau /2)x^{*}(t-\tau /2)e^{-j2\pi \tau f}\,d\tau }$

Cohen's kernel function : ${\displaystyle \Pi (t,f)=\delta _{(0,0)}(t,f)}$

• Spectrogram

${\displaystyle SP_{x}(t,f)=|ST_{x}(t,f)|^{2}=ST_{x}(t,f)\,ST_{x}^{*}(t,f)}$

where ${\displaystyle ST_{x}}$ is the short-time Fourier transform of ${\displaystyle x}$.

${\displaystyle ST_{x}(t,f)=\int _{-\infty }^{\infty }x(\tau )w^{*}(t-\tau )e^{-j2\pi f\tau }\,d\tau }$

Cohen's kernel function : ${\displaystyle \Pi (t,f)=W_{h}(t,f)}$ which is the WD of the window function itself. This can be verified by applying the convolution property of the Wigner distribution function.

The spectrogram cannot produce interference since it is a positive-valued quadratic distribution.

• Modified form I

${\displaystyle W_{x}(t,f)=\int _{-B}^{B}w(\tau )x(t+\tau /2)x^{*}(t-\tau /2)e^{-j2\pi \tau f}\,d\tau }$

Can't solve the cross term problem, however it can solve the problem of 2 components time difference larger than window size B.

• Modified form II

${\displaystyle W_{x}(t,f)=\int _{-B}^{B}w(\eta )X(f+\eta /2)X^{*}(f-\eta /2)e^{j2\pi t\eta }\,d\eta }$

• Modified form III (Pseudo L-Wigner Distribution)

${\displaystyle W_{x}(t,f)=\int _{-\infty }^{\infty }w(\tau )x^{L}(r+\tau /2L){\overline {x^{*L}(t-\tau /2L)}}e^{-j2\pi \tau f}\,d\tau }$

Where L is any integer greater than 0

Increase L can reduce the influence of cross term (however it can't eliminate completely )

For example, for L=2, the dominant third term is divided by 4 ( which is equivalent to 12dB ).

This gives a significant improvement over the Wigner Distribution.

Properties of L-Wigner Distribution:

1. The L-Wigner Distribution is always real.
2. If the signal is time shifted ${\displaystyle x(t-t0)}$, then its LWD is time shifted as well, ${\displaystyle LWD:W_{x}(t-t0,f)}$
3. The LWD of a modulated signal ${\displaystyle x(t)\exp(j\omega _{0}t)}$is shifted in frequency ${\displaystyle LWD:W_{x}(t,f-f0)}$
4. Is the signal ${\displaystyle x(t)}$is time limited, i.e.,${\displaystyle x(t)=0}$ ${\displaystyle for\left\vert t\right\vert >T,}$ then the L-Wigner distribution is time limited, ${\displaystyle LWD:W_{x}(t,f)=0}$ ${\displaystyle for\left\vert t\right\vert >T}$
5. If the signal ${\displaystyle x(t)}$is band limited with ${\displaystyle f_{m}}$(${\displaystyle F(f)=0}$${\displaystyle for\left\vert f\right\vert >f_{m}}$), then ${\displaystyle LWD:W_{x}(t,f)}$is limited in the frequency domain by ${\displaystyle f_{m}}$as well.
6. Integral of L-Wigner distribution over frequency is equal to the generalized signal power: ${\displaystyle \int _{-\infty }^{\infty }W_{x}(t,f)df=\left\vert x(t)\right\vert ^{2L}}$
7. Integral of ${\displaystyle LWD:W_{x}(t,f)}$over time and frequency is equal to the ${\displaystyle 2L^{th}}$power of the ${\displaystyle 2L^{th}}$norm of signal ${\displaystyle x(t)}$:

${\displaystyle \int _{-\infty }^{\infty }\int _{-\infty }^{\infty }W_{x}(t,f)dtdf=\int _{-\infty }^{\infty }\left\vert x(t)\right\vert ^{2L}dt=\lVert x(t)\rVert _{2L}^{2L}}$

1. The integral over time is:

${\displaystyle \int _{-\infty }^{\infty }W_{x}(t,f)dt=\left\vert F_{L}(f)\right\vert ^{2}=\left\vert \underbrace {F(L_{f})*F(L_{f})*\cdots *F(L_{f})} _{Ltimes}\right\vert ^{2}}$

1. For a large value of ${\displaystyle L(L\rightarrow \infty )}$We may neglect all values of ${\displaystyle LWD:W_{x}(t,f)}$, Comparing them to the one at the points ${\displaystyle (t_{m},f_{m})}$, where the distribution reaches its essential supremum:

${\displaystyle \lim _{L\to \infty }(W_{x}(t,f)/W_{x}(t_{m},f_{m}))={\begin{cases}0,&{\text{if }}f\neq f_{m}{\text{ or }}t\neq t_{m}{\text{ }}\\1,&{\text{if }}f=f_{m}{\text{ and }}t=t_{m}\end{cases}}}$

• Modified form IV (Polynomial Wigner Distribution Function)

${\displaystyle W_{x}(t,f)=\int _{-B}^{B}[\textstyle \prod _{l=1}^{q/2}\displaystyle x(t+d_{l}\tau )x^{*}(t-d_{-l}\tau )]e^{-j2\pi \tau f}\,d\tau }$

When ${\displaystyle q=2}$ and ${\displaystyle d_{l}=d_{-l}=0.5}$, it becomes the original Wigner distribution function.

It can avoid the cross term when the order of phase of the exponential function is no larger than ${\displaystyle q/2+1}$

However the cross term between two components cannot be removed.

${\displaystyle d_{l}}$should be chosen properly such that

${\displaystyle \textstyle \prod _{l=1}^{q/2}\displaystyle x(t+d_{l}\tau )x^{*}(t-d_{-l}\tau )=\exp {\big (}j2\pi \textstyle \sum _{n=1}^{q/2+1}na_{n}t^{n-1}\tau \displaystyle {\big )}}$

${\displaystyle W_{x}(t,f)=\int _{-\infty }^{\infty }\exp {\Bigl (}-j2\pi (f-\sum _{n=1}^{q/2+1}na_{n}t^{n-1})\tau {\Bigr )}d\tau }$

${\displaystyle \cong \delta {\bigl (}f-\sum _{n=1}^{q/2+1}na_{n}t^{n-1}{\bigr )}}$

If ${\displaystyle x(t)=\exp {\bigl (}j2\pi \sum _{n=1}^{q/2+1}a_{n}t^{n}{\bigr )}}$

when ${\displaystyle q=2}$ , ${\displaystyle x(t+d_{l}\tau )x^{*}(t-d_{-l}\tau )=\exp {\bigl (}j2\pi \sum _{n=1}^{q/2+1}na_{n}t^{n-1}\tau {\bigr )}}$

${\displaystyle a_{2}(t+d_{l}\tau )^{2}+a_{1}(t+d_{l}\tau )-a_{2}(t-d_{-l}\tau )^{2}-a_{1}(t-d_{-l}\tau )=2a_{2}t\tau +a_{1}\tau }$

${\displaystyle \Longrightarrow d_{l}+d_{-l}=1,d_{l}-d_{-l}=0}$

${\displaystyle \Longrightarrow d_{l}=d_{-l}=1/2}$

• Pseudo Wigner distribution

${\displaystyle PW_{x}(t,f)=\int _{-\infty }^{\infty }w(\tau /2)w^{*}(-\tau /2)x(t+\tau /2)x^{*}(t-\tau /2)e^{-j2\pi \tau \,f}\,d\tau }$

Cohen's kernel function : ${\displaystyle \Pi (t,f)=\delta _{0}(t)\,W_{h}(t,f)}$ which is concentred on the frequency axis.

Note that the pseudo Wigner can also be written as the Fourier transform of the “spectral-correlation” of the STFT

${\displaystyle PW_{x}(t,f)=\int _{-\infty }^{\infty }ST_{x}(t,f+\nu /2)ST_{x}^{*}(t,f-\nu /2)e^{j2\pi \nu \,t}\,d\nu }$

• Smoothed pseudo Wigner distribution :

In the pseudo Wigner the time windowing acts as a frequency direction smoothing. Therefore, it suppresses the Wigner distribution interference components that oscillate in the frequency direction. Time direction smoothing can be implemented by a time-convolution of the PWD with a lowpass function ${\displaystyle q}$ :

${\displaystyle SPW_{x}(t,f)=[q\,\ast \,PW_{x}(.,f)](t)=\int _{-\infty }^{\infty }q(t-u)\int _{-\infty }^{\infty }w(\tau /2)w^{*}(-\tau /2)x(u+\tau /2)x^{*}(u-\tau /2)e^{-j2\pi \tau \,f}\,d\tau \,du}$

Cohen's kernel function : ${\displaystyle \Pi (t,f)=q(t)\,W(f)}$ where ${\displaystyle W}$ is the Fourier transform of the window ${\displaystyle w}$.

Thus the kernel corresponding to the smoothed pseudo Wigner distribution has a separable form. Note that even if the SPWD and the S-Method both smoothes the WD in the time domain, they are not equivalent in general.

• S-method

${\displaystyle SM(t,f)=\int _{-\infty }^{\infty }ST_{x}(t,f+\nu /2)ST_{x}^{*}(t,f-\nu /2)G(\nu )e^{j2\pi \nu \,t}\,d\nu }$

Cohen's kernel function : ${\displaystyle \Pi (t,f)=g(t)\,W_{h}(t,f)}$

The S-method limits the range of the integral of the PWD with a low-pass windowing function ${\displaystyle g(t)}$ of Fourier transform ${\displaystyle G(f)}$. This results in the cross-term removal, without blurring the auto-terms that are well-concentred along the frequency axis. The S-method strikes a balance in smoothing between the pseudo-Wigner distribution ${\displaystyle PW_{x}}$ [${\displaystyle g(t)=1}$] and the power spectrogram ${\displaystyle SP_{x}}$ [${\displaystyle g(t)=\delta _{0}(t)}$].

Note that in the original 1994 paper, Stankovic defines the S-methode with a modulated version of the short-time Fourier transform :

${\displaystyle SM(t,f)=\int _{-\infty }^{\infty }{\tilde {ST}}_{x}(t,f+\nu ){\tilde {ST}}_{x}^{*}(t,f-\nu )P(\nu )\,d\nu }$

where

${\displaystyle {\tilde {ST}}_{x}(t,f)=\int _{-\infty }^{\infty }x(t+\tau )w^{*}(\tau )e^{-j2\pi f\tau }\,d\tau \quad =ST_{x}(t,f)\,e^{j2\pi ft}}$

Even in this case we still have

${\displaystyle \Pi (t,f)=p(2t)\,W_{h}(t,f)}$

• Time-frequency representation
• Wigner distribution function
• Bilinear time–frequency distribution
• Short-time Fourier transform
• Gabor transform

• P. Gonçalves and R. Baraniuk, “Pseudo Affine Wigner Distributions : Definition and Kernel Formulation”, IEEE Transactions on Signal Processing, vol. 46, no. 6, Jun. 1998
• L. Stankovic, “A Method for Time-Frequency Analysis”, IEEE Transactions on Signal Processing, vol. 42, no. 1, Jan. 1994
• L. J. Stankovic, S. Stankovic, and E. Fakultet, “An analysis of instantaneous frequency representation using time frequency distributions-generalized Wigner distribution,” IEEE Trans. on Signal Processing, pp. 549-552, vol. 43, no. 2, Feb. 1995