我正在尝试使用 python 绘制离散传递函数 H(z) = (0.5 - 0.50z)/(z-1.0) 但我不断收到绘图错误。我试图在 scipy.signal 中使用 dstep 函数。连续绘图工作正常,但我在行 plt.plot(t,y,'k-',label=name)
中遇到错误.我将它用于连续传递函数并且工作正常。我认为其中一个论点一定是不正确的
from scipy import signal
import matplotlib.pyplot as plt
name = 'Discrete Transfer Function'
num = [-0.50, 0.50]
den = [1.0, -1.0]
sys3 = signal.TransferFunction(num, den, dt = 0.001)
t,y = signal.dstep(sys3)
plt.figure(3)
plt.plot(t,y,'k-',label=name)
plt.legend(loc='best')
plt.show()
最佳答案
y
是一个numpy数组的元组,所以它需要被索引为y[0]
并用 ravel()
压平(似乎 matplotlib 可以在不展平的情况下计算出来):
plt.plot(t, y[0], 'k-', label=name)
关于python - 使用 scipy.signal.dstep 绘制离散传递函数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/67248802/