📄 sturanus.pas
字号:
L4 := 0.00000113855 * cos(3.14159265360 + 0.00000000000 * Tau)
+ 0.00000005599 * cos(4.57882424420 + 74.78159856700 * Tau)
+ 0.00000003203 * cos(0.34623003207 + 11.04570026400 * Tau)
+ 0.00000001217 * cos(3.42199121830 + 56.62235130300 * Tau);
L5 := 0.00000000000;
Result := (L0 + L1*Tau + L2*Tau2 + L3*Tau3 + L4*Tau4 + L5*Tau5);
end;
{---------------------------------------------------------------------------}
function GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
var
B0, B1,
B2, B3,
B4, B5 : Double;
begin
B0 := 0.01346277639 * cos(2.61877810550 + 74.78159856700 * Tau)
+ 0.00062341405 * cos(5.08111175860 + 149.56319713000 * Tau)
+ 0.00061601203 * cos(3.14159265360 + 0.00000000000 * Tau)
+ 0.00009963744 * cos(1.61603876360 + 76.26607127600 * Tau)
+ 0.00009926151 * cos(0.57630387917 + 73.29712585900 * Tau)
+ 0.00003259455 * cos(1.26119385960 + 224.34479570000 * Tau)
+ 0.00002972318 * cos(2.24367035540 + 1.48447270830 * Tau)
+ 0.00002010257 * cos(6.05550401090 + 148.07872443000 * Tau)
+ 0.00001522172 * cos(0.27960386377 + 63.73589830300 * Tau)
+ 0.00000924055 * cos(4.03822927850 + 151.04766984000 * Tau)
+ 0.00000760624 * cos(6.14000431920 + 71.81265315100 * Tau)
+ 0.00000522309 * cos(3.32085194770 + 138.51749687000 * Tau)
+ 0.00000462630 * cos(0.74256727574 + 85.82729883100 * Tau)
+ 0.00000436843 * cos(3.38082524320 + 529.69096509000 * Tau)
+ 0.00000434625 * cos(0.34065281858 + 77.75054398400 * Tau)
+ 0.00000430668 * cos(3.55445034850 + 213.29909544000 * Tau)
+ 0.00000420265 * cos(5.21279984790 + 11.04570026400 * Tau)
+ 0.00000244698 * cos(0.78795150326 + 2.96894541660 * Tau)
+ 0.00000232649 * cos(2.25716421380 + 222.86032299000 * Tau)
+ 0.00000215838 * cos(1.59121704940 + 38.13303563800 * Tau)
+ 0.00000179935 * cos(3.72487952670 + 299.12639427000 * Tau)
+ 0.00000174895 * cos(1.23550262210 + 146.59425172000 * Tau)
+ 0.00000173667 * cos(1.93654269130 + 380.12776796000 * Tau)
+ 0.00000160368 * cos(5.33635436460 + 111.43016150000 * Tau)
+ 0.00000144064 * cos(5.96239326410 + 35.16409022100 * Tau)
+ 0.00000116363 * cos(5.73877190010 + 70.84944530400 * Tau)
+ 0.00000106441 * cos(0.94103112994 + 70.32818044200 * Tau)
+ 0.00000102049 * cos(2.61876256510 + 78.71375183000 * Tau);
B1 := 0.00206366162 * cos(4.12394311410 + 74.78159856700 * Tau)
+ 0.00008563230 * cos(0.33819986165 + 149.56319713000 * Tau)
+ 0.00001725703 * cos(2.12193159900 + 73.29712585900 * Tau)
+ 0.00001374449 * cos(0.00000000000 + 0.00000000000 * Tau)
+ 0.00001368860 * cos(3.06861722050 + 76.26607127600 * Tau)
+ 0.00000450639 * cos(3.77656180980 + 1.48447270830 * Tau)
+ 0.00000399847 * cos(2.84767037790 + 224.34479570000 * Tau)
+ 0.00000307214 * cos(1.25456766740 + 148.07872443000 * Tau)
+ 0.00000154336 * cos(3.78575467750 + 63.73589830300 * Tau)
+ 0.00000112432 * cos(5.57299891500 + 151.04766984000 * Tau)
+ 0.00000110888 * cos(5.32888676460 + 138.51749687000 * Tau)
+ 0.00000083493 * cos(3.59152795560 + 71.81265315100 * Tau)
+ 0.00000055573 * cos(3.40135416350 + 85.82729883100 * Tau)
+ 0.00000053690 * cos(1.70455769940 + 77.75054398400 * Tau)
+ 0.00000041912 * cos(1.21476607430 + 11.04570026400 * Tau)
+ 0.00000041377 * cos(4.45476669140 + 78.71375183000 * Tau)
+ 0.00000031959 * cos(3.77446207750 + 222.86032299000 * Tau)
+ 0.00000030297 * cos(2.56371683640 + 2.96894541660 * Tau)
+ 0.00000026977 * cos(5.33695500290 + 213.29909544000 * Tau)
+ 0.00000026222 * cos(0.41620628369 + 380.12776796000 * Tau);
B2 := 0.00009211656 * cos(5.80044305790 + 74.78159856700 * Tau)
+ 0.00000556926 * cos(0.00000000000 + 0.00000000000 * Tau)
+ 0.00000286265 * cos(2.17729776350 + 149.56319713000 * Tau)
+ 0.00000094969 * cos(3.84237569810 + 73.29712585900 * Tau)
+ 0.00000045419 * cos(4.87822046060 + 76.26607127600 * Tau)
+ 0.00000020107 * cos(5.46264485370 + 1.48447270830 * Tau)
+ 0.00000014793 * cos(0.87983715652 + 138.51749687000 * Tau)
+ 0.00000014261 * cos(2.84517742690 + 148.07872443000 * Tau)
+ 0.00000013963 * cos(5.07234043990 + 63.73589830300 * Tau)
+ 0.00000010122 * cos(5.00290894860 + 224.34479570000 * Tau)
+ 0.00000008299 * cos(6.26655615200 + 78.71375183000 * Tau);
B3 := 0.00000267832 * cos(1.25097888290 + 74.78159856700 * Tau)
+ 0.00000011048 * cos(3.14159265360 + 0.00000000000 * Tau)
+ 0.00000006154 * cos(4.00663614490 + 149.56319713000 * Tau)
+ 0.00000003361 * cos(5.77804694940 + 73.29712585900 * Tau);
B4 := 0.00000005719 * cos(2.85499529310 + 74.78159856700 * Tau);
B5 := 0.00000000000;
Result := (B0 + B1*Tau + B2*Tau2 + B3*Tau3 + B4*Tau4 + B5*Tau5);
end;
{---------------------------------------------------------------------------}
function GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5 : Double) : Double;
var
R0, R1,
R2, R3,
R4, R5 : Double;
begin
R0 := 19.21264847900 * cos(0.00000000000 + 0.00000000000 * Tau)
+ 0.88784984055 * cos(5.60377526990 + 74.78159856700 * Tau)
+ 0.03440835545 * cos(0.32836098991 + 73.29712585900 * Tau)
+ 0.02055653495 * cos(1.78295170030 + 149.56319713000 * Tau)
+ 0.00649321851 * cos(4.52247298120 + 76.26607127600 * Tau)
+ 0.00602248144 * cos(3.86003820460 + 63.73589830300 * Tau)
+ 0.00496404171 * cos(1.40139934720 + 454.90936653000 * Tau)
+ 0.00338525522 * cos(1.58002682950 + 138.51749687000 * Tau)
+ 0.00243508222 * cos(1.57086595070 + 71.81265315100 * Tau)
+ 0.00190521915 * cos(1.99809364500 + 1.48447270830 * Tau)
+ 0.00161858251 * cos(2.79137863470 + 148.07872443000 * Tau)
+ 0.00143705902 * cos(1.38368574480 + 11.04570026400 * Tau)
+ 0.00093192359 * cos(0.17437193645 + 36.64856292900 * Tau)
+ 0.00089805842 * cos(3.66105366330 + 109.94568879000 * Tau)
+ 0.00071424265 * cos(4.24509327400 + 224.34479570000 * Tau)
+ 0.00046677322 * cos(1.39976563940 + 35.16409022100 * Tau)
+ 0.00039025681 * cos(3.36234710690 + 277.03499374000 * Tau)
+ 0.00039009624 * cos(1.66971128870 + 70.84944530400 * Tau)
+ 0.00036755160 * cos(3.88648934740 + 146.59425172000 * Tau)
+ 0.00030348875 * cos(0.70100446346 + 151.04766984000 * Tau)
+ 0.00029156264 * cos(3.18056174560 + 77.75054398400 * Tau)
+ 0.00025785805 * cos(3.78537741500 + 85.82729883100 * Tau)
+ 0.00025620360 * cos(5.25656292800 + 380.12776796000 * Tau)
+ 0.00022637152 * cos(0.72519137745 + 529.69096509000 * Tau)
+ 0.00020473163 * cos(2.79639811630 + 70.32818044200 * Tau)
+ 0.00020471584 * cos(1.55588961500 + 202.25339517000 * Tau)
+ 0.00017900561 * cos(0.55455488605 + 2.96894541660 * Tau)
+ 0.00015502809 * cos(5.35405037600 + 38.13303563800 * Tau)
+ 0.00014701566 * cos(4.90434406650 + 108.46121608000 * Tau)
+ 0.00012896507 * cos(2.62154018240 + 111.43016150000 * Tau)
+ 0.00012328151 * cos(5.96039150920 + 127.47179661000 * Tau)
+ 0.00011959355 * cos(1.75044072170 + 984.60033162000 * Tau)
+ 0.00011852996 * cos(0.99342814582 + 52.69019803900 * Tau)
+ 0.00011696085 * cos(3.29825599110 + 3.93215326310 * Tau)
+ 0.00011494701 * cos(0.43774027872 + 65.22037101200 * Tau)
+ 0.00010792699 * cos(1.42104858470 + 213.29909544000 * Tau)
+ 0.00009111446 * cos(4.99638600050 + 62.25142559500 * Tau)
+ 0.00008420550 * cos(5.25350716620 + 222.86032299000 * Tau)
+ 0.00008402147 * cos(5.03877516490 + 415.55249061000 * Tau)
+ 0.00007449125 * cos(0.79491905956 + 351.81659231000 * Tau)
+ 0.00007329454 * cos(3.97277527840 + 183.24281465000 * Tau)
+ 0.00006046370 * cos(5.67960948360 + 78.71375183000 * Tau)
+ 0.00005524133 * cos(3.11499484160 + 9.56122755560 * Tau)
+ 0.00005444878 * cos(5.10575635360 + 145.10977901000 * Tau)
+ 0.00005238103 * cos(2.62960141800 + 33.67961751300 * Tau)
+ 0.00004079167 * cos(3.22064788670 + 340.77089205000 * Tau)
+ 0.00003919476 * cos(4.25015288870 + 39.61750834600 * Tau)
+ 0.00003801606 * cos(6.10985558500 + 184.72728736000 * Tau)
+ 0.00003781219 * cos(3.45840272870 + 456.39383924000 * Tau)
+ 0.00003686787 * cos(2.48718116540 + 453.42489382000 * Tau)
+ 0.00003101743 * cos(4.14031063900 + 219.89137758000 * Tau)
+ 0.00002962641 * cos(0.82977991995 + 56.62235130300 * Tau)
+ 0.00002942239 * cos(0.42393808854 + 299.12639427000 * Tau)
+ 0.00002940492 * cos(2.14637460320 + 137.03302416000 * Tau)
+ 0.00002937799 * cos(3.67657450930 + 140.00196958000 * Tau)
+ 0.00002865128 * cos(0.30996903761 + 12.53017297200 * Tau)
+ 0.00002538032 * cos(4.85457831990 + 131.40394987000 * Tau)
+ 0.00002363550 * cos(0.44253328372 + 554.06998748000 * Tau)
+ 0.00002182572 * cos(2.94040431640 + 305.34616939000 * Tau);
R1 := 0.01479896370 * cos(3.67205705320 + 74.78159856700 * Tau)
+ 0.00071212085 * cos(6.22601006670 + 63.73589830300 * Tau)
+ 0.00068626972 * cos(6.13411265050 + 149.56319713000 * Tau)
+ 0.00024059649 * cos(3.14159265360 + 0.00000000000 * Tau)
+ 0.00021468152 * cos(2.60176704270 + 76.26607127600 * Tau)
+ 0.00020857262 * cos(5.24625494220 + 11.04570026400 * Tau)
+ 0.00011405346 * cos(0.01848461561 + 70.84944530400 * Tau)
+ 0.00007496775 * cos(0.42360033283 + 73.29712585900 * Tau)
+ 0.00004243800 * cos(1.41692350370 + 85.82729883100 * Tau)
+ 0.00003926694 * cos(3.15513991320 + 71.81265315100 * Tau)
+ 0.00003578446 * cos(2.31160668310 + 224.34479570000 * Tau)
+ 0.00003505936 * cos(2.58354048850 + 138.51749687000 * Tau)
+ 0.00003228835 * cos(5.25499602900 + 3.93215326310 * Tau)
+ 0.00003060010 * cos(0.15321893225 + 1.48447270830 * Tau)
+ 0.00002564251 * cos(0.98076846352 + 148.07872443000 * Tau)
+ 0.00002429445 * cos(3.99440122470 + 52.69019803900 * Tau)
+ 0.00001644719 * cos(2.65349313120 + 127.47179661000 * Tau)
+ 0.00001583766 * cos(1.43045619200 + 78.71375183000 * Tau)
+ 0.00001508028 * cos(5.05996325430 + 151.04766984000 * Tau)
+ 0.00001489525 * cos(2.67559167320 + 56.62235130300 * Tau)
+ 0.00001413112 * cos(4.57461892060 + 202.25339517000 * Tau)
+ 0.00001403237 * cos(1.36985349740 + 77.75054398400 * Tau)
+ 0.00001228220 * cos(1.04703640150 + 62.25142559500 * Tau)
+ 0.00001032731 * cos(0.26459059027 + 131.40394987000 * Tau)
+ 0.00000992085 * cos(2.17168865910 + 65.22037101200 * Tau)
+ 0.00000861867 * cos(5.05530802220 + 351.81659231000 * Tau)
+ 0.00000744445 * cos(3.07640148940 + 35.16409022100 * Tau)
+ 0.00000687470 * cos(2.49912565670 + 77.96299230500 * Tau)
+ 0.00000646851 * cos(4.47290422910 + 70.32818044200 * Tau)
+ 0.00000623602 * cos(0.86253073820 + 9.56122755560 * Tau)
+ 0.00000604362 * cos(0.90717667985 + 984.60033162000 * Tau)
+ 0.00000574710 * cos(3.23070708460 + 447.79581953000 * Tau)
+ 0.00000561839 * cos(2.71778158980 + 462.02291353000 * Tau)
+ 0.00000530364 * cos(5.91655309050 + 213.29909544000 * Tau)
+ 0.00000527794 * cos(5.15136007080 + 2.96894541660 * Tau);
R2 := 0.00022439904 * cos(0.69953118760 + 74.78159856700 * Tau)
+ 0.00004727037 * cos(1.69901641490 + 63.73589830300 * Tau)
+ 0.00001681903 * cos(4.64833551730 + 70.84944530400 * Tau)
+ 0.00001649559 * cos(3.09660078980 + 11.04570026400 * Tau)
+ 0.00001433755 * cos(3.52119917950 + 149.56319713000 * Tau)
+ 0.00000770188 * cos(0.00000000000 + 0.00000000000 * Tau)
+ 0.00000500429 * cos(6.17229032220 + 76.26607127600 * Tau)
+ 0.00000461009 * cos(0.76676632849 + 3.93215326310 * Tau)
+ 0.00000390371 * cos(4.49605283500 + 56.62235130300 * Tau)
+ 0.00000389945 * cos(5.52673426380 + 85.82729883100 * Tau)
+ 0.00000292097 * cos(0.20389012095 + 52.69019803900 * Tau)
+ 0.00000286579 * cos(3.53357683270 + 73.29712585900 * Tau)
+ 0.00000272898 * cos(3.84707823650 + 138.51749687000 * Tau)
+ 0.00000219674 * cos(1.96418942890 + 131.40394987000 * Tau)
+ 0.00000215788 * cos(0.84812474187 + 77.96299230500 * Tau)
+ 0.00000205449 * cos(3.24758017120 + 78.71375183000 * Tau)
+ 0.00000148554 * cos(4.89840863840 + 127.47179661000 * Tau)
+ 0.00000128834 * cos(2.08146849520 + 3.18139373770 * Tau);
R3 := 0.00001164382 * cos(4.73453291600 + 74.78159856700 * Tau)
+ 0.00000212367 * cos(3.34255735000 + 63.73589830300 * Tau)
+ 0.00000196408 * cos(2.98004616320 + 70.84944530400 * Tau)
+ 0.00000104527 * cos(0.95807937648 + 11.04570026400 * Tau)
+ 0.00000072540 * cos(0.99701907912 + 149.56319713000 * Tau)
+ 0.00000071681 * cos(0.02528455665 + 56.62235130300 * Tau)
+ 0.00000054875 * cos(2.59436811270 + 3.93215326310 * Tau)
+ 0.00000036377 * cos(5.65035573020 + 77.96299230500 * Tau)
+ 0.00000034029 * cos(3.81553325640 + 76.26607127600 * Tau)
+ 0.00000032081 * cos(3.59825177840 + 131.40394987000 * Tau);
R4 := 0.00000052996 * cos(3.00838033090 + 74.78159856700 * Tau)
+ 0.00000009887 * cos(1.91399083600 + 56.62235130300 * Tau);
R5 := 0.00000000000;
Result := (R0 + R1*Tau + R2*Tau2 + R3*Tau3 + R4*Tau4 + R5*Tau5);
end;
{---------------------------------------------------------------------------}
function ComputeUranus(JD : Double) : TStEclipticalCord;
var
Tau,
Tau2,
Tau3,
Tau4,
Tau5 : Double;
begin
Tau := (JD - 2451545.0) / 365250.0;
Tau2 := sqr(Tau);
Tau3 := Tau * Tau2;
Tau4 := sqr(Tau2);
Tau5 := Tau2 * Tau3;
Result.L0 := GetLongitude(Tau, Tau2, Tau3, Tau4, Tau5);
Result.B0 := GetLatitude(Tau, Tau2, Tau3, Tau4, Tau5);
Result.R0 := GetRadiusVector(Tau, Tau2, Tau3, Tau4, Tau5);
end;
end.
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -