python - 热门创建绘图以可视化结构的模式

Question

我有一个用电线连接的具有 4 个相等质量的摆锤。它们是等距的。我找到了系统的特征值和特征向量。我可以通过哪种方式可视化系统的模式号 1、2 和 3？

这是我的代码：

import numpy as np
from scipy import linalg
from numpy.linalg import eig
import matplotlib.pyplot as plt
# parameters 

L = 0.5 # length
Nx = 4 # number of masses
M = 0.1 # total mass
g= 9.81 # gravity constant
I0 = 0.001 # impulse

Nev = 3 # number of eigenvalues

fs=100 #sampling frequency
t_max=5 #maximum time for the calculation of the response

#Useful variables

dx=L/Nx #x spacing
m=M/Nx #mass of each node

x=[dx,dx,dx,dx]

#matrices
#mass matrix
m_vec=np.diag(np.ones(len(x)))
M=m*m_vec
M1=np.linalg.inv(M)

#Stiffness Matrix
T1=np.array([4,3,2,1])
T2=np.array([3,2,1,0])
T=T1*m*g
Tplus=T2*m*g

K=np.diag(T)
Kplus=np.diag(Tplus)
D1=np.diag(-T, k=1)
K1=D1[1:5,1:5]
D2=np.diag(-T, k=-1)
K2=D2[1:5,1:5]
Ktot=(K+K1+K2+Kplus)/dx

#ProductMatrices

P=np.dot(M1,Ktot)
(w,v)=eig(P)

谢谢