实验3连续信号的频域分析.docVIP

实验三 连续信号的频域分析 1.实验目的： （1）掌握周期信号分解和合成的方法。 （2）掌握非周期信号频谱分析方法。 2.实验原理 （1）周期信号的分析与合成 周期为T的周期信号的傅里叶级数有三角型与指数型两种形式，分别表示为： 和 式子中，各系数的计算公式和相互关系如下： ， 傅里叶级数表明周期信号可以分解为正弦信号或虚指数的线性组合。由三角型的系数可画出周期信号的单边幅度谱和相位谱，由指数型的系数可画出周期信号的双边谱，它们都是离散谱。 上述系数可以用MATLAB函数quad或quadl计算。它们的调用格式为： y = quad(FUN,A,B), 和 y = quadl(FUN,A,B) 其中，FUN是被积函数名或函数句柄；A和B分别是积分区间的下限和上限。比如，系数可以如下计算： ，(T为信号周期) 例1 求周期为4，幅度为1、脉冲宽度为2的对称矩形脉冲信号的三角型傅里叶级数。 参考： T = 4; A = 1; tao = 2; w0 = 2*pi/T; N = 6; f = @(n) (2*quad (@(t)(A*rectpuls(t,tao).*cos(n*w0*t)),-T/2,T/2)/T); for k = 1:N an(k) = f(k); End fun = @(t)(A*rectpuls(t,tao)); a0 = quad (fun,-T/2,T/2)/T; an = [a0,an]; 例2 画出例1中前6次谐波合成的信号波形，要求画2个周期。 t = -T:0.01:T; N=6; n=(0:N).’; ft = an*cos(w0*n*t); plot(t,ft); （2）非周期信号的频谱分析 非周期的频谱是连续谱，MATLAB通常用数值方法求其频谱，一种方法是利用DFT变换，另一种方法是利用函数quad和quadl。调用格式参见上文。 例3求门宽为2的门函数的傅里叶变换，画出其幅度谱。 参考： w1 = -2*pi;w2 = -w1; t1 = -tao/2;t2 = -t1; N = 500; wk = linspace(w1,w2,N); F = zeros(1,N); Fw = @(w)(quad(@(t)(rectpuls(t,tao).*exp(-j*w*t)),t1,t2)); for k = 1:N F(k) = Fw(wk(k)); End Plot(wk,abs(F)); 3.实验内容及其步骤 本次实验需要完成以下内容： 复习有关傅里叶级数和傅里叶变换的有关内容。 求周期为5，幅度为2、脉冲宽度为2的对称矩形脉冲信号的三角型傅里叶级数， 画出其幅度谱（取前7次谐波）。2）将前7次谐波叠加，画出合成的波形（画至少2个周期）。 提示：参考例1，例2. 求单边指数信号的傅里叶变换，画出其在区间[-,]之间的幅度谱。 4. 实验用MATLAB函数介绍 在实验过程中用到的MATLAB指令（函数）有： rectpuls, cos,quad,quadl, plot,linspace,zeros,abs,exp. 5.思考练习题 （1）分析、对比实验内容1中合成波形和原周期信号波形之间的不同，简要说明原因，增多谐波个数，合成波形有什么变化？ （2）对实验过程中所涉及的问题进行分析，试编写和修改相应的程序，得出最终正确的结果和波形图，并对实验报告进行整理分析。 6.实验报告要求： （1）明确实验目的以及实验的原理。 （2）通过实验总结利用MATLAB分析信号频谱的一般方法。 （3）完成思考题的内容，对实验结果及其波形图进行分析对比，总结主要结论。