동력학-공기저항시포물선운동궤적

소개글

동력학-공기저항시포물선운동궤적에 대한 보고서 자료입니다.

목차

1. 계산식의 유도
2. 프로그래밍

본문내용

axis()
변수 범위 지정
hold on
전의 그래프와 함께 고정시켜 표시
legend()
범례의 정의
end
프로그래밍 종료
이제 이들을 조합해보면
c=[0 0.002]; m=100; g=9.81; dt=0.1;
c, m, g, dt 값을 정의
x=zeros(1000);y=zeros(1000);vx=zeros(1000);vy=zeros(1000);
x, y, vx, vy를 0이 1000개인 벡터로 만든다.
x(1)=0; y(1)=0; vx(1)=100; vy(1)=100; k=1;
x, y, vx, vy, k의 초기 값 지정
A=char('B-', 'G-');
그래프 색 지정( 파랑, 초록 순)
for i=1:2
c가 1부터 2까지 반복실행
for k=1:2000
k가 1부터 2000까지 반복실행
vx(k+1)=vx(k)-(c(i)/m*sqrt(vx(k)^2+vy(k)^2)*vx(k))*dt;
유한차분법에 의한 vx 정의
vy(k+1)=vy(k)-(g+c(i)/m*sqrt(vx(k)^2+vy(k)^2)*vy(k))*dt;
유한차분법에 의한 vy 정의
x(k+1)=x(k)+vx(k)*dt;
유한차분법에 의한 x 정의
y(k+1)=y(k)+vy(k)*dt;
유한차분법에 의한 y 정의
end
k 반복문 종료
if i==1
조건 : c가 첫 번째 값일 때
plot(x,y, A(i,:))
x, y 그래프를 A정의대로 그린다.
title('운동하는 물체의 궤적')
제목 : 운동하는 물체의 궤적
ylabel('y(t)')
y축 라벨 'y(t)'
xlabel('x(t)')
x축 라벨 'x(t)'
axis([0 2400 0 700])
범위지정(x:0~2400, y:0~700)
hold on
그래프를 연이어 표시한다.
else
조건이 아닐 때 : 두 번 째 c값
plot(x,y, A(i,:))
x, y 그래프를 A 정의대로 그림
end
조건문종료
legend('c=0','c=0.002')
범례 정의(c=0, 0.002일 때)
end
c 반복문 종료
<매트랩으로 프로그래밍 한 모습>
④ 결과를 도출한다.
→ 위 과정대로 프로그래밍 했을 때 도출되는 그래프의 형상은 다음과 같다.
c=0의 그래프는 공기저항을 무시했을 때, c=0.002는 공기저항 존재시의 궤적의 모습이다.
공기 저항이 존재할 때 덜 멀리 날아갔음을 알 수 있다.
그리고 구간 사이의 시간 차이인 를 더욱 작게 할수록 그래프의 형태는 더 정확해질 것이다.

키워드

매트랩,   동력학,   포물선운동,   공기저항,   항력계수