[공학]Matlab을 이용하여 화성과 지구의 궤도 그리기

소개글

6가지 전통적인 요소인 장반경, eccentricity,perigee 통과시간, ascending node의 적경(right ascension), 궤도면 경사각, 그리고 perigee의 위치를 통해 행성들의 궤도를 그릴 수 있는,
matlab을 이용하여 알고리즘 짜기

목차

없음

본문내용

그 동안 강의를 통해 배운 내용으로 화성과 지구의 날에 따른 거리를 구하는 프로그램을 짜보았다. 거리를 구하려면 먼저 태양을 초점으로 하는 화성과 지구 각각의 날짜에 따른 위치를 알아야한다. 이는 화성과 지구의 궤도를 구해서 날짜에 따른 를 구한다는 말이 된다.
궤도를 구하기 위해서 전통적인 6가지 요소를 사용하게 되는데, 그 첫 번째 요소로 장반경 와 eccentricity , 그리고 perigee 통과시간 가 있다. 이 세 요소가 주어지면 궤도면에서의 위성의 위치를 알 수 있다. 두 번째 요소 그룹은 궤도면을 표현하는데 사용하는데, 이는 ascending node의 적경(right ascension) , 궤도면 경사각 , 그리고 perigee의 위치를 나타내는 이다. 이에 대한 자세한 내용은 모 기관의 홈페이지에서 찾았는데 그 값은 다음 표에 나타내었다.

clear

%%Initial values
T=1:input(`T =`);
t_e=100.46457166*pi/180; t_m=-4.55343205*pi/180 ; %Mean longitude
a_e=1.496*10^8; a_m=2.2792e+008; %Semi-major axis
e_eccent=0.0167175; m_eccent=0.0933865; %% Orbit eccentricity
mu=1.32715*10^11; %%Gravitational parameter

W_e=102.93768193*pi/180; W_m=-23.94362959*pi/180;
%Longitude of the ascending node
w_e=102.93768193*pi/180; w_m=-73.5031682*pi/180;

o_e=0; o_m=49.55953891*pi/180; %Argument of perihelion
i_e=0; i_m=1.84969142*pi/180; %Orbit inclination

t=T*3600*24;

%Compute the mean anomaly
M_e=t_e-W_e+sqrt(mu/a_e^3)*t;
M_m=t_m-W_m+sqrt(mu/a_m^3)*t;

참고 자료

없음