Fourier Cosine Transform

In mathematics, the Fourier sine and cosine transforms are forms of the Fourier transform that do not use complex numbers or require negative frequency. They are the forms originally used by Joseph Fourier and are still preferred in some applications, such as signal processing or statistics.

Definition

The Fourier sine transform of , sometimes denoted by either $^s$ or $_s (f)$, is
$^s(\xi) = \int_^\infty f(t)\sin(2\pi \xi t) \,dt.$

If means time, then is frequency in cycles per unit time, but in the abstract, they can be any pair of variables which are dual to each other.

This transform is necessarily an odd function of frequency, i.e. for all :
$^s(-\xi) = - ^s(\xi).$

The numerical factors in the Fourier transforms are defined uniquely only by their product. Here, in order that the Fourier inversion formula not have any numerical factor, the factor of 2 appears because the sine function has norm of $\tfrac.$

The Fourier cosine transform of , sometimes denoted by either $^c$ or $_c (f)$, is
$^c(\xi) = \int_^\infty f(t)\cos(2\pi \xi t) \,dt.$

It is necessarily an even function of frequency, i.e. for all :
$^c(\xi) = ^c(-\xi).$Since positive frequencies can fully express the transform, the non-trivial concept of negative frequency needed in the regular Fourier transform can be avoided.

Simplification to avoid negative t

Some authors only define the cosine transform for even functions of , in which case its sine transform is zero. Since cosine is also even, a simpler formula can be used,
$^c(\xi)=2 \int_0^\infty f(t)\cos(2\pi \xi t) \,dt.$

Similarly, if is an odd function, then the cosine transform is zero and the sine transform can be simplified to
$^s(\xi) =2\int_0^\infty f(t)\sin(2\pi \xi t) \,dt.$

Other conventions

Just like the Fourier transform takes the form of different equations with different constant factors (see ), other authors also define the cosine transform as
$^c(\xi)=\sqrt \int_0^\infty f(t)\cos(2\pi \xi t) \,dt.$
and sine as
$^s(\xi) =\sqrt \int_0^\infty f(t)\sin(2\pi \xi t) \,dt,$
or, the cosine transform as $F_c(\alpha) = \frac \int_0^\infty f(x) \cos(\alpha x) \, dx$ and the sine transform as $F_s(\alpha) = \frac \int_0^\infty f(x) \sin(\alpha x) \, dx$ using $\alpha$ as the transformation variable. And while is typically used to represent the time domain, is often used alternatively, particularly when representing frequencies in a spatial domain.

Fourier inversion

The original function can be recovered from its transform under the usual hypotheses, that and both of its transforms should be absolutely integrable. For more details on the different hypotheses, see Fourier inversion theorem.

The inversion formula is
$f(t) = \int _^\infty ^c(\xi) \cos (2\pi \xi t) \,d\xi + \int _^\infty ^s(\xi) \sin (2\pi \xi t) \,d\xi,$

which has the advantage that all quantities are real. Using the addition formula for cosine, this can be rewritten as
$f(t) = \int _^\infty \int_^\infty f(\tau) \cos (2\pi\xi (\tau-t) ) \,d\tau\,d\xi.$

If the original function is an even function, then the sine transform is zero; if is an odd function, then the cosine transform is zero. In either case, the inversion formula simplifies.

Relation with complex exponentials

The form of the Fourier transform used more often today is
\begin
\hat(\xi) &= \int_^\infty f(t) e^\,dt \\
&= \int_^\infty f(t) \left(\cos (2\pi\xi t) - i\,\sin (2\pi\xi t)\right) dt && \text \\
&= \left (\int_^\infty f(t)\cos(2\pi \xi t) \,dt \right ) - i \left (\int_^\infty f(t)\sin(2\pi \xi t) \,dt \right )\\ pt&= ^c (\xi) - i ^s (\xi)
\end{align}

Numerical Evaluation

Using standard methods of numerical evaluation for Fourier integrals, such as Gaussian or tanh-sinh quadrature, is likely to lead to completely incorrect results, as the quadrature sum is (for most integrands of interest) highly ill-conditioned.
Special numerical methods which exploit the structure of the oscillation are required, an example of which is Ooura's method for Fourier integralsTakuya Ooura, Masatake Mori, ''A robust double exponential formula for Fourier-type integrals'', Journal of computational and applied mathematics 112.1-2 (1999): 229-241. This method attempts to evaluate the integrand at locations which asymptotically approach the zeros of the oscillation (either the sine or cosine), quickly reducing the magnitude of positive and negative terms which are summed.

See also

* Discrete cosine transform
* Discrete sine transform

References

* Whittaker, Edmund, and James Watson, ''A Course in Modern Analysis'', Fourth Edition, Cambridge Univ. Press, 1927, pp. 189, 211

Integral transforms
Fourier analysis
Mathematical physics