2024华中杯A-1

import math

# 北纬30°35′转换为十进制度数
latitude = 30.58333333

# 大气层外层每月的太阳能辐射强度（I0）
monthly_i0 = {
    1: 1405,
    2: 1394,
    3: 1378,
    4: 1353,
    5: 1334,
    6: 1316,
    7: 1308,
    8: 1315,
    9: 1330,
    10: 1350,
    11: 1372,
    12: 1392
}

# 计算太阳高度角
def calculate_sun_elevation_angle(latitude, declination):
    return 90 - math.fabs(latitude - declination)

# 估算太阳赤纬角
def estimate_declination(day_of_year):
    declination = 23.45 * math.sin(math.radians(360 * (284 + day_of_year) / 365))
    return declination

# 计算太阳直射辐射的函数
def calculate_radiation(i0, panel_tilt, sun_elevation):
    angle_difference = abs(panel_tilt - sun_elevation)
    cosine_factor = math.cos(math.radians(angle_difference))
    if cosine_factor < 0:
        return 0, 0
    max_intensity = i0 * cosine_factor
    total_energy = max_intensity * effective_sunshine_hours
    return max_intensity, total_energy

# 每天有效日照时间
effective_sunshine_hours = 6

# 光伏板的倾角
tilt_angles = [20, 40, 60]

# 计算并存储每月15日的数据
monthly_data = []
for month in range(1, 13):
    day_of_year = (month - 1) * 30 + 15  # 估算每月15日是年中的第几天
    declination = estimate_declination(day_of_year)
    sun_elevation = calculate_sun_elevation_angle(latitude, declination)
    i0 = monthly_i0[month]  # 获取当前月份的大气层外层太阳能辐射强度
    radiation_data = []
    for tilt in tilt_angles:
        max_intensity, total_energy = calculate_radiation(i0, tilt, sun_elevation)
        radiation_data.append((max_intensity, total_energy))
    monthly_data.append(radiation_data)

# 输出结果
for month, data in enumerate(monthly_data, start=1):
    print(f"2025年{month}月15日的数据：")
    for tilt, (max_intensity, total_energy) in zip(tilt_angles, data):
        if max_intensity > 0:
            print(f"倾角 {tilt}°时的最大太阳直射强度：{max_intensity:.2f} W/m^2, "
                  f"太阳直射辐射总能量：{total_energy:.2f} Wh/m^2")
    print()