1 files changed, 71 insertions, 0 deletions

diff --git a/src/calendar.h b/src/calendar.h
new file mode 100644
index 0000000..c87531d
--- /dev/null
+++ b/src/calendar.h
@@ -0,0 +1,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