求时角的Python程序

求时角的Python程序一般涉及到天文学计算，是通过计算某个天体在地球上的位置与观测者位置的差距来计算得到的。以下是一个简单的求时角的Python程序：

import math

# 经度、纬度和观测时间
longitude = 116.3
latitude = 39.9
time = 14.0

# 计算儒略日
def calculate_jd(year, month, day):
    if month <= 2:
        year -= 1
        month += 12
    a = int(year / 100)
    b = 2 - a + int(a / 4)
    jd = int(365.25 * (year + 4716)) + int(30.6001 * (month + 1)) + day + b - 1524.5
    return jd

# 计算黄道坐标
def calculate_ecliptic_coordinate(jd):
    d = jd - 2451545.0
    w = 282.9404 + 4.70935E-5 * d
    e = 0.016709 - 1.151E-9 * d
    m = 356.0470 + 0.9856002585 * d
    l = w + m
    dl = (1.9146 - 0.004817 * math.sin(math.radians(m)) - 0.000014 * math.sin(math.radians(2 * m))) * math.sin(math.radians(l)) + (0.019993 - 0.000101 * math.sin(math.radians(m))) * math.sin(math.radians(2 * l)) + 0.000290 * math.sin(math.radians(3 * l))
    sl = l + dl
    sra = math.atan2(math.cos(math.radians(sl)) * math.sin(math.radians(23.4393)), math.cos(math.radians(sl)))
    dec = math.asin(math.sin(math.radians(sl)) * math.sin(math.radians(23.4393)))
    return sra, dec

# 计算赤道坐标
def calculate_equatorial_coordinate(sra, dec, jd):
    t = (jd - 2451545.0) / 36525.0
    eps = 23.4393 - 0.0130041667 * t
    ra = sra
    dec = dec
    ha = math.radians((time - longitude / 15.0) * 15.0 - ra)
    alt = math.asin(math.sin(dec) * math.sin(math.radians(latitude)) + math.cos(dec) * math.cos(math.radians(latitude)) * math.cos(ha))
    az = math.acos((math.sin(dec) - math.sin(alt) * math.sin(math.radians(latitude))) / (math.cos(alt) * math.cos(math.radians(latitude))))
    if math.sin(ha) > 0:
        az = math.pi * 2 - az
    return alt, az

# 主程序
year = 2021
month = 10
day = 1
jd = calculate_jd(year, month, day)
sra, dec = calculate_ecliptic_coordinate(jd)
alt, az = calculate_equatorial_coordinate(sra, dec, jd)
ha = (time - longitude / 15.0) * 15.0 - sra
ha_hour = ha / 15.0
print("时角为：%.2f" % ha_hour)

注意：这个程序只是一个简单的示例，实际应用中还需要考虑更多的因素，比如大气折射等。如果你需要更精确的计算，建议使用专业的天文学计算软件或库。