Transformation between distributions in time–frequency analysis: Difference between revisions

Revision as of 22:35, 10 May 2012

Latest revision as of 07:11, 18 January 2022

(30 intermediate revisions by 16 users not shown)

Line 1:

In the field of [[time–frequency analysis]], several signal formulations are used to represent the signal in a joint time–frequency ___domain.<ref>L. Cohen, "Time–Frequency Analysis," ''Prentice-Hall'', New York, 1995. {{ISBN|978-0135945322}}</ref>

{{Expert-subject|Mathematics|date=January 2010}}

In the field of [[time–frequency analysis]], the goal is to define signal formulations that are used for representing the signal in a joint time–frequency ___domain (see also [[time–frequency representation]]s<ref>B. Boashash, “Time-Frequency Concepts”, Chapter 1, pp. 3–28, in B. Boashash, ed,, Time-Frequency Signal Analysis & Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003; ISBN 008044335.</ref>). There are several methods and transforms called "time-frequency distributions" (TFDs).<ref>B. Boashash, “Heuristic Formulation of Time-Frequency Distributions”, Chapter 2, pp. 29–58, in B. Boashash, editor, Time-Frequency Signal Analysis and Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003; ISBN 008044335.</ref> The most useful and used methods form a class referred to as "quadratic" or [[bilinear time–frequency distribution]]s. A core member of this class is the [[Wigner–Ville distribution]] (WVD),<ref>B. Boashash, "Note on the Use of the Wigner Distribution for Time Frequency Signal Analysis", IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 36, No. 9, pp. 1518–1521, Sept. 1988</ref> as all other TFDs can be written as a smoothed version of the WVD. Another popular member of this class is the [[spectrogram]] which is the square of the magnitude of the [[short-time Fourier transform]] (STFT). The spectrogram has the advantage of being positive and is easy to interpret, but has disadvantages like being irreversible which means that once the spectrogram of a signal is computed, the original signal can't be extracted from the spectrogram. The theory and methodology for defining a TFD that verifies certain desirable properties is given in the "Theory of Quadratic TFDs" <ref>B. Boashash, “Theory of Quadratic TFDs”, Chapter 3, pp. 59–82, in B. Boashash, editor, Time-Frequency Signal Analysis & Processing: A Comprehensive Reference, Elsevier, Oxford, 2003; ISBN 0-08-044335-4.</ref>.

The scope of this article is to outline some elements of the procedure to transform one distribution into another. The method used to transform a distribution is borrowed from [[quantum mechanics]], even though the subject matter of the article is "signal processing". Noting that a signal can recovered from a particular distribution under certain conditions, given a certain TFD ρ1(t,f) representing the signal in a joint time–frequency ___domain, another different TFD ρ2(t,f) of the same signal can be obtained to calculate any other distribution, by simple smoothing or filtering; some of these relationships are shown below. A full treatment of the question can be given from a signal processing perspective.<ref>B. Boashash, editor, “Time-Frequency Signal Analysis and Processing – A Comprehensive Reference”, Elsevier Science, Oxford, 2003; ISBN 0-08-044335-4</ref>.

There are several methods and transforms called "time-frequency distributions" (TFDs), whose interconnections were organized by Leon Cohen.<ref>L. Cohen, "Generalized phase-space distribution functions," ''J. Math. Phys.'', '''7''' (1966) pp.&nbsp;781–786, [https://dx.doi.org/10.1063/1.1931206 doi:10.1063/1.1931206] </ref>

==General class==

<ref>L. Cohen, "Quantization Problem and Variational Principle in the Phase Space Formulation of Quantum Mechanics," ''J. Math. Phys.'', '''7''' pp.&nbsp;1863–1866, 1976.</ref><ref name="Janssen">A. J. E. M. Janssen, "On the locus and spread of pseudo-density functions in the time frequency plane," ''Philips Journal of Research'', vol. 37, pp.&nbsp;79–110, 1982.</ref><ref>E. Sejdić, I. Djurović, J. Jiang, “Time-frequency feature representation using energy concentration: An overview of recent advances,” ''Digital Signal Processing'', vol. 19, no. 1, pp. 153-183, January 2009.</ref>

In the field of [[time–frequency analysis]], the goal is to define signal formulations that are used for representing the signal in a joint time–frequency ___domain (see also [[time–frequency representation]]s<ref>B. Boashash, “Time-Frequency Concepts”, Chapter 1, pp. 3–28, in B. Boashash, ed,, Time-Frequency Signal Analysis & Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003; ISBN 008044335.</ref>). There are several methods and transforms called "time-frequency distributions" (TFDs).<ref>B. Boashash, “Heuristic Formulation of Time-Frequency Distributions”, Chapter 2, pp. 29–58, in B. Boashash, editor, Time-Frequency Signal Analysis and Processing: A Comprehensive Reference, Elsevier Science, Oxford, 2003; ISBN 008044335.</ref> The most useful and usedpopular methods form a class referred to as "quadratic" or [[bilinear time–frequency distribution]]s. A core member of this class is the [[Wigner–Ville distribution]] (WVD),<ref>B. Boashash, "Note on the Use of the Wigner Distribution for Time Frequency Signal Analysis", IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. 36, No. 9, pp. 1518–1521, Sept. 1988</ref> as all other TFDs can be written as a smoothed versionor convolved versions of the WVD. Another popular member of this class is the [[spectrogram]] which is the square of the magnitude of the [[short-time Fourier transform]] (STFT). The spectrogram has the advantage of being positive and is easy to interpret, but also has disadvantages, like being irreversible, which means that once the spectrogram of a signal is computed, the original signal can't be extracted from the spectrogram. The theory and methodology for defining a TFD that verifies certain desirable properties is given in the "Theory of Quadratic TFDs" .<ref>B. Boashash, “Theory of Quadratic TFDs”, Chapter 3, pp. 59–82, in B. Boashash, editor, Time-Frequency Signal Analysis & Processing: A Comprehensive Reference, Elsevier, Oxford, 2003; {{ISBN |0-08-044335-4}}.</ref>.

The scope of this article is to outlineillustrate some elements of the procedure to transform one distribution into another. The method used to transform a distribution is borrowed from the [[phase space formulation]] of [[quantum mechanics]], even though the subject matter of thethis article is "signal processing". Noting that a signal can be recovered from a particular distribution under certain conditions, given a certain TFD ρ1''ρ''₁(''t'',''f'') representing the signal in a joint time–frequency ___domain, another, different, TFD ρ2''ρ''₂(''t'',''f'') of the same signal can be obtained to calculate any other distribution, by simple smoothing or filtering; some of these relationships are shown below. A full treatment of the question can be given fromin aCohen's signal processing perspective.<ref>B. Boashash, editor, “Time-Frequency Signal Analysis and Processing – A Comprehensive Reference”, Elsevier Science, Oxford, 2003; ISBN 0-08-044335-4</ref>book.

If we use the variable ''ω''=2''πf'', then, borrowing the notations used in the field of quantum mechanics, we can show that time–frequency representation, such as [[Wigner distribution function]] (WDF) and other [[bilinear time–frequency distribution]]s, can be expressed as

==General class==

: <math>C(t,\omega) = \dfrac{1}{4\pi^2}\iiint s^*(u-\dfrac{1}{2}\tau)s(u+\dfrac{1}{2}\tau)\phi(\theta,\tau)e^{-j\theta t-j\tau\omega+j\theta u}\, du\,d\tau\,d\theta ,</math> (1)

If we use the variable {{math|1=''ω'' = 2''πf''}}, then, borrowing the notations used in the field of quantum mechanics, we can show that time–frequency representation, such as [[Wigner distribution function]] (WDF) and other [[bilinear time–frequency distribution]]s, can be expressed as

where <math>\phi(\theta,\tau)</math> is a two dimensional function called the kernel, which determines the distribution and its properties (for a signal processing terminology and treatment of this question, the reader is referred to the references already cited in the introduction).

: {{NumBlk||<math display="block">C(t,\omega) = \dfrac{1}{4\pi^2}\iiint s^*\left(u-\dfrac{1}{2}\tau\right) s\left(u+\dfrac{1}{2}\tau\right)\phi(\theta,\tau)e^{-j\theta t-j\tau\omega+j\theta u}\, du\,d\tau\,d\theta ,</math> (|{{EquationRef|1)}}}}

where <math>\phi(\theta,\tau)</math> is a two dimensional function called the kernel, which determines the distribution and its properties (for a signal processing terminology and treatment of this question, the reader is referred to the references already cited in the introduction).

For theThe kernel of the [[Wigner distribution function]] (WDF) is one. However, it is no particular significance should be attached to that, since it is possible to write the general form so that the kernel of any distribution is one, in which case the kernel of the [[Wigner distribution function]] (WDF) would be something else.

==Characteristic function formulation==

The characteristic function is the double [[Fourier transform]] of the distribution. By inspection of Eq. ({{EquationNote|1}}), we can obtain that

: {{NumBlk||<math display="block">C(t,\omega) = \dfrac{1}{4\pi^2}\iint M(\theta,\tau)e^{-j\theta t-j\tau\omega}\, d\theta\,d\tau</math> (|{{EquationRef|2)}}}}

where

: <math>C(t,\omega) = \dfrac{1}{4\pi^2}\iint M(\theta,\tau)e^{-j\theta t-j\tau\omega}\, d\theta\,d\tau</math> (2)

{{NumBlk||<math display="block">\begin{alignat}{2}

M(\theta,\tau) & = \phi(\theta,\tau)\int s^*\left(u-\dfrac{1}{2}\tau\right)s\left(u+\dfrac{1}{2}\tau\right)e^{j\theta u}\,du \\

where

& = \phi(\theta,\tau)A(\theta,\tau) \\

: <math>\beginend{alignat}</math> | {2{EquationRef|3}}}}

M(\theta,\tau) & = \phi(\theta,\tau)\int s^*(u-\dfrac{1}{2}\tau)s(u+\dfrac{1}{2}\tau)e^{j\theta u}\,du \\

& = \phi(\theta,\tau)A(\theta,\tau) \\

\end{alignat}</math> (3)

and where <math>A(\theta,\tau)</math> is the symmetrical ambiguity function. The characteristic function may be appropriately called the generalized ambiguity function.

Line 30 ⟶ 29:

To obtain that relationship suppose that there are two distributions, <math>C_1</math> and <math>C_2</math>, with corresponding kernels, <math>\phi_1</math> and <math>\phi_2</math>. Their characteristic functions are

: {{NumBlk||<math display="block">M_1(\phi,\tau) = \phi_1(\theta,\tau)\int s^*\left(u-\dfractfrac{1\tau}{2}\tauright)s\left(u+\dfractfrac{1\tau}{2}\tauright)e^{j\theta u}\, du</math> | ({{EquationRef|4)}}}}

: {{NumBlk||<math display="block">M_2(\phi,\tau) = \phi_2(\theta,\tau)\int s^*\left(u-\dfractfrac{1\tau}{2}\tauright)s\left(u+\dfractfrac{1\tau}{2}\tauright)e^{j\theta u}\, du</math> | ({{EquationRef|5)}}}}

Divide one equation by the other to obtain

: <math>M_1(\phi,\tau) = \phi_1(\theta,\tau)\int s^*(u-\dfrac{1}{2}\tau)s(u+\dfrac{1}{2}\tau)e^{j\theta u}\, du</math> (4)

: {{NumBlk||<math display="block">M_1(\phi,\tau) = \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}M_2(\phi,\tau)</math> (|{{EquationRef|6)}}}}

: <math>M_2(\phi,\tau) = \phi_2(\theta,\tau)\int s^*(u-\dfrac{1}{2}\tau)s(u+\dfrac{1}{2}\tau)e^{j\theta u}\, du</math> (5)

Divide one equation by the other to obtain

: <math>M_1(\phi,\tau) = \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}M_2(\phi,\tau)</math> (6)

This is an important relationship because it connects the characteristic functions. For the division to be proper the kernel cannot to be zero in a finite region.

To obtain the relationship between the distributions take the double [[Fourier transform]] of both sides and use Eq. ({{EquationNote|2}})

: {{NumBlk||<math display="block">C_1(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}M_2(\theta,\tau)e^{-j\theta t-j\tau\omega}\, d\theta\,d\tau</math> (|{{EquationRef|7)}}}}

: <math>C_1(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}M_2(\theta,\tau)e^{-j\theta t-j\tau\omega}\, d\theta\,d\tau</math> (7)

Now express <math>M_2</math> in terms of <math>C_2</math> to obtain

: {{NumBlk||<math display="block">C_1(t,\omega) = \dfrac{1}{4\pi^2}\iiiint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}C_2(t,\omega^')e^{j\theta(t^'-t)+j\tau(\omega^'-\omega)}\, d\theta\,d\tau\,dt^'\,d\omega^'</math> (|{{EquationRef|8)}}}}

: <math>C_1(t,\omega) = \dfrac{1}{4\pi^2}\iiiint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}C_2(t,\omega^')e^{j\theta(t^'-t)+j\tau(\omega^'-\omega)}\, d\theta\,d\tau\,dt^'\,d\omega^'</math> (8)

This relationship can be written as

: {{NumBlk||<math display="block">C_1(t,\omega) = \iint g_{12}(t^'-t,\omega^'-\omega)C_2(t,\omega^')\,dt^'\,d\omega^'</math> (|{{EquationRef|9)}}}}

with

: <math>C_1(t,\omega) = \iint g_{12}(t^'-t,\omega^'-\omega)C_2(t,\omega^')\,dt^'\,d\omega^'</math> (9)

: {{NumBlk||<math display="block">g_{12}(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}e^{j\theta t+j\tau\omega}\, d\theta\, d\tau</math> (|{{EquationRef|10)}}}}

with

: <math>g_{12}(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{\phi_1(\theta,\tau)}{\phi_2(\theta,\tau)}e^{j\theta t+j\tau\omega}\, d\theta\, d\tau</math> (10)

==Relation of the spectrogram to other bilinear representations==

Now we specialize to the case where one transform from an arbitrary representation to the spectrogram. In Eq. ({{EquationNote|9}}), both <math>C_1</math> to be the spectrogram and <math>C_2</math> to be arbitrary are set. In addition, to simplify notation, <math>\phi_{SP} = \phi_1</math>, <math>\phi = \phi_2</math>, and <math>g_{SP} = g_{12}</math> are set and written as

: {{NumBlk||<math display="block">C_{SP}(t,\omega) = \iint g_{SP} \left (t^'-t,\omega^'-\omega \right)C \left (t,\omega^' \right )\,dt^'\,d\omega^'</math> (|{{EquationRef|11)}}}}

: <math>C_{SP}(t,\omega) = \iint g_{SP}(t^'-t,\omega^'-\omega)C(t,\omega^')\,dt^'\,d\omega^'</math> (11)

The kernel for the spectrogram with window, <math>h(t)</math>, is <math>A_h(-\theta,\tau)</math> and therefore

: <math display="block">\begin{alignat}{3align}

g_{SP}(t,\omega) & = \dfrac{1}{4\pi^2}\iint \dfrac{A_h(-\theta,\tau)}{\phi(\theta,\tau)}e^{j\theta t+j\tau\omega}\, d\theta\,d\tau \\

& = \dfrac{1}{4\pi^2}\iiint \dfrac{1}{\phi(\theta,\tau)}h^*(u-\dfractfrac{1\tau}{2}\tau)h(u+\dfractfrac{1\tau}{2}\tau)e^{j\theta t+j\tau\omega-j\theta u}\, du\,d\tau\,d\theta \\

& = \dfrac{1}{4\pi^2}\iiint h^*(u-\dfractfrac{1\tau}{2}\tau)h(u+\dfractfrac{1\tau}{2}\tau)\dfrac{\phi(\theta,\tau)}{\phi(\theta,\tau)\phi(-\theta,\tau)}e^{-j\theta t+j\tau\omega+j\theta u}\, du\,d\tau\,d\theta \\

\end{alignatalign}</math> (12)

If takingwe theonly consider kernels for which <math>\phi(-\theta,\tau)\phi(\theta,\tau) = 1</math>, <math>g_{SP}(t,\omega)</math>holds is just the distribution of the window function, except that it is evaluated at <math>-\omega</math>. Therefore,then

<math display="block">g_{SP}(t,\omega) = \dfrac{1}{4\pi^2}\iiint h^*(u-\tfrac{\tau}{2})h(u+\tfrac{\tau}{2}) \phi(\theta,\tau) e^{-j\theta t+j\tau\omega+j\theta u}\, du\,d\tau\,d\theta = C_h(t,-\omega)</math>

and therefore

: <math>g_{SP}(t,\omega) = C_h(t,-\omega)</math> (13)

: <math display="block">C_{SP}(t,\omega) = \iint C_s(t^',\omega^')C_h(t^'-t,\omega^'-\omega)\,dt^'\,d\omega^'</math> (14)

for kernels that satisfy <math>\phi(-\theta,\tau)\phi(\theta,\tau) = 1</math>

and

: <math>C_{SP}(t,\omega) = \iint C_s(t^',\omega^')C_h(t^'-t,\omega^'-\omega)\,dt^'\,d\omega^'</math> (14)

for kernels that satisfy <math>\phi(-\theta,\tau)\phi(\theta,\tau) = 1</math>

This was shown by Janssen[4]. For the case where <math>\phi(-\theta,\tau)\phi(\theta,\tau)</math> does not equal one, then

: <math>C_{SP}(t,\omega) = \iiiint G(t^{''},\omega^{''})C_s(t^',\omega^')C_h(t^{''}+t^'-t,-\omega^{''}+\omega-\omega^')\,dt^'\,dt^{''}\,d\omega^\,d\omega^{''}</math> (15)

This was shown by Janssen[4].<ref Forname="Janssen"/> the case whereWhen <math>\phi(-\theta,\tau)\phi(\theta,\tau)</math> does not equal one, then

: <math display="block">C_{SP}(t,\omega) = \iiiint G(t^{''},\omega^{''})C_s(t^',\omega^')C_h(t^{''}+t^'-t,-\omega^{''}+\omega-\omega^')\,dt^'\,dt^{''}\,d\omega^'\,d\omega^{''}</math> (15)

where

: <math display="block">G(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{e^{-j\theta t-j\tau\omega}}{\phi(\theta,\tau)\phi(-\theta,\tau)}\, d\theta\,d\tau</math> (16)

: <math>G(t,\omega) = \dfrac{1}{4\pi^2}\iint \dfrac{e^{-j\theta t-j\tau\omega}}{\phi(\theta,\tau)\phi(-\theta,\tau)}\, d\theta\,d\tau</math> (16)

==References==

{{Reflist}}

{{DEFAULTSORT:Transformation between distributions in time-frequency analysis}}

[1] L. Cohen, "Time–Frequency Analaysis," ''Prentice-Hall'', New York, 1995.

[2] L. Cohen, "Generalized phase-space distribution functions," ''Jour. Math. Phys.'', vol.7, pp.&nbsp;781–786, 1966.

[3] L. Cohen, "Quantization Problem and Variational Principle in the Phase Space Formulation of Quantum Mechanics," ''Jour. Math. Phys.'', vol.7, pp.&nbsp;1863–1866, 1976.

[4] A. J. E. M. Janssen, "On the locus and spread of pseudo-density functions in the time frequency plane," ''Philips Journal of Research'', vol. 37, pp.&nbsp;79–110, 1982.

[[Category:Time–frequency analysis]]