dsp课程报告【DOC精选】.doc

计算机科学技术系 课程设计报告 课 程 数字信号处理 专 业 通信工程 班 级 2009级本科一班 学生姓名1 陈博俊 学号2009010954 学生姓名2 江 英 学号2009010944 学生姓名3 万永辉 学号2009010917 学生姓名4 王燕龙 学号2009010933 二〇一一年十二月目录 目录 I 1.使用双线性变换设计IIR数字滤波器 - 1 - 1.2设计要求 - 1 - 1.2设计过程 - 1 - 1.3 MATLAB工作脚本 - 4 - 1.4运行结果 - 7 - 2. 使用Blackman窗法设计FIR高通滤波器 - 9 - 2.1设计要求 - 9 - 2.2设计过程 - 10 - 2.3 MATLAB工作脚本 - 10 - 2.3运行结果 - 12 - 参考文献 - 15 - 摘要：通过本次课程设计，主要训练和培养学生综合应用所学数字信号处理课程的滤波器设计的相关知识，独立完成滤波器的设计和MATLAB的仿真实现的能力。通过分析所给规定指标，筛选合适的滤波器，自己编写滤波器的功能函数，用双线性变换法设计IIR数字低通滤波器，以及使用Blackman窗法设计FIR高通滤波器。并用MATLAB进行仿真。 关键词：滤波器；双线性变换；IIR数字滤波器；Blackman窗法；FIR高通滤波器。1.使用双线性变换设计IIR数字滤波器 1.2设计要求 用双线性变换设计巴特沃兹滤波器，技术指标是： 用绝对指标表示为： 1.2设计过程 双线性变换法基本方法是应用下面公式： 由于双线性变换时将s平面的虚轴映射为z平面的单位园，所以产生非线性压缩，需要对离散时间滤波器的频率特性做预畸变，即确定?p、?s，公式如下： 或 对于本设计，对给出的离散时间的频域特性进行变换： 或 在采用MATLAB设计滤波器的时候，经常使用的指标是对应（dB）指标要求Rp和Rs，其中： Rp 是以dB计的通带波纹； Rs是以dB计的阻带衰减。 在本设计中计算方法为： 计算Butterworth低通滤波器的阶数N和归一化的特征频率Wc根据对Buttord的理解自己编写。部分代码： ksp=((10^(0.1*as)-1)/(10^(0.1*rp)-1))^(1/2); osp=ws/wp; N=ceil(log10(ksp)/log10(osp)) wc=wp*(10^(0.1*rp)-1)^(-1/(2*N)); 通过以上的代码实现MATLAB中buttord函数的功能。以上求的的滤波器阶数，和归一化特征频率Wc可用来确定模拟滤波器，根据以下公式计算： （此处?c以上所求的归一化特征频率Wc）① k=0,1,......,N-1 ② ③ 综合②③式，得到： k=0,1,......,N-1 ④ ⑤ 通过butter函数求的A序列为分母的系数，B为分子，所以，此函数可用以下代码实现： syms s; h=1; for k=0:(N-1) s0=wc*exp(j*pi*(1/2+(2*k+1)/(2*N))); h=(s-s0)*h; end A=real(sym2poly(h)) B=wc^N 以上步骤的意义为设计了一个低通模拟滤波器，如果将上式⑤中参数s改成z则设计相应的数字滤波器。我们知道从s变换z的公式为； ⑥ 将⑥式代入①中： 以此求的H(z)的分母，分子系数a,b,此部分可用代码实现为： syms z; h=1; for k=0:(N-1) z0=wc*exp(j*pi*(1/2+(2*k+1)/(2*N)))-2*fs; h=((wc*exp(j*pi*(1/2+(2*k+1)/(2*N)))+2*fs)*z+z0)*h; end b=fliplr(real(sym2poly(h)))*(-1)^N; l=B*(z+1)^N; a=fliplr(sym2poly(l)); 到这就求的了低通数字滤波器的各项参数，通过filter函数即可让输入信号的高频部分滤除。 1.3 MATLAB工作脚本 %%%%