使用matlab代码计算太阳高度角
使用matlab计算太阳高度角
输入的变量为
year(年份)month(月份)day(日期)time(24小时制时间)latitude(纬度)
输出
hangle(即太阳高度角)
代码
function[hangle]=sunshade1(year,month,day,time,latitude)
monthdays=[31 28 31 30 31 30 31 31 30 31 30 31];
latitude=deg2rad(latitude)
for i=1:month-1day=day+monthdays(i);%积日
end
N0=79.6764+0.2422*(year-1985)-floor((year-1985)/4);%floor是向下取整函数
sunangle=2*pi*(day-N0)/365.2422;//日角
degangle=0.3723+23.2567*sin(sunangle)+0.1149*sin(2*sunangle)-0.1712*sin(3*sunangle)-0.758*cos(sunangle)+0.3656*cos(2*sunangle)+0.0201*cos(3*sunangle)//太阳赤纬
degangle=deg2rad(degangle);//角度转弧度
timeangle=deg2rad((time-12)*15);//时角
hangle=asin(sin(latitude)*sin(degangle)+cos(latitude)*cos(degangle)*cos(timeangle));
hangle=rad2deg(hangle)
end
参考太阳高度角定义
使用matlab代码计算太阳高度角相关推荐
- 最大流matlab代码,计算最小费用最大流MATLAB源代码.doc
计算最小费用最大流MATLAB源代码,文件名为mp_mc.mfunction[Mm,mc,Mmr]=mp_mc(a,c)A=a; %各路径最大承载流量矩阵C=c; %各路径花费矩阵Mm=0; %初始可 ...
- 一个计算水文质量变化产生的负荷变形的Matlab代码
地球是一个由大气圈.水圈.岩石圈组成的复杂系统,其内部发生众多的地球动力学过程,比如地震活动.火山喷发.水文负荷等,均涉及到地球在内外力作用下的变形.地球的负荷变形有着相当长的研究历史,在1911年, ...
- matlab中服从高斯分布的矩阵_一些张量的计算步骤matlab代码
Matlab工具箱Tensor Toolbox,安装链接如下 CSDN-专业IT技术社区-登录blog.csdn.net . , 这是我们定义的张量 , matlab代码: %% %工具包:Tens ...
- ahp层次分析法matlab代码_(案例)AHP层次决策分析Matlab编码计算
"层次聚类分析Matlab编码计算" 01 - AHP层次决策分析计算函数 求判断矩阵最大特征根和归一化特征向量: function [maxEigVal,w] = maxEigV ...
- 单频信号的相位谱计算与误差修正-附Matlab代码
一.问题描述 我们在实际处理时经常遇到只有一个正弦信号的情况,其频率为 f 0 {{f}_{0}} f0,在谱分析以后,除了在频率为 f 0 {{f}_{0}} f0处有相位数值外,其他频率处都有 ...
- matlab圆周率计算,matlab代码求圆周率的简单算法
说起圆周率的算法很多人都会想起一大堆的无穷级数等各种表达式,但是这样的算法需要比较高的数学推理水品,而且对于很多初学者而言很难理解.程序员不能仅仅是机械的写程序,必须要真正的理解程序中每个代码的意义, ...
- 信号相角位移量的计算与信号位移计算-附Matlab代码
一.初始相角的位移量 在信号处理中正弦信号经常表示为 x ( n ) = A cos ( 2 π f 0 n / f s + θ ) x\left( n \right)=A\cos (2\pi { ...
- matlab 计算变异系数,变异系数法求权重matlab 代码
利用matlab编程,很容易根据变异系数法,求得多指标综合评价模型的权重. 代码如果有不懂的地方,可以联系我. 变异系数法求权重matlab 代码 clear;clc; [data1,header1] ...
- matlab 计算变异系数,变异系数法求权重matlab代码
<变异系数法求权重matlab代码>由会员分享,可在线阅读,更多相关<变异系数法求权重matlab代码(1页珍藏版)>请在读根文库上搜索. 1.变异系数 法求权重 matlab ...
- matlab三相短路电流计算程序_电力系统三相短路计算的MATLAB代码.doc
电力系统三相短路计算的MATLAB代码 电力系统三相短路计算 main.m clear tim %打开文件 [dfile,pathname]=uigetfile('*.m','Select Data ...
最新文章
- 中国移动用户能不能用WCDMA网?(世界杯与通信2)
- linux shell里面nohup的用法
- python flask 跨域问题 解决方法
- 不从事编程、学python有用吗-为什么你觉得C语言什么都不能做,学了没用?不可能的...
- 对CORS OPTIONS预检请求的一些思考
- python编码和解码_Python中的编码与解码(转)
- php中奖概率算法,可用于刮刮卡,大转盘等抽奖算法
- poj2031(prim)
- weblogic部署方式
- jq+ajax前端上传多张图片_史上最轻量的前端框架-VanillaJS
- OSG 实现跟随节点的相机(转)
- Javascript实现鼠标框选元素后拖拽被框选的元素
- HDU 5009 Paint Pearls
- Linux下如何进入Redis客户端
- layer 在标题中加点击监听事件
- Android UI切图命名规范
- (五)A股市盈率下降近三成
- 一篇对于了解我自己,挖掘我自己,从而成长的文章
- 练习聚合函数的综合练习题
- android 版本分布 最新,安卓手机系统版本分布:Android 9.0占比达34%!
热门文章
- linux小白家教学一
- 人工智能+建筑,会产生什么?
- Mybatis+Servlet+jsp
- 使用Fiddler修改百度的Logo
- 霍夫斯特拉大学计算机科学专业排名,霍夫斯特拉大学专业排名一览(含历年专业排名信息,USNEWS美国大学排名版)...
- PS CS6移动工具中图层与组设置快捷键切换(ps动作录制的妙用)
- 1_绪论 分布式机器学习
- 我想请教一下我这个错误的原因以及怎么解决,求各位大佬帮忙,感激不尽。
- 互联网经济催生了一些新职业,带来新机遇!
- [Python人工智能] 十九.Keras搭建循环神经网络分类案例及RNN原理详解