数字信号处理实验报告吉布斯效应.doc

《数字信号处理》 实验报告 实 验 名 称： 吉布斯效应 专 业 班 级 ： XXXXX 学 生 姓 名 ： XXXXX 学 号 ： XXXXXXXX 指 导 教 师 ： XXXXXX 实 验 时 间： 2011年6月15日Sc（n）=1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n)-2/7/pi*cos(7*pi*n)… 实验结果 （1）由该程序得到的波形图 DeT = 1/128*4; n = -1:DeT:1; nTao = -1/2:DeT:1/2; sTao(1:length(nTao)) = ones(1, length(nTao)); s(1:length(n)) = [zeros(1, 1/2*length(nTao)) sTao zeros(1, 1/2*length(nTao))]; plot(n, s) （2）下图为N=1，N=3，N=5，N=7的图形。 （3）下图为N=9，N=27的图形。 （4）该图为N=1000的波形图。 5 结论 当n值越大时，越接近理论值。 6 心得体会: 本次课程设计中，通过运用matlab应用程序对吉布斯效应实验过程有了进一步的了解，通过改变不同的参数值，由所对应的波形可反映出吉布斯效应。 7 实验代码： %实验设计程序如下： function WavConstruction clear clc close all DeT = 1/128*4; n = -1:DeT:1; nTao = -1/2:DeT:1/2; sTao(1:length(nTao)) = ones(1, length(nTao)); s(1:length(n)) = [zeros(1, 1/2*length(nTao)) sTao zeros(1, 1/2*length(nTao))]; plot(n, s) xlim([-1.05 1.05]) ylim([-0.2 1.21]) figure sc1 = 1/2 + 2/pi*cos(pi*n); subplot(2,2,1) plot(n, sc1) title(N = 1) subplot(2,2,2) sc2 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n); plot(n, sc2) title(N = 3) subplot(2,2,3) sc3 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n); plot(n, sc3) title(N = 5) subplot(2,2,4) sc4 =1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n) - 2/pi/7*cos(7*pi*n); plot(n, sc4) title(N = 7) figure subplot(2,1,1) sc5 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n) - 2/pi/7*cos(7*pi*n) + 2/pi/9*cos(9*pi*n); plot(n, sc5) title(N = 9) sc6 = sc5-2/pi/11*cos(11*pi*n)+2/pi/13*cos(13*pi*n)-2/pi/15*cos(15*pi*n)+2/pi/17*cos(17*pi*n)-2/pi/19*cos(19*pi*n) ... +2/pi/21*cos(21*pi*n)-2/pi/23*cos(23*pi*n)+2/pi/25*cos(25*pi*n)-2/pi/27*cos(27*pi*n); subplot(2,1,2) plot(n, sc6) title(N = 27) figure N_Order = 1000; scN_Order = 1/2; for i_Index = 1:1:N_Order scN_Order = scN_Order - (-1)^i_Index*2/pi/(2*i_Index-1)*cos((2*i_Index-1)*pi*n); end plot(n, scN_Order) title([N = ,int2str(N_Order)]) grid