第二周(3)第五章有限冲击响应滤波器课件.pptVIP

5.6频率采样 * 加窗设计法需要从时域变换到频域，为了避免此转换，可采用更直接的频率采样设计法。它直接指定滤波器响应H(ej?)的各项与长度为N的Hd(ej?)，在围绕单位圆的统一空间频率相对应，并可通过过滤带中的一些自由项减少超调，使得两者之间的误差按某种准则为最小。 频率采样法得到的滤波器一般比用窗设计得到的滤波器更加有效。 * 超调和误差变小 * 频 率采样法从频域出发，把给定的理想频率响应Hd(ejw) 加以等间隔采样得到Hd(k) 即 * 1.窗函数法与频率抽样法的区别 1、窗函数法是从时域出发，把理想的hd(n) 用一定形状的窗函数截取成有限长 的 h(n)，以此h(n)求近似理想的hd(n)，这样得到的频率响应H(ejw)逼近于所求的 理想频率响应Hd(ejw) 。 * 2、频率抽样法则从频域出发，把给定的理想频率响应Hd(ejw)加以等间隔抽样然后以此Hd(k)作为实际FIR滤波器的频率特性的抽样值H(k)，知道H(k)后，由DFT定义，可用频域的这N个抽样值H(k)求唯一确定的有限长序列h(n)，利用这N个频域抽样值 H(k)同样可得FIR滤波器的系统函数H(z)及频率响应。 * 谢 谢! Thank you! * 放大电路的电压放大倍数与频率的关系称为幅频特性,输出信号与输入信号的相位差与频率之间的关系称为相频特性。两者统称频率特性。 * L=2,3,4,5的滤波器权系数 * * 5.3陷波滤波器 * 普通生物医学信号处理问题包括从信号中移去特殊频率或者频率范围（如60/50Hz的工频干扰）的噪声，以便让更高或者更低的频率可无衰减地通过，完成这一任务的滤波器称为陷波（带阻）滤波器。 为了从信号中完全去除指定频率噪声，一种简单的方法是在单位圆上对应于该频率的位置设置一个零点。例如，用180Hz的频率采样，则位于2π/3处设置的零点，可以抑制信号中60Hz工频干扰。 差分方程： * 滤波器零点： 其幅值和相位响应为 信号流图和零极点分布图 * * 该滤波器的慢收敛性使得高于60Hz的信号也被极大消弱 5.4求导数 * 对于数字微分，由于频率响应的周期性，要得到导函数的响应是不可能的。这里所讨论的是在计算的复杂性与真实导数的接近程度和高频噪声的删除之间折中考虑的方法。 几种不同的微分算法：两点差分、三点中心差分、最小平方多项式拟合的信号流图和零极点分布图如下。 * 求导数-两点差分 差分方程： 最简单的方法 * 其幅值响应与理想响应非常接近，但由于在fs/2处未趋向于零，故将大大增强高频噪声，所以经常用一个低通滤波器紧随其后，以消除高频噪声的影响。 求导数-三点中心差分 * 差分方程 * 对于fs/10以上的频率信号，其求导差别较大，但由于在fs/2处趋向于零，故该滤波器具有一些内在的平滑作用。 求导数-五点最小平方多项式差分 * 差分方程 这种滤波器与抛物线平滑滤波器相似，是取抛物线中点的多项式斜率作为导数值。具有四个零点z=1，-1，-0.25?j0.968。 * 由于抛物线拟合的平滑特性大大消弱了高频段的响应，故响应仅在低频段，与实际求导相近。 二次求导 * 差分方程 传递函数 4个零点,幅值放大系数1 5.5窗的设计 * 通常希望频率响应Hd(?)（连续函数）的离散FFT反变换的结果就是预期单位脉冲序列hd(k)（离散函数）,但此序列有无限项，物理上无法实现，故采用加窗的技术。窗设计的目的在于选择有限项的实际h(k),使其频率响应H(ej?)从某种意义上说与Hd(?)相近。 要使实际频率响应与预期频率响应之间的平均方差最小，直接截断hd(k) (脉冲响应h(k)应选择hd(k)的前N项)，即可达到这一目标。当然，这样的截断将导致频率响应急剧转折时的大超调(Gibb现象)。 可看成希望的单位脉冲序列与一矩形窗函数相乘而得。 * 1、设计方法 w(n)：窗函数序列 要选择合适的形状和长度 * 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 其理想单位抽样响应： 中心点为 α 的偶对称无限长非因果序列 * * 取矩形窗： 则FIR滤波器的单位抽样响应： 按第一类线性相位条件，得 * * 加窗处理后对频率响应的影响： 时域乘积相当于频域卷积 而矩形窗的频率响应： * 理想滤波器的频率响应： 其幅度函数： 则FIR滤波器的频率响应： * * 幅度函数： * 加窗函数的影响： 不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。 在 处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少 改变N只能改变