Ch3.5 FIR滤波器的设计.ppt

Ch3.5: FIR滤波器的设计 FIR滤波器的数学表示 H(z)是z-1的N-1阶多项式, 在z平面有(N-1)个零点, z=0是(N-1)阶(重)极点, 收敛域？稳定系统？ FIR滤波器的特点 FIR的两大特点： FIR filter的h(n)有限长，只有零点，无极点，因而是稳定的； FIR filter在满足一定条件下，做到具有严格的线性相位； 适用于高保真度的信号处理： 图像处理以及数据传输, 语音或音频信号处理, FIR滤波器可用DFT实现。 设计思路: 确定指标: 按照任务要求, 确定滤波器的性能指标, 逼近: 用一个因果稳定的离散LTI系统的系统函数来逼近这一性能要求, 实际的实现: 软件: 利用有限精度算法来实现这个系统函数, 硬件: 专用数字滤波器器件, 软硬件的结合, 设计的方法 对理想滤波器频率特性作某种近似, 窗函数设计法 从时域出发，用窗函数把理想的hd(n)截取成有限长的h(n)去逼近理想的hd(n) 频率抽取法 从频域出发，把理想的Hd(j?)抽样成H(k), 进而得到滤波器的H(z) 数值逼近法(最优等波段设计法) 插值法，最小平方近似法，一致逼近法 滤波器的指标 幅度平方响应(幅频响应) 反应了信号x(n)通过该系统后各频率成分衰减的情况 相位响应(相频响应) 反映信号x(n)中各成分通过系统后在时间上的位移情况 群延迟 反映滤波器的平均延迟, 定义为相位对角频率的导数的负值, 一个理想的离散时间系统, 具有: 所希望的幅频响应, 线性相位, 通带内群延迟特性是常数, Linear phase 无失真传输条件，信号按比例放大和时间上有一定的延迟。y(n)=Ax(n-k) 线性系统不失真的条件: 幅频特性为一常数， 具有线性相位，why? Suppose arg[H(?)]=-k?, 输入和输出关系如何？ An example A signal having two sinusoids is filtered with phase β1 and β2. What is the condition when all sinusoids have the same time delay L? FIR滤波器线性相位的条件 Question: condition for undistorted transfer? 如果FIR滤波器的h(n)为实数, 而且满足以下任意条件, 滤波器就具有准确的线性相位, 第一种：偶对称，h(n)=h(N-1-n)， φ(ω)=-(N-1) ω /2 第二种：奇对称，h(n)=-h(N-1-n), φ(ω)=-(N-1) ω /2 +pi/2 对称中心在n=(N-1)/2处, 四种线性相位FIR滤波器 线性相位的特性 h(n)偶对称的H(ω) h(n)奇对称的H(ω) Four different filters h(n), evenly symmetric, odd N, h(n), evenly symmetric, even N, h(n), oddly symmetric, odd N, h(n), oddly symmetric, even N, The above four filters have different frequency response and are thus good for different purposes. 零点位置 零点具有一定的对称性, h(n)为实数的零点构象的四件套: 窗函数设计法(傅立叶级数法) 设计方法: 理想的Hd(ej?)-- hd(n), FIR的h(n)逼近hd(n), 低通滤波器: 理想的矩形频率特性Hd(ej?), hd(n)为无限长的, 非因果系统， 用有限长的窗口函数序列w(n)来截断hd(n), h(n)=w(n)?hd(n), H(ej?)=Hd(ej?)*W(ej?) 窗函数的设计思想 根据给定的滤波器指标, 选择滤波器的长度N和窗函数w(n), w(n)应具有最窄的主瓣和最小的旁瓣, 矩形窗 矩形窗的特点 窗函数为有限序列N, 其响应有一个主瓣, 宽度与1/N成正比, 2*pi/N, 卷积使理想响应变得轮廓模糊, 主瓣在H(ej?)中产生过渡带, 其宽度由主瓣的宽度决定, 主瓣越宽, 过渡带越宽, 旁瓣产生波动, 在通带和阻带中有类似的形状, 即通带和阻带的波动系数相同, 对窗函数的要求 两项要求： 窄的主瓣，以获得较陡的过渡带, 减小窗谱的最大旁瓣的相对幅度，即能量集中于主瓣，使肩峰和波纹减少, 二者不能同时满足, 选择窗函数 增加主瓣宽度以换取对旁瓣的抑制，即得到平坦的幅度响应和较小的阻带波纹， 矩形窗的两个缺点 第一旁瓣的衰减为13dB,