双层频率多元反射超表面及设计方法

1.本发明涉及反射超表面技术领域，主要涉及一种双层频率多元反射超表面及设计方法。


背景技术：

2.随着通信系统的快速发展，多功能微波器件被广泛地应用在信号传输、成像系统等集成设备中，但传统多功能微波器件因其存在体积大、损耗高、效率低等问题，并不符合微波器件的鲁棒性和实用性。相比较而言，具有亚波长厚度的超表面具有对电磁波的幅度、相位和极化良好的调控能力。对于传统光学透镜依靠距离传输获得相位积累，超表面通过与入射电磁波的共振耦合即可获得突变相位，展现出超表面对电磁波强大的操控能力。因为超表面优越的电磁调控能力，使其在电磁波异常反射/折射、雷达散射截面减缩(rcs)、全息、聚焦和涡旋光束发生器等方面都取得了重要进展。
3.为提高设备的集成度和紧凑性，研究人员广泛开展关于多功能超表面器件的研究。目前主要通过电磁波对超表面激励信息(如频率、极化、方向和位置)来实现多功能。其中，频率是电磁波携带的重要信息，频率复用使得器件具有高效的频谱利用率，使频率多功能超表面在电磁调控方面得到广泛应用。其中设计频率多功能超表面的核心在于依靠超表面上工作于不同频段的非对称各向异性单元对电磁波幅度、极化、相位的调控。因其单元具有对电磁波良好的调控能力使得频率多元超表面被广泛应用，其中频率多元相位调控超表面因具有良好的传输特性和全空间电磁波调控能力被广泛应用于多功能微波集成器件中，如涡旋光束和全息图的多频多通道超表面，消色差的全息超表面、宽带自旋解耦的无畸变双相超表面。虽然频率多功能超表面能在多个频带实现功能的集成以应对不同的工作环境，但大多数频率多功能超表面单元都是依靠多个金属谐振器在单层介质上拼接或在多层介质上利用空间复用(透、反射集成)等方法来实现多功能。这种方法由于通道与通道之间不可避免会存在串扰，使效率降低。为避免串扰，这就要求频率多功能超表面的每个工作频段不能相距太近。


技术实现要素：

4.本发明提出一种双层频率多元反射超表面，其特征在于，所述双层频率多元反射超表面通过高效互补谐振器单元的周期性排列构成，所述高效互补谐振器单元包括两层介质板、在所述两层介质板上构建互补形式的双c形开槽谐振器和双c形金属谐振器以及金属地板；
5.所述双层频率多元反射超表面在低频比的条件下实现模式数l＝3的聚焦oam波束和零阶贝塞尔波束。
6.更进一步地，所述高效互补谐振器单元包括依结构顺序设置的双c形开槽谐振器、第一层介质板、双c形金属谐振器、第二层介质板和金属地板。
7.更进一步地，所述第一层介质板厚度h1＝1.5mm，第二层介质板厚度h2＝1.5mm，介
质板材料为f4b；
8.所述双c形开槽谐振器、双c形金属谐振器和金属地板采用的材料为铜，其电导率为σ＝5.8
×
107s/m。
9.更进一步地，所述高效互补谐振器单元周期p＝10.2mm；
10.所述双c形开槽谐振器结构参数为：外环外径r1＝4.75mm、外环内径r2＝4.35mm、中间槽宽为w1＝0.3mm、内金属环宽w2＝0.4mm，外环与内环连接部分金属宽度g1＝0.9mm；
11.所述双c形金属谐振器结构参数为：金属环内径宽度r3＝3.0mm、金属环宽度w3＝0.8mm，金属环缝隙宽度g2＝0.3mm。
12.更进一步地，所述高效互补谐振器单元在两个谐振频率f1和f2处均实现100％交叉圆极化波转化的独立几何相位调控，其中，双c形开槽谐振器工作于频率f1，负责该频率处的几何相位功能调控，双c形金属谐振器工作于频率f2，负责该频率处的几何相位功能调控，即在谐振频率f1和f2处线极化波入射条件下所述高效互补谐振器单元满足：
13.|r
xx
|＝|r
yy
|＝1
[0014][0015]
其中，r
xx
为在x极化的反射系数，r
yy
为在y极化的反射系数，为在y极化的反射系数，为线极化波同极化反射相位。
[0016]
更进一步地，所述双层频率多元反射超表面在工作频率f1＝9.2ghz处实现模式数l＝3的聚焦oam波束，利用相位叠加原理将其分解为聚焦相位部分与涡旋相位部分，对于聚焦相位而言，所述高效互补谐振器单元满足如下相位分布：
[0017][0018]
其中，f0＝214mm为其焦距，p为单元周期，λ为入射电磁波的波长，是任意相位常量，m为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目；
[0019]
对于涡旋波束的相位分布，单元对入射波进行的相位补偿需要满足如下公式：
[0020][0021]
其中，l＝3为oam波束的模式数；
[0022]
依据电磁波的相位叠加原理，将涡旋相位与聚焦相位进行相位叠加，对于双层频率多元反射超表面的最终单元相位分布满足如下公式：
[0023][0024]
表示单元在谐振频率f1＝9.2ghz处最终的反射相位。
[0025]
更进一步地，所述双层频率多元反射超表面在工作频率f2＝11.2ghz处实现零阶贝塞尔波束，所述高效互补谐振器单元在谐振频率f2的补偿相位满足如下公式：
[0026][0027]
其中，p为单元的周期，λ为入射电磁波的波长，β＝30
°
为贝塞尔波束的衍射半角，m为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目。
[0028]
还提出了一种双层频率多元反射超表面设计方法，其特征在于，所述双层频率多元反射超表面设计方法包括以下步骤：
[0029]
步骤1，设计双层频率多元反射超表面所需的高效互补谐振器单元，使高效互补谐振器单元在两个谐振频率处实现100％交叉圆极化波的相位调控；
[0030]
步骤2，通过高效互补谐振器单元非周期性排列构成了频率多元反射超表面；
[0031]
步骤3，对双层频率多元反射超表面性能分析，确认双层频率多元反射超表面实现预期功能；
[0032]
在步骤2中，双层频率多元反射超表面最终构建了在以f1＝9.2ghz与f2＝11.2ghz为工作频率的，分别实现了模式数l＝3的聚焦oam波束和零阶贝塞尔波束的频率复用多功能超表面。
[0033]
更进一步地，在步骤2中，所述双层频率多元反射超表面在工作频率f1＝9.2ghz处实现模式数l＝3的聚焦oam波束，利用相位叠加原理将其分解为聚焦相位部分与涡旋相位部分，对于聚焦相位而言，所述高效互补谐振器单元满足如下相位分布：
[0034][0035]
其中，f0＝214mm为其焦距，p为单元周期，λ为入射电磁波的波长，是任意相位常量，m为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目；
[0036]
对于涡旋波束的相位分布，单元对入射波进行的相位补偿需要满足如下公式：
[0037][0038]
其中，l＝3为oam波束的模式数；
[0039]
依据电磁波的相位叠加原理，将涡旋相位与聚焦相位进行相位叠加，对于双层频率多元反射超表面的最终单元相位分布满足如下公式：
[0040][0041]
表示单元在谐振频率f1＝9.2ghz处最终的反射相位。
[0042]
更进一步地，在步骤2中，所述双层频率多元反射超表面在工作频率f2＝11.2ghz处实现零阶贝塞尔波束，所述高效互补谐振器单元在谐振频率f2的补偿相位满足如下公式：
[0043][0044]
其中，p为单元的周期，λ为入射电磁波的波长，β＝30
°
为贝塞尔波束的衍射半角，m为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目。
[0045]
本发明的优点在于：
[0046]
本发明提供的频率多元反射超表面相比于传统的频率多功能超表面，该超表面在两个谐振频率处具有良好的独立性，且严格满足几何相位原理。
[0047]
传统的频率多功能超表面，其频段所占频谱范围宽，相对频谱利用率低，结构复杂不易加工。本发明提供的频率多元反射超表面所占频谱宽度窄、频比低、结构简单和易于加
工。
附图说明
[0048]
图1为双层频率多元反射超表面的结构图；
[0049]
图2为双层频率多元反射超表面的单元布局和结构参数图；
[0050]
图3为f1与f2谐振频率处x、y极化波入射下互补谐振器单元的电磁响应；
[0051]
图4为f1与f2谐振频率处双c形开槽谐振器、双c形金属谐振器和金属地板的表面电流分布正视图与侧视图；
[0052]
图5为不同α1(α2)的双c形开槽谐振器(双c形金属谐振器)随α2(α1)变化的反射幅度与相位；
[0053]
图6为不同α1(α2)的双c形开槽谐振器(双c形金属谐振器)随频率变化的反射幅度与相位；
[0054]
图7为模式数l＝3的聚焦oam波束单元相位分布图；
[0055]
图8为零阶贝塞尔波束单元相位分布图；
[0056]
图9为零阶贝塞尔波束成波原理图；
[0057]
图10模式数l＝3的聚焦oam波束在焦平面f0＝214mm处振幅、相位仿真结果图；
[0058]
图11零阶贝塞尔波束不同位置|e
lcp
|^2归一化分布图；
[0059]
图12为贝塞尔波束在yoz面不同位置的|e
lcp
|^2归一化分布曲线图。
具体实施方式
[0060]
下面结合附图对本发明的技术方案进行更详细的说明，本发明包括但不仅限于下述实施例。
[0061]
本发明提供了一种双层频率多元反射超表面，其超表面主要通过高效互补谐振器单元的周期性排列构成。高效互补谐振器单元由在两层介质板上构建互补形式的双c形开槽谐振器和双c形金属谐振器以及金属地板构成，由于互补形谐振器具有很高的q值，所以使由单元构成的超表面能在低频比的条件下实现相位的独立调控。
[0062]
为进行验证并探索应用，设计工作于f1＝9.2ghz和f2＝11.2ghz的多功能超表面，并在工作频率f1与f2处分别实现了模式数l＝3的聚焦oam波束和零阶贝塞尔波束。通过在f1与f2工作频率处对超表面进行仿真测试，得出仿真结果与理论结果吻合良好，表明超表面在工作频率处具有良好的波前相位调控能力。与以往报道的多功能超表面相比，本专利所提出器件的双频工作频比仅为1.2，且器件在f1、f2处效率高达86.1％和93％。
[0063]
如附图1-2所示，该高效互补谐振器单元主要包括五部分：双c形开槽谐振器-第一层介质板-双c形金属谐振器-第二层介质板-金属地板的结构顺序组成。其中第一层介质板厚度h1＝1.5mm，第二层介质板厚度h2＝1.5mm，介质板材料为f4b(εr＝2.65，tanδ＝0.001)，该材料具有正切损耗小，色散小，制造难度低等优点。对于互补形金属谐振器以及金属地板，采用的材料为铜，其电导率为σ＝5.8
×
107s/m。由于互补形谐振器具有很高的q值，所以单元能在低频比的条件下实现相位的独立调控，避免因频带窄而造成通道间的串扰。为了单元能在两个谐振频率处有最大的反射幅度，通过对单元中谐振器各项参数的优化，最终得出如下结果。其中单元周期p＝10.2mm，双c形开槽谐振器结构参数为：外环外径r1＝
4.75mm、外环内径r2＝4.35mm、中间槽宽为w1＝0.3mm、内金属环宽w2＝0.4mm。其中外环与内环连接部分金属宽度g1＝0.9mm，双c形金属谐振器结构参数为：金属环内径宽度r3＝3.0mm、金属环宽度w3＝0.8mm。其中金属环缝隙宽度g2＝0.3mm。对于谐振器旋转角度而言，双c形开槽谐振器以所在介质板平面几何中心为原点在介质板平面按逆时针旋转角度为α1，双c形金属谐振器以所在介质板平面几何中心为原点在介质层平面按逆时针旋转角度为α2。
[0064]
根据所需功能对单元进行周期性排布以达到在f1＝9.2ghz与f2＝11.2ghz工作频率处分别实现模式数l＝3的聚焦oam波束和零阶贝塞尔波束。
[0065]
对于多功能超表面不同工作频率处功能的实现，需要选择合适的单元。单元对于线极化电磁波的调控，一般采用改变金属贴片的尺寸使反射电磁波发生传输相位的突变从而来实现对电磁波的控制。对于圆极化波而言，应用几何相位原理，通过旋转对应谐振频率f1与f2的谐振器来实现圆极化反射电磁波相位的调控。为进一步解释单元在圆极化波垂直入射下的高效工作原理，引入线极化波反射模式下的琼斯矩阵其中r
xx
，r
yy
为x、y极化波复反射系数，其中包含了反射波的振幅值与相位值如：
[0066]
由于几何相位原理是通过对单元中谐振器的旋转来使单元对应谐振频率的反射波产生相位突变，所以当谐振器旋转过一个固定角度时与之对应的坐标系也会发生改变。当谐振器逆时针旋转θ角时，旋转后的uov坐标系与xoy坐标系会有如下关系。
[0067][0068]
其中定义为旋转矩阵r(θ)，当谐振器绕中心轴左旋则θ》0，相反右旋θ《0。这里以右旋圆极化波垂直入射为例，因圆极化波能分解为两个相互正交、幅值相等且相位相差90
°
的线极化波，故右旋圆极化波可表示为
[0069]
其中，eo表示入射波的振幅。
[0070]
经谐振器旋转后，反射波在uov坐标系下表示为：
[0071][0072]
其中r
uu
为在u方向的振幅，r
vv
为在v方向的振幅，而则为在u方向引入的相移，而φ
vv
则为在v方向引入的相移，其中θ是旋转角度，由此在圆极化波激励下可以引起单元几何相位的产生。
[0073]
因为结构为反射型超表面，所以由电磁传播理论可知：r
uu
＝r
xx
，r
vv
＝r
yy
，由此反射波可以表示为：
[0074][0075]
[0076][0077][0078]
由以上结果可知反射波存在左旋圆极化波和右旋圆极化波，将其拆分为
[0079][0080][0081]
其中，为右旋圆极化波，为左旋圆极化波；
[0082]
从结果可以看出左旋圆极化部分携带了e
j2θ
相位因子，而当相位因子，而当时结果可以表示为：
[0083][0084][0085]
所以当右旋圆极化波垂直照射到超表面时，反射左旋圆极化波相对于入射波而言发生了极化转换且存在2θ的几何相位变化。高效互补谐振器单元在两个谐振频率f1和f2处均实现100％交叉圆极化波转化的独立几何相位调控，其中，双c形开槽谐振器工作于频率f1，负责该频率处的几何相位功能调控，双c形金属谐振器工作于频率f2，负责该频率处的几何相位功能调控，即在谐振频率f1和f2处线极化波入射条件下所述高效互补谐振器单元满足：
[0086]
|r
xx
|＝|r
yy
|＝1
[0087][0088]
其中，r
xx
为在x极化的反射系数，r
yy
为在y极化的反射系数，为在y极化的反射系数，为线极化波同极化反射相位，为实现圆极化波100％交叉极化转换需要单元满足以上条件。
[0089]
为验证设计单元满足圆极化波交叉极化反射条件，通过时域有限差分法(finite difference time domain,fdtd)进行仿真计算。仿真过程中，以x、y正交线极化波垂直照射单元，沿x、y方向的两个边界设置为周期边界条件，仿真结果如附图3所示。可以看出，x、y极化波入射下，谐振频率f1处主要为双c形开槽谐振器与金属地板参与谐振，而f2谐振频率处主要为双c形金属谐振器与金属地板参与谐振。在不同谐振频率f1与f2的反射场中r
xx
，r
yy
的反射幅度均大于0.92，均可视为接近于1。与此同时，x、y极化波的反射相位与在8～10ghz和10～12ghz的频段范围内相位差值保持均匀变化，并在f1与f2处稳定保持在180
°
左右，满足式(6)x、y线极化反射波激发下的幅度和相位条件，保证了圆极化波照射下单元的
高效率交叉极化反射。
[0090]
为了展示单元在谐振频率处具有良好的反射效果，通过仿真软件cst对单元进行仿真。通过对反射模式下单元的电流分布，如附图4所示，可以得出在谐振频率f1＝9.2ghz圆极化电磁波垂直入射时参与调控圆极化波的双c形开槽谐振器内存在闭合的谐振电流，且与底层金属地板在yoz平面内产生相反的电流分布从而发生了磁谐振现象。当谐振频率f2＝11.2ghz圆极化电磁波垂直入射时参与调控圆极化波的双c形金属谐振器内存在闭合的谐振电流，且与底层金属地板在xoz平面内产生相反的电流分布从而发生了磁谐振现象。
[0091]
为证明圆极化波入射下，高效互补形谐振器单元能在f1＝9.2ghz与f2＝11.2ghz处能实现几何相位的独立调控且交叉极化反射波能实现2π的几何相位覆盖，以右旋圆极化波垂直入射为例对单元进行仿真。设双c形开槽谐振器从初始位置逆时针旋转α1，双c形金属谐振器逆时针旋转α2，旋转步进为30
°
，仿真结果整理如附图5-6所示。
[0092]
如附图5所示，在谐振频率f1＝9.2ghz处，随着双c形金属谐振器旋转角度的改变，不同α1的双c形开槽谐振器的反射幅度趋近于1且反射相位几乎不随α2变化，同时其反射相位差与旋转角α1满足2倍几何相位关系。同理在谐振频率f2＝11.2ghz处，随着双c形开槽谐振器旋转角度的改变，不同α2的双c形金属谐振器的反射幅度与相位可以得到相同结论。通过对双c形开槽谐振器与双c形金属谐振器在谐振频率f1与f2处不同转角下的反射相位误差进行计算，得出其均小于8
°
，因此可忽略，证明单元能在不同两个谐振频率处具有良好的独立相控能力。
[0093]
对于互补谐振器单元的几何相位能否在f1与f2处达到2π的几何相位覆盖，从附图6可以得出。在谐振频率f1和f2处，随着双c形开槽谐振器的转角α1与双c形金属谐振器的转角α2从0
°
变化到180
°
，f1与f2处反射波的几何相位成功实现2π的几何相位覆盖，不同旋转角度的谐振器在工作频率处彼此之间的相位变化满足2倍转角关系，符合几何相位原理，且反射幅度在工作频率处高达0.98。以上结果充分表明高效互补谐振器单元有独立调控不同频率圆极化电磁波的能力，且双c形开槽谐振器与双c形金属谐振器在谐振频率f1与f2处彼此之间互不影响，这为设计多功能超表面打下坚实基础。
[0094]
本发明提出的双层频率多元反射超表面在所设计的单元在不同两个谐振频率处具有良好的相互独立性，且在不同两个谐振频率处单元均实现了2π的几何相位覆盖的前提下，通过对已知单元进行非周期性排列，最终构建了在以f1与f2为工作频率的，分别实现了模式数l＝3的聚焦oam波束和零阶贝塞尔波束的频率复用多功能超表面。
[0095]
对于多功能超表面在工作频率f1＝9.2ghz处实现模式数l＝3的聚焦oam涡旋波束，利用相位叠加原理将其分解为聚焦相位部分与涡旋相位部分，如附图7所示。
[0096]
对于聚焦相位而言，需要单元满足如下相位分布。
[0097][0098]
其中f0＝214mm为其焦距，p为单元周期，λ为入射电磁波的波长，是任意相位常量，m(n)为沿x(y)方向距原点的单元数目。只要使超表面的每个单元的反射相位满足公式，即可实现在预期焦平面上的聚焦效果。
[0099]
同理对于涡旋波束的相位分布，单元对入射波进行的相位补偿需要满足如下公
式：
[0100][0101]
其中l＝3为涡旋波束的模式数。依据电磁波的相位叠加原理，将涡旋相位与聚焦相位进行相位叠加，就可以实现在谐振频率处的目标功能。所以对于目标反射电磁波的最终相位分布满足如下公式。
[0102][0103]
本发明提出的双层频率多元反射超表面在工作频率f2＝11.2ghz处实现零阶贝塞尔波束，需要使超表面单元在谐振频率f2的补偿相位满足如附图8所示分布，其中单元具体相位值如下。
[0104][0105]
式中p为单元的周期，λ为入射电磁波的波长，β＝30
°
为贝塞尔波束的衍射半角。根据贝塞尔波束成波原理图，如附图9所示，可知零阶贝塞尔波束主要由与z轴正方向夹角为β的反射平行波在无衍射距离内叠加形成，因此通过图中几何关系可知最大无衍射距离z
max
为如下所示。
[0106][0107]
其中d与反射超表面的周期相等，通过计算可得z
max
＝309mm。从式中可以看出该方法生成的贝塞尔波束能量主要集中于半径为d/4的圆内，其中能量最高处在z
max
/2处，且波束能量随着距离超表面位置的变化呈现出能量密度先升高后降低的趋势。
[0108]
本发明提出的双层频率多元反射超表面通过cst-matlab的联合仿真对单元进行布阵与验证，包括以下过程：
[0109]
首先通过在全波仿真软件cst中对单元谐振器，介质板，金属地板等结构依据参数进行建模，并设置双c形开槽谐振器与双c形金属谐振器以所在介质板平面几何中心为圆心逆时针旋转角度分别为α1与α2。
[0110]
其次针对不同工作频率f1与f2与超表面所在工作频率处功能的相位分布，通过matlab对单元中双c形开槽谐振器与双c形金属谐振器进行旋转。通过几何相位调节机制，调整单元谐振器旋转角度并对单元模型进行非周期性排列最终实现了多功能超表面的建模。
[0111]
最终在通过仿真软件cst中以右旋圆极化平面波垂直入射超表面对其进行时域仿真。在工作频率f1处，通过对仿真数据进行整理，得到在焦平面f0＝214mm处的电场实部，相位分布图如附图10所示。图中电场实部图显示出电场明显聚焦于焦点处且波束具有与模式数l＝3相等的螺旋臂数目，而电场相位图显示出波束具有与模式数l数量一致的1080
°
相位变化，这与焦平面处理论聚焦oam波束的成波特性相吻合。针对近场仿真超表面的模式数l＝3的聚焦oam波束效率，采用如下方法进行计算：
[0112]
[0113]
其中定义为完全反射波功率与入射波功率之比，为焦平面处焦点功率与焦平面总功率之比，其中焦点功率为以焦点为圆心，主瓣半功率宽度为半径包围的圆功率积分，其积分区域如附图10黑色虚线区域。经计算可得最终计算得到聚焦oam波束的仿真效率为ηf＝87.5％。其效率未能达到100％的主要原因为部分电磁波因衍射未参与聚焦而使得聚焦效率较低。
[0114]
在工作频率f2处，通过对仿真数据进行整理，得到了yoz面和z1＝75、z2＝150、z3＝225和z4＝309mm四个位置处的|e
lcp
|^2归一化分布，如附图11所示。图中显示在波束传播方向的能量较为集中，在不同位置横截面的中心处都具有最大能量强度，且横向能量在偏离中心位置时逐渐振荡减小，符合贝塞尔函数曲线的振荡趋势。
[0115]
通过贝塞尔波束能量集中效率来对超表面的成波效果进行验证，选取z1＝75、z2＝150、z3＝225和z4＝309mm四个观察面的|e
lcp
|^2能量分布并通过下式进行计算。
[0116][0117]
其中，p2是超表面上方同一位置处同一平面内无衍射波束区域(平面中心为圆心，半径为d/4的圆形区域)，p1为同一平面内超表面口径面积(边长为d的方形区域，与超表面等大)；电磁波的能量大小由坡印廷矢量决定，即s＝e
×
h，其中s为坡印廷矢量，e为电场强度，h为磁场强度。
[0118]
最终得出在z1＝75、z2＝150、z3＝225和z4＝309mm四个位置处的纵向贝塞尔波束能量集中效率依次为82％，92％，85％和73％，这与贝塞尔波束能量随传输距离先略微增大而后略微减小的无衍射传输特性相符。为进一步证明器件的优越性能，提取z1＝75、z2＝150、z3＝225和z4＝309mm四个位置处x＝0的|e
lcp
|^2归一化强度分布，如附图12所示。由图可知，主瓣中|e
lcp
|^2强度随位置z的距离增大出现先增大后减小的趋势，且不同位置的副瓣也具有相同现象。这与零阶贝塞尔波束的成波特性完全相符。综上所述，通过对双层频率多元反射超表面模型仿真，其较好地实现了在f1＝9.2ghz与f2＝11.2ghz处模式数l＝3的聚焦oam波束与零阶贝塞尔波束。
[0119]
本发明还提出一种双层频率多元反射超表面设计方法，具体包括以下步骤：
[0120]
步骤1，设计双层频率多元反射超表面所需的高效互补谐振器单元，使高效互补谐振器单元在两个谐振频率处实现100％交叉圆极化波的相位调控；
[0121]
为实现双层频率多元反射超表面在不同工作频率f1与f2处波束的形成，设计了在不同谐振频率f1与f2处能实现圆极化波交叉极化高反射与独立2π相位调控的亚波长尺寸单元。其单元结构主要由三层金属和两层介质板组成，在第一层介质板上印刷有双c形开槽谐振器，其厚度h1＝1.5mm，第二层介质板上印刷有双c形金属谐振器，其介质厚度h2＝1.5mm，最底层印刷有金属地板，所有金属厚度为h＝0.036mm。对于谐振器与金属地板，采用的材料为铜，其电导率为σ＝5.8
×
107s/m，其目的为提高电磁响应效果，减少超表面的插入损耗。对于谐振器而言，双c形开槽谐振器外环外径r1＝4.75mm、外环内径r2＝4.35mm、中间槽宽为w1＝0.3mm和内金属环宽w2＝0.4mm，其中外环与内环连接部分金属宽度g1＝0.9mm。双c形金属谐振器结构参数为：金属环内径宽度r3＝3.0mm、金属环宽度w3＝0.8mm，其中金属环缝隙
宽度g2＝0.3mm。为了验证单元能在两个谐振频率处实现100％交叉圆极化波的相位调控，以几何相位理论为基础对单元结构条件进行理论分析并得出如下结论：
[0122]
高效互补谐振器单元在两个谐振频率f1和f2处均实现100％交叉圆极化波转化的独立几何相位调控，其中，双c形开槽谐振器工作于频率f1，负责该频率处的几何相位功能调控，双c形金属谐振器工作于频率f2，负责该频率处的几何相位功能调控，即在谐振频率f1和f2处线极化波入射条件下所述高效互补谐振器单元满足：
[0123]
|r
xx
|＝|r
yy
|＝1
[0124][0125]
其中，r
xx
为在x极化的反射系数，r
yy
为在y极化的反射系数，为在y极化的反射系数，为线极化波同极化反射相位。
[0126]
通过时域有限差分法(finite difference time domain,fdtd)以x，y正交线极化波垂直入射单元进行仿真验证，最终得出结果与要求高度吻合以此认为单元可以在不同两个谐振频率处实现高度的圆极化波交叉极化的调控。
[0127]
步骤2，通过高效互补谐振器单元非周期性排列构成了频率多元反射超表面；
[0128]
基于单元具有独立调控不同频率圆极化电磁波的能力，双层频率多元反射超表面最终构建了在以f1＝9.2ghz与f2＝11.2ghz为工作频率的，分别实现了模式数l＝3的聚焦oam波束和零阶贝塞尔波束的频率复用多功能超表面，其主要过程如下所示。
[0129]
为在f1＝9.2ghz处超表面实现模式数l＝3的聚焦oam波束，根据相位叠加原理将其分解为聚焦相位部分与涡旋相位部分。
[0130]
其中聚焦相位分布满足公式：
[0131][0132]
其中，f0＝214mm为其焦距，p为单元周期，λ为入射电磁波的波长，是任意相位常量，m为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目；
[0133]
对于涡旋波束的相位分布，单元对入射波进行的相位补偿需要满足如下公式：
[0134][0135]
其中，l＝3为oam波束的模式数；
[0136]
依据电磁波的相位叠加原理，将涡旋相位与聚焦相位进行相位叠加，对于双层频率多元反射超表面的最终单元相位分布满足如下公式：
[0137][0138]
表示单元在谐振频率f1＝9.2ghz处最终的反射相位。
[0139]
双层频率多元反射超表面在工作频率f2＝11.2ghz处实现零阶贝塞尔波束，所述高效互补谐振器单元在谐振频率f2的补偿相位满足如下公式：
[0140][0141]
其中，p为单元的周期，λ为入射电磁波的波长，β＝30
°
为贝塞尔波束的衍射半角，m
为沿x方向距原点的单元数目，n为沿y方向距原点的单元数目。
[0142]
为使超表面能更好地应用于工程实际，可以通过改变阵列排布来控制聚焦oam波束的焦距f0和模式数l，以及贝塞尔波束的衍射半角β来实现对贝塞尔波束最大无衍射距离的自由操控。
[0143]
通过matlab计算不同位置的单元在不同两个谐振频率处所需要的补偿相位，利用cst-matlab进行联合布阵。因为在不同两个谐振频率处单元谐振器的补偿相位与旋转角度满足2倍关系，所以将matlab在f1与f2处计算得到的双c形开槽谐振器与双c形金属谐振器旋转角度通过invoke函数撰写成cst中的vb命令，通过vb命令对超表面进行建模。在cst中对建模得到的超表面进行时域仿真，这里将边界条件均设置为open(add space)条件，将端口设置为以右旋圆极化平面波垂直入射。
[0144]
步骤3，对双层频率多元反射超表面性能分析，确认双层频率多元反射超表面实现预期功能；
[0145]
为对不同工作频率f1与f2处超表面的性能分析，计算了两类波束的效率，对于聚焦oam波束采用如下方式来进行计算：
[0146][0147]
其中定义为完全反射波功率与入射波功率之比，为焦平面处焦点功率与焦平面总功率之比，其中焦点功率为以焦点为圆心，主瓣半功率宽度为半径包围的圆功率积分，其积分区域如附图10黑色虚线区域。经计算可得最终计算得到聚焦oam波束的仿真效率为ηf＝87.5％。其效率未能达到100％的主要原因为部分电磁波因衍射未参与聚焦而使得聚焦效率较低。
[0148]
对于零阶贝塞尔波束，通过计算波束能量集中效率对超表面进行分析。
[0149][0150]
其中，p2是超表面上方同一位置处同一平面内无衍射波束区域(平面中心为圆心，半径为d/4的圆形区域)，p1为同一平面内超表面口径面积(边长为d的方形区域，与超表面等大)；电磁波的能量大小由坡印廷矢量决定，即s＝e
×
h，其中s为坡印廷矢量，e为电场强度，h为磁场强度。
[0151]
在工作频率f2处，通过对零阶贝塞尔波束的仿真数据进行整理，取z1＝75、z2＝150、z3＝225和z4＝309mm四个位置处yoz面的|e
lcp
|^2能量分布数据并进行计算。最终得出在z1＝75、z2＝150、z3＝225和z4＝309mm四个位置处的纵向贝塞尔波束能量集中效率依次为82％，92％，85％和73％，这与贝塞尔波束能量随传输距离先略微增大而后略微减小的无衍射传输特性相符。
[0152]
综上所述，多功能超表面仿真较好地实现了在f1＝9.2ghz与f2＝11.2ghz处模式数l＝3的聚焦oam波束与零阶贝塞尔波束，验证了设计双层频率多元反射超表面能在两个谐振频率处很好的实现预期功能，且在谐振频率f1，f2处的实验效率高达86.1％和93％，证明了多功能超表面的优越性能，同时也为oam波束产生，多功能集成器件设计提供了新途径。
[0153]
本发明不仅局限于上述具体实施方式，本领域一般技术人员根据实施例和附图公开内容，可以采用其它多种具体实施方式实施本发明，因此，凡是采用本发明的设计结构和思路，做一些简单的变换或更改的设计，都落入本发明保护的范围。