scipy.signal.dfreqresp¶
- scipy.signal.dfreqresp(system, w=None, n=10000, whole=False)[源代码]¶
计算离散时间系统的频率响应。
- 参数
- 系统 :一个
dlti
类或描述系统的元组。对象的一个实例 下面给出了元组中的元素数量和解释:
1(实例为
dlti
)2(分子,分母,DT)
3(零点、极点、增益、DT)
4(A、B、C、D、DT)
- wARRAY_LIKE,可选
频率数组(以弧度/样本为单位)。此数组中的每个值都会计算幅值和相位数据。如果没有给出,将会计算出一个合理的设置。
- n整型,可选
在以下情况下要计算的频率点数 w 是不会被给予的。这个 n 频率以对数间隔分布在选定的区间内,以包括系统的极点和零点的影响。
- whole布尔值,可选
通常,如果没有给出‘w’,频率的计算范围是从0到奈奎斯特频率,π弧度/样本(单位圆的上半部分)。如果 whole 为True,则计算0到2*pi弧度/样本的频率。
- 系统 :一个
- 退货
- w一维ndarray
频率阵列 [radians/sample]
- H一维ndarray
复数量级值数组
注意事项
如果传入(num,den)作为
system
,分子和分母的系数应按指数降序指定(例如z^2 + 3z + 5
将表示为[1, 3, 5]
)。0.18.0 新版功能.
示例
生成传递函数的奈奎斯特图
>>> from scipy import signal >>> import matplotlib.pyplot as plt
构造传递函数 \(H(z) = \frac{{1}}{{z^2 + 2z + 3}}\) 采样时间为0.05秒时:
>>> sys = signal.TransferFunction([1], [1, 2, 3], dt=0.05)
>>> w, H = signal.dfreqresp(sys)
>>> plt.figure() >>> plt.plot(H.real, H.imag, "b") >>> plt.plot(H.real, -H.imag, "r") >>> plt.show()