子空间迭代法求解矩阵的前n阶固有振型和固有频率
%子空间迭代法%其实质是不断改善李籽法中的假设阵型质量，使得李兹法的假设阵型所张成的子空间不断接近前R阶主阵型
function =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;
=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 指什么？
我有两个问题：
1，刚度矩阵是奇异的怎么办？
2，求逆运算会不会很不好啊？
页:
[1]