481
2016-03-20 01:54:50
1
clc;
clear all;
data = 0.5:0.5:10;
for i = 1:1:20
tf = 20;
t=0;
v=0;
h=data(i);
while(1)
if (t+data(i))>tf
h=tf-t;
end
dvdt=9.81-(0.3/65)*v*abs(v);
v=v+dvdt*h;
t=t+h;
if t>=tf
break
end
end
data(i) = v;
end
plot(data);
대충만들어서 인덱스가 다르긴한데 뭐 디버깅용으로 보는데는 문제없을듯