把一团“云”当成一个“点”去校准，会发生什么？

先别急着看公式。假设你拍到的是一团边缘柔和的云，但修图软件被告知：“画面中央应该是一颗针尖大小、亮度恰好为 1 的星。”软件又被允许调节每个镜头的亮度和相位。它会怎么做？它会把“云不像星”的差异，错误地解释成镜头参数不对，然后用一套错误的镜头校正去处理原图。最后，云真的被修得更像一颗星。

01｜论文到底在问什么？

这不是纯粹的玩具问题。射电自校准经常用一组 CLEAN 点源分量近似天空；若校准源其实略有分辨，或扩展辐射没有被充分建模，就会出现类似的不匹配。论文作者也提到，HERA-19 调试中曾以点源作为银河系中心的启动模型，不过真实流程随后还会通过自校准继续完善模型。

02｜零基础：干涉仪究竟测到了什么？

单个天线接收天空电波；两根天线组成一条基线（baseline），比较两路信号的相关性。每条基线在某个时刻测得一个复数，这就是可见度（visibility）。它同时包含振幅和相位，可理解为天空亮度分布在某个空间频率上的一笔采样。

随着地球自转，同一条基线会在 uv 平面上画出一条轨迹。把很多基线、很多时刻的可见度组合起来，再做傅里叶反变换，才得到我们熟悉的天空图像。因此，图像里的一个伪结构，往往来自很多基线上的系统偏差共同叠加。

03｜校准：本来是修正仪器，为什么会改造天空？

每根天线都有复增益 g：振幅部分像音量旋钮，相位部分像时间延迟。方向无关校准假设一根天线在整个视场只有一个增益，然后寻找一组增益，使模型可见度尽量接近观测可见度。

如果 ℳ 正确，G 主要吸收仪器和传播路径的误差。如果 ℳ 错了，最小二乘求解器仍然要把 ℛ 和 ℳ 拉近，于是 G 可能同时吸收两类东西：真正的仪器误差，以及“真实天空与错误模型之间的差”。

校准后还要用求得的 G 去“除掉”观测中的增益，得到校正可见度矩阵 𝒮。第二个等号是原文给出的逐元素写法：𝒢 = gg^H 是可见度增益矩阵，⊙ 是 Hadamard（逐元素）乘积，𝒢^⊙-1 是它的 Hadamard 逆。问题在于：若 G 已经混进了天空模型误差，这个除法就不再只是修仪器，而是在对天空做一场有偏的重塑。

04｜作者怎么把看不见的误差“拆出来”？

论文把两种方法叠在一起：扰动分析负责给出近似解析式，外推（extrapolation）负责生成每条基线的、没有 uv 采样旁瓣污染的伪影图，并作为经验基准验证理论。

方法 A：扰动分析

作者假设估计增益只比真实增益多一个较小偏差 Δg，把校准目标在真实解附近展开，并保留一阶项。这样可以得到增益矩阵的近似解析表达式。它的价值不在于逐像素完美预测，而在于揭示误差由哪些阵列几何项、源宽度项和基线比例项控制。

方法 B：外推

东西向阵列的不同 uv 轨迹可以看成同一参考轨迹的缩放。作者把某条基线的响应外推到整个规则 uv 平面，再做 FFT。这样得到的图不受稀疏采样产生的常规旁瓣干扰，更容易看清“校准本身制造了什么”。外推是逐基线进行的。

三组结果，不要混淆

05｜实验设置：作者究竟模拟了什么？

σ = 0.02° 大约是该配置在 1.45 GHz 分辨率的十倍，所以真实源确实被充分分辨，绝不是一个“差不多可以当点源”的对象。这种刻意放大的模型失配，便于作者看清系统误差的结构。

06｜五个核心结果

07｜全图谱：14 幅图逐张读懂

Figure 1WSRT 几何矩阵 Φ：之后所有“例外基线”的根源

怎么看：每个格子是基线 uv 轨迹相对参考轨迹的整数缩放因子 |φ_pq|。靠近对角线的数通常较小，对应较短基线；右侧两列的特殊数值结构后来会解释 Figure 13 在约 1224 m 后的误差反弹。作用：把阵列物理几何转成式 (28)-(38) 中可计算的缩放系数。

Figure 2所有基线的增益与逆增益响应

怎么看：上三角是增益，下三角是逆增益；绿、蓝、红分别代表 S1、S2、S3。三色主瓣通常重合，说明解析推导抓住了中心结构；旁瓣有时分离，说明高阶扰动与倒数近似会影响细节。作用：第一次全阵列验证式 (28)、(29) 与纯外推结果。

Figure 3两个代表性基线的放大图

怎么看：基线 (0-13) 的理论、精确倒数与外推结果吻合较好；基线 (9-10) 的曲线差异更明显。作用：提醒读者“主瓣总体吻合”并不等于每条基线的旁瓣都能被一阶理论准确预测。

Figure 4全部基线的校正可见度与校准后天空

怎么看：上三角给出校正可见度曲线，下三角给出相应天空切片。短基线常出现更高、更窄的中心响应，也更容易出现三种方法不一致。作用：直接把“增益误差”连接到最终科学图像中的源形态改变。

Figure 5五类基线案例：理论从准确到失效

怎么看：五列对应论文选出的代表基线，上排是可见度、下排是天空响应。由左到右可以看到解析结果从较好复现 S3，到低估、过估或无法复现。作用：为 Figure 6(c) 和 Table 2 的五级准确度分类提供具体例子。

Figure 6理论系数统计与基线准确度分类

怎么看：(a) 与 (b) 描述 K̃_pq,rs^-2 的近似均值和标准差，(c) 把基线按理论预测质量分为 0-4 类。作用：把“某些基线理论不准”从观察提升为可分类、可与阵列几何对应的结论。

Figure 7S3、S2、S1 的峰值矩阵

怎么看：三幅矩阵分别给出 S3、S2、S1 的归一化峰值并在 2 处截断。总体上 |φ| 越小，峰值越大，但有受几何矩阵影响的例外。作用：建立“峰值可作为理论是否可靠的代理量”，并支持短基线更受系统误差影响。

Figure 8校准后源的归一化 FWHM

怎么看：S1 和 S2 多数基线的宽度接近原源；S3 中短基线更明显地给出小于 1 的 FWHM。作用：与 Figure 7 的峰值一起证明源被压窄、变得更点状。

Figure 9每条基线的估计总流量

怎么看：色标从 1 到 2，超过 2 的结果被截断。S3 对多数基线接近 1 Jy；少数基线的峰值增加没有被宽度下降完全抵消。作用：验证式 (38)-(39)：总流量被校准拉向模型中的 1 Jy 标尺。

Figure 10天线数 N 对峰值、宽度、流量的影响

怎么看：实线、虚线、点线分别是峰值、FWHM、流量；红、蓝、绿对应 S1、S2、S3。随着 N 从 4 增到 14，没有统一单调趋势。作用：说明“伪影亮度强依赖天线数”不适用于这个单个、中心、部分建模源的特殊场景。

Figure A1两种逆增益近似为什么选择式 (29)

怎么看：红线是纯外推真值，蓝线是式 (28) 的数值精确倒数，绿线是式 (29)，青线是式 (A9)，洋红线是中间量 f。A9 在原点附近更贴近蓝线，但式 (29) 更能覆盖全局形状且便于解析傅里叶变换。作用：解释正文为何采用较简单的一阶倒数近似。

原文两张表与算法框

Table 1实验参数与外推单元测试

不仅列出 WSRT、A、σ、频率、时角、500 时间步和 14 根天线，还用 b₀₁=[0,-696]^T 验证外推值与传统 StEFCal 都为 3.78×10^-17。

Table 2五级基线分类规则

用 S2 与 S3 的峰值区间给每条基线赋 0-4 标签，对应 Figure 6(c) 的颜色类别。

Algorithm 1逐基线外推伪代码

在规则 uv 网格逐点构造 ℛ 与 ℳ、调用任意校准器，并把结果写入目标基线的外推增益网格。

08｜进阶：从公式看“为什么会叠成扩展伪源”

真实天空被写成一个归一化圆高斯：

A 是积分后的总流量，σ 控制宽度。原文的模型天空是 I_ℳ(s) = δ(s)，即视场中心的 1 Jy 点源（论文式 16）。

对增益做一阶扰动并利用东西向阵列的几何缩放后，每条基线 pq 的外推增益近似成为：

式 (28) 是一个常数项，加上一串“在 uv 域仍为高斯”的项。傅里叶反变换后，常数变成中心 δ 点源，高斯仍变成高斯。原文完整式子是：

其中 σ_pq,rs = σφ_rs/φ_pq；φ_pq 是基线 pq 的 uv 轨迹相对于参考轨迹的缩放因子。

校正可见度等于观测可见度乘以增益的倒数。图像域里，乘法对应卷积，所以真实高斯源会与“中心点源 + 一群同心高斯幽灵”卷积，最终仍是一群高斯的组合：

其中 K_N = (2N - 1)/N，K_pq,rs = √(φ_pq² + φ_rs²)/φ_pq，分别对应原文式 (29) 与式 (36)。式 (37) 近似的是每条基线 pq 在校准后的天空亮度分布，因此不同基线得到的 I_𝒮pq(s) 并不完全相同。

它同时解释了三件事：

完整公式图谱：正文 (1)-(47) 与附录 (A1)-(A9)

公式 (1)-(4)：校准问题、校正操作与真实观测

式号	名称	在推导中的作用
(1)	方向无关校准目标	寻找天线增益矩阵 G，使观测 ℛ 与模型预测 GℳGH 的 Frobenius 距离最小。
(2)	Hadamard 形式	把天线增益外积 ggH 记为可见度增益矩阵 𝒢，将矩阵乘法改写为逐元素乘法。
(3)	校正可见度	用估计增益的逆作用于观测，得到校正矩阵 𝒮；这是模型误差重新进入图像的入口。
(4)	真实观测方程	把观测写成真实增益 Gt 作用在真实天空可见度 𝒳 上。

公式 (5)-(10)：外推几何与扰动展开

式号	名称	在推导中的作用
(5)	真实天空与模型的联合目标	把式 (4) 代入校准目标，显式暴露真实天空 𝒳、模型 ℳ、真实增益与估计增益之间的关系。
(6)	东西向阵列几何矩阵 Φ	用带符号的 φpq 编码每条基线相对参考 uv 轨迹的缩放。
(7)	轨迹线性变换	在东西向阵列条件下，把一般变换 Ppq 简化为 φpqb。
(8)	增益扰动模型	令 g=gt+Δg，并展开可见度增益矩阵 𝒢。
(9)	扰动后的校准目标	把式 (8) 代入式 (5)，使未知量变成 Δg。
(10)	扰动项展开	把天空差异与一阶、二阶增益扰动项分离，为忽略 ΔgΔgH 做准备。

公式 (11)-(26)：线性化求解与高斯真天空/点源模型

式号	名称	在推导中的作用
(11)	一阶线性化目标	定义 𝒟=𝒳-ℳ 并忽略二阶小量，只保留可闭式求解的一阶扰动。
(12)	Δg 闭式解	用分块矩阵逆把 Δg 与其共轭写成由 E、F 和真实增益决定的近似解。
(13)	分块矩阵 A-F	定义式 (12) 所需的六个矩阵块，压缩后续代数。
(14)	真实天空高斯源	规定真实天空是中心、总流量 A、标准差 σ 的二维圆高斯。
(15)	高斯归一化常数 B	令 B=2σ²π，使式 (14) 对天空积分后得到总流量 A。
(16)	点源校准模型	把模型天空设为中心 δ 函数，即 1 Jy 点源。
(17)	点源模型可见度矩阵	中心单位点源在所有基线上产生全 1 的模型矩阵 ℳ=11H。
(18)	真实高斯的逐基线可见度	二维高斯的傅里叶变换仍为高斯，衰减由基线长度 bpqTbpq 决定。
(19)	分块矩阵简化	利用 ℳ=11H 把 A-D 四个块化成 N、I 与 11H 的组合。
(20)	E gt	把各天线连接基线上的高斯可见度求和，得到式 (12) 的右端第一部分。
(21)	F gt	给出共轭/列方向对应的第二部分；对当前对称实数场景与式 (20) 结构相似。
(22)	分块矩阵伪逆	给出式 (19) 的 Moore-Penrose 伪逆，使 Δg 可显式计算。
(23)	可见度增益的一阶近似	取 gt=1 并忽略 ΔgΔgH，得到 11H 加两个一阶项。
(24)	逐基线增益的高斯和	证明 gpq 是常数项加多条基线高斯项的线性组合。
(25)	常数系数	规定 cpq,0=1/N。
(26)	逐项系数	依据 rs 与目标基线 pq 是否共享天线，给 cpq,rs 五种取值。

公式 (27)-(38)：索引集合、增益幽灵、逆增益与校准后天空

式号	名称	在推导中的作用
(27)	五个索引集合	严格定义式 (26) 的 I1-I5，消除“共享哪根天线”的歧义。
(28)	外推增益	把式 (24) 映射到参考 uv 平面；φrs²/φpq² 控制每个高斯项的宽度。
(29)	逆增益近似	以一阶 Maclaurin 近似构造 gpq-1；附录解释为何选它。
(30)	真实 uv 轨迹上的逆增益	把式 (29) 采样回基线 pq 的轨迹，供实际校正可见度使用。
(31)	增益伪影图	反傅里叶变换后得到中心 δ 项加一组同心高斯；这是扩展 ghost 的解析来源。
(32)	逆增益伪影图	同样得到中心点项与同心高斯和，但高斯项符号相反并带 A、B 的倒数。
(33)	观测逐基线可见度	在真实增益设为 1 时，把真实观测写成 A 乘 uv 高斯。
(34)	外推观测可见度	映射到目标基线参考坐标后，基线依赖被消去。
(35)	校正可见度	将式 (29) 与式 (34) 相乘；A 被抵消，留下不同宽度的高斯组合。
(36)	Kpq,rs	把 φpq 与 φrs 合并成控制校正后高斯宽度的比例因子。
(37)	校准后天空亮度	式 (35) 的反傅里叶变换；直接显示真实源与同心 ghost 卷积后的最终形态。
(38)	积分流量与中心峰值	给出每条基线校准后源的总流量约为 1，以及中心峰值的解析表达式。

公式 (39)：用峰值和 FWHM 估计高斯总流量

(39)

二维圆高斯流量公式

F=πP·FWHM²/(4ln2)。作者用它从峰值 P 和宽度估计 Figure 9 的总流量；若曲线含负翼，直接按单高斯估计会与数值积分略有差别。

公式 (40)-(43)：为什么外推能复现传统校准

(40)	外推观测矩阵	重写 Algorithm 1 第 9 行，突出 Prs∘Ppq-1 的坐标映射。
(41)	待化简的复合变换	单独抽出几何映射，便于把 b 替换为真实轨迹 bpq。
(42)	复合变换回到 brs	证明目标点位于真实 pq 轨迹时，所有 rs 项恰好回到各自真实 uv 坐标。
(43)	恢复普通观测矩阵	式 (40) 因而退化为传统校准使用的 ℛ，解释 Table 1 单元测试为何完全相等。

公式 (44)-(47)：真实源趋向点源时的四个极限

(44)	观测天空极限	B→0 时，真实高斯趋向 Aδ(s)。
(45)	增益伪影极限	增益图同样趋向 Aδ(s)，在可见度域表现为由 A 缩放的常数增益。
(46)	逆增益极限	逆增益图趋向 A-1δ(s)。
(47)	校正天空极限	A 与 A-1 相消，最终只剩单位点源 δ(s)；若 A=1，ghost 完全消失。

附录公式 (A1)-(A9)：逆增益近似的两条路线

式号	名称	在推导中的作用
(A1)	中间量 fpq	把增益按 A 归一化并减 1，使逆增益可写成 1/(1+f) 的幂级数问题。
(A2)	f 的弱上下界	证明其上界小于 1、下界大于 -2；并指出经验上大部分区域满足 |f|<1。
(A3)	逆增益的精确 f 表达	把目标写为 [A(1+f)]-1，为 Maclaurin 展开做准备。
(A4)	第一条路线的起点	重复精确倒数形式，随后在 |f|≪1 的主瓣区域截断级数。
(A5)	一阶 Maclaurin	使用 1/(1+f)≈1-f。
(A6)	用 g 表示的一阶倒数	把 f 换回 g，得到 2A-1-A-2g；代入式 (28) 即正文式 (29)。
(A7)	第二条路线：二维二阶展开	直接在 u=v=0 对逆增益做二阶 Maclaurin；一阶导数和混合项在对称性下为零。
(A8)	二阶系数 zpq	把二阶曲率写成 A、B、c 与 φ 比例的求和。
(A9)	局部二次近似	得到 A-1+zpqbTb；原点附近更准，但全局形状和可解析变换不如式 (29) 合用。

09｜这篇论文没有证明什么？

好的论文解读不仅要讲结论，还要标出边界。作者的理论和实验都很干净，但现实观测更脏、更复杂。

10｜读论文时最容易卡住的词

Visibility

可见度。两天线信号相关得到的复数采样，是天空亮度傅里叶变换上的一点。

Baseline

基线。两根天线之间的矢量间距；长度和方向决定采样哪一种空间尺度。

Gain

增益。每根天线对信号振幅与相位的复数响应；校准就是估计它。

Sky model

天空模型。算法事先相信天空中有哪些源、位置在哪、亮度和形态如何。

Ghost

幽灵源。由校准系统误差产生的伪天体结构，不是天空中真正存在的辐射。

Extended emission

扩展发射。角尺度大于仪器点扩散函数、不能视作点源的辐射结构。

Extrapolation

外推。利用阵列几何把各基线轨迹映射到参考轨迹，在完整 uv 平面构造逐基线响应。

Geometry matrix Φ

几何矩阵。对东西向阵列，记录各基线 uv 轨迹相对参考轨迹的缩放因子。

Perturbation

扰动。把估计增益写成真实增益加一个小偏差，再通过一阶展开得到可解释的近似式。

11｜用三个问题检查你是否真的懂了