Plan Omega to Alpha - POA

포물선 운동 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

포물선운동 (중력 가속도 포함) hold off; g=9.81; v=input('초기 속도를 입력하시오.'); theta=input('초기 각도를 입력하시오.'); dt=input('시간 간격을 입력하시오'); x=0; y=0; X(1)=0; Y(1)=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); while y>=0        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);         X(k+1)=x;     Y(k+1)=y;          end; z=k+1; q=t2/dt     for p=1:1:z;    x1=X(p);    y1=Y(p);    plot(0.0)    hold on    axis([0,X(k+1),0,Y(round(q+10))]);  ...

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

포물선운동 계산 (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; for k=0:1:t1/dt;        plot(x,y,'-ro')     time= dt*k;         dxdt=dxdt;     dydt=dydt-g*dt;         x=x+(dxdt*dt);     y=y+(dydt*dt);    hold on       title([num2str(time),'초'])    pause(0.1) end