下载文档原格式
日出日落时间的计算以及中国常见的日出日落时间以地球中心为原点O,赤道所在平面为XY平面,东经120度指向西经60度为Y轴正方向.球心指向北极为Z 轴正方向.有了Y轴与Z轴就可定X轴的方向(从东经30度指向西经150度)球面方程:X^2 + Y^2 + Z^2 = 1 (设地球直径为1)日出日落时刻圈方程:Y^2 + Z'^2 = 1 (Z'以Z轴作坐标变换,见下面)Z'=Z*sin(β+90) (β为太直射点纬度)求纬度α度时日出时刻.先解出纬度为α度时的X,Y坐标.X=si nα*sin(β+90)*cos(β+90)/(cosβ*cosβ)Y=-SQRT(1-X^2-sinα*sinα) (SQRT为平方根)有了XY坐标,求反正切,得出一个角度值(由于在XY平面,0度在X轴正向,实际的东经120度在Y轴负方向上,即270度角.所以要换算一下,才能得出经度差)实算一下:代入的纬度为30.15度,夏至日时,太阳直射点纬度为23.4333度.算出X=-0.21835,Y=-0.83578,反正切得出-104.64度.计算时假设在Y轴负方向上(即270度或-90度).两者之间相差14.64度,换成时间就是58.56分钟.(计算出的14.64度的含义是指,夏至日那天,当赤道上(北纬0度)东经120度的地方看到日出时,北纬30.15度,东经(120-14.64)度的地方也正好看到日出.)(换句话说:当赤道上东经120度的地方看到日出时,北纬30.15度东经120度的地方日出已经过去58.56分钟了.由于赤道上是昼夜等分的(假设太阳是个点光源),即日出时刻一定在6:00.那么同一经度的北纬30.15度地方,日出时间是5:01:26左右.东经120度10分.比120度还早了40秒钟.所以日出时间为5:00:46)查寿星万年历,在夏至日的日出时间为4:58:07,日落时间19:04:07.实测数据2009年5月1日星期五所有时间为时间 (任意地点日月升落时刻查询)省会城市北京116.46E 39.92N (-14分9秒)日出：05时15分(-1分15秒/日) 日落：19时09分(+1分0秒/日) 长春125.35E 43.88N (+21分23秒)日出：04时31分(-1分25秒/日) 日落：18时41分(+1分11秒/日) 长沙113E 28.21N(-28分0秒)日出：05时48分(-49秒/日)日落：19时03分(+35秒/日)成都104.06E 30.67N (-1时3分45秒) 日出：06时20分(-53秒/日)日落：19时42分(+39秒/日)重庆106.54E 29.59N (-53分50秒)日出：06时12分(-51秒/日)日落：19时31分(+38秒/日)福州119.3E 26.08N (-2分48秒)日出：05时26分(-46秒/日)日落：18时34分(+31秒/日)贵阳106.71E 26.57N (-53分9秒)日出：06时15分(-46秒/日)日落：19时26分(+32秒/日)广州113.23E 23.16N (-27分4秒)日出：05时54分(-40秒/日)日落：18时55分(+27秒/日)126.63E 45.75N (+26分31秒) 日出：04时21分(-1分30秒/日) 日落：18时40分(+1分17秒/日) 海口110.35E 20.02N (-38分36秒) 日出：06时09分(-35秒/日)日落：19时02分(+22秒/日)杭州120.19E 30.26N (+45秒)日出：05时17分(-53秒/日)日落：18时37分(+39秒/日)合肥117.27E 31.86N (-10分55秒) 日出：05时26分(-56秒/日)日落：18时51分(+42秒/日)呼和浩特111.65E 40.82N (-33分23秒)日出：05时33分(-1分16秒/日) 日落：19时30分(+1分2秒/日)济南117E 36.65N(-12分0秒)日出：05时19分(-1分7秒/日)日落：19时00分(+53秒/日)昆明102.73E 25.04N (-1时9分4秒) 日出：06时33分(-44秒/日)日落：19时39分(+29秒/日)拉萨91.11E 29.97N (-1时55分33秒) 日出：07时13分(-53秒/日)日落：20时33分(+38秒/日)兰州103.73E 36.03N (-1时5分4秒) 日出：06时13分(-1分4秒/日)日落：19时52分(+50秒/日)南昌115.89E 28.68N (-16分26秒)日出：05时36分(-50秒/日)日落：18时52分(+36秒/日)南京118.78E 32.04N(-4分52秒)日出：05时19分(-57秒/日)日落：18时45分(+43秒/日)南宁108.33E 22.84N (-46分40秒) 日出：06时14分(-39秒/日)日落：19时14分(+26秒/日)澳门113.33E 22.13N (-26分40秒) 日出：05时55分(-38秒/日)日落：18时53分(+25秒/日)上海121.48E 31.22N (+5分55秒)日出：05时10分(-55秒/日)日落：18时33分(+41秒/日)阳123.38E 41.8N (+13分31秒) 日出：04时43分(-1分20秒/日) 日落：18时44分(+1分5秒/日)114.48E 38.03N (-22分4秒)日出：05时27分(-1分10秒/日) 日落：19时13分(+56秒/日)台北121.5E 25.05N (+6分0秒)日出：05时19分(-44秒/日)日落：18时24分(+29秒/日)太原112.53E 37.87N (-29分52秒)日出：05时35分(-1分9秒/日) 日落：19时20分(+55秒/日)天津117.2E 39.13N (-11分11秒)日出：05时14分(-1分13秒/日) 日落：19时04分(+58秒/日)武汉114.31E 30.52N (-22分45秒)日出：05时40分(-53秒/日)日落：19时01分(+40秒/日)乌鲁木齐87.68E 43.77N (-2时9分16秒) 日出：07时02分(-1分25秒/日) 日落：21时12分(+1分11秒/日) 西安108.95E 34.27N (-44分11秒)日出：05时55分(-1分0秒/日) 日落：19时28分(+47秒/日)西宁101.74E 36.56N (-1时13分2秒) 日出：06时20分(-1分6秒/日) 日落：20时01分(+51秒/日)香港114.1E 22.2N (-23分36秒) 日出：05时52分(-38秒/日)日落：18时50分(+25秒/日)银川106.27E 38.47N (-54分55秒) 日出：05时58分(-1分9秒/日) 日落：19时47分(+56秒/日)州113.65E 34.76N (-25分23秒) 日出：05时36分(-1分2秒/日) 日落：19时10分(+48秒/日)珠海113.52E 22.3N (-25分55秒) 日出：05时54分(-38秒/日)日落：18时52分(+25秒/日)深圳114.07E 22.62N (-23分43秒) 日出：05时52分(-39秒/日)日落：18时50分(+26秒/日)。
C语言存储日期和时间的方法在C语言中，存储日期和时间有多种方法。
下面将介绍三种常用的方法。
1.使用结构体：C语言中的结构体是一种自定义的数据类型，它可以用来将不同类型的数据组合在一起。
我们可以定义一个包含年、月、日、时、分、秒的结构体来存储日期和时间。
```cstruct DateTimeint year;int month;int day;int hour;int minute;int second;};```然后，我们可以通过创建结构体变量来存储具体的日期和时间值。
```cstruct DateTime dt;dt.year = 2024;dt.month = 9;dt.day = 30;dt.hour = 13;dt.minute = 45;dt.second = 0;```2. 使用時間戳（Timestamp）：```c#include <stdio.h>int maireturn 0;```通过计算时间戳，我们可以将其与一些固定时间点进行比较，从而得到具体的日期和时间。
```c#include <stdio.h>int maireturn 0;```3.使用字符串：在C语言中，我们可以用字符串来表示日期和时间，并使用相关的库函数进行转换。
常用的字符表示日期和时间的格式有"YYYY-MM-DDHH:MM:SS"和"MM/DD/YYYYHH:MM:SS"等。
```c#include <stdio.h>int maiint year, month, day, hour, minute, second;printf("年：%d\n", year);printf("月：%d\n", month);printf("日：%d\n", day);printf("时：%d\n", hour);printf("分：%d\n", minute);printf("秒：%d\n", second);return 0;```总结：存储日期和时间的方法有多种，我们可以使用结构体、时间戳或字符串来表示和存储。
c++日出日落计算公式详细（实用版）目录1.C++语言简介2.日出日落计算公式原理3.C++实现日出日落计算公式的代码示例正文【1.C++语言简介】C++是一种通用的高级编程语言，具有高效、灵活、可移植性强等特点。
C++是 C 语言的增强版，支持面向对象编程、泛型编程等现代编程特性，广泛应用于系统编程、应用开发、游戏开发等领域。
【2.日出日落计算公式原理】日出日落计算公式主要用于计算一天中日出和日落的时间。
其原理基于地球自转和公转的运动规律，以及地理经纬度和时区等参数。
具体公式如下：日出时间 = 12 - (当前经度 - 经度 1) / 15 + 当日前一天时间日落时间 = 12 + (当前经度 - 经度 2) / 15 + 当日前一天时间其中，当前经度是指所求地点的经度，经度 1 和经度 2 分别为日出和日落时太阳所在的经度。
【3.C++实现日出日落计算公式的代码示例】以下是一个简单的 C++程序，实现日出日落计算公式的代码示例：```cpp#include <iostream>#include <cmath>#include <iomanip>using namespace std;int main() {double current_longitude, longitude1, longitude2,current_time;cout << "请输入当前地点的经度 (示例:121.47): ";cin >> current_longitude;cout << "请输入日出时太阳所在的经度 (示例:121.47): ";cin >> longitude1;cout << "请输入日落时太阳所在的经度 (示例:121.47): ";cin >> longitude2;cout << "请输入当前时间 (示例:2021-10-01 00:00:00): ";cin >> current_time;double time_difference = fabs(current_time - 24 * 60 * 60);double sunrise_time = 12 - (current_longitude - longitude1) / 15 + time_difference;double sunset_time = 12 + (current_longitude - longitude2)/ 15 + time_difference;cout << "日出时间：" << setw(2) << setprecision(2) << sunrise_time << " 时" << setw(2) << setprecision(2) << sunrise_time- int(sunrise_time) * 60 << " 分" << endl;cout << "日落时间：" << setw(2) << setprecision(2) << sunset_time << " 时" << setw(2) << setprecision(2) << sunset_time- int(sunset_time) * 60 << " 分" << endl;return 0;}```该程序通过接收用户输入的经度、日出日落时太阳所在的经度以及当前时间，计算出日出和日落的时间，并输出结果。
time.h是C/C++中的日期和时间头文件。
用于需要时间方面的函数。
下面分享time.h头文件中几个常用函数的用法：/********time()函数********/函数原型time_t time(time_t *t);函数说明time_t 是long int 类型。
此函数会返回从公元1970年1月1日的UTC时间从0时0 分0秒算起到现在所经过的秒数。
如果t是空指针，直接返回当前时间。
如果t不是空指针，返回当前时间的同时，将返回值赋予t指向的内存空间。
函数返回值成功则返回秒数，失败则返回((time_t)-1)值，错误原因存于errno中。
示例程序#include <stdio.h>#include <time.h>int main(void){time_t sec = time((time_t*)NULL);printf("1970-01-01至现在的秒数为：%d\n", sec);return 0;}输出结果1970-01-01至现在的秒数为：1542377079/********gmtime()函数********/函数原型struct tm *gmtime(const time_t *timep);函数说明gmtime()将参数timep所指的time_t结构中的信息转换成真实世界所使用的时间日期表示方法，然后将结果由结构tm返回。
结构tm的定义为：struct tm{int tm_sec; //代表目前秒数, 正常范围为0-59, 但允许至61 秒int tm_min; //代表目前分数, 范围0-59int tm_hour; //从午夜算起的时数, 范围为0-23int tm_mday; //目前月份的日数, 范围01-31int tm_mon; //代表目前月份, 从一月算起, 范围从0-11int tm_year; //从1900 年算起至今的年数int tm_wday; //一星期的日数, 从星期一算起, 范围为0-6int tm_yday; //从今年1 月1 日算起至今的天数, 范围为0-365int tm_isdst; //日光节约时间的旗标};函数返回值返回结构tm代表目前UTC时间。
昼长夜长日出日落时间的计算公式一、引言昼长夜长日出日落时间的计算公式，是用来预测在任何给定日期和地点，太阳升起和落下的时间的数学公式。
这个公式是基于地球自转和公转的规律，通过计算太阳的位置和地球上观测点的经度、纬度等因素，可以得出准确的日出和日落时间。
二、日出日落的背后科学原理地球自转是指地球绕自身轴线旋转一周所需的时间，约为24小时。
地球公转是指地球绕太阳旋转一周所需的时间，约为365.25天。
由于地球的轴倾斜，所以地球在公转过程中，不同地点的昼夜长度会有所变化。
三、日出时间的计算公式要计算日出时间，需要考虑观测点的经度、纬度以及日期的因素。
根据计算公式，日出时间可以通过以下步骤得出：1. 首先，确定观测点的经度和纬度。
经度表示观测点在东西方向上的位置，纬度表示观测点在南北方向上的位置。
2. 然后，计算出所在时区的标准子午线时间（即当地的标准时间）。
3. 接下来，根据观测点的经度和标准子午线时间，计算出当地的太阳时，即太阳在当地的时间。
4. 最后，通过将太阳时转换为钟表时间，即可得出日出时间。
四、日落时间的计算公式计算日落时间的方法与计算日出时间类似，只是在计算太阳时时，需要考虑地球自转的方向。
具体步骤如下：1. 确定观测点的经度和纬度。
2. 计算出所在时区的标准子午线时间。
3. 根据观测点的经度和标准子午线时间，计算太阳时。
4. 将太阳时转换为钟表时间，即可得出日落时间。
五、实际应用和注意事项昼长夜长日出日落时间的计算公式在日常生活中有广泛的应用，特别是在农业、航海、天文学等领域。
通过准确预测日出和日落的时间，人们可以合理安排农作物的种植和收割时间，帮助船舶和飞机进行导航，以及进行天文观测等。
然而，需要注意的是，这个计算公式是基于理想情况下的理论推导，实际情况可能会受到地形、气候、大气折射等因素的影响，导致计算结果与实际有一定差异。
因此，在实际应用中，需要综合考虑其他因素，并结合观测数据进行修正，以获得更准确的日出和日落时间。
首先，我们根据当前的日期确定太阳在黄道中的位置。
我们知道，太阳大约在3月20日或3月21日直射赤道，而其运动周期是365天。
所以现在太阳在黄道中
的位置是[注1]
黄道平面与赤道平面的夹角被称为黄赤夹角，大约为23.5度。
根据立体几何的知识，太阳当前直射的纬度B
满足
s。
根据这个式子，就可以用计算器的反三角函数求得B
s
接下来，需要计算的是目前你所在纬度的昼长。
由于黄道、赤道和你所在的纬圈均不同在一个平面，所以在纬圈上，太阳照射的部分并不等于纬圈的一半。
这里令2θ为其超出的角度，则有下式
其中B为你当前的纬度。
解出θ后，可以求得昼长等于
等等，这里我们忽略了两点：1. 阳光在大气中的折射会使得我们能见到地平线下方约37′射来的光线。
2. 太阳不是点光源，它是一个约16′圆盘。
考虑这两点，θ还有一个修正量[注2]
在这个基础上，重新计算昼长T。
到这里，最难的问题已经解决了，还需要把当地时间，根据经度L换算到北京时间。
也就是当地正午（太阳最高时）的北京时间为
和日落时间t
那么日出时间t
r
s
由于忽略了地球公转的椭圆轨道（近似为圆），忽略了地球表面的起伏（近似为球）等因素，以上的计算过程和实际值具有一定误差。
注1：为了方便，也可以夏至（或秋分，冬至）作为起算点。
不过需要对应的加上。
注2：这是一个近似的做法，其依据是太阳日出的过程并不长，可以把弧段近似看作线段。
但这一点只能在中低纬度地区得到保证。
对精确解法有兴趣的同学，可以试着动手算一算。
精美文档 1。
怎样用经纬度计算日出日落的时间下面是一种随经纬度变化的日出日落时间计算方法,我成功运用在一智能路灯控制器中,希望对需要的朋友有帮助。
已知：日出日落时太阳的位置h＝-0.833°,要计算地的地理位置,经度Long,纬度G1at,时区zone,UTo为上次计算的日出日落时间,第一次计算时UTo＝180°。
(1)先计算出从格林威治时间公元2000年1月1日到计算日天数day;(2)计算从格林威治时间公元2000年1月1日到计算日的世纪数t,则t＝(day+UTo／360)／36525;(3)计算太阳的平黄径L=280.460+36000.770某t;(4)计算太阳的平近点角G＝357.528+35999.050某t(5)计算太阳的黄道经度λ＝L+1.915某inG+0.020某in(2G);(6)计算地球的倾角ε＝23.4393-0.0130某t;(7)计算太阳的偏差δ＝arcin(inε某inλ);(8)计算格林威治时间的太阳时间角GHA：GHA=UTo-180-1.915某inG-0.020某in(2G)+2.466某in(2λ)-0.053某in(4λ)(9)计算修正值e：e=arco{[inh-in(Glat)in(δ)]/co(Glat)co(δ)}(10)计算新的日出日落时间UT＝UTo-(GHA+Long±e);其中“+”表示计算日出时间,“-”表示计算日落时间;(11)比较UTo和UT之差的绝对值,如果大于0.1°即0.007小时,把UT作为新的日出日落时间值,重新从第(2)步开始进行迭代计算,如果UTo 和UT之差的绝对值小于0.007小时,则UT即为所求的格林威治日出日落时间;(12)上面的计算以度为单位,即180°=12小时,因此需要转化为以小时表示的时间,再加上所在的时区数Zone,即要计算地的日出日落时间为T=UT/15+Zone上面的计算日出日落时间方法适用于小于北纬60°和南纬60°之间的区域,如果计算位置为西半球时,经度Long为负数。
c++日出日落计算公式详细C++是一种高级编程语言，常用于开发各种类型的应用程序。
C++提供了丰富的库和函数，使得编写复杂的计算公式变得更加简单和高效。
在计算日出日落时间方面，可以利用一些基本的天文学公式来进行计算。
日出和日落的计算是基于地球的自转和公转运动的。
地球绕太阳公转的轨道是一个椭圆轨道，而地球自转的轴倾斜23.5度。
这些因素都会影响日出和日落时间的计算。
以下是一个基本的日出日落计算公式的例子：```cpp#include <iostream>#include <cmath>int main() {int year = 2022; // 设置年份int month = 1; // 设置月份int day = 1; // 设置日期// 计算一年中的第几天int N = 275 * month / 9 - ((month + 9) / 12) * (1 + (year - 4 * (year / 4 + 1)) / 3) + day - 30;// 计算地球与太阳之间的平均太阳时角double M = 0.9856 * N - 3.289;// 计算太阳的真太阳时角double L = fmod(M + 1.916 * sin(M_PI * M / 180) + 0.020 * sin(2 * M_PI * M / 180) + 282.634, 360);// 计算太阳的赤纬角double declination = 23.45 * sin(M_PI * (N + 284) / 365);// 计算地球的时差double eqTime = 4 * ((L - 0.06571 * N - 19.44) / 60);// 计算日出时间double sunrise = 12 - (eqTime + 4 * atan(tan(M_PI * declination / 180)) / 15);// 计算日落时间double sunset = 12 - (eqTime - 4 * atan(tan(M_PI * declination / 180)) / 15);std::cout << "日出时间：" << sunrise << "小时" << std::endl; std::cout << "日落时间：" << sunset << "小时" << std::endl;return 0;}```上述代码给出了一个简单的日出日落时间计算的例子。
c++日出日落计算公式详细摘要：一、引言二、C++日出日落计算公式的原理1.地球自转与公转的关系2.太阳光线与地球表面的夹角3.计算日出日落时间的基本公式三、C++实现日出日落计算公式的步骤1.导入相关库和模块2.获取经纬度和海拔高度3.计算当地时区4.计算时角5.计算太阳赤纬6.计算日出日落时间四、C++代码实现1.初始化相关变量2.获取经纬度和海拔高度3.计算当地时区4.计算时角5.计算太阳赤纬6.计算日出日落时间7.输出结果五、结论正文：一、引言在实际生活中，我们有时需要知道具体的日出日落时间，这对于安排户外活动、摄影等有很大的帮助。
C++日出日落计算公式可以帮助我们快速准确地计算当地的日出日落时间。
本文将详细介绍C++日出日落计算公式的原理和实现方法。
二、C++日出日落计算公式的原理要了解C++日出日落计算公式，首先需要了解地球自转与公转的关系以及太阳光线与地球表面的夹角。
在此基础上，我们可以通过计算时角、太阳赤纬等参数，推导出日出日落时间的基本公式。
1.地球自转与公转的关系地球绕太阳公转的同时，自身也在自转。
地球自转使得地球上的每个地点都能经历一次白天和一次黑夜。
地球公转则决定了不同季节太阳光线照射地球的角度不同。
2.太阳光线与地球表面的夹角太阳光线与地球表面的夹角会影响地球上的日照时间。
当太阳光线垂直照射地球某一点时，该点的日出日落时间最短；而当太阳光线平行照射地球时，该点的日出日落时间最长。
3.计算日出日落时间的基本公式日出日落时间可以通过以下公式计算：Δt = 12 * (1 - sin(δ)) / cos(lat)其中，Δt 表示日出日落时间（小时），δ表示太阳赤纬（弧度），lat 表示纬度（弧度）。
三、C++实现日出日落计算公式的步骤要实现C++日出日落计算公式，需要按照以下步骤进行：1.导入相关库和模块我们需要使用C++标准库中的<iostream>、<cmath>等库，以及一些第三方库，如<sunrise-sunset>等。
日出日落时间计算程序(C语言)//日出日落时间计算C语言程序#define PI 3.1415926#include<math.h>#include<iostream>using namespace std;intdays_of_month_1[]={31,28,31,30,31,30,31,3 1,30,31,30,31};intdays_of_month_2[]={31,29,31,30,31,30,31,3 1,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；若非闰年,返回0int 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日到计算日天数dayslong double t_century(int days, long double UTo){return ((long double)days+UTo/360)/36525;}//求格林威治时间公元2000年1月1日到计算日的世纪数tlong 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(P I/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_lon gitude*PI/180)-0.053*sin(4*ecliptic_longi tude*PI/180));}//求格林威治时间的太阳时间角GHAlong double e(long double h, long double glat, long double sun_deviation){return180/PI*acos((sin(h*PI/180)-sin(glat*PI/180)*sin(sun_deviation*PI/180))/(cos(glat*P I/180)*cos(sun_deviation*PI/180)));}//求修正值elong 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(da ys(year,month,date),UTo)),ecliptic_longit ude(L_sun(t_century(days(year,month,date) ,UTo)),G_sun(t_century(days(year,month,da te),UTo)))),glong,e(h,glat,sun_deviation( earth_tilt(t_century(days(year,month,date ),UTo)),ecliptic_longitude(L_sun(t_centur y(days(year,month,date),UTo)),G_sun(t_cen tury(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(day s(year,month,date),UTo)),ecliptic_longitu de(L_sun(t_century(days(year,month,date), UTo)),G_sun(t_century(days(year,month,dat e),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_cent ury(days(year,month,date),UTo))))));result_set(UT,UTo,glong,glat,year,month,d ate);}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)(ri se/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+Z one(glong))))<<" .\n";else cout<<"The time at which the sun rises is "<<(int)(rise/15+Zone(glong))<<":"<<(int) (60*(rise/15+Zone(glong)-(int)(rise/15+Zo ne(glong))))<<" .\n";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))))<<" .\n";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";}//打印结果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_su n(t_century(days(year,month,date),UTo)),e cliptic_longitude(L_sun(t_century(days(ye ar,month,date),UTo)),G_sun(t_century(days (year,month,date),UTo)))),glong,e(h,glat, sun_deviation(earth_tilt(t_century(days(y ear,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)),ecli ptic_longitude(L_sun(t_century(days(year, month,date),UTo)),G_sun(t_century(days(ye ar,month,date),UTo)))),glong,e(h,glat,sun _deviation(earth_tilt(t_century(days(year ,month,date),UTo)),ecliptic_longitude(L_s un(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;}。
