TMSCx算法编程概要.ppt

FIR滤波器 IIR滤波器 快速傅里叶变换(FFT) 定标与溢出处理 7.1 FIR滤波器 数字滤波器是DSP的基本应用，有2种基本类型： 有限冲激响应滤波器FIR 无限冲激滤波器IIR 一般来说，如果需要线性相位则选择用FIR滤波器，对于相位要求不敏感的场合可以选用IIR滤波器。 本节主要讨论FIR滤波器的DSP实现方法，有关IIR滤波器的实现将在下一节中介绍。 7.1.1 FIR滤波器的基本结构 一个FIR滤波器的输出序列和输入序列之间的关系，满足差分方程： 传递函数为 7.1.2 滤波器设计方法 设计要求：运用双线性变换法基于MATLAB设计一个IIR带通滤波器，其中带通的中心频率为ωp0=0.5π，;通带截止频率ωp1=0.4π,ωp2=0.6π;通带最大衰减αp=3dB;阻带最小衰减αs=15dB;阻带截止频率ωs2=0.7π (1)根据任务,确定性能指标:在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标: 带通滤波器的阻带边界频率关于中心频率ωp0几何对称，因此ws1=wp0- (ws2-wp0)=0.3π 通带截止频率wc1=0.4π,wc2=0.6π;阻带截止频率wr1=0.3π，wr2=0.7π；阻带最小衰减αs=3dB和通带最大衰减αp=15dB; （2）用Ω=2/T*tan(w/2)对带通数字滤波器H(z)的数字边界频率预畸变，得到带通模拟滤波器H(s)的边界频率主要是通带截止频率ωp1,ωp2;阻带截止频率ωs1，ωs2的转换。 为了计算简便，对双线性变换法一般T=2s 通带截止频率 wc1=(2/T)*tan(wp1/2)=tan(0.4π/2)=0.7265 wc2=(2/T)*tan(wp2/2)=tan(0.6π/2)=1.3764 阻带截止频率wr1=(2/T)*tan(ws1/2)=tan(0.3π/2)=0.5095 wr2=(2/T)*tan(ws2/2)=tan(0.7π/2)=1.9626 阻带最小衰减αs=3dB和通带最大衰减αp=15dB; （3）运用低通到带通频率变换公式λ=(((Ω^2)-(Ω0^2))/(B*Ω))将模拟带通滤波器指标转换为模拟低通滤波器指标。 B=wc2-wc1=0.6499 normwr1=(((wr1^2)-(w0^2))/(B*wr1))=2.236 normwr2=(((wr2^2)-(w0^2))/(B*wr2))=2.236 normwc1=(((wc1^2)-(w0^2))/(B*wc1))=1 normwc2=(((wc2^2)-(w0^2))/(B*wc2))=1 得出，normwc=1，normwr=2.236 模拟低通滤波器指标：normwc=1，normwr=2.236，αp=3dB，αs=15Db (4)设计模拟低通原型滤波器。用模拟低通滤波器设计方法得到模拟低通滤波器的传输函数Ha(s);借助巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等。 (5)调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。 (6)利用双线性变换法将模拟带通滤波器Ha(s)转换成数字带通滤波器H(z) FDATOOL辅助设计滤波器系数 1.2 打开MATLAB的FDATool MATLAB集成了一套功能强大的滤波器设计工具FDATool（Filter Design Analysis Tool），可以完成多种滤波器的设计、分析和性能评估。 单击MATLAB主窗口下方的“Start”按钮，如图B.1所示，选择菜单“ToolBox” →“Filter Design” →“Filter Design Analysis Tool（FDATool）”命令，打开FDATool，如图B.2所示。 滤波器类型（Filer Type）为低通（Low Pass） ● 设计方法（Design Method）为FIR，采用窗函数法（Window） ● 滤波器阶数（Filter order）定制为15 ● 窗口类型为Kaiser，Beta为0.5 ● FS为48kHz，FC为10.8kHz 最后单击Design Filter图标，让MATLAB计算FIR滤波器系数并作相关分析。 1.4 滤波器分析 计算完FIR滤波器系数以后，往往需要对设计好的FIR滤波器进行相关的性能分析，以便了解该滤波器是否满足设计要求。分析操作步骤如下： 选择FDATool的菜单“Analysis”→“Magnitude Response”，启动幅频响应分析如图B.4所示，x轴为频率，y轴为幅度