matlab求cos角,科学网—MATLAB求太阳高度角的小程序

参考链接：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)

转载本文请联系原作者获取授权，同时请注明本文来自张乐乐科学网博客。

链接地址：http://blog.sciencenet.cn/blog-1103122-935807.html

上一篇：MODIS积雪产品代码含义

下一篇：Noah-MP模式中的一个小错误

matlab求cos角,科学网—MATLAB求太阳高度角的小程序 - 张乐乐的博文相关推荐

matlab画波传播,科学网—MATLAB绘制波包曲线 - 李金磊的博文
根据 Born 和 Wolf 的 Principles of Optics 一书中1.3节的描述,以及 Fig1.5 ,可绘制出波包曲线如下: 图1 波包曲线图具体程序为: t=-5:0.01:5 ...
matlab做误差棒图,科学网—Matlab误差棒、图例、颜色配比，很有用的代码 - 周锋的博文...
话不多说,直接贴代码和效果图,会用之人自然会觉得有用. 效果图: 代码: clear; clc; load('ALGO_GPS.neu'); % 加载数据 data = ALGO_GPS; dyr ...
matlab箭头大小固定,科学网—Matlab 粗箭头绘制 - 肖鑫的博文
matlab官方工具箱介绍是一个很神奇的地方,原贡献者和评论中的贡献者非常多能帮助解决很多疑难问题,例如上一篇博文便是来自于官方解答最初得到arrow.m是来自叶师兄的工具包,但是这个版本非常老了 ...
matlab 画qq图,科学网—[转载]R语言绘制QQ图 - 刘朋的博文
R语言绘制QQ图实例1: #############加载数据 data R R=apply(R,2,as.numeric) #R语言将字符串矩阵转化为数值型矩阵,apply()函数里面的第2个值,如 ...
matlab stract结构_科学网—[MATLAB]方便快捷读取结构体里数个结构体内的数据 - 胡振东的博文...
clc;clear;close all cd F:\01_DATA\Data_process\Glide load('Glide0.mat') % Glide0.mat里的Glide结构体有G1,G2 ...
matlab对数收益直方图,科学网—MATLAB中绘制数据直方图的新函数histogram2 - 王福昌的博文...
MATLAB中有命令hist3() 可以绘制直方图,竖坐标是频数,这与一些教科书中用纵轴表示频率的做法不一致,有些时候不便于使用.当然,使用者可以自己编写定制能够在纵轴绘出频率的直方图.在MATLAB ...
matlab使用load指令,科学网—matlab常用方法 - 陈超的博文
之前将数组或者矩阵保存为一个mat格式的文件,在进行load命令读取时: s1=load('qiyipuzong.mat'): 得到的s1是struct类型的数据,而我想要的是一个矩阵或者数组. 经过 ...
matlab画波传播,科学网-MATLAB绘制波包曲线-李金磊的博文
根据 Born 和 Wolf 的 Principles of Optics 一书中1.3节的描述,以及 Fig1.5 ,可绘制出波包曲线如下: 图1 波包曲线图具体程序为: t=-5:0.01:5 ...
matlab nan 无色_科学网—Matlab 关于NaN值的填充 - 肖鑫的博文
以前认为很麻烦的事往往一行命令就能搞定,代码如下 figure('position',[100 100 600 500],'PaperPositionMode', 'auto') m_proj('Sa ...

matlab求cos角,科学网—MATLAB求太阳高度角的小程序 - 张乐乐的博文

matlab求cos角,科学网—MATLAB求太阳高度角的小程序 - 张乐乐的博文相关推荐

最新文章

热门文章