太阳角度计算c语言程序

1、精选优质文档-倾情为你奉上*/#include "main.h"#include "cracker.h"#include <math.h>/*-*Function name: leap* Descriptions: 判断当年是否为闰年*input parameters: pTime -日期时间*output parameters: leap（闰年返回1，平年返回0）*-*/ int leap(PCF8563_Time *pTime) int leap; leap=(pTime->Year%4=0&&pTime->Y

2、ear%100!=0)|(pTime->Year%400=0); return leap;/*-*Function name: ：DayAddCalculate* Descriptions: 计算当天为一年中的第几天*input parameters: pTime -日期时间*output parameters: DayAdd-积日*-*/INT32U DayAddCalculate(PCF8563_Time *pTime ) INT32U DayAdd; switch(pTime->Month) case 1: DayAdd=(pTime->Day);break; case

3、 2: DayAdd=31+pTime->Day;break; case 3: DayAdd=31+28+pTime->Day;break; case 4: DayAdd=31+28+31+pTime->Day;break; case 5: DayAdd=31+28+31+30+pTime->Day;break; case 6: DayAdd=31+28+31+30+31+pTime->Day;break; case 7: DayAdd=31+28+31+30+31+30+pTime->Day;break; case 8: DayAdd=31+28+31+3

4、0+31+30+31+pTime->Day;break; case 9: DayAdd=31+28+31+30+31+30+31+31+pTime->Day;break; case 10: DayAdd=31+28+31+30+31+30+31+31+30+pTime->Day;break; case 11: DayAdd=31+28+31+30+31+30+31+31+30+31+pTime->Day;break; case 12: DayAdd=31+28+31+30+31+30+31+31+30+31+30+pTime->Day;break; default

5、: break; /*cout<<"月份输入错误！*/ if(leap(pTime)&&(pTime->Month>=3) DayAdd=DayAdd+1; return DayAdd;/*-*Function name: ：DecAngleCalculate* Descriptions: 计算太阳的赤纬角*input parameters: pTime -日期时间*output parameters:赤纬角DecAngle*-*/FP64 DecAngleCalculate(PCF8563_Time *pTime ) INT32U DayAd

6、d; DayAdd=DayAddCalculate(pTime); return (23.45*sin(2*PI*(284+DayAdd)/365)*PI/180;/*-*Function name: ：HourAngleCalculate* Descriptions: 计算当前时角*input parameters: pTime -日期时间*output parameters:HourAngle*/FP64 HourAngleCalculate(PCF8563_Time *pTime ) return (12-(pTime->Hour)-(pTime->Minute)/60.0-

7、(pTime->Second)/3600.0)*PI/12;/*-*Function name: ：EleAngleCalculate* Descriptions: 计算当前太阳的高度角*input parameters: pTime -日期时间*output parameters:EleAngle*-*/FP64 EleAngleCalculate(PCF8563_Time *pTime ) FP64 DecAngle; FP64 HourAngle; DecAngle=DecAngleCalculate(pTime); HourAngle=HourAngleCalculate(pTi

8、me); return asin(sin(LAT)*sin(DecAngle)+cos(LAT )*cos(DecAngle)*cos(HourAngle);/*-*Function name: ：AziAngleCalculate* Descriptions: 计算当前太阳的方位角*input parameters: pTime -日期时间*output parameters:AziAngle*-FP64 AziAngleCalculate(PCF8563_Time *pTime ) FP64 DecAngle; FP64 EleAngle; DecAngle=DecAngleCalculate(pTime); EleAngle=EleAngleCalculate(pTime ); if(pTime->Hour>12) return acos(sin(EleAngle)*sin(LAT )-sin(DecAngle)/(cos(EleAngle)*cos(LAT ); else re