Астрономический анализ хронологии

ПРИЛОЖЕНИЕ 5. Компьютерная программа геометрического метода датировки звездных конфигураций по собственным движениям с учетом систематической ошибки каталога

!!!!!!!!!!!! Написать описание программы, ее процедур, входных данных и выдаваемого результата.

Компьютерная программа геометрического метода датировки звездных конфигураций по собственным движениям с учетом систематической ошибки каталога =====================================

program perebor; {программа написана на языке Pascal под Delphi4.0}

uses Math;

const

nstar1 = 300; {ограничение числа звезд в конфигурации}

pi = 3.1415926536; {значение константы }

deltaGM = 5; {размах перебора gamma вокруг _stat при поиске

оптимального поворота (в минутах)}

deltaBM = 30; {размах перебора beta вокруг нуля при поиске

оптимального поворота (в минутах)}

gstepM = 1.0; {шаг поиска оптимальной точки по gamma (в минутах)}

bstepM = 1.0; {шаг поиска оптимальной точки по beta (в минутах)}

eps = 30; {окрестность захвата для подсчета хорошо

приблизившихся по широте звезд (в минутах)}

d8 = 900000; {наибольшее допустимое расстояние от звезды до

ближайшей из 8-ми именных }

type

cr1=record

nb : integer;

a,d,va,vd,l,b,cb,sb,Mbs5,Malm : real;

obozn : string;

end;

var

co : array[1..nstar1] of cr1;

ah,am,asec,dg,dm,ds,va,vd,lg,lm,bg,bm,e,ce,se,

lx,clx,slx,bx,cbx,sbx,ly,cly,sly,by,cby,sby,

e1,se1,ce1,ft,ps,mg,maxb1,maxb2,angle,cangle,sangle,

x,y,gr,deltl,ymin,ymax,gstep,bstep,cgstep,

sgstep,cbstep,sbstep,bmax,gamma0,beta0,dl0,dist0,

cminmax,cc,fmax,fminmax,fx,y1,dist1,dBm,dBmm,

deltaG,cdeltaG,sdeltaG,deltaB,cdeltaB,sdeltaB,

cGstat,sGstat,xd1,xd2,d8rad,epsrad : real;

stt,stm,stf : array [1..nstar1,1..6] of real;

Gstat : array [1..30] of real; {значения gamma_stat

полученные из статистической процедуры оценивания}

zv,zvv : array [1..nstar1] of integer;

id : array [1..nstar1] of integer; {признак удержания

звезды за счет близости к 8-ми зведному ядру:

0 - отбрасывается, 1 - удерживается}

agamt,cgamt,sgamt,abett,cbett,sbett : real;

nb,i,j,t,t1,t2,nstar,Ngamma,Nbeta,Ng0,Nb0,ig,ib,Nstep,

Iok,Itek,NBmm,NBm,jj,jj1,i8 : integer;

f,f1,f2 : text;

konec : char;

{************************************************************}

{* vvod *}

{************************************************************}

procedure vvod;

var i : integer;

Mbs5,Malm : real;

ob : string;

begin

assign(f1,’result.txt’);

rewrite(f1);

assign(f2,’sig-max.txt’);

rewrite(f2);

writeln(f1,’ *** Program perebor.pas ***’);

writeln(’ *** Program perebor.pas ***’);

assign(f,’fast.txt’); {fast.txt - ВХОДНОЙ ФАЙЛ С ДАННЫМИ О ЗВЕЗДАХ}

reset(f);

{*********** ЧТЕНИЕ ДАННЫХ **************}

nstar:=0;

while not eof(f) do

begin {while}

nstar:= nstar+1;

i:=nstar;

readln(f,nb,ah,am,asec,dg,dm,ds,Mbs5,va,vd,lg,lm,bg,bm,Malm,ob);

{++++++ строение строки данных в файле fast.txt ++++++++}

{ nb - номер звезды по BS5, }

{ ah - прямое восхождение (часы), }

{ am - прямое восхождение (часовые минуты) ЗНАК ОПУЩЕН, }

{ asec - прямое восхождение (часовые секунды) ЗНАК ОПУЩЕН, }

{ dg - склонение (градусы), }

{ dm - склонение (дуговые минуты), ЗНАК ОПУЩЕН }

{ ds - склонение (дуговые секунды), ЗНАК ОПУЩЕН }

{ va - скорость собственного движения в прямом восхождении,}

{ приведенная к экватору ("/год), }

{ vd - скорость собственного движения в склонении }

{ ("/год), }

{ lg - долгота по Альмагесту (градусы), }

{ lm - долгота по Альмагесту (минуты), НЕОТРИЦАТЕЛЬНА }

{ bg - широта по Альмагесту (градусы), }

{ bm - широта по Альмагесту (минуты) ЗНАК ОПУЩЕН }

{ Mbs5 - величина (яркость) по BS5 }

{ Malm - величина (яркость)по Альмагесту }

{ ob - современное обозначение звезды }

if (ah<0) then

begin

am:= -am;

asec:=-asec;

end;

if (dg<0) then

begin

dm:= -dm;

ds:=- ds;

end;

if (bg<0) then bm:= -bm;

co[i].nb:=nb;

co[i].a:=pi*(ah+am/60+asec/3600)/12;

co[i].d:=pi*(dg+dm/60+ds/3600)/180;

co[i].va:=va*pi/6480.0; {перевод скоростей собств. движения: }

co[i].vd:=vd*pi/6480.0; {секунды/год->радианы/100лет }

co[i].l:=pi*(lg+lm/60)/180;

co[i].b:=pi*(bg+bm/60)/180;

co[i].Malm:=Malm;

co[i].Mbs5:=Mbs5;

co[i].obozn:=ob;

co[i].cb:=cos(co[i].b);

co[i].sb:=sin(co[i].b);

if co[i].cb <> 0 then

co[i].va:=co[i].va/co[i].cb;{теперь скорость НЕ приведена к экватору}

writeln(f1,nb:4,’ ’,ah:4:0,’ ’,am:6:2,’ ’,

dg:4:0,’ ’,dm:6:2,’ ’,

lg:4:0,’ ’,lm:4:0,’ ’,bg:4:0,’ ’,bm:4:0,’ ’,

Malm:3:1,’ ’,Mbs5:3:1,’ ’,ob);

writeln(nb:4,’ ’,ah:4:0,’ ’,am:6:2,’ ’,

dg:4:0,’ ’,dm:6:2,’ ’,

lg:4:0,’ ’,lm:4:0,’ ’,bg:4:0,’ ’,bm:4:0,’ ’,

Malm:3:1,’ ’,Mbs5:3:1,’ ’,ob);

end; {while}

writeln(’nstar= ’,nstar);

writeln(f1,’FAST.TXT: nstar= ’,nstar);

writeln(f1);

{for i:=1 to nstar do

writeln(f1,co[i].nb:4:0,’ ’,co[i].a:7:5,’ ’,co[i].d:7:5,

’ ’,co[i].l:7:5,’ ’,co[i].b:7:5); }

writeln(’VVOD’ );

end; {vvod}

{************************************************************}

{* TURN *}

{************************************************************}

procedure turn;

{lx (clx, slx) - долгота (cos, sin) до поворота,

bx (cbx, sbx) - широта (cos, sin) до поворота,

ly (cly, sly) - долгота (cos, sin) после поворота,

by (cby, sby) - широта (cos, sin) после поворота,

angle (cangle,sangle) - угол (cos,sin) поворота}

var

c,x,y : real;

begin {turn}

sby:= -slx*cbx*sangle + sbx*cangle;

cby:= sqrt(1 - sqr(sby));

if sby=1 then by:= pi/2

else by:= arctan(sby/cby);

c:= cbx*clx;

if c = 0 then

begin

if cbx*cangle+slx*sbx*sangle > 0 then ly := lx

else ly:=lx-pi;

if cbx = 0 then ly:= pi/2;

end

else {если c не равно нулю}

begin

ly:= (slx*cbx*cangle + sbx*sangle)/c;

ly:= arctan(ly);

if ly < 0 then ly:= ly + pi;

{if ly > pi then writeln(’!!!!!!!!!!’); }

{–––––––––––––––––––––––––}

{Если звезда находится в сферическом круге, построенном как на

диаметре на дуге длины angle, соединяющем старый и новый полюс,

то модуль разности ее старой и новой долгот ближе к pi, чем к

нулю. Если же она находится вне этого круга, то модуль разности

долгот ближе к нулю, чем к pi}

y:=pi/2 - bx;

x:=angle*cos(lx+pi/2); {Для скорости счета взята оценка.

На самом деле:

angle*cos(lx+pi/2) <= x <= angle }

if y>x then

begin

if abs(abs(lx-ly)-pi)<pi/2 then ly:=ly+pi;

end

else

begin

if abs(lx-ly)<pi/2 then ly:=ly+pi;

end;

{–––––––––––––––––––––––––}

end; {если c не равно нулю}

cly:= cos(ly);

sly:= sin(ly);

if ly > 2*pi then ly:=ly-2*pi;

if ly < 0 then ly:=ly+2*pi;

end; {turn}

{****************************************************************}

{* PERESCHET NA VREMYA Т *}

{****************************************************************}

procedure pereschet;

var i: integer;

z,zz: real;

{выдает: stt[i,1] = l

stt[i,2] = cos(l)

stt[i,3] = sin(l)

stt[i,4] = b

stt[i,5] = cos(b)

stt[i,6] = sin(b)

где l,b - эклиптикальные координаты звезды

в эпоху t (с учетом собственного движения)}

begin {pereschet}

for i:= 1 to nstar do

begin {for i}

lx := co[i].a + t1*co[i].va;

clx:= cos(lx);

slx:= sin(lx);

bx := co[i].d + t1*co[i].vd;

sbx:= sin(bx);

cbx:= sqrt(1 - sqr(sbx));

cangle:= ce;

sangle:= se;

angle:=e;

turn;

bx := by;

cbx:= cby;

sbx:= sby;

lx:= ly - ft;

if lx < 0 then lx:= lx + 2*pi;

clx:= cos(lx);

slx:= sin(lx);

cangle:= ce1;

sangle:= se1;

angle:=e1;

turn;

stt[i,4]:= by;

stt[i,5]:= cby;

stt[i,6]:= sby;

lx:= ly + ft + ps;

if lx > 2*pi then lx:= lx - 2*pi;

if lx <= -2*pi then lx:= lx + 2*pi;

if lx > 2*pi then lx:= lx - 2*pi;

if lx <= -2*pi then lx:= lx + 2*pi;

stt[i,1]:= lx;

stt[i,2]:= cos(lx);

stt[i,3]:= sin(lx);

{–––––––––––––––––

zz:=mg/60;

z:= (stt[i,1]-co[i].l)*zz;

writeln(f1,co[i].nb:4,’ ’,’L= ’,lx*zz:5:3,’; B= ’,by*zz:5:3);

writeln(co[i].nb:4,’ ’,’L= ’,lx*zz:5:3,’; B= ’,by*zz:5:3);

if abs(z)> 20 then

begin

writeln(f1,’dL=’,z:10:1,’(gr); i= ’,co[i].nb,’ L-alm=’,co[i].l*zz:6:2,

’ B-alm=’,co[i].b*zz:6:2);

writeln(”,’dL=’,z:10:1,’(gr); i= ’,co[i].nb,’ L-alm=’,co[i].l*zz:6:2,

’ B-alm=’,co[i].b*zz:6:2);

end;

z:= (stt[i,4]-co[i].b)*mg;

if abs(z)> 300 then

begin

writeln(f1,’ ’,’dB= ’,z:10:1,’(min); i= ’,i);

writeln(’ ’,’dB= ’,z:10:1,’(min); i= ’,i);

end;

–––––––––––––––––– }

end; {for i}

end; {pereschet}

{************************************************************}

{* DIST (расстояние между точками сферы в радианах) *}

{************************************************************}

function dist(L1:real;B1:real;L2:real;B2:real) : real;

{L1,B1 - долгота и широта первой точки,

L2,B2 - долгота и широта второй точки}

var

X1,X2,Y1,Y2,Z1,Z2,DE,DSIN,DTAN : real;

begin {dist}

X1 := COS(B1)*COS(L1);

Y1 := COS(B1)*SIN(L1);

Z1 := SIN(B1);

X2 := COS(B2)*COS(L2);

Y2 := COS(B2)*SIN(L2);

Z2 := SIN(B2);

DE:=SQRT(SQR(X1-X2)+SQR(Y1-Y2)+SQR(Z1-Z2));

DSIN:= DE/2;

DTAN:=DSIN/SQRT(1.0-SQR(DSIN));

Result:= 2.0*ARCTAN(DTAN);

end;{dist}

{**************************************************************}

{ MAIN PROGRAM }

{**************************************************************}

begin {program}

{***********************************}

vvod; {ввод данных для звезд из файла fast.txt}

{***********************************}

mg:= 180.0*60.0/pi; {к-т пересчета минут в радианы и наоборот}

e:=pi*(23+27/60+8.26/3600)/180; {угол накл. экл. к экватору для t=0}

se:=sin(e);

ce:=cos(e);

d8rad:=d8/mg;

epsrad:=eps/mg;

{––––––––––––––––––––––}

for i:=1 to nstar do

begin

xd1:=10;

for i8:=1 to 8 do {8 звезд информативного

ядра должны стоять вначале!}

begin

xd2:=dist(co[i8].a,co[i8].d,co[i].a,co[i].d);

{ writeln(f1,co[i].nb,’ dist (min) = ’,xd2*mg:4:1); }

if xd2 < xd1 then xd1:=xd2;

end;

xd2:=xd1*mg/60;

{ writeln(f1,co[i].nb,’ dist (grad) = ’,xd2:4:1); }

if xd1 < d8rad then id[i]:=1 else id[i]:=0;

end;

{––––––––––––––––––––––}

gstep:=gstepM/mg;

bstep:=bstepM/mg;

cgstep:= cos(gstep);

sgstep:= sin(gstep);

cbstep:= cos(bstep);

sbstep:= sin(bstep);

deltaG:=deltaGM/mg;

cdeltaG:= cos(deltaG);

sdeltaG:= sin(deltaG);

deltaB:=deltaBM/mg;

cdeltaB:= cos(deltaB);

sdeltaB:= sin(deltaB);

Ngamma:=Trunc(deltaG/gstep); {число шагов по gamma в одну сторону}

Nbeta:= Trunc(deltaB/bstep); {число шагов по beta в одну сторону}

Gstat[1]:= 30.5/mg;

Gstat[2]:= 29.5/mg;

Gstat[3]:= 28.5/mg;

Gstat[4]:= 27.5/mg;

Gstat[5]:= 27.0/mg;

Gstat[6]:= 26.0/mg;

Gstat[7]:= 25.2/mg;

Gstat[8]:= 24.4/mg;

Gstat[9]:= 23.5/mg;

Gstat[10]:= 22.6/mg;

Gstat[11]:= 21.8/mg;

Gstat[12]:= 21.0/mg;

Gstat[13]:= 20.4/mg;

Gstat[14]:= 19.5/mg;

Gstat[15]:= 18.8/mg;

Gstat[16]:= 18.0/mg;

Gstat[17]:= 17.2/mg;

Gstat[18]:= 16.4/mg;

Gstat[19]:= 15.8/mg;

Gstat[20]:= 15.0/mg;

Gstat[21]:= 14.4/mg;

Gstat[22]:= 13.8/mg;

Gstat[23]:= 13.1/mg;

Gstat[24]:= 12.5/mg;

Gstat[25]:= 12.0/mg;

Gstat[26]:= 11.5/mg;

Gstat[27]:= 11.1/mg;

Gstat[28]:= 10.8/mg;

Gstat[29]:= 10.5/mg;

Gstat[30]:= 10.2/mg;

writeln(f2,’ t ’,’sigma ’,’maxB’,’ N-in-eps’);

for t:=1 to 30 do {цикл по времени в прошлое с шагом 1=столетие}

begin {for t}

{ writeln(f1,’T = ’,t:2);

writeln(f1); }

writeln(’T = ’,t:2);

writeln;

t1:=-t;

e1:=(pi/648000.0)*(47.070559+(-0.033769+0.00005*t1)*t1)*t1;

se1:=sin(e1);

ce1:=cos(e1);

ft:=(pi/180.0)*(174+52/60.0 -t1*870.0798/3600.0+t1*t1*0.024578/3600.0);

ps:=(pi/648000.0)*(5026.872+(1.131358+0.000102*t1)*t1)*t1;

{******************************************************}

pereschet; {пересчет координат звезд на эпоху t }

{******************************************************}

cGstat:=cos(Gstat[t]);

sGstat:=sin(Gstat[t]);

angle:= Gstat[t]-deltaG;

cangle:= cdeltaG*cGstat+sdeltaG*sGstat;

sangle:= sGstat*cdeltaG -sdeltaG*cGstat;

{текущий угол поворота по gamma - косинус и синус

устанавливается в начале равным Gstat[t]-deltaG}

{cgamt,sgamt - косинус и синус накопленного угла поворота по gamma}

{cbett,sbett - косинус и синус накопленного угла поворота по beta}

bmax:=1; {заготовка для минимума по поворотам максимальной

широтной невязки по звездам}

dBmm:=1; {заготовка для минимума по поворотам средней

широтной невязки по звездам}

Nbmm:=0; {заготовка для максимума по подкруткам числа

звезд, попавших в eps’ - окрестность альмагестовской}

gamma0:=0; {заготовка для оптимального поворота по gamma}

beta0:=0; {заготовка для оптимального поворота по beta}

dl0:=0; {заготовка для размаха по долготе при минимаксе по широте}

dist0:=0; {заготовка для невязки по дуге при минимаксе по широте}

for ig:=-Ngamma to Ngamma do

begin {for ig - поворот вдоль}

{ writeln(’ig = ’,ig); }

i:=1;

while (i <= nstar) do

begin {while i<=nstar}

lx := stt[i,1];

clx:= stt[i,2];

slx:= stt[i,3];

bx := stt[i,4];

cbx:= stt[i,5];

sbx:= stt[i,6];

turn;

if ly > 3.0*pi/2.0 then x:= ly-2.0*pi else x:=ly;

stm[i,1]:= x+pi/2;

stm[i,2]:= -sly;

stm[i,3]:= cly;

stm[i,4]:= by;

stm[i,5]:= cby;

stm[i,6]:= sby;

i:=i+1;

end; {while i<=nstar}

agamt:=angle;

cgamt:=cangle;

sgamt:=sangle; {запоминаем накопленный угол поворота по gamma,

чтобы вернуться к нему после цикла поворотов

поперек}

angle:= -deltaB;

cangle:= cdeltaB;

sangle:= -sdeltaB; {в начале цикла поворотов по beta устанавли-

ваем угол поворота равным -deltaB}

for ib:= -Nbeta to Nbeta do

begin {for ib - поворот поперек}

i:=1;

maxb1:=0.0;

ymin:=7.0;

ymax:=-7.0;

dBm:=0;

Nbm :=0;

while (i <= nstar) do

begin {while i<=nstar}

lx := stm[i,1];

clx:= stm[i,2];

slx:= stm[i,3];

bx := stm[i,4];

cbx:= stm[i,5];

sbx:= stm[i,6];

turn;

stf[i,2]:=by;

stf[i,3]:=cby;

if ly < pi/2 then y:=ly + 2*pi else y:= ly;

stf[i,1]:=ly - pi/2;

y:= y - pi/2 - co[i].l;

if y < -pi then y:=y+2*pi

else if y> pi then y:= y-2*pi;

if y < -pi then y:= y+2*pi

else if y>pi then y:=y-2*pi;

y1:=y*cby;

if abs(y1)>0.5 then

begin

writeln(f1,’dL*cosB=’,y1:10:5,’(rad); N(BS5)=’,co[i].nb:4);

writeln(’dL*cosB=’,y1:10:5,’(rad); N(BS5)=’,co[i].nb:4);

writeln(f1,’cosB=’,cby:10:5);

writeln(’cosB=’,cby:10:5);

x:=mg/60;

writeln(f1,’by=’,by*x:9:2,’ ly=’,ly*x:9:2);

writeln(’by=’,by*x:9:2,’ ly=’,ly*x:9:2);

writeln(f1,’L-alm=’,co[i].l*x:9:2,’ B-alm=’,co[i].b*x:9:2);

writeln(’L-alm=’,co[i].l*x:9:2,’ B-alm=’,co[i].b*x:9:2);

readln(konec);

end;

stf[i,4]:=y;

if y < ymin then ymin:= y;

if y > ymax then ymax:= y;

{–––-1-й вариант: ядро из 8 звезд всегда удерживается –––––-}

maxb2:= abs(by - co[i].b);

if (id[i]=1) and (maxb2 < epsrad) then

begin

dBm:=dBm+sqr(maxb2);

NBm:=NBm+1;

zv[NBm]:=i;

end;

if maxb2 > maxb1 then

begin

maxb1:= maxb2;

Itek:=i

end;

i:= i+1;

end; {while i<=nstar}

dBm:=sqrt(dBm/NBm);

{–––––2-й вариант: ядро не выделено при удержании –––––}

{ maxb2:= abs(by - co[i].b);

dBm:=dBm+sqr(maxb2);

if maxb2*mg<eps then NBm:=NBm+1;

if maxb2 > maxb1 then

begin

maxb1:= maxb2;

Itek:=i

end;

i:= i+1;

end; {while i<=nstar}

{ dBm:=sqrt(dBm/nstar); }

{––––––––––-конец 2-х вариантов–––––––-}

{====================================================

deltL:=(ymin+ymax)/2; }{- старый расчет оптимальной

подкрутки}

{Уточненный расчет оптимальной подкрутки по долготе:

ищем максимум по C минимума по i величины

cos(B)*[abs(dL(i) - C],

где B - максимум из широты Альмагеста и расчетной широты,

dL(i) - разница между расчетной и альмагестовской

долготой для i-той звезды.

Полученное C дает величину отимальной подкрутки deltL }

x:=0.01;

y:=ymax-ymin;

Nstep:=Trunc(y/x);

cminmax:=ymin;

cc:=ymin;

fminmax:=7;

for i:=1 to Nstep do

begin

cc:=cc+x;

fmax:=0;

for j:=1 to nstar do

begin

fx:=Min(stf[j,3],co[j].cb);

fx:=fx*abs(stf[j,4]-cc);

if fx > fmax then fmax:=fx;

end;

if fmax < fminmax then

begin

fminmax:=fmax;

deltL:=cc;

end;

{====================================================}

{if (maxb1 < bmax) then }

if (dBm < dBmm) then { <- выбран один из трех вариантов }

{ if (NBm > NBmm) then }

begin

bmax:=maxb1;

Iok:=Itek;

Ng0:=ig;

Nb0:=ib;

dBmm:=dBm;

NBmm:=NBm;

for jj:=1 to NBm do

begin

zvv[jj]:=zv[jj];

end;

gr:=0.0;

for i:=1 to nstar do

begin

x:= (stf[i,4]-deltL)*Min(stf[i,3],co[i].cb);

x:= sqr(x);

y:=sqr(stf[i,2] - co[i].b);

x:=sqrt(x+y);

if x > gr then gr:=x;

end;

dist0:=gr;

end; {if maxb1<bmax, if dBm < dBmm или if (NBm > NBmm)}

abett:=angle;

cbett:=cangle;

sbett:=sangle;

angle:=angle+bstep;

cangle:= cbett*cbstep - sbett*sbstep;

sangle:= sbett*cbstep + cbett*sbstep;

end; {for ib - поворот поперек}

angle:= agamt+gstep;

cangle:= cgamt*cgstep - sgamt*sgstep;

sangle:= sgamt*cgstep + cgamt*sgstep;

end; {for ig - поворот вдоль}

{**************************************************************}

{вывод результатов в файл и на печать}

gamma0:= (Ng0*gstep+Gstat[t])*mg;

beta0:=Nb0*bstep*mg;

bmax:=bmax*mg;

dist0:=dist0*mg;

dBmm:=dBmm*mg;

t2:=1900-t*100;

writeln(f1,’==================================================’);

writeln(f1,’Max distance to inf. kernel allowed = ’,d8,’(min)’);

writeln(f1,’ ’,t2:2,’ ’,bmax:4:1,’ (’,co[Iok].nb:4,

’) ’,gamma0:4:1,’ ’,beta0:4:1,’ ’,dist0:4:1);

writeln(f1,’ sigma=’,dBmm:4:1,’ Nstars= ’,NBmm,’(’,eps,’min close)’);

for jj:=1 to NBmm do

begin

jj1:=zvv[jj];

writeln(f1,co[jj1].nb,’ ’,co[jj1].obozn);

end;

writeln(f2,t2:2,’ ’,dBmm:4:1,’ ’,bmax:4:1,’ ’,NBmm);

writeln(’*** T = ’,t2:2,’ ***’);

writeln(’Max distance to inf. kernel allowed = ’,d8,’(min)’);

writeln(’dBmax=’,bmax:4:1,’ i=’,co[Iok].nb:4,’ gamma=’,gamma0:4:1,

’ beta= ’,beta0:4:1,’ dist=’,dist0:4:1);

writeln(’sigma=’,dBmm:4:1,’; Nstars (’,eps,’ min close)=’,NBmm);

end; {for t}

close(f1);

close(f2);

writeln(’Enter any character’);

readln(konec);

end.

===========================================

Входные данные для программы PERESCHET (файл FAST.TXT)

1 столбец – номер звезды в каталоге ярких звезд BS4, BS5;

2 столбец – прямое восхождение RA 1900 по BS5: часы;

3 столбец – прямое восхождение RA 1900 по BS5: минуты;

4 столбец – прямое восхождение RA 1900 по BS5: секунды;

5 столбец – склонение DEC 1900 по BS5: градусы;

6 столбец – склонение DEC 1900 по BS5: минуты;

7 столбец – склонение DEC 1900 по BS5: секунды;

8 столбец – звездная величина по BS5;

9 столбец – скорость собственного движения в RA1900,

приведенная к экватору (по BS4);

10 столбец – скорость собственного движения в DEC1900,

приведенная к экватору (по BS4);

11 столбец – долгота в Альмагесте;

12 столбец – широта в Альмагесте;

13 столбец – яркость в Альмагесте;

14 столбец – современное название звезды по BS5.

1. Файл данных: 8 звезд информативного ядра Альмагеста.

5340 14 11 06.0 +19 42 11 -0.04 -1.098 -1.999 177 00 +31 30 1. 16Alp Boo

1708 05 09 18.0 +45 53 47 0.08 +0.080 -0.423 55 00 +22 30 1. 13Alp Aur

3982 10 03 02.8 +12 27 22 1.35 -0.249 +0.003 122 30 0 10 1. 32Alp Leo

2943 07 34 04.0 +05 28 53 0.38 -0.706 -1.029 89 10 -16 10 1. 10Alp CMi

5056 13 19 55.4 -10 38 22 0.98 -0.043 -0.033 176 40 -2 0 1. 67Alp Vir

6134 16 23 16.4 -26 12 36 0.96 -0.007 -0.023 222 40 -4 0 2. 21Alp Sco

7001 18 33 33.1 +38 41 26 0.03 +0.200 +0.285 257 20 62 0 1. 3Alp Lyr

3449 08 37 29.9 +21 49 42 4.66 -0.103 -0.043 100 20 2 40 3.7 43Gam Cnc

2. Файл данных: именные звезды из A, ZodA, B, ZodB, M,

быстрые (>=0.1"/год по RA1900 или DEC1900) и изолированные на небе

среди звезд сравнимой с ними яркости, что приводит к однозначности их

отождествления в каталоге альмагеста. В начало списка добавлено

информативное ядро Альмагеста из 8-ми звезд.