基于MATLAB的谐波分析FFT概要

《基于MATLAB的谐波分析FFT概要》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、目录（1）Matlab6.5以上版本软件；1绪论11公式分析及计算21.1傅里叶变换的原理21.2傅里叶变换的证明31.3周期信号的分解……………………………………………………………………………..31.4方波的分解…………………………………………………………………………………...52建模与仿真72.1建模72.2仿真83仿真结果分析104小结11参考文献1314武汉理工大学《基础技能强化训练》课程设计说明书14武汉理工大学《基础技能强化训练》课程设计说明书绪论方波是一种非正弦曲线的波形，通常会于电

2、子和讯号处理时出现。由于一般电子零件只有“高(1)”和“低(0)”两个值，方波就自然产生，所以理想方波只有“高”和“低”这两个值。电流的波形为矩形的电流即为方波电流。不论时间轴上下是不是对称的，只要是矩形就可叫方波，必要时，可加“对称”，“不对称”加以说明。而在现实世界，方波只有有限的带宽。因为方波可以快速从一个值转至另一个(即0→1或1→0)，所以方波就用作时钟讯号来准确地触发同步电路。但是如果用频率定义域来表示方波，就会出然一连串的谐波。所以方波可用相应频率的基波及其奇次谐波合成。在电路信号系统的

3、分析中，随着电路规模的加大，微分方程的阶数以及联立后所得的方程的个数也随之加大，加上电器元件的多样化，这些都给解题运算分析电路系统带来了一定的困难。传统的计算机编程语言，如FORTRAN、C语言等，虽然都可以帮助计算，但在处理高阶微分方程和大规模的联立方程组的问题时大量的时间和精力都花在矩阵处理和图形的生成分析等繁琐易错的细节上。而MATLAB凭借其强大的矩阵运算能力、简便的绘图功能、可视化的仿真环境以及丰富的算法工具箱，已成为科研和工程技术人员的有力开发工具。利用MATLAB不仅可以简单快速的求解电

4、路方程，同时，MAYLAB提供的Simulink工具还可以直接建立电路模型，随意改变模型的参数，并且还可以快速得到仿真拟结果，进一步省去了编程的步骤。MATLAB具有数值计算功能；图形处理及可视化功能；可视化建模及动态仿真功能等等。它给用户带来的是最直观，最简洁的程序开发环境。它的语言简洁紧凑，使用方便灵活，程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。同时，它的运算符也很丰富。由于MATLAB是用C语言编写的，MATLAB提供了和C语言几乎一样多的运算符，

5、灵活使用MATLAB的运算符将使程序变得极为简短。它的程序的可移植性，基本上不做修改就可以在各种型号的计算机和操作系统上运行。本文应用MATLAB来验证定理：方波可用相应频率的基波及其奇次谐波合成。14武汉理工大学《基础技能强化训练》课程设计说明书1公式分析及计算1.1傅里叶变换的原理任何具有性质周期为T的波函数都可以表示为三角函数所构成的级数之和，即：（1）其中：t为时间，为角频率。=（T为周期），第一项为直流分量。图1方波所谓周期性函数的傅里叶变换（Fouriertransform）就是将周期性函

6、数张凯成直流分量，基波和所有n次谐波的叠加。图1所示的方波可以写成函数形式：h(0≤t＜)-h(-≤t＜0)在这里，h为常数2。很明显，此方波为奇函数，并且它没有常数项，同时，它是一个周期为T的函数，所以我们可以用傅里叶级数来表示这个函数。14武汉理工大学《基础技能强化训练》课程设计说明书（2）我们把它展开，可以得到：（3）1.2傅里叶变换的证明下面，我们要从数学角度来证明为什么公式（3）能成立。由于这是一个奇函数,常数项可以用积分函数计算出来：所以其常数项不存在，即=0，下面开始计算与:由上式可知，

7、当n=2，4，6，时，；当n=1，3，5，时，。然后我们将与都带入公式（1），就可以得到公式（3）：如果我们取的项越多，就会越逼近原本的方波函数。1.3周期信号的分解14武汉理工大学《基础技能强化训练》课程设计说明书周期信号是定义在（-，）区间，每隔一定的时间T，按相同规律重复变化的信号，它可表示为f(t)=f(t+mT)式中m为任意整数。时间T称为该信号的重复周期，简称为周期。需要指出的是，只有当周期信号满足狄里赫利条件时，才能展开为傅里叶级数。狄里赫利条件是：1）函数在任意有限区间内连续，或只有有

8、限个第一类间断点（当t从左或右趋于这个间断点时，函数有有限的左极限和右极限）2）在一周期内，函数有有限个极大值或极小值。设有周期信号f(t),它的周期是T，角频率=2F=，它可分解为(2-1)上式中系数a,b称为傅里叶系数。为简便，式积分区间（t,t+T）取为（-，）或（0，T）。考虑到正、余弦函数的正交条件，可得傅里叶系数a=(t)cos(nt)dt,n=0,1,2,……(2-2)b=(t)sin(nt)dt,n=0,1,2,……(2-3)式中T为函数