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