ch2_3系统传递函数绪论.ppt

模拟频率f 和数字频率ω 数字滤波器的形状设计可以不依赖采样频率，但采样频率将影响滤波器的带宽。 “基础”P193 例7.18-7.20 这三个例子都说明了： 采样频率对滤波器性能（带宽）的影响。 例2.7.4 已知H(z)=1-z-N，试定性画出系统的幅频特性。 解： 零点： H(z)的极点为z=0，这是一个N阶极点，它不影响系统的频响。零点有N个，由分子多项式的根决定 N个零点等间隔分布在单位圆上，设N=8，极零点分布如图2.7.5所示。当ω从零变化到2π时，每遇到一个零点，幅度为零，在两个零点的中间幅度最大，形成峰值。幅度谷值点频率为：ωk=(2π/N)k,k=0,1,2,…(N-1)。一般将具有如图2.7.5所示的幅度特性的滤波器称为梳状滤波器。 图2.7.5 梳状滤波器的极零点分布及幅度特性 例2.7.5 利用几何法分析矩形序列的幅频特性。 解： 零点： 极点： 设N=8，z=1处的极点零点相互抵消。这样极零点分布及其幅频特性如图2.7.6所示。 阶零点 图2.7.6 N=8矩形序列极零点分布及幅度特性 低通滤波器 由于零点并不能完全组合成复数共轭对，该滤波器系数将是复数。这对滤波器实现不利。 对带通滤波器的设计思路：在单位圆上设定零点的方法来控制滤波器频率响应。若将零极点抵消位置移至某一频率，则可获相应BPF。 实系数带通滤波器 带通滤波器 MATLAB IMPLEMENTATION EXAMPLE1 Given a causal system a. Find H(z) and sketch its pole-zero plot. b. Plot c. Determine the impulse response h(n). y(n)=0.9y(n-1)+x(n) Solution: a. From We will use Matlab to illustrate the use of the zplane function: b=1; a=[1,-0.9]; zplane(b,a); y(n)=0.9y(n-1)+x(n) MATLAB IMPLEMENTATION b. Matlab provides a function called freqz, the simplest form of this function is invoked by: [H, w] = freqz ( b, a, N) which returns the N-point frequency vector w and the N-point complex frequency response vector H of the system, given its numerator and denominator coefficients in vectors b and a. The frequency response is evaluated at N points equally spaced around the upper half of the unit circle. (0 to π) Let us take 100 points along the upper half of the unit circle: % z-平面函数: % a) b) b = 1; a = [1, -0.9]; zplane(b,a); title(极点-零点图); text(0.85,-0.1,0.9);text(0.01,-0.1,0);pause, % b) [H,w]=freqz(b,a,100); magH = abs(H);phaH = angle(H); subplot(2,1,1),plot(w/pi,magH);grid xlabel();ylabel(幅度); title(幅度响应) subplot(2,1,2);plot(w/pi,phaH/pi);grid xlabel(频率 (单位:pi));ylabel(相位 (单位:pi)); title(相位响应), MATLAB IMPLEMENTATION c. Matlab provides a function called filter, which is invoked b