import numpy as np
import matplotlib.pyplot as plt
import scipy.signal as sig
# Parameters defining the system
c = 4 # Damping constant
k = 2 # Stiffness of the spring
m = 20 # Mass
F = 5 # Force
Ft = np.ones(610)*F
# Simulation Parameters
t=np.linspace(1,100,1000)
# System matrices
A = np.array(
0, 1,0,0], [-k/m, -c/m,0,c/m],[0,0,0,1],[0,c/m,k/m,c/m)
B =np.array(
0,0],[1/m,0],[0,0],[0,1/m)
C =np.array(
1, 0,0,0],[0,0,1,0)
D=np.array(
0],[0)
sys = sig.StateSpace(A, B, C, D)
print(sys)
t, y, x = sig.lsim(sys, Ft, t)
plt.plot(t, x)
#plt.plot(t, y)
plt.title('Simulation of Mass-Spring-Damper System')
plt.xlabel('t')
plt.ylabel('x(t)')
plt.grid()
plt.show()
this is the error : ValueError: The input arrays have incompatible shapes.