數(shù)字信號處理實驗報告--吉布斯效應(yīng)

1、數(shù)字信號處理實驗報告實 驗 名 稱：吉布斯效應(yīng)專 業(yè) 班 級 ： xxxxx 學(xué) 生 姓 名 ： xxxxx 學(xué)號 ： xxxxxxxx 指 導(dǎo) 教 師 ： xxxxxx 實 驗 時 間： 2011年 6 月 15 日1 實驗?zāi)康模候炞C吉布斯效應(yīng)2 運(yùn)行環(huán)境：pc機(jī)，windows2000及其以上， matlab 軟件。3 實驗設(shè)計：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)4實驗結(jié)果（1）由該程序得到的波形圖det = 1/128*4; n = -1:det

2、: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 結(jié)論當(dāng) n 值越大時，越接近理論值。6 心得體會 : 本次課程設(shè)計中，通過運(yùn)用matlab 應(yīng)用程序?qū)妓剐?yīng)實驗過程有了進(jìn)一步的了解，通過改變不同的參

3、數(shù)值，由所對應(yīng)的波形可反映出吉布斯效應(yīng)。7 實驗代碼：% 實驗設(shè)計程序如下：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) figur

4、e 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/

5、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/1

6、3*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 on 參考文獻(xiàn)參考文獻(xiàn)格式如下： (1 書籍2 文章例