作业帮 > 数学 > 作业

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

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/09 09:41:03
Matlab求解矩阵方程det(K-w^2*M)=0
如题,其中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了...是哪里出了问题呢?
没有问题的,带回去A的行列式是0的,不过是由于是数值计算,结果和0是差一点点,但基本上都是0,这是数值计算带来的误差.
再问: 谢谢你的回答,难道我代入过程错了? w=8.8465315330125002672255483959155 %得到的一个解 det([K-w^2*M]) 结果是ans = 2.4641e+009 这是为什么呢?是不是我的w没有真正代入呢?
再答: 你带入的方法不对,w是一个数,det([K-w.^2*M])