2. 星载雷达误差源分类与建模
星载雷达的测量误差并非单一来源,而是由载荷硬件、信号传播环境、卫星平台运动以及随机噪声共同叠加形成。若把所有偏差笼统当作“噪声”处理,标定与补偿就会失去针对性。因此,误差源分类是校准工作的第一步:先判断误差属于哪一类,再决定用确定性模型修正、用统计模型抑制,还是通过平台参数与传播模型联合补偿。
一、四类误差的判别依据
工程上通常按“可重复性”和“物理来源”两个维度划分。系统误差在相同工作条件下重复出现,具有确定规律,例如接收通道幅相不一致、天线阵元位置偏差、定时延迟固定偏移;随机误差每次观测不同,只能用统计量描述,例如热噪声、量化噪声、相位噪声。传播误差来自电磁波穿过大气与电离层时的折射、衰减和时延变化;平台误差则来自卫星姿态、轨道、速度与姿态稳定度引起的几何和相位偏差。
判别时可用一个简单准则:若误差随温度、频率或入射角呈稳定函数关系,优先归为系统误差;若其均值为零、方差有限且与观测时刻无确定关系,归为随机误差;若误差与仰角、电离层总电子含量(TEC)或气象参数强相关,归为传播误差;若误差与姿态角、轨道位置或平台振动频谱相关,归为平台误差。
二、数学建模基础
设雷达第 k 次观测的测量值为 zk,真值为 xk,则误差可写为:
z_k = x_k + δ_sys(k) + δ_rand(k) + δ_prop(k) + δ_plat(k)
其中:
δ_sys —— 系统误差,可建模为确定性函数,如 δ_sys = a·sinθ + b·f + c
δ_rand —— 随机误差,通常假设为零均值高斯白噪声 N(0, σ²)
δ_prop —— 传播误差,电离层项 δ_ion ∝ TEC / f²,对流层项与仰角相关
δ_plat —— 平台误差,由姿态角误差 Δφ、Δθ 与轨道位置误差 Δr 耦合产生
系统误差建模的关键是找到可观测的自变量。例如通道幅相不一致可写成频率 f 和温度 T 的二元多项式;天线阵元位置偏差可写成阵元编号 n 的傅里叶级数。随机误差建模则需确定其概率分布,工程上常用零均值高斯模型,并用 Allan 方差或功率谱密度描述其时间相关性。
传播误差中,电离层引起的群时延近似为:
Δτ_ion ≈ 40.28 × TEC / f² (单位:m,TEC 单位 TECU,1 TECU = 10¹⁶ el/m²)
示例:f = 9.6 GHz,TEC = 20 TECU
Δτ_ion ≈ 40.28 × 20 / (9.6×10⁹)² ≈ 8.7×10⁻¹⁸ s
对应距离误差 ≈ c·Δτ ≈ 2.6 mm
该示例说明:在 X 波段、中等电离层条件下,单频电离层距离误差约毫米级;若频率降至 L 波段,误差将放大约 20 倍,必须用双频或模型校正。对流层误差则常用 Saastamoinen 模型或 Hopfield 模型,以地面气压、温度、湿度和仰角为输入。
三、平台误差的耦合关系
平台误差不是简单的加性项,而是通过几何关系耦合进观测方程。设卫星姿态角误差为 Δφ、Δθ、Δψ,天线波束指向在星体坐标系中的单位矢量为 u,则指向偏差可近似为:
Δu ≈ [Δψ, -Δφ, Δθ]ᵀ × u
若 u = [0, 0, 1]ᵀ(正下方波束),则:
Δu ≈ [Δψ, -Δφ, 0]ᵀ
即方位向偏差由偏航角误差 Δψ 主导,俯仰向偏差由横滚角误差 Δφ 主导。
轨道位置误差 Δr 则直接改变斜距,对测高雷达影响尤为明显。平台振动通常表现为 10~200 Hz 的机械频谱,可通过加速度计测量并建立传递函数,再在信号处理中做相位补偿。
四、分类建模的工程流程
实际工作中,建议按以下顺序推进:先通过地面测试分离系统误差,建立查找表或多项式系数;再利用长时间静态观测估计随机误差的方差与相关时间;然后结合电离层 TEC 地图与气象数据计算传播误差修正量;最后用星敏、陀螺和 GPS 轨道数据计算平台误差。四类误差修正后,用残差均方根(RMS)判断校准效果:若残差接近随机误差理论下限,说明系统误差、传播误差和平台误差已基本消除。
需要强调的是,四类误差之间并非完全独立。例如平台振动会调制传播路径,使传播误差出现周期性分量;温度变化既影响系统误差,也改变随机噪声功率。因此,高精度校准通常采用联合估计:把系统误差系数、传播模型参数和平台姿态误差一起放入状态向量,用卡尔曼滤波或最小二乘统一求解。这样得到的误差模型才具有物理一致性和工程可用性。