mkeig.m

来自「利用matalab编制的转子动力学计算软件。其中模型利用有限元建模」· M 代码 · 共 20 行

function [V,D]=mkeig(K,M)
%[Vectors, EigenValues]=mkeig(K,M)
%
% [V, s]=mkeig(K,M)
%	compute the mass-normalised eigenvectors and eigenvalues of s*M*V=M*V
%
% s=mkeig(K,M)
%	compute the  eigenvalues of s*M*V=M*V
%
% by : I. Bucher 1.Nov.1993
%

 [V,D]=eig(K,M);
 [D j]=sort(diag(D));
 V=V(:,j);
 mr=diag(V.'*M*V);
 n=length(mr);
 V=V*diag(ones(n,1)./sqrt(mr(:)));
%