[EE261学习笔记] 9.分布傅里叶变换（又称广义傅里叶变换）及δ函数

在进行傅里叶变换的时候，我们发现，因为傅里叶变换的积分需要被积函数黎曼可积，即仅当$f$可积时，才能计算出对应的$\mathcal{F}f$，然而在实际情况下会有很多函数是不可积的，比如之前提到的三角函数、矩形函数等等。因此就需要一个更普遍的傅里叶变换。

---

建立分布傅里叶变换的主要步骤如下（不仅适用于傅里叶变换推广到广义，也适用于其他）

1. 找到一个$\text{测试函数}$，我们将其记为 $\phi$ ，这个函数是你认为的可以用来解决你所遇到问题的一类最优函数

2. 我们引入一个定义，$\text{分布}$（又称广义函数），我们将其记为 $T$。首先，分布可以包含一切函数；其次，分布是一个更广泛的定义，它记录定义域内所有点的位置。之前提到的三角函数、矩形函数等，都可以看作一个分布。分布$T$是一个作用于测试函数$\phi$的$\text{线性泛函}$，即给定某$T$，$T(\phi )$将给出某个数。同时由线性可知$T({\phi }_1+{\phi }_2)=T({\phi }_1)+T({\phi }_2)$，$T(\alpha \phi )=\alpha T(\phi )$

3. 连续性。如果测试函数序列收敛${\phi }_n\to \phi$，则$T({\phi }_n)\to T(\phi )$。注意，函数序列的收敛是很难的，但是数字构成的序列收敛却很容易做

---

依据前面介绍的从一般情况到广义情况的步骤，我们开始讲解分布傅里叶变换：

1. 在傅里叶变换中，速降函数（schwartz function）就是我们需要的最优解。
   所谓速降函数，指的是这样一类函数$\phi (x)$
   a)$\phi (x)$无限可积，即$\frac{d^n}{d{x}^n}\phi (x)$存在
   b)它的任意阶导数下降的速度都快于$x$的任意次幂，即$|x{|}^n\left|\frac{d^n}{d{x}^n}\phi (x)\right|\to 0$，注意$m$和$n$是两个无关的变量
   速降函数的性质：
   a)设$\mathcal{S}$为速降函数的集合，若 $\phi \in \mathcal{S}$，则 $\mathcal{F}\phi \in \mathcal{S}$，${\mathcal{F}}^{-1}\phi \in \mathcal{S}$
   b)若 $\phi \in \mathcal{S}$，则 $\mathcal{F}{\mathcal{F}}^{-1}\phi ={\mathcal{F}}^{-1}\mathcal{F}\phi =\phi$
   这两个性质说明了：$\phi$的傅里叶变换、逆变换都是速降函数，都可以使积分收敛

2. 我们定义$T$与$\phi$的$\text{配对}$为$⟨T,\phi ⟩$。
   在分布傅里叶变换中
   $$⟨T,\phi ⟩={\int }_{-\infty }^{\infty }T(x)\phi (x)dx$$
   不难证明，配对是线性运算
   同时，为保证配对的收敛性，$T$必须是缓增分布（与速降函数相对应），实际上，在现实中运用分布傅里叶的情况基本都满足配对收敛

3. 该配对也满足连续性，具体证明比较麻烦，略。

---

接下来我们先给出分布傅里叶变换中配对的几个常见公式：

$$\begin{array}{rl}⟨\mathcal{F}T,\phi ⟩& ={\int }_{-\infty }^{\infty }\mathcal{F}T(y)\phi (y)dy\\ & ={\int }_{-\infty }^{\infty }\left({\int }_{-\infty }^{\infty }{e}^{-2\pi ixy}T(x)dx\right)\phi (y)dy\\ & ={\int }_{-\infty }^{\infty }\left({\int }_{-\infty }^{\infty }{e}^{-2\pi ixy}\phi (y)dy\right)T(x)dx\\ & ={\int }_{-\infty }^{\infty }\mathcal{F}\phi (x)T(x)dx\\ & =⟨T,\mathcal{F}\phi ⟩\end{array}$$

注意，在积分的过程中，$\mathcal{F}T(y)$本来不一定可以写成傅里叶积分的形式，但是我们假设它可以，从而推导出这个公式，再用这个公式来解决很多之前无法解决的问题。

同理可以证明以下常用的分布傅里叶变换中的配对公式：

$$⟨\mathcal{F}T,\phi ⟩=⟨T,\mathcal{F}\phi ⟩$$

$$⟨{\mathcal{F}}^{-1}T,\phi ⟩=⟨T,{\mathcal{F}}^{-1}\phi ⟩$$

$$⟨\mathcal{F}{\mathcal{F}}^{-1}T,\phi ⟩=⟨{\mathcal{F}}^{-1}\mathcal{F}T,\phi ⟩=⟨T,\phi ⟩$$

---

接下来我们运用这些公式来解决几个普通傅里叶变换无法解决的傅里叶变换问题

1. $\delta$ 的傅里叶变换
   $\delta$ 函数的相关定义见重侧δ函数
   第1.2例计算中均需要用到的重侧 $\delta$ 函数性质3
   $$\begin{array}{rl}\mathcal{F}\delta =⟨\mathcal{F}\delta ,\phi ⟩=⟨\delta ,\mathcal{F}\phi ⟩& =\mathcal{F}\phi (0)\\ & ={\int }_{-\infty }^{\infty }{e}^{-2\pi i0x}\phi (x)dx\\ & ={\int }_{-\infty }^{\infty }1\phi (x)dx\\ & =⟨1,\phi ⟩\\ & =1\end{array}$$
   这是一个非常直观的时域频域的联系，$\delta$ 函数表示集中于一点，而分布 $1$ 表示在区间内均匀分布。充分体现了时域的压缩就等于频域的拉伸。

2. ${\delta }_a$ 的傅里叶变换
   $$\begin{array}{rl}\mathcal{F}{\delta }_a=⟨\mathcal{F}{\delta }_a,\phi ⟩=⟨{\delta }_a,\mathcal{F}\phi ⟩& =\mathcal{F}\phi (a)\\ & ={\int }_{-\infty }^{\infty }{e}^{-2\pi iax}\phi (x)dx\\ & =⟨e^{-2\pi iax},\phi ⟩\\ & =e^{-2\pi iax}\end{array}$$

3. $e^{2\pi iax}$ 的傅里叶变换
   运用普通的傅里叶变换，积分不收敛，因此我们需要用分布傅里叶积分来计算，将$e^{-2\pi iax}$ 看作分布，于是就有：
   $$\begin{array}{rl}\mathcal{F}\left(e^{2\pi iax}\right)=⟨\mathcal{F}\left(e^{2\pi iax}\right),\phi ⟩& =⟨e^{2\pi iax},\mathcal{F}\phi ⟩\\ & ={\int }_{-\infty }^{\infty }{e}^{2\pi iax}\mathcal{F}\phi (x)dx\end{array}$$
   注意到上式相当于 $\phi$ 的傅里叶变换的逆变换，最后得到的是 $\phi$ 在 $a$ 点的值，即$\phi (a)$，而根据1中的计算，$\phi (a)$ 又等于 $⟨\delta ,\phi ⟩$，因此
   $$\mathcal{F}\left(e^{2\pi iax}\right)=\phi (a)=⟨{\delta }_a,\phi ⟩={\delta }_a$$
   我们将第2.3.两个例子结合在一起，可以得到
   $$\mathcal{F}\left(\mathcal{F}\left(e^{2\pi iax}\right)\right)=e^{-2\pi iax}$$
   这验证了我们在之前文章中提到的
   $$\mathcal{F}\mathcal{F}f=f^{-}$$
   也从侧面证明了分布傅里叶变换与普通傅里叶变换的统一性

4. $\sin 2\pi ax$ 与 $\cos 2\pi ax$ 的傅里叶变换
   首先，将函数用欧拉公式的变形（之前的文章中(2)(3)式）写出，有
   $$\sin 2\pi ax=\frac{1}{2i}\left(e^{2\pi iax}-e^{-2\pi iax}\right)$$
   $$\cos 2\pi ax=\frac{1}{2}\left(e^{2\pi iax}+e^{-2\pi iax}\right)$$
   然后利用例3，可以直接得出：
   $$\mathcal{F}(\sin 2\pi ax)=\frac{1}{2i}\left({\delta }_a-{\delta }_{-a}\right)$$
   $$\mathcal{F}(\cos 2\pi ax)=\frac{1}{2}\left({\delta }_a+{\delta }_{-a}\right)$$