<< Chapter < Page | Chapter >> Page > |
In this lab, we learn how to compute the continuous-time Fourier transform (CTFT), normally referred to as Fourier transform, numerically and examine its properties. Also, we explore noise cancellation and amplitude modulation as applications of Fourier transform.
The continuous-time Fourier transform (CTFT) (commonly known as Fourier transform) of an aperiodic signal $x(t)$ is given by
The signal $x(t)$ can be recovered from $X(\omega )$ via this inverse transform equation
Some of the properties of CTFT are listed in [link] .
Properties | Time domain | Frequency domain |
Time shift | $x(t-{t}_{0})$ | $X(\omega ){e}^{-{\mathrm{j\omega t}}_{0}}$ |
Time scaling | $x(\text{at})$ | $\frac{1}{\mid a\mid}X\left(\frac{\omega}{a}\right)$ |
Linearity | ${a}_{1}{x}_{1}(t)+{a}_{2}{x}_{2}(t)$ | ${a}_{1}{X}_{1}(\omega )+{a}_{2}{X}_{2}(\omega )$ |
Time convolution | $x(t)\ast h(t)$ | $X(\omega )H(\omega )$ |
Frequency convolution | $x(t)h(t)$ | $X(\omega )\ast H(\omega )$ |
Refer to signals and systems textbooks [link] - [link] for more theoretical details on this transform.
Assuming that the signal $x(t)$ is zero for $t\text{<0}$ and $t\ge T$ , we can approximate the CTFT integration in Equation (1) as follows:
where
$T=\mathrm{N\tau}$ and N is an integer. For sufficiently small
$\tau $ , the above summation provides a close approximation to the CTFT integral. The summation
$\sum _{n=0}^{N-1}x(\mathrm{n\tau}){e}^{-\mathrm{j\omega n\tau}}$ is widely used in digital signal processing (DSP), and both LabVIEW MathScript and LabVIEW have a built-in function for it called
fft
. In a .m file, if N samples
$x(\mathrm{n\tau})$ are stored in a vector
$x$ , then the function call
>>xw=tau*fft
(x)
calculates
where
with N assumed to be even. The
fft
function returns the positive frequency samples before the negative frequency samples. To place the frequency samples in the right order, use the function
fftshift
as indicated below:
>>xw=fftshift(tau*fft
(x
)
)
Note that
$X(\omega )$ is a vector (actually, a complex vector) of dimension N.
$X(\omega )$ is complex in general despite the fact that
$x(t)$ is real-valued. The magnitude of
$X(\omega )$ can be computed using the function
abs
and the phase of
$X(\omega )$ using the function
angle
.
Notification Switch
Would you like to follow the 'An interactive approach to signals and systems laboratory' conversation and receive update notifications?