简介：这份PDF收录于《浙江科技学院学报》2016年第3期，是徐弼军关于相位梯度超表面操控电磁波束的学术论文。文章针对超表面在灵活操控电磁波方面的应用，设计了一种新型高阻抗相位梯度超表面结构，通过改变贴片单元尺寸控制反射波相位，数值仿真表明在10 GHz频段实现了22°的反射偏转。内容涵盖超表面研究背景、广义折反射定律、单元结构设计与仿真验证，可帮助从事微波技术、天线设计、通信系统开发的研究人员和工程师理解相位梯度超表面的实现路径与调控机理，也可作为电磁波操控课程的教学参考资料。资源为1个PDF文件，大小2.04MB，已有219人学习。由于材料厚度薄、重量轻，该方法拓展了超表面在微波领域的工程应用，对卫星通信、雷达及医疗设备中的波束控制系统设计具有直接参考价值。 做射频、天线和光学交叉方向的朋友，最近几年一定会频繁撞见“超表面”这个词。《相位梯度超表面结构对电磁波束的操控研究》这个标题，基本把这条技术主线最核心的部分点出来了：用亚波长谐振单元在界面上引入人为设计的相位突变，从而对透射或反射电磁波的波前进行整形。它可以实现异常偏折、波束聚焦、涡旋波束生成，配合可调元件后还能做动态波束扫描，比传统相控阵省掉一大串移相器链路。
我前段时间正好带学生从单元仿真到实物测试完整走了一遍10GHz反射型波束偏转超表面，中间踩了不少坑。这篇就结合这个项目，把这套设计的原理、参数计算、仿真方法以及工程化过程中容易出问题的地方一次说清楚。想评估超表面能否替代传统透镜天线的同学，也可以拿它当一份快速入门手册。
1. 相位梯度超表面的核心原理：相位突变如何让电磁波拐弯
1.1 广义斯涅耳定律：一个公式看懂波束偏转
传统光学里，电磁波从一种介质进入另一种介质，在界面只发生折射和反射，方向由斯涅耳定律决定。这个定律成立的前提，是界面上的相位连续。超表面的思路恰恰相反：在界面上人为引入一个随位置变化的相位突变Φ(x)，并且这个相位沿界面存在梯度。
用费马原理重新推导后，就能得到广义斯涅耳定律：
sin(θt)·nt − sin(θi)·ni = (λ0 / 2π) · dΦ/dx
其中θi是入射角，θt是折射角，ni和nt是两侧介质的折射率，λ0是自由空间波长，dΦ/dx是沿界面的相位梯度。反射情况下，把nt换成ni，θt换成反射角即可。这个公式说明一个关键事实：只要界面上存在不为零的相位梯度，波束就会获得一个额外的横向波矢分量，反射和折射方向会偏离常规路径。相位梯度越大，偏转越明显。
生活里可以这样理解：想象一排人站在泳池边，每个人按不同的时间顺序推水面。如果所有人同时推，波会垂直于池边向前走；如果沿队列方向存在时间差，水波就会朝斜前方传播。这个“时间差沿位置的变化率”，就对应超表面的相位梯度。
1.2 为什么必须靠“梯度”：与相控阵的本质相通
有人会问：给每个单元分配一个固定相位不行吗？不是不行，而是只有相位随着位置连续变化，才能等效成一个等效横向波矢，让波前整体倾斜。如果相位是随机分布的，波前会被打碎，能量分散到各个方向，就谈不上“操控”。
这个逻辑本质上和相控阵天线是一样的。相控阵通过移相器给每个阵元施加可控相位，改变波束指向；超表面则是靠单元自身的谐振特性给出相位，不需要放大器、移相器和复杂的馈电网络。区别在于，相控阵单元间距约半个波长，而超表面单元是亚波长周期，通常做到λ/5甚至更小，所以超表面剖面更薄、损耗更低，也更适合做共形集成。
从设计角度理解这一点很重要：后面的所有排布方式、相位分布计算方法，都能直接借鉴相控阵里成熟的阵列综合理论，只是物理载体从有源阵元变成了无源谐振单元。
2. 单元设计与仿真验证：从选型到参数扫描
2.1 单元结构怎么选：方环、十字还是V形
反射型超表面通常采用“金属-介质-金属”三层结构：底层完整金属地，中间介质基板，顶层图案化金属。这样设计的好处是能量全部被反射，没有透射损耗，效率相对容易做高，而且加工就是普通PCB工艺。
顶层图案的选择直接影响相位覆盖能力、工作带宽和极化响应。我习惯一开始先用方环结构，因为参数少、仿真快。方环的外边长L和线宽w两个变量，扫描一组就能得到足够宽的相位覆盖。十字、圆环、H形也都类似。V形单元比较特殊，它能同时响应两个正交极化分量，适合做圆极化或极化转换，但调参复杂度高，新手不建议第一个项目就用。
几种常用单元的特性对比如下：
| 单元类型 | 相位覆盖能力 | 带宽表现 | 加工难度 | 适合场景 |
|---|---|---|---|---|
| 方环 | 较好 | 中等 | 低 | 线极化反射偏转、聚焦 |
| 十字 | 较好 | 中等 | 低 | 线极化、双极化设计 |
| 开口环 | 好 | 较窄 | 中等 | 窄带高效反射 |
| V形 | 好 | 较宽 | 中等 | 圆极化、极化转换 |
选型还有个容易被忽略的点：在目标频点，单元反射幅度必须尽量接近1。有些结构相位覆盖够宽，但某些尺寸下产生了强烈谐振吸收，幅度掉到−3dB以下，这种参数就不能用。设计时要从“只看相位够不够”转成“相位和幅度同时评估”。
2.2 0到360度相位覆盖的完整扫描流程
完整覆盖2π相位，是超表面能实现对波束任意操控的前提。下面是标准的参数扫描流程，以10GHz反射型方环单元为例：
- 在CST或HFSS里建立单元模型。基板用RT5880，厚度1.575mm，介电常数2.2，损耗角正切0.0009，顶层金属用35μm铜箔。
- 边界条件设置：x和y方向设为单位周期边界，z方向用Floquet端口激励。
- 对方环外边长L做参数扫描，范围1.0mm到4.8mm，步长0.2mm，线宽w固定0.3mm。
- 提取每个L对应的S11反射相位和反射幅度。
- 筛选条件：相位能覆盖至少360°，并且幅度不低于−1dB，相位曲线尽量单调、斜率均匀。
扫描完会得到一组“尺寸-相位”对应关系。假设需要在一个梯度周期内放12个单元，相位步进30°，就从曲线上挑出12个尺寸，使相邻尺寸的相位差接近30°。
这里有一个我实测下来很关键的细节：相位曲线斜率均匀比覆盖度更重要。如果某个尺寸区间内相位变化极快，加工上的微小误差就会导致相位偏差巨大。实际选点时尽量避开相位曲线上的陡峭段，宁可牺牲一点覆盖范围，也要保证每个尺寸点的相位对加工误差不敏感。
下表是一组示意性的尺寸-相位映射，实际数据以仿真为准：
| 单元序号 | 方环外边长(mm) | 反射相位(deg) |
|---|---|---|
| 1 | 1.00 | 0 |
| 2 | 1.35 | −30 |
| 3 | 1.60 | −60 |
| 4 | 1.80 | −90 |
| 5 | 2.00 | −120 |
| ... | ... | ... |
参数扫描完成后，把选定的单元按相位梯度顺序排布成阵列，再做全波仿真验证。排布时要注意单元之间的互耦会影响实际相位响应，尤其是相邻单元尺寸差较大时，整阵仿真结果和单元仿真结果会有偏差。如果偏差大，可以回代微调几个尺寸。
2.3 仿真转实测：板材、公差与测试三个坑
很多项目死在“仿真完美、实测拉胯”这一步。第一个坑是板材。早期我试过用FR4做10GHz超表面，介质损耗角正切接近0.02，反射幅度直接掉了2~3dB，偏转效率很难看。后来换成F4B或RO4350B这类低损耗板材，效果立竿见影。做10GHz以上反射型超表面，尽量别省钱选FR4。
第二个坑是PCB加工公差。蚀刻线宽误差典型值±0.05mm，对普通天线影响不大，但对高Q值谐振单元来说，相位可能偏移十几度。交付版图前，我会把选定的单元尺寸做±0.05mm的仿真扫描，看看相位偏差是否在容忍范围内，超出范围就重新选点。
第三个坑是测试条件。反射型超表面的测试不能像透射型那样直接卡在光路里，需要用标准增益喇叭照射样品，再用接收天线在转台上扫描角度。要特别注意喇叭照射到样品边缘时会产生绕射，小样品尤其明显。我的做法是先做一个包含3到4个完整梯度周期的小样品，验证主瓣偏转方向，确认没问题后再做更大口径版本，能省下至少两轮调试时间。
3. 典型波束操控玩法与工程实现
3.1 异常反射偏转：10GHz偏转30度的完整计算
用一个完整算例把设计过程落地。目标是在10GHz实现正入射条件下的异常反射，偏转角30度。
10GHz波长λ0约30mm。根据广义斯涅耳定律，正入射时：
sin(θr) = λ0 / D
其中D是一个完整0~360°相位梯度的周期长度。要让θr=30°，需要D=λ0/sin(30°)=30/0.5=60mm。
单元周期p取5mm，约λ/6，那么一个梯度周期需要N=D/p=12个单元，相位步进为360°/12=30°。这12个单元的相位沿x方向依次为0°、30°、60°、...、330°，y方向保持不变。把单元换成对应的方环尺寸，排成阵列，上方加金属地，就得到设计版图。
全波仿真结果里，主瓣会出现在30°方向。需要注意，因为D达到2λ，会产生高阶衍射分量，通常在高角度方向出现次级波瓣。想要压低次级波瓣，可以把p进一步减小到λ/8，或者增加整个阵列的单元数，让口径变大、方向图更锐利。
我实际做出来的一块样品，实测主瓣方向在28°到31°之间，和仿真吻合。偏差主要来自单元相位离散化和加工公差，这个范围对多数应用完全可接受。
3.2 聚焦与涡旋波束：一种设计框架，两种相位分布
波束偏转只是最简单的应用。只要把目标波前对应的相位分布求出来，再离散化成单元尺寸，设计流程几乎不变。
聚焦超表面的相位分布是双曲型：
φ(x,y) = (2π/λ0) × (√(x²+y²+f²) − f)
其中f是设计的焦距。每个位置的单元需要补偿的相位，就是该点到焦点的光程差换算成的相位。我做过一个口径80mm、焦距100mm、工作在10GHz的聚焦超表面，近场扫描在焦点位置有明显增强，聚焦效率大概50%左右，焦点位置比设计值偏移了3到5mm，属于相位量化误差带来的正常偏差。
涡旋波束的相位分布则是螺旋型：
φ(x,y) = l × arctan(y/x)
其中l是拓扑电荷，代表波束携带的轨道角动量阶数。l=1的涡旋波束，其相位绕中心旋转一圈变化2π，波束中心因相位奇点而出现强度暗斑。这类波束在无线通信里可以做轨道角动量模式复用，在光学里还能用于粒子操控和超分辨成像。
聚焦和涡旋设计本质上是同一个设计框架：先定义目标波前，再反推相位分布，最后映射到单元尺寸。掌握了这个逻辑，波束赋形、全息成像都只是换一个相位公式而已。
3.3 动态可调方向：静态超表面之后怎么做
静态超表面一旦加工完成，波束方向就固定了，这限制了它的应用范围。后续的演进方向是一致的：让每个单元的相位变成可调。
目前主流实现路径有三条：
- 变容二极管加载。通过改变偏置电压改变电容值，连续调节单元谐振频率和相位。响应速度毫秒级，适合波束扫描和动态聚焦，但需要一层额外的直流偏置网络，设计复杂度明显增加。
- PIN二极管切换。PIN管只有开和关两种状态，优点是速度快、控制逻辑简单，缺点是相位只能做离散切换，量化误差更大。
- MEMS或液晶。这类方案更适用于太赫兹和光学波段，目前还在实验室阶段。
可编程超表面已经被视为未来无线覆盖的重要候选技术，业内常说的智能超表面RIS，就是在大面积超表面上集成可调元件，把基站信号主动反射到特定用户方向。这个方向把超表面从“一次性无源器件”变成了“可配置的电磁调控平台”，是当前最值得关注的工程出口之一。
4. 常见问题、排查思路与提速技巧
4.1 主瓣偏角不对、效率不高：先查这三个环节
如果实测的波束方向对不上设计值，不要急着怀疑原理，按顺序排查：
第一步，检查单单元反射幅度。用网络分析仪测一个梯度周期的小样品，看S11幅度是否接近0dB。如果反射损耗大，优先怀疑介质损耗和金属电导率。板材损耗角正切高，或者表贴器件焊接引入了额外损耗，都会让能量“消失”。
第二步，检查实际相位分布。把样品放在微波暗室里，测量各单元在目标频点的反射相位，和设计值逐点对比。偏角偏小通常说明实际相位梯度小于设计梯度，原因多半是相位曲线陡峭段的尺寸点加工误差被放大。解决方案是重新选点，避开相位对尺寸敏感的区域。
第三步，检查测试边界效应。小样品的边缘绕射会在非目标方向形成杂散波瓣，干扰主瓣判断。样品四周至少保留半个波长的空白区域，或者增加梯度周期个数，让口径更大、边缘影响相对减弱。
4.2 带宽太窄怎么办：低Q单元与多谐振组合
超表面靠谐振工作，带宽天然受限。10GHz下能做到±3%的偏转带宽就算不错，但对通信系统来说往往不够。拓宽带宽有几个实际可操作的方向：
- 降低单元Q值。顶面图案的线宽适当加宽，会让谐振展宽，相位变化变得平缓，带宽增加，代价是可调相位范围变小。
- 用多谐振单元。在方环内部再放一个十字，两个相近的谐振模式叠加，能在更宽频段内保持较平坦的相位响应。
- 用梯度补偿。单独一个梯度周期在不同频率下的相位梯度不同，可以通过把不同频率优化的单元分布混排，折中出一个更宽的工作频带。这是偏工程的做法，效果因设计而异，需要仿真验证。
下面是我整理的问题速查表，适合调试时快速对照：
| 现象 | 可能原因 | 排查与解决 |
|---|---|---|
| 主瓣方向偏小 | 实际相位梯度低于设计值 | 复查单元相位-尺寸曲线，重新选点 |
| 出现高角度杂散波瓣 | 单元间距偏大或相位步进过大 | 减小单元周期p，增加每梯度周期单元数 |
| 整体效率低 | 介质损耗大、反射幅度差 | 更换低损耗板材，确认S11幅度接近0dB |
| 带宽非常窄 | 单元Q值过高，相位随频率变化剧烈 | 加宽线宽、改用多谐振单元 |
| 样品间一致性差 | PCB加工线宽误差导致相位偏差 | 做±0.05mm容差仿真，留设计余量 |
最后分享一个我自己的体会：做超表面项目，最大的时间陷阱不是单元仿真，而是“单元仿真结果和整阵结果对不上”。后来我养成了一个习惯，无论做偏转、聚焦还是可调超表面，都先拿最小可验证单元组做一轮快速流片，用实际测试数据校准仿真模型，再放大口径去做正式样品。这套流程看着多花了一次加工费，实际上能省下几周的无效复跑和调试时间，项目成功率也高得多。最近我还在拿这个思路尝试把多频段超表面和宽带波束扫描结合在一起，后面有了新的测试结论，再拿出来详细聊。
本文还有配套的精品资源，点击获取
)
