Plan Omega to Alpha - POA

각도에 따른 포물선 운동 (공기 마찰력) 공기의 마찰력을 속도의 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

감소 운동 (시간에 따라서 마찰력에 의해 속도가 감소하는 운동) 감소하는 속도는 현재의 속도의 10 %의 비율로 감소한다. hold off x=0; dxdt=input('초기 속도를 입력하시오'); dt=0.1; t=10; for k=0:1:t/dt;     time=dt*k;     dxdt=dxdt-0.1*dxdt*dt;     x=x+dxdt*dt;     plot(x,0,'-ro');     pause(0.1)     title([num2str(time),'초'])     hold on    end; hold of f