利用日期、经纬度求日出日落时间 C语言程序代码

作者&投稿：藤倩（若有异议请与网页底部的电邮联系）

谁有利用经纬度计算日出日落时间的计算公式，我在~

我能给你一个粗略的计算公式
日出时间计算公式：
24*(180+时区*15-经度-acos(-tan(-23.4*cos(360*(日期序列数+9)/365))*tan(纬度))/360
我国时区为东8区，时区=8
经度、纬度采用角度制，东经、北纬为正，西经、南纬为负
日期序列数为当天在这一年中的序列，如2月11日就是42
因为计算机一般采用弧度制，上面公式可变化为：
24*(180+时区*15-经度-acos(-tan(-23.4*cos(2*π*(日期序列数+9)/365)*π/180)*tan(纬度*π/180))*180/π)/360
计算结果是一个小于24的数值，如6.69，表示6：41
日落时间计算公式：24*(1+(时区*15-经度)/180)-日出时间

我能给你一个粗略的计算公式
日出时间计算公式：
24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360
我国时区为东8区，时区=8
经度、纬度采用角度制，东经、北纬为正，西经、南纬为负
日期序列数为当天在这一年中的序列，如2月11日就是42
因为计算机一般采用弧度制，上面公式可变化为：
24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(2*π*(日期序列数+9)/365)*π/180)*TAN(纬度*π/180))*180/π)/360
计算结果是一个小于24的数值，如6.69，表示6：41

日落时间计算公式：24*(1+(时区*15-经度)/180)-日出时间

#define PI 3.1415926 

 #include<math.h> 

 #include<iostream> 

 using namespace std; 


 int days_of_month_1[]={31,28,31,30,31,30,31,31,30,31,30,31}; 

 int days_of_month_2[]={31,29,31,30,31,30,31,31,30,31,30,31}; 

 long double h=-0.833; 

 //定义全局变量 


void input_date(int c[]){ 

     int i; 

     cout<<"Enter the date (form: 2009 03 10):"<<endl; 

     for(i=0;i<3;i++){ 

         cin>>c[i]; 

     } 

 } 

 //输入日期 


void input_glat(int c[]){ 

     int i; 

     cout<<"Enter the degree of latitude(range: 0°- 60°,form: 40 40 40 (means 40°40′40″)):"<<endl; 

     for(i=0;i<3;i++){ 

         cin>>c[i]; 

     } 

 } 

 //输入纬度 


void input_glong(int c[]){ 

     int i; 

     cout<<"Enter the degree of longitude(west is negativ,form: 40 40 40 (means 40°40′40″)):"<<endl; 

     for(i=0;i<3;i++){ 

         cin>>c[i]; 

     } 

 } 

 //输入经度 


int leap_year(int year){ 

     if(((year%400==0) || (year%100!=0) && (year%4==0))) return 1; 

     else return 0; 

 } 

 //判断是否为闰年:若为闰年，返回1；若非闰年，返回0 


 int days(int year, int month, int date){ 

     int i,a=0; 

     for(i=2000;i<year;i++){ 

         if(leap_year(i)) a=a+366; 

         else a=a+365; 

     } 

     if(leap_year(year)){ 

         for(i=0;i<month-1;i++){ 

             a=a+days_of_month_2[i]; 

         } 

     } 

     else { 

         for(i=0;i<month-1;i++){ 

             a=a+days_of_month_1[i]; 

         } 

     } 

     a=a+date; 

     return a; 

 } 

 //求从格林威治时间公元2000年1月1日到计算日天数days 


 long double t_century(int days, long double UTo){ 

     return ((long double)days+UTo/360)/36525; 

 } 

 //求格林威治时间公元2000年1月1日到计算日的世纪数t 


 long double L_sun(long double t_century){ 

     return (280.460+36000.770*t_century); 

 } 

 //求太阳的平黄径 


long double G_sun(long double t_century){ 

     return (357.528+35999.050*t_century); 

 } 

 //求太阳的平近点角 


long double ecliptic_longitude(long double L_sun,long double G_sun){ 

     return (L_sun+1.915*sin(G_sun*PI/180)+0.02*sin(2*G_sun*PI/180)); 

 } 

 //求黄道经度 


long double earth_tilt(long double t_century){ 

     return (23.4393-0.0130*t_century); 

 } 

 //求地球倾角 


long double sun_deviation(long double earth_tilt, long double ecliptic_longitude){ 

     return (180/PI*asin(sin(PI/180*earth_tilt)*sin(PI/180*ecliptic_longitude))); 

 } 

 //求太阳偏差 


long double GHA(long double UTo, long double G_sun, long double ecliptic_longitude){ 

     return (UTo-180-1.915*sin(G_sun*PI/180)-0.02*sin(2*G_sun*PI/180)+2.466*sin(2*ecliptic_longitude*PI/180)-0.053*sin(4*ecliptic_longitude*PI/180)); 

 } 

 //求格林威治时间的太阳时间角GHA 


 long double e(long double h, long double glat, long double sun_deviation){ 

     return 180/PI*acos((sin(h*PI/180)-sin(glat*PI/180)*sin(sun_deviation*PI/180))/(cos(glat*PI/180)*cos(sun_deviation*PI/180))); 

 } 

 //求修正值e 


 long double UT_rise(long double UTo, long double GHA, long double glong, long double e){ 

     return (UTo-(GHA+glong+e)); 

 } 

 //求日出时间 


long double UT_set(long double UTo, long double GHA, long double glong, long double e){ 

     return (UTo-(GHA+glong-e)); 

 } 

 //求日落时间 


long double result_rise(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date){ 

     long double d; 

     if(UT>=UTo) d=UT-UTo; 

     else d=UTo-UT; 

     if(d>=0.1) { 

         UTo=UT; 

         UT=UT_rise(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))))); 

         result_rise(UT,UTo,glong,glat,year,month,date); 

     } 

     return UT; 

 } 

 //判断并返回结果（日出） 


long double result_set(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date){ 

     long double d; 

     if(UT>=UTo) d=UT-UTo; 

     else d=UTo-UT; 

     if(d>=0.1){ 

         UTo=UT; 

         UT=UT_set(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))))); 

         result_set(UT,UTo,glong,glat,year,month,date); 

     } 

     return UT; 

 } 

 //判断并返回结果（日落） 


int Zone(long double glong){ 

     if(glong>=0) return (int)((int)(glong/15.0)+1); 

     else return (int)((int)(glong/15.0)-1); 

 } 

 //求时区 


void output(long double rise, long double set, long double glong){ 

     if((int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))<10) 

         cout<<"The time at which the sun rises is "<<(int)(rise/15+Zone(glong))<<":0"<<(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))<<" .
"; 

     else cout<<"The time at which the sun rises is "<<(int)(rise/15+Zone(glong))<<":"<<(int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))<<" .
"; 

     if((int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))<10) 

         cout<<"The time at which the sun sets is "<<(int)(set/15+Zone(glong))<<": "<<(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))<<" .
"; 

     else cout<<"The time at which the sun sets is "<<(int)(set/15+Zone(glong))<<":"<<(int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))<<" .
"; 

 } 

 //打印结果 

         

int main(){ 

     long double UTo=180.0; 

     int year,month,date; 

     long double glat,glong; 

     int c[3]; 

     input_date(c); 

     year=c[0]; 

     month=c[1]; 

     date=c[2]; 

     input_glat(c); 

     glat=c[0]+c[1]/60+c[2]/3600; 

     input_glong(c); 

     glong=c[0]+c[1]/60+c[2]/3600; 

     long double rise,set; 

     rise=result_rise(UT_rise(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))))),UTo,glong,glat,year,month,date); 

     set=result_set(UT_set(UTo,GHA(UTo,G_sun(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))),glong,e(h,glat,sun_deviation(earth_tilt(t_century(days(year,month,date),UTo)),ecliptic_longitude(L_sun(t_century(days(year,month,date),UTo)),G_sun(t_century(days(year,month,date),UTo)))))),UTo,glong,glat,year,month,date); 

     output(rise,set,glong); 

     system("pause"); 

     return 0; 

 }

求太阳直射点的经纬度
太阳直射点的经纬度可以根据特殊日期来确定。1. 根据特殊日期（二分二至日）计算。当北半球夏至日（6月22日前后）时，太阳直射点的纬度为23°26”N；当北半球冬至日（12月22日前后）时，太阳直射点的纬度为23°26”S；当北半球春分日（3月21日前后）和秋分日（9月23日前后）时，太阳直射点的...

java根据经纬度算出日出日落时间
那么我们回家拿个灯照一个球形体，并模拟球形赤道线与太阳的夹角在变化（季节变化）就能知道怎么计算了，当然还需要把这个计算的逻辑转换成为数学方法；这里面我们需要计算的就是(对于某个季度，比如假设太阳照在北纬23.5 度时）对于每一个纬度圈，太阳能照到的弧长是多少，照到和照不到的也就是对应着...

求根据经纬度和日期时间计算太阳高度及太阳东西方向公式。
太阳赤纬（与太阳直射点纬度相等）以δ表示，观测地地理纬度用φ表示（太阳赤纬与地理纬度都是北纬为正，南纬为负）正午太阳高度角h=90°-|φ-δ|

太阳直射南北回归线的时间是什么时候?
回归线是太阳每年在地球上直射来回移动的分界线，太阳直射北回归线的这一天是北半球的夏至日，同时也是南半球的冬至日。南北回归线的经纬度。南北回归线就是南纬和北纬纬度为23°26′的那条纬线（北回归线：23°26′N，南回归线23°26′S）。南回归线是太阳在南半球能够直射到的最远位置，大约在南纬...

不同经度纬度地区位于同一晨线上判断日期以及原因
在古代航海中关于经纬度的判断一直是个很重要的问题,我给你解释下几种方法:1:关于纬度:纬度测量一直是相对来说比较容易的,主要依靠北极星的高度来测定纬度,在北半球非常适用,北极星高度算高度作为古代航海的纬度测量方法中最精确最方便的。还有就是依靠太阳中天高度测定当地纬度法也是很容易的。2:关于经度:经度的测量...

北京时间2015年10月30日15时求此时华盛顿北纬77度的日期和区时
华盛顿在北纬77°？？？怎么我了解到的是北纬39°呢？？再说了计算时间，要个纬度来干吗？华盛顿经纬度：38°53'N，77°02'W；时间计算只需要经度，这里取 77°W；77°W所在时区：77°W÷15°\/区=西5区……2°W，余数（2°W）小于半时区跨度（7.5°\/半区），舍去，即华盛顿属于西5区。北...

举个列子,太阳升起60经纬度,太阳落山时还是不是60?
地球经度是记录和判断不同地区时间早晚的依据。地区不同，经纬度就不同，经度不同，日出时刻就不同，纬度不同，一般昼夜长短就不同。日出时间计算公式：24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)\/365))*TAN(纬度))\/360 我国时区为东8区,时区=8。

地理题,关于测量当地经纬度的题目
根据图，影子最短时为2m，和竿长相等，太阳高度角为45°，当天为冬至日，太阳直射23°26'S，当地纬度是90°-45°-23°26'=21°34'N。影子最短即正午时间为11:56，在北京时间12:00前4分钟，当地经度就是在120°E以东1°（经度每差1°时间差4分钟），即121°E。立杆无影即太阳直射，当地在...

时间为12 月22 日前后,此时,太阳直射点的经纬度是45°E、23.5°S。北...
从图上可以看出是南半球（逆时针自转）。南极圈内为极昼，所以太阳直射南回归线，时间为12月22日，北半球的节气是冬至，太阳直射点的经纬度是45°E、23.5°S。即45°E时间为12点，北京时间（120°E）是17：00(东经45度属于东三区，所以东三区为十二时。而北京在东八区，所以加上五个小时)...

日分界线的经纬度是多少
是东、西经180度，从南极到北极贯穿所有纬度线。但实际上并不是南北一条直线，原因是为了照顾处在这条经度线上的国家，避免产生一个国家同时存在两个日期的现象。位于这条线东侧的时间和西侧的时间相同，日期相差一天。如果东侧是六月三十日零点，西侧就是七月一日零点。具体到每一个国家，其日期变更...

昌乐县19176126562： 高中地理... 知道经纬度如何求日出日落时间呢 - ？
宁琴云南： 日出时间计算公式: 24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360 我国时区为东8区,时区=8 经度、纬度采用角度制,东经、北纬为正,西经、南纬为负日期序列数为当天在这一年中的序列,...

昌乐县19176126562： ..知道经纬度如何求日出日落时间呢 - ？
宁琴云南：[答案] 日出时间计算公式:24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360我国时区为东8区,时区=8经度、纬度采用角度制,东经、北纬为正,西经、南纬为负日期序列数为当天在这一年中的序...

昌乐县19176126562： 请问怎样根据经纬度计算日出日日落时间 - ？
宁琴云南：[答案] 日出时间计算公式: 24*(180+时区*15-经度-ACOS(-TAN(-23.4*COS(360*(日期序列数+9)/365))*TAN(纬度))/360 我国时区为东8区,时区=8 经度、纬度采用角度制,东经、北纬为正,西经、南纬为负日期序列数为当天在这一年中的序列,如2月...

昌乐县19176126562： 利用日期、经纬度求日出日落时间 C语言程序代码 - ？
宁琴云南：#define PI 3.1415926 #include #include using namespace std;int days_of_month_1[]={31,28,31,30,31,30,31,31,30,31,30,31}; int days_of_month_2[]={31,29,31,30,31,30,31,31,30,31,30,31}; long double h=-0.833; //定义全局变...

昌乐县19176126562： 请问,如何通过经纬度算出当地的日出日落时间 - ？
宁琴云南： 很复杂的公式很复杂的公式,我以前试图用单片机完成,没有成功,可以用近似的方法倒是可以,但是有点误差. 纬度算出日出日落时间的算法?应该还有日期吧!要不然日出日落可不一样!是不是要做太阳能的控制装置?呵呵你应该到一些天文网站找相关的资料还有说比如知道北京的日出日落, 拿这个已知的跟北京的经度纬度做比较看

昌乐县19176126562： 求根据经纬度和日期时间计算太阳高度及太阳东西方向公式.还有就是能不能计算出当天日出日落时间?公式是什么? - ？
宁琴云南：[答案] 太阳赤纬(与太阳直射点纬度相等)以δ表示,观测地地理纬度用φ表示(太阳赤纬与地理纬度都是北纬为正,南纬为负) 正午太阳高度角h=90°-|φ-δ|

昌乐县19176126562： 已知一个点的经纬求日出日落时间？
宁琴云南：首先,要引入晨昏线的概念(关于晨昏线详情参阅资料 http://baike.baidu.com/view/218583.htm) 因为赤道上是终年昼夜平分的,也就是早上6点太阳升起,晚上6点太阳落下.晨昏线只有在春分秋分的时候是与经线相重合的,其余的时间都是与经线相切的其次,如果你知道一点的位置,看该点的纬线与晨线(或昏线这里用晨线举例)相交的那点的经度与此时赤道上为6时的那根经线之间的经度数之差,再算出时差,就可以推出此时该点所在的纬线的日出时间,日落时间则是12-日出时间=日落时间了如果有图,解释起来会简单些,可惜好像粘不上图~

昌乐县19176126562： 求日出、日落时间计算公式已知经纬度(单位:角度制的度),日期(1月1日记为1,没过一天加一,不考虑闰年)假设地球为正球体,不考虑大气折射,近... - ？
宁琴云南：[答案] 春秋分为全球昼夜等分3.21 日出时间为:6时日落时间为18时即昼长=日落时间-日出时间 18-6=12个小时日出时间=12-12/2=6日落时间=12+12/2=18函数如下:日出时间:12-昼长/2日落时间:12+昼长/2 其实记住一点就ok...

昌乐县19176126562： 知道日期和纬度怎样求当天日出时间? - ？
宁琴云南： 你好,这个条件不足,不能求出准确时间.但是可以肯定是楼上的说法不正确.应该是在六点之后.因为12月17日,已经过了秋分日(9月22日),也就是说太阳直射点已经到了南纬21.5度附近,北纬24度当天最大太阳高度角也只能是44.5度左右,昼短夜长,所以日出应是在六点之后,日落在18点之前.希望对你有所帮助.

昌乐县19176126562： 求根据日期、经纬度计算太阳起落时间的EXCEL表格 - ？
宁琴云南： 众所周知的:当太阳升起的太阳的位置,并设置H = -0.833°,你要计算的地理位置,经度纬度G1at龙,时区,区反对派最后的日出和日落时间计算,首先计算反对派= 180°. (1)第一种计算法GMT 2000年1月1日AD的数量来计算的天数天; (2...

你可能想看的相关专题

星空见康网

利用日期、经纬度求日出日落时间 C语言程序代码

你可能想看的相关专题