要求:
(1)用java编写
(2)输入是 经纬度 和 日期
(3)输出是日出时间和日落时间
急求代码 感谢
2条回答 默认 最新
- qq_35660280 2017-06-18 03:29关注
java写的我找不到 但是我找到了c语言的 要改成java 应该不难 我c语言版本的也能用 有需要的自己改一下吧
下面我分享一下c语言版本的 我测试 误差在几分钟之内 很好用
#include
#include
#include
#include "math.h"
#include "time.h"define M_PI 3.14159265358979323846
static double RAD = 180.0 * 3600 / M_PI;
static double richu;
static double midDayTime;
static double dawnTime;
static double jd;
static double wd;
static double jd_degrees;
static double jd_minutes;
static double jd_seconds;
static double wd_degrees;
static double wd_minutes;
static double wd_seconds;
/*************************- 儒略日的计算 *
- @param y 年 *
- @param M 月 *
- @param d 日 *
- @param h 小时 *
- @param m 分 *
- @param s秒 *
-
@return int
***************************/
static double timeToDouble(int y, int M, double d)
{
// double A=0;
double B = 0;
double jd = 0;//设Y为给定年份,M为月份,D为该月日期(可以带小数)。
//若M > 2,Y和M不变,若 M =1或2,以Y–1代Y,以M+12代M,换句话说,如果日期在1月或2月,则被看作是在前一年的13月或14月。
//对格里高利历有 :A = INT(Y/100) B = 2 - A + INT(A/4)
//对儒略历,取 B = 0
//JD = INT(365.25(Y+4716))+INT(30.6001(M+1))+D+B-1524.5 (7.1)
B = -13;
jd = floor(365.25 * (y + 4716)) + floor(30.60001 * (M + 1)) + B + d - 1524.5;
return jd;
}
static void doubleToStr(double time, char *str)
{
double t;
int h, m, s;t = time + 0.5; t = (t - (int)t) * 24; h = (int)t; t = (t - h) * 60; m = (int)t; t = (t - m) * 60; s = (int)t; sprintf(str, "%02d:%02d:%02d", h, m, s);
}
/****************************
- @param t 儒略世纪数 *
- @return 太阳黄经 *****************************/ static double sunHJ(double t) { double j; t = t + (32.0 * (t + 1.8) * (t + 1.8) - 20) / 86400.0 / 36525.0; // 儒略世纪年数,力学时 j = 48950621.66 + 6283319653.318 * t + 53 * t * t - 994 + 334166 * cos(4.669257 + 628.307585 * t) + 3489 * cos(4.6261 + 1256.61517 * t) + 2060.6 * cos(2.67823 + 628.307585 * t) * t; return (j / 10000000); }
static double mod(double num1, double num2)
{
num2 = fabs(num2);
// 只是取决于Num1的符号
return num1 >= 0 ? (num1 - (floor(num1 / num2)) * num2) : ((floor(fabs(num1) / num2)) * num2 - fabs(num1));
}
/********************************- 保证角度∈(-π,π) *
- @param ag
- @return ag ***********************************/ static double degree(double ag) { ag = mod(ag, 2 * M_PI); if (ag <= -M_PI){ ag = ag + 2 * M_PI; } else if (ag>M_PI){ ag = ag - 2 * M_PI; } return ag; }
/***********************************
*- @param date 儒略日平午 *
- @param lo 地理经度 *
- @param la 地理纬度 *
- @param tz 时区 *
-
@return 太阳升起时间
*************************************/
double sunRiseTime(double date, double lo, double la, double tz)
{
double t, j, sinJ, cosJ, gst, E, a, D, cosH0, cosH1, H0, H1, H;
date = date - tz;
// 太阳黄经以及它的正余弦值
t = date / 36525;
j = sunHJ(t);
// 太阳黄经以及它的正余弦值
sinJ = sin(j);
cosJ = cos(j);
// 其中2*M_PI*(0.7790572732640 + 1.00273781191135448*jd)恒星时(子午圈位置)
gst = 2 * M_PI * (0.779057273264 + 1.00273781191135 * date) + (0.014506 + 4612.15739966 * t + 1.39667721 * t * t) / RAD;
E = (84381.406 - 46.836769 * t) / RAD; // 黄赤交角
a = atan2(sinJ * cos(E), cosJ);// '太阳赤经
D = asin(sin(E) * sinJ); // 太阳赤纬
cosH0 = (sin(-50 * 60 / RAD) - sin(la) * sin(D)) / (cos(la) * cos(D)); // 日出的时角计算,地平线下50分
cosH1 = (sin(-6 * 3600 / RAD) - sin(la) * sin(D)) / (cos(la) * cos(D)); // 天亮的时角计算,地平线下6度,若为航海请改为地平线下12度
// 严格应当区分极昼极夜,本程序不算
if (cosH0 >= 1 || cosH0 <= -1){
return -0.5;// 极昼
}
H0 = -acos(cosH0); // 升点时角(日出)若去掉负号 就是降点时角,也可以利用中天和升点计算
H1 = -acos(cosH1);
H = gst - lo - a; // 太阳时角
midDayTime = date - degree(H) / M_PI / 2 + tz; // 中天时间
dawnTime = date - degree(H - H1) / M_PI / 2 + tz;// 天亮时间
return date - degree(H - H0) / M_PI / 2 + tz; // 日出时间,函数返回值
}
void input_gJD(int c[]){int i;
printf("Enter the degree of longitude(form: 40 40 (means 40°40′40″)):");
for (i = 0; i<3; i++){scanf("%d", &c[i]);
}
for (i = 0; i < 3; i++){
printf("< %d>", c[i]);
};
printf("\n");
}
//输入经度
void input_gWD(int c[]){
int i; printf("Enter the degree of latitude(form: 40 40 40 (means 40°40′40″)):"); for (i = 0; i<3; i++){ scanf("%d",&c[i]); } for (i = 0; i < 3; i++){ printf( "< %d>",c[i]); }; printf("\n");
}
//输入经度
void input_date(int c[]){int i; printf("Enter the date (form: 2009 03 10):"); for (i = 0; i<3; i++){ scanf("%d",&c[i]); }
}
//输入日期
int main()
{time_t curtime; struct tm * timeinfo; //int i, j, k; char timestr[20]; int year; int month; int day; //int hour; //int min; //int sec; int tz; int c[3]; //FILE *fp; //上海东经121度29分,北纬31度14分 //jd_degrees = 121; //jd_seconds = 28; //wd_degrees = 31; //wd_seconds = 14; //乌鲁木齐东经87度35分,北纬43度48分 //jd_degrees=87; //jd_seconds=35; //wd_degrees=43; //wd_seconds=48; //广州 E113°16' N23°06' input_gJD(c); jd_degrees = c[0]; jd_minutes = c[1]; jd_seconds = c[2]; input_gWD(c); wd_degrees = c[0]; wd_minutes = c[1]; wd_seconds = c[2]; tz = 8; //得到当前时间 time(&curtime); timeinfo = localtime(&curtime); //printf("The current date/time is: %s", asctime(timeinfo)); //step 1 jd = -(jd_degrees + jd_minutes / 60+jd_seconds/3600) / 180 * M_PI; wd = (wd_degrees +wd_minutes/60+ wd_seconds / 3600) / 180 * M_PI; //printf("将输出经纬度的语句注释掉\n"); //printf("jd=%f\r\n", jd); //printf("wd=%f\r\n", wd);; //step 2 input_date(c); year = c[0]; month = c[1]; day = c[2]; //year = timeinfo->tm_year + 1900; //month = timeinfo->tm_mon + 1; //day = timeinfo->tm_mday; //hour = timeinfo->tm_hour; //min = timeinfo->tm_min;
// sec = timeinfo->tm_sec;
richu = timeToDouble(year, month, day) - 2451544.5;
//printf("mjd:%f\r\n", richu);
printf("%d年%d月%d日\r\n", year, month, day);
for (int i = 0; i < 10; i++){
//时区在此处运用,可以根据经纬度算出 ,也是整型
richu = sunRiseTime(richu, jd, wd, tz / 24.0);// 逐步逼近法算10次
}doubleToStr(richu, timestr); printf("日出时间 %s\r\n", timestr); doubleToStr(midDayTime + midDayTime - richu, timestr); printf("日落时间 %s\r\n", timestr); //printf("richu:%f,midDayTime:%f\r\n", richu, midDayTime); /*doubleToStr(midDayTime, timestr); printf("中天时间 %s\r\n", timestr); doubleToStr(dawnTime, timestr); printf("天亮时间 %s\r\n", timestr); doubleToStr(midDayTime + midDayTime - dawnTime, timestr); printf("天黑时间 %s\r\n", timestr); fp = fopen("richu.txt", "w"); if (fp == NULL){ printf("file open error\r\n"); } /* //char fst[]={"日出 日落\n"}; char str[50]; char m2012[12]={31,29,31,30,31,30,31,31,30,31,30,31}; //2012 //fputs(fst,fp); for(i=0;i<12;i++){ for(j=1;j<=m2012[i];j++){ year=2012; month=i+1; day=j; richu = timeToDouble(year,month,day*1.0) - 2451544.5; for (k = 0; k < 10; k++){ richu = sunRiseTime(richu, jd, wd, tz/24.0);// 逐步逼近法算10次 } sprintf(str,"%d-%02d-%02d ",year,month,(int)day); doubleToStr(richu,timestr); strcat(str,timestr); strcat(str," "); doubleToStr(midDayTime + midDayTime - richu, timestr); strcat(str,timestr); strcat(str,"\n"); fputs(str,fp); } } */ //fclose(fp); //printf("Hello world!\n"); return 0;
}
本回答被题主选为最佳回答 , 对您是否有帮助呢?解决 无用评论 打赏 举报
悬赏问题
- ¥15 Unet采样阶段的res_samples问题
- ¥60 Python+pygame坦克大战游戏开发实验报告
- ¥15 R语言regionNames()和demomap()无法选中中文地区的问题
- ¥15 Open GL ES 的使用
- ¥15 我如果只想表示节点的结构信息,使用GCN方法不进行训练可以吗
- ¥15 QT6将音频采样数据转PCM
- ¥15 下面三个文件分别是OFDM波形的数据,我的思路公式和我写的成像算法代码,有没有人能帮我改一改,如何解决?
- ¥15 Ubuntu打开gazebo模型调不出来,如何解决?
- ¥100 有chang请一位会arm和dsp的朋友解读一个工程
- ¥15 查询优化:A表100000行,B表2000 行,内存页大小只有20页,运行时3页,设计两个表等值连接的最简单的算法