问题在这里:
r = I + k + np.square(k) * ((1 -c)/(s**2))
np.square(k)对矩阵的每个元素求平方。您想要
np.matmul(k,k)或
k @ k矩阵乘以本身。
我还将实施该
s=0答案的注释中提到的附带情况(尤其是),否则在很多情况下您都会出错。

问题在这里:
r = I + k + np.square(k) * ((1 -c)/(s**2))
np.square(k)对矩阵的每个元素求平方。您想要
np.matmul(k,k)或
k @ k矩阵乘以本身。
我还将实施该
s=0答案的注释中提到的附带情况(尤其是),否则在很多情况下您都会出错。