In this module, we will derive an expansion for arbitrary discrete-time functions, and in doing so, derive the Discrete Time Fourier Transform (DTFT).

Since complex exponentials are eigenfunctions of linear time-invariant (LTI) systems, calculating the output of an LTI system $\mathscr{H}$ given ${e}^{i\omega n}$ as an input amounts to simple multiplication, where ${\omega}_{0}=\frac{2\pi k}{N}$ , and where $H\left[k\right]\in \mathbb{C}$ is the eigenvalue corresponding to k. As shown in the figure, a simple exponential input would yield the output

$$y\left[n\right]=H\left[k\right]{e}^{i\omega n}$$

1

Using this and the fact that $\mathscr{H}$ is linear, calculating $y\left[n\right]$ for combinations of complex exponentials is also straightforward.

$${c}_{1}{e}^{i{\omega}_{1}n}+{c}_{2}{e}^{i{\omega}_{2}n}\to {c}_{1}H\left[{k}_{1}\right]{e}^{i{\omega}_{1}n}+{c}_{2}H\left[{k}_{2}\right]{e}^{i{\omega}_{1}n}$$ $$\sum _{l}^{}{c}_{l}{e}^{i{\omega}_{\mathrm{l}}n}\to \sum _{l}^{}{c}_{l}H\left[{k}_{l}\right]{e}^{i{\omega}_{\mathrm{l}}n}$$

The action of $H$ on an input such
as those in the two equations above is easy to explain.
** $\mathscr{H}$ independently
scales** each exponential component
${e}^{i{\omega}_{\mathrm{l}}n}$
by a different complex number
$H\left[{k}_{l}\right]\in \mathbb{C}$. As such, if we can write a function
$y\left[n\right]$
as a combination of complex exponentials it allows us to easily calculate the output of a system.

Now, we will look to use the power of complex exponentials to see how we may represent arbitrary signals in terms of a set of simpler functions by superposition of a number of complex exponentials. Below we will present the Discrete-Time Fourier Transform (DTFT). Because the DTFT deals with nonperiodic signals, we must find a way to include all real frequencies in the general equations. For the DTFT we simply utilize summation over all real numbers rather than summation over integers in order to express the aperiodic signals.

It can be demonstrated that an arbitrary Discrete Time-periodic function $f\left[n\right]$ can be written as a linear combination of harmonic complex sinusoids

$$f\left[n\right]=\sum _{k=0}^{N-1}{c}_{k}{e}^{i{\omega}_{0}kn}$$

2

The
${c}_{n}$ - called the Fourier coefficients -
tell us "how much" of the sinusoid
${e}^{j{\omega}_{0}kn}$ is in
$f\left[n\right]$.
The formula shows
$f\left[n\right]$ as a sum of complex exponentials, each of which is easily processed by an
LTI system (since it is an eigenfunction of
**every** LTI system). Mathematically,
it tells us that the set of
complex exponentials
$\left\{\forall k,k\in \mathbb{Z}:\left({e}^{j{\omega}_{0}kn}\right)\right\}$ form a basis for the space of N-periodic discrete
time functions.

Now, in order to take this useful tool and apply it to arbitrary non-periodic signals, we will have to delve deeper into the use of the superposition principle. Let ${s}_{T}\left(t\right)$ be a periodic signal having period $T$. We want to consider what happens to this signal's spectrum as the period goes to infinity. We denote the spectrum for any assumed value of the period by ${c}_{n}\left(T\right)$. We calculate the spectrum according to the Fourier formula for a periodic signal, known as the Fourier Series (for more on this derivation, see the section on Fourier Series.)

$${c}_{n}=\frac{1}{T}{\int}_{0}^{T}\phantom{\rule{-0.166667em}{0ex}}s\left(t\right)exp(-\u0131{\omega}_{0}t)\phantom{\rule{0.166667em}{0ex}}dt$$

3

$${S}_{T}\left(f\right)\equiv T{c}_{n}=\frac{1}{T}{\int}_{0}^{T}\phantom{\rule{-0.166667em}{0ex}}({S}_{T}\left(f\right)exp(\u0131{\omega}_{0}t)\phantom{\rule{0.166667em}{0ex}}dt$$

4

$${s}_{T}\left(t\right)=\sum _{-\infty}^{\infty}\phantom{\rule{-0.166667em}{0ex}}f\left(t\right)exp(\u0131{\omega}_{0}t)\phantom{\rule{0.166667em}{0ex}}\frac{1}{T})$$

5

$$\underset{T\to \infty}{lim}{s}_{T}\left(t\right)\equiv s\left(t\right)=\underset{-\infty}{\overset{\infty}{\int}}\phantom{\rule{-0.166667em}{0ex}}S\left(f\right)exp(\u0131{\omega}_{0}t)\phantom{\rule{0.166667em}{0ex}}df$$

6

$$S\left(f\right)=\underset{-\infty}{\overset{\infty}{\int}}\phantom{\rule{-0.166667em}{0ex}}s\left(t\right)exp(-\u0131{\omega}_{0}t)\phantom{\rule{0.166667em}{0ex}}dt$$

7

Discrete-Time Fourier Transform

$$\mathcal{F}\left[\Omega \right]=\sum _{n=-\infty}^{\infty}f\left[n\right]{e}^{-(i\Omega n)}$$

8

Inverse DTFT

$$f\left(n\right)=\frac{1}{2\pi}{\int}_{0}^{2\pi}\mathcal{F}\left(\Omega \right){e}^{i\Omega n}d\Omega $$

9

Warning:

It is not uncommon to see
the above formula written slightly different. One of the
most common differences is the way
that the exponential is written. The above equations use the radial
frequency variable $\Omega $ in the exponential, where
$\Omega =2\pi f$, but it is also common to include the more explicit expression,
$i2\pi ft$, in the exponential.
Click here for an
overview of the notation used in Connexion's DSP modules.
Because complex exponentials are eigenfunctions of LTI systems, it is often useful to represent signals using a set of complex exponentials as a basis. The discrete time Fourier transform synthesis formula expresses a discrete time, aperiodic function as the infinite sum of continuous frequency complex exponentials.

$$\mathcal{F}\left[\Omega \right]=\sum _{n=-\infty}^{\infty}f\left[n\right]{e}^{-(i\Omega n)}$$

10

$$f\left(n\right)=\frac{1}{2\pi}{\int}_{0}^{2\pi}\mathcal{F}\left(\Omega \right){e}^{i\Omega n}d\Omega $$

11