032 Matlab 각도에 따른 포물선 운동 (공기 마찰력)

각도에 따른 포물선 운동 (공기 마찰력)

공기의 마찰력을 속도의 0.1 % 만큼 감소할때를 고려

hold off;

g=9.81;
v=input('초기 속도를 입력하시오.');
theta=input('초기 각도를 입력하시오.');

dt=input('시간 간격을 입력하시오');
x=0;
y=0;
dxdt=v*cos(theta*(pi/180));
dydt=v*sin(theta*(pi/180));

t1=(2*v*sin(theta*(pi/180)))/(g+0.1*v);
for k=0:1:t1/dt;
   
    plot(x,y,'-ro')
    time= dt*k;
   
    dxdtdt=-0.1*dxdt;
    dydtdt=-0.1*dydt;
   
    dxdt=dxdt+dxdtdt*dt;
    dydt=dydt+-g*dt+dydtdt*dt;
   
    x=x+(dxdt*dt);
    y=y+(dydt*dt);
    hold on
  
    title([num2str(time),'초'])
    pause(0.1)
end