#include"math.h"

intGetDayOfYear(intYear,intMonth,intDay)
{
//功能：給定年月日得到此日在一年中的天數(shù)，1月1日為第1天，分閏年和平年
//作者：汪幫主2009-6-6wzj08@mails.jlu.edu.cn

intDayOfYear;
intFebDays;

if((Year%4==0&&Year%100!=0)||(Year%400==0))
{
FebDays=29;
}
else
{
FebDays=28;
}

DayOfYear=Day;
if(Month>1)DayOfYear=DayOfYear+31;
if(Month>2)DayOfYear=DayOfYear+FebDays;
if(Month>3)DayOfYear=DayOfYear+31;
if(Month>4)DayOfYear=DayOfYear+30;
if(Month>5)DayOfYear=DayOfYear+31;
if(Month>6)DayOfYear=DayOfYear+30;
if(Month>7)DayOfYear=DayOfYear+31;
if(Month>8)DayOfYear=DayOfYear+31;
if(Month>9)DayOfYear=DayOfYear+30;
if(Month>10)DayOfYear=DayOfYear+31;
if(Month>11)DayOfYear=DayOfYear+30;
returnDayOfYear;
}

voidGetSolarAngle(intyear,intmonth,intday,inthour,intmin,doublesec,doublelat,doublelon,intTimeZone,intN,double&AzimuthAngle,double&ZenithAngle)
{
//功能：計算太陽天頂角和太陽方位角
//作者：汪幫主2009-6-6wzj08@mails.jlu.edu.cn
doubleN0,sitar,ED,dLon,Et,Ho,viewAng;
constdoublePI=3.14159265;
doubledTimeAngle,gtdt,latitudeArc,latitude,HeightAngleArc;
doubleAzimuthAngleArc,HANoon,dTA,TimeNoon,CosAzimuthAngle,HeightAngle;

N0=79.6764+0.2422*(year-1985)-floor((year-1985)/4.0);

sitar=2*PI*(N-N0)/365.2422;
ED=0.3723+23.2567*sin(sitar)+0.1149*sin(2*sitar)-0.1712*sin(3*sitar)-0.758*cos(sitar)+0.3656*cos(2*sitar)+0.0201*cos(3*sitar);
ED=ED*PI/180;//ED本身有符號，無需判斷正負。
//ED=(m_northOrSouth==0)?ED:(-ED);

dLon=0.0;
//地球上某一點與其所在時區(qū)中心的經(jīng)度差
if(lon>=0)
{
if(TimeZone==-1)
{
dLon=lon-(floor((lon*10-75)/150)+1)*15.0;
}
else
{
dLon=lon-TimeZone*15.0;
}
}
else
{if(TimeZone==-1)
{
dLon=(floor((lon*10-75)/150)+1)*15.0-lon;
}
else
{
dLon=TimeZone*15.0-lon;
}
}

Et=0.0028-1.9857*sin(sitar)+9.9059*sin(2*sitar)-7.0924*cos(sitar)-0.6882*cos(2*sitar);//視差

gtdt=hour+min/60.0+sec/3600.0+dLon/15;//地方時
gtdt=gtdt+Et/60.0;
dTimeAngle=15.0*(gtdt-12);
dTimeAngle=dTimeAngle*PI/180;
latitudeArc=lat*PI/180;

HeightAngleArc=asin(sin(latitudeArc)*sin(ED)+cos(latitudeArc)*cos(ED)*cos(dTimeAngle));

CosAzimuthAngle=(sin(HeightAngleArc)*sin(latitudeArc)-sin(ED))/cos(HeightAngleArc)/cos(latitudeArc);
AzimuthAngleArc=acos(CosAzimuthAngle);
HeightAngle=HeightAngleArc*180/PI;
AzimuthAngle=AzimuthAngleArc*180/PI;

if(dTimeAngle<0)
{
AzimuthAngle=180-AzimuthAngle;
}
else
{
AzimuthAngle=180+AzimuthAngle;
}
ZenithAngle=90-HeightAngle;
}

總結(jié)

以上是生活随笔為你收集整理的太阳方位角的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。