源重构近远场变换的探头补偿算法
A Probe Correction Technique for Near Field to Far Field Transform Based on Sources Reconstruction Method
通讯作者: 董欣欣,女,硕士生;E-mail:15091542878@163.com.
责任编辑: 王一凡
收稿日期: 2021-06-23
基金资助: |
|
Received: 2021-06-23
作者简介 About authors
袁浩波(1980-),男,湖北省天门市人,副教授,从事计算电磁学、电磁兼容方面的研究.
源重构近远场变换方法在天线测量、天线诊断以及成像等方面的应用越来越多.该方法将待测天线的近场分布等效为包围待测天线的一个平面上的磁流作用,从而进行近远场变换,但现有源重构方法由于没有合适的探头补偿算法而精度较低.对此,首先根据互易定理给出等效磁流和探头接收信号之间的积分方程;然后用矩量法将其转化为矩阵方程,并求解出磁流分布;最后得到待测天线的方向图.波导探头测量喇叭天线的例子表明:所提方法得到的方向图的均方根误差仅为现有无补偿方法方向图均方根误差的1/3,该探头补偿算法准确有效.
关键词:
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:
本文引用格式
袁浩波, 董欣欣, 张瑞雪, 鲁榄埔, 陈曦.
YUAN Haobo, DONG Xinxin, ZHANG Ruixue, LU Lanpu, CHEN Xi.
SRM经常用于平面扫描的近远场变换,它最大的优点是对扫描面没有任何限制.比如该方法用于电磁场全息成像时[12],会将扫描面扩展为一个包围待测天线的长方体盒子,从而避免扫描面截断的问题.因此,SRM不需要考虑可信角,可以测量宽波束天线,还可以得到天线的后向方向图.
最合理的SRM应该采用长方体包围面,并配合完善的探头补偿算法.作为初步的尝试,本文设计了一种基于互易定理的探头补偿算法,它适用于任意扫描面的情况.由于现有硬件设备不能进行长方体包围面扫描,所以只讨论平面扫描的情况.研究旨在摸索提高SRM精度的可能途径,从而进一步改善天线测量、天线诊断以及电磁成像等实际应用的性能.
1 互易定理
式中:s为S平面面积;EA为待测天线在S平面上的电场强度;HB和EB分别为探头天线在S平面上的磁场强度和电场强度;HA为待测天线在S平面上的磁场强度;
图1
根据等效原理,假定在天线A的口面附近建立一个无穷大平面S',然后将S'的左侧包围区域全部替换成理想磁导体,而右侧不变.那么S'上的等效磁流M(r')可以在S'右侧区域产生与原问题中EA和HA完全相同的场.这两种场可以根据电场积分算子K和磁场积分算子L写成等效磁流的表达式:
式中:s'为平面S'面积;η和k分别为自由空间波阻抗和波数;j=
式中:sB为区域SB面积.上式右侧是个双重面积分,计算量较大.
2 矩量法求解积分方程
采用类似于文献[5]中的矩量法求解式 (6).首先仅考虑有限大区域上的磁流分布和有限大的测量面,即对理论上无限大的磁流面S'和测量面S进行截断.
其次,将S'剖分成I个均匀的小矩形面片,每个面片尺寸为D'x×D'y,并将S'上的磁流分布用脉冲基函数展开.第i(1≤i≤I)个面片S'i上的磁流为M(r')=Mix
再次,将S剖分成N个均匀的小矩形面片,每个面片尺寸为Dx×Dy.同时,第n次测量的采样点刚好处于第n(1≤n≤N)个小矩形面片Sn的中心rn=(xn, yn, zn).探头沿着两个极化方向接收的信号分别记为Pnx和Pny.相对应的,尽管假定探头的场分布始终不变,但是当探头位置移动到rn,将探头场的切向分量记为
最后,式(6)右侧的积分区域SB的几何中心位置rn由第n个采样点确定.SB往往覆盖S面上若干个矩形面片,其具体尺寸由探头口面场的分布确定.比如对于波导探头来说,其电磁场在口面之外快速衰减到几乎没有,因此SB一般取口面向外扩展1个波长即可.继续扩大该尺寸,则算法消耗的时间快速增加,但对于整个算法精度的提高没有明显的作用.
采用上述离散过程可将式(6)展开为
写成矩阵方程:
其中矩阵元素为
求解式(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λ.平面S和S'中心均与喇叭天线口面中心对齐.探头天线口面附近进行互易定理积分的区域为SB.本文算法中系统矩阵的规模为 3 040×1 800,将探头补偿的平面波谱算法[1]得到的方向图作为对比的参考结果.
图2
首先取互易定理的积分面积sB=9λ2,测试不同奇异值截断时的方向图.如图3和图4所示,本文方法求解式(9)时奇异值σ的截断主要影响方向图左右两侧,即可信角的附近范围.图中:θ为待求场点的矢径与z'轴的夹角;Eunit为归一化电场强度的幅度,Eunit=20lg(
图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
图5
图6
最后,在本文方法中取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
图8
4 结语
本文提出一种对源重构近远场变换进行探头补偿的算法,通过喇叭天线的例子验证该方法的有效性.此外,该方法对扫描面和测量探头都没有限制.后续工作中,首先将把扫描面扩展为长方体盒子的形式,以便解决平面近场测量中对于扫描面的截断问题.其次将尝试采用多探头测量系统以提高测量效率,并通过本文方法进行补偿,以实现高精度的测量.
参考文献
The UCLA bi-polar planar-near-field antenna-measurement and diagnostics range
[J].
A new method for avoiding the truncation error in near-field antennas measurements
[J].DOI:10.1109/TAP.2006.882151 URL [本文引用: 1]
Reduction of truncation errors in planar, cylindrical, and partial spherical near-field antenna measurements
[J].
Planar near-field to far-field transformation using an equivalent magnetic current approach
[J].DOI:10.1109/8.202712 URL [本文引用: 3]
A source reconstruction technique using eigenmode currents
[C]
An inverse fast multipole method for geometry reconstruction using scattered field information
[J].DOI:10.1109/TAP.2012.2196950 URL [本文引用: 1]
Fast antenna characterization using the sources reconstruction method on graphics processors
[J].DOI:10.2528/PIER11121408 URL [本文引用: 2]
Source reconstruction method based on machine learning algorithms
[C]
Improved source reconstruction technique for antenna diagnostics
[DB/OL].(
The sources reconstruction method for antenna diagnostics and imaging applications
[M]
Accuracy estimations of microwave holography from planar near-field measurements
[C]
Multiplicatively regularized source reconstruction method for phaseless planar near-field antenna measurements
[J].DOI:10.1109/TAP.2017.2670518 URL [本文引用: 1]
An iterative approach for EMI source reconstruction based on phaseless and single-plane near-field scanning
[J].DOI:10.1109/TEMC.2017.2756902 URL [本文引用: 1]
Stable source reconstruction from a finite number of measurements in the multi-frequency inverse source problem
[J].DOI:10.1088/1361-6420/aaba83 URL [本文引用: 1]
A novel method for source reconstruction based on spherical wave expansion and NF-FF transformation using amplitude-only data
[J].DOI:10.1109/TAP.2019.2908036 URL [本文引用: 1]
A reaction theorem and its application to antenna impedance calculations
[J].DOI:10.1109/TAP.1961.1145068 URL [本文引用: 1]
/
〈 |
|
〉 |
