*Pure mathematics is much more than an armory of tools and techniques for the applied mathematician. On the other hand, the pure mathematician has ever been grateful to applied mathematics for stimulus and inspiration. From the vibrations of the violin string they have drawn enchanting harmonies of Fourier Series, and to study the triode valve they have invented a whole theory of non-linear oscillations.*

**George Frederick James Temple** In *100 Years of Mathematics: a Personal Viewpoint* (1981).

The **Fourier Series** is a very important mathematics tool discovered by Jean-Baptiste Joseph Fourier in the 18th century. The Fourier series is used in many important areas of science and engineering. They are used to give an analytical approximate description of complex periodic function or series of data. In this blog, I am going to give a short introduction to it.

Let start with some definition. A function is said to be periodic with period T (T being a nonzero constant) if it if it is defined for all real and if there is some positive number, , such that

An example of periodic functions are the *sinusoidal function* (or *sinusoidal oscillation* or *sinusoidal signal*)} given by the function where

- : Amplitude. It gives how high the graph rises above the t-axis at its maximum values.
- : phase lag. Position of the maximum of the function.
- : the time delay or time shift. It tell you how far the graph of the has been shifted along the t-axis.
- : angular frequency. It gives the number of complete oscillation of in a time interval of length .
- : frequency of . It gives the number of complete oscillation of in a time interval of length 1.
- : period.Time required for one complete oscillation.

If the function is a periodic function with period The function in the interval can be represented by a trigonometric series as:

namely as a summation of sinusoidal functions of different frequency and amplitude that allows to approximate any periodic function with any level of accuracy.

In order to achieve this result we need to calculate the coefficients and in the series. But for that, we need to recall of properties of these functions and of their integrals.

### Integral of Orthogonal Functions

This means that for :

for :

This means for all value of that

This means that for :

for :

### Determination of

We start to integrate both the members of the Fourier series from to

For the integrals formulas that we have derived before, we obtain

namely,

is equal to the average value of the function

### Determination of

We start by multiplying both the members of the Fourier series by and then integrate from to

Shall we consider the integrals that we get on the right-hand-side:

First term:

The second integral for is equal to:

The third integral for the orthonormality of the function is always equal to:

Therefore:

namely, for :

### Determination of

We start by multiplying both the members of the Fourier series by and then integrate from to

Shall we consider the integrals that we get on the right-hand-side:

*First term:*

*Second term:*

*Third term:*

Therefore:

namely, for :

**Summary**

** The coefficients of the Fourier series **

By replacing m with n we can finally express the coefficient of the series in the final form.

with

or in the same way we can integrate in the interval with the same result:

For a given generic interval

with

and the Fourier series is give by

### Example 1

In this example, we will calculate the Fourier coefficient of a the *square wave *function*. *This function is of paramount importance in digital electronics since it used in the electronic signal to represent unit of information (bit) as on/off state.

Mathematically, the square wave is represented by functions as

The plot of the extended periodic function is given in the following figure.

We now calculate the Fourier coefficients starting from the term

The coefficients are given by

.

And finally

or

The Fourier series is

.

In this way, for example, the first 3 harmonics of the series are

### Example 2

In this second example the square wave signal is defined as a odd function such that and mathematically, the square wave is represented by functions as

The plot of the extended periodic function is given in the following figure.

We now calculate the Fourier coefficients starting from the term

The coefficients are given by

.

And finally

or

The Fourier series is

.

The plot of the first 100 terms of the series is given in the following animation.

The overshoting oscillations in correspondence of the discontinuity of the square wave is called *Gibbs phenomenon *by* W. Gibbs that *pointed it out in a letter to Nature in 1899. The oscillation overshoot the function of a fixed amount that can be calculated. In addition, the function oscillate with an amplitude that is decreating going far from the discontinuity point.

### Fourier Series using Complex Notation

Let be a piecewise linear function on (This means that may possess a finite number of finite discontinuities on the interval). Then can be expanded in a

Using Euler’s identities:

where ,

the Fourier series of can be written in a more compact from as

where

### The Discrete Fourier Transform (DFT)

IN CONSTRUCTION

## Appendix I

*A program in Awk language to calculate the DFT.*

This program read a file containing a list of values representing a function or a signal or a generic set of data to analyze using the harmonical analysis. Then it calculates the Fourier coefficients and finally reconstructs the Fourier series using a given number of terms (n).

```
[code language="cpp"]
#==========================================================
#
# NAME: dft.awk
#
#==========================================================
# DESCRIPTION: Calculate the fourier series of a signal
# INPUT: The program read a two column file containing the
# signal to elaborate
# OUTPUT: The output contain the reconstructed signal
# using the first n Fourier components.
#==========================================================
# DEVELOPED USING: gawk
# FUNCTIONS:
# ReadData() : read a two column data file from stdin
# CalculateDFT(): Calculate the DFT.
# Reconstruct() : Reconstruct/filter the signal using the
# first n components of the DFT, and print on the
# stdout thereconstructed signal.
#==========================================================
# Author: Danilo Roccatano
# Version: 1.0
# Copyright (C): 2016 Danilo Roccatano
#==========================================================
#
#
function ReadData() {
#
# Read data from standard input
#
np=0
while (getline > 0 ) {
if (NF >0) {
P[np]=$2
np++
}
}
np--
}
function CalculateDFT() {
#
# Perfrom the discrete fourier transform of the data
#
PI=4*atan2(1,1)
for (j=1;j<=n;j++) {
A[j]=0
B[j]=0
for (i=0;i<=np-1;i++) {
x=2*PI*i*j/np
A[j]+=P[i]*cos(x)
B[j]+=P[i]*sin(x)
}
A[j]/=(np/2)
B[j]/=(np/2)
}
}
function CalcA0(){
#
# Calculate the average value of A0
#
A0=0
for (i=0;i<=np;i++) {
A0+=P[i]
}
A0=A0/np
}
function Reconstruct() {
for (i=0;i<=np;i++){
M[i]=0
for (j=1;j<=np;j++) {
F=2*PI*j*i/np
M[i]=M[i]+A[j]*cos(F)+B[j]*sin(F)
}
print i,M[i]+A0
}
}
BEGIN {
#
# n:= number of armonics to consider
#
n=30
ReadData()
CalcA0()
CalculateDFT()
#
# Reconstruct the signal using n harmonics
#
Reconstruct()
}
[/code]
```

### Appendix II

*A program in tcl/tk language for the harmonic analysis of a drawn signal*

This program was inspired to a similar program in BASIC language in a book that I read many years ago [1]. Using the mouse, you can draw a function in the upper canvas, calculate numerically the Fourier coefficients and plot a number of terms to your choice of the discrete Fourier series.

```
[code language="bash"]
#! /bin/sh
# the next line restarts using wish \
exec wish "$0" "$@"
#==========================================================
# NAME: FourierV3.0.tcl
#==========================================================
# DESCRIPTION: Perfrom the harmonic analysis of a drawn
# signal.
#==========================================================
# DEVELOPED USING: tcl/tk
#==========================================================
# Author: Danilo Roccatano (c) 2018
# Version: 3.0
#==========================================================
#
#
package require Tk
package require xyplot
package require Plotchart
wm title . ""
array set Yv {}
array set Yvt {}
array set A {}
array set B {}
array set P {}
set width 1000
set height 200
set midX [expr { $width/2 }]
set midY [expr { $height/ 2 }]
set imidY [expr { 1.0/ $midY }]
set xsta 0
set xend 0
set x0 0
set y0 0
set np $width
set np1 [expr $np-1]
set np2 [expr $np/2]
set inp2 [expr 1.0/$np2 ]
set n 10
set A0 0
set pi 3.1415926535897931
set K 2*$pi/$np
set COLORS {
navy
cyan
green chartreuse
gold
chocolate
red
magenta
purple
blue
turquoise
aquamarine
yellow
sienna
peru
burlywood
wheat
tan
firebrick
brown
salmon
orange
coral
tomato
pink
maroon
violet
orchid
lavender
}
set nCOLORS [llength $COLORS]
set icol 0
puts "$nCOLORS"
global xsta xend width height midY K inp2 imidY A0 COLORS icol
for { set i 1 } { $i <= $width+10 } { incr i 1 } {
set Yv($i) [expr { $midY } ]
set P($i) 0
set Yvt($i) 0
set A($i) 0
set B($i) 0
}
proc DrawAxis {w} {
global width height midY
$w create line 0 $midY $width $midY -tags "Xaxis" -fill black -width 1
}
proc redraw {w ar1 ar2 ar3 x y} {
upvar $ar1 Yv
upvar $ar2 Yvt
upvar $ar3 P
global xsta xend width K midY imidY A0 np
set x0 0
set y0 0
set icol 0
$w delete "all"
DrawAxis .c
set A0 0
for { set x $xsta } { $x <= $xend } { incr x 1 } {
set ix [expr int($x)]
set Yv($ix) 0
set Yv($ix) $Yvt($x)
set Yvt($ix) 0
}
set A0 [expr $A0/$np]
set y0 0
for { set x 1 } { $x <= $width } { incr x 1 } {
set ix [expr int($x)]
if {$Yv($ix) != 0 } {
$w create line $x0 $y0 $x $Yv($ix) -fill red -width 1
set x0 $x
set y0 $Yv($ix)
}
if {$Yv($ix) == 0 } {
set Yv($ix) $y0
}
set P($ix) [expr -($Yv($ix)-$midY)*$imidY]
}
}
proc down {w Yv Yvt x y} {
global xsta x0 y0
set ::ID [$w create line $x $y $x $y -fill black]
set xsta $x
set x0 $x
set y0 $y
}
proc move {w ar x y} {
upvar $ar Yvt
global xend xsta x0 y0
set x [$w canvasx $x]
set y [$w canvasx $y]
$w create line $x0 $y0 $x $y -fill black
set xend $x
set ix [expr int($x)]
set Yvt($ix) $y
set x0 $x
set y0 $y
}
proc DFT {w w1 w2 n ar1 ar2 ar3 ar4} {
#
# Perfrom the discrete fourier transform of the data
#
upvar $ar1 A
upvar $ar2 B
upvar $ar3 Yv
upvar $ar4 P
global pi np np1 np2 inp2 K midY A0 icol
puts "Calculating the coefficients ... "
for { set j 1} {$j <= $n} {incr j 1} {
set A($j) 0
set B($j) 0
for { set i 1} {$i <= $np} {incr i 1} {
set x [expr $K*($i-1)*$j]
set A($j) [expr $A($j)+$P($i)*cos($x)]
set B($j) [expr $B($j)+$P($i)*sin($x)]
}
set A($j) [expr $A($j)*$inp2]
set B($j) [expr $B($j)*$inp2]
}
PlotPSpect $w1 $w2 $n A B
puts "Reconstructing the signal using $n harmonics"
set icol [expr $icol +1]
if {$icol > 29} {set icol 0}
RecSignal $w $n A B
}
proc RecSignal {w n ar1 ar2 } {
upvar $ar1 A
upvar $ar2 B
set x0 0
set y0 0
global pi np midY K A0 COLORS icol
for { set i 1} {$i <= [expr $np+1]} {incr i 1} {
set M($i) 0
for { set j 1} {$j <= $n} {incr j 1} {
set x [expr $K*($i-1)*$j]
set M($i) [expr $M($i)+$A($j)*cos($x)+$B($j)*sin($x)]
}
set M($i) [expr $midY-$M($i)*$midY]
$w create line $x0 $y0 $i $M($i) -fill [lindex $COLORS $icol] -width 1
set x0 $i
set y0 $M($i)
}
}
proc ClrCanvas {w } {
$w delete "all"
set icol 0
DrawAxis $w
}
proc PlotPSpect {w1 w2 n ar1 ar2} {
upvar $ar1 A
upvar $ar2 B
set xval {}
global pi
$w1 delete all
$w2 delete all
set p [::Plotchart::createHistogram $w1 {0 30 2} {0 1.0 0.1}]
set p1 [::Plotchart::createXYPlot $w2 {0 30 2} {-180. 180 30 }]
$p title "Power Spectrum"
$p xtext "# harmonics"
$p dataconfig data -fillcolour cyan -width 2 -colour blue
$p1 dataconfig data -fillcolour cyan -width 2 -colour red -type symbol
$p1 title "Phase angle"
$p1 xtext "# harmonics"
for { set i 1 } { $i <= $n } { incr i } {
$p plot data $i [expr sqrt($A($i)*$A($i)+$B($i)*$B($i))]
$p1 plot data $i [expr -atan2($B($i),$A($i))*180./$pi]
puts "$i [expr -atan2($B($i),$A($i))]"
# puts "$A($i) $B($i)"
}
}
#
#-- Building the UI
set colors { black red orange }
label .ltitle \
-background magenta -padx 64 -relief raised -text {Harmonic Analysis (c) Danilo Roccatano 2002-2018}
set xyp [xyplot .xyp -xformat "%5.0f" -yformat "%5.4f" -title "Power Spectrum" -background gray90]
button .b1 -text FOURIER -command { DFT .c .c1 .c2 $n A B Yv P}
scale .scale -orient horizontal -length 300 -from 1 -to 30 \
-tickinterval 3 -variable n
label .lscale -text " Harmonics"
label .descr -justify left -background pink -padx 32 -relief raised -wraplength 800 -text {Calculation of the discrete Fourier Harmonics of a generic curve drawn in the above canvas. Use the mouse pointer to draw a curve. Everytime you click on the canvas the new curve will be merged with the existent one. Once you you finish, select the number f harmonics to consider for the curve reconstruction (defauls 10) and then press FOURIER button. The program calculate also the power spectrum and the phase angle for each harmonic term of the series. HAVE FUN!}
button .b2 -text "Clear All" -padx 20 -pady 20 -command { ClrCanvas .c }
canvas .c -relief raised -borderwidth 1 -width 1000 -height 200 \
-background white -relief raised -border 2
pack [canvas .c1 -width 400 -height 300 -bg white]
pack [canvas .c2 -width 400 -height 300 -bg white]
grid .ltitle -row 0 -columnspan 4 -sticky news
grid .descr -row 1 -columnspan 4 -sticky news
grid .c -row 2 -columnspan 4 -sticky news
grid .c1 -row 3 -columnspan 2 -columnspan 2 -sticky news
grid .c2 -row 3 -column 3 -columnspan 2 -sticky news
grid .b1 -row 4 -column 0 -sticky w
grid .b2 -row 4 -column 1
grid .lscale -row 4 -column 3 -sticky e
grid .scale -row 4 -column 3 -sticky nw
#-- The current mode is retrieved at runtime from the global ''Mode'' variable:
bind .c <1> {down %W Yv Yvt %x %y}
bind .c <ButtonRelease-1> {redraw %W Yv Yvt P %x %y }
bind .c <B1-Motion> {move %W Yvt %x %y}
bind .c <3> {%W delete current}
DrawAxis .c
[/code]
```