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

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

用法:

scipy.signal.freqresp(system, w=None, n=10000)#

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

參數 ：：

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

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

• 1 (instance of lti)

• 2 (num, den)

• 3 (zeros, poles, gain)

• 4 (A, B, C, D)

w： 數組，可選

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

n： 整數，可選

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

返回 ：：

w： 一維ndarray

頻率陣列 [rad/s]

H： 一維ndarray

複數幅度值數組

注意：

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

例子：

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

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

構造傳遞函數\(H(s) = \frac{5}{(s-1)^3}\)：

>>> s1 = signal.ZerosPolesGain([], [1, 1, 1], [5])

>>> w, H = signal.freqresp(s1)

>>> plt.figure()
>>> plt.plot(H.real, H.imag, "b")
>>> plt.plot(H.real, -H.imag, "r")
>>> plt.show()