Тригонометричні ефемериди планет Сонячної системи

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;

Image1.Canvas.Brush.Color:=clBlack;

Image1.Canvas.Brush.Style:=bsSolid;