武科大信号处理综合课程设计技术分析.doc

二○一四～二○一五学年第 一 学期 电子信息工程系 课程设计报告书 班 级： 电子信息工程 学 号： 姓 名： 课程名称： 信号处理综合课程设计 学时学分： 2 周 1 学 分 指导教师： 杨 永 立 二○一五 年 一月 信号谱分析和滤波器设计 一、设计任务和要求 有一个含有5Hz、15Hz和000Hz的混合正弦波信号，试设计滤波，。 要求： ① 确定采样频率、采样长度，对进行采样； ② 分别用IIR滤波器和FIR滤波器； ③ 分析并画出分离前后信号的频谱，画出滤波器的频率响应，评价分离的效果，评价设计的滤波器是否满足要求。 二、原理分析和设计 混合正弦波信号 =sin(2π500t)+sin(2πt)+sin(2πt) 其包含500Hz、1500Hz和5000Hz频率成分，如图1所示。 图1 混合正弦波信号及频谱图 2、由于最高频率为5000Hz，故采样频率应大于10000Hz，实际设计中考虑到余量，采样频率可以取为20000～30000Hz即可。为显示出更好的滤波效果，本设计采用30000Hz频率进行采样。 三、IIR滤波器的设计 （1） 第一步，为分离出500Hz的正弦波成分，可以考虑采用低通滤波器，屏蔽掉1500Hz和5000Hz频率成分。在这里，我采用了巴特沃斯低通滤波器的设计方法，MATLAB中提供了用于设计巴特沃斯原型低通滤波器的函数，其调用格式如下： [N,Wn]=buttord(wp,ws,rp,rs,’s’):其中，wp和ws分别是通带边界频率和阻带边界频率，单位为rad/s；rp和rs分别为通带最大衰减和阻带最小衰减；返回的参数N和Wn分别为滤波器的阶数和3dB截止频率。 为实现500Hz频率分量的分离，可将通带边界频率设置为wp1=2π600rad/s，阻带边界频率ws1=2π×1500rad/s，通带最大衰减rp=1dB，阻带最小衰减rs=60dB。并采用冲激响应不变法将模拟滤波器转换成数字滤波器。巴特沃斯低通滤波器频率响应图如下： 图2 巴特沃斯低通滤波器频率响应 同时将横轴归一化为频率Hz。混合正弦波经巴特沃斯低通滤波器后便可分理出500Hz正弦波信号，其效果如图3所示： 图3 500Hz频率分量的时域和频谱图 （2） 第二步，为分离出1500Hz的正弦波成分，可以考虑采用带通滤波器，屏蔽掉500Hz和5000Hz频率成分。 设置带通滤波器通带下限频率和上限频率分别为1300Hz、1700Hz，中心频率为1500Hz，带宽B=400Hz；阻带下限频率和上限频率分别为600Hz、4000Hz；通带最大衰减为1dB，阻带最小衰减为100dB，采用巴特沃斯带通滤波器设计。并采用冲激响应不变法将模拟滤波器转换成数字滤波器。巴特沃斯带通滤波器频率响应图如下： 图4 带通滤波器频率响应 同时将横轴归一化为频率Hz。混合正弦波经巴特沃斯带通滤波器后便可分理出1500Hz正弦波信号，其效果如图5所示： 图5 1500Hz频率分量的时域和频谱图 （3） 第三步，为分离出5000Hz的正弦波成分，可以考虑采用高通滤波器，屏蔽掉500Hz和1500Hz频率成分。 在高通滤波器设计中，将通带截止频率设置为5000Hz，阻带截止频率为4000Hz，并采用双线性变换法设计数字滤波器；通带最大衰减为1dB，阻带最小衰减为60dB。巴特沃斯高通滤波器频率响应图如下： 图6 高通滤波器频率响应 同时将横轴归一化为频率Hz。混合正弦波经巴特沃斯高通滤波器后便可分理出5000Hz正弦波信号，其效果如图7所示： 图7 5000Hz频率分量的时域和频谱图 四、FIR滤波器的设计 MATLAB提供了函数w=hamming（N）’ftype’,Window（N））2的向量；设计低通和高通滤波器时，无需’ftype’，Window表示设计滤波器所采用的窗函数类型，Window的长度为N，若Window缺省，则fir1默认使用哈明窗；h对应设计好的滤波器的系数h(n)。 （1） 滤除500Hz频率成分，设计FIR低通滤波器，取截止频率wp1= 0.0191π，窗宽N=180。FIR低通滤波器幅频和相频特性如图8所示，分离效果如图9所示。 图8 FIR低通滤波器幅频和相频特性 图9 500Hz频率分量的时域和频谱图 （2）滤除1500Hz频率成分，设计FIR带通滤波器，取截止频率wp2=[0.066,0.16];，窗宽N=200。FIR带通滤波器幅频和相频特性如图10所示，分离效