三角波-锯齿波-方波合成与分解MATLAB实现

连续周期信号的傅立叶级数分析及其MATLAB实现I连续周期信号的傅立叶级数分析及其MATLAB实现摘要现在，MATLAB已经发展成为适合多学科的大型软件，在世界各高校，MALAB已经成为线性代数、数值分析、数理统计、优化方法、自动控制、数字信号处理、动态系统仿真等高级课程的基本应用工具。本次课程设计则在深入研究连续时间信号傅里叶级数分析理论知识的基础上，利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，通过MATLAB编程进行图形功能仿真，从而实现连续时间周期信号频域分析的仿真波形，包括以下内容：用MATLAB实现周期信号的傅里叶级数分解与综合的波形；用MATLAB实现周期信号的单边频谱及双边频谱的波形与分析；用MATLAB实现典型周期信号的频谱的波形。关键词MATLAB图形处理傅里叶级数周期信号频谱连续周期信号的傅立叶级数分析及其MATLAB实现目录摘要..................................................................................................................................I1MATLAB..............................................................................................................................11.1简介..........................................................................................................................11.2特点.......................................................................................................................1２连续周期信号的傅立叶级数...............................................................................32.1连续时间周期信号的分解.....................................................................................31.1.１三角形式的傅里叶级数..............................................................................32.１.２指数形式的傅里叶级数............................................................................43连续周期信号的傅立叶级数分解与合成的MATLAB实现..............................53.1方波.........................................................................................................................53.2锯齿波......................................................................................................................93.3三角波...................................................................................................................134连续周期信号的频谱分析的MATLAB实现.........................................................21参考文献............................................................................................................................25附录................................................................................................................................26连续周期信号的傅立叶级数分析及其MATLAB实现11MATLAB1.1简介MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。MATLAB是一种科学计算软件，主要适用于矩阵计算及控制和信息处理领域的分析设计。它使用方便，输入简捷，运算效率高，内容丰富，并且很容易由用户自行扩展。Matlab当前已成为美国和其他发达国家在大学教学和科学研究中最常用而必不可少的工具MATLAB需要较多的高等数学的知识，要随着年级的增加，随着matlab在专业技术上应用才能逐渐深入掌握MATLAB在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。1.2特点第一、语言简洁，编程效率高因为MATLAB定义了专门用于矩阵运算的运算符，使得矩阵运算很简单。对于这些运算，几条语句即可代替数十条甚至上百条C语言及Fortran语言的功能。第二、交互性能好，使用方便在MATLAB的命令窗口中，输入一条命令，立即就能看到该命令的执行结果，体现了良好的交互性。交互方式减少了编程和调试程序的工作量，给使用者带来了极大的方便。因为不像使用C语言和Fortran语言那样，首先编写源程序，然后对其进行编译、连接、待形成可执行文件后，方可运行程序，得出结果。第三、强大的绘图功能，便于数据可视化连续周期信号的傅立叶级数分析及其MATLAB实现2MATLAB提供了两个层次的图形命令：一种是对图形句柄进行的低级图形命令，另一种是建立在低级图形命令之上的高级图形命令。利用MATLAB的高级图形命令可以轻而易举地绘制二维、三维乃至四维图形，并可进行图形和坐标的标识、视角和光照设计、色彩精细控制等等。第四、学科众多，领域广泛的应用工具箱其工具箱分为两大类：功能性工具箱和学科性工具箱。功能性工具箱主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。学科性工具箱专业性比较强，如控制系统工具箱、信号处理工具箱、神经网络工具箱、最优化工具箱、金融工具箱等，用户可以直接利用这些工具箱进行相关领域的科学研究。第五、开放性好，易于扩充除内部函数外，MATLAB的其他文件都是公开的、可读可改的源文件，体现了MATLAB的开放性特点。用户可修改源文件和加入自己的文件，甚至可以构造自己的工具箱。连续周期信号的傅立叶级数分析及其MATLAB实现3T2２连续周期信号的傅立叶级数频域分析法即傅里叶分析法，它是变换域分析法的基石。其中，傅里叶级数是变换域分析法的理论基础，傅里叶变换作为频域分析法的重要数学工具，具有明确的物理意义，在不同的领域得到广泛的应用。2.1连续时间周期信号的分解以高等数学的知识，任何周期为T的周期函数，在满足狄里赫利条件时，则该周期信号可以展开成傅里叶级数。傅里叶级数有三角形式和指数形式两种。1.1.１三角形式的傅里叶级数式中系数、称为傅里叶系数，可由下式求得2/)(0atf其中，为基波频率,为n次谐波频率。如果将式中同频率的正弦和余弦分量合并，则三角形式的傅里叶级数可表示为：上式中nanb222222012cos2sinTTTTTTnnaftdtTaftntdtTbftntdtT22n2101tcos2nnnAfAnt1,2,3,n0022,1,2,arctannnnnnnAaAabnab2324连续周期信号的傅立叶级数分析及其MATLAB实现4可以看出，傅里叶系数和都是或的函数，其中和是或的偶函数，即有；而和是或的奇函数，即有。2.１.２指数形式的傅里叶级数根据欧拉公式：并考虑和奇偶性可将改写为指数形式的傅里叶级数：,0,1,2,3,jntnnftFen即周期信号可分解为一系列不同频率的虚指数信号之和，式中nF称为傅里叶复系数，可由下式求得：221TTjntnFftedtT00cos,1,2,sinnnnnnnaAaAnbAnanbnnnanAnnnananbnnnnbnb][21)cos()()(nntnjtnjneetn)62(27)82()52(连续周期信号的傅立叶级数分析及其MATLAB实现53连续周期信号的傅立叶级数分解与合成的MATLAB实现3.1方波周期方波的傅里叶级数展开式：（3-1）MATLAB实现程序：clc;clearall;closeall;n=input('n=')T0=0.01;A=1;t1=[0:0.0001:0.05];z=square(2*pi*100*t);x=0;pi=3.1415926;w0=2*pi/T0;fori=1:nfw(i)=(2*i-1)*w0;a(i)=4*A/(pi*(2*i-1));y(i,:)=a(i)*sin(fw(i)*t);fai(i)=0;x=x+y(i,:);figure(1)a=1:100subplot(n,1,i);plot(t,y(i,:));xlabel('时间/s');ylabel('振幅')title([‘方波’,num2str(a(i)),'次谐波'])end;figure(2);subplot(2,1,1);连续周期信号的傅立叶级数分析及其MATLAB实现6plot(t,z);;xlabel('时间/s');axis([-0,0.05,-2,2])ylabel('振幅')title('方波')subplot(2,1,2);plot(t,x);xlabel('时间/s');ylabel('振幅');title([num2str(a(i)),'次谐波合成图'])生成图形：图3-1方波的谐波00.0050.010.0150.020.0250.030.0350.040.0450.05-202时间/s振幅方波1次谐波00.0050.010.0150.020.0250.030.0350.040.0450.05-0.500.5时间/s振幅方波2次谐波00.0050.010.0150.020.0250.030.0350.040.0450.05-0.500.5时间/s振幅方波3次谐波00.0050.010.0150.020.0250.030.0350.040.0450.05-0.200.2时间/s振幅方波4次谐波00.00