sun_angle文档
sun_angle给出了太阳方位角和仰角地球上任何时间在任何地点。这个函数是改编自SolarAzEl达林c . Koblick函数。
内容
语法
[az, el] = sun_angle (t,纬度,经度)[az, el] = sun_angle (t,纬度,经度,h)
描述
(阿兹,el) = sun_angle (t,纬度,经度)让太阳方位角和仰角度在指定的地理位置和时间t在UTC。输入t可以datenum、datestr或datetime格式。
(阿兹,el) = sun_angle (t,纬度,经度,h)指定在米的海拔高度。如果没有指定高度,高度会自动通过CDT功能决定的topo_interp。
示例1:时间序列在一个位置
在我写这篇文章时,我坐在咖啡店帕萨迪纳市加州(34.16 n, 118.13 w)。现在确定太阳的方位角和仰角,在接下来的10天。
t = linspace(现在,现在+ 10,10000);%一系列时间10天t = t + 8/24;%将太平洋标准时间转换为UTC(阿兹,el) = sun_angle (t、34.16、-118.13);
使用subsubplot绘制太阳的角度:
subsubplot(2, 1, 1)情节(t, az)轴紧ylabel的方位(度)datetick (“x”,“keeplimits”)subsubplot(2, 1, 2)情节(t, el)轴紧ylabel的海拔(度)datetick (“x”,“keeplimits”)
示例2:全球网格
使用cdtgrid创建一个quarter-degree全球网格,得到全球各地太阳角度UTC罢工的午夜5月27日,2019:
(纬度、经度)= cdtgrid (1/4);(阿兹,el) = sun_angle (“1984年5月27日,就是”纬度,经度);
使用imagescn绘制太阳方位并设置colormapcmocean阶段。为上下文使用添加政治界限边界。
图imagescn(经度、纬度、az) cb = colorbar;ylabel (cb、“太阳方位角(度)”360年)caxis ([0]) cmocean阶段边界(“国家”,“颜色”rgb (“灰色”))
覆盖太阳高度角的轮廓。当然,太阳仰角只是物理海拔0以上太阳在哪里,所以只有包括积极的轮廓:
持有在轮廓(经度、纬度、el [0 0),“k”)%固体零线轮廓(经度、纬度、el、0:10:90凯西:”)文本(0,0,{“年代夜间;“到处都是南”;“固体的轮廓”},…“水平的”,“中心”)
示例3:微小的海拔高度的影响
太阳的角是受海拔高度的影响。默认情况下,sun_angle使用topo_interp确定地表海拔高程如果没有指定仰角。但是我们考虑一个假想的塔,高10000米在伦敦的中心(51.51 n, 0.123 w)。什么太阳方位角和仰角沿着这塔7月14日下午7点45 2019 ?
t =2019年7月14日晚上7点的;lat = 51.51;朗= -0.123;z = 0:10000;%沿着塔高度(米)(阿兹,el) = sun_angle (t,纬度,经度,z);图次要情节(1、2、1)情节(az-az (1), z)%情节相对于基地轴紧包含‘\δ方位(度)ylabel的高度(米)次要情节(1、2、2)情节(el-el (1), z)%情节相对于基地轴紧包含‘\三角高程(度)
注意x轴上的科学记数法。这告诉我们高程的影响是微乎其微的几千米的规模。
作者信息
这个函数是由达林c Koblick使用这里描述的公式:http://stjarnhimlen.se/comp/tutorial.html # 5。这是改编的气候数据为Matlab工具箱2019年由乍得格林。