山东大学信息学院信号和系统第五次实验报告5.docVIP

山东大学·信号与系统·实验报告 PAGE \* MERGEFORMAT PAGE \* MERGEFORMAT - 2 - 信息科学与工程学院 2017－2018学年第二学期 实 验 报 告 课程名称： 信号与系统实验 实验名称： 周期信号的傅里叶级数与频谱 实验报告 【实验源程序及图形】 一：画出频谱图 n=-20:20; tau=0.00004;T=0.0002;w1=2*pi/T; x=n*tau/T; fn=5*tau/T*sinc(x); plot(311); stem(n*w1,fn),grid on title(\tau=40us,T=0.2ms) 二：观测吉布斯现象 画出原始方波 t=-0.0003:0.0000001:0.0003;ft=2.5+square(10000*pi*(t+0.00002),20)*5/2; plot(t,ft),grid on axis([-0.0003,0.0003,-1,6]) t=-0.0003:0.0000001:0.0003; n=1:1:1; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) title([最大谐波次数=1]) n=1:1:3; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) t=-0.0003:0.0000001:0.0003; n=1:1:7; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) title([最大谐波次数=7]) n=1:1:11; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) t=-0.0003:0.0000001:0.0003; n=1:1:17; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) title([最大谐波次数=17]) n=1:1:57; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) t=-0.0003:0.0000001:0.0003; n=1:1:79; b=10*sin(n*pi/5)./(n*pi); x=b*cos(10000*pi*n*t)+1; plot(t,x),grid on axis([-0.0003,0.0003,-1,6]) title([最大谐波次数=79]) 【实验分析与讨论】 在实验前，需要先计算傅里叶展开式，确定对应系数，以及展开形式， 作出图像后，用计算结果与之比对，看是否符合 观察到当最大谐波次数足够大时，图像非常逼近方波，只是在间断点处存在过冲，这种过冲大小大约为跳变值的9%，即为吉布斯现象