数字信号处理课后答案第6章(高西全丁美玉第三版).pptVIP

　　%程序ex612.m 　　wp=［0.25， 0.45］； ws=［0.15， 0.55］； Rp=3； 　　As=15； %设置带通数字滤波器指标参数　　［N， wc］=buttord(wp， ws， Rp， As)； 　　　　%计算带通滤波器阶数N和3 dB截止频率Wc　　［B， A］=butter(N， wc)； 　　　　 %计算带通滤波器系统函数分子分母多项式系数向　　　　　量A和B　　myplot(B， A)； 　　　%调用自编绘图函数myplot绘制带通滤波器的损耗函　　　　数曲线 题12解图 　　13*． 设计巴特沃斯数字带通滤波器， 要求通带范围为0.25π rad≤w ≤0.45π rad, 通带最大衰减为3 dB， 阻带范围为0≤w≤0.15π rad和0.55π rad≤w ≤π rad， 阻带最小衰减为40 dB。 调用MATLAB工具箱函数buttord和butter设计， 并显示数字滤波器系统函数H(z)的系数， 绘制数字滤波器的损耗函数和相频特性曲线。 这种设计对应于脉冲响应不变法还是双线性变换法？ 　　解： 调用函数buttord和butter设计巴特沃斯数字带通滤波器程序ex613.m如下：  　　%程序ex613.m 　　wp=［0.25， 0.45］； ws=［0.15， 0.55］； 　　rp=3； rs=40；  　　［N， wc］=buttord(wp， ws， rp， rs)；  　　［B， A］=butter(N， wc) 　　clf； mpplot(B， A， rs) 　　由图可见， 该数字滤波器近似为低通滤波器。 且T越小， 滤波器频率混叠越小， 滤波特性越好（即选择性越好）。 反之，T越大， 极点　　　　　　　离单位圆越远， 选择性越差， 而且频率混叠越严重， ω=π附近衰减越小， 使数字滤波器频响特性不能模拟原模拟滤波器的频响特性。  　　7． 假设某模拟滤波器Ha(s)是一个低通滤波器， 又知 　　　　　　　， 数字滤波器H(z)的通带中心位于下面哪种情况？并说明原因。  　　（1） ω=0（低通）。  　　（2） ω=π（高通）。  　　（3） 除0或π以外的某一频率（带通）。 　　解： 方法一 按题意可写出 故 即 原模拟低通滤波器以Ω=0为通带中心， 由上式可知， Ω=0时， 对应于ω=π， 故答案为(2)。 　　方法二 找出对应于Ω=0的数字频率ω的对应值即可。 　　令z=1， 对应于ejω=1， 应有ω=0， 则 对应的不是模拟低通中心频率， 所以， 答案（1）即ω=0（低通）不对。  　　令z=－1， 对应于ejω=－1， 应有ω=π， 则 即将Ω=0映射到ω=π处， 所以答案为(2)。  　　方法三 直接根据双线性变换法设计公式及模拟滤波器由低通到高通频率变换公式求解。 双线性变换设计公式为 当T=2时， 则H(z)亦为低通。 如果将Ha(s)变换为高通滤波器： ， 这时，如果Ha(s)为低通， 则可将Hah(s)用双线性变换法变成数字高通； 这正是题中所给变换关系，所以数字滤波器 通带 中心位于ω=π， 故答案（2）正确。 　　8． 题8图是由RC组成的模拟滤波器， 写出其系统函数Ha(s)， 并选用一种合适的转换方法， 将Ha(s)转换成数字滤波器H(z)， 最后画出网络结构图。 题 8 图 解: 模拟RC滤波网络的频率响应函数为 显然， Ha(jΩ)具有高通特性， 用脉冲响应不变法必然会产生严重的频率混叠失真。 所以应选用双线性变换法。 将Ha(jΩ)中的jΩ用s代替， 可得到RC滤波网络的系统函数： 用双线性变换法设计公式, 可得 H(z)的结构图如题8解图所示。 题8解图 　　由图可见， 在模拟域由一个R和一个C组成的RC滤波网络， 用双线性变换法转换成数字滤波器后， 用两个乘法器、 两个加法器和一个单位延迟器实现其数字滤波功能。 也可用软件实现该数字滤波功能。 由滤波器差分方程编写程序较容易。 为此， 由H(z)求出差分方程。 编程序实现差分方程中的计算， 即可实现对输入信号序列x(n)的高通滤波。  　　9． 设计低通数字滤波器， 要求通带内频率低于0.2π rad时， 容许幅度误差在1 dB之内； 频率在0.3π到π之间的阻带衰减大于10 dB。 试采用巴特沃斯型模拟滤波器进行设计， 用脉冲响应不变法进行转换， 采样间隔T=1 ms。 　　解: 本题要求用巴特沃斯型模拟滤波