(新)第2章.pptVIP

2.3.2　周期序列的傅里叶变换表示式 （2.6.7） （2.6.8） (2.6.9) 系统的频响特性由(2.6.8)式和(2.6.9)式确定。当频率ω从0变化到2π时，这些向量的终点B沿单位圆逆时针旋转一周，按照(2.6.8)式和(2.6.9)式，分别估算出系统的幅频特性和相频特性。例如图2.6.2表示了具有一个零点和两个极点的频率特性。 图2.6.2　频响的几何表示法 　　按照(2.6.8)式，知道零极点的分布后，可以很容易地确定零极点位置对系统特性的影响。当B点转到极点附近时，极点相量长度最短，因而幅度特性可能出现峰值，且极点愈靠近单位圆，极点相量长度愈短，峰值愈高愈尖锐。如果极点在单位圆上，则幅度特性为∞，系统不稳定。对于零点，情况相反，当B点转到零点附近时，零点相量长度变短，幅度特性将出现谷值，零点愈靠近单位圆，谷值愈接近零。当零点处在单位圆上时，谷值为零。总结以上结论：极点位置主要影响频响的峰值位置及尖锐程度，零点位置主要影响频响的谷点位置及形状。 　　这种通过零极点位置分布分析系统频响的几何方法为我们提供了一个直观的概念，对于分析和设计系统是十分有用的。基于这种概念，可以用零极点累试法设计简单滤波器。 　　下面介绍用MATLAB计算零、极点及频率响应曲线。首先介绍MATLAB工具箱中两个函数zplane和freqz的功能和调用格式。 　　zplane　绘制H(z)的零、极点图。 　　zplane(z, p)　绘制出列向量z中的零点（以符号“○”表示）和列向量p中的极点（以符号“×”表示），同时画出参考单位圆，并在多阶零点和极点的右上角标出其阶数。如果z和p为矩阵，则zplane以不同的颜色分别绘出z和p各列中的零点和极点。 　　zplane(B, A)绘制出系统函数H(z)的零极点图。其中B和A为系统函数H(z) = B(z)/A(z)的分子和分母多项式系数向量。假设系统函数H(z)用下式表示: 则 B=B(1)　B(2)　B(3)　　B(M+1)］，A=［A(1)　A(2)　A(3)　　　　A(N+1)］ 　　freqz　计算数字滤波器H(z)的频率响应。 　　H=freqz(B, A, w)　计算由向量w指定的数字频率点上数字滤波器H(z)的频率响应H(ejw)，结果存于H向量中。B和A仍为H(z)的分子和分母多项式系数向量（同上）。 ［H, w］= freqz(B, A, M)　计算出M个频率点上的频率响应，存放在H向量中，M个频率存放在向量w中。freqz函数自动将这M个频率点均匀设置在频率范围［0, π］上。 　　［H, w］ = freqz(B, A, M, ′whole′)　自动将M个频率点均匀设置在频率范围［0, 2π］上。 　　当然，还可以由频率响应向量H得到各采样频点上的幅频响应函数和相频响应函数；再调用plot绘制其曲线图。 |H(ejω)|=abs(H) φ(ω)=angle(H) 式中，abs函数的功能是对复数求模，对实数求绝对值；angle函数的功能是求复数的相角。 　　freqz(B, A)　自动选取512个频率点计算。不带输出向量的freqz函数将自动绘出固定格式的幅频响应和相频响应曲线。所谓固定格式, 是指频率范围为［0, π］，频率和相位是线性坐标，幅频响应为对数坐标。　其他几种调用格式可用命令help查阅。 　　【例2.6.2】　已知H(z)=z－1，分析其频率特性。 　　解　由H(z)=z－1，可知极点为z=0，幅频特性|H(ejω)|=1, 相频特性φ(ω)=－ω，频响特性如图 2.6.3所示。用几何方法也容易确定，当ω=0转到ω=2π时，极点向量的长度始终为1。由该例可以得到结论：处于原点处的零点或极点，由于零点向量长度或者极点向量长度始终为1，因此原点处的零极点不影响系统的幅频响应特性, 但对相频特性有贡献。 图2.6.3　H(z)=z－1的频响特性 W(z)的收敛域为|a||z|≤∞；被积函数υ平面上的收敛域为max(|a|, 0)|υ|min(|a－1|, |z|)，υ平面上极点：a、a－1和z，c内极点：z=a。 令 　　10． 帕斯维尔（Parseval）定理 　　设X(z)=ZT［x(n)］　　Rx－|z|Rx+ 　　　Y(z)=ZT［x(n)］　　Rx－|z|Rx+ 　　　Rx－Ry－1, Rx+Ry+1 　　那么 （2.5.27） υ平面上，c所在的收敛域为 　　利用复卷积定理可以证明上面的重要的帕斯维尔定理。　证明　令　　　　　w(n)=x(n)y*(n) 按照（2.5.24）式得到： 按照（2.5.25）式，Rx－Ry－|z|Rx+Ry+; 按照假设，z=1在收敛域中，将z=1代入W(z)中，则