系统建模与设计报告课程设计-基于方波信号的分解与合成.doc

湖南文理学院系统建模与设计报告 专业班级： 电信10101班 学生姓名： 学生学号： 指导教师： 设计时间： 2012-12-20 基于matlab的信号分解与合成一、课程设计目的1、学习MATLAB 软件的使用.2、使学生掌握利用工具软件来实现信号系统基本概念、基本原理的方法。二、基本要求 掌握数字信号处理的基本概念、基本理论和基本方法； 学会 MATLAB 的使用，掌握 MATLAB 的程序设计方法； 学会用 MATLAB 对信号进行分析和处理； 信号的各参数需由键盘输入，输入不同参数即可得不同的x(t) 和x(n)； 撰写课程设计论文，用数字信号处理基本理论分析结果。三、设计方法与步骤1、信号分解为正交函数 设有个函数在区间构成一个正交函数空间。将任一函数用这个正交函数的线性组合来近似，可表示为 （1-1）这里的问题是：如何选择才能得到最佳近似。显然，应选取个系数使实际函数与近似函数之间误差在区间内为最小。这里“误差最小”不是平均误差最小，因为在平均误差最小甚至等于零的情况下，也可能有较大的正误差和负误差在平均过程中相互抵消，以致不能正确反映两函数的近似程度。通常选择误差均方值（或称为方均值）最小，这时，可以认为已经得到了最好的近似。误差的均方值也称为均方误差，用符号表示 （1-2）在中，为球的使均方误差最小的第个系数，必须使即 （1-3）展开上式的被积函数，注意到有序号不同的正交函数相乘的各项，其积分均为零，而且所有不包含的各项对求导也等于零。这样，式（1-3）中只有两项不为零，它可以写为交换微分与积分次序，得于是可求得 （1-4）式中 （1-5）这就是满足最小均方误差的条件下，式（1-1）中各系数的表达式。此时，能获得最佳近似。 当按式（1-4）选取系数时，将代入到式（1-2），可以得到最佳近似条件下的均方误差为 考虑到，得 （1-6）利用上式可直接求得在给定项数的条件下的最小均方误差。 有均方误差的定义式（1-4）可见，由于函数平方后再积分，因而不可能为负，即恒有0。由式（1-5）可见，在用正交函数去近似（或逼近）时，所取的项数越多，即愈大，则均方误差愈小。当时，。由式（1-6）可得，如，则有 （1-7）式（1-7）称为帕斯瓦尔（Parseval）方程。 如果信号是电压或电流，那么，式（1-7）等号左端就是在区间信号的能量，等号右边是在区间信号各正交分量的能量之和。式（1-7）表明：在区间信号所含能量恒等于此信号在完备函数集中各正交分量能量的总和。与此相反，如果信号在正交函数中的各正交分量能量的总和小于信号本身的能量，这时式（1-7）不成立，该正交函数集不完备。 这样，当时，均方误差，式（1-1）可写成 （1-8）即函数在区间可分解为无穷多项正交函数之和。2、周期信号的信号分解与合成2.1周期连续信号的特点周期连续信号有如下特点：（1）满足，m是整数，是周期。从波形上看，有一个时间跨度为的基本波形，其余的是该基本波形经平移的整数倍后的重新拷贝。（2）在一个周期内的积分，其值与积分的起点和终点无关，即有 （3）将周期信号展开成傅里叶级数具有的以下显著优点是：三角函数和指数函数是自然界中最常见、最基本的函数。三角函数和复指数函数是间谐函数，用它们表示时间信号，就自然地建立了时间和频率这两个基本物理量之间的关系。间谐信号较其他信号更容易产生和处理。三角函数（或指数函数）信号通过线性时不变系统后，仍为三角函数（或指数函数），其重复频率不变，只是幅度和相位有变化。线性时不变系统对三角函数（或指数函数）信号的响应可以很方便地求的。很多系统（例如滤波器、信息传输等）的特性主要是由其频域特性来描述的，因此常常更需要知道的并不是这些系统的冲激响应，而是其冲激响应所对应的频率特性。时域中的卷积运算在频域会转化为乘积运算，从而找到了计算卷积的一种新方法，这可使时域中难以实现的卷积运算求解便于实现。周期信号当满足狄里赫利（Dirichlet）条件时可以展开成傅里叶级数。傅里叶级数分三角形式和指数形式两种。狄里赫利条件如下：（1）在一个周期内，是绝对可积，即。（2）在一个周期内，的最大值和最小值的数目是有限个。（3）在一个周期内，只有有限个间断点，而且在这些间断点上，函数值必须是有限个。2.2 周期为T的信号的三角形式的傅里叶级数表示的一般形式设有周期信号，它的周期为T，角频率，则的三角傅里叶级数表示的一般形式为 （2-1） 其中 可以写成更紧凑的和式为：式（2-1）中的系数、称为傅里叶系数，为在函数中的分量（相对大小）；为在函数中的分量，它可由式（1-4）求得。为简便，式（1-4）的积分区间取为或。考虑到正、余弦函数的正交条件，由式（1-4），可得傅立叶系数 （2-2） 周期信号也可分解为一系列余弦信号，即： （2-3）其中 式（2-3）表明，任何满足狄里赫利条件的周期函数都可以分解为直流和许多余弦（或正弦）分量。其中第一项是常数项，它是周期信号中所包含的直流分量；式中第二项称为基波或一次谐波，它的角频率和原周期信号相同，是基波振幅，是基波初相角；式中第三项称为二次谐波，它的频率是基波频率的二倍，是二次谐波的振幅，是其初相角。以此类推，还有四三次、四次、谐波。一般而言，称为次谐波，是次谐波的振幅，是其初相角。式（2-3）表明，周期信号可以分解为各次谐波分量。 式(2-2)表示周期信号可以分解成直流分量和各次谐波分量的叠加，用直流分量和各次谐波分量代替原来的周期信号，原则上应该是无穷多项的叠加，实际应用中只取其中的前项，产生的误差函数用en（t）来表示 （2-4）另外一个衡量误差大小的函数为方均误差【10】： （2-5）3、 MATLAB的仿真实现MATLAB是目前世界上最流行的、应用最广泛的工程计算和仿真软件，它将计算、可视化和编程等功能同时集于一个易于开发的环境。MATLAB是Matrix Laboratory的缩写，是一个包含众多工程计算和仿真的庞大系统。MATLAB是一个交互式开发系统，其基本数据要素是矩阵。语法规则简单，适应于专业科技人员的思维方式和书写习惯；它用解释方式工作，编写程序和运行同步，键入程序立即得出结果，因此人机交互更加简洁和智能化；而且MATLAB可适用于多种平台，随着计算机软、硬件的更新而及时升级,shide编程和调试效率大大提高3.1基于MATLAB方波信号的分解与合成现以周期为T、幅值为1的方波信号为例 方波信号的分解与合成 图1 周期为T的方波图由式（2-2）可得 考虑到，可得 将它们代入（2-1）式，可得图1所示的方波信号的傅立叶级数展开式为 它只含一、三、五、奇次谐波分量。周期为T=1的可分解为 方波信号的分解仿真由周期T=1为例：图2为周期为T=1的方波信号，经傅立叶级数分解以后而得到的基波到七次谐波的仿真图，左上角为基波图，它是一个非常正规的正弦波，幅值在1到1.5之间，要高于原方波的幅值。而且它的角频率与原方波信号相同。右上角为三次谐波图，其也是正弦波，明显，其幅值降到了0.5以下，但是三次谐波的频率是基波的1.5倍。其它图形依次为五次谐波，七次谐波。 图2 周期为T=1方波信号的分解图 方波信号的合成仿真图3为方波信号分解以后取有限次谐波的合成波形。左上方图是单独的基波，是正弦波，波身较为平滑，波峰和波谷尖锐。右上方是基波和三次谐波叠加而成的波，大体仍是正弦的形式，但是波身已经比单独的基波较为陡峭，波峰和波谷出现波动，已经趋向方波，有了方波的雏形。以下依次叠加起五次谐波，七次谐波的波形。图3 周期为T=1方波信号的合成图 图4 偶次谐波与奇次谐波的对比由图4可以看出，由于原方波信号经傅立叶级数分解后，偶次谐波不存在，所以在图中只能观察到奇次谐波。 方波信号的频谱图 图5 方波信号的频谱图 图5为周期信号的频谱图，在频谱图中，=1时，信号的幅值在1.2到1.4之间，=2时，信号的幅值为0，=3时，幅值在0.2到0.4之间， =4、5、6、7、时，幅值有起伏，但总体趋势是呈下降趋势。 方波信号的误差分析 方波信号的误差分析表3-1 方波信号前七项合成的误差分析前N之和基波基波+三次谐波基波+三次谐波+五次谐波基波+三次谐波+五次谐波+七次谐波en（t）0.99800.99600.99400.9920图6 方波的误差分析图由图6和表3-1知道，在信号合成时，其叠加的谐波次数越多，将产生的误差值将越小，说明，合成波形越加的向原三角波形靠近。3.2结论1、 由图5可见，周期信号的频谱图有以下特点：2、 （1）离散性。频谱图中的变量为，由于n只能是整数（单边频谱中是正整数），因而谱线是离散的而非连续的，谱线的间隔是，所以周期信号的频谱是离散频谱。（2）谐波性。由于n只取整数，因而谱线在频谱轴上的位置是基频的整数倍。（3）收敛性。幅度谱中各谱线的高度尽管不一定岁随谐波次数的增高作单调的减小，中间有可能有起伏，但总的趋势是随n的增高而减小的，当n为时，高度趋于零。 二、由图2可以得出，任何周期信号都可以由一系列的正弦（或余弦）波组成，随着谐波次数的增大，谐波的幅值越来越小，频率越来越大。三、由图3可以得到，合成波形所包含的谐波分量越多时，除间断点附近外，它越接近与原波形信号。在间断点附近，随着所含有的谐波次数的增加，合成波形的波身越陡峭，波峰越靠近间断点，但尖峰幅度并未明显减小。在傅立叶级数的项数取得很大时，间断点出尖峰下的面积非常小以趋近于零，因而在均方的意义上合成波形同原波形的真值之间没有区别。附录：（程序）方波的分解：t=-3*pi:pi/100000:3*pi;f=square(2*pi*t,50);f1=4*sin(2*t*pi)/pi;f2=4*sin(6*t*pi)/(pi*3);f3=4*sin(10*t*pi)/(pi*5);f4=4*sin(14*t*pi)/(pi*7);subplot(221),plot(t,f1);hold onplot(t,f,r-);grid on;axis(-2 2 -1.5 1.5);subplot(222),plot(t,f2);hold onplot(t,f,r-);grid on;axis(-2 2 -1.5 1.5);subplot(223),plot(t,f3);hold onplot(t,f,r-);grid on;axis(-2 2 -1.5 1.5);subplot(224),plot(t,f4);hold onplot(t,f,r-);grid on;axis(-2 2 -1.5 1.5);方波信号的合成：t=-3*pi:pi/100000:3*pi;f1=4*sin(2*t*pi)/pi;f2=4*sin(2*t*pi)/pi+4*sin(6*t*pi)/(pi*3);f3=4*sin(2*t*pi)/pi+4*sin(6*t*pi)/(pi*3)+4*sin(10*t*pi)/(pi*5);f4=4*sin(2*t*pi)/pi+4*sin(6*t*pi)/(pi*3)+4*sin(10*t*pi)/(pi*5)+4*sin(14*t*pi)/(pi*7);subplot(221),plot(t,f1);hold onplot(t,f,r-);axis(-2 2 -1.5 1.5);grid on;subplot(222),plot(t,f2);hold onplot(t,f,r-);axis(-2 2 -1.5 1.5);grid on;subplot(223),plot(t,f3);hold onplot(t,f,r-);axis(-2 2 -1.5 1.5);grid on;subplot(224),plot(t,f4);hold onplot(t,f,r-);axis(-2 2 -1.5 1.5);grid on;方波信号的奇偶次谐波的对比：t=-3*pi:pi/100000:3*pi;f1=4*sin(2*t*pi)/pi+4*sin(6*t*pi)/(pi*3)+4*sin(10*t*pi)/(pi*5)+4*sin(14*t*pi)/(pi*7);subplot(211),plot(t,f1);axis(-2 2 -1.5 1.5);grid on;t=-3*pi:pi/100:3*pi;f2=0;subplot(212),plot(t,f2);误差分析：e1=f-f1;e2=f-f2;e3=f-f3;e4=f-f4;subplot(221),plot(t,e1);axis(-2 2 -1.5 1.5);grid on;subplot(222),plot(t,e2);axis(-2 2 -1.5 1.5);grid on;subplot(223),plot(t,e3);axis(-2 2 -1.5 1.5);grid on;subplot(224),plot(t,e4);axis(-2 2 -1.5 1.5);grid on;方波信号的频谱图：N=7;n=1:N;for i=1:2:NC(i)=4/(pi*(2*i-1) ;end;stem(n,C) ;4、 总结与体会 首先。这次系统建模能够得以顺利完成，并非我一人之功劳，是指导老师和帮助过我的同学对我的教诲、帮助和鼓励的结果。我要在这里对他们表示深深的谢意！通过这次课程设计，让我对MATLAB这个软件有了很丰富的了解，掌握了MATLAB一些基本的功能和技巧，并且能在设计时充分利用。这对我以后的学习有很大的帮助，尤其作为电子专业的我更加需要掌握这个软件。MATLAB的功能很强大，在这次实习中感受到了它对信号的处理功能极其强大，而且直观的图示和简单方便的操作更让人易懂。但我知道这只是冰山一角，需要学习的东西还很多，所以以后