Lab 11-Fourier Transform
Lab 11-Fourier Transform
[Spring 2012]
In other words, the Fourier transform of a signal x(t ) is a signal X () . An alternative way of
writing 11.1 is (11.2)
F
x (t )
X ( )
X ( ) F {x(t )}
x (t )e
j t
dt
From 11.3 it is clear that X () is a complex function of . In case the Fourier transform of x(t ) has
to be expressed in the frequency domain f , then substituting by 2 f in 11.3, yields 11.4 given
as
X ( f ) F {x(t )}
x (t )e
j 2 t
dt
In order to return from the frequency domain back to the time domain the inverse Fourier
transform is implemented. The inverse Fourier transform is denoted by the symbol F 1{} ; that is,
one can write (11.5) as
x(t ) F 1{ X ( )}
or alternatively (11.6) as
1
[FOURIER TRANSFORM]
[Spring 2012]
F
X ()
x (t )
X ( )e
j t
x(t ) F 1{ X ( f )}
X ( f )e
j 2 ft
df
The cyclic frequency is measured in rad/s, while the frequency f is measured in Hertz. The
Fourier transform of a signal is called (frequency) spectrum.
Results
ans=pi^(1/2)*exp(1/4*w^2)
ans= exp(1/4*w^2)*pi^(1/2)
Comments
The Fourier transform of
2 4
.
x(t ) is X () e
The result is verified
according to 11.3.
Example:
Compute the inverse Fourier transform of the function X () 1 (1 j)
Commands
Comments
X=1/(1+j*w);
ans=exp(-x)*heaviside(x)
The inverse Fourier transform
Ifourier(X)
of the function is expressed
with x as the independent
variable. The result is e x u ( x) .
However, it is more appropriate to use t as independent variable . This is accomplished by using
the syntax ifourier(X,t)
2
Results
[FOURIER TRANSFORM]
[Spring 2012]
Commands
Results
X=1/(1+j*w);
Ifourier(X,t)
ans=exp(-t)*heaviside(t)
Comments
The inverse Fourier transform
of the function is expressed
with t as the independent
variable. The result is e t u (t ) .
The command fourier should be executed as fourier(x,w). Using this syntax is optimal
since it allows the computation of constant functions.
Commands
Results
Comments
?? Function fourier is not defined for The simple syntax is not
values of class double
appropriate for computing the
Fourier transform of constant
functions.
ans=2*pi*dirac(w)
Using complex syntax allows
us to compute the Fourier
transform (if it exits) of any
function.
ans=2*pi*dirac(s)
Changing the variable in which
the result is given.
ans=exp(-n)*heaviside(n)
A variable change is also
possible
at
ifourier
command.
x=1;
fourier(x)
fourier(x,w)
syms s
fourier(x,s)
syms n
X=1/(1+j*w);
Ifourier(X,n)
Frequency Domain
Commands
X ( )
1
syms t w w0 t0
1
u (t )
2 (t )
(1 j) (t )
(t t 0 )
e j t0
x=dirac(t);
fourier(x,w)
fourier(1,w)
X=1/(j*w)+pi*dirac(w);
ifourier(X,t)
X=dirac(t-t0);
fourier(x,w)
Results
ans=1
ans=2*pi*dirac(w)
ans=heaviside(t)
ans=exp(-i*t0*w)
[FOURIER TRANSFORM]
[Spring 2012]
e j 0t
2 ( 0 )
cos( 0t )
( 0 ) ( 0 )
sin(0t )
j ( 0 ) j ( 0 )
e at u (t )
(a) 0
1 ( j a )
te at u (t )
1 ( j a ) 2
(a) 0
(t n 1 (n 1)!)e at u (t )1 ( j a)2
(a) 0
X=2*pi*dirac(w-wo);
ifourier(X,t)
X=pi*dirac(w-w0)+dirac(w+w0);
ifourier(X,t)
X=(pi/j)*dirac(w-w0)-dirac(w+w0);
x=ifourier(X,t)
a=8;
x=exp(-a*t)*heaviside(t);
X=fourier(x,w)
x=(t)*exp(-a*t)*heaviside(t);
fourier(x,w)
ans=exp(i*t0*w)
n=4;
X=1/(j*w+a)^n;
Ifourier(X,t)
ans=1/6*t^3*exp(8*t)*heaviside(t)
ans=cos(w0*t)
ans=sin(w0*t)
X=1/(8+i*w)
ans=1/(8+i*w)^2
Tasks
Task 01: Plot the Fourier transform of signal x (t ) cos(t ) .
Task 02: Plot the Fourier transform of signal x(t ) sin( t ) / ( t ) .
Task 03: Plot the inverse Fourier transform of signal X () sin / .
Task 04: Suppose that a signal x(t ) is given by x (t ) te3t . Compute the Fourier transform X () of
the signal x(t ) and plot for 20 20 rad/s.
a. The magnitude of X ()
b. The angle of X ()
c. The real part of X ()
d. The imaginary part of X ()