Matlab求解矩阵方程det（K-w^2*M）=0

Question

如题，其中K,M已知。clear
M=[46.512 0 0 0 0;0 92.55 0 0 0; 0 0 92.55 0 0; 0 0 0 92.55 0; 0 0 0 0 92.55]
k=74070
K=k*[1 -1 0 0 0; -1 2 -1 0 0; 0 -1 2 -1 0; 0 0 -1 2 -1; 0 0 0 -1 2]
syms w
A=K-w^2*M
[w]=solve(det(A));w

之后用得到的w再带回去，A的行列式却不等于0了...求指点，是哪里出了问题呢？

Answer 1

没有问题的，带回去A的行列式是0的，不过是由于是数值计算，结果和0是差一点点，但基本上都是0，这是数值计算带来的误差。

追问

谢谢你的回答，难道我代入过程错了？
w=8.8465315330125002672255483959155 %得到的一个解
det([K-w^2*M])

结果是ans = 2.4641e+009
这是为什么呢？是不是我的w没有真正代入呢？

追答

你带入的方法不对，w是一个数，det([K-w.^2*M])

Answer 2

没有问题的，带回去A的行列式是0的，不过是由于是数值计算，结果和0是差一点点，但基本上都是0，这是数值计算带来的误差。 但是0AW 就算了。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。