Matlab中zplane函数用法

1
2
3
4
5
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。

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

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

1
2
3
4
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弧度/样本。

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

1
zplane(b,a)

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
[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

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