034 Matlab 포물선 운동

포물선 운동

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));
k=0;
t1=(2*v*sin(theta*(pi/180)))/(g+0.1*v);
t2=(v*sin(theta*(pi/180)))/(g+0.1*v);
plot(0,0,'-ro')
while y>=0
   
    plot(0,0)
    hold on
   
    k=k+1;   
   
    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);
    plot(x,y,'-ro')  
    title([num2str(time),'초'])
    pause(0.1)
    hold off
end