- 145 名前：ご冗談でしょう？名無しさん mailto:sage [2012/01/23(月) 23:30:13.17 ID:???]
- n=length(W);
W=W*1000;
K=K*10^5;
F=F*1000;
dt=T(1,2)-T(1,1);
tt=(length(T)-1)*dt;
m1=W(1,1)/9.8;
m2=W(1,2)/9.8;
m=[m1 0;0 m2];
k1=K(1,1);
k2=K(1,2);
k=[k1+k2 -k2;-k2 k2];
%固有値の計算
ww=m\k;
[w,e]=eig(ww);
www=sqrt(e);
disp(w)
disp(www)
dooo=k\(P*F);
%x''=m-1*(F-kx)により以下のようにx''を求める。
ao=m\(P*F-k*dooo);
f=P*F;
%ここでは2つあるうちの2つ目の固有円振動数を使用
aa=f*cos(www(2,2)*T);
