User Tools

Site Tools


math121a-f23:october_4_wednesday

October 4 (Wednesday)

What's coming in the second part of this course?

  • Fourier transform. Given f(x)f(x), we want to express f(x)=eikxg(k)dkf(x) = \int e^{ikx} g(k) dk, since eikxe^{ikx} is easy to deal with.
  • Laplace transform. Given f(t)f(t), on t[0,)t \in [0, \infty), we want to write f(t)=c+iReptF(p)dpf(t) = \int_{c + i \R} e^{pt} F(p) dp.
  • Gamma function, Stirling Formula

What is Fourier Transformation

When we solve equation of the form df(x)dx=λf(x), \frac{df(x)}{dx} = \lambda f(x), the general soluition is f(x)=ceλx. f(x) = c e^{\lambda x}. We say eλxe^{\lambda x} is an eigenfunction for the operator d/dxd/dx with eigenvalue λ\lambda, here λ\lambda can be any complex number.

  • For λ\lambda purely imaginary, λ=ik\lambda = i k, the function eikxe^{i k x}, which has constant size eikx=1|e^{ikx}| = 1.
  • For λ\lambda purely real, λR\lambda \in \R, the function eλxe^{\lambda x}, which is real, and has exponential growth (if λ>0\lambda > 0) or exponential decay.
  • For λ\lambda general complex number, say λ=a+ib\lambda = a + i b for a,bRa,b\in\R, then eλx=eaxeibxe^{\lambda x} = e^{ax} e^{ibx} has both oscillation factor eibxe^{ibx} and exponential growth/decay eaxe^{ax}.

If you give me a function f(x)f(x), we can try to write it as a linear combination of these easy to understand pieces eλxe^{\lambda x}. Say f(x)=Cc(λ)eλxdλ f(x) = \int_{C} c(\lambda) e^{\lambda x} d\lambda The benefit for doing this? If a differential operator d/dxd/dx walk to it, we would have (d/dx)f(x)=λc(λ)(d/dx)eλxdλ=λc(λ)λeλxdλ. (d/dx) f(x) = \int_{\lambda} c(\lambda) (d/dx) e^{\lambda x} d\lambda = \int_{\lambda} c(\lambda) \lambda e^{\lambda x} d\lambda. Of course, you would immediately complain:

  • what is this integration contour CC? Is it along the real axis of λ\lambda? Is it along the imaginary axis of λ\lambda?
  • Why we can switch the order of integration and differentiation?

Let's keep these questions in our head, we will return to those.

So, what is Fourier transformation? If you go on wiki, or open some textbook, you see the following definition: given a function f(x)f(x) (with some condition on it), we can define F(p):=xRf(x)eipxdx. F(p) := \int_{x \in \R} f(x) e^{-ipx} dx. This is called the Fourier transformation.

Let's throw in some input and see what we get.

Example 1

f(x)=cosh(x)=(ex+ex)/2. f(x) = \cosh(x) = (e^{x} + e^{-x} )/2. This function grows to infinity when x±x \to \pm \infty, indeed, when x+x \to +\infty, exe^x dominate, and xx \to -\infty we have exe^{-x} dominate. Either way, you blow up. Good luck when integrating Rf(x)dx\int_\R f(x) dx.

But, let's try to do the Fourier transformation anyway. We get F(p)=xf(x)eipxdx=(1/2)Re(1ip)xdx+(1/2)Re(1ip)xdx. F(p) = \int_x f(x) e^{-ipx} dx = (1/2) \int_\R e^{(1-ip)x} dx + (1/2) \int_\R e^{(-1-ip)x} dx. OK, we cannot continue.

Lesson 1

f(x)f(x) need to have sufficient 'decay' near infinity for the Fourier transformation to make sense. More precisely, a sufficient condition is that f(x)f(x) is 'absolutely integrable', which means Rf(x)dx<. \int_\R |f(x)|dx < \infty. Given this condition, we know F(p)=Reipxf(x)dxReipxf(x)dx=Rf(x)dx<. |F(p)| = | \int_\R e^{-ipx} f(x) dx | \leq \int_\R |e^{-ipx} f(x)| dx = \int_\R |f(x)| dx < \infty. so we are safe.

Example 2

The Gaussian. f(x)=ex2.f(x) = e^{-x^2}. It decays pretty nicely in both directions.

Since we learned our lesson, we need to do a check first Rex2dx=π. \int_\R e^{-x^2} dx = \sqrt{\pi}. Great, it is finite.

Now, let's Fourier transform. F(p)=eipxx2dx=Re(xip/2)2p2/4dx=ep2/4uip/2+Reu2du F(p) = \int e^{-ipx - x^2} dx = \int_\R e^{- (x-ip/2)^2 - p^2/4} dx = e^{-p^2/4} \int_{u \in -ip/2+\R} e^{-u^2} du where we used the new variable u=ip/2+xu = -ip/2+x.

We can move the contour for uu as long as the integral is convergent. Recall from last Wednesday's lecture, eu2e^{-u^2} decays fast to 00 when u|u| \to \infty in the sectors arg(u)(π/4,π/4)\arg(u) \in (-\pi/4, \pi/4) and arg(u)π+(π/4,π/4)\arg(u) \in \pi + (-\pi/4, \pi/4). So, we are going to move the integration contour of uu from the shifted real line ip/2+R-ip/2+\R back to R\R. We get uip/2+Reu2du=uReu2du=π. \int_{u \in -ip/2+\R} e^{-u^2} du = \int_{u \in \R} e^{-u^2} du = \sqrt{\pi}. Thus, we get F(p)=πep2/4. F(p) = \sqrt{\pi} e^{-p^2/4}.

Success!

Example 3

Let f(x)=1/(1+x2)f(x) = 1/(1+x^2). Let's check that f(x)f(x) is absolutely integrable first Rf(x)dx=CRf(z)dz=2πiResz=if(z)=2πi12i=π. \int_\R f(x) dx = \int_{C_R} f(z) dz = 2\pi i Res_{z=i} f(z) = 2\pi i \frac{1}{2i } = \pi. where C+,RC_{+,R} is the contour consist of two part

  • [R,R][-R, R]
  • a semi-circle in the upper half plane of radius RR.

Equivalently, we can choose a different contour C,RC_{-,R}, which also consist of two part

  • [R,R][-R, R]
  • a semi-circle in the lower half plane of radius RR.

Note that when we go around C,RC_{-,R} it goes around the singularity of f(z)f(z) at z=iz=-i in the 'negative' (clockwise) direction, hence we get C,Rf(z)dz=(2πi)Resz=if(z)=2πi12i=π. \int_{C_{-,R}} f(z) dz = (-2\pi i) Res_{z=-i} f(z) = -2\pi i \frac{1}{-2i} = \pi. Whew! the two minus sign cancels, we get the same answer.

Now, let's compute the little variation F(p)=Reipxf(x)dx. F(p) = \int_\R e^{-ipx} f(x) dx. We first take the 'truncation' approximation, which recovers the original limit under the limit RR \to \infty. F(p)=limRRReipxf(x)dx. F(p) = \lim_{R \to \infty} \int_{-R}^R e^{-ipx} f(x) dx. Then, we try to 'close the contour', by connecting the point RR to R-R along some way. Due to this exponential factor eipze^{-ipz}, it matters whether we close the contour in the upper half-plane, or the lower one. Indeed, we have eipz=eRe(ip(x+iy))=eRe(ipx+py)=epy=epIm(z). |e^{-ipz}| = e^{Re(-ip(x+iy))} = e^{Re(-ipx + py)} = e^{py} = e^{p Im(z)}. So,

  • if p0p \geq 0, to avoid eipz|e^{-ipz}| blow up to infinity, we need to keep Im(z)Im(z) bounded from above on the contour.
  • if p0p \leq 0, to avoid eipz|e^{-ipz}| blow up to infinity, we need to keep Im(z)Im(z) bounded from below on the contour.

Thus, for p0p \geq 0, we can complete the by the lower-half-semicircle C,RC_{-,R}. We get (p0)F(p)=C,Reipz1+z2dz=(2πi)Resz=ieipz1+z2=πep. (p\geq 0) F(p) = \int_{C_{-,R}} \frac{e^{-ipz}}{1+z^2} dz = (-2\pi i) Res_{z=-i} \frac{e^{-ipz}}{1+z^2} = \pi e^{-p}. Similarly, for p0p \geq 0, we get (p0)F(p)=C+,Reipz1+z2dz=(2πi)Resz=ieipz1+z2=πep. (p\leq 0) F(p) = \int_{C_{+,R}} \frac{e^{-ipz}}{1+z^2} dz = (2\pi i) Res_{z=i} \frac{e^{-ipz}}{1+z^2} = \pi e^{p}. We can write the result in a cool way as F(p)=πep. F(p) = \pi e^{-|p|}. So when p|p| \to \infty, we see F(p)0F(p) \to 0 very quickly.

Exercise

warm up:Does the function f(x)=1f(x)=1 admits Fourier transformation? Why? How about f(x)=1/(1+x)f(x) = 1 / (1+|x|)?

Find the Fourier transformation of the following functions. f(x)={10<x<10else f(x) = \begin{cases} 1 & 0<x<1 \cr 0 & \text{else} \end{cases}

f(x)={1+x1<x<01x0<x<10else f(x) = \begin{cases} 1+x & -1<x<0 \cr 1-x & 0<x<1 \cr 0 & \text{else} \end{cases}

math121a-f23/october_4_wednesday.txt · Last modified: 2023/10/05 22:17 by pzhou