数字信号处理报告(3)

可知： 同理：

带通滤波器的单位冲击响应为： h(n)=(Sinω2n- Sinω1n)/ (πn) 所以： y(n)?

m?0?h(m)x(n?m)

N?1四、实验结果及分析：

以下各图输入序列均为128个点，输出序列均为256个点，滤波器窗函数取样点的数目N均为61。 低通滤波

（1） 单位冲击序列输入：

输入序列1x(n)0.500204080100n加矩形窗低通滤波后频谱601201401幅度/dB0.50050150200n加矩形窗低通滤波后输出序列1002503000.2y(n)0-0.2050100150n200250300

输入序列1x(n)0.500204080100n加海明窗低通滤波后频谱601201400.5幅度/dB0050150200n加海明窗低通滤波后输出序列1002503000.1y(n)0-0.1050100150n200250300

加矩形窗后响应0幅度/dB-10-20-3000.10.20.30.40.50.6归一化频率/p加海明窗后响应0.70.80.91-5-10幅度/dB-15-20-2500.10.20.30.40.50.6归一化频率/p0.70.80.91

主要代码：

n=0:127;N=30;x(1)=1;x(2:128)=0;

h=sin(0.2*pi*(n-30))./(pi*(n-30+eps));y=conv(x,h);Y=abs(fft(y)); subplot(3,1,1); stem(n,x); xlabel('n'); ylabel('x(n)'); title('输入序列');

subplot(3,1,2); stem(0:254,Y) ;xlabel('n'); ylabel('幅度/dB'); title('加矩形窗低通滤波后频谱'); subplot(3,1,3); stem(0:254,y); xlabel('n'); ylabel('y(n)'); title('加矩形窗低通滤波后输出序列');

hd=sin(0.2*pi*(n-30))./(pi*(n-30+eps));

w=0.5-0.5*cos(2*pi*n/128); h=hd.*w; y=conv(x,h); Y=abs(fft(y)); subplot(3,1,1) ;stem(n,x); xlabel('n'); ylabel('x(n)'); title('输入序列');

subplot(3,1,2); stem(0:254,Y); xlabel('n'); ylabel('幅度/dB'); title('加海明窗低通滤波后频谱'); subplot(3,1,3); stem(0:254,y); xlabel('n'); ylabel('y(n)'); title('加海明窗低通滤波后输出序列'); n=0:127;h=sin(0.2*pi*(n-30))./(pi*(n-30+eps));[h1,w1]=freqz(h,1); subplot(2,1,1); plot(w1/pi,20*log10(abs(h1)));