先贴在这了，后面应该用得着

    1. #define PI 3.1415926 
    1.  
    1.  #include<math.h> 
    1.  
    1.  #include<iostream> 
    1.  
    1.  using namespace std; 
    1.  
    1.  
    1.  int days_of_month_1[]={31,28,31,30,31,30,31,31,30,31,30,31}; 
    1.  
    1.  int days_of_month_2[]={31,29,31,30,31,30,31,31,30,31,30,31}; 
    1.  
    1.  long double h=-0.833; 
    1.  
    1.  //定义全局变量 
    1.  
    1.  
    1. void input_date(int c[]){ 
    1.  
    1.      int i; 
    1.  
    1.      cout<<"Enter the date (form: 2009 03 10):"<<endl; 
    1.  
    1.      for(i=0;i<3;i++){ 
    1.  
    1.          cin>>c[i]; 
    1.  
    1.      } 
    1.  
    1.  } 
    1.  
    1.  //输入日期 
    1.  
    1.  
    1. void input_glat(int c[]){ 
    1.  
    1.      int i; 
    1.  
    1.      cout<<"Enter the degree of latitude(range: 0°- 60°,form: 40 40 40 (means 40°40′40″)):"<<endl; 
    1.  
    1.      for(i=0;i<3;i++){ 
    1.  
    1.          cin>>c[i]; 
    1.  
    1.      } 
    1.  
    1.  } 
    1.  
    1.  //输入纬度 
    1.  
    1.  
    1. void input_glong(int c[]){ 
    1.  
    1.      int i; 
    1.  
    1.      cout<<"Enter the degree of longitude(west is negativ,form: 40 40 40 (means 40°40′40″)):"<<endl; 
    1.  
    1.      for(i=0;i<3;i++){ 
    1.  
    1.          cin>>c[i]; 
    1.  
    1.      } 
    1.  
    1.  } 
    1.  
    1.  //输入经度 
    1.  
    1.  
    1. int leap_year(int year){ 
    1.  
    1.      if(((year%400==0) || (year%100!=0) && (year%4==0))) return 1; 
    1.  
    1.      else return 0; 
    1.  
    1.  } 
    1.  
    1.  //判断是否为闰年:若为闰年，返回1；若非闰年，返回0 
    1.  
    1.  
    1.  int days(int year, int month, int date){ 
    1.  
    1.      int i,a=0; 
    1.  
    1.      for(i=2000;i<year;i++){ 
    1.  
    1.          if(leap_year(i)) a=a+366; 
    1.  
    1.          else a=a+365; 
    1.  
    1.      } 
    1.  
    1.      if(leap_year(year)){ 
    1.  
    1.          for(i=0;i<month-1;i++){ 
    1.  
    1.              a=a+days_of_month_2[i]; 
    1.  
    1.          } 
    1.  
    1.      } 
    1.  
    1.      else { 
    1.  
    1.          for(i=0;i<month-1;i++){ 
    1.  
    1.              a=a+days_of_month_1[i]; 
    1.  
    1.          } 
    1.  
    1.      } 
    1.  
    1.      a=a+date; 
    1.  
    1.      return a; 
    1.  
    1.  } 
    1.  
    1.  //求从格林威治时间公元2000年1月1日到计算日天数days 
    1.  
    1.  
    1.  long double t_century(int days, long double UTo){ 
    1.  
    1.      return ((long double)days+UTo/360)/36525; 
    1.  
    1.  } 
    1.  
    1.  //求格林威治时间公元2000年1月1日到计算日的世纪数t 
    1.  
    1.  
    1.  long double L_sun(long double t_century){ 
    1.  
    1.      return (280.460+36000.770*t_century); 
    1.  
    1.  } 
    1.  
    1.  //求太阳的平黄径 
    1.  
    1.  
    1. long double G_sun(long double t_century){ 
    1.  
    1.      return (357.528+35999.050*t_century); 
    1.  
    1.  } 
    1.  
    1.  //求太阳的平近点角 
    1.  
    1.  
    1. long double ecliptic_longitude(long double L_sun,long double G_sun){ 
    1.  
    1.      return (L_sun+1.915*sin(G_sun*PI/180)+0.02*sin(2*G_sun*PI/180)); 
    1.  
    1.  } 
    1.  
    1.  //求黄道经度 
    1.  
    1.  
    1. long double earth_tilt(long double t_century){ 
    1.  
    1.      return (23.4393-0.0130*t_century); 
    1.  
    1.  } 
    1.  
    1.  //求地球倾角 
    1.  
    1.  
    1. long double sun_deviation(long double earth_tilt, long double ecliptic_longitude){ 
    1.  
    1.      return (180/PI*asin(sin(PI/180*earth_tilt)*sin(PI/180*ecliptic_longitude))); 
    1.  
    1.  } 
    1.  
    1.  //求太阳偏差 
    1.  
    1.  
    1. long double GHA(long double UTo, long double G_sun, long double ecliptic_longitude){ 
    1.  
    1.      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)); 
    1.  
    1.  } 
    1.  
    1.  //求格林威治时间的太阳时间角GHA 
    1.  
    1.  
    1.  long double e(long double h, long double glat, long double sun_deviation){ 
    1.  
    1.      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))); 
    1.  
    1.  } 
    1.  
    1.  //求修正值e 
    1.  
    1.  
    1.  long double UT_rise(long double UTo, long double GHA, long double glong, long double e){ 
    1.  
    1.      return (UTo-(GHA+glong+e)); 
    1.  
    1.  } 
    1.  
    1.  //求日出时间 
    1.  
    1.  
    1. long double UT_set(long double UTo, long double GHA, long double glong, long double e){ 
    1.  
    1.      return (UTo-(GHA+glong-e)); 
    1.  
    1.  } 
    1.  
    1.  //求日落时间 
    1.  
    1.  
    1. long double result_rise(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date){ 
    1.  
    1.      long double d; 
    1.  
    1.      if(UT>=UTo) d=UT-UTo; 
    1.  
    1.      else d=UTo-UT; 
    1.  
    1.      if(d>=0.1) { 
    1.  
    1.          UTo=UT; 
    1.  
    1.          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)))))); 
    1.  
    1.          result_rise(UT,UTo,glong,glat,year,month,date); 
    1.  
    1.      } 
    1.  
    1.      return UT; 
    1.  
    1.  } 
    1.  
    1.  //判断并返回结果（日出） 
    1.  
    1.  
    1. long double result_set(long double UT, long double UTo, long double glong, long double glat, int year, int month, int date){ 
    1.  
    1.      long double d; 
    1.  
    1.      if(UT>=UTo) d=UT-UTo; 
    1.  
    1.      else d=UTo-UT; 
    1.  
    1.      if(d>=0.1){ 
    1.  
    1.          UTo=UT; 
    1.  
    1.          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)))))); 
    1.  
    1.          result_set(UT,UTo,glong,glat,year,month,date); 
    1.  
    1.      } 
    1.  
    1.      return UT; 
    1.  
    1.  } 
    1.  
    1.  //判断并返回结果（日落） 
    1.  
    1.  
    1. int Zone(long double glong){ 
    1.  
    1.      if(glong>=0) return (int)((int)(glong/15.0)+1); 
    1.  
    1.      else return (int)((int)(glong/15.0)-1); 
    1.  
    1.  } 
    1.  
    1.  //求时区 
    1.  
    1.  
    1. void output(long double rise, long double set, long double glong){ 
    1.  
    1.      if((int)(60*(rise/15+Zone(glong)-(int)(rise/15+Zone(glong))))<10) 
    1.  
    1.          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))))<<" .\n"; 
    1.  
    1.      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))))<<" .\n"; 
    1.  
    1.      if((int)(60*(set/15+Zone(glong)-(int)(set/15+Zone(glong))))<10) 
    1.  
    1.          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))))<<" .\n"; 
    1.  
    1.      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))))<<" .\n"; 
    1.  
    1.  } 
    1.  
    1.  //打印结果 
    1.  
    1.           
    1.  
    1. int main(){ 
    1.  
    1.      long double UTo=180.0; 
    1.  
    1.      int year,month,date; 
    1.  
    1.      long double glat,glong; 
    1.  
    1.      int c[3]; 
    1.  
    1.      input_date(c); 
    1.  
    1.      year=c[0]; 
    1.  
    1.      month=c[1]; 
    1.  
    1.      date=c[2]; 
    1.  
    1.      input_glat(c); 
    1.  
    1.      glat=c[0]+c[1]/60+c[2]/3600; 
    1.  
    1.      input_glong(c); 
    1.  
    1.      glong=c[0]+c[1]/60+c[2]/3600; 
    1.  
    1.      long double rise,set; 
    1.  
    1.      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); 
    1.  
    1.      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); 
    1.  
    1.      output(rise,set,glong); 
    1.  
    1.      system("pause"); 
    1.  
    1.      return 0; 
    1.  
    1.  }

利用日期、经纬度求日出日落时间 C语言程序代码（zz）的更多相关文章

  1. 利用OD破解一个简单的C语言程序

    最近在学习汇编(看的是王爽老师的<汇编语言(第三版)>),然后想尝试使用OD(Ollydbg)软件破解一个简单的C语言程序练练手. 环境: C语言编译环境:VC++6.0 系统:在Wind ...

  2. 编写一个求圆面积的C语言程序

    #include<stdio.h>   //文件包含//#define PI 3.14     //宏定义//void main()         { float r,s; scanf( ...

  3. C++/C语言程序代码

    //-----------------------------------1 #include <stdio.h> #include<stdlib.h> void main() ...

  4. C#可用的日出日落时间类

    一个现成代码的公共类库,复制下来作为一个类文件就可以调用了.一般不需要了解实现过程,各种数学公式太麻烦. 调用方法: SunTimeResult result = SunTimes.GetSunTim ...

  5. java获取日出日落时间

    import java.math.BigDecimal; import java.text.ParseException; import java.text.SimpleDateFormat; imp ...

  6. 指针直接赋值为整型AND利用宏定义求结构体成员偏移量

    首先我们要更正一个很熟悉的概念,那就是指针不仅仅是“地址”,指针还有一个很重要的特性,那就是“类型”. 指针初始化时,“=”的右操作数; 除外,该语句表示指针为空): 所以 ; 这样的代码是不允许的. ...

  7. Kotlin入门(18)利用单例对象获取时间

    前面介绍了,使用扩展函数可以很方便地扩充数组Array的处理功能,例如交换两个数组元素.求数组的最大元素等等.那么除了数组之外,日期和时间的相关操作,也是很常见的,比如获取当前日期,获取当前时间.获取 ...

  8. SQLServer 日期函数大全 SQLServer 时间函数大全

    原文地址:https://www.cnblogs.com/zhangpengnike/p/6122588.html 一.统计语句 1.--统计当前[>当天00点以后的数据] SELECT * F ...

  9. Js 日期转换函数(UTC时间转换及日期想加减)

    IOS上Js日期转换中new Date("yyyy-mm-dd")不能正常工作,必须使用new Date("yyyy/MM/dd"); 日期相加减: Date. ...

随机推荐

  1. 【模块化编程】理解requireJS-实现一个简单的模块加载器

    在前文中我们不止一次强调过模块化编程的重要性,以及其可以解决的问题: ① 解决单文件变量命名冲突问题 ② 解决前端多人协作问题 ③ 解决文件依赖问题 ④ 按需加载(这个说法其实很假了) ⑤ ..... ...

  2. [deviceone开发]-do_SegmentView和do_SlideView联动的示例

    一.简介 示例展示do_SegmentView和do_SlideView联动的使用,这二个组件很常用,而且这个组合也非常常用,类似网易新闻的效果,上面滑动带动下面的slideview滑动,反过来也是. ...

  3. [JS]笔记13之Date对象

    -->获取与设置时间的方法-->使用Date对象制作相应的效果 1.设置时间创建一个时间对象 new Date(time); 设置时间 time 从1970年1月1日至几种格式:new D ...

  4. SharePoint 2013 母版页取消和HTML页关联

    前言:在新版本的SharePoint 2013上,有新的功能可以通过HTML导入母版页,然后HTML和Master页面相关联,更改HTML页的时候,Master会自动同步修改,然而,有些时候我们不需要 ...

  5. Android自定义控件4--优酷菜单的菜单键及细节补充

    在上篇文章中实现了优酷菜单执行动画,本文接着完善已经实现的动画功能 本文地址:http://www.cnblogs.com/wuyudong/p/5915958.html ,转载请注明源地址. 已经实 ...

  6. Core Data的一些常见用法

    一.简介 Core Data是一个纯粹的面向对象框架,其本质就是一个ORM(对象关系映射:Object Relational Mapping),能以面向对象的方式操作SQLite数据库.在实际开发中绝 ...

  7. 使用FragmentTabHost+TabLayout+ViewPager实现双层嵌套Tab

    大多数应用程序都会在底部使用3~5个Tab对应用程序的主要功能进行划分,对于一些信息量非常大的应用程序,还需要在每个Tab下继续划分子Tab对信息进行分类显示. 本文实现采用FragmentTabHo ...

  8. 适配ios10（iTunes找不到构建版本）

    前两天上架App遇到一个比较神奇的问题,打包好的项目使用Application Loader上传成功,但是在iTunes里面却找不到构建版本,App的活动页面也没有相应的版本. 之前了解IOS10对用 ...

  9. -bash: .bash_profile: command not found

    今天有一同事安装了ORACLE后,在切换账号时遇到错误提示"-bash: .bash_profile: command not found".如下所示 [root@GLETestL ...

  10. SQL性能优化：如何定位网络性能问题

    一同事跟我反馈他遇到了一个SQL性能问题,他说全表只有69条记录,客户端执行耗费了两分多钟,这不科学呀.要我分析一下原因并解决.我按照类似表结构,构造了一个案例,测试截图如下所示 这个表有13800K ...