傅里叶变换:只需要异想天开?
By 苏剑林 | 2014-04-25 | 44632位读者 |在对数学或物理进行事后分析,往往会发现一些奇怪的现象,也有可能得到一些更为深刻有趣的结果。比如本文所要谈及的傅里叶变换,可以由一种“异想天开”的思路得来。
洛朗展式 #
我们知道,在原点处形态良好的函数,可以展开为泰勒级数
$$f(x)=\sum_{n=0}^{\infty}a_n x^n$$
我们发现,上面的幂都是正的,为什么不能包含$x$的负数次幂呢?比如$\frac{\sin z}{z^2}$展开为
$$\frac{1}{z}-\frac{z}{6}+\frac{z^3}{120}\dots$$
显然也是一件合理的事情。于是,结合复变函数,我们得到解析函数的洛朗展式
$$f(z)=\sum_{n=-\infty}^{+\infty}a_n z^n$$
这是函数的双边展开。其中
$$a_n=\frac{1}{2\pi i}\int_{\gamma} \frac{f(z)}{z^{n+1}}dz$$
$\gamma$是积分曲线$|z|=\rho,\rho>0$。这个公式基于下面一个很显然的事实($\alpha\in\mathbb{Z}$)
$$\int_{\gamma} z^{\alpha}dz=\left\{ {\begin{array}{\cdot {20}{c}}
{2\pi i,\alpha=-1;}\\
{0,\alpha\neq -1.}
\end{array}} \right.$$
半整数幂级数 #
可是,还有一个问题,为什么不能包含分数次幂呢?诸如$\sqrt{z}$的函数如何展开呢?这说明,洛朗展式还是可以拓展的。为了举例说明,我们不妨把“半整数幂”也加进去,考虑:
$$f(z)=\sum_{n=-\infty}^{+\infty}a_n z^{n/2}$$
如何求各项系数?我们将它转化为熟悉的洛朗展式即可。设$z^{1/2}=\xi$,则
$$f(\xi^2)=\sum_{n=-\infty}^{+\infty}a_n \xi^n$$
洛朗展式我们知道怎么求系数了,即
$$\begin{aligned}
a_n=\frac{1}{2\pi i}\int_{\gamma}\frac{f(\xi^2)}{\xi^{n+1}}d\xi &=\frac{1}{2\pi i}\int_{2\gamma}\frac{f(z)}{z^{(n+1)/2}}dz^{1/2}\\
&=\frac{1}{4\pi i}\int_{2\gamma}\frac{f(z)}{z^{n/2+1}}dz
\end{aligned}$$
$\gamma$是以原点为圆心的一个圆(逆时针绕原点一周),$n\gamma$则表示逆时针绕原点$n$周。新级数是洛朗展式的推广之一。
一劳永逸 #
可是,疑问又来了,为什么不能有1/3整数幂呢?为什么不能有无理数次幂呢?这类“异想天开”可以是无穷尽的。于是乎,为了避免被进一步提问,我们干脆把所有实数幂都考虑进去,这也就是本节的“一劳永逸”的意思。函数
$$f(z)=\int_{-\infty}^{+\infty}a(x)z^x dx$$
这包含了上面的两种形式。(离散的$a_n$对应带有狄拉克函数$\delta(x)$的$a(x)$。)。下面我们来推导$a(x)$的表达式。我们把求和离散化,从原点出发,以$\Delta x$为步长进行求和
$$f(z)=\lim_{\Delta x \to 0} \sum_{n=-\infty}^{+\infty}a(n\Delta x)\Delta x z^{n\Delta x}$$
设$z^{\Delta x}=u$,代入
$$f(\xi^{1/\Delta x})=\sum_{n=-\infty}^{+\infty}a(n\Delta x)\Delta x \xi^{n}$$
同样用洛朗展式的方法求系数,得
$$\begin{aligned}a(n\Delta x)\Delta x=\frac{1}{2\pi i}\int_{\gamma}\frac{f(\xi^{1/\Delta x})}{\xi^{n+1}}d\xi &=\frac{1}{2\pi i}\int_{\left(\frac{1}{\Delta x}\right)\gamma}\frac{f(z)}{z^{(n+1)\Delta x}}dz^{\Delta x}\\
&=\frac{1}{2\pi i}\int_{\left(\frac{1}{\Delta x}\right)\gamma}\frac{\Delta x f(z)}{z^{n\Delta x+1}}dz
\end{aligned}$$
此即
$$\begin{aligned}a(x)=\frac{1}{2\pi i}\int_{\left(\frac{1}{\Delta x}\right)\gamma}\frac{f(z)}{z^{x+1}}dz
\end{aligned}$$
积分路径为绕原点的圆无穷多圈。由于变换之间可逆,这成为了两个函数之间的相互变换。
傅里叶变换 #
现在我们换个熟悉的记号,设$z=e^{-i\omega}$,记$a(x)$换成记号$ f(x)$,对应的$f(z)$则换成$F(\omega)$,有
$$F(\omega)=\int_{-\infty}^{+\infty}f(x)e^{-i\omega x} dx$$
以及
$$f(x)=-\frac{1}{2\pi i}\int_{-\infty}^{+\infty} \frac{F(\omega)}{e^{-i\omega(x+1)}}de^{-i\omega}=\frac{1}{2\pi }\int_{-\infty}^{+\infty} F(\omega) e^{i\omega x}d\omega$$
我们取的积分曲线$\gamma$是逆时针的,但是$e^{-i\omega x}$是顺时针的,所以要加个负号。这就得到了傅里叶变换及其逆变换。
稍微总结 #
读者最大的疑问应该是为什么积分限是负无穷到正无穷,而不是0到正无穷?这个确实不容易解释清楚,严格推导傅里叶逆变换还是需要《数学物理方法》上面的步骤。当然,勉强来说,从本文的角度,还是可以给积分限一个解释的(解释而不是证明)。
如果积分上限选取适当,那么积分下限时可以随意选取的,比如积分下限取0,那么积分上限就必须取$2\pi N,N\in \mathbb{Z}$且$N\to \infty$。但是,趋于无穷时,加上个限制($N\in\mathbb{Z}$而不是$N\in \mathbb{R}$)总是不舒服的,由于正无穷具有不确定性,我们干脆在下限也来个不确定性(负无穷),这样子,由于没有确定起点,也就没有必要给终点加上个限制了。于是积分上下限就可以写成正负无穷,至于区域无穷的部分的瑕疵,在积分的时候自然地被抹平了。(PS:这段话相当含糊,但是为了使文章不至于陷入繁琐的技术细节中,只能说到这个份上。事实上,基于数学分析中严格的瑕积分理论,可以给上面这段话一个证明。这主要用到诸如$\int_0^{\infty} f(x)\cos(\omega x)dx$之类的积分技巧。希望深究的朋友,不妨尝试一下?)
事实上,上面的推导基本上都是不严格的,本文的目的是想通过这个思路给大家一种对傅里叶变换的“比较自然”的认识,并且提供大家将两种看似不相关的东西(傅里叶变换和洛朗展式)联系起来的思路。通过复变函数,可以将很多数学分支巧妙地联系起来。这似乎是设计者的一个法则:和谐统一。如果真是这样子的话,我们将有更深层次的理由去探索数学、探索科学,那就是——为了美!
转载到请包括本文地址:https://spaces.ac.cn/archives/2555
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Apr. 25, 2014). 《傅里叶变换:只需要异想天开? 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/2555
@online{kexuefm-2555,
title={傅里叶变换:只需要异想天开?},
author={苏剑林},
year={2014},
month={Apr},
url={\url{https://spaces.ac.cn/archives/2555}},
}
January 5th, 2015
楼主说的不就是Z变换吗?只不过连续化了而已
谢谢这位读者的提醒。我之前只是就像文中所说的那样研究下去得到的结果,并不知道有专门研究这个变换的内容。
我抽时间读读Z变换的资料吧,有收获再和大家分享。
July 16th, 2015
一旦连续化,洛朗展开就成Fourier变换了,当初怎看不出来,楼主辛苦啦。
ps,那个,一劳永逸里,文字第六行,设z,,,,,(有问题)
January 18th, 2023
苏神,你能不能顺手推一个向量值函数的傅里叶变换方法?现在流行的多维傅里叶变换,都是定义域是多维的,然后信号值域还是一维的,如果能有信号值域是向量值的版本是不是数学上就更加完美了?
多维傅立叶变换不就是一维的平凡推广?
目前教材上的多维傅里叶变换,被变换的函数定义域是多维的,值域仍然是一个scalar,比如图像领域中的2D-FFT。。。
我想要的是向量值函数的傅里叶变换,例如$\pmb{x(t)}=(\cos t,\sin t)$,它的自变量是t,应变量是$(\cos t,\sin t)$,求这种信号的傅里叶变换。
$$\int \boldsymbol{f}(\boldsymbol{x})e^{i\boldsymbol{\omega}\cdot\boldsymbol{x}}d\boldsymbol{x}$$
这个运算对于向量值函数也是可以正常操作的。对于你举的例子$(\cos t,\sin t)$,我理解就是每个分量独立进行傅立叶变换。
November 3rd, 2023
确实可以,我试了一下,如果加入了类似复数向量模长这种概念,好像他也满足很多经典的傅里叶变换的性质,甚至于帕塞瓦尔定理也能轻易地推广开来。事实上好像是可以定义$\Bbb R^m \rightarrow \Bbb R^n$的傅里叶变换的。
January 13th, 2025
工程中使用傅里叶变换的时候其实有很多限制,有时候确实有很多缺陷。比如一个信号是$f(t)+f(t-t_0)$,他的傅里叶变换是$F(\omega)+F(\omega)e^{-\text j\omega}$,就是有时移效应,导致它的频谱和功率谱都不是线性叠加的。在有些情况下确实比较难弄,比如我想通过某种统计量描述信号的形状,或者发现曲线某种拓扑上的规律的时候,往往希望不要将位移叠加到相位上去,而是在谱值上能直观线性叠加,尽管能量最后是满足帕塞瓦尔公式是线性叠加的,但是遗憾的是谱不是。
这个问题还不能更换其他正交基来解决。本质上是,如果把一个信号看成是一个无穷长度的向量,那么时移相当于在向量的前面几个维数上补了零,相当于做了一个换轴的正交变换。这个正交变换后的向量在原来的基的投影肯定会变,只不过傅里叶基刚好是$\sin$和$\cos$,所以引入了所谓的相位,由于同一个位移对应不同的频率时移动的相位不一样,所以出现了所谓的$e^{-\text j\omega}$这个东西。换其他的正交基,本质上还是投影的问题,一样不会变。
那可不可以把所谓的位移丢出去,变成另一个参数呢?比如加窗傅里叶变换(盖波变换),短时傅里叶变换、小波变换等等,基本搞的就是那一套东西。以小波变换为例,确实有$W_{f(t-\tau)}(a,b)=W_{f(t)}(a,b-\tau)$,但实际上还是取决于小波基是个什么东西,最后总是另一个空间矢量的叠加,似乎在这个角度不怎么适合分析形状?而且不管是哪套,本质上等于叠加上去的加窗傅里叶变换,且不说这类变换本身受到时频测不准原理的限制,就单单说把信号的某一个点加窗然后抠下来做一个变换这件事情来说,本质上是将一个信号$f(t),t\in(-\infty,+\infty)$变成了$(-\infty,a)$和$(b,+\infty)$为$0$,$[a,b]$上为$f(t)$的东西,就算是分析信号,这个窗还在两个边缘引入了两个阶跃信号在里面,可以认为每一个区间信号分析都是不精确的。非常好玩的一点就是,如果你单独去看某一个区间的信号做分析,按照这个说法,就是不精确的,但是连在一起看,又变得精确了,因为叠加在一起刚好正负抵消,就是不能抽出来看。。。似乎一抽出来就是不准的?这是不是也能看成是“观察不准原理”?
January 13th, 2025
瞬时频率的思路还要求信号是可导的,又无法描述复杂的信号结构,好像还挺麻烦的。。。
January 13th, 2025
瞬时频率的思路还要求信号是可导的,又无法描述复杂的信号结构,好像还挺麻烦的。。。