时空编码超表面在射频计算中的创新应用

1. 项目概述：基于时空编码超表面的射频计算革命

在传统电子系统中，射频信号处理长期受限于数字信号处理器（DSP）的架构瓶颈。每次信号转换都需要经过模数/数模转换（ADC/DAC）模块，这不仅导致硬件冗余和能耗增加，还引入了不可避免的传输延迟。针对这一行业痛点，我们团队开发了基于时空编码超表面（STCM）的可编程射频计算系统，实现了信号处理从"电路域"到"电磁波域"的范式转移。

这个系统的核心突破在于：通过超表面原子单元的时空编码策略，直接在电磁波反射过程中完成傅里叶变换、卷积运算等关键信号处理操作。就像用光学棱镜分解白光一样，我们的STCM能够实时"分解"和"重组"射频信号的频谱特征，而无需传统数字处理链条的介入。在5GHz工作频率下的实测表明，系统对线性调频信号（LFM）的频谱分析误差小于1%，速度测量精度达到±0.75m/s，距离分辨率优于300米。

2. 系统架构与工作原理

2.1 时空编码超表面的物理实现

我们采用的STCM样品为8×12阵列的2比特可重构超表面，每个单元采用三层结构设计（图3a）：

• 顶层：六边形贴片与直流馈线组成的辐射层，集成两个PIN二极管实现00/01/10/11四种状态切换
• 中层：介电常数2.65、损耗角正切0.001的F4B介质基板
• 底层：金属接地板

关键设计参数包括单元周期P=20mm，贴片尺寸L1=6.62mm，通过正交极化方向的金属化过孔设计，使系统在±40°入射角范围内保持稳定的反射特性。实测显示在5GHz频点，超表面能提供0°、90°、180°、270°四种离散相位响应，反射损耗均小于3dB。

注意事项：在实际组装时，二极管偏置线路的走线需要做λ/4阻抗变换处理，避免射频信号泄漏到控制电路。我们采用蛇形走线设计，既满足直流馈电需求，又保持了单元间的电磁隔离。

2.2 等效反射系数调制技术

传统2比特超表面仅能提供四种固定相位状态，远不能满足射频计算对连续反射系数的需求。我们创新性地提出空间编码调制法，通过快速切换相邻单元的编码状态（>1kHz刷新率），在时域上合成等效的连续反射系数Γ_equ。如图3d所示，这种方法能将可用反射系数扩展到整个IQ平面内半径为A_ref的圆域。

具体实现算法如下：

1. 将目标复数反射系数Γ_target分解为四个基本状态的线性组合
2. 根据分解权重动态分配相邻单元的编码模式
3. 通过空间平均效应实现等效调制

实测表明，该方法可将相位分辨率提升至1°，幅度调节步进达到0.1dB，完全满足傅里叶变换和卷积运算的精度要求。

3. 核心算法实现

3.1 电磁波域傅里叶变换

傅里叶变换的物理实现本质上是将数学基底函数映射到超表面的时变反射系数上。当入射信号s(t)与STCM相互作用时，反射信号可表示为：

r(t) = Γ(t)·s(t) // 时域相乘

通过设计Γ(t)为一系列正交的复指数函数，即可在反射波中合成输入信号的频谱。以32点FFT为例，STCM的反射系数矩阵构建为：

Γ_FFT = [1, 1, 1, ..., 1 1, e^(-j2π/N), e^(-j4π/N), ..., e^(-j2π(N-1)/N) ... 1, e^(-j2π(N-1)/N), ..., e^(-j2π(N-1)^2/N)]

图2的实验结果验证了该方法的有效性：对5kHz单音信号进行32点和64点傅里叶变换时，频谱峰值位置误差小于0.1%，旁瓣抑制比达到25dB以上。值得注意的是，通过简单地延长采样窗口（1.6ms→3.2ms），系统频率分辨率可从625Hz提升至312.5Hz，展现出优异的可重构性。

3.2 实时卷积运算

卷积核的物理实现依赖于STCM对时反信号的精确复现。对于匹配滤波应用，我们将期望信号s(t)的共轭时反s*(-τ)编码为超表面的反射系数序列。当回波信号y(t)入射时，反射过程即完成卷积积分：

y(t) ⊗ s(t) = ∫y(τ)·s*(τ-t)dτ

图5展示了两种典型卷积核的实现效果：

1. 汉宁窗sinc函数低通滤波器：截止频率25kHz，带外抑制>30dB
2. 13位巴克码匹配滤波器：主旁瓣比达到22.3dB，接近理论极限值22.8dB

实测中发现，卷积运算对时序同步精度要求极高。我们采用FPGA产生纳秒级精度的控制时序，确保编码序列与输入信号的严格对齐，将时间抖动控制在100ps以内。

4. 雷达应用验证

4.1 速度测量模式

在连续波雷达架构下，多普勒频移与目标速度满足：

f_d = 2v/λ

STCM通过实时傅里叶变换直接提取f_d。如图6所示，当采样点数从500增至1000时：

• 速度分辨率从±1.5m/s提升至±0.75m/s
• 对200m/s目标的测量误差从1m/s降至0.5m/s
• 处理延迟从5ms降低到2.8ms

4.2 距离测量模式

采用脉冲压缩技术时，距离分辨率取决于信号带宽：

ΔR = c/(2B)

我们测试了两种波形：

1. LFM脉冲（500kHz带宽）：实测距离分辨率310m，与理论值300m吻合
2. Gold码序列：在12km距离上测距误差<50m，验证了方法的普适性

特别值得注意的是，系统支持波形参数的软件定义重构。如图6f所示，只需通过FPGA重新加载编码序列，即可在20μs内完成从LFM到Gold码的切换，无需任何硬件改动。

5. 工程实践中的关键挑战

5.1 单元间耦合抑制

密集阵列中单元互耦会导致反射系数失真。我们通过三项措施解决：

1. 采用六边形栅格排列，比矩形栅格降低15%的耦合强度
2. 在单元间添加电磁带隙（EBG）结构，表面波抑制比提升8dB
3. 引入预失真补偿算法，将互耦引起的相位误差从±12°降至±2°

5.2 动态范围优化

系统动态范围受两个因素制约：

• 接收机灵敏度：-110dBm
• STCM非线性失真：三阶交调点+25dBm

通过自适应功率控制算法，在保持最佳信噪比的同时，将带内失真分量抑制在-50dBc以下。实测显示，当输入信号功率在-30dBm至-10dBm范围内变化时，傅里叶变换的频谱纯度保持稳定。

6. 性能对比与优势分析

与传统DSP方案相比，STCM射频计算系统展现出显著优势：

这种架构特别适合机载雷达等对SWaP（Size, Weight and Power）要求严苛的场景。在无人机载合成孔径雷达（SAR）测试中，系统将信号处理模块重量从1.2kg降至0.3kg，同时功耗降低4倍。

7. 未来演进方向

当前系统还存在两个主要局限：

1. 工作带宽受限（4.8-5.2GHz），下一步将开发倍频程超宽带超表面
2. 计算功能较为基础，正在研究神经网络赋能的智能超表面计算架构

我们已经验证了通过级联多个STCM实现连续傅里叶变换-滤波-逆变换的处理链条可行性。这为构建全模拟信号处理系统奠定了基础，有望将雷达信号处理能效比再提升一个数量级。