control的官方网站

函数示例

强迫响应forced_response

import numpy as np
import os
import sys
import control as ctrl
import matplotlib.pyplot as pltdef lim_x(x, lim=0):res = 0if x > lim:res = 1else:res = 0return resif __name__ == "__main__":print(__file__)fc = 20Ts = 50*10e-6T = 1/np.pi/2/fcomega = 2*np.pi*5/2t = np.arange(0, 2, Ts)y = np.sin(omega*t)ys = [lim_x(k) for k in y]plt.close("all")plt.figure("响应函数")lowf = ctrl.tf([1], [T, 1])lowf_z = ctrl.c2d(lowf, Ts)tout, yout = ctrl.forced_response(lowf_z, T=t, U=ys, X0=0)plt.plot(tout, yout)plt.plot(t, ys)plt.legend(["input", "output"])plt.xlabel("time/s")plt.show()# print(lowf, lowf_z)