4.2 频谱后端：通道化与数据流水线
4.2 频谱后端：通道化与数据流水线
本节摘要：数字化之后的信号进入频谱后端：多相滤波器组把宽带切成成千上万条窄通道，逐通道累计成时间片，再按科学目标分流——谱线模式攒通道分辨率，脉冲星模式攒时间分辨率并在流水线里完成消色散与折叠。通道宽、时间分辨率与数据量构成三角账，本节用可复算的演算把这笔账算清。
从一维波形到二维数据立方
上一节排好了计划，本节看数据本身长什么样。天线输出经接收机数字化后是一维电压流，但科学问题几乎都问在二维上：频率（谱线在哪个频率、什么形状）与时间（何时到达、如何变化）。频谱后端的工作就是把一维流切成二维网格：横轴频率、纵轴时间，每个格子里是平均功率。切法不唯一，切得不合适，第三章攒下的灵敏度与第二章讲的色散修正都会在流水线里漏掉。这一站是全章工程味最浓的一站，却是读懂数据格式的必经之路。
图 4-2：观测数据流水线：从电压流到科学产品
通道化：切多窄才算够
通道宽度的账有两个方向。往窄切：谱线观测要分辨线形，速度分辨率换算成频率分辨率是速度比光速乘静止频率；往宽切：每条通道的噪声随通道变窄而变大（平方根规律又来了），太窄的通道在带宽定标后会留下平台状起伏。最优通道宽通常取谱线线宽的三分之一到五分之一——足以刻画线形，又不至于让通道噪声失控。
# 通道账：频率分辨率、速度分辨率与色散拖尾 import math c = 299792.458 # km/s bw = 400e6 for nch in (1024, 4096, 8192, 16384): dnu = bw / nch dv = c * dnu / 1420.40575e6 # 21cm 处的速度分辨率 print(f"{nch:6d} 通道 -> 通道宽 {dnu/1e3:7.1f} kHz，21cm 速度分辨率 {dv:5.1f} km/s") # 色散拖尾：单通道内因色散造成的脉冲展宽 # 拖尾（毫秒）= 8.3 x DM x 通道宽（GHz）x 频率（GHz）的三次方倒数 DM = 300.0 nch = 8192 dnu_GHz = bw / nch / 1e9 nu = 1.3 smear_ms = 8.3 * DM * dnu_GHz * nu**-3 print(f"DM={DM}、1.3GHz、通道宽 {dnu_GHz*1e6:.0f} kHz -> 通道内拖尾 {smear_ms:.2f} ms")
输出（真实运行结果）：
1024 通道 -> 通道宽 390.6 kHz，21cm 速度分辨率 82.5 km/s 4096 通道 -> 通道宽 97.7 kHz，21cm 速度分辨率 20.6 km/s 8192 通道 -> 通道宽 48.8 kHz，21cm 速度分辨率 10.3 km/s 16384 通道 -> 通道宽 24.4 kHz，21cm 速度分辨率 5.2 km/s DM=300、1.3GHz、通道宽 49 kHz -> 通道内拖尾 0.06 ms
（两段演算一次运行输出。）第一段账告诉谱线观测员：分辨十公里每秒的星系谱线至少要八千通道。第二段账告诉脉冲星观测员：色散量三百的信号在千赫级通道内每通道拖尾不足零点一毫秒，比典型毫秒采样还短——通道够细，色散就只表现为通道间的时差，可在处理时对齐；通道太粗，拖尾直接糊在通道内部，神仙难救。
消色散与折叠：脉冲星流水线的心脏
第二章算过色散延迟：DM 三十的信号，零点四与一点四吉赫之间差七百多毫秒。流水线的对策分两步。消色散：把各通道的时间轴按色散公式逐条平移对齐，使脉冲在所有通道同时到达，再逐时间片求和——色散量未知时（巡天情形）就对一列候选色散量逐个试验，看哪个值让脉冲最尖。折叠：已知周期的脉冲星，把长时间数据按周期相位折起来累加，噪声按平方根稀释、脉冲按幅度累加，几万圈下来轮廓浮现。这一对操作就是脉冲星天文的数据引擎，第六章第五节的毫秒脉冲星钟表学全靠它供数。
脉冲星巡天流水线（简化） 电压流 -> 通道化 -> 检波落盘（毫秒片） -> 对一列候选色散量逐个消色散 -> 对每个候选做快速折叠与周期搜索 -> 候选打分（信噪比、色散一致性、重复性） -> 人工或机器复核 -> 新脉冲星档案
数据量的三角账
通道数、时间分辨率、位深三者相乘就是数据率，上一节算过一条四百兆赫链路八比特每秒上百吉字节。流水线的所有"聪明"都在跟这笔账搏斗：相干消色散保真但费算力，非相干消色散便宜但留拖尾；巡天模式先粗筛再精测，把窄通道深采样只花在候选上。理解了三角账，就读懂了为什么各台站的后端配置单长得千差万别——不是技术偏好，是各自科学目标在这笔账上的不同解。
⚠️ 易错点：把"通道多"当"数据好"。通道超过科学需要，只放大存储与处理成本，还让每通道信噪比变差、定标曲线起伏；通道数是配平项，不是荣誉项。
后端配置速查
把三角账的常见解法列成速查表，写提案时对号入座：
科学目标 通道策略 时间策略 数据率压力 银河氢巡天 宽带加万级通道 秒级积分 中 谱线缺源搜索 中带宽加窄通道 分钟级积分 低 脉冲星到达时间 中带宽加微秒采样 毫秒级折叠 高 暂现源巡天 宽视场加毫秒采样 实时触发落盘 极高
表里的每一行都对应后端的一种固化模式，切换模式通常要重新分配通道化资源——这正是观测计划与后端配置必须一起设计的原因：提案里的通道数不是随手填的，它牵动整条流水线的排班。
通道泄漏与平台起伏
后端还有两个隐蔽的坑。通道泄漏：多相滤波器组的通道响应带旁瓣，强谱线旁边的通道会沾到邻居的亮度，定标后表现为谱线两侧的"裙边"；平台起伏：定标后的带通在通道间留有百分之几的纹波残余，宽而平的真实谱容易被纹波伪装出结构。对策分别是给通道化加窗（压旁瓣）与提高定标循环的密度——两个坑都属于"处理得当就消失、处理不当就写进论文"的类型，写数据处理说明时值得主动交代。
本节要点回顾
- 二维网格：后端把一维电压流切成频率乘时间的功率网格，一切科学分析都发生在这张网格上；
- 通道宽三帕斯卡：通道窄于线宽三分之一到五分之一，兼顾线形与通道噪声；
- 速度换算：速度分辨率等于光速乘通道宽除以观测频率，21cm 处四十千赫约每秒十公里；
- 色散两步走：通道够细则拖尾不糊，通道间时差在流水线里对齐；巡天对候选色散量逐个试；
- 折叠放大信噪比：已知周期按相位折叠，噪声平方根稀释、脉冲线性累加；
- 三角账定配置：通道、时间分辨率、位深相乘即数据率，配置单是科学目标的账面解。
数据组织完毕，下一站把仪器自己留在数据里的指纹剥掉：带宽响应、流量标度与相位增益的校准。