这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用。
傅里叶变换的符号在不同的书籍可能有不同的写法:
如正变换的符号:$\mathcal{F} f(s)$,$\hat{f}(s)$,$F(s)$
如反变换的符号:$\mathcal{F}^{-1}f(t)$,$\check{f}(t)$,$f(t)$
傅里叶变换的公式也没有统一的写法:
本课程采用的是如下公式
$\mathcal{F} f(s) = \displaystyle{\int_{-\infty}^{\infty}e^{-\pi ist}f(t)dt }$
另外有些书本的写法是
$\mathcal{F} f(s) = \displaystyle{\int_{-\infty}^{\infty}e^{-ist}f(t)dt }$
这是由于采用不同的周期而导致的,但是尽管写法不同,但表示的都是同样的意思。
高斯函数的归一化(积分为1)式子如下:
$f(t) = e^{-\pi t^2}$
高斯函数图像如下:
对高斯函数进行积分过程如下:
由于高斯函数的变量$t$是在幂的位置上,而且是二次方,因此无法直接用$dt$对其进行积分计算。下面采用极坐标方法
$\begin{align*}
\left(\displaystyle{\int_{-\infty}^{\infty}e^{-\pi t^2}dt}\right)^2
&=\displaystyle{\int_{-\infty}^{\infty}e^{-\pi x^2}dx\times \int_{-\infty}^{\infty}e^{-\pi y^2}dy}\\
&=\displaystyle{\iint_{-\infty}^{\infty}e^{-\pi(x^2+y^2)}dxdy}\\
&=\int_0^{2\pi}\int_{0}^{\infty}e^{-\pi r^2}rdrd\theta\\
&=2\pi\int_{0}^{\infty}e^{-\pi r^2}rdr\\
&=2\pi\int_{0}^{\infty}e^{-\pi r^2}d(\frac{1}{2}r^2)\\
&=\frac{2\pi}{\pi}\times\frac{1}{2}\int_0^{\infty}e^{-\pi r^2}d\pi r^2\\
&=\int_0^{\infty}e^{-s}ds\\
&=\left. -e^{-s}\right|_0^{\infty}\\
&=0-(-1)\\
&=1
\end{align*}$
那么该高斯函数的积分为
$\displaystyle{\int_{-\infty}^{\infty}e^{-\pi t^2}dt = \sqrt{1} = 1 }$
下面对高斯函数进行傅里叶变换
$\begin{align*}
F(s)=\mathcal{F} f(s)
&=\int_{-\infty}^{\infty}e^{-2\pi ist}e^{-\pi t^2}dt
\end{align*}$
这也是一个非常难以积分的项,我们需要采用其他巧妙的方法:微分
$\begin{align*}
F'(s)=\mathcal{F} f'(s)
&=\int_{-\infty}^{\infty}\frac{d(e^{-2\pi ist})}{ds}e^{-\pi t^2}dt\\
&=\int_{-\infty}^{\infty}-2\pi ite^{-2\pi ist}e^{-\pi t^2}dt\\
&=i\int_{-\infty}^{\infty}e^{-2\pi ist}(-2\pi te^{-\pi t^2})dt\\
&=i\left(\left. e^{-2\pi ist}e^{-\pi t^2}\right|_{-\infty}^{\infty}-\int_{-\infty}^{\infty}e^{-\pi t^2}(-2\pi ise^{-2\pi ist})dt\right)\\
&=-2\pi s\int_{-\infty}^{\infty}e^{-2\pi ist}e^{-\pi t^2}dt\qquad eliminate\ \left. e^{-2\pi ist}e^{-\pi t^2}\right|_{-\infty}^{\infty}\ because\ |e^{-2\pi ist}|=1,\lim_{t\to\infty}e^{-\pi t^2}=0\\
&=-2\pi sF(s)
\end{align*}$
求偏微分方程,得
$F(s) = F(0)e^{-\pi s^2} = \displaystyle{\int_{-\infty}^{\infty}e^{-\pi t^2}dt\times e^{-\pi s^2} } = e^{-\pi s^2}$
也就是说归化为1的高斯函数的傅里叶变换还是归化为1的高斯函数
这是一个新的定义,目的是为了方便式子的表达,定义如下
令$f^{-}(t) = f(-t)$
$f^{-}(t)$即为$f(t)$的反转
回顾一下傅里叶变换:
$F(s) = \mathcal{F} f(s) = \displaystyle{\int_{-\infty}^{\infty}e^{-2\pi ist}f(t)dt }$
当取值为$-s$时,
$F(-s) = \mathcal{F} f(-s) = \displaystyle{\int_{-\infty}^{\infty}e^{2\pi ist}f(t)dt } = \mathcal{F}_{-1}f(s)$
一般来说,$f(t)$是时域,$F(s)$是频域,$f(t)$通过傅里叶变换得到$F(s)$,$F(s)$通过逆变换得到$f(t)$。不过上面的式子是对$f(t)$进行傅里叶逆变换,在这里,我们并不需要分析这个等式所表示的含义,而是把傅里叶变换当作工具使用。
把反转信号引入傅里叶变换的对偶性中,得$\mathcal{F} f(-s) = (\mathcal{F} f)^{-}(s)$,而且上面对偶性讨论已得出结论:$\mathcal{F} f(-s) = \mathcal{F}^{-1}f(s)$,即有
$(\mathcal{F} f)^{-}(s) = \mathcal{F}^{-1}f(s)$
$(\mathcal{F} f)^{-} = \mathcal{F}^{-}f$
函数的傅里叶变换的反转等于对该函数进行傅里叶逆变换。
如果对$f^{-}(t)$进行傅里叶变换会得到什么结果呢?
$\begin{align*}
\mathcal{F}(f^{-}(s))
&= \int_{-\infty}^{\infty}e^{-2\pi ist}f(-t)dt\\
&= \int_{+\infty}^{-\infty}e^{-2\pi is(-u)}f(u)d(-t) \qquad let \ u=-t\\
&= \int_{-\infty}^{\infty}e^{2\pi isu}f(u)du\\
&= \mathcal{F}^{-1}f(s)
\end{align*}$
即,
$\mathcal{F}(f^{-}) = \mathcal{F}^{-}f$
函数的反转的傅里叶变换等于对该函数进行傅里叶逆变换。
把对偶定理1与对偶定理2结合起来,得
$(\mathcal{F} f)^{-} = \mathcal{F}(f^{-})$
函数的傅里叶变换的反转等于对该函数反转的傅里叶变换
对函数进行两次傅里叶变换
$\mathcal{F}\mathcal{F} f = \mathcal{F}(\mathcal{F} f) = \mathcal{F} (\mathcal{F}^{-}(f^{-})) = f^{-}$
函数连续进行两次傅里叶变换等于该函数的反转。
对偶定理的目的是为了方便计算,如:
求$sinc$函数的傅里叶变换。
$sinc = \frac{sin \pi s}{\pi s}$
由上一节课我们知道$\pi$函数经过傅里叶变换后得到$sinc$函数,那么我们就运用傅里叶变换的对偶定理能进行如下计算
$\mathcal{F} sinc = \mathcal{F}\mathcal{F} \pi = \pi^{-} = \pi$