Function calculate_direct_rate

backend/orbit.py:309–348

直连链路速率 (Mbps)，由链路预算引擎驱动（含大气衰减）。 替代原 base_rate × exp(-d/k) 经验公式， 使用 FSPL/EIRP/G-T/C-N0 推导物理上准确的可达速率； 通过仰角与地面站降雨率叠加 ITU-R P.838/P.618 大气衰减。

(sat_pos, gs, data_type=None)

def calculate_direct_rate(sat_pos, gs, data_type=None):
    """直连链路速率 (Mbps)，由链路预算引擎驱动（含大气衰减）。

    替代原 base_rate × exp(-d/k) 经验公式，
    使用 FSPL/EIRP/G-T/C-N0 推导物理上准确的可达速率；
    通过仰角与地面站降雨率叠加 ITU-R P.838/P.618 大气衰减。
    """
    gs_lat = gs["lat"]
    gs_lon = gs["lon"]
    gs_alt_m = gs.get("alt", 0.0)
    sat_alt_m = sat_pos.get("alt", 500e3)

    # 精确计算仰角（用于大气衰减路径折算）
    look = compute_look_angles(
        sat_lat=sat_pos["lat"],
        sat_lon=sat_pos["lon"],
        sat_alt_m=sat_alt_m,
        gs_lat=gs_lat,
        gs_lon=gs_lon,
        gs_alt_m=gs_alt_m,
    )
    elevation_deg = look["elevation_deg"]
    dist_km = look["slant_range_km"] if look["slant_range_km"] > 0 else _slant_range_km(
        sat_pos, {"lat": gs_lat, "lon": gs_lon, "alt": gs_alt_m}
    )

    antenna_type = gs.get("antenna_type", "Ka")
    # 地面站可配置降雨率场景（默认晴天）
    rainfall_rate_mm_h = gs.get("rainfall_rate_mm_h", 0.0)
    gs_altitude_km = gs_alt_m / 1000.0

    result = link_budget_direct(
        distance_km=dist_km,
        antenna_type=antenna_type,
        data_type=data_type,
        elevation_deg=elevation_deg,
        rainfall_rate_mm_h=rainfall_rate_mm_h,
        gs_altitude_km=gs_altitude_km,
    )
    return max(result.achievable_rate_mbps, 5.0)


def calculate_relay_rate(sat_pos, geo_pos, gs, data_type=None):