估计太阳8字曲线使用行星的星历表和ECI AER转换
这个例子展示了如何估计8字曲线。8字曲线的曲线代表太阳的角偏移的变化从它的意思是天球上的位置相对于一个特定的地理位置在地球表面。在这个例子中,8字曲线估计相对于在格林威治皇家天文台,联合王国。评估后,图8字曲线的例子。
这个示例使用的数据,您可以下载使用aeroDataPackage函数。
确定每年的日期来计算太阳的8字曲线
指定的日期来计算8字曲线。在这个例子中,这些日期范围从1月1日,2014年到2014年12月31日中午UTC。
1:365 dv = datetime(2014年1日,12日,0,0);dvUTC = [dv。年的dv。月的dv。天的dv。小时的dv。分钟' dv.Second '];
计算太阳的位置
使用planetEphemeris函数来计算太阳的位置。在这个例子中:
的tdbjuliandate函数计算动态重心的朱利安日期时间(TDB)。
的tdbjuliandate函数需要地面时间(TT)。
陆地从UTC时间在秒的计算需要的差异协调世界时(UTC)及国际原子时(TAI)。
2014年,这种差异(dAT)是35秒。
近似地球时间(secTT)是dAT + 32.184秒。
陆地的时间年、月、日、小时、分钟和秒是包含在dvTT数组中。
dAT = 35;secTT = dAT + 32.184;dvTT = dv + secTT / 86400;
估计动态重心的朱利安日期时间基于陆地的时间使用dvTT数组中。
jdTDB = tdbjuliandate ([dvTT。年的dvTT。月的dvTT。天的dvTT。小时的dvTT。分钟' dvTT.Second ']);
这些日期的确定太阳的位置:
posSun = planetEphemeris (jdTDB,“地球”,“太阳”)* 1000;
计算UTC和主要区别世界时(不一致)每
计算UTC的区别和不一致,每deltaUT1,使用修改后的朱利安日期UTC。
mjdUTC = mjuliandate (dvUTC);dUT1 = deltaUT1 (mjdUTC);
计算天体中间杆的极地运动和位移(CIP)
CIP的极地运动和位移计算使用修改后的朱利安日期UTC。
点= polarMotion (mjdUTC);dCIP = deltaCIP (mjdUTC);
指定的几何位置在格林威治皇家天文台,联合王国
指定的位置的几何位置估计8字曲线。在本例中,这个位置是纬度、经度和海拔在格林威治皇家天文台(51.48度、0.0015度,0米高度)。
LLAGreenwich = (51.48, -0.0015, 0);爱尔兰= eci2aer (posSun、dvUTC repmat (LLAGreenwich、长度(jdTDB), 1),…“deltaAT”,dAT *(长度(jdTDB), 1),“deltaUT1”dUT1,…“PolarMotion”点,“dCIP”dCIP);
年内指定天的8字曲线照片你想阴谋
8字曲线,感兴趣的你可以画出天年内8字曲线。这个示例图:
2014年每个月的第一天。
夏至、冬至日有关。
春天和秋天春分和秋分。
获得2014年每个月的第一天:
aerFirstMonth = aer (dvUTC (:, 3) = = 1:);
让冬至、春分和秋分(2014年3/20,6/21,9/22,12/21):
solsticeEquinox =[爱尔兰(dvUTC (:, 2) = = 3 & dvUTC (:, 3) = = 20, 1:2);…爱尔兰(dvUTC (:, 2) = = 6 & dvUTC(:, 3) = = 21日1:2);…爱尔兰(dvUTC (:, 2) = = 9 & dvUTC(:, 3) = = 22日1:2);…爱尔兰(dvUTC (:, 2) = = 12 & dvUTC(:, 3) = = 21日1:2)];
阴谋的结果
图8字曲线。全年8字曲线,情节点,每月的第一天,秋分),二至点。
firstDays = (12: 1:1) +“/”+ 1;f =图;情节(aer(: 1)、爱尔兰(:,2),“。”,…solsticeEquinox (: 1), solsticeEquinox (:, 2),“ks”,…aerFirstMonth (: 1), aerFirstMonth (:, 2),“柯”,…“MarkerSize”8“MarkerFaceColor”,“k”);标题(的8字曲线观察到格林威治天文台);包含(的方位(度));ylabel (的海拔(度));轴([176185、10、70])文本(aerFirstMonth (: 1) +。1,aerFirstMonth (:, 2) + 1.2, firstDays,“颜色”,“k”,“HorizontalAlignment”,“左”);文本(solsticeEquinox (1,1) +。2,solsticeEquinox -1.5 (1、2),“春分”,“颜色”,“k”,“HorizontalAlignment”,“左”);文本(solsticeEquinox (2, 1), solsticeEquinox (2, 2) + 2.5,“夏至”,“颜色”,“k”,“HorizontalAlignment”,“左”);文本(solsticeEquinox (3,1) +。1,solsticeEquinox (3 2) + 1.2,“秋分”,“颜色”,“k”,“HorizontalAlignment”,“左”);文本(solsticeEquinox (4,1) +。1,solsticeEquinox (4,2) -2.5,“冬至”,“颜色”,“k”,“HorizontalAlignment”,“左”);
另请参阅
planetEphemeris
|aeroDataPackage
|deltaCIP
|爱尔兰
|eci2aer
|polarMotion
|tdbjuliandate