Python SciPy signal.dfreqresp用法及代碼示例

本文簡要介紹 python 語言中 scipy.signal.dfreqresp 的用法。

用法:

scipy.signal.dfreqresp(system, w=None, n=10000, whole=False)#

計算discrete-time 係統的頻率響應。

參數 ：：

system： dlti 類的實例或說明係統的元組。

下麵給出元組中元素的數量和解釋：

• 1 (instance of dlti)

• 2 (numerator, denominator, dt)

• 3 (zeros, poles, gain, dt)

• 4 (A, B, C, D, dt)

w： 數組，可選

頻率數組(以弧度/樣本為單位)。為該數組中的每個值計算幅度和相位數據。如果沒有給出合理的設置，將被計算。

n： 整數，可選

如果未給出 w，則要計算的頻點數。 n 個頻率以對數方式隔開，間隔選擇包括係統極點和零點的影響。

whole： 布爾型，可選

通常，如果沒有給出‘w’，則從 0 到奈奎斯特頻率計算頻率，pi 弧度/樣本(upper-half of unit-circle)。如果整體為真，則計算從 0 到 2*pi 弧度/樣本的頻率。

返回 ：：

w： 一維ndarray

頻率陣列 [弧度/樣本]

H： 一維ndarray

複數幅度值數組

注意：

如果 (num, den) 為 system 傳入，則分子和分母的係數應按 index 降序指定(例如 z^2 + 3z + 5 將表示為 [1, 3, 5] )。

例子：

生成傳遞函數的奈奎斯特圖

>>> from scipy import signal
>>> import matplotlib.pyplot as plt

構造采樣時間為 0.05 秒的傳遞函數\(H(z) = \frac{1}{z^2 + 2z + 3}\)：

>>> 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()