FIR最小方差法带阻滤波器.docVIP

南华大学 数字信号处理课程设计 学院： 电气工程学院 学生姓名： 王明明 学号： 20094470149 题目编号： 0904 设计题目： 基于最小方差法的数字频变带阻滤波器 指导老师： 陈忠泽 2013年1月 1．课程设计的内容和要求（包括原始数据、技术要求、工作要求等）： 1.设计内容： 根据自己在班里的学号49（最后两位）查表一得到一个四位数，0904，由该四位数索引表二确定待设计数字滤波器的类型：最小方差FIR带阻（数字频带变换）滤波器 2.滤波器的设计指标： ⑴通带下截止频率 ⑵阻带下截止频率 ⑶阻带上截止频率 ⑷通带上截止频率 ⑸通带最大衰减 ⑹阻带最小衰减 其中，—你的学号的最后两位 3. 滤波器的初始设计通过手工计算完成； 4. 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器结构进行分析）； 5. 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 6. 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 7. 课程设计结束时提交设计说明书。  2．对课程设计成果的要求〔包括图表（或实物）等硬件要求〕： 滤波器的初始设计通过手工计算完成； 在计算机辅助计算基础上分析滤波器结构对其性能指标的影响（至少选择两种以上合适的滤波器结构进行分析）； 在计算机辅助计算基础上分析滤波器参数的字长对其性能指标的影响； 以上各项要有理论分析和推导、原程序以及表示计算结果的图表； 课程设计结束时提交设计说明书。 3．主要参考文献： [1]高息全 丁美玉.《数字信号处理》[M].西安：西安电子科技大学出版社，2008.8 [2]陈怀琛.《数字信号处理教程——MATLAB释义与实现》[M].北京：电子工业出版社，2004.12 [3]张德丰.《详解MATLAB数字信号处理》[M].北京：电子工业出版社，2010.6 [4]飞思科技产品研发中心.《MATLAB7辅助信号处理技术与应用》[M].北京：电子工业出版社，2005.3 4．课程设计工作进度计划： 序号 起 迄 日 期 工 作 内 容 1 2012.12.26-2013.12.31 接到题目，搜集资料 2 2012.12.31-2013.1.3 整理资料，构思设计方案 3 2013.1.3-2013.1.5 手工计算进行滤波器的初步设计 4 2013.1.5-2013.1.7 完善初步设计，学习Matlab软件操作 5 2013.1.8-2013.1.12 通过Matlab软件分析设计内容，逐步落实课题目标 6 2013.1.13-2013.1.14 上交课程设计，并做细节修改并完成设计 主指导教师 日期： 年 月 日 一 手工计算完成最小方差FIR带阻滤波器初始设计 带阻滤波器的设计指标： =49 通带下截止频率, =0.53289πrad*0.6=0.319734πrad=1.004474rad/sample 带下截止频率, =0.79934πrad*0.6=0.479604πrad=1.506720rad/sample 带上截止频率, =1.86512πrad*0.6=1.11907πrad=3.515662rad/sample ⑷通带上截止频率, =2.13156πrad*0.6=1.27893πrad=4.017877rad/sample ⑸通带最大衰减, ⑹阻带最小衰减 由得到性能指标在MATLAB中的常用形式 再将其除以采样频率Fs转换为归一化频率： 通带下截止频率：0.159867πrad/sample 阻带下截止频率：0.239802πrad/sample 阻带上截止频率：0.55954πrad/sample 通带上截止频率：0.63947πrad/sample 将ap=1dB，as=60dB带入公式 ξ1=10^(ap/20)-1)/(10^(ap/20)+1， ξ2=10^(-as/20) 得 ξ1=0.057501128 ξ2=0.001 由凯泽提出逼近n的公式 求得 Ws-Wp=（0.479604π-0.319734π）+（1.27893π-1.11907π）=0.318730π Δ =0.159365 N=56 3．在Matlab中利用REM