文档介绍:实验三周期信号的频谱分析一、实验目的 1 、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法; 2、观察截短傅里叶级数而产生的“ Gibbs 现象”, 了解其特点以及产生的原因; 3 、掌握各种典型的连续时间非周期信号的频谱特征。二、原理说明: 1 、连续时间周期信号的傅里叶级数分析任何一个周期为 T 1 的正弦周期信号, 只要满足狄利克利条件, 就可以展开成傅里叶级数。其中三角傅里叶级数为: ?????? 1 000 )] sin( ) cos( [)( k kktkbtkaatx?? 或:?????? 1 00) cos( )( k kktkcatx?? 其中 1 02T ???,称为信号的基本频率( Fundamental frequency ),kkbaa,和, 0 分别是信号)(tx 的直流分量、余弦分量幅度和正弦分量幅度, kkc?、为合并同频率项之后各正弦谐波分量的幅度和初相位, 它们都是频率 0?k 的函数, 绘制出它们与0?k 之间的图像, 称为信号的频谱图( 简称“频谱”),kc -0?k 图像为幅度谱, k?-0?k 图像为相位谱。三角形式傅里叶级数表明, 如果一个周期信号 x(t) , 满足狄里克利条件, 那么, 它就可以被看作是由很多不同频率的互为谐波关系( harmonically related ) 的正弦信号所组成,其中每一个不同频率的正弦信号称为正弦谐波分量(ponent) ,其幅度( amplitude )为 kc 。也可以反过来理解三角傅里叶级数:用无限多个正弦谐波分量可以合成一个任意的非正弦周期信号。指数形式的傅里叶级数为:?????? k t jkkeatx)( ? 其中, ka 为指数形式的傅里叶级数的系数,按如下公式计算: ???? 2/2/ 1)( 1 TT t jk kdtetxT a ? 指数形式的傅里叶级数告诉我们, 如果一个周期信号 x(t) , 满足狄里克利条件, 那么, 它就可以被看作是由很多不同频率的互为谐波关系( harmonically related ) 的周期复指数信号所组成, 其中每一个不同频率的周期复指数信号称为基本频率分量, 其复幅度( complex amplitude )为ka 。这里“复幅度( complex amplitude )”指的是 ka 通常是复数。上面的傅里叶级数的合成式说明, 我们可以用无穷多个不同频率的周期复指数信号来合成任意一个周期信号。然而, 用计算机( 或任何其它设备) 合成一个周期信号, 显然不可能做到用无限多个谐波来合成, 只能取这些有限个谐波分量来近似合成。假设谐波项数为 N ,则上面的和成式为: ???? NNk t jkkeatx)( ? 显然, N 越大,所选项数越多,有限项级数合成的结果越逼近原信号 x(t) 。本实验可以比较直观地了解傅里叶级数的物理意义, 并观察到级数中各频率分量对波形的影响包括“ Gibbs ”现象: 即信号在不连续点附近存在一个幅度大约为 9% 的过冲, 且所选谐波次数越多, 过冲点越向不连续点靠近。这一现象在观察周期矩形波信号和周期锯齿波信号时可以看得很清楚。 2 、傅里叶级数的 MATLAB 计算设周期信号 x(t) 的基本周期为 T 1 ,且满足狄里克利条件,则其傅里叶级数的系数可由式 计算得到。式 重写如下: ???? 2/2/ 1)( 1 TT t jk kdtetxT a ?基本频率为: 1 02T ???对周期信号进行分析时, 我们往往只需对其在一个周期内进行分析即可, 通常选择主周期( Principle period ) 。假定 x 1 (t) 是 x(t) 中的主周期,则???? 2/2/ 11)( 1 TT t jk kdtetxT a ?计算机不能计算无穷多个系数,所以我们假设需要计算的谐波次数为 N ,则总的系数个数为 2N+1 个。在确定了时间范围和时间变化的步长即 T 1和 dt 之后, 对某一个系数,上述系数的积分公式可以近似为: ??????? n t jkn TT t jk kTdtetxdtetxT a 1 2/2/ 11/)()( 1 ??1 21/],,[ )]( ),( ),([Tdteeetxtxtx t jkt jkt jkM???????????对于全部需要的 2N+1 个系数, 上面的计算可以按照矩阵运算实现。 MATLAB 实现系数计算的程序如下: dt= ; T= 2;t= -T/2:dt:T/2; w0= 2*pi/T; x1= input( ‘ Type in the periodic signal x(t) over one period x1(t)= ’);