from sippy import system_identification as sysid
openLoopInput = [0, 0, 0, 0, 0, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100,]
m = sysid(openLoopOutput, openLoopInput, 'FIR', FIR_orders=[10,1], tsample=1)

<TransferFunction>: sys[0]
Inputs (1): ['u[0]']
Outputs (1): ['y[0]']

0.1744 z^9 + 0.1744 z^8 + 0.1744 z^7 + 0.1744 z^6 + 0.1744 z^5 + 0.011 z^4 + 0.007 z^3 + 0.003 z^2
--------------------------------------------------------------------------------------------------
                                               z^11

dt = 1.0

simInput = [0, 0, 0, 0, 0, 120, 120, 120, 120, 120, 120, 120,]
dataLen = len(simInput)
simout, _, _ = control.matlab.lsim(m.G, simInput, np.linspace(0, dataLen-1, dataLen))

control/timeresp.py:1083: UserWarning: return_x specified for a transfer function system. Internal conversion to state space used; results may meaningless.