- 407 名前：ご冗談でしょう？名無しさん [2012/01/28(土) 22:06:27.15 ID:Qjvpgfb+]
- >>404
返信ありがとうございます。したことというと、
最初の初期条件の時にT=0としているためプログラムを以下のように
書き換えました。t=dt:dt:ttにしました。0から始めると、初期条件を
書き換えてしまうので。
W=W*1000;
F=F*1000;
dt=T(1,2)-T(1,1);
tt=(length(T)-1)*dt;
m1=W/9.8;
m=m1;
k=K;
k=k*10^5;
ww=k/m;
www=sqrt(ww);
dooo=(P*F)/k;
aoo=(P*F-k*dooo)/m;
A=m*6/dt^2+k;
B=m*6/dt^2;
C=m*6/dt;
D=2*m;
xx1=[];
xx2=[];
xx3=[];
aa=[];
f=P*F;
aa=f*cos(www*T);
xx1(1,1)=dooo;
xx2(1,1)=vo;
xx3(1,1)=aoo;
