Controllable Manipulation of Wi-Fi Signals Using Tunable Metasurface
DOI: 10.12000/JR21012
-
Abstract: In this paper, we propose the utilization of a programmable metasurface for flexibly manipulating ambient Wi-Fi signals. First, we propose a new and efficient optimization algorithm CWGS (Complex Weighted Gerchberg-Saxton), which is based on an electromagnetic scattering model of the metasurface. The proposed algorithm quickly redesigns the complex amplitude distribution of the Wi-Fi field bounced off the programmable metasurface to enhance the Wi-Fi signals at desired locations significantly. Second, we fabricated a large-scale programmable metasurface that operates at the 2.4 GHz frequency band. We conducted several experiments using the fabricated metasurface to verify the proposed optimization algorithm’s feasibility and effectiveness. Both the theoretical and experimental results show that the programmable metasurface can dynamically boost Wi-Fi signals at multiple locations. Besides, we experimentally verified that using the developed strategy could improve the Wi-Fi signals by 23.5 dB. The results of our work improve the usability and practicality of the programmable metasurface in real-world applications and pave the way for wireless communications, future smart homes, and other applications.摘要: 可编程超表面是由可调谐的单元在二维平面上组成的平面阵列,具有任意、动态操控电磁波波前的能力,是微波领域前沿研究方向之一。尽管目前基于可编程超表面的电磁调控研究已取得丰硕进展,但是现有技术都需要采用专用发射源主动馈电超表面,这不仅增加了实际系统的复杂性和成本,而且在一定程度上限制了这些技术在现实环境中的应用。因此,该文提出了一种利用可编程超表面灵活调控周围Wi-Fi信号的方法,并通过理论和实验证明了其对Wi-Fi信号的优越调控性能。首先,该文首次提出一种基于可编程超表面散射模型的高效优化算法CWGS,该算法可以重新设计可编程超表面散射场的复振幅分布,从而实现指定位置处的Wi-Fi信号显著增强。其次,该文制作了一款工作频率为2.4 GHz的大规模可编程超表面,并基于该超表面实验验证了优化算法应用于Wi-Fi信号增强调控的可行性和有效性。理论和实验结果均证明,可编程超表面可以在多个位置处实现Wi-Fi信号的动态增强。实验结果表明,经可编程超表面调控后的Wi-Fi信号强度提高了23.5 dB。该文提出的方法提高了可编程超表面在实际应用中的可用性和实用性,有望为无线通信、未来智能家居等领域开辟新的道路。
-
关键词:
- 商用Wi-Fi 信号 /
- 可编码超表面 /
- 空间能量分配
-
1. 引言
合成孔径雷达(Synthetic Aperture Radar, SAR)可以实现全天时、全天候的高分辨率成像[1],星载SAR系统为全球地理遥感、地形测绘、形变监测等应用提供了丰富的高分辨数据。分布式星载SAR系统是一种特殊模式的合成孔径雷达系统,雷达发射机和接收机位于不同平台,这种收发分置的特点带来了很多优势,比如灵活的基线配置,多样的接收模式,可以提供更精细的地物分类信息等,具有单基情况无可比拟的优势。但是收发分置的特点也带来了很多新的挑战,尤其是时间同步、相位同步和波束同步等三大同步问题[2,3]。由于收发雷达不共频率源,导致了时间同步和相位同步问题[4],时间同步问题是指,发射雷达和接收雷达使用不同的晶振产生脉冲信号,导致了触发脉冲不一致,进而导致信号发射时刻与接收窗起始时刻不一致,使距离向的录取回波减少,时间同步误差会在回波中引入相位误差,并且晶振偏差会随时间累积,对回波数据产生较大影响,严重降低雷达图像质量[5,6]。为达到时间同步的要求,Tandem-X SAR系统采用频率一致性指标为10–8量级的晶振,并采用跳跃PRI的方法,确保600 s内,时间同步误差小于5 μs[7]。相位同步问题是指,由于晶振频率的差异,导致了发射信号调制频率和接收信号解调频率不一致,造成解调回波中存在沿方位时间变化的相位误差,该误差在成像后传递到图像域,导致方位向位置偏移、主瓣展宽、积分旁瓣比恶化等问题,并且相位误差最终会传递到干涉相位,造成干涉测高精度损失。TanDEM-X SAR系统采用脉冲交替对传的同步方案解决相位同步的问题。波束同步问题是指,由于收发天线的控制精度不同,使得收发天线无法在同一时刻完全覆盖地面同一区域,导致双基合成天线方向图的增益降低,同时接收天线无法获得全部的目标散射信号,造成有效观测幅宽损失。为保证系统系能,TanDEM-X SAR系统的波束指向精度控制在0.01°,最差情况下导致的天线增益损失小于0.1 dB。
为了消除相位同步带来的影响,可以建立专门的同步链路,文献[8]提出了卫星间连续波对传的同步方法,文献[9]和文献[10]进一步提出了脉冲交替对传的方法,文献[11]提出一种基于信号交换的调频连续波合成孔径雷达相位同步方法。无论采用哪种同步链路方式,目的都是获得相位误差信息,最后在进行数据处理时提取并补偿相位误差,这也是解决分布式卫星SAR相位同步的优选方法[12]。文献[13]建立了一种时间同步误差模型,分析了时间同步误差对成像和干涉测高性能的影响,但是没有分析时间同步和相位同步的联系。时间同步和相位同步来自同一误差源,都是由于收发平台采用了独立的晶振系统来校准系统时钟[14],所以两者之间存在一定的关系。波束同步的研究主要集中在波束同步控制方法,以及波束同步误差对雷达成像性能的影响,文献[15]分析了几种可行的飞行控制策略,同时分析了波束同步对收发平台和姿态控制精度的要求,文献[16]分析了系统波束控制精度需要满足的要求。但是,在不同的卫星姿态控制精度和波束指向精度情况下,波束同步误差对分布式星载SAR系统的综合影响还需要进一步分析。
本文结构安排如下。第2节介绍了时间同步误差模型和波束同步误差模型。第3节基于误差模型分析了时间同步和相位同步导致的相位误差之间的关系,并分析晶振稳定度需要满足的条件。第4节分析了波束同步误差对合成天线方向图、双基图像信噪比、多普勒去相干以及波束覆盖偏差的具体影响。第5节,进行仿真验证,为分布式SAR系统设计提供参考。
2. 误差模型
2.1 时间同步误差模型
收发雷达之间的时间同步误差模型可以表示为常数项、线性项和随机项3部分[13],随机项由频率源的短期稳定性决定,由于随机的特点,该项影响结果不具备积累性,一般只有在频率源稳定度较低或者所需相干积累时间较长时有较大影响[17–20],因此该模型可以表达为:
ε(t)=Δt+βt (1) 其中,
Δt 为常数项,表示初始状态主辅星间雷达累积的定时信号之差;βt 为线性项,由主辅星晶振频率偏差导致,β 可以表示为:β=ΔPRI1−ΔPRI2PRI=Δf1f1−Δf2f2≈Δf1−Δf2f0 (2) 其中,f0表示标称频率,
Δf1 和Δf2 分别表示发射雷达和接收雷达的频率偏移。当存在同步误差时,在以发射信号为参考的时间坐标系中,不同脉冲的回波在接收雷达采样窗中的位置会发生偏移[21]。2.2 波束同步误差模型
卫星姿态误差会影响卫星和地面目标的波束传播距离,同时影响回波多普勒频率,如图1所示,X轴从星下点指向卫星方向,Y轴指向方位向方向,Z轴指向距离向方向,卫星姿态误差可以分解为偏航角、俯仰角和横滚角,其中Rs表示斜距,R0表示最近斜距,P0表示姿态导引后期望的波束中心地面落点,P表示波束中心实际落点,
θy 表示偏航角,θp 表示俯仰角,θr 表示横滚角。可见,偏航角和俯仰角分别导致波束中心线绕X轴和Z轴旋转,进而偏离零多普勒面,均会导致斜视角。横滚角导致波束中心线绕Y轴旋转,所以横滚角引起波束在距离向移动,在距离向产生波束指向误差。根据几何关系可以得到,偏航角和俯仰角引起的斜视角分别为:
θsy=arctan(tanθy⋅sinθl) (3) θsp=arctan(tanθp⋅cosθl) (4) 其中,
θl 为下视角。横滚角引起的距离向波束指向误差为:θsr=θr (5) 3. 时间同步误差分析
主星发射的线性调频脉冲信号可以表示为:
s(τ,t)=rect(τTp)exp(jπkτ2)exp(jπf0τ) (6) 其中,
rect[⋅] 为矩形窗函数,τ 为距离向快时间,t为方位向慢时间,Tp为脉冲宽度,f0为发射信号载频,k为调频率。将时间同步误差模型代入发射信号,经辅星雷达接收解调后含有时间同步误差的回波信号为:sD(τ,t)=rect(τ−Rt/c−ε(t)Tp)⋅rect(t−tc)⋅exp(jπk(τ−Rt/c−ε(t))2)⋅exp(−j2πRt/λ)exp(−j2πf0ε(t)) (7) 其中,c为光速,
λ 为波长,tc为波束中心偏离时间,Rt为雷达和卫星之间的波束传播距离,R(t)=√R12+V21(t−t01)2+√R22+V22(t−t02)2 (8) 其中,R1和R2分别表示收发最近斜距,t01和t02分别表示波束中心穿越时刻,V1和V2分别表示收发卫星的速度。为简化分析,下文中,假设两星同时穿越波束中心时刻且以该时刻为方位时间零点,即t01=t02=0。
以主星发射信号,辅星接收信号为例,采用距离多普勒聚焦算法,推导时间同步误差对SAR成像的影响,假设主星雷达晶振为标准频率源,辅星雷达晶振含有误差项,即时间同步误差全部表现在辅星获得的图像上。
3.1 固定时间误差分析
辅星雷达只存在固定时间误差的情况下,即
ε(t)=Δt (9) 成像处理之后的辅星图像将存在时间同步误差项,具体可以表示为:
σ(τ,t)=sinc(B(τ−(R1+R2)/c−Δt))⋅sinc(Bd⋅t)⋅exp(−j2πR1+R2λ)⋅exp(j2πfact)⋅exp(−j2πf0Δt) (10) 其中,Bd是多普勒带宽,fac是多普勒中心频率。由式(10)可以得到,固定时间同步误差在辅星图像中引入距离向主瓣偏移和固定相位误差,固定相位误差在经过干涉处理之后会产生常量干涉相位误差。主瓣偏移量和干涉相位误差分别为:
Δr=Δt⋅c2 (11) φe=−2πf0Δt (12) 常量相位误差会影响绝对测高精度。
3.2 时间同步与相位同步关系
辅星雷达只存在线性时间误差的情况下,即
ε(t)=βt (13) 对式(7)中带有线性时间误差的相位项进一步推导可以得到
φlinear=−2πβf0t=−2πΔf1−Δf2f0f0t=−2π(Δf1−Δf2)t (14) 由于晶振频率存在偏移,t时刻主星的晶振频率为
f1=f0+Δf1 ,辅星的晶振频率为f2=f0+Δf2 ,根据文献[9]中对脉冲对传相位同步方案的分析,辅星需要补偿的相位偏差为:φc(t)=2π(Δf2−Δf1)t+π(Δf2−Δf1)(τ12+τsys)−πfdτsys+12(nφ2(t+τsys)+nφ2(t+τ)−nφ1(t)−nφ1(t+τ+τsys))+12(φSNR1(t)−φSNR2(t))+Δφsys(t)+Δφant(t) (15) 其中,fd是多普勒频率,
τsys 是同步脉冲的发射延迟,nφ1 和nφ2 分别代表晶振的相位噪声,φSNR1 和φSNR2 分别代表主星和辅星的接收机噪声,Δφsys(t) 代表通道链路的相位噪声,Δφant(t) 代表天线方向图引入的相位噪声。式(15)的第1项
2π(Δf2−Δf1) ,等于式(7)中的线性偏差φlinear ,即随方位时间线性变化的相位误差部分可以通过相位同步方案补偿,补偿后的回波信号经过距离徙动校正后在方位频域的表达形式为:SR−com(τ,fa)=sinc(B(τ−R1+R2c+βfaKa))⋅Wa(fa−fac)⋅exp(−j2π(R1+R2)/λ)⋅exp(jπf2aKa) (16) 其中,Ka为方位向调频率。式(16)中
β(fa/faKaKa) 项会造成距离徙动校正误差,使方位信号弥散在若干距离单元上,导致主瓣展宽。如果要求距离徙动校正误差造成的主瓣展宽小于2%,那么要求未被校正的徙动量小于半个分辨率单元[22],即max(βKa⋅fa)<0.5⋅0.886B (17) 根据时频关系可以得到
βKa⋅Ka⋅Tsar<0.5⋅0.886B⇒β<0.443B⋅Tsar (18) 其中Tsar为合成孔径时间。
满足式(18)的情况下,经方位压缩后的SAR信号为:
σ(τ,t)=sinc(B(τ−R1+R2c))⋅sinc(Bd⋅t)⋅exp(−j2πR1+R2λ) ⋅exp(j2πfact) (19) 此时主要导致主瓣展宽,严重情况可能导致散焦。
4. 波束同步误差分析
理想情况下,SAR系统处于匀速运动状态,但是在实际工作情况下,卫星的运动受到控制系统精度的影响,会产生姿态误差,所以必须保证分布式系统的姿态控制精度[23]。波束同步误差可以分为两个方面,一方面是卫星姿态引起的误差,另一方面是波束指向误差。波束同步误差会导致天线照射范围不一致,导致信噪比损失,多普勒带宽偏移,方位向天线增益降低等问题[24]。同时,由于照射区域的偏差,还会导致距离向重叠幅宽减小。
4.1 对方位向影响
本节首先引入文献[7]中介绍的基于椭圆轨道的双星一发双收工作模式下的姿态导引方式,在此基础上,分析波束同步误差对多普勒去相干、天线增益以及双基图像信噪比损失的影响。
双星有两种姿态导引方式[7],一种是只针对主星进行零多普勒姿态导引,此时,场景内的回波多普勒中心频率不为0。另一种是主辅双星均进行零多普勒姿态导引,此时辅星与主星照射的共同区域等效为正侧视,即场景内多普勒中心频率为0。图2和图3给出了在主辅星顺轨基线从0到3 km情况下,不同姿态导引方式下方位向合成天线方向图的情况,其它仿真参数如表1所示。
图2中,由于仅主星进行零多普勒姿态导引,导致合成后的方位天线方向图指向等效于斜视,即场景回波的多普勒中心频率不为0,产生方位多普勒去相干,导致相干系数发生变化,如图4所示。
考虑到顺轨基线和波束同步误差的综合影响,相干系数的计算方法可以表示为:
γdc=1−|ΔfDop1+ΔfDop2Bd| (20) ΔfDop1 表示由顺轨基线导致的多普勒频率偏移,计算公式为[7]:ΔfDop1≈v⋅dalongλ⋅r (21) 其中,dalong为基线长度,r为平均斜距。
ΔfDop2 表示由波束指向误差和卫星姿态误差导致的多普勒频率偏移,ΔfDop2=2v⋅sin(θsq+Δθa)λ−2v⋅sin(θsq)λ (22) Δθa=θerr+θsy+θsp (23) 其中,
θsq 表示斜视角,Δθa 表示方位向波束同步误差,包括方位向波束指向误差以及偏航角和俯仰角对斜视角的影响部分。图3中,由于对主辅双星均进行零多普勒姿态导引,合成后的方位天线方向图等效于正侧视,此时对多普勒去相干影响很小。考虑到主辅双星之间顺轨基线的影响,合成后的天线方向图幅度会有一定的下降,从仿真中可以看出,在顺轨基线小于3 km的情况下,导致的幅度衰减小于0.05 dB,可以忽略。在不同顺轨基线情况下,主辅星同时进行零多普勒导引的合成天线方向图对频谱加权的计算公式为[7]:
Wbistat(fDop;ΔfDop)=sinc(l2v(fDop−ΔfDop))⋅sinc(l2v(fDop+ΔfDop)) (24) 其中,l为天线长度,v为卫星速度,
ΔfDop 是由顺轨基线导致的多普勒频率偏移。在此基础上,加入波束指向误差和卫星姿态误差的影响,顺轨基线和波束同步误差的综合加权公式可以表示为:
Wbistat(fDop;ΔfDop)=sinc(l2v(fDop−ΔfDop1−ΔfDop2))⋅sinc(l2v(fDop+ΔfDop1+ΔfDop2)) (25) 由于波束同步误差,导致成像后辅星图像信噪比降低,信噪比损失量为[25]:
SNRloss=∫θbw1/2−θbw1/2(sinc2(l⋅sinφ/λ)⋅sinc2(l⋅sin(φ−Δθa)/λ))dφ∫θbw1/2−θbw1/2(sinc2(l⋅sinφ/λ)⋅sinc2(l⋅sinφ/λ))dφ (26) 其中,
φ 为方位向离轴角,θbw1 为方位向波束宽度。4.2 对距离向影响
针对距离向,波束控制误差主要影响主、辅星的重叠幅宽,当主、辅星距离向同时存在相反方向的指向偏差时,距离向的重叠区域(图5和图6中黄色区域)会大幅减小,影响后续的干涉处理。根据几何模型,可以推导出覆盖区域偏差的计算公式:
ΔW=(H+Re)cos(θl+θbw2/2−Δθr)−√(H+Re)2cos2(θl+θbw2/2−Δθr)−2⋅H⋅Re−(H+Re)cos(θl+θbw2/2+Δθr)+√(H+Re)2cos2(θl+θbw2/2+Δθr)−2⋅H⋅Re (27) 其中,H为卫星高度,Re为地球半径,
θl 为下视角,θbw2 为俯仰向波束宽度,Δθr 为距离向波束同步误差。5. 仿真结果
仿真参数如表1所示。
表 1 仿真参数表Table 1. Simulation parameters参数 数值 载频(GHz) 1.30 卫星高度(km) 700 方位向天线长度(m) 6 脉冲宽度(μs) 80 距离向带宽(MHz) 100 PRF(Hz) 1600 卫星有效速度(m/s) 7500 距离向波束指向精度(°) 0.06 方位向波束指向精度(°) 0.05 姿态控制精度(°) 0.02 顺轨基线(km) 3 5.1 时间同步误差
仿真中加入参数为
Δt=100 ns 的固定时间误差,点目标分析结果如图7所示,通过计算可以得到理论上的距离向位置偏移为Δr=Δt⋅c/2=15 m ,仿真得到的对应偏差为15 m,如图7(b)所示。时间同步误差会引起主瓣展宽,如式(18)所示,仿真参数下,当
β<1.2×10−9 时,引起的展宽可以忽略。当β=5×10−8 时,距离向和方位向剖面图如图8所示,由仿真结果可以看到明显的主瓣展宽。5.2 波束同步误差
本节的仿真给出了当顺轨基线为3 km时,波束指向误差和卫星姿态误差对系统性能的影响。图9和图10分别比较了方位向波束指向精度在0.01°到0.05°情况下,对合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.05°波束指向误差造成方位天线增益下降小于0.05 dB,造成相干系数减小0.02。
图11和图12分别给出了卫星姿态控制精度在0.005°到0.020°的情况下,对方位向合成天线方向图及相干系数的影响,仿真结果表明,与无误差情况对比,0.020°的姿态控制误差造成方位天线增益下降小于0.02 dB,造成相干系数减小0.01。
不同波束指向和卫星姿态误差情况下,对合成天线方向图及相干系数的综合影响结果如图13和图14所示。波束指向精度为0.05°,卫星姿态控制精度0.02°的情况下,和无误差情况对比,造成方位天线增益下降小于0.1 dB,造成相干系数减小0.03。
不同波束指向和卫星姿态误差情况下,辅星图像信噪比损失情况如表2所示,0.05°的波束指向误差和0.02°的姿态误差导致的辅星图像信噪比损失小于0.1 dB。
表 2 波束同步误差造成的辅星图像信噪比损失(dB)Table 2. Influence of beam synchronization error on SNR (dB)波束指向精度(°) 姿态控制精度(°) 0.01 0.02 0.01 0.021 0.025 0.02 0.028 0.032 0.03 0.036 0.039 0.04 0.046 0.048 0.05 0.055 0.057 0°~0.06°的波束指向误差对距离向重叠幅宽的影响如图15所示。横滚角会在距离向产生波束控制误差,对距离向覆盖偏差产生影响,0.02°的姿态控制误差对重叠幅宽的影响如图16所示。横滚角从0°步进至0.02°,距离向波束指向精度从0°步进至0.06°情况下,两者对距离向重叠幅宽的综合影响如图17所示。上述结果表明,0.02°的姿态控制误差对距离向重叠幅宽的影响为368 m, 0.06°的波束指向误差对距离向重叠幅宽的影响为1100 m,由于波束指向和姿态控制导致的误差均为随机误差,所以综合影响采用几何平均数的形式,两者综合之后对距离向重叠幅宽的影响为1164 m。考虑到轨道控制精度的影响,在重轨工作模式下,两次测量中卫星的位置会在一定范围内偏移,在卫星姿态误差和波束指向误差基础上,系统设计中应综合考虑卫星位置偏移的影响。
由上述分析可知,在仿真精度情况下,波束同步精度对多普勒去相干、天线增益以及双基图像信噪比损失的影响较小,但对距离向重叠幅宽影响较大,因此需要在波位设计中保留一定的幅宽余量,或者根据幅宽需求,对波束同步精度提出要求。不同卫星姿态控制精度和波束指向精度对重叠幅宽的综合影响总结在表3可供系统设计参考。针对不同的卫星系统,可以根据系统参数结合式(27)求得对波束同步误差的限制。
表 3 波束同步误差对重叠幅宽的影响(m)Table 3. Influence of beam synchronization error to overlapping area (m)波束指向精度(°) 姿态控制精度(°) 0.01 0.02 0.01 260 411 0.02 411 520 0.03 582 663 0.04 759 823 0.05 938 991 0.06 1119 1164 6. 结论
本文分析了分布式SAR系统中时间同步误差和波束同步误差的影响,根据固定项和线性项的时间同步误差模型,推导时间同步和相位同步的关系,并分析线性项需要满足的条件。波束同步部分,本文分析了波束指向误差和卫星姿态误差对方位向合成天线方向图、多普勒去相干以及双基图像信噪比损失情况的影响,并根据几何关系推导距离向重叠幅宽的损失量。根据仿真参数,给出实验结果,以便在系统设计中充分考虑同步问题的影响,提出指标要求,为系统设计提供指导。
-
Figure 4. Results comparison between prediction and measurement
The first row are the coding sequences. The second to fourth rows are the normalized spatial intensity distributions predicted respectively at z=0.984 m, z=1.257 m, and z=1.531 m using the proposed method. The fifth row are the measured normalized spatial intensity distributions at z=1.531 m using the near-field scanning technology
Table 1. The focus efficiency for six different cases
Cases
(dB)Center
(%)Left
(%)Right
(%)Up
(%)Down
(%)Double
(%)–3 25.36 30.73 40.24 28.51 35.27 21.27 –6 37.42 45.19 59.32 41.81 51.96 38.65 Table 2. Comparing with the existing works dealing with EM waves manipulation using metasurface
Ref. f0 (GHz) Source type Working mode of metasurface Optimization
goal of algorithmWorking
regionThe shape
of energy
allocationReconfig-urability Maximum
focusing
efficiency[6] 6.9 Dedicated source Active Field’s amplitude pattern Near-field Point-like YES – [7] 10 Dedicated source Active Field’s amplitude pattern Near-field Point-like NO 57.8% [8] 5.75 Dedicated source Active Field’s amplitude pattern Near-field Point-like YES – [9] 10 Dedicated source Active Field’s amplitude pattern Near-field Point-like YES – [10] 20 Dedicated source Active Field’s amplitude pattern Near-field Point-like NO 52.3% [11] 10 Dedicated source Active Field’s amplitude pattern Near-field Point-like NO – [12] 5.8 Dedicated source Active Field’s amplitude pattern Near-field Point-like NO – [31] 14.3 Dedicated source Active Field’s amplitude pattern Far-field Single beam YES – [23] 5.8 Dedicated source Active Field’s amplitude pattern Near-field Point-like NO – [24] 2.3 Dedicated source Active Field’s amplitude pattern Far-field Single beam YES – [25] 10 Dedicated source Active Field’s amplitude pattern Far-field Single beam YES Proposed 2.4 Commodity wireless device Passive Field’s complex amplitude pattern Not limited Not limited YES 59.32% -
[1] YANG Huanhuan, CAO Xiangyu, YANG Fan, et al. A programmable metasurface with dynamic polarization, scattering and focusing control[J]. Scientific Reports, 2016, 6: 35692. doi: 10.1038/srep35692 [2] LI Lianlin, CUI Tiejun, JI Wei, et al. Electromagnetic reprogrammable coding-metasurface holograms[J]. Nature Communications, 2017, 8(1): 197. doi: 10.1038/s41467-017-00164-9 [3] LAROUCHE S, TSAI Y J, TYLER T, et al. Infrared metamaterial phase holograms[J]. Nature Materials, 2012, 11(5): 450–454. doi: 10.1038/nmat3278 [4] LI Lianlin, RUAN Hengxin, LIU Che, et al. Machine-learning reprogrammable metasurface imager[J]. Nature Communications, 2019, 10(1): 1082. doi: 10.1038/s41467-019-09103-2 [5] LI Lianlin and CUI Tiejun. Information metamaterials – from effective media to real-time information processing systems[J]. Nanophotonics, 2019, 8(5): 703–724. doi: 10.1515/nanoph-2019-0006 [6] CHEN Ke, FENG Yijun, MONTICONE F, et al. A reconfigurable active huygens’ metalens[J]. Advanced Materials, 2017, 29(17): 1606422. doi: 10.1002/adma.201606422 [7] WANG Yue, GUAN Chunsheng, DING Xumin, et al. Multi-focus hologram utilizing Pancharatnam-Berry phase elements based metamirror[J]. Optics Letters, 2019, 44(9): 2189–2192. doi: 10.1364/OL.44.002189 [8] WANG Zehao, LIAO Dashuang, ZHANG Ting, et al. Metasurface-based focus-tunable mirror[J]. Optics Express, 2019, 27(21): 30332–30339. doi: 10.1364/OE.27.030332 [9] ZHANG Kuang, YUAN Yueyi, DING Xumin, et al. High-efficiency metalenses with switchable functionalities in microwave region[J]. ACS Applied Materials & Interfaces, 2019, 11(31): 28423–28430. [10] GOWDA V R, IMANI M F, SLEASMAN T, et al. Focusing microwaves in the fresnel zone with a cavity-backed holographic metasurface[J]. IEEE Access, 2018, 6: 12815–12824. doi: 10.1109/ACCESS.2018.2802379 [11] ZHANG Pei, LI Long, ZHANG Xuanming, et al. Design, measurement and analysis of near-field focusing reflective metasurface for dual-polarization and multi-focus wireless power transfer[J]. IEEE Access, 2019, 7: 110387–110399. doi: 10.1109/ACCESS.2019.2934135 [12] YU Shixing, LIU Haixia, and LI Long. Design of near-field focused metasurface for high-efficient wireless power transfer with multifocus characteristics[J]. IEEE Transactions on Industrial Electronics, 2019, 66(5): 3993–4002. doi: 10.1109/TIE.2018.2815991 [13] ZHAO Jie, YANG Xi, DAI Junyan, et al. Programmable time-domain digital-coding metasurface for non-linear harmonic manipulation and new wireless communication systems[J]. National Science Review, 2019, 6(2): 231–238. doi: 10.1093/nsr/nwy135 [14] ZHANG Xu, GRAJAL J, VAZQUEZ-ROY J L, et al. Two-dimensional MoS2-enabled flexible rectenna for Wi-Fi-band wireless energy harvesting[J]. Nature, 2019, 566(7744): 368–372. doi: 10.1038/s41586-019-0892-1 [15] WANG Guanhua, ZOU Yongpan, ZHOU Zimu, et al. We can hear you with Wi-Fi![J]. IEEE Transactions on Mobile Computing, 2016, 15(11): 2907–2920. doi: 10.1109/TMC.2016.2517630 [16] ALI K, LIU A X, WANG Wei, et al. Keystroke recognition using WiFi signals[C]. The 21st Annual International Conference on Mobile Computing and Networking, Paris, France, 2015: 90–102. [17] WU Xuangou, CHU Zhaobin, YANG Panlong, et al. TW-see: Human activity recognition through the wall with commodity Wi-Fi devices[J]. IEEE Transactions on Vehicular Technology, 2019, 68(1): 306–319. doi: 10.1109/TVT.2018.2878754 [18] XU Qinyi, CHEN Yan, WANG Beibei, et al. Radio biometrics: Human recognition through a wall[J]. IEEE Transactions on Information Forensics and Security, 2017, 12(5): 1141–1155. doi: 10.1109/TIFS.2016.2647224 [19] HOLL P M and REINHARD F. Holography of Wi-Fi radiation[J]. Physical Review Letters, 2017, 118(18): 183901. doi: 10.1103/PhysRevLett.118.183901 [20] ZHONG Wei, HE W K, WANG Longgang, et al. Through-the-wall imaging using Wi-Fi signals[C]. The 2018 12th International Symposium on Antennas, Propagation and EM Theory, Hangzhou, China, 2018: 1–3. [21] LI Lianlin, SHUANG Ya, MA Qian, et al. Intelligent metasurface imager and recognizer[J]. Light: Science & Applications, 2019, 8: 97. [22] TANG Wankai, DAI Junyan, CHEN Mingzheng, et al. Programmable metasurface-based RF chain-free 8PSK wireless transmitter[J]. Electronics Letters, 2019, 55(7): 417–420. doi: 10.1049/el.2019.0400 [23] GOWDA V R, YURDUSEVEN O, LIPWORTH G, et al. Wireless power transfer in the radiative near field[J]. IEEE Antennas and Wireless Propagation Letters, 2016, 15: 1865–1868. doi: 10.1109/LAWP.2016.2542138 [24] DAI Linglong, WANG Bichai, WANG Min, et al. Reconfigurable intelligent surface-based wireless communications: Antenna design, prototyping, and experimental results[J]. IEEE Access, 2020, 8: 45913–45923. doi: 10.1109/ACCESS.2020.2977772 [25] RATNI B, DE LUSTRAC A, PIAU G P, et al. Reconfigurable meta-mirror for wavefronts control: Applications to microwave antennas[J]. Optics Express, 2018, 26(3): 2613–2624. doi: 10.1364/OE.26.002613 [26] XU Hexiu, SUN Shulin, TANG Shiwei, et al. Dynamical control on helicity of electromagnetic waves by tunable metasurfaces[J]. Scientific Reports, 2016, 6: 27503. doi: 10.1038/srep27503 [27] WANG Di, YIN Lizheng, HUANG Tiejun, et al. Design of a 1 bit broadband space-time-coding digital metasurface element[J]. IEEE Antennas and Wireless Propagation Letters, 2020, 19(4): 611–615. doi: 10.1109/LAWP.2020.2973424 [28] YANG Xue, XU Shenheng, YANG Fan, et al. A novel 2-bit reconfigurable reflectarray element for both linear and circular polarizations[C]. 2017 IEEE International Symposium on Antennas and Propagation & USNC/URSI National Radio Science Meeting, San Diego, USA, 2017: 2083–2084. [29] KAMODA H, IWASAKI T, TSUMOCHI J, et al. 60-GHz electronically reconfigurable large reflectarray using single-bit phase shifters[J]. IEEE Transactions on Antennas and Propagation, 2011, 59(7): 2524–2531. doi: 10.1109/TAP.2011.2152338 [30] LI Long, YANG Yang, and LIANG Changhong. A wide-angle polarization-insensitive ultra-thin metamaterial absorber with three resonant modes[J]. Journal of Applied Physics, 2011, 110(6): 063702. doi: 10.1063/1.3638118 [31] YANG Huanhuan, YANG Fan, CAO Xiangyu, et al. A 1600-element dual-frequency electronically reconfigurable reflectarray at X/Ku-Band[J]. IEEE Transactions on Antennas and Propagation, 2017, 65(6): 3024–3032. doi: 10.1109/TAP.2017.2694703 [32] WANG Ke, ZHAO Jie, CHENG Qiang, et al. Broadband and broad-angle low-scattering metasurface based on hybrid optimization algorithm[J]. Scientific Reports, 2014, 4: 5935. [33] GUAN Chunsheng, WANG Zhuochao, DING Xumin, et al. Coding Huygens’ metasurface for enhanced quality holographic imaging[J]. Optics Express, 2019, 27(5): 7108–7119. doi: 10.1364/OE.27.007108 -