Tools Profile Research Publications Activities Tools etc.  Matlab用 便利ソース  固有振動数等を計算 （ファンクションファイル）  function [wn,v]=eigsort(k,m); [vtem,d]=eig(k,m); Omega=sqrt(diag(d)); [wn,isort]=sort(Omega); il=length(wn);     for i=1:il     v(:,i)=vtem(:,isort(i));     end  Matlab用 便利ソース  状態方程式表現に変換 （ファンクションファイル）  function [A,B]=de2ss(M,C,K,H) n=length(M); [Hc Hr]=size(H); A=[zeros(n,n),eye(n);-inv(M)*K,-inv(M)*C]; B=[zeros(n,Hr);inv(M)*H];