这个例子显示了如何估算太阳的安息。分析是表示相对于地球表面上的特定地理位置从天体球体上的平均位置的角度偏移的变化的曲线。在这个例子中,Analemma相对于英国格林威治的皇家天文台估计。估计后,示例绘制了analemma。
此示例使用您可以使用的数据下载Aerodatapackage.命令。
指定计算Analemma的日期。在这个例子中,这些日期为2014年1月1日至12月31日在中午UTC的年度。
DV = DateTime(2014,1,1:365,12,0,0);dvutc = [dv.year'dv.month'dv.day'dv.hour'dh.minute'dv.second'];
使用Planetephemeris函数来计算Sun的位置。在这个例子中:
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的地面时间来估算动态重心时间的Julian日期:
jdtdb = tdbjuliandate([dvtt.year'dvtt.month'dvtt.day'dvttt.hour'dvtt.minute'dvttt.second']);
确定太阳的位置:
Hostun = Planetephemeris(JDTDB,'地球'那'太阳')* 1000;
计算UTC和UT1,Deltaut1之间的差异,使用UTC的修改后的朱利安日期。
mjdutc = mjuliandate(dvutc);dut1 = deltaut1(mjdutc);
计算CIP的极性运动和位移使用修改的UTC的朱利安日期。
PM = Polarmotion(MJDUTC);dcip = deltacip(mjdutc);
指定用于估计安疗的位置的地理位置位置。在这个例子中,这个位置是格林威治皇家天文台的纬度,经度和高度(北北,0.0015度西,0米高度)。
llagreenwich = [51.48,-0.0015,0];AER = ECI2AER(POSSUN,DVUTC,REPMAT(LLAGRENWICH,长度(JDTDB),1),......'deltaat',dat * on(长度(jdtdb),1),'deltaut1',dut1,......'polarmotion',下午,'dcip',dcip);
在analemma上,您可以在分析年内绘制兴趣日。此示例绘图:
2014年每个月的第一天。
夏季和冬季索尔斯特。
春天和瀑布。
获得2014年每个月的第一天:
AERFIRSTMONTH = AER(DVUTC(:,3)== 1,:);
获得索尔斯特和昼夜计划(2014年为3/20,6 / 21,9 / 22,122,121):
solsticeequinox = [AER(DVUTC(:,2)== 3,DVUTC(:3)== 20),1)AER(和(DVUTC(:,2)== 3,DVUTC(:,3)== 20),2);......AER(和(DVUTC(:,2)== 6,DVUTC(:,3)== 21),1)AER(和(DVUTC(:,2)== 6,DVUTC(:,3)== 21),2);......AER(和(DVUTC(:,2)== 9,DVUTC(:,3)== 22),1)AER(和(DVUTC(:,2)== 9,DVUTC(:,3)== 22),2);......AER(和(DVUTC(:,2)== 12,DVUTC(:3)== 21),1)AER(和(DVUTC(:,2)== 12,DVUTC(:,3)== 21),2)];
绘制analemma。沿着安排,整年的情节点,本月的第一天,昼夜平分和索斯。
为了II = 12:-1:1小星期{II} = [NUM2STR(II)'/'num2str(1)];结尾f =数字;绘图(AER(:,1),AER(:,2),'。'那......Solsticeequinox(:,1),溶精(:,2),'ks'那......Aerfirstmonth(:,1),Aerfirstmonth(:,2),'ko'那......'Markersize',8,'markerfacecolor'那'K');标题(在Greenwich观测所观察'Analemma');Xlabel('azimuth [deg]');ylabel('海拔[deg]');轴([176,185,10,70])文本(Aerfirstmonth(:,1)+。1,Aerfirstmonth(:,2)+1.2,Firstdays,'颜色'那'K'那'水平对齐'那'剩下');文字(溶精(1,1)+。2,Solsticeequinox(1,2)-1.5,'春分'那'颜色'那'K'那'水平对齐'那'剩下');文字(溶精(2,1),溶酶体喹氏素(2,2)+2.5,'夏天'那'颜色'那'K'那'水平对齐'那'剩下');文字(溶精(3,1)+。1,Solsticeequinox(3,2)+1.2,'跌落股份'那'颜色'那'K'那'水平对齐'那'剩下');文字(溶精(4,1)+。1,溶酶体喹氏蛋杆(4,2)-2.5,'冬至'那'颜色'那'K'那'水平对齐'那'剩下');