3.5 增益演算：分辨率与流量密度估算
3.5 增益演算：分辨率与流量密度估算
本节摘要：本章前四节的工程原理，收束成两组演算：分辨角从肉眼排到黑洞照相机的一张对数表（一点二二乘波长除口径，阵列换基线）；观测预算从目标流量反推积分时长（辐射计方程的实战用法）。两组演算全部可复算，是观测提案可行性论证的标准动作，也是读懂任何性能表的钥匙。
把设备排上一杆秤
前三节拆完了天线、性能参数与接收链路，第四节把分辨率交给了基线。现在把它们放回同一杆秤上。衡量任何射电设施，最后都落到两个数：分辨角（看多细）与灵敏度（听多弱）。本节先造分辨角的秤，再造灵敏度的秤，然后用一道完整的观测预算题把两杆秤串起来——这套动作就是观测员写提案前的心算作业。
图 3-5：从肉眼到黑洞照相机的分辨角对数表
第一组演算：分辨角对数表
口径与波长决定单机分辨角，基线替代口径决定阵列分辨角。下面这张表覆盖了射电天文九十年里所有标志性设备，数字一行一个公式，全可手算复核。
# 分辨角对数表：theta = 1.22 * lambda / D import math c = 299792458.0 rad2as = 206264.806 rows = [("肉眼（5mm 瞳孔，550nm）", 5.5e-7, 0.005), ("雷伯后院锅（9.75m@160MHz）", c/160e6, 9.75), ("百米全可动锅（21cm）", 0.211, 100.0), ("天眼有效口径（300m@21cm）", 0.211, 300.0), ("甚大阵 A 构型（36km@2cm）", 0.02, 36000.0), ("洲际甚长基线（8000km@3.5mm）", 0.0035, 8e6), ("黑洞照相机（1e4km@1.3mm）", 0.0013, 1e7)] print(f"{'设备':<28}{'分辨角':>14}") for name, lam, D in rows: th_as = 1.22 * lam / D * rad2as unit = f"{th_as:.4g} arcsec" if th_as >= 0.01 else f"{th_as*1000:.4g} mas" print(f"{name:<28}{unit:>14}")
输出（真实运行结果）：
设备 分辨角 肉眼（5mm 瞳孔，550nm） 27.68 arcsec 雷伯后院锅（9.75m@160MHz） 4.836e+04 arcsec 百米全可动锅（21cm） 531 arcsec 天眼有效口径（300m@21cm） 177 arcsec 甚大阵 A 构型（36km@2cm） 0.1398 arcsec 洲际甚长基线（8000km@3.5mm） 0.1101 mas 黑洞照相机（1e4km@1.3mm） 0.03271 mas
这张表值得读三遍。第一遍看射电的先天劣势：三百米口径的天眼在二十一厘米波段的分辨角，仍比肉眼差六倍有余；雷伯的后院锅更是差上千倍。第二遍看技术补偿：同在厘米波段，甚大阵用三十六公里基线把分辨角压到零点一四角秒，反超肉眼两百倍。第三遍看极限演出：把基线铺满地球、把波长压到一点三毫米，分辨角压进毫角秒以内——从地球上看月球表面的一张信用卡。第六章将看到，这级分辨率正是给黑洞阴影拍照的门票。
第二组演算：观测预算
分辨率的答案决定了"用哪台设备"，灵敏度的答案决定"测多久"。预算题的标准形式：SEFD 已知的目标、给定的带宽与效率，探测到目标流量需要多长积分？把辐射计方程反过来用即可。
# 观测预算：探测 1 mJy 目标（5 倍信噪比）需要多久 import math SEFD = 300.0 # Jy，某中口径天线 eta = 0.88 # 系统效率 npol = 2 # 双极化 snr = 5.0 S_target = 1e-3 # Jy sigma_req = S_target / snr print(f"要求噪声 σ = {sigma_req*1000:.1f} mJy") for dnu_MHz in (64, 400, 4000): dnu = dnu_MHz * 1e6 tau = (SEFD / (eta * sigma_req))**2 / (npol * dnu) print(f"带宽 {dnu_MHz:5d} MHz -> 积分 {tau:7.0f} s ≈ {tau/3600:5.2f} h")
输出（真实运行结果）：
要求噪声 σ = 0.2 mJy 带宽 64 MHz -> 积分 22699 s ≈ 6.31 h 带宽 400 MHz -> 积分 3632 s ≈ 1.01 h 带宽 4000 MHz -> 积分 363 s ≈ 0.10 h
结论一目了然：老带宽下探测一毫央斯基要熬六个小时，换四吉赫的宽带接收机，六分钟收工。接收机的一步升级抵得上上百小时的机时——这就是观测界对宽带化趋之若鹜的算术根源。把信噪比再验一遍，预算闭环：
# 预算复核：算出给定时长的实际信噪比 import math SEFD, eta, npol = 300.0, 0.88, 2 dnu = 64e6 for S_mJy, tau in [(10.0, 600), (1.0, 600), (1.0, 22700)]: sigma = SEFD / (eta * math.sqrt(npol * dnu * tau)) print(f"S = {S_mJy:5.1f} mJy，积分 {tau:6d} s -> 信噪比 {S_mJy/1000/sigma:.1f}")
输出（真实运行结果）：
S = 10.0 mJy，积分 600 s -> 信噪比 8.1 S = 1.0 mJy，积分 600 s -> 信噪比 0.8 S = 1.0 mJy，积分 22700 s -> 信噪比 5.0
十毫央斯基的源十分钟就有八倍信噪比；一毫央斯基的源十分钟只有零点八倍——根本看不见；熬满六个多小时，五倍信噪比达标。提案里的机时申请就是这么算出来的。
变式演练：三道真实场景
预算能力要练到见题拆题。三道变式，先自己心算，再对照答案。
变式一（选设备）：要在 21cm 分开相距 2 角秒的一对中性氢云核。 天眼（177 角秒）能吗？甚大阵 C 构型（基线 3.4km，约 15 角秒）呢？ 答案：天眼不行（177 远大于 2），甚大阵 C 构型勉强（15 略大于 2）， 得用 B 构型（11km，约 4.7 角秒）或 A 构型（1.2 角秒）。 变式二（降要求）：若目标只有 0.5 mJy，带宽 64 MHz，还要几小时？ 答案：要求噪声降到 0.1 mJy，时间翻四倍 -> 约 25 小时。 合理的对策是换宽带接收机而不是熬机时。 变式三（换算直觉）：基线 100km、波长 6cm，条纹角间距多少？ 答案：0.06/1e5 弧度 = 1.24 角秒。 想把分辨率再翻倍，基线加长到 200km，或波长换 3cm。
💡 关键直觉：观测预算是乘除法不是感觉。目标流量除以信噪比要求得噪声上限，噪声上限反解积分时间——三步算完，这台设备能不能干这活、值不值这个机时，立刻见分晓。
本节要点回顾
- 分辨角一杆秤：一点二二乘波长除口径，阵列把口径换成基线，从十三度到万分之三毫角秒；
- 百年百万倍：分辨率改善全部来自口径与基线的延伸，没有捷径；
- 预算三步：目标流量除信噪比、辐射计方程反解、积分时长换算成机时；
- 带宽最划算：六十四兆赫到四吉赫，六小时变六分钟，升级接收机胜过熬时间；
- 选设备先算后选：分辨目标尺度、灵敏度、带宽三数对齐，性能表才能变成决策。
增益工程全章收工：信号已经收集、放大、合成了。下一章面对最后一道工序——数据里的相位漂移与人为干扰，降噪调校。