/**
* @covers Marando\Meeus\Nutation\Nutation::trueObliquity
* @see p.148, ex. 22.a
*/
public function testTrueObliquity()
{
$ε = Nutation::trueObliquity(AstroDate::parse('1987-Apr-10'));
$this->assertEquals(23, $ε->d, 'd');
$this->assertEquals(26, $ε->m, 'm');
$this->assertEquals(36.85, $ε->s, 's', 0.01);
}
/**
* Creates a new Cartesian vector instance
*
* @param Frame $frame Reference frame
* @param Epoch|AstroDate $epoch Observation epoch
* @param Distance $x x position
* @param Distance $y y position
* @param Distance $z z position
* @param Velocity $vx x velocity
* @param Velocity $vy y velocity
* @param Velocity $vz z velocity
*/
public function __construct(Frame $frame, $epoch, Distance $x, Distance $y, Distance $z, Velocity $vx = null, Velocity $vy = null, Velocity $vz = null)
{
// Set reference frame and observation epoch
$this->frame = $frame;
$this->epoch = $epoch instanceof Epoch ? $epoch : $epoch->toEpoch();
// Set position components
$this->x = $x;
$this->y = $y;
$this->z = $z;
// Set velocity components
$this->vx = $vx;
$this->vy = $vy;
$this->vz = $vz;
// Set default rounding and units
$this->decimalPlaces = 13;
}
/**
* @covers Marando\AstroCoord\Frame::FK4
*/
public function testFK4()
{
$equinoxes = [Epoch::B1950(), Epoch::B1900(), Epoch::B(1948), Epoch::dt(AstroDate::now())];
foreach ($equinoxes as $e) {
$equinox = $e;
$fk4 = Frame::FK4($e);
$this->assertEquals('FK4', $fk4->name);
$this->assertEquals($equinox, $fk4->equinox);
}
}
public function testToEclip()
{
// Earth -> Venus @ 2015-Nov-08 23:20:34.000 UT
$dt = AstroDate::jd(2457335.47261574);
$x = Distance::au(-0.7956853147170494);
$y = Distance::au(-0.00807301690301796);
$z = Distance::au(0.01392567642390632);
$c = new Cartesian(Frame::ICRF(), $dt->toEpoch(), $x, $y, $z);
$equat = $c->toEquat();
$eclip = $equat->toEclip();
$prec = Angle::arcmin(13)->deg;
$this->assertEquals(180.3510658, $eclip->lon->deg, 'λ', $prec);
$this->assertEquals(1.1487465, $eclip->lat->deg, 'β', $prec);
}
/**
* Calculates the nutations in longitude (Δψ) and latitude (Δε) for a date
*
* @param AstroDate $date
* @return static
*
* @see Meeus, Jean. Astronomical Algorithms. Richmond, Virg.: Willmann-Bell,
* 2009. 143-147. Print.
*/
public static function find(AstroDate $date)
{
// Time factor
$t = ($date->copy()->toUTC()->jd - 2451545.0) / 36525;
/*
* Required terms
* ----------------
* D = mean elongation of the Moon from the Sun
* M = mean anomaly of the Earth
* M´ = mean anomaly of the Moon
* F = argument of latitude for Moon
* Ω = long. of asc. node for Moon's mean ecl. orbit (mean equi. of date)
*/
$coeffD = [297.85036, 445267.1148, -0.0019142, 1 / 189474];
$coeffM = [357.52772, 35999.05034, -0.0001603, -1 / 300000];
$coeffM´ = [134.96298, 477198.867398, 0.0086972, 1 / 5620];
$coeffF = [93.27191000000001, 483202.017538, -0.0036825, 1 / 327270];
$coeffΩ = [125.04452, -1934.136261, 0.0020708, 1.0 / 450000];
// Calculate the terms
$D = Angle::deg(static::Horner($t, $coeffD))->norm();
$M = Angle::deg(static::Horner($t, $coeffM))->norm();
$M´ = Angle::deg(static::Horner($t, $coeffM´))->norm();
$F = Angle::deg(static::Horner($t, $coeffF))->norm();
$Ω = Angle::deg(static::Horner($t, $coeffΩ))->norm();
// Nutation coefficient terms
$nutationTerms = static::NutationTerms();
// Evaluate the nutation terms
$Δψ = 0;
$Δε = 0;
for ($i = 0; $i < count($nutationTerms); $i++) {
$row = $nutationTerms[$i];
$arg = 0 + $row[0] * $D->rad + $row[1] * $M->rad + $row[2] * $M´->rad + $row[3] * $F->rad + $row[4] * $Ω->rad;
$Δψ += ($row[5] + $row[6] * $t) * sin($arg) / 10000.0 / Time::SEC_IN_HOUR;
$Δε += ($row[7] + $row[8] * $t) * cos($arg) / 10000.0 / Time::SEC_IN_HOUR;
}
// Store as angles
$Δψ = Angle::deg($Δψ);
$Δε = Angle::deg($Δε);
// Return the nutation
return new Nutation($Δψ, $Δε);
}
/**
* @covers Marando\AstroDate\AstroDate::equinoxAutumn
*/
public function testEquinoxAutumn()
{
$sep = [[1996, 22, 18, 1, 8], [1997, 22, 23, 56, 49], [1998, 23, 5, 38, 15], [1999, 23, 11, 32, 34], [2000, 22, 17, 28, 40], [2001, 22, 23, 5, 32], [2002, 23, 4, 56, 28], [2003, 23, 10, 47, 53], [2004, 22, 16, 30, 54], [2005, 22, 22, 24, 14]];
foreach ($sep as $e) {
$equinox = AstroDate::equinoxAutumn($e[0]);
$expected = new AstroDate($e[0], 9, $e[1], $e[2], $e[3], $e[4]);
$this->assertEquals($expected->toJD(), $equinox->toJD(), $e[0], 0.001);
}
}
public function test()
{
return;
$frame = Frame::ICRF();
$epoch = Epoch::J(2000);
echo $epoch->toDate()->toTDB();
echo "\n";
echo "\n";
exit;
$ra = \Marando\Units\Time::hms(11, 16, 46.6);
$dec = \Marando\Units\Angle::dms(5, 45, 32.5);
$dist = Distance::au(5.8);
echo "\n" . ($eq = new Equat($frame, $epoch, $ra, $dec, $dist));
echo "\n" . ($eq = (new Equat($frame, $epoch, $ra, $dec, $dist))->apparent());
echo "\n";
return;
$frame = Frame::ICRF();
$epoch = AstroDate::jd(2457335.472615741)->toEpoch();
$ra = \Marando\Units\Time::hms(11, 16, 46.6);
$dec = \Marando\Units\Angle::dms(5, 45, 32.5);
$dist = Distance::au(5.8);
echo "\n" . ($eq = new Equat($frame, $epoch, $ra, $dec, $dist));
echo "\n" . ($eq = (new Equat($frame, $epoch, $ra, $dec, $dist))->apparent());
$eq = new Equat($frame, $epoch, $ra, $dec, $dist);
$eq->topo = Geo::deg(27, -82);
echo "\n" . $eq->toHoriz();
echo "\n" . $eq->apparent();
echo "\n" . ($h = $eq->toHoriz());
echo "\n" . $h->az->deg . "\t" . $h->alt->deg;
echo "\n" . $eq->apparent(\Marando\Units\Pressure::mbar(100), \Marando\Units\Temperature::F(90), 0.85);
echo "\n" . $eq->toHoriz();
return;
$ra = \Marando\Units\Time::hours(14.424354);
$dec = \Marando\Units\Angle::deg(33.54366);
$dist = Distance::au(1.5);
$geo = Geo::deg(27, -82);
$e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist);
$e->topo = $geo;
echo "\n\n" . $e;
echo "\n" . $e->apparent();
$e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist);
$e->topo = $geo;
echo "\n" . $e->apparent();
$e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist);
$e->topo = $geo;
echo "\n" . $e->toHoriz();
return;
$e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist);
echo "\n\n" . $e;
echo "\n" . $e->apparent($geo);
echo "\n" . $e->apparent()->toHoriz($geo);
$e = new Equat(Frame::ICRF(), Epoch::jd(2455586), $ra, $dec, $dist);
echo "\n\n" . $e;
echo "\n" . $e->apparent();
echo "\n" . $e->apparent()->toHoriz();
$e = new Equat(Frame::ICRF(), Epoch::jd(2451586), $ra, $dec, $dist);
$e->topo = $geo;
echo "\n\n" . $e;
echo "\n" . $e->apparent();
echo "\n" . $e->apparent()->toHoriz();
return;
// Position of Mercury
$x = Distance::au(+1.18052679326447);
$y = Distance::au(-0.3650485652522116);
$z = Distance::au(-0.2123422968928603);
// Velocity of Mercury
$vx = Velocity::aud(0.02052320553396918);
$vy = Velocity::aud(0.02828864020900079);
$vz = Velocity::aud(0.01145246106968459);
// Frame & epoch
$frame = Frame::ICRF();
$epoch = AstroDate::parse('2015-Mar-20')->toEpoch();
$c = new Cartesian($frame, $epoch, $x, $y, $z, $vx, $vy, $vz);
echo "\n\n";
echo $c;
echo "\n\n";
echo $c->setUnit('km km/d');
echo "\n\n";
echo $c->setUnit('km km/s');
echo "\n\n";
echo $c->toEquat();
echo "\n\n";
echo $c->toEquat()->apparent();
echo "\n\n";
echo $c->toEquat()->apparent(Geo::deg(27, 278));
echo "\n\n";
echo $c->toEquat()->toHoriz(Geo::deg(27, 278));
echo "\n\n";
}
public function test()
{
return;
// Converts to J2000.0
echo "\n" . AstroDate::parse('2000-Jan-1 12:00:00 TT')->toEpoch();
// Isn't J2000.0 so represented as date
echo "\n" . ($e = AstroDate::parse('2000-Jan-1 14:00:00 TT')->toEpoch());
echo "\n\n" . AstroDate::now()->toEpoch();
$d = AstroDate::now()->setTimezone('est');
echo "\n\n" . $d->sinceMidnight();
echo "\n\n";
return;
echo "\n" . ($d = AstroDate::now()->setTimezone(TimeZone::UT(-7)));
echo "\n" . TimeScale::TT()->abr;
echo "\n" . TimeScale::TT()->name;
echo "\n" . Epoch::J2000();
echo "\n" . Epoch::J1900();
echo "\n" . Epoch::B1950();
echo "\n" . Epoch::B1900();
echo "\n" . Epoch::B(1954.4423);
echo "\n" . Epoch::jd(2455200.5);
echo "\n\n" . Epoch::J2000()->toDate();
echo "\n" . Epoch::J1900()->toDate();
echo "\n" . Epoch::B1950()->toDate();
echo "\n" . Epoch::B1900()->toDate();
echo "\n" . Epoch::B(1954.4423)->toDate();
echo "\n\n" . AstroDate::now();
echo "\n" . ($d = AstroDate::now()->setTimezone(TimeZone::parse('EST')));
echo "\n" . $d->toTDB();
echo "\n" . $d->toUTC();
echo "\n" . ($d = AstroDate::now()->setTimezone('MST'));
echo "\n" . ($d = AstroDate::now()->setTimezone('PST'));
return;
echo "\n" . AstroDate::parse('2015-Dec-10 6:00');
echo "\n" . AstroDate::parse('2015-Dec-10 6:00')->setTimezone(TimeZone::parse('EST'));
echo "\n" . AstroDate::now();
echo "\n" . AstroDate::now()->setTimezone(TimeZone::parse('EST'));
echo "\n" . AstroDate::now()->setTimezone(TimeZone::parse('PST'));
echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(6));
echo "\n" . TimeZone::parse('est')->offset;
echo "\n" . TimeZone::parse('est')->offset(2451545.5);
echo "\n" . TimeZone::parse('est')->offset(2451589.5);
return;
echo "\n" . TimeZone::UT(-2);
echo "\n" . TimeZone::parse('EST');
echo "\n" . AstroDate::now()->setTimezone(TimeZone::EST())->format(AstroDate::FORMAT_GOOGLE);
echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(-10.5))->format(AstroDate::FORMAT_GOOGLE);
echo "\n" . AstroDate::now()->setTimezone(TimeZone::UT(-10))->format(AstroDate::FORMAT_GOOGLE);
return;
$t = timezone_abbreviations_list();
$s = timezone_identifiers_list();
$e = timezone_name_from_abbr('est');
var_dump($t);
return;
$d = AstroDate::now()->setTimezone(TimeZone::EST());
echo "\n" . $d->format('r Y-M-c h:i:s.u A T');
echo "\n" . $d->sub(Time::days(15))->format(AstroDate::FORMAT_GENERIC);
echo "\n" . $d->sidereal('a');
echo "\n" . $d->sidereal('m');
echo "\n" . $d->format(AstroDate::FORMAT_JPL_FRAC);
echo "\n" . $d->sidereal('a', Angle::deg(-82.47));
echo "\n" . $d->sidereal('m', Angle::deg(-82.47));
echo "\n" . $d->sinceMidnight();
echo "\n" . $d->untilMidnight();
/**
* TODO:
*
* - Figure out formula for DST start and end
* - Figure out formula for Week # of year
*
*
*/
return;
echo "\n" . $d->format(DateTime::RSS);
echo "\n" . $d;
echo "\n" . $d->format(AstroDate::FORMAT_GENERIC);
echo "\n" . $d->format(DateTime::ISO8601);
return;
$str = '2016-Nov-14 17:07:07.120';
echo "\n" . ($d = AstroDate::parse($str));
echo "\n" . ($d = AstroDate::parse($str)->setTimezone(TimeZone::EST()));
echo "\n" . $d->format('Z T P O e u s i h H G g a A y Y n M m F W z D Y-m-d j');
$str = '2016-Nov-14 17:07:07.120 TAI';
echo "\n" . ($d = AstroDate::parse($str));
echo "\n" . $d->format('O e u s i h H G g a A y Y n M m F W z D Y-m-d j');
$str = '2015-Nov-1 17:07:07.120 UTC';
$d = AstroDate::parse($str);
echo "\n" . $d->format('z w S N l L D Y-m-d j');
return;
echo "\n\n" . ($str = '2015-Nov-16 17:07:07.120 TT');
echo "\n" . AstroDate::parse($str);
echo "\n\n" . ($str = '-1950-1-16 17:07:07 UTC');
echo "\n" . ($d = AstroDate::parse($str));
echo "\n" . $d->format('Y-m-d');
return;
echo "\n" . ($d = AstroDate::now());
echo "\n" . $d->toUT1();
echo "\n" . $d->toTAI();
echo "\n\n" . ($d = AstroDate::now(TimeZone::EST()));
echo "\n" . $d->toUT1();
echo "\n" . $d->toTAI();
return;
$d = new AstroDate(2017, 11, 15, 7, 0, 0);
echo "\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::EST());
//var_dump($d);
$d = new AstroDate(2017, 6, 15, 7, 0, 0);
echo "\n\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::EST());
//var_dump($d);
$d = new AstroDate(2017, 11, 15, 7, 0, 0);
echo "\n\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::EST());
//var_dump($d);
$d = new AstroDate(2017, 6, 15, 6, 0, 0);
echo "\n\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::EST());
//var_dump($d);
return;
$d = new AstroDate(2017, 11, 15, 0, 0, 0);
echo "\n" . $d;
echo "\n" . $d->jd();
echo "\n" . $d->jd(12);
echo "\n" . $d->toTT();
echo "\n" . $d->toJD();
echo "\n" . $d->jd(12);
echo "\n";
$d = new AstroDate(2015, 11, 15, 20, 23, 18, TimeZone::EST(), TimeScale::TT());
echo "\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::UTC());
echo "\n" . $d->toUT1();
echo "\n" . $d->setTimezone(TimeZone::EST());
echo "\n" . $d->jd();
echo "\n" . $d->jd(12);
echo "\n" . $d->mjd();
echo "\n" . $d->mjd(12);
echo "\n" . $d->monthName();
echo "\n" . $d->mjd();
echo "\n\n" . $d->toTAI();
echo "\n" . $d->year;
echo "\n" . $d->month;
echo "\n" . $d->day;
echo "\n" . $d->hour;
echo "\n" . $d->min;
echo "\n" . $d->sec;
echo "\n" . $d->micro;
echo "\n" . $d->timezone;
echo "\n" . $d->timescale;
$d->year = 2020;
echo "\n" . $d;
$d->month = 12.13432;
echo "\n" . $d;
$d->sec = 12;
echo "\n" . $d;
echo "\n" . $d->add(Time::sec(0.5));
echo "\n" . $d->setDateTime(2017, 9, 4, 18, 34, 34.234);
var_dump($d->isLeapYear());
echo "\n" . $d->setDateTime(2016, 9, 4, 18, 34, 34.234);
var_dump($d->isLeapYear());
echo "\n" . $d->dayName();
$a = new AstroDate(2016, 12, 25, 19, 10, 2);
$b = new AstroDate(2015, 12, 25, 19, 10, 2);
echo "\n" . $a->diff($b);
echo "\n" . ($a = new AstroDate(2016, 12, 31, 19, 10, 2));
echo "\n" . ($b = new AstroDate(2015, 12, 31, 19, 10, 2));
echo "\n" . $a->dayOfYear();
echo "\n" . $b->dayOfYear();
echo "\n\n" . $b;
echo "\n" . $b->toTDB();
echo "\n" . $b->setTimezone(TimeZone::UTC());
echo "\n" . $b->toTDB();
echo "\n" . $b->setTimezone(TimeZone::EST());
echo "\n" . $b->toTDB();
echo "\n\n\n";
$d = new AstroDate(2015, 11, 15, 20, 23, 18.454334);
echo "\n" . $d;
echo "\n" . $d->setTimezone(TimeZone::EST());
echo "\n" . $d->toUTC();
echo "\n" . $d->toTAI();
echo "\n" . $d->setTimezone(TimeZone::EST());
echo "\n" . $d->toTT();
return;
$d = new AstroDate(2015, 11, 15, 20, 23, 18.454334);
echo "\n" . $d;
$d = new AstroDate(2015, 11, 15, 20, 23, 18.454334, TimeZone::EST());
echo "\n" . $d;
echo "\n" . $d->setDate(2020, 10, 2);
echo "\n" . $d->setTime(23, 59, 1);
echo "\n" . $d->setTimezone(TimeZone::UTC());
echo "\n" . $d->setTimezone(TimeZone::EST());
}
public function testCreateWithAstroDate()
{
$x = Distance::au(+9.32);
$y = Distance::au(-5.25);
$z = Distance::au(-6.89);
$c = new Cartesian(Frame::ICRF(), AstroDate::now(), $x, $y, $z);
}
/**
* Finds the difference between this AstroDate and another
* @param AstroDate $b
* @return static
*/
public function diff(AstroDate $b)
{
$prec = 12;
$jd1 = $this->toJD($prec);
$jd2 = $b->toJD($prec);
$days = bcsub($jd1, $jd2, $prec);
return Time::days(-1 * $days);
}
