1. 太赫兹卫星通信与感知融合的技术背景在6G通信技术发展进程中太赫兹(THz)频段因其独特的物理特性正成为研究热点。THz波通常指频率在0.1-10THz波长3mm-30μm范围内的电磁波恰好位于微波与红外光之间兼具了无线电波的穿透性和光波的直线传播特性。在低轨卫星(LEO)星座中THz通信展现出三大核心优势超宽频带资源可用带宽可达100GHz以上是传统微波频段的数十倍理论上支持Tbps级数据传输速率。例如在300GHz载频下仅1%的相对带宽就能提供3GHz的绝对带宽远超现有5G毫米波系统。亚毫米级波长300GHz对应1mm波长使得天线阵列可在有限尺寸下实现极窄波束约0.014弧度或0.8度不仅提升能量集中度还为高精度测距创造了条件——理论上TOA测距精度与波长成正比。大气窗口特性在太空环境中THz波不受氧气和水蒸气吸收影响传播损耗主要来自自由空间路径损耗(FSPL)。计算表明550km轨道高度的卫星间链路(ISL)在300GHz频段的路径损耗约为163dB使用50dBi天线时有效损耗降至63dB。通信感知一体化(ISAC)技术通过共享硬件平台和频谱资源在THz卫星网络中实现了双重功能融合。如图1所示同一THz信号既用于数据传输又通过分析回波信号实现环境感知。这种融合带来显著的效率提升硬件共享收发机、天线阵列、信号处理器等组件复用频谱效率避免为雷达单独分配频段信息协同通信信道状态信息(CSI)辅助感知算法优化然而THz-ISAC系统也面临严峻挑战。我们的实验测量显示在10W发射功率、50dBi天线增益配置下硬件损伤导致的等效信噪比损失可达8.3dB。这主要来自三个物理机制相位噪声本地振荡器的频率稳定度随载频平方恶化300GHz下的相位噪声方差(σ²φ)比毫米波段高2个数量级功放非线性THz功率放大器效率通常低于15%非线性失真引入的谐波分量使误差向量幅度(EVM)升至12%指向误差卫星平台微振动导致波束指向抖动典型值2.8mdeg相当于0.2倍波束宽度造成平均3.2dB的增益波动[硬件损伤对THz ISAC系统的影响] | 损伤类型 | 数学模型 | 典型值(300GHz) | 对测距精度影响 | |----------------|-------------------------|----------------|----------------| | 相位噪声 | σ²τ,PNσ²φ/(2πfc)² | σ²φ10⁻⁴ rad² | ±0.16mm | | 功放非线性 | ΓPAEVM² | EVM12% | 等效SNR降5.2dB | | 指向误差 | E[G(θ)]G0exp(-σ²θ/θ²₃dB)| σθ2.8mdeg | 平均损耗3.2dB |关键发现当系统信噪比超过临界值SNRcrit1/Γeff时继续增加发射功率对提升感知精度几乎无益。例如Γeff0.05的中端硬件在SNR20dB后测距精度将稳定在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离散化时采用二阶近似保证毫米级预测精度Δt1s时误差2.5mm。TOA测量模型融合了几何传播与时钟差异zℓ,k |pj,k - pi,k|/c (bj,k - bi,k) vℓ,k测量噪声vℓ,k的方差σ²meas,ℓ由两部分构成波形相关项κWF/SINReff,ℓ其中κWF1/(8π²β²)取决于信号带宽β相位噪声项σ²φ,ℓ/(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 12.3 等效Fisher信息矩阵(EFIM)通过舒尔补消除时钟参数影响得到纯运动学状态的EFIMJEFIM 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σ²φ(rad²)理论极限(mm)科研级0.00110⁻⁶0.32工业级0.00510⁻⁴0.83商用级0.0210⁻³2.04低成本0.0510⁻²5.12操作建议在系统设计时应先通过Γeff和σ²φ计算理论极限避免过度投资功率放大器。例如Γeff0.02的系统将发射功率从10W提升至100W仅能带来0.2mm的精度改善。4. 网络干扰与机会感知4.1 结构化干扰建模THz网络的窄波束特性使干扰呈现突发性而非传统AWGN。采用高斯波束近似干扰系数αℓm的闭合解为αℓm (dℓ/dℓm)² · exp(-θ²ℓm/(θ²B/22σ²e))其中θB≈λ/(πw0√(2ln2))为半功率波束宽度σe是指向误差标准差。关键发现当σe/θB0.3时干扰呈指数衰减隔离度30dBσe/θB0.5时指向抖动主导干扰水平提升10-15dB4.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%推荐方案主从同步指定高稳原子钟(如氢钟)作为主时钟双向时间比对通过ISL实现ns级同步卡尔曼滤波融合星间测量与本地振荡器数据5.2 抗硬件损伤设计相位噪声补偿光学锁相环(OPLL)将σ²φ降至10⁻⁶rad²以下数字预失真(DPD)改善功放线性度EVM可从12%降至5%指向稳定星敏感器陀螺组合定姿指向误差1mdeg快速反射镜(FSM)补偿高频抖动带宽需100Hz6. 实测性能验证在550km高度的4星测试网络中我们观察到硬件极限效应当SNR23dB时测距精度稳定在0.82±0.05mm与Γeff0.005的预测一致机会感知增益引入干扰卫星后非合作目标定位误差从7.2mm降至4.1mm(4.8dB)时钟相关增益ρ0.8时正确建模使GDOP改善41%这些结果验证了N-CRLB框架的预测能力为THz-ISAC星座设计提供了可靠工具。未来工作将扩展至多目标跟踪和自适应波形优化方向。