您在這裡

Derivation of the Fourier Transform

2 六月, 2016 - 15:17

Fourier series clearly open the frequency domain as an interesting and useful way of determining how circuits and systems respond to periodic input signals. Can we use similar techniques for nonperiodic signals? What is the response of the filter to a single pulse? Addressing these issues requires us to find the Fourier spectrum of all signals, both periodic and nonperiodic ones. We need a Definition for the Fourier spectrum of a signal, periodic or not. This spectrum is calculated by what is known as the Fourier transform.

Let sT (t) be a periodic signal having period T. We want to consider what happens to this signal's spectrum as we let the period become longer and longer. We denote the spectrum for any assumed value of the period by ck (T). We calculate the spectrum according to the familiar formula

c_k(T)=\frac{1}{T}\int_{-\left ( \frac{T}{2} \right )}^{\frac{T}{2}}s_T(t)e^{-\left ( \frac{j2\pi kt}{T} \right )}dt

where we have used a symmetric placement of the integration interval about the origin for subsequent deriva tional convenience. Let f be a fxed frequency equaling Tk ; we vary the frequency index k proportionally as we increase the period. Define

S_T(f)\equiv Tc_k(T)=\int_{-\left ( \frac{T}{2} \right )}^{\frac{T}{2}}s_T(t)e^{-\left ( {j2\pi ft} \right )}dt

making the corresponding Fourier series

s_T(t)=\sum_{k=-\infty }^{\infty }\left ( S_T(f)e^{j2\pi ft}\frac{1}{T} \right )

As the period increases, the spectral lines become closer together, becoming a continuum. Therefore,

media/image472.png

with

S(f)=\int_{-\infty }^{\infty }s(t)e^{-(j2\pi ft)}dt

S (f) is the Fourier transform of s (t) (the Fourier transform is symbolically denoted by the uppercase version of the signal's symbol) and is defined for any signal for which the integral ((4.33)) converges.

Example 4.4

Let's calculate the Fourier transform of the pulse signal (Section 2.2.5: Pulse), p(t).

P(f)=\int_{-\infty }^{\infty }p(t)e^{-(j2\pi ft)}dt=\int_{0}^{\triangle}e^{-(j2\pi ft)}dt=\frac{1}{-(j2\pi f)}(e^{-(j2\pif\triangle)}-1)

P(f)=e^{-(j\pi f\triangle)}\frac{sin(\pi f\triangle)}{\pi f}

Note how closely this result resembles the expression for Fourier series coefficients of the Figure 4.10.

media/image476.png
Figure 4.11 Spectrum
The upper plot shows the magnitude of the Fourier series spectrum for the case of T =1 with the Fourier transform of p(t) shown as a dashed line. For the bottom panel, we expanded the period to T =5, keeping the pulse's duration fixed at 0.2, and computed its Fourier series coefficients.  

Figure 4.11 (Spectrum) shows how increasing the period does indeed lead to a continuum of coefficients, and that the Fourier transform does correspond to what the continuum becomes. The quantity \frac{sin(t)}{t}has a special name, the sinc (pronounced "sink") function, and is denoted by sinc (t). Thus, the magnitude of the pulse's Fourier transform equals |ΔsincfΔ) |.

The Fourier transform relates a signal's time and frequency domain representations to each other. The direct Fourier transform (or simply the Fourier transform) calculates a signal's frequency domain representation from its time-domain variant ((4.34)). The inverse Fourier transform finds the time-domain representation from the frequency domain. Rather than explicitly writing the required integral, we often symbolically express these transform calculations as F (s) and F1 (S), respectively.

\begin{align*} F(s)&=S(f)\\ &=\int_{-\infty }^{\infty }s(t)e^{-(j2\pi ft)}dt\\ \end{align*}    

\begin{align*} F^{-1}(s)&=s(f)\\ &=\int_{-\infty }^{\infty }S(t)e^{+j2\pi ft}df\\ \end{align*}

We must have s (t)= F1 (F (s (t))) and S (f)= F(F1 (S (f)), and these results are indeed valid with minor exceptions.

Note: Recall that the Fourier series for a square wave gives a value for the signal at the discontinuities equal to the average value of the jump. This value may difer from how the signal is defined in the time domain, but being unequal at a point is indeed minor.

Showing that you "get back to where you started" is difficult from an analytic viewpoint, and we won't try here. Note that the direct and inverse transforms differ only in the sign of the exponent.

Exercise 4.8.1

The differing exponent signs means that some curious results occur when we use the wrong sign. What is F (S (f))? In other words, use the wrong exponent sign in evaluating the inverse Fourier transform.

Properties of the Fourier transform and some useful transform pairs are provided in the accompanying tables (Table 4.1 and Table 4.2). Especially important among these properties is Parseval's Theorem, which states that power computed in either domain equals the power in the other.

\int_{-\infty }^{\infty }s^2(t)dt=\int_{-\infty }^{\infty }(|S(f)|)^2df

Of practical importance is the conjugate symmetry property: When s (t) is real-valued, the spectrum at negative frequencies equals the complex conjugate of the spectrum at the corresponding positive frequencies. Consequently, we need only plot the positive frequency portion of the spectrum (we can easily determine the remainder of the spectrum).

Exercise 4.8.2

How many Fourier transform operations need to be applied to get the original signal back:
F(...(F(s)))=s(t)?

Note that the mathematical relationships between the time domain and frequency domain versions of the same signal are termed transforms. We are transforming (in the nontechnical meaning of the word) a signal from one representation to another. We express Fourier transform pairs as (s (t) ↔ S (f)). A signal's time and frequency domain representations are uniquely related to each other. A signal thus "exists" in both the time and frequency domains, with the Fourier transform bridging between the two. We can define an information carrying signal in either the time or frequency domains; it behooves the wise engineer to use the simpler of the two.

A common misunderstanding is that while a signal exists in both the time and frequency domains, a single formula expressing a signal must contain only time or frequency: Both cannot be present simultaneously. This situation mirrors what happens with complex amplitudes in circuits: As we reveal how communications systems work and are designed, we will define signals entirely in the frequency domain without explicitly finding their time domain variants. This idea is shown in another module (Section 4.6) where we define Fourier series coefficients according to letter to be transmitted. Thus, a signal, though most familiarly defined in the time-domain, really can be defined equally as well (and sometimes more easily) in the frequency domain. For example, impedances depend on frequency and the time variable cannot appear.

We will learn (Section 4.9) that finding a linear, time-invariant system's output in the time domain can be most easily calculated by determining the input signal's spectrum, performing a simple calculation in the frequency domain, and inverse transforming the result. Furthermore, understanding communications and information processing systems requires a thorough understanding of signal structure and of how systems work in both the time and frequency domains.

The only difficulty in calculating the Fourier transform of any signal occurs when we have periodic signals (in either domain). Realizing that the Fourier series is a special case of the Fourier transform, we simply calculate the Fourier series coefficients instead, and plot them along with the spectra of nonperiodic signals on the same frequency axis.

Table 4.1 Short Table of Fourier Transform Pairs
s(t) S(f)
e^{-(at)}u\left ( t \right ) \frac{1}{j2\pi f+a}
e^{\left ( -a \right )\left | t \right |} \frac{2a}{4\pi ^{2}f^{2}+a^{2}}
p(t)=\left\{\begin{matrix} 1\if\left | t \right | < \frac{\Delta }{2}& & \\ 0\ if \left | t \right | > \frac{\Delta }{2}& & \end{matrix}\right. \frac{sin\left ( \pi f\Delta \right )}{\pi f}
\frac{sin\left ( 2\pi Wt \right )}{\pi t} S(f)=\left\{\begin{matrix} 1\if\left | f \right | < W& & \\ 0\if\left | f \right |> W& & \end{matrix}\right.
 
Table 4.2 Fourier Transform Properties
  Time-Domain Frequency Domain
Linearity a_{1}s_{1}\left ( t \right )+a_{2}s_{2}\left ( t \right ) a_{1}S_{1}\left ( f \right )+a_{2}S_{2}\left ( f \right )
Conjugate Symmetry s\left ( t \right )\in \mathbb{R} S\left ( f \right )=S\left ( -f \right )^{*}
Even Symmetry s\left ( t \right )=s\left ( -t \right ) S\left ( f \right )=S\left ( -f\right )
Odd Symmetry s\left ( t \right )=-\left ( s\left ( -t \right ) \right ) S\left ( f \right )=-\left ( S\left ( -f \right ) \right )
Scale Change s\left ( at \right ) \frac{1}{\left | a \right |}S\left ( \frac{f}{a} \right )
Time Delay s\left ( t-\tau \right ) e^{-\left ( j2\pi f\tau \right ) }S\left ( f \right )
Complex Modulation e^{j2\pi f_{0}t}s\left ( t \right ) S\left ( f-f_{0} \right )
Amplitude Modulation by Cosine s\left ( t \right )cos\left ( 2\pi f_{0}t \right ) \frac{S\left ( f-f_{0} \right )+S\left ( f+f_{0} \right )}{2}
Amplitude Modulation by Sine s\left ( t \right )sin\left ( 2\pi f_{0}t \right ) \frac{S\left ( f-f_{0} \right )-S\left ( f+f_{0} \right )}{2j}
Differentiation \frac{d}{dt}s\left ( t \right ) j2\pi fS\left ( f \right )
Multiplication by t \int_{-\infty }^{t}s\left ( \alpha \right )d\alpha \frac{1}{j2\pi f}S\left ( f \right )\if\ S\left ( 0 \right )=0
Area ts\left ( t \right ) \frac{1}{-\left ( j2\pi \right )}\frac{d}{df}S\left ( f \right )
Value at Origin \int_{-\infty }^{\infty }s\left ( t \right )dt S\left ( 0 \right )
Parseval"s Theorem \int_{-\infty }^{\infty }\left ( \left | s\left ( t \right ) \right | \right )^{2}dt \int_{-\infty }^{\infty }\left ( \left | S\left ( f \right ) \right | \right )^{2}df
 

Example 4.5

In communications, a very important operation on a signal s (t) is to amplitude modulate it. Using this operation more as an example rather than elaborating the communications aspects here, we want to compute the Fourier transform the spectrum of
(1+s(t))cos(2\pi f_ct)
Thus,
(1+s(t))cos(2\pi f_c t)=cos(2\pi f_ct)+s(t)cos(2\pi f_ct)

For the spectrum of cos (2πfct), we use the Fourier series. Its period is \frac{1}{f_c}, and its only nonzero Fourier coefficients are c\pm 1=\frac{1}{2}. The second term is not periodic unless s (t) has the same period as the sinusoid. Using Euler's relation, the spectrum of the second term can be derived as

s(t)cos(2\pi f_c t)=\int_{-\infty }^{\infty }S(f)e^{j2\pi ft}df\cos(2\pi f_ct)

Using Euler's relation for the cosine,

\begin{align*} \left ( s(t)cos(2\pi f_c t) \right )&=\frac{1}{2}\int_{-\infty }^{\infty }S(f)e^{j2\pi (f+f_c)^t}df+\frac{1}{2}\int_{-\infty }^{\infty }S(f)e^{j2\pi (f+f_c)^t}df\\ \left ( s(t)cos(2\pi f_c t) \right )&=\frac{1}{2}\int_{-\infty }^{\infty }S(f-f_c)e^{j2\pi ft}df+\frac{1}{2}\int_{-\infty }^{\infty }S(f+f_c)e^{j2\pi ft}df\\ (s(t)cos&(2\pi f_c t))=\int_{-\infty }^{\infty }\frac{S(f-f_c)+S(f+f_c)}{2}e^{j2\pi ft}df\\ \end{align*}

Exploiting the uniqueness property of the Fourier transform, we have

F(s(t)cos(2\pi f_c t))=\frac{S(f-f_c)+S(f+f_c)}{2}

This component of the spectrum consists of the original signal's spectrum delayed and advanced in frequency. The spectrum of the amplitude modulated signal is shown in Figure 4.12.

media/image492.png
Figure 4.12 Spectrum of the amplitude modulated signal
A signal which has a triangular shaped spectrum is shown in the top plot. Its highest frequency the largest frequency containing power is W Hz. Once amplitude modulated, the resulting spectrum has "lines" corresponding to the Fourier series components at ± (fc) and the original triangular spectrum shifted to components at ± (fc) and scaled by .  

Note how in this figure the signal s (t) is defined in the frequency domain. To find its time domain representation, we simply use the inverse Fourier transform.

Exercise 4.8.3

What is the signal s (t) that corresponds to the spectrum shown in the upper panel of Figure 4.12?

Exercise 4.8.4

What is the power in x (t), the amplitude-modulated signal? Try the calculation in both the time and frequency domains.

In this example, we call the signal s (t) a baseband signal because its power is contained at low frequencies. Signals such as speech and the Dow Jones averages are baseband signals. The baseband signal's bandwidth equals W , the highest frequency at which it has power. Since x (t)'s spectrum is confned to a frequency band not close to the origin (we assume fc» W ) , we have a bandpass signal. The bandwidth of a bandpasssignal is not its highest frequency, but the range of positive frequencies where the signal has power. Thus, in this example, the bandwidth is 2W Hz. Why a signal's bandwidth should depend on its spectral shape will become clear once we develop communications systems.