||
参考链接:http://bbs.06climate.com/forum.php?mod=viewthread&tid=36366
代码部分:
function HSI=calHSI(year,month,day,hour,min,lat,lon)
%HSI太阳高度角(单位:°)
%year年
%month月
%day日
%hour小时
%min分钟
%lat纬度(单位:°)
%lon经度(单位:°)
d=datenum(year,month,day)-datenum(year,1,1)+1;
T=2*pi*(d-1)/365; %单位:弧度
SD=(0.006918-0.399912*cos(T)+0.070257*sin(T)-0.006758*cos(2*T)+0.000907*sin(2*T)-0.002697*cos(3*T)+0.00148*sin(3*T));
ST=hour+min/60+(lon-120)/15;
w=15*(ST-12);
HSI=asind(sind(lat)*sin(SD)+cosd(lat)*cos(SD)*cosd(w));
end
另外,在excel中对2011/2/25 14:00:00格式的数据提取年、月、日、小时、分钟的方法:
=YEAR(A1)
=MONTH(A1)
=DAY(A1)
=HOUR(A1)
=MINUTE(A1)
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 06:36
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社