理论课第7章yyxMATLAB在数字信号处理中的应用剖析.ppt

* 例7.24 系统函数的计算机推导 数字滤波器的网络结构图实际上也是一种信号流图。因此【例6.20】中介绍的方法和公式同样可以用来求离散域的数字滤波器的系统函数。不同的地方仅仅在于节点方程中出现了作为系数的符号变量z?1，它将出现在系数矩阵中。MATLAB是不能处理上标变量的，因此在程序中设q=z ?1，在计算完成后再人工地把结果中的q恢复为z?1。 * 例7.24的结构图与方程 * 例7.24的方程的矩阵形式 由此可以求出系统函数 * 例7.24解出的系统函数 程序运行的结果 如果加入一个激励x(n)?A?(n)，则得出 * 7.5 FIR数字滤波器设计 滤波器的特性指标 用绝对值δ1,δ2表示; 用分贝Rp,Rs表示 * （1）窗函数法设计FIR滤波器 先根据?c和N求出相应的理想滤波器单位脉冲响应hd(n)。 第二步要选择合适的窗函数w(n)来截取hd(n)的适当长度（即阶数），以保证实现要求的阻带衰减；最后得到FIR滤波器单位脉冲响应h(n)=hd(n).*w(n),即其系数。 * （2）等波纹最佳一致逼近法 （2）等波纹最佳一致逼近法: 信号处理工具箱采用remez算法实现线性相位FIR数字滤波器的等波纹最佳一致逼近设计。其优点是，设计指标相同时，使滤波器阶数最低；或阶数相同时，使通带最平坦，阻带最小衰减最大；通带和阻带均为等波纹形式，最适合设计片段常数特性的滤波器。其调用格式如下： ? b = remez(N, f, m, w, ftype) 其中N由remezord函数求出: ? [N, fo, mo, w] = remezord(f, m, dev, Fs) 输入变元dev为各逼近频段允许的波纹振幅。 remez函数可直接调用remezord返回的参数如下： ? b=remez(N, fo, mo, w) * 例7.25 窗函数法设计数字滤波器 分别用矩形窗和Hamming窗设计线性相位FIR低通滤波器。要求通带截止频率?c = ?/4，单位脉冲响应h(n)的长度N = 21。绘出h(n)及其幅频响应特性曲线。 先求出相应的理想滤波器（本例应为理想低通）单位脉冲响应hd(n)，再根据阻带最小衰减选择合适的窗函数w(n)，最后得到FIR滤波器单位脉冲响应h(n)=hd(n).*w(n)。 * 例7.25 窗函数法设计数字滤波器 本题中，?c = ?/4，N = 21，所以线性相位理想低通滤波器的单位脉冲响应为: 为了满足线性相位FIR滤波器条件h(n) = h(N-1-n)，要求? = (N-1)/2 = 10。 信号处理工具箱中有窗生成函数boxcar，hamming，hanning和blackman等。 * 例7.25 窗函数法设计数字滤波器 对两种窗函数的设计结果分别如右图7.25-1和图7.25-2所示。 * 工具箱设计函数fir1和fir2 MATLAB提供了基于窗函数法的FIR滤波器设计函数fir1和fir2,其功能及用法如下。 ? fir1功能：标准频率响应形状。 格式：b=fir1(N, wc, ‘ftype’, window)。 当wc=[wc1，wc2]时，是的带通滤波器。 当ftype=high时，设计高通FIR滤波器； 当ftype=stop时，设计带阻FIR滤波器。 ? fir2功能：任意频率响应形状。 格式：b = fir2(N, f, m, window) * 例7.26 窗函数法设计带通滤波器 使用fir1函数b = fir1(N, wc, window) 编程 参数?c为行向量?c = [?lp/?，?hp/?] 根据阻带最小衰减Rs = 60dB选择窗函数类型和阶次。可以查上面列出的“窗函数设计滤波器时的阶数选择表”。选blackman窗，其滤波器阻带最小衰减可达到74dB，其窗口长度M由过渡带宽度B = 0.15? 决定，Blackman窗设计的滤波器过渡带宽度为12?/M，故M取80。因M = N+1，所以滤波器阶数N = 79。 * 例7.27 用remez函数低通滤波器 解：先由题意计算设计参数 f = [1/4,5/16]，m = [1,0]; dev的计算稍复杂一些，由于 所以 有了这几个参数就可以调用remezord和remez函数了. * 例7.27 用remez函数低通滤波器 横线为-3dB，两条竖线分别位于频率?/4和5?/16。显然，通带指标稍有富裕，过渡带宽度和阻带最小衰减刚好满足指标要求。 程序输出的幅频特性 * 例7.28 remez函数设计高通滤波器 观察等波纹逼近法中加权系数w(?)及滤波器阶数N的作用和影响。期望逼近的滤波器通带为[3?/4，?]，阻带为[0，23?/