Fourier分析基础(一)——Fourier级数

前言

傅立叶分析的作用是把一个函数变成一堆三角函数的和的形式,也就是分解。首先引入的是傅立叶级数,Fourier级数的作用是把函数变为可数无限个三角函数的和,而且这些三角函数的频率都是某个基频的整数倍。如果这个基频无限趋近于0,那么在极限的情况下这函数的参数(频率)就连续了,将连续时域函数映射到连续的频域函数的变换就是标准的傅立叶变换。

由于工程采集的信号大多都是离散的,把时域离散化以后不可能在得到连续的频域函数,所以在频域上也不连续了,这种离散时域序列到离散频域序列的变换称之为离散傅立叶变换(DFT),然后有人开发出了快速计算的快速傅立叶变换(FFT)。

以上介绍的每一种Fourier变换都有其逆变换。

Fourier级数

考虑一下,假设存在2个序列clip_image002[6]clip_image004[4],还有一个数字clip_image006[6]

现在有一个时域上变化的函数clip_image008[4]。这个函数可以表达为如下的形式:

clip_image010[4]

这就是傅立叶级数,傅立叶变换最基础的形式。上式中和式的形式或许并不直观,如果画出一部分或许会直观一些。

如图所示,是从圆频率为1rad/s~10rad/s的正弦信号的合集。

clip_image011[4]

clip_image013[4]空间上的Fourier级数展开

考虑傅立叶级数,首先考虑周期函数在clip_image015[4]上的展开,但是在展开之前,需要做几个计算和证明。

计算clip_image017[4]clip_image019[4]就不计算了,这个是一样的)。

积分过程略,得到。

clip_image021[4]

这样就求出了clip_image023[4]clip_image025[6]上的范数,就是clip_image027[4]

那么正交基不妨使用clip_image029[4]构成。

然而我并不知道这个集合是不是正交的基,需要证明啊!

下面证明如下的积分关系成立:

clip_image031[4]

clip_image033[4]

clip_image035[4]

其中,clip_image037[4],可以的很轻松得到:2,4,5的正确性,关键在于1,3,6的证明上。

证明了以上六个等式,也就等同证明了clip_image038[4]clip_image025[7]空间的一组正交基。

首先我们都造:

clip_image040[4]

clip_image042[4]

辣么两个加一下就得到:

clip_image044[4]

减一下就得到:

clip_image046[4]

clip_image048[4]的时候,上面那俩货clip_image050[4]上积分是0可以用肉眼看出。

看不出的……呵呵

最后证明1式,首先假设clip_image052[4],很轻松就证明了,然后clip_image054[4]的时候,函数是奇函数,所以在相对原点对称的区间上积分是0。

证完收工。

现在我们搞到了这么一组正交基:

clip_image056[4]

下面就要用它分解函数了,也就是计算clip_image058[4]这些系数。

clip_image060[4]

这是个求平均值的操作应该很好理解。

clip_image062[4]

clip_image064[4]

以上两个公式应该也没有理解上的困难,但是有两个细节,第一个,为什么clip_image066[4]的系数使用和直流量(或者叫平均值)一样的clip_image068[4],而不是使用clip_image070[4],因为这个clip_image006[7]本质上是和基clip_image072[4]做内积的时候,clip_image074[4]时候的特殊情况。

其次,为何前面的系数是clip_image076[6],其实是为了匹配形式:

clip_image077[4]

我们的三角函数基都是带clip_image079[8]的,如果需要不带,那么系数就要乘以clip_image079[9],本来求系数很干脆:
clip_image081[4]

clip_image083[4]

这样表达式就要写成:

clip_image085[4]

如果让新的系数clip_image087[4]clip_image089[4]成为系数,那么就可以得到直接乘sin和cos的简洁形式,这个时候积分式变为:

clip_image091[4]

clip_image093[4]

两个clip_image079[10]相乘以后就得到了clip_image076[7]

clip_image095[4]空间上的Fourier级数展开

如果我们把上面那个宽度为clip_image097[4]的区域平移会发生什么?

答:毛线都不会发生。

那么我们需要改变什么呢?

答:改一下积分的上下限就可以了。以clip_image002[7]的求法为例,原来的式子需要改为:

clip_image099[4]

其它的以此类推都一样,没有任何变化。

之所以敢这么搞,是因为推导出正交基的六个等式的对错没有因为增加了这个clip_image101[4]发生变化,其引理自然也不变化。

clip_image103[4]空间上的Fourier级数展开

这个变化好像就在问:如果区间拉伸了怎么办?

函数的定义域(区间)拉伸了,那么我们把基也一起拉伸就OK了。

那么,问题来了:怎么拉伸呢?

很简单,先直接抻长了:

clip_image104[4]

然后保证范数是1就可以了。

范数的计算和上面范数的计算几乎是一毛一样的,所以在此略过直接说结果:

clip_image106[4]空间上满足:

clip_image108[4]

的参数的求取公式如下:

clip_image110[4]

clip_image112[4]

clip_image114[4]

Fourier级数的最终自由形态:clip_image117[4]空间上的级数展开

有了上面两份推导,现在可以正式得出更加一般的展开:

clip_image119[4]

这个时候,区间的长度是clip_image121[4],等同于刚才clip_image123[4],那么带入得到从参数换为函数的变换如上所示。

这个时候参数的求取是这样的:

clip_image125[4]

clip_image127[4]

clip_image129[4]

这样就结束了,我们得到了一种在任意有限区间对函数进行展开的方法,这样我们就可以开发出一些更加奇葩的东西,比如对一个函数使用Fourier基进行分段拟合。

Trick:

在对函数进行人肉傅立叶级数展开的时候,如果区间关于原点对称,首先观察函数的形式,如果是奇函数,那么cos的部分(含直流分量)就可以不看了:一定没有~

如果是偶函数,那么sin的部分就可以不用看了。

大道至简:Fourier级数的复指数形式

由著名的欧拉公式:clip_image132[4]

其中clip_image134[4],在许多场合一般使用j代替i,但是实际上一般定义clip_image136[4],但是由于clip_image138[4]的两个共轭的结果就是clip_image140[4],所以上述关系一般用不到,只是在做一些说明的时候为了方便使用j代替-i。

我们不妨设clip_image143[4],那么就可以得到Fourier的复指数形式:

clip_image145[4]

从三个公式变成一个公式,果然身心都愉悦了好多呢,但是这么变换为什么是有效的呢?

以下证明分两步走,第一步,证明为什么欧拉公式形状这么诡异但是是对的。第二,为什么能用复数吧sin和cos结合起来。

关于欧拉公式的证明:

首先把指数函数和两个三角函数进行Taylor展开:

clip_image147[4]

clip_image149[4]

clip_image151[4]

然后令clip_image153[4]可以推导出下式:

clip_image155[7]

就这么愉快的证明了。

那么为什么复数可以用来把两个函数糅合到一起呢?

因为在复平面上,实数轴和虚数轴是垂直的。

而且,我们以三角函数过零,且导数大于0的点为相位的0点,那么在clip_image157[7]的时候,sin的相位是0,而cos的相位是90°,也就是clip_image159[7]

那么,如果sin的相位是-90°,那么cos的相位就是0°,这是一个相对的关系,并无绝对。我们继续分析,其实clip_image161[10]是一条三维空间的螺旋线:

clip_image162[7]

除了有在复平面上能直接表示相位的良好性质以外,clip_image161[11]还是一个合格的复周期函数呢~那么我们不妨抛弃从旧版傅立叶到新版傅立叶的映射,重构整个证明吧!

现在我们新建规范正交基:clip_image165[7],我们来证明它是规范正交的。

根据内积的定义好像clip_image167[7],但是在复数下这是错误的!

复数下面其实是:

clip_image169[7]

这样就可以用共轭消除虚数部分了,得到的范数才是真实的。

话说指数函数真的是一个专业对抗求导积分500年的函数,它的积分更容易,在此不证明了,直接给出结论:它的范数是clip_image171[10],所以上面给出的正交基是规范的。

那么还不知道正交不正交呢(然而我就急着叫他规范正交基了)。

clip_image173[7]

clip_image175[7]的时候可以得到函数积分为0,如果clip_image177[7]的时候就反常了,然而并不是很难解,直接clip_image179[7]常数项积分谁都会,得到clip_image171[11],就是上面的结论。(一不小心还是证明了不想写的东西)

其实这个时候傅立叶级数已经脱胎换骨了,从根本上换了一套基,也就从原来的实数领域推广成了更为广阔的复数领域,但是形式却更加的简洁了。

原文地址:https://www.cnblogs.com/TsingJyujing/p/4951596.html