summaryrefslogtreecommitdiffstats
path: root/src/calendar.h
blob: c87531d185e9516b88d053aca2db9dc7676a9c29 (plain) (blame)
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
#ifndef _CALENDAR_H_
#define _CALENDAR_H_


#include <iostream>
#include <string>
#include <ctime>
#include <vector>
#include <fstream>
#include <math.h>
using namespace std;

typedef struct tm Date;

extern double pi;
extern double delta;

class Julian
{
    private:
        static double d[23][5];

        double dt_ext(int y,double jsd);

        //计算力学时与世界时之差,传入年份
        double delta_t(int y);
    public:
        // 计算儒略日
        double getJulianDay(time_t time);

        // 计算儒略千年数
        double getJulianKiloYear(time_t time);

        //儒略千年数转时间戳
        time_t kiloYearToTime(double t,int year);
};


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 nution(double T);

        //光行差修正,返回弧度制
        double aberration(double R);

        // 获取地日运行参数,L为地球日心黄经,B为地球日心黄纬,R为地日距离
        vector<vector<double>> get_parameters(double t);

    public:

        double sun_longitude(double t);
};

#endif