Signals occur in a wide range of physical phenomenon. They might be human speech, blood pressure variations with time, seismic waves, radar and sonar signals, pictures or images, stress and strain signals in a building structure, stock market prices, a city's population, or temperature across a plate. These signals are often modeled or represented by a real or complex valued mathematical function of one or more variables. For example, speech is modeled by a function representing air pressure varying with time. The function is acting as a mathematical analogy to the speech signal and, therefore, is called an analog signal. For these signals, the independent variable is time and it changes continuously so that the term continuous-time signal is also used. In our discussion, we talk of the mathematical function as the signal even though it is really a model or representation of the physical signal.
The description of signals in terms of their sinusoidal frequency content has proven to be one of the most powerful tools of continuous and discrete-time signal description, analysis, and processing. For that reason, we will start the discussion of signals with a development of Fourier transform methods. We will first review the continuous-time methods of the Fourier series (FS), the Fourier transform or integral (FT), and the Laplace transform (LT). Next the discrete-time methods will be developed in more detail with the discrete Fourier transform (DFT) applied to finite length signals followed by the discrete-time Fourier transform (DTFT) for infinitely long signals and ending with the Z-transform which allows the powerful tools of complex variable theory to be applied.
More recently, a new tool has been developed for the analysis of signals. Wavelets and wavelet transforms 9, 1, 5, 16, 15 are another more flexible expansion system that also can describe continuous and discrete-time, finite or infinite duration signals. We will very briefly introduce the ideas behind wavelet-based signal analysis.
The problem of expanding a finite length signal in a trigonometric series was posed and studied in the late 1700's by renowned mathematicians such as Bernoulli, d'Alembert, Euler, Lagrange, and Gauss. Indeed, what we now call the Fourier series and the formulas for the coefficients were used by Euler in 1780. However, it was the presentation in 1807 and the paper in 1822 by Fourier stating that an arbitrary function could be represented by a series of sines and cosines that brought the problem to everyone's attention and started serious theoretical investigations and practical applications that continue to this day 8, 3, 11, 10, 7, 12. The theoretical work has been at the center of analysis and the practical applications have been of major significance in virtually every field of quantitative science and technology. For these reasons and others, the Fourier series is worth our serious attention in a study of signal processing.
We assume that the signal x(t) to be analyzed is well described by a real or complex valued function of a real variable t defined over a finite interval {0≤t≤T}. The trigonometric series expansion of x(t) is given by
where xk(t)=cos(2πkt/T) and yk(t)=sin(2πkt/T) are the basis functions for the expansion. The energy or power in an electrical, mechanical, etc. system is a function of the square of voltage, current, velocity, pressure, etc. For this reason, the natural setting for a representation of signals is the Hilbert space of L2[0,T]. This modern formulation of the problem is developed in 6, 11. The sinusoidal basis functions in the trigonometric expansion form a complete orthogonal set in L2[0,T]. The orthogonality is easily seen from inner products
and
where δ(t) is the Kronecker delta function with δ(0)=1 and δ(k≠0)=0. Because of this, the kth coefficients in the series can be found by taking the inner product of x(t) with the kth basis functions. This gives for the coefficients
and
where T is the time interval of interest or the period of a periodic signal. Because of the orthogonality of the basis functions, a finite Fourier series formed by truncating the infinite series is an optimal least squared error approximation to x(t). If the finite series is defined by
the squared error is
which is minimized over all a(k) and b(k) by Equation 1.4 and Equation 1.5. This is an extraordinarily important property.
It follows that if x(t)∈L2[0,T], then the series converges to x(t) in the sense that ε→0 as N→∞6, 11. The question of point-wise convergence is more difficult. A sufficient condition that is adequate for most application states: If f(x) is bounded, is piece-wise continuous, and has no more than a finite number of maxima over an interval, the Fourier series converges point-wise to f(x) at all points of continuity and to the arithmetic mean at points of discontinuities. If f(x) is continuous, the series converges uniformly at all points 11, 8, 3.
A useful condition 6, 11 states that if x(t) and its derivatives through the qth derivative are defined and have bounded variation, the Fourier coefficients a(k) and b(k) asymptotically drop off at least as fast as as k→∞. This ties global rates of convergence of the coefficients to local smoothness conditions of the function.
The form of the Fourier series using both sines and cosines makes determination of the peak value or of the location of a particular frequency term difficult. A different form that explicitly gives the peak value of the sinusoid of that frequency and the location or phase shift of that sinusoid is given by
and, using Euler's relation and the usual electrical engineering notation of ,
the complex exponential form is obtained as
where
The coefficient equation is
The coefficients in these three forms are related by
and
It is easier to evaluate a signal in terms of c(k) or d(k) and θ(k) than in terms of a(k) and b(k). The first two are polar representation of a complex value and the last is rectangular. The exponential form is easier to work with mathematically.
Although the function to be expanded is defined only over a specific finite region, the series converges to a function that is defined over the real line and is periodic. It is equal to the original function over the region of definition and is a periodic extension outside of the region. Indeed, one could artificially extend the given function at the outset and then the expansion would converge everywhere.
It can be very helpful to develop a geometric view of the Fourier series where x(t) is considered to be a vector and the basis functions are the coordinate or basis vectors. The coefficients become the projections of x(t) on the coordinates. The ideas of a measure of distance, size, and orthogonality are important and the definition of error is easy to picture. This is done in 6, 11, 17 using Hilbert space methods.
The properties of the Fourier series are important in applying it to signal analysis and to interpreting it. The main properties are given here using the notation that the Fourier series of a real valued function x(t) over {0≤t≤T} is given by F{x(t)}=c(k) and denotes the periodic extensions of x(t).
Linear: F{x+y}=F{x}+F{y}
Idea of superposition. Also scalability: F{ax}=aF{x}
Extensions of x(t):
is periodic.
Even and Odd Parts: x(t)=u(t)+jv(t) and
u | v | A | B | |C| | θ |
even | 0 | even | 0 | even | 0 |
odd | 0 | 0 | odd | even | 0 |
0 | even | 0 | even | even | π/2 |
0 | odd | odd | 0 | even | π/2 |
Convolution: If continuous cyclic convolution is defined by
then
Multiplication: If discrete convolution is defined by
then
This property is the inverse of property 4 and vice versa.
Parseval:
This property says the energy calculated in the time domain is the same as
that calculated in the frequency (or Fourier) domain.
Shift:
A shift in the time domain results in a linear phase shift in the frequency domain.
Modulate:
Modulation in the time domain results in a shift in the frequency domain. This
property is the inverse of property 7.
Orthogonality of basis functions:
Orthogonality allows the calculation of coefficients using inner products in Equation 1.4 and Equation 1.5. It also allows Parseval's Theorem in property 6. A relaxed version of orthogonality is called “tight frames" and is important in over-specified systems, especially in wavelets.
An example of the Fourier series is the expansion of a square wave signal with period 2π. The expansion is
Because x(t) is odd, there are no cosine terms (all a(k)=0) and,
because of its symmetries, there are no even harmonics (even k terms are
zero). The function is well defined and bounded; its derivative is not,
therefore, the coefficients drop off as