Реферат: Тригонометричні ефемериди планет Сонячної системи
while EG>=360 do
begin
EG:=EG-360
end;
ER1:=(Eao*(1-sqr(Eeo)))/(1+Eeo*cos((EG-Eap)*PI/180)); //формула 2
ER:=ER1;
Until DR<=ED; //обчислення до моменту спостереження
end;
end.
unit Mars;
interface
Uses KonstPLN,Eagth,Math;
Const AD5=32852; //номер дня 9 січня 1990р. з 1 березня 1900р.
AG5=235.0; AR5=1.538; //початкове положення Марса 9 січня 1990р.
Procedure Mars1(var AD,AG,AR,x,AR0,ASD,AH,AR1,AGD,RM,EG:real);
Procedure Mars2(var AAD1,ARS,AZE,AZP,ABD,AV,AYD,AYDS,AAD0,AAD,AA,AAM,AAS,AID1,AID2,AID,AB,
ABM,AF, ER,AG,EG,AR,AAH,ABG:real;var R0:integer);
Procedure Mars3(var AT1,ATK,ATKL,ATKLM,AP,AT2,AT3,ATG,ATS,ATSL,ATSLM,ATZ,ATZL,ATZLM,AZ,
AA,EG,k,LB,AYD,AV,FI,AB,ATKLH,ATSLH,ATZLH:real;var n:integer);
implementation
Procedure Mars1;
begin
AD:=AD5; AG:=AG5; AR:=AR5; x:=0;
AR0:=Avo*86400/AO*N0; //формула 4 шлях проходження за N4 дні
ASD:=(PI*sqr(Aao)*sqrt(1-sqr(Aeo))/Atd)*N0;//формула 3 площа еліпса орбіти за N4 дні
Repeat
AD:=AD+N0; //2 закон Кеплера
AH:=2*ASD/AR; //радіус-вектор через N4 дні
AR1:=sqrt(sqr(AH)+sqr(AR-sqrt(abs(sqr(AR0)-sqr(AH)))));//формула 5
AGD:=arctan((1/AR1*AH)/sqrt(1-sqr(1/AR1*AH)))*180/PI;//формула 6
AG:=AG+AGD; //кут проходу за N4 дні
While AG>=360 do
begin
AG:=AG-360
end;
AR1:=(Aao*(1-sqr(Aeo)))/(1+Aeo*cos((AG-Aap)*PI/180));//формула 2
AR:=AR1;
x:=x+1 //кількість циклів обчислення
Until RM<=AD; //обчислення до моменту спостереження
end;//Mars1
Procedure Mars2;
begin
AAD1:=arctan((ER*sin((AG-EG)*PI/180))/(AR-ER*cos((AG-EG)*PI/180)))*180/PI;//зміщення прямого
ARS:=sin((AG-Avy)*PI/180)*sin(Aie*PI/180); //піднесення, формула 7
AZE:=AR*ARS;
AZP:=sqrt(sqr(ER*sin((EG-AG)*PI/180))+sqr(AR*sqrt(1-sqr(ARS))-ER*cos((EG-AG)*PI/180)));
ABD:=arctan(AZE/AZP)*180/PI; //зміщення схилення, формула 8
AV:=sqrt(sqr(AZE)+sqr(AZP));//лінійна відстань між Землею і Марсом, формула 9
AYD:=2*arctan(Are/(AV*AO))*180/PI;//кутовий діаметр Марса, формула 10
AYDS:=AYD*3600;
AAD0:=arctan(sqrt(1-sqr((sqr(AV)+sqr(AR)-sqr(ER))/(2*AV*AR)))/((sqr(AV)+sqr(AR)-sqr(ER))/
(2*AV*AR)))*180/PI;
AAD:=sqrt(abs(sqr(AAD0)-sqr(ABD)));
If AAD1<0 then
begin
AAD:=-AAD;
end;
AA:=(AG+AAD)/15; //пряме піднесення Марса
If AA>=24 then
begin
AA:=AA-24;
end;
If AA<0 then
begin
AA:=AA+24;
end;
AAH:=int(AA);
AAM:=frac(AA)*60;
AAS:=frac(AAM)*60;
AID1:=sin((AA*15)*PI/180)*(Eri-0.00013011111*(R0-1900));//схилення точки екліптики
AID2:=arctan((sin((180-AA*15)*PI/180)*sin((Eri-0.00013011111*(R0-1900))*PI/180))/
(sqrt(1-sqr(sin((180-AA*15)*PI/180)*sin((Eri-0.00013011111*(R0-1900))*PI/180)))))*
180/PI; //те саме другим методом
AID:=(AID1+AID2)/2; //їх середнє значення
AB:=AID+ABD; //Схилення Марса
ABG:=int(AB);
ABM:=frac(AB)*60;
AF:=(1/2)+(cos(AAD0*PI/180)/2); //фаза диска Марса, формула 11
end;//Mars2
Procedure Mars3;
begin
AT1:=AA-(EG/15+k);
If AT1<0 then
begin
AT1:=AT1+24;
end;
ATK:=AT1/1.002738; //час кульмінації Марса, формула 14
ATKL:=ATK-LB/15+n; //місцевий час кульмінації
if ATKL<0 then
begin
ATKL:=ATKL+24;
end;
ATKLH:=int(ATKL); ATKLM:=frac(ATKL)*60;
AP:=90+0.59+(AYD/2)-(arctan(Ere/(AV*AO))*180/PI);
AT2:=(cos(AP*PI/180)-sin(FI*PI/180)*sin(AB*PI/180))/(cos(FI*PI/180)*cos(AB*PI/180));
AT3:=arctan(sqrt(1-sqr(AT2))/AT2)*180/PI;//годинний кут сходу і заходу Марса
if AT3<=0 then //формула 12
begin
AT3:=AT3+180;
end;
ATG:=AT3/15;
ATS:=ATK-ATG; //час сходу Марса
if ATS<0 then
begin
ATS:=ATS+24;
end;
ATSL:=ATS-LB/15+n; //місцевий час сходу
if ATSL<0 then
begin
ATSL:=ATSL+24;
end;
ATSLH:=int(ATSL); ATSLM:=frac(ATSL)*60;
ATZ:=ATK+ATG; //час заходу Марса
if ATZ>=24 then
begin
ATZ:=ATZ-24;
end;
ATZL:=ATZ-LB/15+n; //місцевий час заходу
if ATZL<0 then
begin
ATZL:=ATZL+24;
end;
ATZLH:=int(ATZL); ATZLM:=frac(ATZL)*60;
AZ:=arctan((cos(AB*PI/180)*sin(AT3*PI/180)/sin(AP*PI/180))/sqrt(1-sqr((cos(AB*PI/180)*
sin(AT3*PI/180))/sin(AP*PI/180))))*180/PI; //азимут Марса при сході і заході
if AB>0 then //формула 13
begin
AZ:=180-AZ;
end;
end;//Mars3
end.
unit Common;
interface
var
ZPL,ZD,ZM,ZR,ZPL1,ZD1,ZM1,ZR1 : integer;
ZAA,ZAB,ZAA1,ZAB1 :real;
implementation
end.
unit Unit3;
interface
uses
Windows, Messages, SysUtils, Variants, Classes, Graphics, Controls, Forms,
Dialogs, ExtCtrls, StdCtrls, Buttons, Common,Math,Unit4;
type
TForm3 = class(TForm)
Image1: TImage;
SpeedButton1: TSpeedButton;
SpeedButton2: TSpeedButton;
Button1: TButton;
Label1: TLabel;
-----------------------
Label31: TLabel;
StaticText1: TStaticText;
procedure FormCreate(Sender: TObject);
procedure Button1Click(Sender: TObject);
procedure SpeedButton1Click(Sender: TObject);
procedure SpeedButton2Click(Sender: TObject);
private
{ Private declarations }
public
{ Public declarations }
end;
Mas = array[1..311] of real;
var
Form3: TForm3;
ZA,ZB,ZP : Mas;
ZAG,ZBG : integer;
ZAM,ZBM,x,y,x3,y3 : real;
x4,y4,x5,y5,x0,y0,x1,y1,x2,y2 :variant;
implementation
{$R *.dfm}
procedure TForm3.FormCreate(Sender: TObject);
var i:integer;
begin
SpeedButton2.visible:=false;
Label1.Visible:=false;
Label2.Visible:=false;
Label3.Visible:=false;
Label4.Visible:=false;
Label5.Visible:=false;
Label6.Visible:=false;
Страницы: 1, 2, 3, 4, 5, 6, 7, 8
