基于MATLAB的谐波分析FFT

目录

（1）Matlab6.5以上版本软件； ..........................................................错误!未定义书签。绪论. (1)

1 公式分析及计算 (2)

1.1傅里叶变换的原理 (2)

1.2傅里叶变换的证明 (3)

1.3 周期信号的分解 (3)

1.4 方波的分解 (5)

2 建模与仿真 (7)

2.1建模 (7)

2.2仿真 (8)

3 仿真结果分析 (10)

4 小结 (11)

参考文献 (13)

绪论

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

在电路信号系统的分析中，随着电路规模的加大，微分方程的阶数以及联立后所得的方程的个数也随之加大，加上电器元件的多样化，这些都给解题运算分析电路系统带来了一定的困难。传统的计算机编程语言，如FORTRAN、C语言等，虽然都可以帮助计算，但在处理高阶微分方程和大规模的联立方程组的问题时大量的时间和精力都花在矩阵处理和图形的生成分析等繁琐易错的细节上。而MATLAB凭借其强大的矩阵运算能力、简便的绘图功能、可视化的仿真环境以及丰富的算法工具箱，已成为科研和工程技术人员的有力开发工具。利用MATLAB不仅可以简单快速的求解电路方程，同时，MAYLAB提供的Simulink工具还可以直接建立电路模型，随意改变模型的参数，并且还可以快速得到仿真拟结果，进一步省去了编程的步骤。MATLAB具有数值计算功能；图形处理及可视化功能；可视化建模及动态仿真功能等等。它给用户带来的是最直观，最简洁的程序开发环境。它的语言简洁紧凑，使用方便灵活，程序书写形式自由，利用起丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。同时，它的运算符也很丰富。由于MATLAB是用C

语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短。它的程序的可移植性，基本上不做修改就可以在各种型号的计算机和操作系统上运行。

本文应用MATLAB来验证定理：方波可用相应频率的基波及其奇次谐波合成。

1 公式分析及计算

1.1傅里叶变换的原理

任何具有性质周期为T 的波函数()f t 都可以表示为三角函数所构成的级数之和，即：

∑∞

=++=1

0)sin cos (21

)(n n n t n b t n a a t f ωω （1）

其中：t 为时间，ω为角频率。ω=T

π

2（T 为周期），第一项01a 2为直流分量。

图1 方波

所谓周期性函数的傅里叶变换（Fourier transform ）就是将周期性函数张凯成直流分量，基波和所有n 次谐波的叠加。

图1所示的方波可以写成函数形式：

h (0≤t ＜2

T )

()f t =

-h (-2

T ≤t ＜0)

在这里，h 为常数2。很明显，此方波为奇函数，并且它没有常数项，同时，它是一个周期为T 的函数，所以我们可以用傅里叶级数来表示这个函数。

()f t =

∑∞

=--1

])12sin[()1

21

(

4n t n n h

ωπ

（2） 我们把它展开，可以得到：

)7sin 7

1

5sin 513sin 31(sin 4)( ++++=

t t t t h

t f ωωωωπ （3）

1.2傅里叶变换的证明

下面，我们要从数学角度来证明为什么公式（3）能成立。由于这是一个奇函数,常数项0a 可以用积分函数计算出来：

2

202a =()0T

T

f x dt T -=⎰

所以其常数项不存在，即0a =0，下面开始计算n a 与n b :

22

2=()cos()T

T n a f x n t dt T ω-⎰

020

222()cos()cos()T

T h n t dt h n t dt T T ωω-=-+⎰⎰

[][]0

202

22sin()sin()T T h h n t n t T n T n ωωωω-=-+

0=

22

2=()sin()T

T n b f x n t dt T ω-⎰

[][]0202

22cos()cos()T

T h h n t n t T n T n ωωωω-=-+ 2[1(1)]n h

n π

=

-- 由上式可知，当n =2，4，6

，

时，0n b =；当n =1，3，5

，时，4n h

b n π

=

。然后我们将n a 与n b 都带入公式（1），就可以得到公式（3）：

)7sin 7

1

5sin 513sin 31(sin 4)( ++++=

t t t t h

t f ωωωωπ 如果我们取的项越多，就会越逼近原本的方波函数。

1.3 周期信号的分解