基于MATLAB的谐波分析FFT

目录（1）Matlab6.5以上版本软件；..........................................................错误!未定义书签。绪论.....................................................................................................................................11公式分析及计算..............................................................................................................21.1傅里叶变换的原理.........................................................................................................21.2傅里叶变换的证明.........................................................................................................31.3周期信号的分解……………………………………………………………………………..31.4方波的分解…………………………………………………………………………………...52建模与仿真......................................................................................................................72.1建模...............................................................................................................................72.2仿真...............................................................................................................................83仿真结果分析................................................................................................................104小结................................................................................................................................11参考文献...........................................................................................................................13武汉理工大学《基础技能强化训练》课程设计说明书1绪论方波是一种非正弦曲线的波形，通常会于电子和讯号处理时出现。由于一般电子零件只有“高(1)”和“低(0)”两个值，方波就自然产生，所以理想方波只有“高”和“低”这两个值。电流的波形为矩形的电流即为方波电流。不论时间轴上下是不是对称的，只要是矩形就可叫方波，必要时，可加“对称”，“不对称”加以说明。而在现实世界，方波只有有限的带宽。因为方波可以快速从一个值转至另一个(即0→1或1→0)，所以方波就用作时钟讯号来准确地触发同步电路。但是如果用频率定义域来表示方波，就会出然一连串的谐波。所以方波可用相应频率的基波及其奇次谐波合成。在电路信号系统的分析中，随着电路规模的加大，微分方程的阶数以及联立后所得的方程的个数也随之加大，加上电器元件的多样化，这些都给解题运算分析电路系统带来了一定的困难。传统的计算机编程语言，如FORTRAN、C语言等，虽然都可以帮助计算，但在处理高阶微分方程和大规模的联立方程组的问题时大量的时间和精力都花在矩阵处理和图形的生成分析等繁琐易错的细节上。而MATLAB凭借其强大的矩阵运算能力、简便的绘图功能、可视化的仿真环境以及丰富的算法工具箱，已成为科研和工程技术人员的有力开发工具。利用MATLAB不仅可以简单快速的求解电路方程，同时，MAYLAB提供的Simulink工具还可以直接建立电路模型，随意改变模型的参数，并且还可以快速得到仿真拟结果，进一步省去了编程的步骤。MATLAB具有数值计算功能；图形处理及可视化功能；可视化建模及动态仿真功能等等。它给用户带来的是最直观，最简洁的程序开发环境。它的语言简洁紧凑，使用方便灵活，程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。同时，它的运算符也很丰富。由于MATLAB是用C语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短。它的程序的可移植性，基本上不做修改就可以在各种型号的计算机和操作系统上运行。本文应用MATLAB来验证定理：方波可用相应频率的基波及其奇次谐波合成。武汉理工大学《基础技能强化训练》课程设计说明书21公式分析及计算1.1傅里叶变换的原理任何具有性质周期为T的波函数()ft都可以表示为三角函数所构成的级数之和，即：10)sincos(21)(nnntnbtnaatf（1）其中：t为时间，为角频率。=T2（T为周期），第一项01a2为直流分量。图1方波所谓周期性函数的傅里叶变换（Fouriertransform）就是将周期性函数张凯成直流分量，基波和所有n次谐波的叠加。图1所示的方波可以写成函数形式：h(0≤t＜2T)()ft-h(-2T≤t＜0)在这里，h为常数2。很明显，此方波为奇函数，并且它没有常数项，同时，它是一个周期为T的函数，所以我们可以用傅里叶级数来表示这个函数。武汉理工大学《基础技能强化训练》课程设计说明书3()ft1])12sin[()121(4ntnnh（2）我们把它展开，可以得到：)7sin715sin513sin31(sin4)(tttthtf（3）1.2傅里叶变换的证明下面，我们要从数学角度来证明为什么公式（3）能成立。由于这是一个奇函数,常数项0a可以用积分函数计算出来：2202a=()0TTfxdtT所以其常数项不存在，即0a=0，下面开始计算na与nb:222=()cos()TTnafxntdtT020222()cos()cos()TThntdthntdtTT020222sin()sin()TThhntntTnTn0222=()sin()TTnbfxntdtT020222cos()cos()TThhntntTnTn2[1(1)]nhn由上式可知，当n=2，4，6，时，0nb；当n=1，3，5，时，4nhbn。然后我们将na与nb都带入公式（1），就可以得到公式（3）：)7sin715sin513sin31(sin4)(tttthtf如果我们取的项越多，就会越逼近原本的方波函数。1.3周期信号的分解武汉理工大学《基础技能强化训练》课程设计说明书4周期信号是定义在（-，）区间，每隔一定的时间T，按相同规律重复变化的信号，它可表示为f(t)=f(t+mT)式中m为任意整数。时间T称为该信号的重复周期，简称为周期。需要指出的是，只有当周期信号满足狄里赫利条件时，才能展开为傅里叶级数。狄里赫利条件是：1）函数在任意有限区间内连续，或只有有限个第一类间断点（当t从左或右趋于这个间断点时，函数有有限的左极限和右极限）2）在一周期内，函数有有限个极大值或极小值。设有周期信号f(t),它的周期是T，角频率=2F=T2，它可分解为(2-1)上式中系数an,bn称为傅里叶系数。为简便，式积分区间（t0,t0+T）取为（-2T，2T）或（0，T）。考虑到正、余弦函数的正交条件，可得傅里叶系数an=222TTTf(t)cos(nt)dt,n=0,1,2,……(2-2)bn=222TTTf(t)sin(nt)dt,n=0,1,2,……(2-3)式中T为函数f(t)的周期，=T2为角频率，由上述两式，傅里叶系数an和bn都tnbtbtbtnatataatfnnsin2sinsincos2coscos2)(21210武汉理工大学《基础技能强化训练》课程设计说明书5是n的函数，其中an是n的偶函数，即an=an；而bn是n的奇函数，既有bn=-bn将式（2-1）中同频率项合并，可写成如下形式f(t)=)2cos()cos(222110tAtAA……=)cos(210nnntnAA(2-4)式中00aAAn=22ba，n=1,2,……n=-arctan(nnab)如将式（2-4）的形式化为（2-1）的形式，他们系数之间的关系为a00AannnAcos,n=1,2,……bnnnAsin,式（2-4）表明，任何满足狄里赫利条件的周期函数可分解为直流和许多余弦（或正弦）分量。其中第一项20A是常数项，它是周期信号中所包涵的直流分量；式中第二项A)cos(11t称为基波或一次谐波，它的角频率与原周期信号相同，A1是基波振幅，1是基波初相角；式中第三项A)2cos(22t称为二次谐波，它的频率是基波频率的两倍，A2是二次谐波振幅，2是其初相角。以此类推，还有三次、四次、……谐波。一般而言，A)cos(nntn称为n次谐波，An是n次谐波的振幅，n是其初相角。式（2-4）表明，周期函数可以分解为各谐波分量。1.4方波的分解设方波信号f(t)的周期为T，宽度为2T，将其展开为傅里叶级数由式（2-2）和（2-3）可得22)cos()(2TTndttntfTa武汉理工大学《基础技能强化训练》课程设计说明书6=dttnTT)cos()1(202+dttnTT)cos()1(22002sin1220sin12TtnnTTtnnT考虑到T2，可得0na2002sin2sin12TTndttnTdttnTb02cos1220cos12TtnnTTtnnT,5,3,1n,6,4,240cos12nnnn将它们代入到式（2-1），得到信号的傅里叶级数展开式为tnttttfsinn15sin513sin31sin4)(，n=1，3，5，……它只含一、三、五…奇次谐波分量。下图中画出了一个周期的方波组成情况，由图可见，当它包含的谐波分量愈多时，波形就愈接近原来的方波信号tf（图中虚线所示），其均方误差愈小，还可以看出，频率较低的谐波，其振幅较大，他们组成方波的主体，而频率较高的高次谐波振幅较小，它们主要影响波形的细节，波形中所包含的高次谐波愈多，波形的边缘愈陡峭。（a）基波（b）基波“+”三次谐波武汉理工大学《基础技能强化训练》课程设计说明书7（c）基波+三次谐波+五次谐波（d）基波+三次谐波+五次谐波+七次谐波由图中还可以看出，合成波形所包含的谐波，除间断点附近外，它愈接近于