control的官方网站

函数示例

强迫响应forced_response

import numpy as np
import os
import sys
import control as ctrl
import matplotlib.pyplot as plt


def lim_x(x, lim=0):
    res = 0
    if x > lim:
        res = 1
    else:
        res = 0
    return res


if __name__ == "__main__":
    print(__file__)

    fc = 20
    Ts = 50*10e-6
    T = 1/np.pi/2/fc

    omega = 2*np.pi*5/2
    t = 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)