Our goal is to write Fourier series for aperiodic signals. Here's our prototypical example.
f(t):={10 if t∈(−1,1). otherwise.
The guiding intuition is that all signals are technically periodic—it's just that the period of an aperiodic signal is T→∞. Instead of a discrete sum∑k=−∞∞F[k]ejktω0, we'll write f as a continuous integral∫−∞∞F(ω)ejtωdω.
But let's take things slow. To build up to being able to say that f(t) has period T→∞, let's start by considering what happens when we force f(t) to have some finite period T by constructing an alternative fT(t).
fT(t):={10 if (t%T)∈(−1,1). otherwise.
As a function of T, we now determine the array of Fourier coefficients FT for the periodic signal fT(t).
And here's a plot of T⋅FT[k]=(2πk/T)2sin(2πk/T) for T=16 and T=8. (Note that k can be any integer, not just those pictured below.)
Remark. Convince yourself that doubling T has the effect of halving the “height” and doubling the “width” of the raw coefficients FT[k]. If we wanted to “normalize” the graphs for T=16 and T=8 so they looked similar, we might instead graph, say, T⋅FT[k] against 2πk/T…
And so our Fourier series for fT(t) looks like the following:
fT(t)=T1k=−∞∑∞2πk/T2sin(2πk/T)ejt(2πk/T).
There's a reparametrization begging to be made here: ωk:=2πk/T for each k∈Z. Note that Δω:=ωk+1−ωk=2π/T. Now we can send T→∞ and see what happens to the Fourier series.
And so our Fourier series now becomes a Fourier integral!
§ Lecture: The CT Fourier Transform
Naturally, we'd like to generalize: what's the relationship between f(t) and F(ω)=ω2sin(ω), and vice versa?
Theorem. (Fourier Transform) Given a CT signal f(t), define its Fourier transform to be F(ω):=∫−∞∞f(t)e−jωtdt. It then follows that f(t)=2π1∫−∞∞F(ω)ejtωdω.
Proof: Follow the same recipe as before. Let fT(t) denote the extension of f(t) to have period T. Then its Fourier series is:
FT[k]=T1∫Tf(t)e−jt(2πk/T)dt
In our prototypical example, we picked F(ω)=ω2sin(ω) so that F(ω):=limT→∞T⋅FT[k], where ω=2πk/T. Let's follow the recipe and see what happens.
Remark. Comparing and contrasting periodic and aperiodic signals…
All the frequencies in a periodic signal are captured in one period (∫T), whereas aperiodic signals need all of R(∫−∞∞).
Periodic signals have a discrete set of frequencies (∑k=−∞∞), whereas aperiodic signals have a continuous set (∫−∞∞dω).
Importantly, there is now a clearer dual relationship between f(t) and F(ω). If the equation F(ω):=∫−∞∞f(t)e−jωtdt represents a Fourier transform, then the equation f(t)=2π1∫−∞∞F(ω)ejtωdω represents an inverse Fourier transform.
Theorem. (Time-Frequency Duality) If the Fourier transform of f(t) is F(ω), then the Fourier transform of F(t) is 2πf(−ω).
Proof: Just look at the formulas. ■
§ Lecture and Recitation: Properties of the CT Fourier Transform
Theorem. (Time Delay) The FT of f(t−τ) is e−jωτF(ω), and the inverse FT of F(ω−ω0) is ejtω0f(t).
Proof: True by inspecting the relationship F(ω):=∫−∞∞f(t)e−jωtdt. ■
Note that this suggests that our prototypical example was special in having a real-valued Fourier transform. Here's why.
Theorem. (Even / Odd) If f(t) is real and even, F(ω) is real and even. If f(t) is real and odd, F(ω) is imaginary and odd.
Proof: Again, true by inspecting the relationship F(ω):=∫−∞∞f(t)e−jωtdt. Another reason to see why is to apply the same argument as used for periodic signals; even signals are sums of pure cosines, and odd signals are sums of pure sines. ■
More generally, we can say the following:
Theorem. (Hermitian Symmetry) If f(t) is real, then F(−ω)=F(ω) for all ω.
Proof: Again, just inspect the relationship F(ω):=∫−∞∞f(t)e−jωtdt. ■
We can also say that stretching a CT signal in one direction will cause its Fourier transform to stretch in the other direction.
Theorem. (Time Scaling) The FT of f(at) is ∣a∣1F(aω), and the inverse FT of F(aω) is ∣a∣1f(a1t).
Proof: Again, just inspect the relationship F(ω):=∫−∞∞f(t)e−jωtdt. The absolute value signs ∣a∣ are necessary because when a is negative, the u-substitution u=at causes the bounds of the integral to flip.
Alternatively, we can argue heuristically. This relationship says that, say, doubling the width of f(t) causes its FT to shrink in width and double in height.
This should make sense; the signal f(t/2) can be thought of as “playing back the signal f(t) at ×0.5 speed”, so the frequencies prominent in f(t/2) should be half of those prominent in f(t). (And to balance it out, all frequency strengths are doubled.) ■
Here's some intuition as to why “balancing it out” is necessary.
Theorem. (FT Areas) The area under the curve of f(t) is F(0), and the area under the curve of F(ω) is 2πf(0).
Proof: Immediate from plugging in t=0 or ω=0 as needed. ■
Finally, here's a property whose frequency half doesn't have an analogue for Fourier series.
Theorem. (Derivative FTs) The FT of dtdf(t) is jωF(ω), and the inverse FT of dωdF(ω) is −jtf(t).
Proof: This one is trickier. To show the FT of dtdf(t) is jωF(ω), differentiate the synthesis equation.
Consider our prototypical function f(t) from earlier. Imagine compressing f(t) inward by a factor of A, producing functions like:
fA(t):={A0 if t∈(−A1,A1) otherwise.
Then the area under the curve of fA(t) is constantly 2, so FA(0) must be constantly 2. According to “Time Scaling”, in fact, the effect is that FA(ω) “stretches outward” as A tends to infinity, while still maintaining a y-intercept of 2.
At the limit of A→∞, we discover the following (up to a factor of 2).
Theorem. (Dirac Delta FT) Let the Dirac delta functionδ(t) be defined so that δ(t)=0 for all t=0 and ∫−∞∞δ(t)dt=1. Then the FT of δ(t) is the constant function F(ω)=1.
The dual also holds: the FT of the constant function f(t)=1 is F(ω)=2πδ(ω).
Proof: This follows from the heuristic argument from earlier. You can also formally verify this using the FT equations. ■
The Dirac delta function also appears when taking the Fourier transform of a periodic signal.
Theorem. (Periodic FT) Suppose f(t) is a periodic signal with Fourier coefficients F[k], so that f(t)=∑k=−∞∞F[k]ejt(2πk/T). Then the Fourier transform of f(t) is F(ω)=∑k=−∞∞2πF[k]δ(ω−T2πk).
Proof: Nicely, the FT is linear, so it suffices to find the FT of F[k]ejt(2πk/T). This can be done nicely via “Time Delay”.
This should also just heuristically make sense: integrating over a discrete set of Dirac deltas is the same as summing over a discrete set of complex exponentials. ■