Matlab中zplane函数用法

目录

语法

描述

示例

椭圆低通滤波器的极点和零点

传递函数的零点和极点

---

zplane离散时间系统的零极点图

语法

zplane(z,p)
zplane(b,a)
[hz,hp,ht] = zplane(___)
zplane(d)
[vz,vp,vk] = zplane(d)

描述

zplane(z,p)在当前图形窗口中 绘制列向量中指定的零点z和列向量p中指定的极点。符号'o'代表零，符号 'x'代表极。该图包括单位圆作为参考。如果z和p是矩阵，然后 zplane绘制中的列中的极点和零点 z，并p以不同的颜色。

zplane(b,a)，其中b和a是行向量，首先用于roots查找由分子系数b和分母系数 表示的传递函数的零点和极点 a。

[hz,hp,ht] = zplane(___) 将控点向量返回零线hz和极线hp。ht是轴/单位圆线和文本对象的句柄的向量，当存在多个零或极点时，这些对象存在。

zplane(d)找到由数字滤波器表示的传递函数的零点和极点 d。用于根据频率响应规格designfilt生成d。零极点图显示在FVToll中。

[vz,vp,vk] = zplane(d) 返回与数字滤波器相对应的零（向量vz），极点（向量 vp）和增益（标量vk）d。

示例

椭圆低通滤波器的极点和零点

        对于以1000Hz采样的数据，绘制一个四阶椭圆低通数字滤波器的零点和零点，截止频率为200Hz，通带中的波纹为3dB，阻带中的衰减为30dB。

[z,p,k] = ellip(4,3,30,200/500);
zplane(z,p)
grid
title('4th-Order Elliptic Lowpass Digital Filter')

        使用创建相同的过滤器designfilt。使用zplane绘制极点和零点。注意zplanecall的这种语法fvtool。

d = designfilt('lowpassiir','FilterOrder',4,'PassbandFrequency',200, ...
               'PassbandRipple',3,'StopbandAttenuation',30, ...
               'DesignMethod','ellip','SampleRate',1000);
zplane(d)

传递函数的零点和极点

        设计一个阻带衰减为20dB的8阶Chebyshev II型带通滤波器。将阻带边沿频率指定为π/8 rad/sample和 5π/8弧度/样本。

[b,a] = cheby2(8/2,20,[1 5]/8);

        使用zplane绘制传递函数的极点和零点。

zplane(b,a)

        可视化滤波器的零相响应。覆盖单位圆，极点和零位。

[hw,fw] = zerophase(b,a,1024,"whole");

z = roots(b);
p = roots(a);

plot3(cos(fw),sin(fw),hw)
hold on
plot3(cos(fw),sin(fw),zeros(size(fw)),'--')
plot3(real(z),imag(z),zeros(size(z)),'o')
plot3(real(p),imag(p),zeros(size(p)),'x')
hold off
xlabel("Real")
ylabel("Imaginary")
view(35,40)
grid

传递函数系数，指定为行向量。传递函数根据z^-1定义 ：

数字滤波器，指定为digitalFilter对象。使用designfilt生成数字 基于频率响应规格的滤波器。数字滤波器的零点，极点和增益d作为列向量和标量返回。

提示

• 可以覆盖的自动缩放zplane 使用

axis([xmin xmax ymin ymax])

        调用zplane之后。 当一个或多个零或极点的幅度太大而其他零点或极点围绕原点紧密组合且难以区分时，此缩放很有用。