zoukankan      html  css  js  c++  java
  • [傅里叶变换及其应用学习笔记] 七. 傅里叶正(反)变换复习

    这份是本人的学习笔记,课程为网易公开课上的斯坦福大学公开课:傅里叶变换及其应用。

    傅里叶变换没有统一的定义

    符号

    傅里叶变换的符号在不同的书籍可能有不同的写法:

    如正变换的符号:$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^{-2pi ist}f(t)dt }$

    另外有些书本的写法是

    $mathcal{F} f(s) = displaystyle{int_{-infty}^{infty}e^{-ist}f(t)dt }$

    这是由于采用不同的周期而导致的,但是尽管写法不同,但表示的都是同样的意思。

    高斯(Gaussian)函数的傅里叶变换

    高斯函数的归一化(积分为1)式子如下:

    $f(t) = e^{-pi t^2}$

    高斯函数图像如下:

    Fourier 7

    对高斯函数进行积分过程如下:

    由于高斯函数的变量$t$是在幂的位置上,而且是二次方,因此无法直接用$dt$对其进行积分计算。下面采用极坐标方法

    $egin{align*}
    left(displaystyle{int_{-infty}^{infty}e^{-pi t^2}dt} ight)^2
    &=displaystyle{int_{-infty}^{infty}e^{-pi x^2}dx imes int_{-infty}^{infty}e^{-pi y^2}dy}\
    &=displaystyle{iint_{-infty}^{infty}e^{-pi(x^2+y^2)}dxdy}\
    &=int_0^{2pi}int_{0}^{infty}e^{-pi r^2}rdrd heta\
    &=2piint_{0}^{infty}e^{-pi r^2}rdr\
    &=2piint_{0}^{infty}e^{-pi r^2}d(frac{1}{2}r^2)\
    &=frac{2pi}{pi} imesfrac{1}{2}int_0^{infty}e^{-pi r^2}dpi r^2\
    &=int_0^{infty}e^{-s}ds\
    &=left. -e^{-s} ight|_0^{infty}\
    &=0-(-1)\
    &=1
    end{align*}$

    那么该高斯函数的积分为

    $displaystyle{int_{-infty}^{infty}e^{-pi t^2}dt = sqrt{1} = 1 }$

    下面对高斯函数进行傅里叶变换

    $egin{align*}
    F(s)=mathcal{F} f(s)
    &=int_{-infty}^{infty}e^{-2pi ist}e^{-pi t^2}dt
    end{align*}$

    这也是一个非常难以积分的项,我们需要采用其他巧妙的方法:微分

    $egin{align*}
    F'(s)=mathcal{F} f'(s)
    &=int_{-infty}^{infty}frac{d(e^{-2pi ist})}{ds}e^{-pi t^2}dt\
    &=int_{-infty}^{infty}-2pi ite^{-2pi ist}e^{-pi t^2}dt\
    &=iint_{-infty}^{infty}e^{-2pi ist}(-2pi te^{-pi t^2})dt\
    &=ileft(left. e^{-2pi ist}e^{-pi t^2} ight|_{-infty}^{infty}-int_{-infty}^{infty}e^{-pi t^2}(-2pi ise^{-2pi ist})dt ight)\
    &=-2pi sint_{-infty}^{infty}e^{-2pi ist}e^{-pi t^2}dtqquad eliminate  left. e^{-2pi ist}e^{-pi t^2} ight|_{-infty}^{infty} because |e^{-2pi ist}|=1,lim_{t oinfty}e^{-pi t^2}=0\
    &=-2pi sF(s)
    end{align*}$

    求偏微分方程,得

    $F(s) = F(0)e^{-pi s^2} = displaystyle{int_{-infty}^{infty}e^{-pi t^2}dt imes e^{-pi s^2} } = e^{-pi s^2}$

    也就是说归化为1的高斯函数的傅里叶变换还是归化为1的高斯函数

    反转信号(reverse signal)

    这是一个新的定义,目的是为了方便式子的表达,定义如下

    令$f^{-}(t) = f(-t)$

    $f^{-}(t)$即为$f(t)$的反转

    傅里叶变换的对偶性(Fourier Transform Duality)

    回顾一下傅里叶变换:

    $F(s) = mathcal{F} f(s) = displaystyle{int_{-infty}^{infty}e^{-2pi ist}f(t)dt }$

    当取值为$-s$时,

    $F(-s) = mathcal{F} f(-s) = displaystyle{int_{-infty}^{infty}e^{2pi ist}f(t)dt } = mathcal{F}^{-1}f(s)$

    一般来说,$f(t)$是时域,$F(s)$是频域,$f(t)$通过傅里叶变换得到$F(s)$,$F(s)$通过逆变换得到$f(t)$。不过上面的式子是对$f(t)$进行傅里叶逆变换,在这里,我们并不需要分析这个等式所表示的含义,而是把傅里叶变换当作工具使用。

    对偶定理1

    把反转信号引入傅里叶变换的对偶性中,得$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$

    函数的傅里叶变换的反转等于对该函数进行傅里叶逆变换。

    对偶定理2

    如果对$f^{-}(t)$进行傅里叶变换会得到什么结果呢?

    $egin{align*}
    mathcal{F}(f^{-}(s))
    &= int_{-infty}^{infty}e^{-2pi ist}f(-t)dt\
    &= int_{+infty}^{-infty}e^{-2pi is(-u)}f(u)d(-u) qquad let u=-t\
    &= int_{-infty}^{infty}e^{2pi isu}f(u)du\
    &= mathcal{F}^{-1}f(s)
    end{align*}$

    即,

    $mathcal{F}(f^{-}) = mathcal{F}^{-}f$

    函数的反转的傅里叶变换等于对该函数进行傅里叶逆变换。

    对偶定理3

    把对偶定理1与对偶定理2结合起来,得

    $(mathcal{F} f)^{-} = mathcal{F}(f^{-})$

    函数的傅里叶变换的反转等于对该函数反转的傅里叶变换

    对偶定理4

    对函数进行两次傅里叶变换

    $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$

  • 相关阅读:
    Android 入门到精通 (Index)
    负载平衡与冗余备份方案概述
    Android 程序组件交互分析
    复制时保留文件的目录结构
    notepad++中设置tab缩进的宽度
    scws
    php 将字符(包括汉字) 转换成16进制 (apache access log 中文显示16进制码)
    批量修改完整版本
    根据端口号查进程
    php性能优化
  • 原文地址:https://www.cnblogs.com/TaigaCon/p/5001598.html
Copyright © 2011-2022 走看看