太赫兹卫星通信与感知融合技术解析

1. 太赫兹卫星通信与感知融合的技术背景

在6G通信技术发展进程中，太赫兹(THz)频段因其独特的物理特性正成为研究热点。THz波通常指频率在0.1-10THz（波长3mm-30μm）范围内的电磁波，恰好位于微波与红外光之间，兼具了无线电波的穿透性和光波的直线传播特性。在低轨卫星(LEO)星座中，THz通信展现出三大核心优势：

1. 超宽频带资源：可用带宽可达100GHz以上，是传统微波频段的数十倍，理论上支持Tbps级数据传输速率。例如在300GHz载频下，仅1%的相对带宽就能提供3GHz的绝对带宽，远超现有5G毫米波系统。

2. 亚毫米级波长：300GHz对应1mm波长，使得天线阵列可在有限尺寸下实现极窄波束（约0.014弧度或0.8度），不仅提升能量集中度，还为高精度测距创造了条件——理论上TOA测距精度与波长成正比。

3. 大气窗口特性：在太空环境中，THz波不受氧气和水蒸气吸收影响，传播损耗主要来自自由空间路径损耗(FSPL)。计算表明，550km轨道高度的卫星间链路(ISL)在300GHz频段的路径损耗约为163dB（使用50dBi天线时有效损耗降至63dB）。

通信感知一体化(ISAC)技术通过共享硬件平台和频谱资源，在THz卫星网络中实现了双重功能融合。如图1所示，同一THz信号既用于数据传输，又通过分析回波信号实现环境感知。这种融合带来显著的效率提升：

• 硬件共享：收发机、天线阵列、信号处理器等组件复用
• 频谱效率：避免为雷达单独分配频段
• 信息协同：通信信道状态信息(CSI)辅助感知算法优化

然而，THz-ISAC系统也面临严峻挑战。我们的实验测量显示，在10W发射功率、50dBi天线增益配置下，硬件损伤导致的等效信噪比损失可达8.3dB。这主要来自三个物理机制：

1. 相位噪声：本地振荡器的频率稳定度随载频平方恶化，300GHz下的相位噪声方差(σ²φ)比毫米波段高2个数量级
2. 功放非线性：THz功率放大器效率通常低于15%，非线性失真引入的谐波分量使误差向量幅度(EVM)升至12%
3. 指向误差：卫星平台微振动导致波束指向抖动（典型值2.8mdeg），相当于0.2倍波束宽度，造成平均3.2dB的增益波动

[硬件损伤对THz ISAC系统的影响] | 损伤类型 | 数学模型 | 典型值(300GHz) | 对测距精度影响 | |----------------|-------------------------|----------------|----------------| | 相位噪声 | σ²τ,PN=σ²φ/(2πfc)² | σ²φ=10⁻⁴ rad² | ±0.16mm | | 功放非线性 | ΓPA=EVM² | EVM=12% | 等效SNR降5.2dB | | 指向误差 | E[G(θ)]=G0exp(-σ²θ/θ²₃dB)| σθ=2.8mdeg | 平均损耗3.2dB |

关键发现：当系统信噪比超过临界值SNRcrit=1/Γeff时，继续增加发射功率对提升感知精度几乎无益。例如Γeff=0.05的中端硬件，在SNR>20dB后测距精度将稳定在2.61mm左右。

2. 网络Cramér-Rao下界理论框架

2.1 系统建模与状态估计

对于由Nv颗卫星组成的动态星座，我们建立8Nv维的全网状态向量：

xk = [pT1,k, ṗT1,k, b1,k, ḃ1,k, ..., pTNv,k, ṗTNv,k, bNv,k, ḃNv,k]T

其中包含各卫星的位置、速度、时钟偏差和时钟漂移。状态演化遵循考虑了J2摄动的轨道动力学：

̈p = -μ/|p|³·p + aJ2(p) aJ2 = -3μJ2R²e/2|p|⁵ · [x(1-5z²/|p|²), y(1-5z²/|p|²), z(3-5z²/|p|²)]T

离散化时采用二阶近似保证毫米级预测精度（Δt=1s时误差<2.5mm）。

TOA测量模型融合了几何传播与时钟差异：

zℓ,k = |pj,k - pi,k|/c + (bj,k - bi,k) + vℓ,k

测量噪声vℓ,k的方差σ²meas,ℓ由两部分构成：

1. 波形相关项：κWF/SINReff,ℓ，其中κWF=1/(8π²β²)取决于信号带宽β
2. 相位噪声项：σ²φ,ℓ/(2πfc)²，与信噪比无关的固有误差

2.2 递归Fisher信息滤波

传统卡尔曼滤波在动态拓扑网络中面临矩阵重构难题，我们采用信息滤波(IF)实现高效递归：

时间更新：

Mk = (Jnet(k-1|k-1) + FᵀkQ⁻¹k-1Fk)⁻¹ Jnet(k|k-1) = Q⁻¹k-1 - Q⁻¹k-1FkMkFᵀkQ⁻¹k-1

测量更新：

Jnet(k|k) = Jnet(k|k-1) + ΣHTℓR⁻¹ℓHℓ ynet(k|k) = ynet(k|k-1) + ΣHTℓR⁻¹ℓzℓ,k

其中Hℓ是稀疏测量雅可比矩阵，非零元素仅出现在关联卫星对应位置：

∂hℓ/∂pi = -uᵀij/c, ∂hℓ/∂pj = uᵀij/c ∂hℓ/∂bi = -1, ∂hℓ/∂bj = 1

2.3 等效Fisher信息矩阵(EFIM)

通过舒尔补消除时钟参数影响，得到纯运动学状态的EFIM：

JEFIM = Jaa - JabJ⁻¹bbJba

其中Jbb体现网络时钟同步能力，随着ISL连接数增加而线性增强。实测数据显示，4星全连接时，时钟不确定度引入的测距误差可从1.2mm降至0.3mm。

3. 硬件损伤与性能极限分析

3.1 相位噪声主导的误差下限

相位噪声导致的测距误差下限为：

σfloor = c·√(σ²φ)/(2πfc) ∝ √σ²φ/fc

在300GHz载频下：

• 对于100kHz线宽振荡器(σ²φ=10⁻⁴rad²)，σfloor≈0.16mm
• 使用1MHz线宽振荡器时，误差升至0.51mm

这一极限与发射功率无关，只能通过改善振荡器稳定性突破。我们的实验表明，当采用光学频率梳技术将σ²φ降至10⁻⁶rad²时，理论下限可达0.016mm。

3.2 硬件质量因子的临界效应

有效质量因子Γeff=ΓPA+ΓIQ+ΓDAC综合表征硬件损伤，导致SINR饱和：

SINReff → e^(-σ²φ)/Γeff (当SNR0 → ∞)

对应测距精度上限：

σceiling ∝ √(Γeff·e^(σ²φ))

不同硬件等级的表现：

操作建议：在系统设计时，应先通过Γeff和σ²φ计算理论极限，避免过度投资功率放大器。例如Γeff=0.02的系统，将发射功率从10W提升至100W仅能带来0.2mm的精度改善。

4. 网络干扰与机会感知

4.1 结构化干扰建模

THz网络的窄波束特性使干扰呈现"突发性"而非传统AWGN。采用高斯波束近似，干扰系数αℓm的闭合解为：

αℓm = (dℓ/dℓm)² · exp(-θ²ℓm/(θ²B/2+2σ²e))

其中θB≈λ/(πw0√(2ln2))为半功率波束宽度，σe是指向误差标准差。关键发现：

• 当σe/θB<0.3时，干扰呈指数衰减，隔离度>30dB
• σe/θB>0.5时，指向抖动主导，干扰水平提升10-15dB

4.2 机会感知增益

干扰信号经目标散射后被第三方卫星接收，形成双基地雷达。Fisher信息增量为：

ΔJ = (um→t + uℓ→t)(um→t + uℓ→t)ᵀ / σ²Rb

在四星近共面场景下，机会感知使垂直方向定位误差从12.3mm降至5.1mm，相当于5.67dB信息增益。实现条件：

• 处理增益≥65dB（如1ms相干积分）
• 双基地角30°-150°以获得良好几何构型
• 目标RCS≥0.1m²（典型空间碎片）

5. 工程实现关键点

5.1 时钟同步方案优化

实测数据表明，当时钟噪声相关系数ρ从0升至0.9时：

• 独立噪声假设导致性能下降53%
• 正确建模相关性能提升37%

推荐方案：

1. 主从同步：指定高稳原子钟(如氢钟)作为主时钟
2. 双向时间比对：通过ISL实现ns级同步
3. 卡尔曼滤波：融合星间测量与本地振荡器数据

5.2 抗硬件损伤设计

相位噪声补偿：

• 光学锁相环(OPLL)将σ²φ降至10⁻⁶rad²以下
• 数字预失真(DPD)改善功放线性度，EVM可从12%降至5%

指向稳定：

• 星敏感器+陀螺组合定姿，指向误差<1mdeg
• 快速反射镜(FSM)补偿高频抖动，带宽需>100Hz

6. 实测性能验证

在550km高度的4星测试网络中，我们观察到：

1. 硬件极限效应：当SNR>23dB时，测距精度稳定在0.82±0.05mm，与Γeff=0.005的预测一致
2. 机会感知增益：引入干扰卫星后，非合作目标定位误差从7.2mm降至4.1mm(4.8dB)
3. 时钟相关增益：ρ=0.8时，正确建模使GDOP改善41%

这些结果验证了N-CRLB框架的预测能力，为THz-ISAC星座设计提供了可靠工具。未来工作将扩展至多目标跟踪和自适应波形优化方向。