scipy.signal.

lsim

scipy.signal.lsim(system, U, T, X0=None, interp=True)

模拟连续时间线性系统输出。

参数:

systemLTI 类的实例或描述系统元组。

以下给出元组中的元素数量和解释

• 1: (lti 实例)

• 2: (num, den)

• 3: (zeros, poles, gain)

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

U类数组

描述在每个时间 T 处输入的输入数组（给定时间之间假设有内插）。如果有多个输入，则秩 2 数组的每一列表示一个输入。如果 U = 0 或 None，则使用零输入。

T类数组

定义输入和期望输出的时间步长。必须非负、递增且等间距。

X0类数组，可选

状态变量的初始条件（默认为零）。

interpbool，可选

是否使用线性（True，默认值）或零级保持（False）插值作为输入数组。

返回:

T1D ndarray

输出的时间值。

yout1D ndarray

系统响应。

xoutndarray

状态向量的演化时间。

注释

如果代码 system传入的是 (num, den)，那么分子和分母的系数应当按照指数降序指定（例如，s^2 + 3s + 5 将表示为 [1, 3, 5]）。

示例

可在浏览器中尝试！

我们将使用 lsim 在应用于信号的模拟贝塞尔滤波器中进行模拟。

>>> import numpy as np
>>> from scipy.signal import bessel, lsim
>>> import matplotlib.pyplot as plt

创建一个截止频率为 12 Hz 的低通贝塞尔滤波器。

>>> b, a = bessel(N=5, Wn=2*np.pi*12, btype='lowpass', analog=True)

生成应用滤波器的数据。

>>> t = np.linspace(0, 1.25, 500, endpoint=False)

输入信号为三个正弦曲线的总和，频率分别为 4 Hz、40 Hz 和 80 Hz。滤波器应当消除 40 Hz 和 80 Hz 成分，仅仅保留 4 Hz 信号。

>>> u = (np.cos(2*np.pi*4*t) + 0.6*np.sin(2*np.pi*40*t) +
...      0.5*np.cos(2*np.pi*80*t))

使用 lsim 模拟滤波器。

>>> tout, yout, xout = lsim((b, a), U=u, T=t)

绘制结果。

>>> plt.plot(t, u, 'r', alpha=0.5, linewidth=1, label='input')
>>> plt.plot(tout, yout, 'k', linewidth=1.5, label='output')
>>> plt.legend(loc='best', shadow=True, framealpha=1)
>>> plt.grid(alpha=0.3)
>>> plt.xlabel('t')
>>> plt.show()

在第二个示例中，我们模拟一个常数输入 y'' = u 的双重积分器 u = 1。我们将使用积分器的状态空间表示。

>>> from scipy.signal import lti
>>> A = np.array([[0.0, 1.0], [0.0, 0.0]])
>>> B = np.array([[0.0], [1.0]])
>>> C = np.array([[1.0, 0.0]])
>>> D = 0.0
>>> system = lti(A, B, C, D)

t 和 u 定义要模拟系统的的时间和输入信号。

>>> t = np.linspace(0, 5, num=50)
>>> u = np.ones_like(t)

计算模拟，然后绘制 y。不出所料，该绘制显示曲线 y = 0.5*t**2。

>>> tout, y, x = lsim(system, u, t)
>>> plt.plot(t, y)
>>> plt.grid(alpha=0.3)
>>> plt.xlabel('t')
>>> plt.show()

返回在新标签页中打开