Scilab 中的零极点图

Question

在尝试为离散时间全通系统生成极点零图时，该系统在 3/4 处有一个极点，在 4/3 处有一个零，我使用的 scilab 代码不起作用。

我用了他的代码

z=poly(0,'z')
n = [z-4/3];          // Defining Nr polynomial
d = [z-3/4];          // Defining Dr polynomial
h=syslin('d',n./d);  // Defining linear system
plzr(h)

在上面的示例中，极点和零点均位于实轴上，即 r = 3/4 和 4/3，且 theta = 0 度。 假设我们将极点和零分别放置在 3/4 e^j pi/3, 4/3 e^j pi/3 处，那么

z=poly(0,'z')
n = [z-(4/3)*exp(%i*%pi/3)];          // Defining Nr polynomial
d = [z-(3/4)z*exp(%i*%pi/3)];          // Defining Dr polynomial
h=syslin('d',n./d);  // Defining linear system
plzr(h)

它不起作用！有谁可以帮忙吗？？

Answer 1

plzr

仅适用于真实系统，因此您必须添加复杂的共轭极点和零点：

z=poly(0,'z')
n = [z-(4/3)*exp(%i*%pi/3)];          // Defining Nr polynomial
d = [z-(3/4)*exp(%i*%pi/3)];          // Defining Dr polynomial
h=syslin('d',real(n'*n)/real(d'*d));  // Defining linear system
plzr(h)