Subject: О прецессии
Date   : 30 May 1998 16:14 GMT
From   : Anatoly Zaitsev [AZ] (AZ@stel.sebastopol.ua)
To     : All

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

 { коррекция долготы за прецессию }
 { учет собственного движения     }
 t := (dat + utc/24 - 0.5 - 2451545.0)/365.25; { число лет от 2000 }
 x := (50.290966*t+0.000222226*t*t/2)*Rds;     { величина прецессии }
 for i := 1 to sdbm do
 if sdb^[i].lamb <> 0 then begin
   sdb^[i].lamb := sdb^[i].lamb + x + sdb^[i].mL*t*0.01;
   sdb^[i].beta := sdb^[i].beta + sdb^[i].mB*t*0.01;
   { аберрация }
   dl := 20.5*Rds*cos(sdb^[i].lamb - Ls + Pi);
   db := 20.5*Rds*sin(sdb^[i].beta)*sin(sdb^[i].lamb - Ls + Pi);
   { нутация }
   de := ...;
   df := ...;
   sdb^[i].lamb := sdb^[i].lamb + dl + de;
   sdb^[i].beta := sdb^[i].beta + db + df*sin(sdb^[i].lamb);
 end;

Здесь:
dat - юлианский день
utc - всемирное время в часах,
Rds - коэффициент перевода секунд дуги в радианы,
sdbm - размер базы данных звезд,
sdb^[i].lamb - долгота i-ой звезды,
sdb^[i].beta - широта i-ой звезды,
sdb^[i].mL - собственное движение по долготе за 100 лет,
sdb^[i].mB - собственное движение по широте за 100 лет,
Ls - долгота Солнца.