U为主动减震施加的力。

 

通过python control库，可以得到在单位脉冲激励下X1的响应曲线

import control
import matplotlib.pyplot as plt


M1 = 2500
M2 = 320
K1 = 80000
K2 = 500000
b1 = 350
b2 = 15020

s = control.TransferFunction.s

H = (b1*b2*s*s + (b1*K2+b2*K1)*s + K1*K2) / \
    ((M1*s*s + b1*s + K1) * (M2*s*s + (b1+b2)*s + (K1+K2)) - (b1*s+K1)*(b1*s+K1))

t, y = control.impulse_response(H)

plt.plot(t, y)
plt.show()