首页 > 代码库 > 时域和频域变换之---傅里叶级数的数学推导
时域和频域变换之---傅里叶级数的数学推导
废话不多说先列提纲:
0.概述-需求分析-功能描述-受限和缺点改进+知识点预备
1.泰勒级数和傅里叶级数的本质区别,泰勒展开
2. 函数投影和向量正交
3.两个不变函数求导是本身e^x,sinx,cosx也是为什么要傅里叶转换的原因!
4.傅里叶技术推到过程
5.附录参考资料
0.有些时候,尤其是在图像处理中,矩阵运算数据量太大,特征提取量多,此时可以通过时域转频域来减少计算量,而且此转换不会损失数据完整性。
时域转频域的方法有周期函数用傅里叶技术,非周期函数(没有间断点的函数)用傅里叶转换,类似于直方图的分析。
两角和公式
tan(α+β)=(tanα+tanβ)/(1-tanαtanβ)
tan(α-β)=(tanα-tanβ)/(1+tanαtanβ)
cos(α+β)=cosαcosβ-sinαsinβ
cos(α-β)=cosαcosβ+sinαsinβ
sin(α+β)=sinαcosβ+cosαsinβ
sin(α-β)=sinαcosβ -cosαsinβ
积化和差
sinαsinβ =-[cos(α+β)-cos(α-β)] /2
cosαcosβ = [cos(α+β)+cos(α-β)]/2
sinαcosβ = [sin(α+β)+sin(α-β)]/2
cosαsinβ = [sin(α+β)-sin(α-β)]/2
二倍角公式
正弦
sin2A=2sinA·cosA
余弦
1.Cos2a=Cos^2(a)-Sin^2(a)
2.Cos2a=1-2Sin^2(a)
3.Cos2a=2Cos^2(a)-1
即Cos2a=Cos^2(a)-Sin^2(a)=2Cos^2(a)-1=1-2Sin^2(a)
正切
tan2A=(2tanA)/(1-tan^2(A))
平方关系:
sin^2(α)+cos^2(α)=1
1+tan^2(α)=sec^2(α)
1+cot^2(α)=csc^2(α)
平常针对不同条件的常用的两个公式
sin^2(α)+cos^2(α)=1
tan α *cot α=1
(1).
泰勒级数是求导函数组成的变化特征函数和;反应变化剧烈程度
傅里叶级数是频谱叠加的三角函数和;反应变化频率本质属性
(2).
从几何的角度来看。傅里叶告诉我们,f(x) 可以用下面这组由无限多个三角函数(包括常数)组成的“正交基”来展开,
傅里叶级数展开其实只是在做一个动作,那就是把函数“投影”到一系列由三角函数构成的“坐标轴”上。
(3).
正弦和余弦是二阶偏微分方程(含有电容等 元件的方程),而电容是可以隔直流通交流的
(4).
首先,隆重推出傅里叶级数的公式,不过这个东西属于“文物”级别的,诞生于19世纪初,因为傅里叶他老人家生于1768年,死于1830年。
但傅里叶级数在数论、组合数学、信号处理、概率论、统计学、密码学、声学、光学等领域都有着广泛的应用,这不由得让人肃然起敬。一打开《信号与系统》、《锁相环原理》等书籍,动不动就跳出一个“傅里叶级数”或“傅里叶变换”,弄一长串公式,让人云山雾罩。
如下就是傅里叶级数的公式:
不客气地说,这个公式可以说是像“臭婆娘的裹脚布——又臭又长”,而且来历相当蹊跷,不知那个傅里叶什么时候灵光乍现,把一个周期函数f(t)硬生生地写成这么一大堆东西。单看那个①式,就是把周期函数f(t)描述成一个常数系数a0、及1倍ω的sin和cos函数、2倍ω的sin和cos函数等、到n倍ω的sin和cos函数等一系列式子的和,且每项都有不同的系数,即An和Bn,至于这些系数,需要用积分来解得,即②③④式,不过为了积分方便,积分区间一般设为[-π, π],也相当一个周期T的宽度。
能否从数学的角度推导出此公式,以使傅里叶级数来得明白些,让我等能了解它的前世今生呢?下面来详细解释一下此公式的得出过程:
1、把一个周期函数表示成三角级数:
首先,周期函数是客观世界中周期运动的数学表述,如物体挂在弹簧上作简谐振动、单摆振动、无线电电子振荡器的电子振荡等,大多可以表述为:
f(x)=A sin(ωt+ψ)
这里t表示时间,A表示振幅,ω为角频率,ψ为初相(与考察时设置原点位置有关)。
然而,世界上许多周期信号并非正弦函数那么简单,如方波、三角波等。傅叶里就想,能否用一系列的三角函数An sin(nωt+ψ)之和来表示那个较复杂的周期函数f(t)呢?因为正弦函数sin可以说是最简单的周期函数了。于是,傅里叶写出下式:(关于傅里叶推导纯属猜想)
这里,t是变量,其他都是常数。与上面最简单的正弦周期函数相比,5式中多了一个n,且n从1到无穷大。这里f(t)是已知函数,也就是需要分解的原周期函数。从公式5来看,傅里叶是想把一个周期函数表示成许多正弦函数的线性叠加,这许许多多的正弦函数有着不同的幅度分量(即式中An)、有不同的周期或说是频率(是原周期函数的整数倍,即n)、有不同的初相角(即ψ),当然还有一项常数项(即A0)。要命的是,这个n是从1到无穷大,也就是是一个无穷级数。
应该说,傅里叶是一个天才,想得那么复杂。一般人不太会把一个简单的周期函数弄成这么一个复杂的表示式。但傅里叶认为,式子右边一大堆的函数,其实都是最简单的正弦函数,有利于后续的分析和计算。当然,这个式能否成立,关键是级数中的每一项都有一个未知系数,如A0、An等,如果能把这些系数求出来,那么5式就可以成立。当然在5式中,唯一已知的就是原周期函数f(t),那么只需用已知函数f(t)来表达出各项系数,上式就可以成立,也能计算了。
于是乎,傅里叶首先对式5作如下变形:
这样,公式5就可以写成如下公式6的形式:
这个公式6就是通常形式的三角级数,接下来的任务就是要把各项系数an和bn及a0用已知函数f(t)来表达出来。
2、三角函数的正交性:
这是为下一步傅里叶级数展开时所用积分的准备知识。一个三角函数系:1,cosx , sinx , cos2x , sin2x , … , cosnx , sinnx , … 如果这一堆函数(包括常数1)中任何两个不同函数的乘积在区间[-π, π]上的积分等于零,就说三角函数系在区间[-π, π]上正交,即有如下式子:
以上各式在区间[-π, π]的定积分均为0,第1第2式可视为三角函数cos和sin与1相乘的积分;第3-5式则为sin和cos的不同组合相乘的积分式。除了这5个式子外,不可能再有其他的组合了。注意,第4第5两个式中,k不能等于n,否则就不属于“三角函数系中任意两个不同函数”的定义了,变成同一函数的平方了。但第3式中,k与n可以相等,相等时也是二个不同函数。下面通过计算第4式的定积分来验证其正确性,第4式中二函数相乘可以写成:
可见在指定[-π, π]的区间里,该式的定积分为0。其他式也可逐一验证。
3、函数展开成傅里叶级数:
先把傅里叶级数表示为下式,即⑥式:
对⑥式从[-π, π]积分,得:
这就求得了第一个系数a0的表达式,即最上边傅里叶级数公式里的②式。接下来再求an和bn的表达式。用cos(kωt)乘⑥式的二边得:
至此,已经求得傅里叶级数中各系数的表达式,只要这些积分都存在,那么⑥式等号右侧所表示的傅里叶级数就能用来表达原函数f(t)。上述过程就是整个傅里叶级数的推导过程。事实上,如果能够写出⑥式,不难求出各个系数的表达式,关键是人们不会想到一个周期函数竟然可以用一些简单的正弦或余弦函数来表达,且这个表达式是一个无穷级数。这当然就是数学家傅里叶的天才之作了,我等只有拼命理解的份了。
综上,傅里叶级数的产生过程可以分为以下三步:
1、设想可以把一个周期函数f(t)通过最简单的一系列正弦函数来表示,即5式;
2、通过变形后用三角级数(含sin和cos)来表示;
3、通过积分,把各未知系数用f(t)的积分式来表达;
4、最后得到的4个表达式就是傅里叶级数公式。
在电子学中,傅里叶级数是一种频域分析工具,可以理解成一种复杂的周期波分解成直流项、基波(角频率为ω)和各次谐波(角频率为nω)的和,也就是级数中的各项。一般,随着n的增大,各次谐波的能量逐渐衰减,所以一般从级数中取前n项之和就可以很好接近原周期波形。这是傅里叶级数在电子学分析中的重要应用。
当然还有一疑问就是关于周期的取值l和t的不同,有的地方是w=l/2,有的就直接是w=an
:
0.傅里叶就是吧f(x)划分成不同频率三角函数的和
1.用内积法分解出每一个分量的系数
2.接地气的计算系数
在高数课本中是如上所求,但在《信号与系统》一书中傅里叶变换一节中直流分量为A/2,但其中对A和B的求解公式是一样的,这是怎么回事?所求结果肯定是不同的,
(5).
http://www.360doc.com/content/13/0328/12/202378_274443797.shtml;
http://blog.sina.com.cn/s/blog_57ad1bd20100txgs.html
http://www.360doc.com/content/13/0301/19/202378_268716061.shtml(几何公式)
http://blog.renren.com/share/343320656/15540620254(投影+正交)
http://www.360doc.com/content/11/0126/14/2355320_89141366.shtml(泰勒公式)
http://www.360doc.com/content/12/0120/16/99504_180542746.shtml(泰勒级数)
http://www.360doc.com/content/13/0205/02/1489589_264295813.shtml#(非正弦周期 的傅里叶级数)