1、健雄职业技术学院毕业设计（论文）1目录1.绪论 .31.1 数字滤波器定义 31.2 数字滤波器的分类 31.3 数字滤波器的技术指标 32.FIR数字滤波器的条件和特点 52.1 线性相位 FIR 数字滤波器 .52.2 第一类线性相位对 H(N)的约束条件 62.3 线性相位 FIR 滤波器的幅度特性 HG(W)的特点 .72.4 线性相位 FIR 滤波器的零极点分布特点 83.FIR数字滤波器的 MATLAB设计 .93.1 通过窗函数法设计 FIR 滤波器 .93.1.1 窗函数设计原理 93.1.2 典型窗函数介绍 .113.1.3 用窗函数法设计 FIR 录波器的步骤 133.1.
2、4 通过 M 文件实 FIR 数字语音低通滤波器的设计 .133.2 通过等波纹最佳逼近法设计 FIR 数字滤波器 .173.2.1 remez 和 remezord 函数滤波器指标 .173.2.2 通 remez 和 remezord 函数实现 FIR 数字语音低通滤波器的设计 .183.3 通过 FDATOOL实现 FIR 数字滤波器的设计及分析 19雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计24.利用 SIMULINK对 FIR滤波器进行仿真 .234.1SIMULINK模块库简介 .234.2 搭建语音信号 FIR 低通滤波仿真模型 .23总结 2
3、7致谢 28健雄职业技术学院毕业设计（论文）31.绪论当前，通信技术高速发展，业务范围不断扩大，人们对产品的需求迅速增长，滤波器在这些产品电路中就扮演着重要的角色，当然数字滤波器更影响着人们生活的方方面面。数字滤波器既可以是用数字硬件装配成德一台完成给定运算的专用数字计算机，也可将所需的运算编成程序，让通用计算机来执行。MATlAB 是“矩阵实验室” （MATrix LABoratoy）的缩写，它是由美国 Mathworks 公司于 1984 年正式推出的，是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。MATLAB 是功能强大的科学及工程计算软件，它不但具有以矩阵
4、计算为基础的强大数学计算和分析功能，而且还具有丰富的可视化图形表现功能和方便的程序计算能力。MATLAB 的应用领域极为广泛，除数学计算和分析外，还被广泛地应用于自动控制、系统仿真、数字信号处理、图形图像分析、数理统计、人工智能、虚拟现实技术、通信工程、金融系统等领域，因此，MATLAB 是面向 21 世纪的计算机程序设计及科学计算语言。1.1数字滤波器定义所谓数字滤波器，是指输入、输出均为数字信号，通过数字运算处理改变入信号所含频率成分的相对比例，或者滤除某些频率成分的数字器件或程序。因此数字滤波的概念和模拟滤波相同，只是信号的形式和实现滤波的方法不同。正因为数字滤波器通过数值运算实现滤波，
5、所以数字滤波器的处理精度高、稳定、体积小、重量轻、灵活、不存在阻抗匹配问题，可以实现模拟滤波器不能实现的特殊滤波功能。1.2数字滤波器的分类数字滤波器按照不同的分类方法，有许多种类，总体上可分为两大类。一类是经典滤波器，即一般的线性系统滤波器。另一类即所谓的现代滤波器，它的理论建立在随机信号处理的理论基础上，利用了随机信号内部的统计特性对信号进行滤波。从功能上可以分为低通、高通、带通和带阻滤波器，根据数字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计
6、41.3数字滤波器的技术指标1)通带边界频率 pw2)阻带截止频率 s3)通带内允许最大衰减 pmax()20lgjwpHedB0pw4)阻带内允许最小衰减 sax()20lgjwp edH通 带 中阻 带 中健雄职业技术学院毕业设计（论文）52.FIR数字滤波器的条件和特点与 IIR 滤波器相比，FIR 滤波器在保证幅度特性满足技术要求的同时，很容易做到具有严格的线性相位特性。本章用 N 表示 FIR 滤波器单位脉冲响应 h(n)的长度，其系统函数H(z)为：10()()nnHzhzH(z)是 的 N-1 次多项式，所以它在 z 平面内有 N-1 个零点，在原点 z=0 处有 N-1 个重极
7、1z点，因此 H(z)永远稳定。稳定性和相位线性是 FIR 滤波器最为特出的优点。2.1 线性相位 FIR数字滤波器对于长度为 N 的 h(n)，频率响应函数为：（2-1）10()()NjwjwnnHehe（2-2）()j jg式中 为幅度特性， 为相位特性，值得注意的是， 不同于 ，()gw()w()gHw()jwe为实函数，其值可正可负， 总是正。线性相位 FIR 滤波器是指 是 w 的H()jwHe 线性函数，即为常数 (2-3)()如果 满足下式 ()w雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计6是起始相位 (2-4)0()w0严格的说此时 不具有线
8、性相位特性，但以上两种情况都满足群延时是一常数，即： ()dw也称这种情况为线性相位。一般称满足(2-3)式是第一类线性相位，满足(2-4) 式是第二类线性相位。本论文的 FIR 滤波器设计仅采用第一类线性相位。2.2第一类线性相位对 h(n)的约束条件线性相位 FIR 滤波器的时域约束条件是指满足相位线性时，对 h(n)的约束条件。第一类线性相位的 FIR 滤波器的相位函数为 ，联立(2-1) 式和 (2-2)式可得到：()w10()()NjwjwnjgnHeheHe（2-5）10()cosi)()cosin)Ngnhj j 由式（2.1.5）可得： 10()cs()csNgnwH（2-6）
9、10()i()ingnh将（7.2.2）中两式相除可得到： 10()cosiniNnNwh即： 1100()cosin()sincoNNnww移相得： 10()csisic)Nnh健雄职业技术学院毕业设计（论文）7（2-7） 10()sin()0Nnhw函数 h(n)sinw(n- )若能关于(N-1)/2 奇对称，则为（2-7）的一组解，显然 sinw(n- )关于 n= 奇对称，若取 =(N-1)/2，则可取 h(n)关于 n=(N-1)/2 成偶对称，即满足如下条件：()w12N（2-8）1)hnNn0n由以上推论可知：如果要求单位脉冲响应为 h(n),长度为 N 的 FIR 滤波器具有
10、第一类相位特性，只需满足 即可。()2.3线性相位 FIR滤波器的幅度特性 Hg(w)的特点实质上，幅度特性特点就是 FIR 滤波器线性相位的频率域约束条件，将时域约束条件代入（2-1）式中，即可推导出 FIR 线性相位的 Hg(w)的约束条件，当 N 取奇偶不同时，Hg(w)的约束条件不同，本论文仅针对第一类相位特性进行讨论，所以只有两种情况。为了书写方面，引入一下两个参数：， 12N12M表示取不大于 的最大整数，显然，当 为奇数时 = ，当()/M()/ M为偶数是 - = 。2情况 1： ， 为奇数。()1)hnNn0(jwjwjwngnHeehe1 (1)2011()()2 2201
11、120()()2()cos()NMj jnjwNnnjwjjnNMjnhheeh 所以（2-8）10()2()cos()MgnHwhwn雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计8由于 关于 三点偶对称，因此情况 1 可以实现低通、高通、带通、cos()wn0,2带阻滤波器。情况 2： ， 为偶数()1)hNn仿照情况 1 可推导： 100()()()2()cos()MjwjwjwnjgnnHeehehwn （2-9）0()2()cos()Mgn此时 关于 偶对称，关于 奇对称，此时不能实现高通和带阻滤cosw,w波器。2.4线性相位 FIR滤波器的零极点分
12、布特点 10()()NnnHzhz将 代入上式()1hn（2-10）1001(1)(1)()NNnnnmzhzzH由（2-10）可以分析出，若 是 的零点，则 也是其零点，又因为 是zi()z1iz ()hn实序列，所以 的零点必成共轭对出现，所以 和 也是其零点。()Hz *i*()i(z)是 的 N-1 次多项式，所以它在 z 平面内有 N-1 个零点，在原点 z=0 处有 N-1 个1重极点，因此 H(z)永远稳定。健雄职业技术学院毕业设计（论文）93.FIR数字滤波器的 matlab设计3.1 通过窗函数法设计 FIR滤波器3.1.1 窗函数设计原理设希望逼近的滤波器频率响应函数为 ,
13、，其单位脉冲响应为 。()jwdHe()dhn（3-1）()jwjnddnHeh（3-2）1(2cjwnddhe如果由已知的 得出 ，经过 Z 变换可得到滤波器的系统函数，但理想器的幅频()jwde()dn响应在边界频率处往往是不连续的，因而 是无限时宽的，且是非因果序列。例如，()dhn线性相位理想低通滤波器的频率响应函数 为：jwHe(3-3)其单位脉冲响应 为：()dhn（3-4）sin()12cwcjnd wed 由上式看到，理想低通滤波器的单位脉冲响应 是无限长，且是非因果序列。()dhn的波形如图 3-1（a）所示。为了构造一个长度为 N 的第一类线性相位的 FIR 滤波器，()d
14、hn只有截取 一段，并保证截取的一段关于偶 对称。设截取的一段用() (1)/2表示，即：()d(3-5)()()dnhR雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计10式中 为一矩形序列，长度为 N，波形如图 3-1（b）所示。由该图可知，当 取()nR 值为 时，截取的一段 关于 偶对称，保证所设计的滤波器具有线1/2N()hn(1)/2性相位。我们所设计的滤波器单位脉冲响应为 ，长度为 N，其系统函数为 ，() ()Hz10NnnHzhz这样用一个有限序列 代替 ，肯定会引起误差，表现在时域就是吉布斯效应。()hn()d在(3-5)中， 就是起对无限时宽
15、非因果序列 截断作用，所以形象地将()nR()dhn称为窗函数，将 称为对 的加窗处理。()nR()dnh()dhn图 3-1 窗函数设计法的时域波形图加窗后的幅频特性如图 3-2健雄职业技术学院毕业设计（论文）11图 3-2 幅频响应函数曲线3.1.2 典型窗函数介绍1）矩形窗 ()RNwn其频率响应 为：()jwRWe1100()NNj jwnjnRne1()2si/2)(jN()jwRgWe其幅度响应 为：()Rgwsin(/2)()RgN2）汉宁窗升余弦窗 2()0.51cos()(Hn NnwR雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计12其频率响
16、应 为：()jwHnWejnFT1220.5().25()()1NjwRgRgRgWweN 12()NjwHnge其幅度响应 为：W22()0.5().25()()11HngRgRgRgwwWN 3）布莱克曼窗4()0.4250.8cos)(11Bl Nnnn其频率响应 为：jwlWe22()()1144()()11()0.20.5. jwjwj j NNBl RRRjwjNeWee其幅度响应 为：lg()BWlg 220.420.5()()114.()1RRgRgwWwNN 表 3-1 5 种窗函数基本参数类型窗函数的旁瓣峰值 n过渡带宽度 B加窗后滤波器的阻带最小衰减 srectwin 1
17、3 4/N 21bartlet 三角 25 8/N 25hanning 31 8/N 44hamming 41 8/N 53blackman 57 12/N 74健雄职业技术学院毕业设计（论文）133.1.3用窗函数法设计 FIR录波器的步骤根据窗函数设计 FIR 滤波器的步骤如下：（1） 根据阻带衰减指标，选择合适的窗函数类型，再根据过渡带指标，估算出窗口长度 N。具体参数参考表 3.1.2.1.（2） 构造希望逼近的频率响应函数 。()jwdHe很据所设计滤波器（低通、高通、带通、带阻）及通带边界频率和阻带截止频率构造希望逼近的频率响应函数。一般想滤波器的截止频率 取: , 为通cw2ps
18、p带边界频率， 为阻带截止频率。sw（3） 计算 。()dhn1()2jwnddHe（4） 加窗得到设计结果： 。()dh3.1.4通过 M 文件实 FIR数字语音低通滤波器的设计实际设计时一般使用 matlab 工具箱函数。可调用工具箱函数 fir1 实现窗函数设计。Fir1 是利用窗函数法设计线性相位 fir 数字滤波器的工具箱函数，实现线性相位 FIR 数字滤波器的标准窗函数设计。这里所指的“标准” ，是指设计低通、高通、带通和带阻FIR 滤波器时， 分别表示相应的线性相位理想低通、高通、带通和带阻滤波器的()jwdHe频率响 应函数。因而将所设计的滤波器的频率响应函数称为标准频率响应。
19、Fir1 的调用格式及功能：hn=fir1(m,wc),返回 6dB 截止频率为 的 M 阶滤波器系数向量 hn，默认选用哈明cw窗。滤波器单位脉冲响应 与 的关系为：()hn1)0,12.n其中 为对 归一化的数字角频率， 。cwcw当 =wcl,wcu得到的是带通滤波器，其-6dB 通带为 。clwu雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计14hn=fir(M,wc,ftype)，可设计高通和带阻滤波器。当 ftype=high 时，设计高通 FIR 滤波器；当 ftype=stop，且 =wcl,wcu时，设计带阻 FIR 滤波器。cwhn=fir(
20、M,wc,window)，可设计指定窗函数类型的 FIR 滤波器系数向量。缺省时默认为哈明窗。利用 matlab 设计 FIR 语音低通滤波器，技术指标要求：通带 0 3000HZ，阻带f3500HZ 最小衰减 50dB，采样频率为 fs=22050HZ。f实现程序：clear all;fs=22050Hz;f=0:fs-1;wp=3000*2*pi;ws=3500*2*pi; deltw=(ws-wp)/fs;N0=ceil(6.6*pi/deltw);N=N0+mod(N0+1,2);hn=fir1(N-1,(wp+ws)/(2*fs*pi);Hw=fft(hn,fs);abs_Hw=20
21、*log10(abs(Hw);phase=angle(Hw);subplot(3,1,1);stem(0:N-1,hn);title(单位脉冲响应);xlabel(n);ylabel(h(n);subplot(3,1,2);plot(f,abs_Hw);title(幅频特性);xlabel(f);ylabel(20lg(|H(ejw)|);axis(0 11025 -150 50);grid onsubplot(3,1,3);plot(f,phase);title(相频特性);xlabel(f);ylabel(phase);axis(0 11025 -4 4);grid on健雄职业技术学院毕
22、业设计（论文）15图3-3 窗函数法实现的FIR滤波器的响应特性有一个音乐文件，其信号频谱为：3-5 无噪音的音乐频谱插入白噪声后的频谱为：3-5 加噪音的音乐频谱雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计16利用所设计的 FIR 数字语音低通滤波器滤波后的频谱为：3-6 滤波后的音乐频谱实现程序：%-读取文件-y,fs,nbyte=wavread(c:因为爱情);M=length(y);f=(1:M)*(fs/M); %-加噪处理 -noise=0.01*randn(M,2);si=y+noise;Y=fft(y);abs_Yw=1/fs*abs(Y);S
23、i=fft(si);abs_Siw=1/fs*abs(Si);%-滤波器设-wp=3000*2*pi;ws=3500*2*pi; deltw=(ws-wp)/fs;N0=ceil(6.6*pi/deltw);N=N0+mod(N0+1,2);hn=fir1(N-1,(wp+ws)/(2*fs*pi);%-滤波处理-x=fftfilt(hn,si);Xw=fft(x);abs_Xw=1/fs*abs(Xw);%-画图部分-subplot(3,1,1);plot(f,abs_Yw);title(无噪声的信号频谱);xlabel(f/HZ);ylabel(幅度); subplot(3,1,2);pl
24、ot(f, abs_Siw);title(无噪声的信号频谱);xlabel(f/HZ);ylabel(幅度);subplot(3,1,3); plot(f,abs_Xw);title(处理后的信号频谱); xlabel(f/HZ);ylabel(幅度);健雄职业技术学院毕业设计（论文）17%-结束-3.2通过等波纹最佳逼近法设计 FIR数字滤波器等波纹最佳逼近法是一种优化设计方法，使最大误差最小化，并在整个逼近频段上均匀分布。用等文博最佳逼近法设计的 FIR 滤波器的幅频响应在整个通带和阻带都是等波纹的。而且可以分别控制通带和阻带的波纹幅度。最佳逼近是指在滤波器长度给定的条件下，使加权误差最小
25、化。由于这种设计方法使最大误差均与分布，所以设计的滤波器性价比最高。阶数相同时，这种设计是最大逼近误差最小，即通带最大衰减最小，最带最小衰减最大；指标相同时，滤波器的阶数最低。3.2.1 remez和 remezord函数滤波器指标1）remez 和 remezord 函数remez 函数实现线性相位 FIR 数字滤波器的等文波最佳逼近设计。其调用格式为：hn=remez(M,f,m,w)调用结果返回单位脉冲响应 hn。M：滤波器的阶数，f 和 m 给出希望逼近的特性。W 为误差加权量。remezord 函数可根据逼近指标估算出等文波最佳逼近 FIR 滤波器的最低阶数 M，误差加权量 w 和归
26、一化边界频率向量 f。其返回参数供 remez 函数调用。其格式为：M,f0,m0,w=remezord(f,m,rip,Fs)2）滤波器指标低通滤波器设计指标逼近通带：0， ，通带最大衰减： dB；逼近阻带： ， ，阻带最小衰p psw减： dB，remezord 调用参数：s12(,)0,pswfmrip高通滤波器设计指标逼近通带： ， ，通带最大衰减： dB；逼近阻带：0， ，阻带最小衰pwpsw减： dB，remezord 调用参数：s雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计1821(,)1,0,pswfmrip带通滤波器设计指标逼近通带： ， ，
27、通带最大衰减： dB；逼近阻带：0， ， , ，plu pslwsu阻带最小衰减： dB，remezord 调用参数：s21(,)0,1,plusl swf mrip带阻滤波器设计指标逼近通带：0， ， , ，通带最大衰减： dB；逼近阻带： ， ，plpupslwsu阻带最小衰减： dB，remezord 调用参数：s12(,)1,0,pl puslwf mrip3.2.2通 remez和 remezord函数实现 FIR数字语音低通滤波器的设计利用 matlab 设计 FIR 语音低通滤波器，技术指标要求：通带 0 3000HZ 最大衰f减 1dB，阻带 3500HZ 最小衰减 50dB，
28、采样频率为 fs=22050HZ。fFs=22025;fps=3000,3500;m=1,0;rp=1;rs=50;dat1=(10(rp/20)-1)/(10(rp/20)+1);dat2=10(-rs/20);rip=dat1,dat2;M0,f0,m0,w=remezord(fps,m,rip,Fs);M=M0+mod(M0+1,2);hn=remez(M,f0,m0,w);n=0:M;Hw=fft(hn,Fs);abs_Hw=20*log10(abs(Hw);%-画图部分-subplot(2,1,1);stem(n,hn);title(单位脉冲响应 );xlabel(n);ylabel
29、(h(n);subplot(2,1,2);f=0:Fs-1;plot(f,abs_Hw);grid on;title(幅频响应);xlabel(f/HZ);ylabel(20*lg(|Hw|);健雄职业技术学院毕业设计（论文）19图3-7 remez和remezord函数实现的滤波器响应特性分析图3-3和图3-7易得，在满足相同技术指标时，用等文波最佳逼近法设计fIR滤波器的阶数要低于用窗函数设计的阶数。3.3通过 FDATool实现 FIR数字滤波器的设计及分析FDATool（Filter Design and Analysis Tool）是一个强大的数字滤波器设计与分析工具。利用它可方便地
30、以设计出满足各种性能指标的滤波器，并可以查看滤波器各种分析图形。等待滤波器设计满意后，还可以把其系数直接导出为 matlab 变量，文本文件或 C语言头文件。在命令窗中运行 FDATool，启动滤波器分析工具，界面如下图 3.3.1 所示。图 3-8 FDATool 界面雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计20下面利用 FDATool 设计 FIR 语音低通滤波器，技术指标要求：通带 03000HZ 最大衰减 1dB，阻带 3500HZ 最小衰减 50dB，采样频率为f ffs=22050HZ。FDATool 界面设计参数如图 3-9图 3-9 FDA
31、Tool 界面参数设置设置滤波器为低通滤波器 图 3-10 滤波器类型设置用窗函数法设计 FIR 数字滤波器图 3-11 设计方法设置健雄职业技术学院毕业设计（论文）21采用哈明窗函数图 3-12 窗函数类型设置根据过渡带带宽估算滤波器阶数M=N-1，图 3-13 滤波器阶数设置采样频率 ，截止频率sF2pscf图 3-14 滤波器采样频率设置图 3-15、图 3-16、图 3-17 是利用 FDATool 所设计的 FIR 数字语音低通滤波器的响应特性。雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计22图 3-15 单位脉冲响应图 3-16 幅频特性图 3-1
32、7 相频特性由此可见利用 FDATool 工具箱可以很方便的设计我们所需的 FIR 数字滤波器，并可以生成.fda 格式的工程文件，供 simulink 仿真使用。健雄职业技术学院毕业设计（论文）234.利用 Simulink对 FIR滤波器进行仿真Simulink 是用来对动态系统进行建模、仿真和分析的软件包。使用 Simulink 来建模、分析和仿真各种动态系统，将是一种非常轻松的事情。它提供了一种图形化的交互环境，利用鼠标拖动大方法便能迅速地建立起系统框图模型。4.1Simulink模块库简介Simulink 模块库极其丰富，在本论文中不可能作全面介绍，下面只介绍在 4.2 章节中搭建语
33、音信号 FIR 低通滤波仿真模型所需要的模块库。1）信号源模块库脉冲输入（Discrete Impulse）读工作库模块（From Workspace）2）数学运算模块加法器（Add）3）信号处理模块数字滤波器设计（Digital Filter Design）4）输出模块示波器（Scope ）频谱仪（Spectrum ）扬声器（To Audio Device）4.2搭建语音信号 FIR低通滤波仿真模型通过 4.1 章节所介绍的模块库搭建如下仿真模型：雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计24图 4-1 仿真模型图通过示波器 Scope 可观察语音信号的时
34、域波形图 4-2图 4-2 音乐波形图通过频谱仪可观察到 Spectrum Scope1、Spectrum Scope2、Spectrum Scope3、可观察信号加噪前后及处理后的频谱图。健雄职业技术学院毕业设计（论文）25图 4-4 Spectrum Scope1 所示信号频谱图 4-5 Spectrum Scope2 所示信号频谱雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计26图 4-6 Spectrum Scope3 所示信号频谱健雄职业技术学院毕业设计（论文）27总结滤波是信号处理中最基本有极为重要的技术，利用滤波器技术可以从复杂的信号中提取出所需要
35、的信号，抑制不需要的信号。FIR 滤波器的优点FIR 滤波器能严格做到线性相位或群延时为常数，而 IIR 滤波器只能逼近线性相位。FIR 滤波器是全零点型滤波器，总是稳定的，不会因滤波运算的舍入误差而产生极限环振荡现象对同样幅度相应的滤波器，用 FIR 滤波器实现比用 IIR 滤波器实现需要较高的阶数，多达 510 倍。当滤波器的特性要求较高时，用 FIR 滤波器来实现，滤波过程需要较多的计算时间Matlab 软件作为强大的计算软件， Matlab 在滤波器设计方面还提供了编写图形用户界面的功能。所谓图形用户界面, 简称为 GUI, 是由各种图形对象, 如图形窗口菜单按钮、文本框等构建的用户界
36、面，它是用户和计算机之间交流的工具。Matlab 将所有 GUl 支持的用户控件都集成在这个环境中并提供界面外观、属性和行为响应方式的设置方法,随着版本的提高,这种能力还会不断加强。而且具有强大的绘图功能,可以轻松的获得更高质量的曲线图。Matlab 语言进行程序设计和利用信号处理工具箱的 FDATool 工具进行界面设计可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。还可在 MTLAB 环境下的仿真软件 Simulink 对所设计的滤波器进行模拟仿真。雍鹏亮 基于 Matlab/SIMULINK 的 线性相位 FIR 滤波器设计28致谢毕业意
37、味着一个人 一个阶段学习生涯的结束。在大学里，毕业论文是宣告这一事实的标准。从大一到现在3年的学习3年的磨练，在此刻沉淀成一篇毕业论文。通过这次毕业设计，使我受到了一次用专业知识、专业技能分析和解决问题的全面系统的锻炼。我学到了很多书本上学不到的实践知识，尤其是对电子硬件的设计有了更深的理解，了解了设计一个产品从方案的选择到最终投产的一系列的设计阶段，理解了一般的产品调试技术和焊接技术，我真正体会到理论和实践相结合的重要性。再此有无数的感谢！感谢我的家人和朋友，他们的支持和鼓励使我顺利完成学业。感谢教了我3年的老师，那些曾经的岁月，曾经的年华，我们一起走过的路。感谢我的同学们，无论是教室里的切
38、磋，餐桌旁的高谈阔论，还是寝室里的欢声笑语，都曾留下我们美好的足迹。这是我最可宝贵的财富之一。当然，这里最要感谢的是张宏杰老师，他是我毕业设计的指导老师。衷心感谢他在这段时间对我的帮助和对我的指导！张老师的言传身教，我将谨记于心，老师广博的学识和严谨的治学态度使我受益终生！健雄职业技术学院毕业设计（论文）29参考文献1 丁玉美,高西全.数字信号处理M. 3 版. 西安：西安电子科技大学出版社，2008.2刘学. MATLAB/SIMULINK 通信系统建模与仿真M. 北京：电子工业出版社，2011.3唐向宏.MATLAB 在电子信息类课程的应用M.北京：电子工业出版社,2006.4姚剑清.简明数字信号处理M.北京：人民邮电出版社，2008 .5 董长虹等. MATLAB信号处理与应用M.北京：国防工业出版社，2005.6 楼顺天.李博菡.基于MATLAB的系统分析与设计M.西安：西安电子科技大学出版社，1998.7 张葛祥,李 娜. MATLAB仿真技术与应用M.北京：清华大学出版社，2003.8 刘兴钊.数字信号处理M.北京：电子工业出版社，2010.
