概述

希尔伯特空间是一个完备的内积空间，其标准正交函数系，直观来看就是向量空间中基的延伸。其为基于任意正交系上的多项式表示的傅立叶级数和傅立叶变换提供了一种有效的表述方式，而这也是泛函分析的核心概念之一。下文中我们将通过希尔伯特空间的标准正交函数系推导周期函数和有限区间上函数的傅立叶级数表示，并进一步推出傅里叶积分来表示无穷区间的非周期函数，最后引入复数形式的傅立叶积分，引出傅立叶变换。在这一系列推导中，鉴于篇幅，主动略去了一些比较关键的部分，比如可积性及级数收敛性的讨论，有兴趣的读者可以在了解大致原理后，进行细致的理论推导以作补充。为了便于理解希尔伯特空间的概念，引用知乎上面的一段回答：

知乎问答

希尔伯特空间

若无限维酉空间中每个基本序列收敛于V中的元素，则称是完备的。一个完备的无限维酉空间称为希尔伯特空间又称为空间。在维空间中矢量被定义为,在无限维空间中矢量被定义为从变换到的函数。在希尔伯特空间中矢量的加法和乘法定义为函数的加法与函数和数的乘法。

内积

对于,则二者的内积定义为：
度量
- 对于,其长度定义为：
若则二者之间的距离定义为：

直观来看就是两个函数的均方差，就是以均方差来作为空间的距离的度量。
若,则二者的夹角定义为：
正交函数系
- 若非零矢量的内积，则由空间的夹角定义公式可知,此时称矢量正交。
- 若且两两正交，设
  
  则有
  
  由空间长度和正交的定义可推出：
  
  上述积分都是指勒贝格积分有意义。
- 若函数系中任意两个函数相互正交，即
  
  则称这个函数系为正交函数系，若还满足
  
  则称此函数系为标准正交系。
依标准正交函数系的分解

若在空间中给定一个完备的标准正交函数系（即不可能再加一个不恒为零的函数与系中的一切函数正交），则对于任意函数都可根据这个标准正交函数系展开成级数（平均收敛）：

为在这个标准正交函数系上的投影：

很容易证明

代表空间中矢量的长度平方等于该矢量在完备的标准正交系中的矢量上的投影平方和。

傅立叶级数

希尔伯特空间是有限维欧几里得空间的推广，与欧几里得空间相同，希尔伯特空间也是内积空间，也有距离和角的概念，并且不同于欧几里得空间，空间具有完备性：希尔伯特空间内的所有的柯西列会收敛到一点。因此微积分中的大部分概念可无障碍推广至希尔伯特空间中。希尔伯特空间提供了一个很强大理论工具：对于空间中的任意函数都可以由空间中完备的标准正交函数系展开成级数。也就是说，可以通过一个标准正交函数系去逼近一个任意函数（这些函数都是基于希尔伯特空间的）。

三角函数系

基于上述的讨论，接下来讨论如何通过一个标准正交函数系来展开任意函数。基于标准正交函数的定义，直观来看，三角函数系似乎完美的切合，对于三角函数系

由于

又由奇函数的性质直接得到（15）（16）定积分等式成立：

又：

同理得到：

于是得到三角函数系（13）在是一个正交函数系，但是，细心的读者可能发现，三角函数系（13）并不是一个标准正交函数系，为此基于（13）构造三角函数系

显然，三角函数系（20）是一个正交函数系，下面进一步证明该函数系是一个标准正交函数系，则只需证明对于该函数系任意函数有（9）式成立即可：

$\int_{-\pi}^{\pi}(\frac{cos(kx)}{\sqrt{\pi}})^2\mathrm{d}x=\int_{-\pi}^{\pi}\frac{cos^2(kx)}{\pi}\mathrm{d}x\overset{三角降幂公式}{=}\int_{-\pi}^{\pi}\frac{1+cos(2kx)}{2\pi}\mathrm{d}x=(\frac{1}{2\pi}x+\frac{1}{4k\pi}sin(k\pi))｜_{-\pi}^{\pi}=1\\\tag{21}$

$\int_{-\pi}^{\pi}(\frac{sin(kx)}{\sqrt{\pi}})^2\mathrm{d}x=\int_{-\pi}^{\pi}\frac{sin^2(kx)}{\pi}\mathrm{d}x\overset{三角降幂公式}{=}\int_{-\pi}^{\pi}\frac{1-cos(2kx)}{2\pi}\mathrm{d}x=(\frac{1}{2\pi}x-\frac{1}{4k\pi}sin(k\pi))｜_{-\pi}^{\pi}=1\tag{22}$

由（20）（21）（22）三式可得，三角函数系（19）为希尔伯特空间下的标准正交函数系。则任意定义在的函数有：

整理得：
$f(x)=c_0\frac{1}{\sqrt{2\pi}}+\sum_{i=1}^{\infty}a_i\frac{cos(ix)}{\sqrt{\pi}}+\sum_{i=0}^{\infty}b_i\frac{sin(ix)}{\sqrt{\pi}}=c_0\frac{1}{\sqrt{2\pi}}+\sum_{i=1}^{\infty}(a_i\frac{cos(ix)}{\sqrt{\pi}}+b_i\frac{sin(ix)}{\sqrt{\pi}})\tag{24}$
此时，只需确定系数即可得到在在的级数展开形式：

:

对等式（24）两边同时求积分得到：

结合（14）（15）（16）得到：

求得:

:

结合（21）式，对等式（24）两边同乘,得到：

同样对等式（30）两边同时求积分得到：
$\int_{-\pi}^{\pi}f(x)cos(jx)\mathrm{d}x=\color{#F00}{\int_{-\pi}^{\pi}c_0\frac{1}{\sqrt{2\pi}}cos(jx)\mathrm{d}x}+\color{#00F}{\sum_{i=1}^{\infty}\int_{-\pi}^{\pi}a_i\frac{cos(ix)cos(jx)}{\sqrt{\pi}}\mathrm{d}x}+\color{#F00}{\sum_{i=1}^{\infty}\int_{-\pi}^{\pi}b_i\frac{sin(ix)cos(jx)}{\sqrt{\pi}}\mathrm{d}x}\tag{29}$
由（14）（16）式可知，等式（29）中两个红色定积分都为，对于蓝色定积分，由（17）（21）式可知，当且仅当的项积分为，其余项的积分都为，故：
$\int_{-\pi}^{\pi}f(x)cos(jx)\mathrm{d}x=\color{#00F}{\sum_{i=1}^{\infty}\int_{-\pi}^{\pi}a_i\frac{cos(ix)cos(jx)}{\sqrt{\pi}}\mathrm{d}x}=\color{#F00}{\sum_{i=1,i\neq j }^{\infty}\int_{-\pi}^{\pi}a_i\frac{cos(ix)cos(jx)}{\sqrt{\pi}}\mathrm{d}x}+\color{#0F0}{\int_{-\pi}^{\pi}a_j\frac{cos(jx)cos(jx)}{\sqrt{\pi}}\mathrm{d}x}\tag{30}$
将蓝色部分的积分拆分为红色和绿色两部分积分之和（绿色积分为时），显然由（17）式红色部分积分仍旧为，由（21）式绿色部分积分为，进一步得到:

:

同上可得到：

周期与非周期下的傅里叶级数

在上文中，已经求出了系数的表达式，将这些系数代入（24）式，整理得到：
$\begin{equation} \begin{split} f(x)&=c_0\frac{1}{\sqrt{2\pi}}+\sum_{i=1}^{\infty}(a_i\frac{cos(ix)}{\sqrt{\pi}}+b_i\frac{sin(ix)}{\sqrt{\pi}})\\ &=（\frac{1}{\sqrt{2\pi}}\int_{-\pi}^{\pi}f(x)\mathrm{d}x）\frac{1}{\sqrt{2\pi}}+\sum_{i=1}^{\infty}[(\frac{1}{\sqrt{\pi}}\int_{-\pi}^{\pi}f(x)cos(ix)\mathrm{d}x)\frac{cos(ix)}{\sqrt{\pi}}+(\frac{1}{\sqrt{\pi}}\int_{-\pi}^{\pi}f(x)sin(ix)\mathrm{d}x)\frac{sin(ix)}{\sqrt{\pi}}]\\ &=\frac{1}{2\pi}\int_{-\pi}^{\pi}f(x)\mathrm{d}x+\sum_{i=0}^{\infty}[\color{#00F}{\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)cos(ix)\mathrm{d}x} * \color{#F0F}{cos(ix)}+\color{#00F}{\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)sin(ix)\mathrm{d}x} * \color{#F0F}{sin(ix)}] \end{split} \end{equation}\tag{33}$
此时，令
$\begin{equation} \begin{split} a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)\mathrm{d}x\\ a_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)cos(nx)\mathrm{d}x\\ b_n&=\frac{1}{\pi}\int_{-\pi}^{\pi}f(x)sin(nx)\mathrm{d}x \end{split} \end{equation}\tag{34}$
得到：

式（35）显然就是傅立叶级数。级数（35）的收敛性证明，涉及较多泛函分析的内容，在此不做展开，有时间在另开一篇文章说明，有兴趣的读者可以尝试一下证明。细心的读者可能注意到的定义域是，那么如果对于任意周期，级数还成立吗？将定义域拓展到实数域后，级数（35）还成立吗？对此，下面进一步讨论。

周期函数的傅立叶级数

当是一个周期为的周期函数时，那么在区间中，作变换：

则，为周期的周期函数，构造标准正交三角函数系：

进一步得到此时的傅立叶系数：
$\begin{equation} \begin{split} a_0&=\frac{1}{T}\int_{-T}^{T}f(x)\mathrm{d}x\\ a_n&=\frac{1}{T}\int_{-T}^{T}f(x)cos(\frac{n\pi}{T}x)\mathrm{d}x\\ b_n&=\frac{1}{T}\int_{-T}^{T}f(x)sin(\frac{n\pi}{T}x)\mathrm{d}x \end{split} \end{equation}\tag{38}$
则此时的傅立叶级数为：

因此，对于可积的任意周期函数，都能展开为对应的傅立叶级数，并且由于周期函数的特性，当定义域拓展到实数域时也是成立的。

非周期函数的傅立叶积分

在更多情况下，需要处理非周期函数，对此比较直观的处理技巧是，将非周期函数视为周期为的周期函数，即。设是定义在上，并且在定义域绝对可积，是在有限区间上的截取，因为可视为该有限区间上的周期函数，令则得到的傅立叶级数为：

其中：
$\begin{equation} \begin{split} a_0&=\frac{1}{T}\int_{-T}^{T}f(t)\mathrm{d}t\\ a_n&=\frac{1}{T}\int_{-T}^{T}f(t)cos(n\omega t)\mathrm{d}t\\ b_n&=\frac{1}{T}\int_{-T}^{T}f(t)sin(n \omega t)\mathrm{d}t \end{split} \end{equation}\tag{41}$
得到：
$\begin{equation} \begin{split} \xi_T(x)&=\frac{1}{2T}\int_{-T}^{T}f(t)\mathrm{d}t+\sum_{n=1}^{\infty}[\frac{1}{T}\int_{-T}^{T}f(t)cos(n\omega t)\mathrm{d}tcos(n\omega x)+\frac{1}{T}\int_{-T}^{T}f(t)sin(n \omega t)\mathrm{d}tsin(n\omega x)] \\ &=\frac{1}{2T}\int_{-T}^{T}f(t)\mathrm{d}t+\sum_{n=1}^{\infty}\frac{1}{T}\int_{-T}^{T}f(t)[cos(n\omega t)cos(n\omega x)+sin(n \omega t)sin(n\omega x)]\mathrm{d}t \\ &=\frac{1}{2T}\int_{-T}^{T}f(t)\mathrm{d}t+\sum_{n=1}^{\infty}\frac{1}{T}\int_{-T}^{T}f(t)cos[n\omega(x-t)]\mathrm{d}t \end{split} \end{equation}\tag{42}$
于是有：

关（43）式的两个极限，接下来分别进行讨论：

对于第一个极限：

由于在上绝对可积，则：

所以，有：

对于第二个极限，先直接给出结论：

令，有：

于是，得到非周期函数的傅立叶积分表示：

或写为：

其中：

综合来看，周期函数和有限区间上的函数可以用傅立叶级数来表示；而无穷区间的非周期函数，用傅里叶积分表示，对应了频率的连续分布。

傅里叶变换

傅里叶级数的本质是函数在某个函数空间中各个基底的投影和，在上文中我们通过引入希尔伯特空间，构造标准正交三角函数系进而推导出傅立叶级数与傅立叶积分。然而，这一切都是基于实数域推导，那么这种思路在复数域是否也成立呢，答案是显而易见的。下面，我们通过欧拉公式（公式证明见文末）来得到傅立叶积分的复数形式：

观察（47）式，由于：

因为
$\begin{equation} \begin{split} \int_{-\alpha}^{\alpha}f(t)e^{i\lambda(x-t)}\mathrm{d}t&=\int_{-\alpha}^{\alpha}f(t)\{cos[\lambda(x-t)]+isin[\lambda(x-t)]\}\mathrm{d}t\\ &=\int_{-\alpha}^{\alpha}f(t)cos[\lambda(x-t)]\mathrm{d}t+i*0（奇函数在对称区间的积分为零） \end{split} \end{equation}\tag{52}$
则非周期函数的傅立叶积分（47）可改写为：
$\begin{equation} \begin{split} f(x)&=\frac{1}{\pi}\int_{0}^{+\infty}[\int_{-\infty}^{+\infty}f(t)cos[\lambda(x-t)]\mathrm{d}t]\mathrm{d}\lambda\\ &=\frac{1}{\pi}\int_{0}^{+\infty}[\int_{-\infty}^{+\infty}f(t)e^{i\lambda(x-t)}\mathrm{d}t]\mathrm{d}\lambda\\ &=\frac{1}{\pi}\int_{0}^{+\infty}e^{i\lambda x}\int_{-\infty}^{+\infty}f(t)e^{-i\lambda t}\mathrm{d}t\mathrm{d}\lambda\\ &=\frac{1}{2\pi}\int_{-\infty}^{+\infty}e^{i\lambda x}\color{#F00}{\int_{-\infty}^{+\infty}f(t)e^{-i\lambda t}\mathrm{d}t}\mathrm{d}\lambda\\ \end{split} \end{equation}\tag{53}$
上式最后一个等式即为傅立叶积分的复数形式，而红色积分部分就是大名鼎鼎的傅立叶变换也叫像函数，是一个复数表示振幅和相位：

而也称为傅立叶逆变换也叫本函数：

傅里叶变换一词既指变换操作本身（将函数进行傅里叶变换），又指该操作所生成的复数函数（是的傅里叶变换）,需要注意的是，一般情况下傅立叶变换是可逆的。

傅立叶变换的性质

傅立叶级数使用不同频率的三角函数和来表示周期函数和有限区上的函数，而傅立叶积分则是对频率作无穷积分来表示无穷区间上的函数。本质上其实是从不同的角度刻画相同的函数，所以你经常可以听到这样的说法，傅立叶变换是一种线性积分变换，常用于信号时域到频域之间的变换，这里说的时域是从时间的角度描述函数或信号，而频域则是从频率的角度描述函数或信号。本质上傅立叶变换就像化学分析，像分析物质的基本成分一样，确定函数或信号的基本组成。

接下来讨论一下傅立叶变换的一些基本性质，为了方便描述，约定为傅立叶变换的作用算子，即为的傅立叶变换，表示的傅立叶逆变换，并且函数都存在傅立叶变换：

线性性质

两函数之和的傅里叶变换等于各自的傅立叶变换之和：
频移性质

时移特性
帕塞瓦尔定理

若平方可积，则有：
卷积的傅里叶变换

若且在定义域内绝对可积，定义卷积函数：

则有：
$\begin{equation} \begin{split} &\mathscr{F}[f*g]=\mathscr{F}[f]\cdot\mathscr{F}[g] \\ &\mathscr{F}^{-1}[F(\lambda)*G(\lambda)] = 2\pi\mathscr{F}^{-1}[F(\lambda)]\cdot\mathscr{F}^{-1}[G(\lambda)] \end{split} \end{equation}$
傅里叶变换在时域和频域之间搭起来一座桥梁，一些在时域很难解决甚至无法解决的问题，在频域下却可以轻松得到解决，在信号、图像处理还有偏微分方程等领域都有着广泛的应用。

离散傅里叶变换

离散傅里叶变换（Discrete Fourier Transform，缩写为DFT），是傅里叶变换在时域和频域上都呈离散的形式，将信号的时域采样变换为其离散时间傅里叶变换的频域采样。

对于序列

其离散傅里叶变换（DFT）如下：

离散傅里叶变换的逆变换（IDFT）如下：

离散傅里叶变换的应用：

数据压缩

由于人类感官的分辨能力存在极限，因此很多有损压缩算法利用这一点将语音、音频、图像、视频等信号的高频部分除去。高频信号对应于信号的细节，滤除高频信号可以在人类感官可以接受的范围内获得很高的压缩比。这一去除高频分量的处理就是通过离散傅里叶变换完成的。将时域或空域的信号转换到频域，仅储存或传输较低频率上的系数，在解压缩端采用逆变换即可重建信号。

长整数与多项式乘法

目前长整数或多项式乘法最快速的算法是基于离散傅里叶变换的。由于整数（或多项式）乘法是逐位（或逐项）乘累加的形式，因此整数（或多项式）乘积的数字（或系数）可以用乘数数字（或乘式系数）的卷积表示。利用卷积定理，只要将数字（或系数）序列通过离散傅里叶变换变到频域，就可以将逐个乘累加的卷积变为对位的乘法，从而减少计算量，再以一次逆变换便可以得到乘法结果。需要注意整数乘法还有进位的问题。

求解偏微分方程

离散傅里叶变换及其多维形式在偏微分方程的求解中也有应用。此时DFT被看作傅里叶级数的近似。傅里叶级数将函数在复指数上展开，这正是微分算子的特征方程：

因此，通过傅里叶级数的形式，线性常微分方程被转换为代数方程，而后者是很容易求解的。此时得到的结果是偏微分方程解的级数表示，只要通过DFT逆变换即可得到其一般表示，这种方法被称作谱方法或级数解法。

快速傅里叶变换

离散傅里叶变换十分强大，但是计算复杂度较高，对于一个大小为的序列，其离散傅里叶级数的复杂度为，对于一些需要实时计算的场景，不太能满足需求。因此，快速傅里叶变换（FFT）应运而生，快速傅里叶变换是快速计算序列的离散傅里叶变换]叶变换)（DFT）或其逆变换的方法。傅里叶分析将信号从原始域（通常是时间或空间）转换到频域的表示或者逆过来转换。FFT会通过把DFT矩阵分解为稀疏因子之积来快速计算此类变换，因此，它能够将计算DFT的复杂度从降低到。

FFT的本质就是通过不断的把长序列的DFT分解为几个短序列的DFT，并利用单位根的周期性和对称性来减少计算量。FFT算法有很多种，不过大致可以分为两类：

按抽取方法可分为：

时域抽取法（DIT）

频域抽取大（DIF）

按基数可分为：

基2-FFT算法

基4-FFT算法

混合基FFT算法

分裂基FFT算法

Cooley-Tukey算法

Cooley-Tukey算法是最常见的FFT算法。这一方法以分治法为策略递归地将长度为的离散傅里叶变换分解为长度为的个较短序列的离散傅里叶变换，以及与个转因子的复数乘法。

Cooley-Tukey算法最有名的应用，是将序列长为的DFT分割为两个长为的子序列的DFT，因此这一应用只适用于序列长度为2的幂的DFT计算，即基2-FFT。实际上，如同高斯和Cooley与Tukey都指出的那样，Cooley-Tukey算法也可以用于序列长度N 为任意因数分解形式的DFT，即混合基FFT，而且还可以应用于其他诸如分裂基FFT等变种**。尽管Cooley-Tukey算法的基本思路是采用递归的方法进行计算，大多数传统的算法实现都将显式的递归算法改写为非递归的形式。另外，因为Cooley-Tukey算法是将DFT分解为较小长度的多个DFT，因此它可以同任一种其他的DFT算法联合使用。

基2时间抽取法

基2时间抽取算法是Cooley-Tukey算法的一种分支，当序列的点数为（若不满足，可补零），此时的Cooley-Tukey算法称之为基2时间抽取法。由于序列点数为2的整数幂，则可以将序列按序号的奇偶性分为两组：

偶序列
奇序列

即一组由偶数序号组成，另外一组由奇数序号组成（注意数据长度为）

互质因子算法

Winograd算法

拉普拉斯变换

拉普拉斯变换

傅里叶变换是将函数分解到频率不同、幅值恒为1的单位圆上；拉普拉斯变换是将函数分解到频率幅值都在变化的圆上。因为拉普拉斯变换的基有两个变量，因此更灵活，适用范围更广。

关于快速傅里叶变换只是做了简单的介绍和梳理,详细内容等以后有时间再更新。

概念解析

酉空间

设为一个复数域上的线形空间，若在中定义了两个变量的内积（数量积），记作，且满足：

，其中是的共轭

，当且仅当时等号成立

,对任意

则称为酉空间(空间)，又称为内积空间，当为实数域时，此时的内积是可交换的，有限维的实酉空间也就是欧几里德空间。直观来说，酉空间就是将欧几里德空间的内积运算从实数域拓展到复数域。

完备性

一个向量空间具有完备性指空间中的任何柯西序列都收敛在该空间之内。。
柯西列

柯西列就是空间中元素构成的一个序列，并且这个序列在无穷远处两个元素之间的距离趋于零。准确的说，如果空间中有一个序列 ,当的时候，（即二者的距离趋零），则就是一个柯西列，也就是说完备性保证了取序列极限不会跑到空间外面去。一个不完备的例子就是有理数的集合，例如这个集合可以用柯西列的极限去逼近，而这个极限并不在有理数这个集合中，所以有理数集合是不完备的，而实数集合是完备的。
三角恒等式

维基百科：三角恒等式
内积空间

内积空间是线性代数里的基本概念，是增添了一个额外的结构的向量空间。这个额外的结构叫做内积或标量积。内积将一对向量与一个标量连接起来，允许我们严格地谈论向量的“夹角”和“长度”，并进一步谈论向量的正交性。内积空间由欧几里得空间抽象而来（内积是点积的抽象），这是泛函分析讨论的课题。
欧拉公式的证明

欧拉公式（50）的证明方式有两种，第一种是构造函数：

利用拉格朗日中值定理证明即可；第二种则是使用麦克劳林级数分别展开得到：
$e^{ix}=\sum_{i=0}^{\infty}\frac{(ix)^n}{n!}=1+ix-\frac{x^2}{2!}-i\frac{x^3}{3!}+\dots+i^n\frac{x^n}{n!}+\dots（\forall x \in \R）\\ cos(x)=\sum_{i=0}^{\infty}\frac{(-1)^n}{n!}x^{2n}=1-\frac{x^2}{2!}+\frac{x^4}{4!}+\dots+\frac{(-1)^n}{(2n)!}x^{2n}+\dots（\forall x \in \R）\\ isin(x)=i\sum_{i=0}^{\infty}\frac{(-1)^n}{(2n+1)!}x^{2n+1}=ix-i\frac{x^3}{3!}+i\frac{x^5}{5!}+\dots+i\frac{(-1)^n}{(2n+1)!}x^{2n+1}+\dots（\forall x \in \R）$
后两者的麦克劳林级数展开相加刚好等于前者的麦克劳林级数展开。在此只补充拉格朗日中值定理，具体证明过程较简单，不做详细展开。

拉格朗日中值定理:

如果函数，在闭区间上连续，在开区间内可微。则少存在一点，使下面的等式成立:

推论：

若函数的导数恒为零，则为常值函数，即。

傅里叶变换

概述

希尔伯特空间

傅立叶级数

三角函数系

周期与非周期下的傅里叶级数

周期函数的傅立叶级数

非周期函数的傅立叶积分

傅里叶变换

傅立叶变换的性质

离散傅里叶变换

快速傅里叶变换

Cooley-Tukey算法

基2时间抽取法

互质因子算法

Winograd算法

概念解析

你可能感兴趣的:(傅里叶变换)