二十四节气

复制import datetime
import math

# 计算太阳黄经（简化版，误差约±2分钟）
# 基于天文算法，输入年份和节气序号（0=小寒, 23=冬至）
def solar_term_jd(year, term_index):
    # 近似公式：节气日期在1月5日附近，每15°黄经对应约15.2天
    base = datetime.date(year, 1, 1).toordinal()
    # 2000年小寒时刻（儒略日）
    j2000 = 2451545.0
    # 每个节气间隔约15.218天
    days_offset = term_index * 15.2184
    # 修正年差（回归年365.2422天）
    year_offset = (year - 2000) * 365.2422
    jd = j2000 + year_offset + days_offset
    return jd

# 儒略日转公历
def jd_to_date(jd):
    jd = jd + 0.5
    z = int(jd)
    f = jd - z
    a = int((z - 1867216.25) / 36524.25)
    b = z + 1 + a - int(a / 4)
    c = b + 1524
    d = int((c - 122.1) / 365.25)
    e = int(365.25 * d)
    g = int((c - e) / 30.6001)
    day = c - e - int(30.6001 * g) + f
    month = g - 1 if g < 14 else g - 13
    year = d - 4716 if month > 2 else d - 4715
    return datetime.date(year, month, int(day))

# 示例：计算2024年立春（节气索引3）
jd = solar_term_jd(2024, 3)
date = jd_to_date(jd)
print(f"2024年立春约在：{date}")  # 输出近似日期

复制package main

import (
	"fmt"
	"math"
	"time"
)

// 计算节气近似儒略日
func solarTermJD(year, termIndex int) float64 {
	j2000 := 2451545.0
	daysOffset := float64(termIndex) * 15.2184
	yearOffset := float64(year-2000) * 365.2422
	return j2000 + yearOffset + daysOffset
}

// 儒略日转时间
func jdToTime(jd float64) time.Time {
	jd += 0.5
	z := int(jd)
	f := jd - float64(z)
	a := int((float64(z) - 1867216.25) / 36524.25)
	b := z + 1 + a - int(float64(a)/4)
	c := b + 1524
	d := int((float64(c) - 122.1) / 365.25)
	e := int(365.25 * float64(d))
	g := int((float64(c) - float64(e)) / 30.6001)
	day := float64(c) - float64(e) - float64(int(30.6001*float64(g))) + f
	month := g - 1
	if g >= 14 {
		month = g - 13
	}
	year := d - 4716
	if month <= 2 {
		year = d - 4715
	}
	return time.Date(year, time.Month(month), int(day), 0, 0, 0, 0, time.UTC)
}

func main() {
	// 计算2024年立春（索引3）
	jd := solarTermJD(2024, 3)
	t := jdToTime(jd)
	fmt.Printf("2024年立春约在：%s\n", t.Format("2006-01-02"))
}

复制// 节气计算（简化版，误差约±2分钟）
function solarTermJD(year, termIndex) {
    const j2000 = 2451545.0;
    const daysOffset = termIndex * 15.2184;
    const yearOffset = (year - 2000) * 365.2422;
    return j2000 + yearOffset + daysOffset;
}

function jdToDate(jd) {
    jd += 0.5;
    const z = Math.floor(jd);
    const f = jd - z;
    const a = Math.floor((z - 1867216.25) / 36524.25);
    const b = z + 1 + a - Math.floor(a / 4);
    const c = b + 1524;
    const d = Math.floor((c - 122.1) / 365.25);
    const e = Math.floor(365.25 * d);
    const g = Math.floor((c - e) / 30.6001);
    const day = c - e - Math.floor(30.6001 * g) + f;
    const month = g < 14 ? g - 1 : g - 13;
    const year = month > 2 ? d - 4716 : d - 4715;
    return new Date(year, month - 1, Math.floor(day));
}

// 示例：2024年立春（索引3）
const jd = solarTermJD(2024, 3);
const date = jdToDate(jd);
console.log(`2024年立春约在：${date.toISOString().slice(0, 10)}`);