next up previous
Next: About this document ... Up: Time-Frequency Reassignment (Charpentier 1986, Previous: Framework: the Wigner-Ville Distribution

Implementation via STFTs for Audio (Fitz, Haken, and Christensen, ICMC 2000)

For some discrete time domain signal $x(n)$, equations 10 and 11 become
$\displaystyle \hat{t}_{k,n}$ $\textstyle =$ $\displaystyle t_n - \Re\left\{\ensuremath{\frac{X_{th;k,n} X^*_{h;k,n}}{\vert X_{h;k,n}\vert}}\right\}$ (12)
$\displaystyle \hat{\omega}_{k,n}$ $\textstyle =$ $\displaystyle \omega_k +\Im\left\{\ensuremath{\frac{X_{dh;k,n} X^*_{h;k,n}}{\vert X_{h;k,n}\vert}}\right\}$ (13)

Where the authors use a conventional DFT of a windowed signal:
$\displaystyle X_{h;k,n} = \sum_{l=-\infty}^\infty h_{n-l} x_l exp(-j2\pi l k/N)$     (14)


and also define time and frequency derivatives of the window and use them in ``modified'' DFTs:
$\displaystyle X_{dh;k,n}$ $\textstyle =$ $\displaystyle \sum_{l=-\infty}^\infty \ensuremath{\frac{d h_{n-l}}{dt}}x_l exp(-j2\pi l k/N)$ (15)
$\displaystyle X_{th;k,n}$ $\textstyle =$ $\displaystyle \sum_{l=-\infty}^\infty t\cdot h_{n-l} x_l exp(-j2\pi l k/N)$ (16)



(See overhead figures for illustrations of windows.)

NOTES:
next up previous
Next: About this document ... Up: Time-Frequency Reassignment (Charpentier 1986, Previous: Framework: the Wigner-Ville Distribution
Aaron S. Master 2003-02-04