子空间迭代法求解矩阵的前n阶固有振型和固有频率

博大广阔

%子空间迭代法
%其实质是不断改善李籽法中的假设阵型质量，使得李兹法的假设阵型所张成的子空间不断接近前R阶主阵型
function [R,XX]=zmatrix_interation(K,M,number,accuracy)
if nargin<3;warndlg('wrong inpiut');end
if nargin<4; accuracy=0.001;end
D=inv(K)*M;  

%%%
那个
%%%

erro=10^4;
  while erro>=accuracy  
  X1=D*Xold
  a=max(X1);
  for k=1:1:size(X1,2)   
  X1(:,k)=X1(:,k)./a(k);
  end
K1=X1.'*K*X1;
M1=X1.'*M*X1;
[V,rr]=eig(inv(K1)*M1);
Xnew=Xold*V;
a=Xnew-Xold;
erro=norm(a(:,1),2);
Xold=Xnew;
   end
R=rr;
XX=Xold;
end

小胡

您好！我想知道式中的number 和accuracy 指什么？

mxlzhenzhu

我有两个问题：
1，刚度矩阵是奇异的怎么办？
2，求逆运算会不会很不好啊？