西南石油大学实验报告
一实验目的:
学习利用MATLAB进行基于窗函数的FIR数字滤波器的设计。
二实验内容:
利用矩形窗、哈明窗、汉宁窗和布莱克曼窗设计一个FIR低通滤波器,已知ωc=0.25π,N=10
三实验步骤:
1、实验程序
N=10;
M=128;
b1=fir1(N,0.25,boxcar(N+1)); %用矩形窗作为冲激响应的窗函数
b2=fir1(N,0.25,hamming(N+1)); %用哈明窗作为冲激响应的窗函数
b3=fir1(N,0.25,blackman(N+1)); %用布莱克曼窗作为冲激响应的窗函数
b4=fir1(N,0.25,hanning(N+1)); %用汉宁窗作为冲激响应的窗函数
h1=freqz(b1,1,M); %矩形窗对应的频率响应
h2=freqz(b2,1,M); %哈明窗对应的频率响应
h3=freqz(b3,1,M); %布莱克曼窗对应的频率响应
h4=freqz(b4,1,M); %汉宁窗对应的频率响应
f=0:0.5/M:0.5-0.5/M;
plot(f,abs(h1),'-.',f,abs(h2),f,abs(h3),'*',f,abs(h4),':');
legend('¾ØÐδ°','¹þÃ÷´°','BLACKMAN','ººÄþ´°');
grid;
ylabel('magnitude response');
xlabel('w/(2*pi)');
axis([0 0.5 0 1.2]);
set(gca,'XTickMode','manual','XTick',[0,0.25,0.5]);
set(gca,'YTickMode','manual','XTick',[0,0.5,1]);
2、运行结果如下图所示
四、实验结论
在实验后,我对MATLAB软件有了进一步的了解,也在不断的实践中,更多的熟悉了MATLAB的编程,在编程方面一点点的有了进步。
只是,还需要课后多加的去研究MATLAB。
此次实验收获很大。