一种网状反射面天线金属丝网等效电磁参数推演方法
【技术领域】
[0001] 本发明属于雷达技术领域，具体涉及雷达天线领域中的一种网状反射面天线金属 丝网等效电磁参数推演方法。
【背景技术】
[0002] 网状反射面天线广泛应用在通信、雷达、射电天文学、微波通信、卫星通信和跟踪 以及遥感等各个领域。金属丝网作为网状反射面天线的反射面以完成电波发射与接收任 务。天线电性能是评价网状反射面天线性能好坏的重要标准，而金属丝网作为反射材料对 天线电性能的影响至关重要。金属丝网在进行电性能分析过程中，需要建立导线或条带单 元；而导线或条带单元庞大的数目限制了电性能计算的高效进行。通过将金属丝网等效成 某实体介质，在获得该实体介质等效电磁参数的情况下分析实体介质的电性能，将有助于 天线电性能的高效计算。因此有必要对金属丝网进行等效电磁参数的研宄。
[0003] A. Miura 等人在文献 "Spaceborne mesh ref lector antennas with complex weaves: extended PO/periodic-MoM analysis"（IEEETransactions on Antennas and Propagation，Vol. 55, No. 4, pp. 1022-1029, 2007 年）中，公开了一种采用扩展的物理光学 法结合周期矩量法分析网状反射面天线金属丝网电性能的分析方法。采用A. Miura提出的 方法在分析金属丝网电性能时需要建立金属丝网的条带模型，条带模型的增加限制了天线 电性能的高效计算。将金属丝网等效为某一实体介质，通过实体介质等效电磁参数分析电 性能将在保证计算精度的前提下，有效减小计算时间。因此有必要开展针对网状反射面天 线的金属丝网等效电磁参数的研宄。
【发明内容】
[0004] 本发明的目的是克服上述现有技术的不足，提供一种网状反射面天线金属丝网等 效电磁参数推演方法。该方法通过求解非线性方程组获得网状反射面天线金属丝网等效电 磁参数。
[0005] 本发明的技术方案是：一种网状反射面天线金属丝网等效电磁参数推演方法，包 括如下步骤：
[0006] (1)输入金属丝网结构参数和电参数
[0007] 输入用户提供的网状反射面天线金属丝网结构参数和电参数信息，其中结构参数 包括丝网直径、丝网横向尺寸、丝网纵向尺寸，电参数包括工作波长、电磁波入射角、迭代收 敛精度；
[0008] (2)计算金属丝网电磁透射系数
[0009] 根据用户提供的结构参数，采用周期矩量法计算处于用户提供的电参数下的金属 丝网电磁透射系数，获得金属丝网在垂直极化电磁波入射时的反射系数Γ' ±与水平极化 电磁波入射时的反射系数//;
[0010] (3)计算实体介质的电磁透射系数
[0011] 假设存在某一实体介质，其复相对介电常数毛可表示为元其中ε ' ^ 表示实数项，ε " ^表示复数项；复相对磁导率兵可表示为众=A-J//,"，其中μ'表示实 数项，μ" ^表示复数项。通过下式计算实体介质的反射系数：
【主权项】
1. 一种网状反射面天线金属丝网等效电磁参数推演方法，其特征在于，包括如下步 骤： (1) 输入金属丝网结构参数和电参数 输入用户提供的网状反射面天线金属丝网结构参数和电参数信息，其中结构参数包括 丝网直径、丝网横向尺寸、丝网纵向尺寸，电参数包括工作波长、电磁波入射角、迭代收敛精 度； (2) 计算金属丝网电磁透射系数 根据用户提供的结构参数，采用周期矩量法计算处于用户提供的电参数下的金属丝网 电磁透射系数，获得金属丝网在垂直极化电磁波入射时的反射系数Γ' ±与水平极化电磁 波入射时的反射系数Γ ' //; (3) 计算实体介质的电磁透射系数 假设存在某一实体介质，其复相对介电常数元可表示为毛，其中ε' ^表示 实数项，ε "彦示复数项；复相对磁导率总可表示为其中μ'表示实数项， μ" t表示复数项。通过下式计算实体介质的反射系数：
其中，Γ±表示实体介质在垂直极化电磁波入射时的反射系数，Γ#表示实体介质在水 平极化电磁波入射时的反射系数，化表示电磁波入射角，足表示垂直极化入射时的该介质 的波阻抗，之,表示水平极化入射时的该介质的波阻抗，由下式获得tA :
其中，1^=2 31/^为自由空间波数，π为圆周率，λ为工作波长，γ'为相位常数， γ "为衰减常数，两者可由下式得到：
(4) 遵循金属丝网与实体介质的电磁透射系数相等的原则，建立关于电磁参数的非线 性方程组： 遵循金属丝网与实体介质分别在垂直极化电磁波入射与水平极化电磁波入射下的反 射系数对应项相等的原则，建立关于电磁参数的非线性方程组，得到下式 r 丄=r ,丄(7) r "= r ' " (8) 其中，Γ±表示实体介质在垂直极化电磁波入射时的反射系数，Γ' ±为金属丝网在垂 直极化电磁波入射时的反射系数，Γ //表示实体介质在水平极化电磁波入射时的反射系数， Γ ' //为金属丝网在水平极化电磁波入射时的反射系数； (5) 采用模函数极小法构造非线性方程组的模函数 采用模函数最小法构造非线性方程组的模函数，可以得到如下表达式： Φ = 0- 丄-Γ，丄）2+0-"-Γ， "）2 (9) 其中，Φ为模函数的函数值，用以评价非线性方程组的解的逼近程度，Γ±表示实体介 质在垂直极化电磁波入射时的反射系数，Γ ' ±为金属丝网在垂直极化电磁波入射时的反 射系数，Γ#表示实体介质在水平极化电磁波入射时的反射系数，Γ 为金属丝网在水平 极化电磁波入射时的反射系数。 (6) 采用模函数极小法求解实体介质的等效电磁参数。
2. 根据权利要求1所述的一种网状反射面天线金属丝网等效电磁参数推演方法，其特 征在于，步骤（2)所述的周期矩量法是一种求解金属丝网在垂直极化与水平极化电磁波入 射时的反射系数的一种计算方法，该方法的计算公式如下：
其中，€(/)表示反射电场矢量，ρ表示场点位置矢量，π表示圆周率，j表示虚数单位， ω表示圆频率，Utl表示真空磁导率，ε ^表示真空介电常数，RF'):表示处于源点p处的面 电流密度矢量，P表示源点位置矢量，表示源点与场点之间的周期格林函数，▽表 示微分算子，s表示积分面积。
3. 根据权利要求1所述的一种网状反射面天线金属丝网等效电磁参数推演方法，其特 征在于，步骤（5)和（6)中所述的模函数极小法是一种求解非线性方程组的迭代算法，计算 步骤如下： 1) 设迭代收敛精度为ε，迭代次数为k = 1 ; 2) 采用下式计算模函数的梯度列向量
其中，gk表示模函数相对于电磁参数的梯度列向量，/&表示实体介质的电磁透 射系数对电磁参数X的导数，5表示求偏导数运算，X表示电磁参数^ ^ ε " ^ ^ μ" y Γ±表示实体介质在垂直极化电磁波入射时的反射系数，Γ' ±为金属丝网在垂直 极化电磁波入射时的反射系数，Γ#表示实体介质在水平极化电磁波入射时的反射系数， Γ' #为金属丝网在水平极化电磁波入射时的反射系数； 3) 采用下式计算Jacobi矩阵
其中，Jk表示模函数对电磁参数的Jacobi矩阵，δ:表示求偏导数运算，Xp x2、x3、乂4分 别表示电磁参数ε ' p ε " p μ ' p μ " μ Γ±表示实体介质在垂直极化电磁波入射时 的反射系数，Γ#表示实体介质在水平极化电磁波入射时的反射系数； 4) 采用下式计算最优步长因子
其中，Xk表示最优步长因子，gk表示模函数相对于电磁参数的梯度列向量，Jk表示模 函数对电磁参数的Jacobi矩阵，上标T表示求转置运算； 5) 采用下式计算下一次迭代的电磁参数 V= xk-Akgk (14) 其中，Xk表示本次迭代的电磁参数，Xk+1表示下一次迭代的电磁参数，λk表示最优步长 因子，gk表示模函数相对于电磁参数的梯度列向量； 6) 判断非线性方程组迭代过程是否收敛，收敛准则通过下式计算得到 Iv-XkI I 彡 ε (15) 其中，Xk表示本次迭代的电磁参数，X !^表示下一次迭代的电磁参数，ε为迭代收敛精 度； 7) 如果不收敛，转至步骤2);否则，输出结果。
【专利摘要】本发明属于雷达技术领域，具体提供了一种网状反射面天线金属丝网等效电磁参数推演方法，其步骤包括：(1)输入金属丝网结构参数和电参数；(2)计算金属丝网电磁透射系数；(3)计算实体介质的电磁透射系数；(4)遵循金属丝网与实体介质的电磁透射系数相等的原则，建立关于电磁参数的非线性方程组；(5)采用模函数极小法构造非线性方程组的模函数；(6)采用模函数极小法求解实体介质的等效电磁参数。本发明通过求解非线性方程组获得网状反射面天线金属丝网等效电磁参数，弥补了现有方法无法计算金属丝网等效电磁参数的空白。
【IPC分类】G06F17-50
【公开号】CN104636551
【申请号】CN201510061879
【发明人】杨东武, 段宝岩, 张树新, 张逸群, 李申, 杜敬利
【申请人】西安电子科技大学
【公开日】2015年5月20日
【申请日】2015年2月5日