上海交通大学学报, 2022, 56(11): 1541-1546 doi: 10.16183/j.cnki.jsjtu.2021.221

电子信息与电气工程

源重构近远场变换的探头补偿算法

袁浩波, 董欣欣,, 张瑞雪, 鲁榄埔, 陈曦

西安电子科技大学 电子工程学院,西安 710071

A Probe Correction Technique for Near Field to Far Field Transform Based on Sources Reconstruction Method

YUAN Haobo, DONG Xinxin,, ZHANG Ruixue, LU Lanpu, CHEN Xi

School of Electronics Engineering, Xidian University, Xi’an 710071, China

通讯作者: 董欣欣,女,硕士生;E-mail:15091542878@163.com.

责任编辑: 王一凡

收稿日期: 2021-06-23  

基金资助: 国家自然科学基金(61501341)
陕西省自然科学基金(2021JM-129)

Received: 2021-06-23  

作者简介 About authors

袁浩波(1980-),男,湖北省天门市人,副教授,从事计算电磁学、电磁兼容方面的研究.

摘要

源重构近远场变换方法在天线测量、天线诊断以及成像等方面的应用越来越多.该方法将待测天线的近场分布等效为包围待测天线的一个平面上的磁流作用,从而进行近远场变换,但现有源重构方法由于没有合适的探头补偿算法而精度较低.对此,首先根据互易定理给出等效磁流和探头接收信号之间的积分方程;然后用矩量法将其转化为矩阵方程,并求解出磁流分布;最后得到待测天线的方向图.波导探头测量喇叭天线的例子表明:所提方法得到的方向图的均方根误差仅为现有无补偿方法方向图均方根误差的1/3,该探头补偿算法准确有效.

关键词: 近远场变换; 源重构方法; 探头补偿; 互易定理; 矩量法

Abstract

The sources reconstruction method (SRM) is widely applied to antenna measurement, antenna diagnostics and imaging, etc. It is a near field to far field transform, which utilizes near-field data to determine an equivalent magnetic current source over a fictitious flat surface enclosing the antenna under test. But the existing SRM is low in accuracy, since it does not have a proper probe correction algorithm. First, this paper introduces an integral equation connecting the equivalent magnetic current and the signal received by the probe according to the reciprocity theorem. Secondly, the integral equation is converted to a matrix equation by the method of moments and the magnetic current is determined by solving the matrix equation. Finally, the magnetic current leads to the pattern of the antenna under test readily. As an example, a horn is measured with a wave-guide probe and the receiving signals are processed by the proposed method to obtain the pattern, whose root mean square errors (RMSE) are only 1/3 of those of the pattern obtained by the SRM without probe correction. Therefore, the proposed probe correction algorithm is accurate and effective for the SRM.

Keywords: near field to far field transform; sources reconstruction method (SRM); probe correction; reciprocity theorem; method of moment

PDF (1218KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

袁浩波, 董欣欣, 张瑞雪, 鲁榄埔, 陈曦. 源重构近远场变换的探头补偿算法[J]. 上海交通大学学报, 2022, 56(11): 1541-1546 doi:10.16183/j.cnki.jsjtu.2021.221

YUAN Haobo, DONG Xinxin, ZHANG Ruixue, LU Lanpu, CHEN Xi. A Probe Correction Technique for Near Field to Far Field Transform Based on Sources Reconstruction Method[J]. Journal of Shanghai Jiaotong University, 2022, 56(11): 1541-1546 doi:10.16183/j.cnki.jsjtu.2021.221

平面近场测量是天线测量与诊断的最重要技术之一,主要分成两种.主流的技术采用平面波谱展开待测天线和探头天线的近场,然后采用快速傅里叶变换算出待测天线的波谱[1].该方法一般用于测量窄波束天线.由于测量平面不可能无限大,需要进行截断,所以平面近场测量的方向图中存在一个可信角[2],在可信角范围内的结果比较可靠,在此范围之外则不太可靠.即使采用某些特殊的技巧对截断误差进行弥补[3-4],也不能使得可信角范围外的方向图足够准确.此外,平面近场测量难以得到天线的后向方向图.

另一种平面近场测量方法称为源重构方法(SRM)[5-6].首先假定测量得到的接收信号即为待测天线的近区电场;然后采用等效磁流建立积分方程,并转化为矩阵方程进行求解;最后对磁流积分得到远场方向图.该方法在快速多极子、图形处理器(GPU)和机器学习等[7-9]最新技术的支持下,已经广泛应用于天线诊断、电磁成像和无相位测量等[10-16]各个方面.

SRM经常用于平面扫描的近远场变换,它最大的优点是对扫描面没有任何限制.比如该方法用于电磁场全息成像时[12],会将扫描面扩展为一个包围待测天线的长方体盒子,从而避免扫描面截断的问题.因此,SRM不需要考虑可信角,可以测量宽波束天线,还可以得到天线的后向方向图.

然而,SRM存在两个主要缺点.一是采用矩量法求解积分方程时速度较慢,这个问题目前可以通过快速多极子方法[8]或GPU计算[9]进行加速.另一个问题是,SRM还没有合适的探头补偿算法,这制约了算法的使用范围和发展前景.

最合理的SRM应该采用长方体包围面,并配合完善的探头补偿算法.作为初步的尝试,本文设计了一种基于互易定理的探头补偿算法,它适用于任意扫描面的情况.由于现有硬件设备不能进行长方体包围面扫描,所以只讨论平面扫描的情况.研究旨在摸索提高SRM精度的可能途径,从而进一步改善天线测量、天线诊断以及电磁成像等实际应用的性能.

1 互易定理

图1为平面测量系统的结构图.将待测天线记为A,探头天线记为B,探头在扫描平面S上进行测量.图中:r=(x, y, z)为场点坐标;r'=(x', y', z')为源点坐标;R=|r'-r|为源点与场点之间的距离.根据洛伦兹互易定理[17]可知探头天线接收的信号功率为

P(r)=∫s(EA×HB-EB×HA)·ds

式中:sS平面面积;EA为待测天线在S平面上的电场强度;HBEB分别为探头天线在S平面上的磁场强度和电场强度;HA为待测天线在S平面上的磁场强度;n^为探头平面单位法向量.

图1

图1   天线测量系统模型

Fig.1   Antenna measurement system model


根据等效原理,假定在天线A的口面附近建立一个无穷大平面S',然后将S'的左侧包围区域全部替换成理想磁导体,而右侧不变.那么S'上的等效磁流M(r')可以在S'右侧区域产生与原问题中EAHA完全相同的场.这两种场可以根据电场积分算子K和磁场积分算子L写成等效磁流的表达式:

EA(r)=-K(M)=-∫s'G1(R)(r'-rM(r')ds'
HA(r)= 1ηL(M)= 1jkηs'[2G1(R)M(r')+G2(R)(r'-r)×(r'-r)×M(r')]ds'
G1(R)= 1+jkR4πR3e-jkR
G2(R)= 3-k2R2+j3kR4πR5e-jkR

式中:s'为平面S'面积;ηk分别为自由空间波阻抗和波数;j=-1为复数虚部的符号.考虑到探头的辐射场主要分布在其口面附近一个较小区域SB内(见图1),将式(1)的积分区域S进行截断,并改写为

P(r)= sBHB×K(M)-EB×1ηL(M)· n^dsB

式中:sB为区域SB面积.上式右侧是个双重面积分,计算量较大.

2 矩量法求解积分方程

采用类似于文献[5]中的矩量法求解式 (6).首先仅考虑有限大区域上的磁流分布和有限大的测量面,即对理论上无限大的磁流面S'和测量面S进行截断.

其次,将S'剖分成I个均匀的小矩形面片,每个面片尺寸为D'x×D'y,并将S'上的磁流分布用脉冲基函数展开.i(1≤iI)个面片S'i上的磁流为M(r')=Mixx^+Miyy^,其含义是:只要处于面片S'i上,那么磁流为常矢量,其两个分量的数值大小MixMiy是待求的未知系数.

再次,将S剖分成N个均匀的小矩形面片,每个面片尺寸为Dx×Dy.同时,第n次测量的采样点刚好处于第n(1≤nN)个小矩形面片Sn的中心rn=(xn, yn, zn).探头沿着两个极化方向接收的信号分别记为PnxPny.相对应的,尽管假定探头的场分布始终不变,但是当探头位置移动到rn,将探头场的切向分量记为EBnx,HBnx,EBnyHBny.注意到探头场的法向分量对于式(6)右侧的积分值没有贡献,因此不予考虑.

最后,式(6)右侧的积分区域SB的几何中心位置rn由第n个采样点确定.SB往往覆盖S面上若干个矩形面片,其具体尺寸由探头口面场的分布确定.比如对于波导探头来说,其电磁场在口面之外快速衰减到几乎没有,因此SB一般取口面向外扩展1个波长即可.继续扩大该尺寸,则算法消耗的时间快速增加,但对于整个算法精度的提高没有明显的作用.

采用上述离散过程可将式(6)展开为

Pnx= i=1IMxisB(HBnxx^×Kni(x^)- EBnxx^× 1ηLni(x^))· n^dsB+MyisB(HBnxx^×Kni(y^)- EBnxx^× 1ηLni(y^))· n^dsB
Pny= i=1IMxisB(HBnyy^×Kni(x^)- EBnyy^× 1ηLni(x^))· n^dsB+MyisB(HBnyy^×Kni(y^)- EnyBy^× 1ηLni(y^))· n^dsB

写成矩阵方程:

PxPy2N=αβμυ2N×2NMxMy2I

其中矩阵元素为

α= sBHBnxx^×Kni(x^)-EBnxx^×1ηLni(x^)· n^dsB
β= sBHBnxx^×Kni(y^)-EBnxx^×1ηLni(y^)· n^dsB
μ= sBHBnyy^×Kni(x^)-EBnyy^×1ηLni(x^)· n^dsB
υ= sBHBnyy^×Kni(y^)-EBnyy^×1ηLni(y^)· n^dsB

求解式(9)得到磁流分布的展开系数,再通过式(2)即可算出待测天线的方向图.式(9)的系统矩阵是强奇异性的矩阵,求解时很不稳定.首先对系统矩阵进行奇异值分解[10],然后丢弃较小的奇异值以及对应的奇异矢量,最后解出磁流展开的系数.这种解法相当于通过滤波压制了数值误差以及测量数据中的噪声,因此得到的方向图比较光滑.

3 计算结果

图2所示的喇叭天线工作频率为8.15 GHz,对应自由空间的波长λ= 36.81 mm.喇叭的口面尺寸为122.5 mm×84.5 mm.探头天线为BJ70矩形开口波导,其口面大小为34.85 mm×15.80 mm.测量平面S大小为550.94 mm×512.94 mm,网格数设为 38×40,平面S距离喇叭天线口面150 mm,探头天线距离测量平面1 mm.此测量系统的可信角范围为[-55°,55°].等效磁流面S'大小为6λ×6λ,网格数为30×30,距离待测天线口面0.1λ.平面SS'中心均与喇叭天线口面中心对齐.探头天线口面附近进行互易定理积分的区域为SB.本文算法中系统矩阵的规模为 3 040×1 800,将探头补偿的平面波谱算法[1]得到的方向图作为对比的参考结果.

图2

图2   喇叭天线测量系统

Fig.2   Measurement configuration of a horn antenna


首先取互易定理的积分面积sB=9λ2,测试不同奇异值截断时的方向图.如图3图4所示,本文方法求解式(9)时奇异值σ的截断主要影响方向图左右两侧,即可信角的附近范围.图中:θ为待求场点的矢径与z'轴的夹角;Eunit为归一化电场强度的幅度,Eunit=20lg(E/Emax),其中E为电场强度,|E|max 为最大电场幅度.保留的奇异值越多,则方向图左右两侧的抖动越明显,反之,则方向图越光滑.这主要是因为采用脉冲基函数来表示等效磁流导致积分后填充的矩阵十分病态,对误差十分敏感,在保证矩阵特性的前提下保留较少的奇异值能够有效抑制噪声.从总体上看,奇异值的截断对方向图影响不大.图中丢弃σmax×10-2以下奇异值时,仅保留了349个奇异值,此时E面和H面方向图的均方根误差(RMSE)只有0.53 dB和 0.56 dB.

图3

图3   不同奇异值截断时的E面方向图

Fig.3   E-plane pattern for different truncation of the singular value of the system matrix


图4

图4   不同奇异值截断时的H面方向图

Fig.4   H-plane pattern for different truncation of the singular value of the system matrix


其次,求解矩阵方程时丢弃σmax×10-2以下奇异值,测试不同积分面积sB时的方向图.从图5图6可见,积分面积sB越大,则计算结果越接近参考结果.但积分面积增加至一定程度时,计算结果将趋于稳定.为了降低计算量,并保持较好的计算精度,sB一般取略大于波导探头口面的尺寸即可.

图5

图5   不同积分面积sB下的E面方向图

Fig.5   E-plane pattern for different integral areas sB


图6

图6   不同积分面积sB时的H面方向图

Fig.6   H-plane pattern for different integral areas sB


最后,在本文方法中取sB=9λ2,对σmax×10-2以下奇异值进行截断,其结果与无补偿的源重构算法[5]以及探头补偿的平面波谱算法对比.从图7图8可以看出无补偿的源重构算法在磁面上的主瓣明显偏离参考方向图,其RMSE为2.21 dB;在E面上的RMSE也高达1.55 dB.而本文方法在两个主面上的RMSE分别为0.53 dB和0.56 dB.因此本文的探头补偿算法更准确有效.在计算时间方面,无补偿算法平均计算时间为10 s,平面波谱算法的平均计算时间为0.5 s,有补偿算法由于需要在每次探头移动后进行互易定理积分,平均计算时长为 134 s.测量算法运行一次即可得出结果,但缩短计算时间并不是该算法研究的重点.

图7

图7   喇叭天线E面方向图

Fig.7   E-plane pattern of the horn computed by different methods


图8

图8   喇叭天线H面方向图

Fig.8   H-plane pattern of the horn computed by different methods


4 结语

本文提出一种对源重构近远场变换进行探头补偿的算法,通过喇叭天线的例子验证该方法的有效性.此外,该方法对扫描面和测量探头都没有限制.后续工作中,首先将把扫描面扩展为长方体盒子的形式,以便解决平面近场测量中对于扫描面的截断问题.其次将尝试采用多探头测量系统以提高测量效率,并通过本文方法进行补偿,以实现高精度的测量.

参考文献

RAHMAT-SAMII Y, WILLIAMS L I, YACCARINO R G.

The UCLA bi-polar planar-near-field antenna-measurement and diagnostics range

[J]. IEEE Antennas and Propagation Magazine, 1995, 37(6): 16-35.

[本文引用: 2]

GREGSON S, MCCORMICK J, PARINI C, et al. Principles of planar near-field antenna measurements[M]. Stevenage, London, UK: The Institution of Engineering and Technology, 2007.

[本文引用: 1]

BUCCI O M, MIGLIORE M D.

A new method for avoiding the truncation error in near-field antennas measurements

[J]. IEEE Transactions on Antennas and Propagation, 2006, 54(10): 2940-2952.

DOI:10.1109/TAP.2006.882151      URL     [本文引用: 1]

CANO-FÁCILA F J, PIVNENKO S, SIERRA-CASTAÑER M.

Reduction of truncation errors in planar, cylindrical, and partial spherical near-field antenna measurements

[J]. International Journal of Antennas and Propagation, 2012, 2012: 438727.

[本文引用: 1]

PETRE P, SARKAR T K.

Planar near-field to far-field transformation using an equivalent magnetic current approach

[J]. IEEE Transactions on Antennas and Propagation, 1992, 40(11): 1348-1356.

DOI:10.1109/8.202712      URL     [本文引用: 3]

KONNO K, CHEN Q.

A source reconstruction technique using eigenmode currents

[C]∥2019 International Symposium on Antennas and Propagation. Xi’an, China: IEEE, 2019: 1-3.

[本文引用: 1]

ALVAREZ Y, MARTINEZ-LORENZO J Á, LAS-HERAS F, et al.

An inverse fast multipole method for geometry reconstruction using scattered field information

[J]. IEEE Transactions on Antennas and Propagation, 2012, 60(7): 3351-3360.

DOI:10.1109/TAP.2012.2196950      URL     [本文引用: 1]

LOPEZ-FERNANDEZ J A, LOPEZ-PORTUGUES M, ALVAREZ LOPEZ Y, et al.

Fast antenna characterization using the sources reconstruction method on graphics processors

[J]. Progress in Electromagnetics Research, 2012, 126: 185-201.

DOI:10.2528/PIER11121408      URL     [本文引用: 2]

YAO H M, JIANG L J, SHA W E I.

Source reconstruction method based on machine learning algorithms

[C]∥2019 Joint International Symposium on Electromagnetic Compatibility, Sapporo and Asia-Pacific International Symposium on Electromagnetic Compatibility. Sapporo, Japan: IEEE, 2019: 774-777.

[本文引用: 2]

JØRGENSEN E, MEINCKE P, CAPPELLIN C, et al.

Improved source reconstruction technique for antenna diagnostics

[DB/OL].(2010-10-05)[2021-06-03]. https:∥www.researchgate.net/publication/228880223_Improved_source_reconstruction_technique_for_antenna_diagnostics#citations.

URL     [本文引用: 2]

ALVAREZ Y, LAS-HERAS F, GARCIA C.

The sources reconstruction method for antenna diagnostics and imaging applications

[M]∥Solutions and applications of scattering, propagation, radiation and emission of electromagnetic waves. London: IntechOpen, 2012.

[本文引用: 1]

ROSE C A.

Accuracy estimations of microwave holography from planar near-field measurements

[C] ∥AMTA Meeting and Symposium. Philadelphia, PA, USA: AMTA, 2000: 12-18.

[本文引用: 2]

BROWN T, JEFFREY I, MOJABI P.

Multiplicatively regularized source reconstruction method for phaseless planar near-field antenna measurements

[J]. IEEE Transactions on Antennas and Propagation, 2017, 65(4): 2020-2031.

DOI:10.1109/TAP.2017.2670518      URL     [本文引用: 1]

SHU Y F, WEI X C, YANG R, et al.

An iterative approach for EMI source reconstruction based on phaseless and single-plane near-field scanning

[J]. IEEE Transactions on Electromagnetic Compatibility, 2018, 60(4): 937-944.

DOI:10.1109/TEMC.2017.2756902      URL     [本文引用: 1]

KARAMEHMEDOVIĆ M, KIRKEBY A, KNUDSEN K.

Stable source reconstruction from a finite number of measurements in the multi-frequency inverse source problem

[J]. Inverse Problems, 2018, 34(6): 065004.

DOI:10.1088/1361-6420/aaba83      URL     [本文引用: 1]

XIE Z X, ZHANG Y H, HE S Y, et al.

A novel method for source reconstruction based on spherical wave expansion and NF-FF transformation using amplitude-only data

[J]. IEEE Transactions on Antennas and Propagation, 2019, 67(7): 4756-4767.

DOI:10.1109/TAP.2019.2908036      URL     [本文引用: 1]

RICHMOND J.

A reaction theorem and its application to antenna impedance calculations

[J]. IRE Transactions on Antennas and Propagation, 1961, 9(6): 515-520.

DOI:10.1109/TAP.1961.1145068      URL     [本文引用: 1]

/