Но на таких картинках для планет Солнечной системы Вы заметите глазом отклонение от окружности разве что для орбиты Меркурия (и то, пожалуй, с трудом - одна ось будет больше другой на пару миллиметров).
Так и получается. Нарисовал эллипс для Земли, Меркурия и Марса. У меня получилось так :
Все выглядят как окружности. Но на различных картинках в интернете. Рисуют так:
Наверное, я где-то ошибся. Для рисования я использовал программу Scilab. Рисовал так :
(Оффтоп)
Код:
function b=ab(a1,e)
b=sqrt((a1^2)-(e^2)*(a1^2))
endfunction
function b=elipseMy(a1,b1,t)
x1=a1*cos(t)
y1=b1*sin(t)
b=list(x1,y1)
endfunction
function km=ae(val)
km=val*149599300
endfunction
pt=[0:0.001:2*%pi]
////////////////
//солнце только окружность
SunR = (1.392*10^9)/(2*1000)
SunA=SunR
SunB=SunR
SunO=elipseMy(SunA,SunB,t)
//////////////////////
//юпитер
//большая полуось
JupA=ae(5.20288700)
//эксцентриситетом.
JupE=0.048775
//малая полуось
JupB=ab(JupA,JupE)
printf("Jupiter a=%f, b=%f, a/b=%f \n",JupA,JupB,JupA/JupB)
JupO=elipseMy(JupA,JupB,t)
/////////////////////////
//земля
//большая полуось
EaA=ae(1.00000261)
//эксцентриситетом.
EaE=0.01671123
//малая полуось
EaB=ab(EaA,EaE)
EaO=elipseMy(EaA,EaB,t)
printf("Eath a=%f, b=%f, a/b=%f \n",EaA,EaB,EaA/EaB)
/////////////////////////
//Меркурий
//большая полуось
MerA=ae(0.38709927)
//эксцентриситетом.
MerE= 0.20563593
//малая полуось
MerB=ab(MerA,MerE)
MerO=elipseMy(MerA,MerB,t)
printf("Mercury a=%f, b=%f, A/B=%f \n",MerA,MerB,MerA/MerB)
/////////////////////////
//Марс
//большая полуось
MarA=ae(1.52371034)
//эксцентриситетом.
MarE= 0.09339410
//малая полуось
MarB=ab(MarA,MarE)
MarO=elipseMy(MarA,MarB,t)
printf("Mars a=%f, b=%f, A/B=%f \n",MarA,MarB,MarA/MarB)
square(-1,-1,1,1)
a=gca();
a.data_bounds=[-10*10^8,-10*10^8; 10*10^8,10*10^8]
a.box="on"
a.grid=[2,2]
a.auto_margins="on"
plot(SunO(1),SunO(2),"y->")
poly1= a.children.children(1)
poly1.thickness = 5000
plot(JupO(1),JupO(2),"g")
plot(EaO(1),EaO(2),"b")
plot(MerO(1),MerO(2),"m")
plot(MarO(1),MarO(2),"r")
Вывод:
Код:
Jupiter a=778348253.179100, b=777421856.385990, a/b=1.001192
Eath a=149599690.454173, b=149578800.001233, a/b=1.000140
Mercury a=57909779.822511, b=56672164.610569, A/B=1.021838
Mars a=227946000.266762, b=226949698.260825, A/B=1.004390
Данные брал отсюда:
http://www.webcitation.org/6BOqmsTQGИз таблицы 8.10.2 .
Что я не так сделал?