CN104733850B - 基于全息调制的人工张量阻抗表面天线及其实现方法 - Google Patents

基于全息调制的人工张量阻抗表面天线及其实现方法 Download PDF

Info

Publication number
CN104733850B
CN104733850B CN201510185161.7A CN201510185161A CN104733850B CN 104733850 B CN104733850 B CN 104733850B CN 201510185161 A CN201510185161 A CN 201510185161A CN 104733850 B CN104733850 B CN 104733850B
Authority
CN
China
Prior art keywords
mrow
msub
impedance
mtd
mfrac
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510185161.7A
Other languages
English (en)
Other versions
CN104733850A (zh
Inventor
欧阳骏
龙宇
张袁
张凯之
隆锐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510185161.7A priority Critical patent/CN104733850B/zh
Publication of CN104733850A publication Critical patent/CN104733850A/zh
Application granted granted Critical
Publication of CN104733850B publication Critical patent/CN104733850B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerials With Secondary Devices (AREA)

Abstract

本发明公开了一种基于全息调制的人工张量阻抗表面天线及其实现方法,系统包括由多个大小相等的单元晶格人工张量阻抗表面;所述的单元晶格从上至下依次包括金属地板、介质基体和金属贴片;通过改变金属贴片的几何参数控制人工张量阻抗表面的单元晶格所对应的表面阻抗,从而使天线表面上的表面阻抗依照源场和目标场的全息干涉图样规则分布;方法包括以下步骤:S1:标量阻抗提取仿真;S2:等效标量表面阻抗曲线拟合;S3:张量表面阻抗全息调制;S4:天线建模仿真。本发明目的在于利用全息张量表面阻抗调原理获得特定计划方式的目标辐射场的同时,还提出了一种张量表面阻抗的等效标量表面阻抗拟合准则,有效减少张量表面阻抗调制建模所需的数据量。

Description

基于全息调制的人工张量阻抗表面天线及其实现方法
技术领域
本发明涉及一种基于全息调制的人工张量阻抗表面天线及其实现方法。
背景技术
在无线通信技术蓬勃发展的今天,各类无线通信应用已经完全融入了人们日常生活,成为不可或缺的一部分。在军事领域,以雷达为代表的天线系统,对军队作战能力也起着至关重要的作用。而作为发送和接收电磁波的核心部件,天线的电磁性能和物理特性对整个无线通信系统来说是至关重要的。以战机机载天线为例,除了对天线的增益、极化等传统的天线指标有较高或特定要求外,还需要天线具有剖面低、便于共形集成和重量轻等特性。现有的一些天线技术例如反射面天线、平板反射阵天线、阵列天线等等,往往存在诸如剖面较高、空间占用大和馈电复杂等缺点,不利于与载体集成。
作为光学全息技术在微波领域的延伸和拓展,全息天线以及基于全息原理的阻抗调制表面天线在设计方法上沿袭了物理光学全息技术的思路,即通过记录波束干涉信息和源场照射来获得目标漏波辐射场;在加工制作方面,全息阻抗调制表面天线多采用与传统微带天线相同的加工工艺和制作流程,不需要复杂的馈电网络,也不需要空间照射馈源,无论是在加工精度还是在成本控制方面,都能很好地控制。也正是因为这种简单的结构,使得全息阻抗调制表面天线在微波乃至毫米波频段具有很强的潜在优势。并且其重量较轻且剖面低,在载体共形集成方面的应用前景广阔。
通常,标量阻抗调制表面只能将线极化的源场转换成线极化的目标场,要实现对天线极化特性的控制,例如利用线极化表面波实现圆极化的远区辐射场,需要更复杂的张量阻抗调制表面。在现有的技术中,中国专利申请201410343597.X公开了全息人工阻抗表面共形天线,201310323806.X和201310292298.3则分别提出了基于全息表面阻抗调制的锥形波束和多波束天线。但这三种天线无一例外地都使用了标量表面阻抗,其表面阻抗值是不具有方向性的标量数值,可变参量只有金属贴片间距这一项。这就使得这些形式的天线难以实现对具有多个电场分量的目标场进行阻抗调制,其极化方式固定,灵活性不足。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于全息调制的人工张量阻抗表面天线及其实现方法,解决现有技术使用了标量表面阻抗、可变参量只有金属贴片间距这一项使得难以实现对具有多个电场分量的目标场进行阻抗调制,极化方式固定,灵活性不足的问题。
本发明的目的是通过以下技术方案来实现的:基于全息调制的人工张量阻抗表面天线,它包括人工张量阻抗表面,所述的人工张量阻抗表面包括多个大小相等的单元晶格;所述的单元晶格从上至下依次包括金属地板、介质基体和金属贴片;通过改变金属贴片的几何参数控制人工张量阻抗表面的单元晶格所对应的表面阻抗,从而使天线表面上的表面阻抗依照源场和目标场的全息干涉图样规则分布;所述的源场是激励源照射在天线表面上的表面波,所述的目标场是期望得到的平面波辐射场。
基于全息调制的人工张量阻抗表面天线的实现方法,它包括以下步骤:
S1:标量阻抗提取仿真,包括以下子步骤:
S11:在电磁仿真软件Ansoft HFSS中,使用单元晶格和周期边界条件模拟无限大均匀阻抗表面,所述的周期边界条件为主从边界;
S12:对使用不同组合的主从边界的相位差,模拟不同的表面波传播方向θt,从而得到张量表面阻抗的等效标量表面阻抗Ze随表面波传播方向θt变化的等效标量表面阻抗曲线;所述的主从边界的相位差表示从边界相对于主边界的相位延迟;
S2:等效标量表面阻抗曲线拟合:对不同贴片间距g和不同表面波传播方向θt的单元晶格模型进行仿真,得到一系列不同贴片间距g值的等效标量表面阻抗曲线,从而进一步提炼出Zemax和g之间的数学关系,为全息阻抗调制做好准备;其中,将等效标量表面阻抗曲线的两个极轴作为主要的研究对象,针对表面阻抗值和贴片参数g之间的映射提出了一种简化准则:对于开缝金属贴片的阻抗表面来说,在开缝大小固定的情况下,可以认为等效标量表面阻抗曲线的主轴方向θemax与贴片开缝角度θs一致即θemax=θs,且主轴所对应的等效标量表面阻抗最大值Zemax仅与单元晶格之间的贴片间距g有关;
S3:张量表面阻抗全息调制,包括以下子步骤:
S31:构造全息张量表面阻抗调制公式,所述的全息张量表面阻抗调制公式满足张量阻抗矩阵Z是反Hermit纯虚数二阶矩阵;
S32:计算源场的表面电流表达式和目标场的电场表达式;所述的目标场的电场表达式与期望的极化方式相关;
S33:联合步骤S31和步骤S32中得到的三个公式,计算张量表面阻抗Z的各阻抗分量的全息调制公式;
S34:根据各阻抗分量的公式,直接计算得到理论的等效标量表面阻抗曲线,并进而得到的等效标量表面阻抗最大值Zemax和主轴方向θemax
S35:结合步骤S2中得到的Zemax和g之间的数学关系,求得该点处单元晶格所需的贴片间距g和贴片开缝角度θs,完成天线模型的建立;
S4:天线建模仿真:采用单极子作为整个天线的表面波激励源;在天线平面的坐标原点处打孔放置单极子天线。
所述的步骤S12包括以下子步骤:
S121:在单元晶格模型上的四个壁上设置两对主从边界,其主从边界相位差分别设置成φx和φy,可得到:
kx=φx/a,ky=φy/a;
式中,kx和ky分别表示表面波在x方向和y方向传播的波数,a表示表面阻抗单元晶格的横向尺寸;
S122:将表面波在阻抗表面上传播的波矢量kt和与x轴的夹角θt进行表示:
S123:定义阻抗表面的等效折射率n,表示为:
式中,c=3×108m/s自由空间中的光速;ω是在该贴片几何参数下,阻抗表面的本征频率;k0是自由空间中的波数;vp和kt分别是表面波沿着阻抗表面传播的相速和波数;
S124:定义标量表面阻抗:
式中,η0=120π是自由空间中的波阻抗,αz是阻抗表面法向方向即z方向上表面波传播的衰减常数,j是虚数,j的平方为-1;
S125:选取合适的kt,通过φy和φx的不同组合,便能得到在该表面波传播方向上的本征频率ω,进而得到该方向上的标量表面阻抗值:
S126:张量表面阻抗Z与电场E和表面电流J的关系为:
TM模,其电场表达式可以写成:
S127:联合以上方程解得张量表面阻抗的归一化等效标量表面阻抗为:
S128:再得到三组不同的αz/k0之后,又可以将其代入步骤S127中的公式,反过来求得张量表面阻抗矩阵Z,进而再利用步骤S127中的公式得到等效标量表面阻抗曲线。
所述的步骤S3包括以下子步骤:
S31:构造全息张量表面阻抗调制公式,所述的全息张量表面阻抗调制公式满足张量阻抗矩阵Z是反Hermit纯虚数二阶矩阵:
式中,X是等效标量表面阻抗平均值,M是张量表面阻抗调制深度;Erad和Jsurf分别是阻抗表面上的目标场电场和表面电流表达式;上标H表示共轭转置,符号表示Kronecker乘积;
S32:计算源场的表面电流表达式和目标场的电场表达式;所述的目标场的电场表达式与期望的极化方式相关;其中,单极子激励起的表面电流表达式表示为:
式中,t代表平面内点(x,y)的矢量坐标,是该点到坐标原点的距离;
以右旋圆极化天线为例的目标场电场表达式表示为:
Erad=(jcosθ,1,-jsinθ)e-jk·t
其中,k=(sinθ,0,cosθ)是目标场在自由空间中传播的波矢量,θ是目标场主瓣往x方向偏离平面法向方向的角度。
S33:联合步骤S31和步骤S32中得到的三个公式,计算张量表面阻抗Z的各阻抗分量的全息调制公式;
式中γ=kx sinθ-ktr;
S34:根据各阻抗分量的公式,直接计算得到理论的等效标量表面阻抗曲线,并进而得到的等效标量表面阻抗最大值Zemax和主轴方向θemax
S35:结合步骤S2中得到的Zemax和g之间的数学关系,求得该点处单元晶格所需的贴片间距g和贴片开缝角度θs,完成天线模型的建立。
所述的单极子的工作频率与基于全息调制的人工张量阻抗表面天线一致,单极子的长度约为1/4工作波长;所述的单极子天线由SMA进行馈电激励。
所述的由SMA进行馈电激励由用等长度的SMA接头内导体探针直接替代。
所述的右旋圆极化天线实例由51×51个单元晶格构成,其物理尺寸是153mm×153mm,电尺寸6.1λ×6.1λ。
本发明的有益效果是:
(1)本发明目的在于利用全息张量表面阻抗调原理,在实现特定方向上高增益辐射的同时,获得特定计划方式(例如右旋圆极化)的目标辐射场;
(2)本发明提出了一种张量表面阻抗的等效标量表面阻抗拟合准则,有效地减少了张量表面阻抗调制建模所需的数据量,提高了天线设计的可操作性:
本发明还针对这种等效标量表面阻抗曲线进一步归纳总结,强调将等效标量表面阻抗曲线两个极轴作为主要的研究对象,而非直接关注张量阻抗矩阵Z各分量的具体数值,并针对表面阻抗值和贴片参数之间的映射提出了一种简化准则:对于开缝金属贴片的阻抗表面来说,在开缝大小固定的情况下,可以认为等效标量表面阻抗曲线的主轴方向θemax与开缝方向θs一致,且主轴所对应的等效标量表面阻抗最大值Zemax仅与单元晶格之间的贴片间距g有关;
尽管只匹配主轴信息θemax和Zemax并不能使实际单元晶格所对应的阻抗曲线与阻抗调制所需的理论阻抗曲线完全重合,但这种方式将晶格可变参量减少到了两个,大大降低了设计难度;同时θemax和Zemax也是等效标量表面阻抗曲线最主要的两个特征,在满足θemax和Zemax匹配的情况下,阻抗调制效应并不会受到太大影响。
附图说明
图1为HFSS仿真模型中的单元晶格;
图2为单元晶格的主从边界相位差设置;
图3为典型的等效标量表面阻抗曲线;
图4本发明的仿真S11;
图5为本发明的仿真增益;
图6为本发明方法流程图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案:基于全息调制的人工张量阻抗表面天线,它包括人工张量阻抗表面,所述的人工张量阻抗表面包括多个大小相等的单元晶格;所述的单元晶格从上至下依次包括金属地板、介质基体和金属贴片;通过改变金属贴片的几何参数控制人工张量阻抗表面的单元晶格所对应的表面阻抗,从而使天线表面上的表面阻抗依照源场和目标场的全息干涉图样规则分布;所述的源场是激励源照射在天线表面上的表面波,所述的目标场是期望得到的平面波辐射场;采用简单结构的天线(例如单极子天线)作为激励源。
如图6所示,基于全息调制的人工张量阻抗表面天线的实现方法,它包括以下步骤:
S1:标量阻抗提取仿真,包括以下子步骤:
S11:如图1所示,在Ansoft HFSS中使用单元晶格和周期边界条件仿真计算张量阻抗表面的等效标量表面阻抗,其中介质材料采用Rogers RO3006,厚度为1.27mm,晶格横向尺寸边长3mm;设计天线工作频率12GHz,空气腔高度约为1/2工作波长;
S12:对使用不同组合的主从边界的相位差,模拟不同的表面波传播方向θt,从而得到张量表面阻抗的等效标量表面阻抗Ze随表面波传播方向θt变化的等效标量表面阻抗曲线,包括以下子步骤:所述的主从边界的相位差表示从边界相对于主边界的相位延迟;
S121:如图2所示,在单元晶格模型上的四个壁上设置两对主从边界,分别为master x、master y、slave x和slave y,其主从边界相位差分别设置成φx和φy,可得到:
kx=φx/a,ky=φy/a;
式中,kx和ky分别表示表面波在x方向和y方向传播的波数,ky表示,a表示表面阻抗单元晶格的横向尺寸;
S122:将表面波在阻抗表面上传播的波矢量kt和与x轴的夹角θt进行表示:
S123:定义阻抗表面的等效折射率n,表示为:
式中,c=3×108m/s自由空间中的光速;ω是在该贴片几何参数下,阻抗表面的本征频率;k0是自由空间中的波数;vp和kt分别是表面波沿着阻抗表面传播的相速和波数;
S124:定义标量表面阻抗:
式中,η0=120π是自由空间中的波阻抗,αz是阻抗表面法向方向即z方向上表面波传播的衰减常数,j是虚数,j的平方为-1;
S125:选取合适的kt,通过φy和φx的不同组合,便能得到在该表面波传播方向上的本征频率ω,进而得到该方向上的标量表面阻抗值:
强调一下是这个方向上的标量表面阻抗值。因为这时候的标量表面阻抗会随着相位差的不同而变化,并不是一个固定的数值。
S126:张量表面阻抗Z与电场E和表面电流J的关系为:
TM模,其电场表达式可以写成:
S127:联合以上方程解得张量表面阻抗的归一化等效标量表面阻抗为:
所以,张量表面阻抗可由等效标量表面阻抗αz/k0描述,而等效标量表面阻抗与表面波传播方向θt有关,且计算提取方法与标量表面阻抗所采用的步骤S125中的公式相同;
S128:再得到三组不同的αz/k0之后,又可以将其代入步骤S127中的公式,反过来求得张量表面阻抗矩阵Z,进而再利用步骤S127中的公式得到等效标量表面阻抗曲线。
S2:等效标量表面阻抗曲线拟合:对不同贴片间距g和不同表面波传播方向θt的单元晶格模型进行仿真,得到一系列不同贴片间距g值的等效标量表面阻抗曲线,从而进一步提炼出Zemax和g之间的数学关系,为全息阻抗调制做好准备;其中,将等效标量表面阻抗曲线的两个极轴作为主要的研究对象,针对表面阻抗值和贴片参数g之间的映射提出了一种简化准则:对于开缝金属贴片的阻抗表面来说,在开缝大小固定的情况下,可以认为等效标量表面阻抗曲线的主轴方向θemax与贴片开缝角度θs一致即θemax=θs,且主轴所对应的等效标量表面阻抗最大值Zemax仅与单元晶格之间的贴片间距g有关;
如图3所示,图中数据点即为按照上述方法仿真求得的某一金属贴片开缝角度60°单元晶格的等效标量表面阻抗;实曲线为根据数据点计算得到的等效标量表面阻抗曲线;两条虚直线为等效标量表面阻抗的长短极轴。
在本发明所提出的简化规则下,对不同贴片大小的单元晶格进行仿真求得各自的等效标量表面阻抗曲线之后,可以进一步拟合得到等效标量表面阻抗最大值Zemax和贴片间距g之间的数学对应关系;以圆形金属贴片开缝宽度0.2mm的单元晶格模型为例,对应的拟合公式为:
Zemax=80.2218+121.8477g-1-27.2299g-2+2.1952g-3
而主轴的角度则与贴片开缝角度一致,即:
θemax=θs
S3:张量表面阻抗全息调制,包括以下子步骤:
S31:构造全息张量表面阻抗调制公式,所述的全息张量表面阻抗调制公式满足张量阻抗矩阵Z是反Hermit纯虚数二阶矩阵:
式中,X是等效标量表面阻抗平均值,M是张量表面阻抗调制深度;Erad和Jsurf分别是阻抗表面上的目标场电场和表面电流表达式;上标H表示共轭转置,符号表示Kronecker乘积;
S32:计算源场的表面电流表达式和目标场的电场表达式;所述的目标场的电场表达式与期望的极化方式相关;其中,单极子激励起的表面电流表达式表示为:
式中,t代表平面内点(x,y)的矢量坐标,是该点到坐标原点的距离;
以右旋圆极化天线为例的目标场电场表达式表示为:
Erad=(jcosθ,1,-jsinθ)e-jk·t
其中,k=(sinθ,0,cosθ)是目标场在自由空间中传播的波矢量,θ是目标场主瓣往x方向偏离平面法向方向的角度。
S33:联合步骤S31和步骤S32中得到的三个公式,计算张量表面阻抗Z的各阻抗分量的全息调制公式;
式中γ=kxsinθ-ktr;
S34:根据各阻抗分量的公式,直接计算得到理论的等效标量表面阻抗曲线,并进而得到的等效标量表面阻抗最大值Zemax和主轴方向θemax
S35:结合步骤S2中得到的Zemax和g之间的数学关系,求得该点处单元晶格所需的贴片间距g和贴片开缝角度θs,完成天线模型的建立。
S4:天线建模仿真:采用单极子作为整个天线的表面波激励源;在天线平面的坐标原点处打孔放置单极子天线;所述的单极子的工作频率与基于全息调制的人工张量阻抗表面天线一致,单极子的长度约为1/4工作波长;所述的单极子天线由SMA进行馈电激励。所述的单极子天线由SMA进行馈电激励由用等长度的SMA接头内导体探针直接替代。
所述的右旋圆极化天线实例由51×51个单元晶格构成,其物理尺寸是153mm×153mm,电尺寸6.1λ×6.1λ。
图4和图5分别是本发明的仿真得到的回波损耗S11和中心频率处的方向图。
S11表示的是天线的回波损耗,表征了天线的跟输入端口的匹配情况。如果S11很小的话,说明天线匹配良好;否则馈电端口的能量会被天线反射回来,造成性能的降低或者恶化。

Claims (4)

1.基于全息调制的人工张量阻抗表面天线,它包括人工张量阻抗表面和激励源,所述的人工张量阻抗表面包括多个大小相等的单元晶格;所述的单元晶格从上至下依次包括金属地板、介质基体和金属贴片;通过改变金属贴片的几何参数控制人工张量阻抗表面的单元晶格所对应的表面阻抗,从而使天线表面上的表面阻抗依照源场和目标场的全息干涉图样规则分布;所述的源场是激励源照射在天线表面上的表面波,所述的目标场是期望得到的平面波辐射场;其特征在于:
所述的基于全息调制的人工张量阻抗表面天线,在实现特定方向上高增益的同时,获得特定计划方式的目标辐射场;提出了一种张量表面阻抗的等效标量表面阻抗拟合准则,有效地减少了张量表面阻抗调制建模所需的数据量,提高了天线设计的可操作性;包括如下步骤:
S1:标量阻抗提取仿真,包括以下子步骤:
S11:在电磁仿真软件Ansoft HFSS中,使用单元晶格和周期边界条件模拟无限大均匀阻抗表面,所述的周期边界条件为主从边界;
S12:对使用不同组合的主从边界的相位差,模拟不同的表面波传播方向θt,从而得到张量表面阻抗的等效标量表面阻抗Ze随表面波传播方向θt变化的等效标量表面阻抗曲线;所述的主从边界的相位差表示从边界相对于主边界的相位延迟;
S2:等效标量表面阻抗曲线拟合:对不同贴片间距g和不同表面波传播方向θt的单元晶格模型进行仿真,得到一系列不同贴片间距g值的等效标量表面阻抗曲线,从而进一步提炼出Zemax和g之间的数学关系,为全息阻抗调制做好准备;其中,将等效标量表面阻抗曲线的两个极轴作为主要的研究对象,针对表面阻抗值和贴片参数g之间的映射提出了一种简化准则:对于开缝金属贴片的阻抗表面来说,在开缝大小固定的情况下,可以认为等效标量表面阻抗曲线的主轴方向θemax与贴片开缝角度θs一致即θemax=θs,且主轴所对应的等效标量表面阻抗最大值Zemax仅与单元晶格之间的贴片间距g有关;
S3:张量表面阻抗全息调制,包括以下子步骤:
S31:构造全息张量表面阻抗调制公式,所述的全息张量表面阻抗调制公式满足张量阻抗矩阵Z是反Hermit纯虚数二阶矩阵;
S32:计算源场的表面电流表达式和目标场的电场表达式;所述的目标场的电场表达式与期望的极化方式相关;
S33:联合步骤S31和步骤S32中得到的三个公式,计算张量表面阻抗Z的各阻抗分量的全息调制公式;
S34:根据各阻抗分量的公式,直接计算得到理论的等效标量表面阻抗曲线,并进而得到的等效标量表面阻抗最大值Zemax和主轴方向θemax
S35:结合步骤S2中得到的Zemax和g之间的数学关系,求得该点处单元晶格所需的贴片间距g和贴片开缝角度θs,完成天线模型的建立;
S4:天线建模仿真:采用单极子作为整个天线的表面波激励源;在天线平面的坐标原点处打孔放置单极子天线;
所述的步骤S12包括以下子步骤:
S121:在单元晶格模型上的四个壁上设置两对主从边界,其主从边界相位差分别设置成φx和φy,可得到:
kx=φx/a,ky=φy/a;
式中,kx和ky分别表示表面波在x方向和y方向传播的波数,a表示表面阻抗单元晶格的横向尺寸;
S122:将表面波在阻抗表面上传播的波矢量kt和与x轴的夹角θt进行表示:
S123:定义阻抗表面的等效折射率n,表示为:
<mrow> <mi>n</mi> <mo>=</mo> <mfrac> <mi>c</mi> <msub> <mi>v</mi> <mi>p</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <mi>&amp;omega;</mi> <mo>/</mo> <msub> <mi>v</mi> <mi>p</mi> </msub> </mrow> <mrow> <mi>&amp;omega;</mi> <mo>/</mo> <mi>c</mi> </mrow> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>k</mi> <mi>t</mi> </msub> <msub> <mi>k</mi> <mn>0</mn> </msub> </mfrac> <mo>;</mo> </mrow>
式中,c=3×108m/s自由空间中的光速;ω是在该贴片几何参数下,阻抗表面的本征频率;k0是自由空间中的波数;vp和kt分别是表面波沿着阻抗表面传播的相速和波数;
S124:定义标量表面阻抗:
<mrow> <mi>Z</mi> <mo>=</mo> <msub> <mi>j&amp;eta;</mi> <mn>0</mn> </msub> <mfrac> <msub> <mi>&amp;alpha;</mi> <mi>z</mi> </msub> <msub> <mi>k</mi> <mn>0</mn> </msub> </mfrac> <mo>=</mo> <msub> <mi>j&amp;eta;</mi> <mn>0</mn> </msub> <msqrt> <mfrac> <mrow> <msubsup> <mi>k</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mfrac> </msqrt> <mo>=</mo> <msub> <mi>j&amp;eta;</mi> <mn>0</mn> </msub> <msqrt> <mrow> <mfrac> <msubsup> <mi>k</mi> <mi>t</mi> <mn>2</mn> </msubsup> <msubsup> <mi>k</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mfrac> <mo>-</mo> <mn>1</mn> </mrow> </msqrt> <mo>;</mo> </mrow>
式中,η0=120π是自由空间中的波阻抗,αz是阻抗表面法向方向即z方向上表面波传播的衰减常数,j是虚数,j的平方为-1;
S125:选取合适的kt,通过φy和φx的不同组合,便能得到在该表面波传播方向上的本征频率ω,进而得到该方向上的标量表面阻抗值:
<mrow> <mi>Z</mi> <mo>=</mo> <msub> <mi>j&amp;eta;</mi> <mn>0</mn> </msub> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msqrt> <mrow> <msubsup> <mi>&amp;phi;</mi> <mi>x</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;phi;</mi> <mi>y</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mi>c</mi> </mrow> <mrow> <mi>a</mi> <mi>&amp;omega;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>1</mn> </mrow> </msqrt> <mo>;</mo> </mrow>
S126:张量表面阻抗Z与电场E和表面电流J的关系为:
<mrow> <mi>E</mi> <mo>=</mo> <mi>Z</mi> <mo>&amp;CenterDot;</mo> <mi>J</mi> <mo>&amp;DoubleRightArrow;</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>E</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>E</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>x</mi> </mrow> </msub> </mtd> <mtd> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <msub> <mi>J</mi> <mi>x</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>J</mi> <mi>y</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
TM模,其电场表达式可以写成:
<mrow> <msub> <mi>E</mi> <mrow> <mi>T</mi> <mi>M</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <msub> <mi>k</mi> <mn>0</mn> </msub> </mfrac> <mo>&amp;lsqb;</mo> <mo>-</mo> <msub> <mi>e</mi> <mi>z</mi> </msub> <msubsup> <mi>k</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>j&amp;alpha;</mi> <mi>z</mi> </msub> <msub> <mi>k</mi> <mi>t</mi> </msub> <mo>&amp;rsqb;</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jk</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>r</mi> </mrow> </msup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mi>z</mi> </msub> <mi>z</mi> </mrow> </msup> <mo>;</mo> </mrow>
S127:联合以上方程解得张量表面阻抗的归一化等效标量表面阻抗为:
<mrow> <mtable> <mtr> <mtd> <mrow> <mfrac> <msub> <mi>&amp;alpha;</mi> <mi>z</mi> </msub> <msub> <mi>k</mi> <mn>0</mn> </msub> </mfrac> <mo>=</mo> <mo>{</mo> <mo>-</mo> <mi>j</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;eta;</mi> <mn>0</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mrow> <mo>(</mo> <msubsup> <mi>&amp;eta;</mi> <mn>0</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <mo>-</mo> <msubsup> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>4</mn> <msubsup> <mi>&amp;eta;</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>-</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mi>sin</mi> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;times;</mo> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mi>sin</mi> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <msup> <mo>&amp;rsqb;</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> </msup> <mo>}</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>&amp;times;</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <mn>2</mn> <msub> <mi>&amp;eta;</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>-</mo> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mi>sin</mi> <mn>2</mn> <msub> <mi>&amp;theta;</mi> <mi>t</mi> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
S128:再得到三组不同的αz/k0之后,又可以将其代入步骤S127中的公式,反过来求得张量表面阻抗矩阵Z,进而再利用步骤S127中的公式得到等效标量表面阻抗曲线;
所述的步骤S3包括以下子步骤:
S31:构造全息张量表面阻抗调制公式,所述的全息张量表面阻抗调制公式满足张量阻抗矩阵Z是反Hermit纯虚数二阶矩阵:
<mrow> <mi>Z</mi> <mo>=</mo> <mi>j</mi> <mfenced open = "(" close = ")"> <mtable> <mtr> <mtd> <mi>X</mi> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mi>X</mi> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mi>j</mi> <mfrac> <mi>M</mi> <mn>2</mn> </mfrac> <mi>Im</mi> <mrow> <mo>(</mo> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>a</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;CircleTimes;</mo> <msubsup> <mi>J</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>f</mi> </mrow> <mi>H</mi> </msubsup> <mo>-</mo> <msub> <mi>J</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>f</mi> </mrow> </msub> <mo>&amp;CircleTimes;</mo> <msubsup> <mi>E</mi> <mrow> <mi>r</mi> <mi>a</mi> <mi>d</mi> </mrow> <mi>H</mi> </msubsup> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式中,X是等效标量表面阻抗平均值,M是张量表面阻抗调制深度;Erad和Jsurf分别是阻抗表面上的目标场电场和表面电流表达式;上标H表示共轭转置,符号表示Kronecker乘积;
S32:计算源场的表面电流表达式和目标场的电场表达式;所述的目标场的电场表达式与期望的极化方式相关;其中,单极子激励起的表面电流表达式表示为:
<mrow> <msub> <mi>J</mi> <mrow> <mi>s</mi> <mi>u</mi> <mi>r</mi> <mi>f</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mi>r</mi> </mfrac> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>jk</mi> <mi>t</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>t</mi> </mrow> </msup> <mo>;</mo> </mrow>
式中,t代表平面内点(x,y)的矢量坐标,是该点到坐标原点的距离;
以右旋圆极化天线为例的目标场电场表达式表示为:
Erad=(j cosθ,1,-j sinθ)e-jk·t
其中,k=(sinθ,0,cosθ)是目标场在自由空间中传播的波矢量,θ是目标场主瓣往x方向偏离平面法向方向的角度。
S33:联合步骤S31和步骤S32中得到的三个公式,计算张量表面阻抗Z的各阻抗分量的全息调制公式;
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>x</mi> </mrow> </msub> <mo>=</mo> <mi>X</mi> <mo>+</mo> <mi>M</mi> <mfrac> <mi>x</mi> <mi>r</mi> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;theta;</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;gamma;</mi> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>x</mi> <mi>y</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mi>M</mi> <mn>2</mn> </mfrac> <mo>(</mo> <mfrac> <mi>y</mi> <mi>r</mi> </mfrac> <mi>cos</mi> <mi>&amp;theta;</mi> <mi>cos</mi> <mi>&amp;gamma;</mi> <mo>-</mo> <mfrac> <mi>x</mi> <mi>r</mi> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>Z</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <mo>=</mo> <mi>X</mi> <mo>-</mo> <mi>M</mi> <mfrac> <mi>y</mi> <mi>r</mi> </mfrac> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> </mtd> </mtr> </mtable> </mfenced> <mo>;</mo> </mrow>
式中γ=kxsinθ-ktr;
S34:根据各阻抗分量的公式,直接计算得到理论的等效标量表面阻抗曲线,并进而得到的等效标量表面阻抗最大值Zemax和主轴方向θemax
S35:结合步骤S2中得到的Zemax和g之间的数学关系,求得该点处单元晶格所需的贴片间距g和贴片开缝角度θs,完成天线模型的建立。
2.根据权利要求1所述的基于全息调制的人工张量阻抗表面天线的实现方法,其特征在于:所述的单极子的工作频率与基于全息调制的人工张量阻抗表面天线一致,单极子的长度约为1/4工作波长;所述的单极子天线由SMA进行馈电激励。
3.根据权利要求2所述的基于全息调制的人工张量阻抗表面天线的实现方法,其特征在于:所述的由SMA进行馈电激励由用等长度的SMA接头内导体探针直接替代。
4.根据权利要求1所述的基于全息调制的人工张量阻抗表面天线的实现方法,其特征在于:所述的右旋圆极化天线实例由51×51个单元晶格构成,其物理尺寸是153mm×153mm,电尺寸6.1λ×6.1λ。
CN201510185161.7A 2015-04-17 2015-04-17 基于全息调制的人工张量阻抗表面天线及其实现方法 Active CN104733850B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510185161.7A CN104733850B (zh) 2015-04-17 2015-04-17 基于全息调制的人工张量阻抗表面天线及其实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510185161.7A CN104733850B (zh) 2015-04-17 2015-04-17 基于全息调制的人工张量阻抗表面天线及其实现方法

Publications (2)

Publication Number Publication Date
CN104733850A CN104733850A (zh) 2015-06-24
CN104733850B true CN104733850B (zh) 2017-10-10

Family

ID=53457488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510185161.7A Active CN104733850B (zh) 2015-04-17 2015-04-17 基于全息调制的人工张量阻抗表面天线及其实现方法

Country Status (1)

Country Link
CN (1) CN104733850B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106887691B (zh) * 2017-04-28 2019-04-09 电子科技大学 互补地全息调制表面双波束高增益天线
CN107394408A (zh) * 2017-07-24 2017-11-24 电子科技大学 基于全息张量阻抗表面的涡旋电磁波束天线及其实现方法
CN107368664A (zh) * 2017-08-01 2017-11-21 成都大学 一种高增益低交叉极化的全息阻抗调制表面天线设计方法及天线
CN108011186A (zh) * 2017-11-20 2018-05-08 电子科技大学 基于张量阻抗表面的方向图可重构大角度扫描相控阵天线
CN109546314A (zh) * 2018-10-22 2019-03-29 成都大学 高增益全息阻抗调制表面天线设计方法和天线
CN112736483B (zh) * 2020-12-25 2021-09-28 北京大学 一种极化可重构的二维波束扫描全息天线及其实现方法
CN113764894B (zh) * 2021-09-10 2022-10-18 西安电子科技大学 一种三波束独立极化的全息人工阻抗表面天线
CN113991318B (zh) * 2021-10-26 2023-01-20 西安电子科技大学 一种基于全息张量阻抗表面的共形表面波天线及其设计方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7911407B1 (en) * 2008-06-12 2011-03-22 Hrl Laboratories, Llc Method for designing artificial surface impedance structures characterized by an impedance tensor with complex components
US7929147B1 (en) * 2008-05-31 2011-04-19 Hrl Laboratories, Llc Method and system for determining an optimized artificial impedance surface
CN103367926A (zh) * 2013-07-11 2013-10-23 东南大学 一种基于全息阻抗表面的多波束天线设计方法
CN103367884A (zh) * 2013-07-29 2013-10-23 东南大学 一种基于全息原理和阻抗表面的低剖面锥形出射方向图天线
CN104112901A (zh) * 2014-07-18 2014-10-22 电子科技大学 全息人工阻抗表面共形天线

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7929147B1 (en) * 2008-05-31 2011-04-19 Hrl Laboratories, Llc Method and system for determining an optimized artificial impedance surface
US7911407B1 (en) * 2008-06-12 2011-03-22 Hrl Laboratories, Llc Method for designing artificial surface impedance structures characterized by an impedance tensor with complex components
CN103367926A (zh) * 2013-07-11 2013-10-23 东南大学 一种基于全息阻抗表面的多波束天线设计方法
CN103367884A (zh) * 2013-07-29 2013-10-23 东南大学 一种基于全息原理和阻抗表面的低剖面锥形出射方向图天线
CN104112901A (zh) * 2014-07-18 2014-10-22 电子科技大学 全息人工阻抗表面共形天线

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Scalar and Tensor Holographic Artificial Impedance Surfaces";Bryan H. Fong etc.;《IEEE Transactions on Antennas and Propagation》;20101031;第58卷(第10期);第3212页左栏第4段-右栏最后1段、第3213页左栏第2段-3219页右栏第2段,图2、3、7、8、10、11 *
"全息阻抗调制表面天线的分析与设计";张剑;《中国优秀硕士学位论文全文数据库(电子期刊)基础科学辑》;20130630;第2.2节第1段-第2.3.2节最后1段、第3.1节第1段-第3.3.2节最后1段,图3-2、3-4、3-6、3-9、3-18) *

Also Published As

Publication number Publication date
CN104733850A (zh) 2015-06-24

Similar Documents

Publication Publication Date Title
CN104733850B (zh) 基于全息调制的人工张量阻抗表面天线及其实现方法
Budhu et al. A novel and systematic approach to inhomogeneous dielectric lens design based on curved ray geometrical optics and particle swarm optimization
Jiang et al. Wideband transmit arrays based on anisotropic impedance surfaces for circularly polarized single-feed multibeam generation in the Q-band
CN104112901B (zh) 全息人工阻抗表面共形天线
CN103646151B (zh) 平面反射阵天线设计方法
Budhu et al. Design of planar and conformal, passive, lossless metasurfaces that beamform
CN104615836B (zh) 星载微带阵列天线热变形对电性能影响的快速预测方法
CN107394408A (zh) 基于全息张量阻抗表面的涡旋电磁波束天线及其实现方法
CN110011069B (zh) 基于全息超表面的共形表面波天线
CN103367926B (zh) 一种基于全息阻抗表面的多波束天线设计方法
CN203644950U (zh) 基于平板反射阵列的紧缩场天线
CN107404011A (zh) 一种基于特征模理论的低频机载多天线系统及其设计方法
CN113708077A (zh) 一种基于幅度和相位调制的波束赋形超表面结构及方法
CN109670140A (zh) 一种基于阵元特征模式的阵列天线电性能分析方法
CN110729821B (zh) 用于多目标无线能量传输的准无衍射波束形成方法
CN104916922A (zh) 一种面向无线能量传输系统的平面阵列天线布阵方法
Al-Sadoon et al. A New Beamforming Approach Using 60 GHz Antenna Arrays for Multi-Beams 5G Applications
Zhang et al. Electromagnetic-model-driven twin delayed deep deterministic policy gradient algorithm for stealthy conformal array antenna
CN103646126A (zh) 微带阵列聚焦天线的设计方法及微带阵列聚焦天线
Yao et al. Effect of wire space and weaving pattern on performance of microstrip antennas integrated in the three dimensional orthogonal woven composites
CN113991318B (zh) 一种基于全息张量阻抗表面的共形表面波天线及其设计方法
Li et al. Optimization design of skin antenna based on Bayesian optimization
Kim et al. High impedance holographic metasurfaces for conformal and high gain antenna applications
CN103928763A (zh) 基于变换光学的多波束天线
CN102680802A (zh) 一种紧缩场产生装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant