Physics Library
 An open source physics library
Encyclopedia | Forums | Docs | Random  

Equivalence of Real and Complex Fourier Forms

(Topic)

FS20: Equivalence of Real and Complex Fourier Forms

The real Fourier series and the complex Fourier series are not two different theories of periodic functions. They describe the same harmonic components using different coordinates. In FS04–FS06 we extracted coefficients of cosine and sine functions. FS18 rewrote those functions as oppositely rotating complex exponentials, and FS19 derived the complex coefficients by orthogonal projection. Now we will join those two developments into a reversible calculation.

A key conceptual question is this: if a real signal has one cosine-and-sine harmonic at frequency n, why does its complex spectrum have two coefficients, one at +n and one at −n? The answer is not additional physical energy or a new frequency. It is a different bookkeeping system for the very same real oscillation.

The central results we will earn for n ≥ 1 are

      an − ibn             an + ibn
Cn =  --------,    C −n =  --------,
         2                    2
(1)

and

an = Cn + C −n,     bn = i(Cn − C −n).
(2)

We will also show that C0 = a0∕2 and, for real-valued f, that C−n = Cn.

PIC

Figure 1. The real and complex descriptions are two invertible coordinate systems for one harmonic. Neither introduces a new physical oscillation.

1 Begin with the two formulas on the same interval

Let f be a real-valued, 2L-periodic function. Write

𝜃 = πx-.
     L
(3)

We use 𝜃 only as an abbreviation; the period in x remains 2L. The real Fourier form is

             ∑∞
f (x) ∼ a0-+     [a  cos(n𝜃) + b sin(n𝜃)].
        2          n           n
             n=1
(4)

The complex Fourier form is

        ∑∞
f(x) ∼       Cnein𝜃.
       n=− ∞
(5)

In the first sum, the index starts at 1, and there are two real coefficients for each positive harmonic. In the second sum, the index spans all integers and there is one complex coefficient for each signed harmonic. The complex sum also contains n = 0.

When making algebraic transformations, begin with finite sums. This avoids silently exchanging infinite sums or limits. Set

               ∑N
RN  (x) = a0-+     [an cos(n𝜃) + bnsin(n𝜃)],
          2    n=1
(6)

and

          ∑N
QN (x) =       Cnein𝜃.
         n=− N
(7)

Our first task is to choose Cn so that QN(x) = RN(x) identically, for every x, without taking N to infinity.

2 One cosine requires a positive and a negative frequency

Euler’s formula reads

 iα
e  =  cosα + isinα.
(8)

Changing the sign of α gives

 −iα
e    = cos α − isin α.
(9)

Add these equations. The imaginary terms cancel:

eiα + e− iα = 2 cosα.
(10)

Subtract instead. The real parts cancel:

eiα − e−iα = 2isin α.
(11)

Therefore

        eiα-+-e−iα
cos α =      2    ,
(12)

and

       eiα − e−iα
sin α = ----------.
           2i
(13)

These are identities, not approximations.

An important geometric point follows. The phasor ein𝜃 rotates counterclockwise as 𝜃 increases for positive n, while e−in𝜃 rotates in the opposite direction. Their sum has no imaginary component and produces a cosine. Their suitably scaled difference produces a sine.

PIC

Figure 2. Euler’s formula relates complex rotation to real cosine and sine components. A real sinusoid can be assembled from a conjugate pair of counter-rotating complex modes.

3 Derive the conversion from real coefficients to complex coefficients

Take just one real harmonic, not the whole series:

hn (x ) = ancos(n𝜃) + bn sin(n𝜃 ).
(14)

Replace cosine and sine using Equations (12) and (13):

hn(x) = an-
 2( in𝜃   −in𝜃)
 e   + e
+ bn
---
2i(           )
 ein𝜃 − e−in𝜃. (15)

The only algebraic subtlety is that

1
i-=  − i
(16)

because i(−i) = 1. Therefore the sine part can be rewritten as

bn                  ibn      ibn
--(ein𝜃 − e−in𝜃) = −---ein𝜃 + ---e− in𝜃.
2i                   2        2
(17)

Now collect the multiplier of ein𝜃 and the multiplier of e−in𝜃:

hn(x) = (          )
  an-−  ibn-
   2     2ein𝜃
+ (          )
  an-+  ibn-
   2     2e−in𝜃. (18)

Compare this equation with

            in𝜃        −in𝜃
hn (x) = Cne   +  C−ne    .
(19)

The coefficients must be

C  =  an −-ibn-,    n ≥  1,
 n       2
(20)

and

       an-+-ibn
C −n =     2   ,     n ≥ 1.
(21)

What happens to the constant term?

At zero frequency, ei0𝜃 = 1. There is only one zero-frequency mode, not a positive-negative pair. The constant term in Equation (6) is already a0∕2. Thus

      a0
C0 =  --.
      2
(22)

Do not apply Equation (20) at n = 0 while also inventing a separate C−0. That would count the DC component twice.

PIC

Figure 3. The real cosine coefficient contributes equally to +n and −n, while the real sine coefficient contributes oppositely signed imaginary parts.

4 Derive the reverse conversion from complex to real

Now suppose the complex coefficients Cn and C−n are known. Start with the two equations just obtained:

2Cn =  an − ibn,
(23)

2C −n = an + ibn.
(24)

Add them:

2Cn  + 2C −n = 2an,
(25)

so

an = Cn +  C− n.
(26)

Subtract the second equation from the first:

2(Cn − C −n) = − 2ibn.
(27)

Divide by −2i, or multiply by i∕2:

bn = i(Cn  − C −n).
(28)

Finally,

a0 = 2C0.
(29)

The two forms are therefore related by a reversible linear change of coordinates. We did not approximate anything in either direction.

Why real signals have conjugate pairs

For a real-valued function, the Fourier coefficients an and bn are real numbers. Conjugating Equation (20) gives

---
Cn =  an-+-ibn-=  C− n.
         2
(30)

Consequently,

an = 2 Re Cn,
(31)

and

b  = − 2 Im  C .
 n           n
(32)

The minus sign in the formula for bn follows from our e+in𝜃 synthesis convention. Other Fourier sign conventions can reverse it, so the safe approach is to derive the conversion once from Euler’s formula.

We can verify the real-valued nature of a conjugate pair directly:

Cnein𝜃 + C −ne−in𝜃 = C nein𝜃 + C nein𝜃
= 2 Re (      )
 Cnein𝜃. (33)

Thus a conjugate pair is exactly the real harmonic an cos(n𝜃) + bn sin(n𝜃).

5 Worked example 1: reconstruct a finite signal both ways

Reuse the FS19 waveform

f (x) = 2 + 3cosx −  2sin(2x),    − π ≤  x ≤ π.
(34)

From the real form we read

a0 = 4,  a1 =  3,  b2 = − 2,
(35)

with every other nonzero-index coefficient equal to zero. Use our conversion formulas one harmonic at a time.

At zero frequency,

C0 =  a0∕2 = 2.
(36)

At the first harmonic, b1 = 0, giving

C  =  3-−-i0 = 3-,  C   =  3.
  1     2      2      −1   2
(37)

At the second harmonic, a2 = 0 and b2 = −2, giving

      0 −-i(−-2)-             0-+-i(−-2)
C2 =      2     = i,  C −2 =      2     = − i.
(38)

Thus the complex form is the finite expression

f(x) = 2 + 3-eix + 3-e−ix + ie2ix − ie−2ix.
           2      2
(39)

Group the pair at n = ±1:

3(eix + e− ix) = 3cos x.
2
(40)

Group the pair at n = ±2:

i(e2ix − e−2ix) = − 2sin2x.
(41)

We recover Equation (34) exactly.

As an independent check, use the reverse formulas. For n = 2,

a2 = C2 + C− 2 = i − i = 0,
(42)

and

b2 = i(C2 − C −2) = i(2i) = − 2.
(43)

Notice how the sine coefficient is stored as imaginary components with opposite signs in the complex spectrum.

PIC

Figure 4. The same finite signal written in real and complex coefficient form. The zero-frequency real coefficient is a0 = 4 whereas C0 = 2; the other harmonics are divided between paired positive and negative frequencies.

6 Amplitude and phase: interpreting a complex coefficient

We often express one real harmonic as

h  (x) = A  cos(n𝜃 − φ ).
  n       n           n
(44)

Expand the cosine of a difference:

hn(x) = An cos φn cos(n𝜃)
+ An sin φn sin(n𝜃). (45)

Compare with an cos(n𝜃) + bn sin(n𝜃) to identify

an =  An cosφn,     bn = An sin φn.
(46)

Therefore

      ∘ --------
An =    a2n + b2n.
(47)

The phase φn is chosen with the correct quadrant, for example using the two-argument arctangent atan2(bn,an).

Substitute Equation (46) into the complex conversion:

Cn = A
--n
 2(cos φn − isinφn )
= An-
 2e−iφn . (48)

Hence

|C  | = An,     argC   = − φ    (C  ⁄= 0).
   n    2            n      n     n
(49)

The complex spectrum has two conjugate coefficients of magnitude An∕2, but the real sinusoid has amplitude An. It would be wrong to read |Cn| alone as the real amplitude without accounting for the negative-frequency partner.

A quick phase example

Take

            (      )
h(x) = 2 cos x −  π- .
                  6
(50)

The real coefficients are a1 = √ --
  3 and b1 = 1. The complex coefficients are

C1  = e−iπ∕6,    C −1 = eiπ∕6.
(51)

Each has magnitude 1, while the recovered real sinusoid has amplitude 2. The negative argument of C1 is consistent with our cos(x − φ) phase convention.

PIC

Figure 5. A positive-frequency coefficient and its conjugate form symmetric phasors. The coefficient phase records the shift of the real sinusoid; two complex magnitudes An∕2 combine to recover amplitude An.

7 Worked example 2: the odd square wave

Recall from FS07 the odd square wave of amplitude A on [−π,π]:

       {
         − A,  − π < x <  0,
f(x) =
         A,    0 < x <  π.
(52)

It is odd, so a0 = an = 0. Its sine coefficients are

bn =  2A-[1 − (− 1)n].
      πn
(53)

For positive n, the complex coefficient is therefore

C  =  − iA-[1 − (− 1)n].
  n     πn
(54)

For odd positive n this reduces to

       2iA-
Cn = −  πn .
(55)

Its negative-frequency partner is

       2iA
C− n = ----.
       πn
(56)

Both even-harmonic coefficients are zero; the conjugate pair at each odd harmonic recreates the sine term. Notice that |Cn| = 2A∕(πn) whereas the real sine amplitude is bn = 4A∕(πn).

PIC

Figure 6. The real odd square wave has positive coefficients on odd sine harmonics. Its complex spectrum stores them as an antisymmetric imaginary pair at positive and negative indices. The vertical scales reflect the factor of two.

8 Worked example 3: an even rectangular pulse

Define a periodic centered pulse by

       {
         A,   |x | < d,
f(x ) =                        0 < d < L.
         0,   d < |x| ≤ L,
(57)

Because the pulse is even, all bn = 0. Its constant term is

       ∫ d
a =  1-    A dx =  2Ad-.
 0   L  − d         L
(58)

For n ≥ 1 the cosine coefficient is

an = 2A-
 L∫ 0d cos (     )
  nπx-
   L dx
= 2A
---
nπ sin (     )
  nπd
  ----
   L. (59)

Thus the complex coefficients are

C  =  Ad-,
 0    L
(60)

and

                   (     )
              A      n πd
Cn =  C− n = nπ-sin  --L-  ,     n ≥ 1.
(61)

These coefficients are real and even in their index. This is the spectral signature of an even, real waveform about the chosen origin. The sign of a coefficient can still be negative, so a plot of coefficient magnitude alone loses information.

PIC

Figure 7. A centered even pulse and its real, even two-sided complex coefficient sequence. The horizontal waveform coordinate u = 3x∕L is used for the illustrated choice d = L∕3.

9 Time shifts change phase, not coefficient magnitude

For a periodic time signal, replace 𝜃 by ω0t, where ω0 = 2π∕T. Suppose

         ∞
        ∑       inω0t
f (t) =       Cne    .
       n=− ∞
(62)

Shift the entire waveform to the right by t0:

g(t) = f(t − t0)
= ∑ nCneinω0(t−t0)
= ∑ n(          )
 Cne −inω0t0einω0t. (63)

Therefore

C (g) = C (f)e− inω0t0.
  n      n
(64)

Since |e−inω0t0| = 1, the magnitude stays unchanged:

|C (ng)| = |C (nf)|.
(65)

A delay changes the phases by −nω0t0; harmonics at higher n rotate farther.

To see the same transformation in real coefficients, set α = nω0t0 and expand a single shifted harmonic. One finds

 (g)    (f)         (f)
an  =  an  cosα − bn  sin α,
(66)

 (g)    (f)         (f)
bn  =  an sin α + bn cos α.
(67)

These formulas are a rotation in the (an,bn) coefficient plane. In the complex plane, the same rotation is expressed more economically by multiplication by one exponential.

A quick check

Start with f(t) = cos(ω0t) and delay by T∕4. Then α = π∕2, and the shifted waveform is sin(ω0t). The real formulas yield a1(g) = 0 and b 1(g) = 1. The complex formula yields C 1(g) = −i∕2, C−1(g) = i∕2. Both descriptions agree.

10 Finite energy bookkeeping also agrees

FS12 established that orthogonal projection splits the squared norm among retained directions. For a real Fourier partial sum on [−L,L], orthogonality gives

∫ L               La2      ∑N
    |RN  (x)|2 dx = ---0 + L    (a2n + b2n).
 −L                 2      n=1
(68)

Using the complex basis, whose squared norm is 2L, we obtain

∫                       N
   L        2          ∑       2
     |QN  (x )|dx  = 2L      |Cn |.
  −L                  n=−N
(69)

But the finite sums are the same function. The right sides must agree. Indeed, for each positive n,

|Cn|2 + |C −n|2 = a2+  b2
-n---n--
  4 + a2 + b2
-n----n-
   4
=  2   2
an +-bn--
  2. (70)

Multiplying by 2L gives precisely L(an2 + b n2). The DC contribution is 2L|C 0|2 = La 02∕2.

This identity is finite, so it does not yet require a theorem about the convergence of the infinite Fourier series. FS13 develops the infinite-series energy statement, Parseval’s theorem, under the appropriate completeness assumptions.

11 What if the function itself is complex-valued?

Our algebraic identities converting an,bn and C±n do not require f to be real, provided the real-form coefficients are allowed to be complex. However, the symmetry

       ---
C− n = Cn
(71)

does require the signal to be real-valued (in the usual almost-everywhere sense). For example,

f(x) = eix
(72)

has C1 = 1 and C−1 = 0. This is perfectly valid for a complex signal; no conjugate partner is required.

Complex signals are useful mathematical objects, and in electronics the in-phase/quadrature signal is commonly represented as a complex envelope. The physical meaning of this complex notation depends on the system and should not be confused with the automatically real quantity Re{ei𝜃}.

12 Physics example: periodic forcing of a damped oscillator

Consider the linear mass-spring-damper equation

M ¨y + by˙+ ky = F (t).
(73)

Suppose the forcing is a real periodic signal of fundamental angular frequency ω0. Write the complex expansions

        ∑      inω0t            ∑      inω0t
F (t) =    Fne     ,    y (t) =     Yne    .
         n                      n
(74)

For each exponential, differentiation multiplies by inω0; differentiating twice multiplies by −(nω0)2. Consequently the equation for each harmonic becomes

[                     ]
 k − M (nω0 )2 + ibnω0 Yn = Fn.
(75)

For a stable forced response,

               F
Yn =  ----------n----------.
      k − M (nω0 )2 + ibnω0
(76)

A numerical two-harmonic forcing

Let M = 1 kg, b = 2 N s∕m, k = 25 N∕m, and ω0 = 3 rad∕s. Apply

F (t) = 6cos(3t) + 4sin(6t)  (N).
(77)

The nonzero force coefficients are

F1 = F− 1 = 3,    F2 = − 2i,     F−2 = 2i.
(78)

At n = 1, the response is

Y 1 = -----3-----
25 − 9 + 6i = ---3---
16 + 6i
= 48-−-18i-
  292. (79)

At n = 2,

Y 2 =     − 2i
--------------
25 − 36 + 12i =    − 2i
----------
− 11 + 12i
= −-24 +-22i
   265. (80)

The negative-frequency responses are their conjugates because the forcing and the system coefficients are real. Convert the two positive-frequency coefficients back to the real basis:

               24                       9
a1 = 2 Re Y1 = ---,    b1 = − 2Im Y1 =  --,
               73                       73
(81)

                  48                         44
a2 = 2 Re Y2 = − ----,    b2 = − 2Im Y2 =  − ---.
                 265                         265
(82)

Therefore the real steady-state displacement is

y(t) = 24
---
73 cos(3t) + 9
---
73 sin(3t)
−-48-
265 cos(6t) −44--
265 sin(6t), (83)

with displacement measured in meters. Direct substitution into Equation (73) reproduces the original force.

PIC

Figure 8. A two-harmonic periodic force and the corresponding steady displacement. Complex coefficients make each harmonic response algebraic; conjugate pairing converts the result into the real displacement observed in the laboratory.

13 Choosing which representation to use

The real form is often the more intuitive choice when drawing a waveform, separating even and odd contributions, specifying a directly measured cosine or sine amplitude, or teaching elementary harmonic synthesis. The complex form is usually more compact when differentiating harmonics, handling shifts, keeping phase and amplitude together, or solving linear differential equations mode by mode.

Neither form is intrinsically more accurate. Their finite partial sums are identical functions when the coefficients are converted correctly. Any difference in an actual computation comes from truncation, rounding, or using different conventions, not from a mathematical difference between real and complex Fourier series.

Common convention traps

  • Mixing C0 and a0: the constant term is C0 = a0∕2, not a0.
  • Ignoring negative frequencies: for a real sinusoid, two complex exponentials reconstruct one real harmonic.
  • Forgetting the imaginary sign: with synthesis e+in𝜃, b n = −2 Im Cn for positive n.
  • Reading |Cn| as the real amplitude: for n ≥ 1, the real amplitude is 2|Cn|.
  • Treating ein𝜃 as automatically real: a single complex exponential is generally complex; reality requires conjugate pairing.
  • Confusing series with transforms: a periodic signal has discrete harmonic indices. FS26 develops the different continuous-frequency Fourier transform.

14 Self-checks

  1. If f(x) = 5 + 4 cos 3x, find C0,C3,C−3 and all other nonzero coefficients. Answer: C0 = 5, C±3 = 2.
  2. If C2 = 1 − 2i and f is real, what is C−2 and what real harmonic do they represent? Answer: C−2 = 1 + 2i and 2 cos 2x + 4 sin 2x.
  3. What is the amplitude of a real harmonic whose positive-frequency coefficient has modulus 0.4? Answer: 0.8.
  4. Can C−1 = 0 and C1 = 1 describe a real-valued waveform? Answer: not by themselves; they give eix, a complex waveform.
  5. Why does cos x delayed by one quarter-period become sin x? Answer: the positive-frequency coefficient rotates by e−iπ∕2 = −i; the real coefficient vector rotates from (a1,b1) = (1, 0) to (0, 1).

15 What FS20 has established

From FS02’s projection idea to FS19’s complex coefficient formula, we have arrived at a complete dictionary between real and complex representations. Euler’s formula pairs two counter-rotating exponentials into one real harmonic. Their complex coefficients store both the amplitude and phase, and their conjugate symmetry guarantees a real result. The finite projection energy is identical in the two descriptions, and linear physical systems can be solved harmonic by harmonic in complex form before returning to measurable real variables.

This prepares two important directions: FS21 develops Fourier modes as eigenfunctions in separation-of-variables problems, and FS26 explains how discrete harmonics become a continuous-frequency integral for nonperiodic signals.

References

[1]   E. Kreyszig, Advanced Engineering Mathematics, 10th ed., Wiley, 2011.

[2]   R. N. Bracewell, The Fourier Transform and Its Applications, 3rd ed., McGraw-Hill, 2000.

[3]   A. V. Oppenheim and A. S. Willsky, Signals and Systems, 2nd ed., Prentice Hall, 1997.

[4]   G. B. Arfken, H. J. Weber, and F. E. Harris, Mathematical Methods for Physicists, 7th ed., Academic Press, 2012.


"Equivalence of Real and Complex Fourier Forms" is owned by bloftin.
(view preamble)
View style:
Other names:  FS20
Also defines:  real-complex Fourier coefficient conversion, conjugate-pair reconstruction, two-sided Fourier spectrum, one-sided harmonic amplitude, complex coefficient phase, Fourier phase convention
Keywords:  Fourier series, real Fourier series, complex Fourier series, Euler formula, complex coefficients, conjugate symmetry, positive frequency, negative frequency, amplitude spectrum, phase spectrum, orthogonal projection, harmonic forcing

Cross-references: representations, FS02's, vector, Fourier transform, FS26, computation, partial sums, differential equations, displacement, force, angular frequency, theorem, Fourier partial sum, norm, complex plane, periodic, odd square wave, FS07, amplitude, magnitude, phase, negative frequencies, formula, Fourier coefficients, DC component, phasor, identities, Euler's formula, algebraic, period, coordinate systems, oscillation, system, energy, spectrum, frequency, FS19, complex exponentials, FS18, functions, harmonic, periodic functions, complex Fourier series, Fourier series
There are 3 references to this object.

This is version 1 of Equivalence of Real and Complex Fourier Forms, born on 2026-11-11.
Object id is 1580, canonical name is EquivalenceOfRealAndComplexFourierForms.
Accessed 8 times total.

Classification:
Physics Classification: 02.30.Nw (Fourier analysis)
 02.30.Lt (Sequences, series, and summability)
 02.10.Ud (Linear algebra)

Pending Errata and Addenda

None.

Discussion

No messages.

Interact