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介质基板底层金属接地板关键设计参数包括单元周期P20mm贴片尺寸L16.62mm通过正交极化方向的金属化过孔设计使系统在±40°入射角范围内保持稳定的反射特性。实测显示在5GHz频点超表面能提供0°、90°、180°、270°四种离散相位响应反射损耗均小于3dB。注意事项在实际组装时二极管偏置线路的走线需要做λ/4阻抗变换处理避免射频信号泄漏到控制电路。我们采用蛇形走线设计既满足直流馈电需求又保持了单元间的电磁隔离。2.2 等效反射系数调制技术传统2比特超表面仅能提供四种固定相位状态远不能满足射频计算对连续反射系数的需求。我们创新性地提出空间编码调制法通过快速切换相邻单元的编码状态1kHz刷新率在时域上合成等效的连续反射系数Γ_equ。如图3d所示这种方法能将可用反射系数扩展到整个IQ平面内半径为A_ref的圆域。具体实现算法如下将目标复数反射系数Γ_target分解为四个基本状态的线性组合根据分解权重动态分配相邻单元的编码模式通过空间平均效应实现等效调制实测表明该方法可将相位分辨率提升至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展示了两种典型卷积核的实现效果汉宁窗sinc函数低通滤波器截止频率25kHz带外抑制30dB13位巴克码匹配滤波器主旁瓣比达到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.8ms4.2 距离测量模式采用脉冲压缩技术时距离分辨率取决于信号带宽ΔR c/(2B)我们测试了两种波形LFM脉冲500kHz带宽实测距离分辨率310m与理论值300m吻合Gold码序列在12km距离上测距误差50m验证了方法的普适性特别值得注意的是系统支持波形参数的软件定义重构。如图6f所示只需通过FPGA重新加载编码序列即可在20μs内完成从LFM到Gold码的切换无需任何硬件改动。5. 工程实践中的关键挑战5.1 单元间耦合抑制密集阵列中单元互耦会导致反射系数失真。我们通过三项措施解决采用六边形栅格排列比矩形栅格降低15%的耦合强度在单元间添加电磁带隙EBG结构表面波抑制比提升8dB引入预失真补偿算法将互耦引起的相位误差从±12°降至±2°5.2 动态范围优化系统动态范围受两个因素制约接收机灵敏度-110dBmSTCM非线性失真三阶交调点25dBm通过自适应功率控制算法在保持最佳信噪比的同时将带内失真分量抑制在-50dBc以下。实测显示当输入信号功率在-30dBm至-10dBm范围内变化时傅里叶变换的频谱纯度保持稳定。6. 性能对比与优势分析与传统DSP方案相比STCM射频计算系统展现出显著优势指标传统方案STCM方案提升幅度处理延迟2.1ms0.3ms7倍功耗(5GHz CW处理)3.8W0.9W76%硬件复杂度需要ADC/FPGA纯无源结构简化60%可重构性固件升级实时编程响应快100倍这种架构特别适合机载雷达等对SWaPSize, Weight and Power要求严苛的场景。在无人机载合成孔径雷达SAR测试中系统将信号处理模块重量从1.2kg降至0.3kg同时功耗降低4倍。7. 未来演进方向当前系统还存在两个主要局限工作带宽受限4.8-5.2GHz下一步将开发倍频程超宽带超表面计算功能较为基础正在研究神经网络赋能的智能超表面计算架构我们已经验证了通过级联多个STCM实现连续傅里叶变换-滤波-逆变换的处理链条可行性。这为构建全模拟信号处理系统奠定了基础有望将雷达信号处理能效比再提升一个数量级。