本节摘要:校准是把仪器从数据里剥离的过程:带宽定标抚平频响起伏,流量定标把数据绑到绝对标度,增益与相位解消除天线与大气引入的复数误差,自校准用目标自身模型迭代收敛。校准的哲学是承认数据是"天空与仪器的卷积",先把仪器测清楚。本节按校准流程走完全链,并给出绝对流量换算的演算。
上一节的流水线吐出了规整的数据网格,但网格里既有天空,也有仪器:带宽内有滤波器的纹波、各通道增益不齐;天线增益随仰角变形;相位以分钟量级漂移;极化两路存在串扰。读谱的人若不知道这些指纹,会把仪器的纹波当成源的谱特征、把相位漂移当成源的视运动。校准的全部工作是先量出仪器的指纹,再从数据里减去它——前两节的观测计划(定时切定标源)就是为这一节准备的原料。
真实的校准按固定工序走。第一步带宽定标:用噪声二极管或定标源测出各通道的相对响应,抚平频带内的纹波,让每条通道有统一的单位。第二步流量定标:观测特性已知的定标源(如仙后座 A、天鹅座 A 这类标准烛光),把数据绑到绝对流量标度——从此数据单位是央斯基,而不是任意的仪器计数。第三步增益与相位解:对相位定标源的每个时间片解出各天线的复数增益(幅度与相位),插值到目标观测,把目标数据里对应的误差除掉。第四步(干涉阵专属)极化定标:解开两路极化的相位差与泄漏,为偏振测量铺路(本章第六节回收)。
不是亮源都能当定标源。流量定标源要满足三条件:亮(信噪比好)、稳(几年量级不变)、已知(有长期监测的模型频谱)。北天的仙后座 A 因超新星遗迹的自行在缓慢变暗,它的模型连年变率都写进了标度表——定标源是被"供养"的国民资产。相位定标源另有一套资格:离目标近(几度内,大气与仪器误差才能共用)、点状(结构简单才能反解增益)、够亮(几十秒积分就能测出解)。近到什么程度算近?经验上几度以内共用同一条大气路径,相位解的插值才站得住。
# 绝对流量换算:从定标源模型到目标流量 # 标准烛光的模型频谱:S = S0 * (nu/1GHz)^alpha S_casA_1GHz = 2723.0 # Jy,仙后座 A 在 1 GHz 的近似值 alpha_casA = -0.77 # 它的谱指数(还在缓慢变暗,精确模型带年变项) f = 1.42e9 S_model = S_casA_1GHz * (f/1e9)**alpha_casA print(f"仙后座 A 模型流量 @1.42GHz = {S_model:.0f} Jy") # 观测得到的目标与定标源功率计数(同一配置下) P_cal, P_target = 5000.0, 0.9 # 任意单位 S_target = P_target / P_cal * S_model print(f"目标功率是定标源的 {P_target/P_cal*1000:.2f} 千分之一 -> 目标流量 {S_target*1000:.1f} mJy")
输出(真实运行结果):
仙后座 A 模型流量 @1.42GHz = 1657 Jy 目标功率是定标源的 0.18 千分之一 -> 目标流量 298.3 mJy
这就是绝对定标的全部算术:目标与定标源的功率计数同配置同带宽相除,乘上定标源的模型流量。比例式简单,资格严格——定标源的模型流量错了零点几成,下游所有流量、所有光度、所有质量(如第六章用氢线质量公式称出来的氢质量)全部继承这个偏差。
经典校准有天花板:相位解在两次定标切换之间靠插值,插值误差直接糊在图像上。自校准的思路激进:干涉数据对源的亮度分布和天线增益是可交换的——先用手头最好的天图当模型,反解各天线的增益修正;用新增益改进天图;再迭代。天线有几十面时约束远多于自由度,迭代通常在几轮内收敛,动态范围能提升一个数量级以上。
它的风险同样鲜明:如果初始模型错了,自校准会"勤奋地"把模型误差吸进增益解里,越校越像错误模型——俗话叫"校出了你开始时的想象"。防线有三条:初始模型用物理上可信的成分(点源骨架);只解幅度与相位中漂移快的分量;每轮迭代后人工检查解的合理性。自校准是强效药,剂量全靠流程控制。
⚠️ 易错点:把校准当成"让数据好看"的美化步骤。校准的唯一目标是可复现的物理单位与可解释的误差预算;一切让图"更好看"却拿不出定标依据的操作,都该在发表前撤回。
校准做不到完美,成熟观测员给每类残差记账:绝对流量标度百分之几(定标源模型的不确定度)、相位残余几度(切换间隔内插值误差)、带通残余百分之几(纹波未抚平部分)。误差预算的意义在比较:目标物理量若只需要两成的精度,这些残余无关痛痒;若要测百分之几的谱指数,标度残余就成了主要误差源——该多切几次流量定标而不是多熬积分。校准的功夫最终体现在这张预算表上,而不是图像的观感上。
定标误差会沿科学链条逐级放大:流量标度差百分之五,两频点间的谱指数就差百分之五;由谱指数推电子能谱、由能谱推磁场、再由磁场推星系演化的结论,误差层层传递。经验法则是:越靠下游的物理结论,对上游定标的相对误差越敏感。所以高精度项目(精密谱指数、偏振角、微弱信号统计)在计划阶段就要把"定标循环加密"写进提案,而不是事后指望数据处理补救——上游错的部分,下游补不回来。
仪器指纹剥净了。最后一站处理最嚣张的对手:比宇宙源亮百亿倍的人为干扰,以及保护收听权的频谱制度。