武汉理工大学《MATLAB原理与应用》课程设计
图8fdatool工作栏
7.2 FIR滤波器参数设置
fdatool界面总共分两大部分，一部分是design filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。design filter部分主要分为： filter type（滤波器类型）选项，包括lowpass（低通）、highpass（高通）、bandpass（带通）、bandstop（带阻）和特殊的fir滤波器。 design method（设计方法）选项，包括iir滤波器的butterworth（巴特沃思）法、chebyshev type i（切比雪夫i型）法、 chebyshev type ii（切比雪夫ii型） 法、elliptic（椭圆滤波器）法和fir滤波器的equiripple法、least-squares（最小乘方）法、window（窗函数）法。 filter order（滤波器阶数）选项，定义滤波器的阶数，包括specify order（指定阶数）和minimum order（最小阶数）。在specify order中填入所要设计的滤波器的阶数（n阶滤波器，specify order＝n-1），如果选择minimum order则matlab根据所选择的滤波器类型自动使用最小阶数。 frenquency specifications选项，可以详细定义频带的各参数，包括采样频率fs和频带的截止频率。它的具体选项由filter type选项和design method选项决定 在Matlab 命令窗口输入FDATool 命令,按确定调出FDATool 界面。在
17
武汉理工大学《MATLAB原理与应用》课程设计
ResponseType 下选择滤波器的类型为Highpass。在Design Method 下选择设计方法为FIRequiripple,在Filter Order选择 minimum order，在frenquency specifications里输入wstop 0.5和 wpass 0.6，在Magnituide Specifications 里输入astop 50,执行后的图如图所示。
图9 FDATOOL工具箱设计高通FIR滤波器
18
武汉理工大学《MATLAB原理与应用》课程设计
8心得体会
这学期我们学习了《数字信号处理》，这是我们的专业课，是理论方面的指示。正所谓“纸上谈兵终觉浅，觉知此事要躬行”。学习任何知识，仅从理论上去求知，而不去实践、探索是不够的，所以在暨《数字信号处理》之后紧接着来一次MATLAB课程设计是很及时、很必要的。这样不仅能加深我们对数字信号处理的认知，而且还及时、真正的做到了学以致用。
在此过程中，我们通过查找大量资料，请教老师，以及不懈的努力，不仅培养了独立思考、动手操作的能力，在各种其它能力上也都有了提高。更重要的是，在课设中，我们学会了很多学习的方法，而这也是日后最实用的。
不管怎样，这些都是一种锻炼，一种知识的积累，能力的提高。完全可以把这个当作基础东西，只有掌握了这些最基础的，才可以更进一步，取得更好的成绩。很少有人会一步登天，永不言弃才是最重要的。而且，这对于我们的将来也有很大的帮助。
19
武汉理工大学《MATLAB原理与应用》课程设计
参考文献
1.《数字信号处理及其MATLAB实现》. Vinay K.Ingle主编.电子工业出版社.1998年出版
2.《Digital Signal Processing Laboratory Using Matlab》. Sanjit K. Miltra编著. McGraw-Hill出版社.2000年出版
3.《数字信号处理原理与实现》第二版.刘泉主编.电子工业出版社
4.《数字信号处理教程—MATLAB释义与实现》.陈怀琛主编.电子工业出版社.2004年版
5.《信号与系统》.刘泉主编.高等教育出版社.2006年版 6.《MATLAB程序设计》.阮沈青主编.电子工业出版社.2004年版
20
武汉理工大学《MATLAB原理与应用》课程设计
附件：MATLAB程序
N=33;
wc=3*pi/5;T1=0.2;T2=0.7; N1=fix(wc/(2*pi/N)); N2=N-2*N1;
A=[zeros(1,N1),T1,T2,ones(1,N2-4),T2,T1,zeros(1,N1)]; theta=-pi*[1:N]*(N-1)/N; H=A.*exp(j*theta); h=real(ifft(H)); v=1:N; figure(1); plot(v,A,'*'); title('频率采样'); ylabel('H(k)');
axis([0,fix(N*1.1),-0.1,1.1]); figure(2);
stem(v ,h);title('脉冲响应'); ylabel('h(n)'); axis([0,fix(N*1.1),min(h)*1.1,max(h)*1.1]); M=5000; nx=[1:N];
w=linspace(0,pi,M); X=h*exp(-j*nx'*w); figure(3);
plot(w./pi,abs(X),'k');xlabel('\\omega/\\pi');ylabel('|Hd(w)|'); title('幅度响应');axis([0,1,-0.1,1.3]); figure(4);
plot(w./pi,20*log10(abs(X)),'k');title('幅度响应'); xlabel('\\omega/\\pi');ylabel('dB');axis([0,1,-80,10]);
21