信号与系统实验-周期信号的合成与分解

1武汉大学教学实验报告电子信息学院专业年月日实验名称周期信号的合成与分解指导教师姓名年级学号成绩一、预习部分1.实验目的2.实验基本原理3.主要仪器设备（含必要的元器件、工具）一、实验目的1．在理论学习的基础上，通过实验深刻领会周期信号傅里叶级数分解的物理意义。2．理解实际应用中通常采用有限项级数来逼近无限项级数，此时方均误差随项数的增加而减小。3．观察并初步了解Gibbs现象。4．深入理解周期信号的频谱特点，比较不同周期信号频谱的差异。二、实验原理满足Dirichlet条件的周期信号f(t)可以分解成三角函数形式的傅里叶级数，表达式为：式中n为正整数；角频率ω1由周期T1决定：112=T。该式表明：任何2满足Dirichlet条件的周期信号都可以分解成直流分量及许多正弦、余弦分量。这些正弦、余弦分量的频率必定是基频2022141()sin()cos(2)22nEEnftnftn的整数倍。通常把频率为1f的分量称为基波，频率为n1f的分量成为n次谐波。周期信号的频谱只会出现在0,1，21，…，n1，…等离散的频率点上，这种频谱称为离散谱，是周期信号频谱的主要特点。f(t)波形变化越剧烈，所包含的高频分量的比重就越大；变化越平缓，所包含的低频分量的比重就越大。一般来说，将周期信号分解得到的三角函数形式的傅里叶级数的项数是无限的。也就是说，通常只有无穷项的傅里叶级数才能与原函数精确相等。但在实际应用中，显然无法取至无穷多项，而只能采用有限项级数来逼近无穷项级数。而且，所取项数越多，有限项级数就越逼近原函数，原函数与有限项级数间的方均误差就越小，而且低次谐波分量的系数不会因为所取项数的增加而变化。当选取的傅里叶有限级数的项数越多，所合成的波形的峰起就越靠近f(t)的不连续点。当所取得项数N很大时，该峰起值趋于一个常数，约等于总跳变值的9%，这种现象称为Gibbs现象。三、需要掌握的MATLAB函数结果的显示会用到plot函数，请参考MATLAB帮助。3二、实验操作部分1.实验数据、表格及数据处理2.实验操作过程（可用图表示）3.实验结论4四、实验内容1．周期对称方波信号的合成图示方波既是一个奇对称信号，又是一个奇谐信号。根据函数的对称性与傅里叶系数的关系可知，它可以用无穷个奇次谐波分量的傅里叶级数来表示:选取奇对称周期方波的周期T=0.02s，幅度E=6，请采用有限项级数替代无限项级数来逼近该函数。分别取前1、2、5和100项有限级数来近似，编写程序并把结果显示在一幅图中，观察它们逼近方波的过程。MATLAB程序如下：%奇对称方波合成f(t)=2*E/pi*(sum(sin(2*pi*f*t*(2*i+1)/(2*i+1));i=0,1,2¡­¡­)t=0:0.001:0.1;E=6;A=2*E/pi;T=0.02;f=1/T;y=A*sin(2*pi*f*t);subplot(221)plot(t,y);5axis([0,0.1,-4,4]);xlabel('time');ylabel('前1项有限级数');y=A*(sin(2*pi*f*t)+sin(2*pi*f*t*3)/3);subplot(222);plot(t,y);axis([0,0.1,-4,4]);xlabel('time');ylabel('前2项有限级数');y=A*(sin(2*pi*f*t)+sin(2*pi*f*t*3)/3+sin(2*pi*f*t*5)/5+sin(2*pi*f*t*7)/7+sin(2*pi*f*t*9)/9);subplot(223)plot(t,y);axis([0,0.1,-4,4]);xlabel('time');ylabel('前5项有限级数');t=0:0.001:0.1;y=0;N=10;fori=1:N%设置N的值来选项数；y=y+A*(sin((2*i-1)*2*pi*f*t)/(2*i-1));endsubplot(224);plot(t,y);axis([0,0.1,-4,4]);xlabel('time');ylabel('前10项有限级数');6显示结果如图4-2所示。图4-2奇对称方波信号的合成2．观察Gibbs现象分别取前5,6,7和8项有限级数来逼近奇对称方波，观察Gibbs现象。MATLAB程序如下：%观察Gibbs现象%奇对称方波合成，观察Gibbs现象t=0:0.001:0.04;sishu=12/pi;y=0;fori=1:5y=y+sishu*(sin((2*i-1)*100*pi*t)/(2*i-1));endsubplot(221)plot(t,y);axis([0,0.04,-4,4]);xlabel('time');ylabel('前5项有限级数);y=0;7fori=1:6y=y+sishu*(sin((2*i-1)*100*pi*t)/(2*i-1));endsubplot(222);plot(t,y);axis([0,0.04,-4,4]);xlabel('time');ylabel('前6项有限级数');y=0;fori=1:7y=y+sishu*(sin((2*i-1)*100*pi*t)/(2*i-1));endsubplot(223)plot(t,y);axis([0,0.04,-4,4]);xlabel('time');ylabel('前7项有限级数');y=0;fori=1:8y=y+sishu*(sin((2*i-1)*2*pi*f*t)/(2*i-1));endsubplot(224);plot(t,y);axis([0,0.04,-4,4]);xlabel('time');ylabel('前8项有限级数');显示结果如图4-3所示。8图4-3Gibbs现象3．周期对称三角信号的合成设计采用有限项级数逼近偶对称周期三角信号的实验，编制程序并显示结果。4．周期信号的频谱分析奇对称方波信号与偶对称三角信号的频谱，编制程序并显示结果，深入讨论周期信号的频谱特点和两信号频谱的差异。五、实验要求1.输入实验内容1中提供的奇对称方波信号合成的MATLAB程序，生成M文件，编译并运行，观察合成结果。2.输入实验内容2中提供的有限项级数逼近方波信号的MATLAB程序，生成M文件，编译并运行，观察Gibbs现象。3.自行编制完整的MATLAB程序，完成实验内容3中偶对称三角信号的合成。在实验报告中给出程序和显示结果。9该信号的傅里叶级数表示为:2022141()sin()cos(2)22nEEnftnftn选取偶对称周期三角信号T=0.02s，幅度E=6，采用有限项级数替代无限项级数来逼近该函数。分别取前1、2、5和100项有限级数来近似。MATLAB程序如下：%偶对称三角波合成f(t)=E/2+4*E/((2*i-1)*pi)^2*(cos((2*i-1)*2*pi*f*t))E=6;T=0.02;f=1/T;t=-0.1:0.001:0.1;y=E/2;fori=1:1y=y+4*E/((2*i-1)*pi)^2*(cos((2*i-1)*2*pi*f*t));endsubplot(221)plot(t,y);axis([-0.10.106.5]);gridon;xlabel('time');ylabel('前1项有限级数');y=E/2;fori=1:2y=y+4*E/((2*i-1)*pi)^2*(cos((2*i-1)*2*pi*f*t));endsubplot(222)plot(t,y);axis([-0.10.106.5]);gridon;xlabel('time');ylabel('前2项有限级数');y=E/2;fori=1:5y=y+4*E/((2*i-1)*pi)^2*(cos((2*i-1)*2*pi*f*t));endsubplot(223)plot(t,y);axis([-0.10.106.5]);gridon;10xlabel('time');ylabel('前5项有限级数');y=E/2;fori=1:100y=y+4*E/((2*i-1)*pi)^2*(cos((2*i-1)*2*pi*f*t));endsubplot(224)plot(t,y);axis([-0.10.106.5]);gridon;xlabel('time');ylabel('前100项有限级数');显示结果如图4-4所示。11图4-4偶对称三角波信号的合成4.自行编制完整的MATLAB程序，完成实验内容4中奇对称方波信号和偶对称三角波信号的频谱分析。在实验报告中给出程序和显示结果，讨论周期信号的频谱特点和两信号频谱的差异。MATLAB程序如下：%周期性方波的傅立叶级数展开幅度频谱m=1:2:25;E=6;an=2*E./m/pi;subplot(211);stem(m,an,'fill');holdon;m0=0;a0=0;stem(m0,a0,'fill');gridon;xlabel('基频倍数');ylabel('周期性方波幅度频谱');holdoff;%周期性三角波的傅立叶级数展开幅度频谱m=1:2:25;E=6;an=4*E./(pi.*m).^2;subplot(212);stem(m,an,'fill');holdon;m0=0;y0=E/2;stem(m0,y0,'fill');gridon;xlabel('基频倍数');ylabel('周期性三角波幅度频谱');holdoff;显示结果如图4-5所示。12周期信号的频谱具有如下特点：(1)离散性。周期信号的频谱是由不连续的谱线组成，每条谱线代表一个谐波分量。(2)谐波性。频谱中每条谱线只出现在基波频率的整数倍上，由于这里的方波和三角波都是对称的，频谱中只包含基频的奇数倍频率。(3)收敛性。各频率分量的谱线高度表示各次谐波分量的幅值。两信号频谱的差异:由以上周期性方波和三角波信号的频谱分析可知，周期性三角波信号的各次谐波幅值衰减比周期性方波的频谱衰减快得多，这说明三角波的频率结构中低频成分较多，而方波的高频成分比较多，同时偶对称的三角波中含有直流分量而奇对称的方波中没有直流分量。六、思考题1.利用有限项的指数形式的傅里叶级数重复奇对称方波信号的合成。答：其指数形式的傅里叶级数的表示为：程序如下：13t=0:0.001:1;T=0.2;f=1/T;w=2*pi*f;E=6/pi;y=0;forn=1:100y=y+E*(exp(1i*(2*n-1)*w*t-1i*0.5*pi))/(2*n-1);endplot(t,y);gridon;axis([0,1,-4,4]);xlabel('time');2.分析时域信号的间断性与其频谱谐波收敛速率的对应关系14若时域信号间断点较多，则说明其高频分量较多，则谐波收敛速度会变慢。也可以从冲击信号的角度来分享这个问题，当信号出现间断时，例如方波的陡降，在陡降处相当于叠加了一个冲击信号，而冲击信号的频谱是在负无穷到正无穷上均匀分布的，这导致了最终频谱的分布变得更加均匀，高频分量更多，谐波收敛更慢。15三、实验效果分析（包括仪器设备等使用效果）实验分析：1、图形曲线不连续是因为matlab中作图时是取的有限的点，无法做到连续连线，故画出的图形曲线会出现间断或转折等情况。2、作的图形不是完全标准的方波或三角波是因为我们是用有限项傅里叶级数去逼近的，无法到达用无穷项去逼近作图的效果，如果绘图时选取的点数目过少，也会导致图像的转折处出现缓降而变得倾斜。四、教师评语指导教师年月日