1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
#ifndef _CALENDAR_H_
#define _CALENDAR_H_
#include <ctime>
#include <fstream>
#include <iostream>
#include <math.h>
#include <string>
#include <vector>
using namespace std;
// 哈尔滨南岗地方经纬度
#define M_PI 3.14159265358979323846
#define LONGITUDE 126.66854 * M_PI / 180
#define LATITUDE 45.75996 * M_PI / 180
typedef struct tm Date;
string radToDMS(int lowerBound, int upperBound, double radians);
class Julian {
private:
static double d[23][5];
double dt_ext(int y, double jsd);
public:
// 计算儒略日
double getJulianDay(time_t time);
// 计算儒略千年数
double getJulianKiloYear(time_t time);
//儒略千年数转时间戳
time_t kiloYearToTime(double t, int year);
//计算力学时与世界时之差,传入年份
double delta_t(int y);
};
class parameter {
private:
// 计算地球日心黄经
double get_longitude(vector<double> l, double t);
// 计算太阳地心黄纬
double get_latitude(vector<double> b, double t);
double get_R(vector<double> r, double t);
// 转换FK5误差，返回弧度制
double delta_FK5(double L, double B, double T);
//获取章动有关角
vector<double> get_angle(double T);
//章动修正
double nutation(double T);
//光行差修正，返回弧度制
double aberration(double R);
// 获取地日运行参数，L为地球日心黄经，B为地球日心黄纬，R为地日距离
vector<vector<double>> get_parameters(double t);
public:
pair<double, double> sun_longitude(double t);
double get_epsilon(double T);
// 计算时角，返回弧度制
double getHourAngle(double julianKiloYear, double longitude,
double latitude, double alpha);
};
extern Julian julian;
extern parameter p;
extern double epsilon;
#endif
|
