用matlab求解机械振动的固有频率

如果M质量矩阵和K刚度矩阵元素很多，矩阵很大，如何快速求解固有频率w呢
（K-wM）q=0
用eig命令好像求解很慢
大家会做吗，请教一下

试试这个！

代码：

clear all
clc

tic
% Sylvester方程 AX + XB = -C 
% 套用上述方程
A=zeros(1000,1000);
K=rand(1000,1000)*100;
M=rand(1000,1000)*100;
%结果
w=lyap(A,M,-K);
%验证
norm(A*w + w*M + -K)
toc

有用请采纳！