Tools



ProfileResearchPublicationsActivitiesToolsetc.
 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];