 # Complex Fourier Series

In an earlier module, we showed that a square wave could be expressed as a superposition of pulses. As useful as this decomposition was in this example, it does not generalize well to other periodic signals: How can a superposition of pulses equal a smooth signal like a sinusoid? Because of the importance of sinusoids to linear systems, you might wonder whether they could be added together to represent a large number of periodic signals. You would be right and in good company as well. Euler and Gauss in particular worried about this problem, and Jean Baptiste Fourier got the credit even though tough mathematical issues were not settled until later. They worked on what is now known as the Fourier series: representing any periodic signal as a superposition of sinusoids.

But the Fourier series goes well beyond being another signal decomposition method. Rather, the Fourier series begins our journey to appreciate how a signal can be described in either the time-domain or the frequency-domain with no compromise. Let

be a periodic signal with period

. We want to show that periodic signals, even those that have constant-valued segments like a square wave, can be expressed as sum of harmonically related sine waves: sinusoids having frequencies that are integer multiples of the fundamental frequency. Because the signal has period

, the fundamental frequency is

. The complex Fourier series expresses the signal as a superposition of complex exponentials having frequencies

,

.

with

.The real and imaginary parts of the Fourier coefficients

are written in this unusual way for convenience in defining the classic Fourier series. The zeroth coefficient equals the signal's average value and is real- valued for real-valued signals:

. The family of functions

are called basis functions and form the foundation of the Fourier series. No matter what the periodic signal might be, these functions are always present and form the representation's building blocks. They depend on the signal period

, and are indexed by

.

### Key point:

Assuming we know the period, knowing the Fourier coefficients is equivalent to knowing the signal. Thus, it makes no difference if we have a time-domain or a frequency- domain characterization of the signal.

### Exercise

What is the complex Fourier series for a sinusoid?

Because of Euler's relation,

Thus,

, and the other coefficients are zero.

To find the Fourier coefficients, we note the orthogonality property

Assuming for the moment that the complex Fourier series "works," we can find a signal's complex Fourier coefficients, its spectrum, by exploiting the orthogonality properties of harmonically related complex exponentials. Simply multiply each side of Equation by

and integrate over the interval

.

Finding the Fourier series coefficients for the square wave

is very simple. Mathematically, this signal can be expressed as

The expression for the Fourier coefficients has the form

When integrating an expression containing i, treat it just like any other constant.

The two integrals are very similar, one equaling the negative of the other. The final expression becomes

Consequently, the square wave equals a sum of complex exponentials, but only those having frequencies equal to odd multiples of the fundamental frequency

. The coefficients decay slowly as the frequency index

increases. This index corresponds to the

-th harmonic of the signal's period.

A signal's Fourier series spectrum

has interesting properties.

Property:

If

is real,

(real-valued periodic signals have conjugate-symmetric spectra).

This result follows from the integral that calculates the

from the signal. Furthermore, this result means that

: The real part of the Fourier coefficients for real-valued signals is even. Similarly,

: The imaginary parts of the Fourier coefficients have odd symmetry. Consequently, if you are given the Fourier coefficients for positive indices and zero and are told the signal is real-valued, you can find the negative-indexed coefficients, hence the entire spectrum. This kind of symmetry,

is known as conjugate symmetry.

Property:

If

, which says the signal has even symmetry about the origin,

.

Given the previous property for real-valued signals, the Fourier coefficients of even signals are real-valued. A real-valued Fourier expansion amounts to an expansion in terms of only cosines, which is the simplest example of an even signal.

Property:

If

,

which says the signal has odd symmetry,

.

Therefore, the Fourier coefficients are purely imaginary. The square wave is a great example of an odd-symmetric signal.

Property:

The spectral coefficients for a periodic signal delayed by

,

are

where

denotes the spectrum of

. Delaying a signal by

seconds results in a spectrum having a linear phase shift of

in comparison to the spectrum of the undelayed signal. Note that the spectral magnitude is unaffected. Showing this property is easy.

Note that the range of integration extends over a period of the integrand. Consequently, it should not matter how we integrate over a period, which means that

and we have our result.

The complex Fourier series obeys Parseval's Theorem, one of the most important results in signal analysis. This general mathematical result says you can calculate a signal's power in either the time domain or the frequency domain.

## Parseval's Theorem

Average power calculated in the time domain equals the power calculated in the frequency domain.

This result is a (simpler) re-expression of how to calculate a signal's power than with the real-valued Fourier series expression for power.

Let's calculate the Fourier coefficients of the periodic pulse signal shown here.

The pulse width is

,the period

, and the amplitude

. The complex Fourier spectrum of this signal is given by

At this point, simplifying this expression requires knowing an interesting property.

Armed with this result, we can simply express the Fourier series coefficients for our pulse sequence.

Because this signal is real-valued, we find that the coefficients do indeed have conjugate symmetry:

. The periodic pulse signal has neither even nor odd symmetry; consequently, no additional symmetry exists in the spectrum. Because the spectrum is complex valued, to plot it we need to calculate its magnitude and phase.

The function

equals -1 if its argument is negative and zero otherwise. The somewhat complicated expression for the phase results because the sine term can be negative; magnitudes must be positive, leaving the occasional negative values to be accounted for as a phase shift of

.

#### Periodic Pulse Sequence

Also note the presence of a linear phase term (the first term in

is proportional to frequency

). Comparing this term with that predicted from delaying a signal, a delay of

is present in our signal. Advancing the signal by this amount centers the pulse about the origin, leaving an even signal, which in turn means that its spectrum is real-valued. Thus, our calculated spectrum is consistent with the properties of the Fourier spectrum.

### Exercise

What is the value of

? Recalling that this spectral coefficient corresponds to the signal's average value, does your answer make sense?

. This quantity clearly corresponds to the periodic pulse signal's average value.

The phase plot shown in Figure 2 requires some explanation as it does not seem to agree with what Equation suggests. There, the phase has a linear component, with a jump of

every time the sinusoidal term changes sign. We must realize that any integer multiple of

can be added to a phase at each frequency without affecting the value of the complex spectrum. We see that at frequency index 4 the phase is nearly

. The phase at index 5 is undefined because the magnitude is zero in this example. At index 6, the formula suggests that the phase of the linear term should be less than

(more negative). In addition, we expect a shift of

in the phase between indices 4 and 6. Thus, the phase value predicted by the formula is a little less than

without affecting the value of the spectrum at index 6, the result is a slightly negative number as shown. Thus, the formula and the plot do agree. In phase calculations like those made in MATLAB, values are usually confined to the range

by adding some (possibly negative) multiple of

to each phase value. Swipe left and right to change pages.  