1、zplane函数
zplane(z, p) 绘制出列向量z中的零点(以符号"○" 表示)和列向量p中的极点(以符号"×"表示)，同时画出参考单位圆，并在多阶零点和极点的右上角标出其阶数。如果z和p为矩阵，则zplane以不同的颜色分别绘出z和p各列中的零点和极点。
zplane(B, A) 绘制出系统函数H(z)的零极点图。其中B和A为系统函数H(z) = B(z)/A(z)的分子和分母多项式系数向量。 zplane(B, A) 输入的是传递函数模型，函数首先调用root函数以求出它们的零极点。
2、freqz函数
MATLAB提供了专门用于求离散系统频响特性的函数freqz(),调用freqz()的格式有一下两种：
（1）[H,w] = freqz(B,A,N)
B和A分别为离散系统的系统函数的分子、分母多项式的系数向量，N为正整数，返回向量H则包含了离散系统频响在0-PI范围内N个频率等分点的值，向量w则包含0-PI范围内N个频率等分点。调用中若N默认，默认值为512。
（2）[H,w] = freqz(B,A,N,'whole')
该调用格式将计算离散系统在0—PI范内的N个频率等分店的频率响应的值。
因此，可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，即可绘制出系统在 或 范围内的频响曲线。
例子：绘制出系统函数H(z)=(z-0.5)/z的频响曲线。
MATLAB命令如下：
参考文献：http://blog.csdn.net/silenceee/article/details/1361258
（3）[h,f]=freqz(b,1,2048,75000,'whole');
B:系统函数分子系数；
A:系统函数分母系数；
2048:整个单位圆上分为多少个点用以计算频率响应；
'whole':在整个单位圆上计算频率响应,对应的角频率范围是[0,2pi],不加whole只在上半圆计算,对应角频率范围是[0,pi]。
参考文献：http://www.ilovematlab.cn/thread-40158-1-1.html
原文：https://blog.csdn.net/yundanfengqing_nuc/article/details/45366175