CN112798876B - 一种利用插值算法的天线近场测量方法及装置 - Google Patents
一种利用插值算法的天线近场测量方法及装置 Download PDFInfo
- Publication number
- CN112798876B CN112798876B CN202110299350.2A CN202110299350A CN112798876B CN 112798876 B CN112798876 B CN 112798876B CN 202110299350 A CN202110299350 A CN 202110299350A CN 112798876 B CN112798876 B CN 112798876B
- Authority
- CN
- China
- Prior art keywords
- field
- antenna
- interpolation
- electric field
- scanning surface
- 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
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 16
- 230000005684 electric field Effects 0.000 claims abstract description 106
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000010586 diagram Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 42
- 238000004364 calculation method Methods 0.000 claims description 31
- 239000000523 sample Substances 0.000 claims description 17
- 238000012360 testing method Methods 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 9
- 230000005672 electromagnetic field Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 15
- 230000008859 change Effects 0.000 abstract description 4
- 230000009466 transformation Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/10—Radiation diagrams of antennas
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0871—Complete apparatus or systems; circuits, e.g. receivers or amplifiers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/0864—Measuring electromagnetic field characteristics characterised by constructional or functional features
- G01R29/0892—Details related to signal analysis or treatment; presenting results, e.g. displays; measuring specific signal features other than field strength, e.g. polarisation, field modes, phase, envelope, maximum value
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本申请公开了一种利用插值算法的天线近场测量方法。步骤S10:采集并记录待测天线近场扫描面上的电场幅度、相位信息;所述近场扫描面为任意曲面形状,且需同时满足以下三个条件。步骤S20:通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值。步骤S30:通过双线性插值对待测天线近场扫描面上的电场相位进行插值。步骤S40:基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场。步骤S50:基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。上述方法中的近场扫描面不局限于平面、柱面和球面,可以根据实际测量环境变化变成任意曲面。
Description
技术领域
本申请涉及一种天线近场测量(near field measurement)方法,特别是一种利用插值(interpolation)算法的天线近场测量方法。
背景技术
天线测量方法主要包括直接法和间接法,直接法即远场测量方法,间接法即近场测量方法。天线近场测量一般是在微波暗室内进行,利用探头天线扫描获取待测天线近场区辐射场的数据,并通过近远场变换计算天线的远场方向图。
根据扫描面的不同类型,常见的天线近场测量包括平面、球面和柱面近场测量。随着共形天线(conformal antenna)越来越多的应用,很难构成平面、球面或柱面的规则扫描面(也称采样面),此时就出现了能够对任意形状天线进行近场测量的方法。申请公布号为CN106291130A、申请公布日为2017年1月4日的中国发明专利申请《一种任意曲面扫描的近场天线测量方法》就提出了一种适用于任意扫描面的近场天线测量方法。该方法采用双线性插值法对近场电场的实部分量与虚部分量进行插值处理,并根据惠更斯原理计算出远场电场。该方法仍存在一些不足。第一,测量***中直接获取到的近场电场分量通常是幅度、相位值,该方法需要通过一次数学变换将幅度、相位分量转换为实部、虚部分量,这增加了额外的计算量。第二,该方法仅在扫描面(也称采样面)为球形时插值误差较小,改变扫描面形状后电场虚部分量插值结果并不理想,影响后续计算的准确度。第三,若将该方法改为对近场电场的幅度、相位分量进行插值处理,由于电场的幅度值变化趋势相较于相位值十分平缓,通过双线性插值得到的曲线不是平滑的曲线,在数值陡变的位置插值结果与理论结果拟合得不好。
发明内容
本申请所要解决的技术问题是提供一种利用插值算法的天线近场测量方法,能够处理在任意曲面采样的近场数据,并且在仅采集少量数据的情况下仍可以计算出高质量远场方向图。
为解决上述技术问题,本申请提出了一种利用插值算法的天线近场测量方法,包括如下步骤。步骤S10:采集并记录待测天线近场扫描面上的电场幅度、相位信息;所述近场扫描面为任意曲面形状,且需同时满足以下三个条件。条件一:所述近场扫描面必须是位于电磁场没有置零的空间内。条件二:所述近场扫描面必须位于待测天线的辐射近场区,与待测天线的距离大于且小于2D2/λ;其中,D表示待测天线的最大尺寸,λ表示待测天线的中心工作频率的波长。条件三:所述近场扫描面的尺寸要大于截断的等效磁流面,以使所述近场扫描面边缘位置的截断电平在-40dB以下。步骤S20:通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值。步骤S30:通过双线性插值对待测天线近场扫描面上的电场相位进行插值。步骤S40:基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场。步骤S50:基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。上述方法中的近场扫描面不局限于平面、柱面和球面,可以根据实际测量环境变化变成任意曲面。
优选地,所述条件二中,所述近场扫描面与待测天线的距离在3至10个λ的范围内。这是一种优选的距离设置值。
进一步地,所述步骤S10中,将待测天线的位置固定,在近场扫描面上通过坐标定位扫描采样点,利用可移动的探头在待测天线前的近场扫描面上按设定的运动轨迹经过每一个扫描采样点,记录每一个扫描采样点位置所测量的电场的幅度、相位值,直至所有的扫描采样点的数据采集完毕。这是对步骤S10的详细说明。
进一步地,所述步骤S10中,所述近场扫描面为平面、弧面、半球面、半椭球面中的任一种。这是近场扫描面的一些优选形状。
优选地,所述步骤S10中,相邻扫描采样点的间隔小于或等于λ。这是一种优选的间隔设置值。
进一步地,所述步骤S10中,、 代表球坐标系下的两个球坐标方向;可移动探头以θ=0°为扫描的初始位置,在方向上每隔度采集一次扫描面上的电场,获得电场各方向的幅度、相位分量;然后探头在方向上转动Δθ度,重复上述步骤,直至θ=90°时数据采集完毕。这是步骤S10的一种具体实现方式,其中可移动探头在和方向上都是等角采样。
进一步地,所述步骤S20中,基于三次样条插值法的电场幅度插值是:三次样条插值函数S(x)在[x0, xn]内存在n+1个数据节点,数据节点就是指待测天线的近场扫描面上的扫描采样点;在每个分段区间[xi, xi+1]内,S(x)=Si(x)都是一个三次多项式,i=0,1,2,…,n-1。
其中,x代表参与插值计算的数据节点,xi代表参与插值计算的第i个数据节点,满足S(xi)=yi,yi代表第i个数据节点采样得到的电场幅度值;ai、bi、ci、di是第i个分段区间的三次样条插值函数中的待求系数。根据数据节点和边界条件计算出三次样条插值函数的二次微分值mi,进而计算出三次样条插值曲线的系数。
公式二:ai=yi 。
公式四:ci=mi/2。
其中,mi代表第i个数据节点的二次微分值,mi+1代表第i+1个数据节点的二次微分值。利用公式一至公式五完成水平方向上的电场幅度插值计算后,利用上述公式进行仰角方向上的电场幅度插值计算。这是步骤S20的一种具体实现方式。
进一步地,所述步骤S30中,基于双线性插值法的电场相位插值是:设任一插值点(x, y)被四个扫描采样点(x1, y1)、(x1, y2)、(x2, y1)、(x2, y2)所包围,这四个扫描采样点的电场相位值分别为、、、,那么该插值点的电场相位值如下。
进一步地,所述步骤S40中,计算远场电场是:基于惠更斯等效原理,利用一个包含待测天线口径面的无限大电璧上的等效磁流源代替待测天线,利用近场电场数据计算等效磁流源,再利用等效磁流源计算远场电场。
其中,代表在面上的二重积分运算,s’表示截断的等效磁流面,表示梯度算子,用来区分等效磁流面和其他位置的梯度算子,代表观察点坐标,观察点表示不在等效磁流面上的点;代表等效源点坐标,等效源点表示在等效磁流面上离散的等效磁流源所在的点;代表等效磁流源,代表自由空间格林函数。第一次利用上述公式七求解等效磁流源,此时代表待测天线的近场扫描面上经过插值后的扫描采样点的电场,为已知参数。第二次利用上述公式七求解待测天线的远场电场,此时等效磁流源为已知参数。这是步骤S40的一种具体实现方式。
本申请还提出了一种利用插值算法的天线近场测量装置,包括近场测试单元、幅度插值单元、相位插值单元、远场计算单元一和远场计算单元二。所述近场测试单元用于采集并记录待测天线近场扫描面上的电场幅度、相位信息;所述近场扫描面为任意曲面形状,且需同时满足以下三个条件。条件一:所述近场扫描面必须是位于电磁场没有置零的空间内。条件二:所述近场扫描面必须位于待测天线的辐射近场区,与待测天线的距离大于且小于2D2/λ;其中,D表示待测天线的最大尺寸,λ表示待测天线的中心工作频率的波长。条件三:所述近场扫描面的尺寸要大于截断的等效磁流面,以使所述近场扫描面边缘位置的截断电平在-40dB以下。所述幅度插值单元用于通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值。所述相位插值单元用于通过双线性插值对待测天线近场扫描面上的电场相位进行插值。所述远场计算单元一用于基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场。所述远场计算单元二用于基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。上述装置中的近场扫描面不局限于平面、柱面和球面,可以根据实际测量环境变化变成任意曲面。
本申请通过可移动探头对待测天线进行近场扫描,采集近场扫描面上电场的幅度、相位信息。所述近场扫描面不局限于平面、柱面和球面,可以根据实际测量环境变化变成任意曲面。本申请针对近场电场的幅度和相位分别选择合适的插值方法进行插值,利用三次样条插值法插值电场幅度,利用双线性插值法插值电场相位,能够减少实际扫描点的数量,降低对测量设备的要求,提升测量效率。本申请根据惠更斯等效原理,利用插值过的电场数据计算等效源进而计算出待测天线的远场电场,近远场变换简单方便,而且可以处理任意扫描面形状的电场数据。
附图说明
图1是本申请提出的利用插值算法的天线近场测量方法的流程图。
图2是待测天线与探头在笛卡尔坐标系与球坐标系的示意图。
图3是采用本申请所述方法获取的待测天线的远场方向图与理论结果的比较示意图。
图4是本申请提出的利用插值算法的天线近场测量装置的结构示意图。
图中附图标记说明:AUT为待测天线;10为近场测试单元;20为幅度插值单元;30为相位插值单元;40为远场计算单元一;50为远场计算单元二。
具体实施方式
请参阅图1,本申请提出的利用插值算法的天线近场测量方法包括如下步骤。
步骤S10:采集并记录待测天线近场扫描面上的电场幅度、相位信息。所述近场扫描面为任意曲面形状,且需同时满足以下三个条件。
条件一:所述近场扫描面必须是位于电磁场没有置零的空间内。
条件二:所述近场扫描面必须位于待测天线的辐射近场区,与待测天线的距离大于 且小于2D2/λ。其中,D表示待测天线的最大尺寸,λ表示待测天线的中心工作频率的波长。优选地,所述近场扫描面通常选在离待测天线3至10个波长的距离范围内,所述波长是指待测天线的中心工作频率的波长。
条件三:所述近场扫描面的尺寸要大于截断的等效磁流面(equivalent magneticcurrent plane),以使得所述近场扫描面边缘位置的截断电平在-40dB以下。
步骤S20:通过三次样条插值(cubic spline interpolation)对待测天线近场扫描面上的电场幅度进行插值,从而计算丰富电场幅度值。
步骤S30:通过双线性插值(bilinear interpolation)对待测天线近场扫描面上的电场相位进行插值,从而计算丰富电场相位值。
步骤S40:基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理(Huygens–Fresnel principle,也称惠更斯-菲涅耳原理)计算待测天线的远场电场。
步骤S50:基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。
请参阅图2,所述步骤S10例如是将待测天线AUT的位置固定,在近场扫描面上通过坐标定位扫描采样点(也称扫描点、采样点、测量点),利用可移动的探头在待测天线前任一曲面形状的近场扫描面(所述近场扫描面需同时满足前述三个条件)上按设定的运动轨迹经过每一个扫描采样点,记录每一个扫描采样点位置所测量的电场的幅度、相位值,直至所有的扫描采样点的数据采集完毕。所述近场扫描面优选为平面、弧面、半球面、半椭球面等。
作为一个示例,本申请还给出一种近场扫描面上的运动轨迹规划方案,采用球坐标。请参阅图2,、、代表球坐标系下的两个球坐标方向。可移动探头以θ=0°为扫描的初始位置,在方向上每隔度采集一次扫描面上的电场,获得电场各方向的幅度、相位分量;然后探头在方向上转动Δθ度,重复上述步骤,直至θ=90°时数据采集完毕。
所述步骤S20中,基于三次样条插值法的电场幅度插值例如是:三次样条插值函数S(x)在[x0, xn]内存在n+1个数据节点。数据节点就是指待测天线的近场扫描面上的扫描采样点,n+1个数据节点例如分别称为x0、x1、x2、……、xn。其中x0代表第1个数据节点,xn代表第n+1个数据节点。插值点是在每两个相邻的数据节点的空隙中进行插值。在每个分段区间[xi, xi+1]内,S(x)=Si(x)都是一个三次多项式,i=0,1,2,…,n-1。
其中,x代表参与插值计算的数据节点,xi代表参与插值计算的第i个数据节点,满足S(xi)=yi,yi代表第i个数据节点采样得到的电场幅度值。ai、bi、ci、di是第i个分段区间的三次样条插值函数中的待求系数。
根据数据节点和边界条件计算出三次样条插值函数的二次微分值mi,进而计算出三次样条插值曲线的系数,如下所示。
公式二:ai=yi 。
公式四:ci=mi/2。
其中,mi代表第i个数据节点的二次微分值,mi+1代表第i+1个数据节点的二次微分值。
利用公式一至公式五完成水平方向上的电场幅度插值计算后,利用上述公式进行仰角方向上的电场幅度插值计算。
所述步骤S30中,基于双线性插值法的电场相位插值例如是:设任一插值点(x, y)被四个扫描采样点(x1, y1)、(x1, y2)、(x2, y1)、(x2, y2)所包围,这四个扫描采样点的电场相位值分别为、、、,那么该插值点的电场相位值如下。
所述步骤S40中,计算远场电场例如是:基于惠更斯等效原理,利用一个包含待测天线口径面的无限大电璧(electric wall)上的等效磁流源(equivalent magneticcurrent source)代替待测天线,利用近场电场数据计算等效磁流源,再利用等效磁流源计算远场电场,计算公式如下。
其中,代表在面上的二重积分运算,s’表示截断的等效磁流面,表示梯度算子,用来区分等效磁流面和其他位置的梯度算子,代表观察点坐标,观察点表示不在等效磁流面上的点;代表等效源点坐标,等效源点表示在等效磁流面上离散的等效磁流源所在的点;代表等效磁流源,代表自由空间格林函数。
所述步骤S50已有现有技术能够实现,在此不作赘述。
在图2所示的实施例中,待测天线AUT为主辐射方向朝向z轴方向的1×25的半波振子天线阵列,这是一个平面阵列天线,天线单元沿y轴排布。近场扫描面形状为半个椭球(ellipsoid)形,该椭球面的x轴方向半径为5λ,y轴方向半径为10λ,z轴方向半径为5λ。可移动探头在和方向上等角采样(equal angle interval sampling,也称等角度间隔采样),采样间隔Δθ和都为4°。请参阅图3,采用本申请的利用插值算法的天线近场测量方法得的xoy面远场方向图在图3中以曲线表示,其中NF2FF表示近场变换(near field to farfield);理论结果在图3中以圆点表示,可以发现两组计算结果在θ为±30°角域内基本一致,证明了本申请的有效性。
在其他的实施例中,待测天线可以是阵列天线,也可以不是阵列天线。当待测天线能构成平面口径时,采用本申请所述方法的计算准确度较高。
请参阅图4,本申请提出的利用插值算法的天线近场测量装置包括近场测试单元10、幅度插值单元20、相位插值单元30、远场计算单元一40和远场计算单元二50。
所述近场测试单元10用于采集并记录待测天线近场扫描面上的电场幅度、相位信息。所述近场扫描面为任意曲面形状,且需同时满足以下三个条件。
条件一:所述近场扫描面必须是位于电磁场没有置零的空间内。
条件二:所述近场扫描面必须位于待测天线的辐射近场区,与待测天线的距离大于 且小于2D2/λ。其中,D表示待测天线的最大尺寸,λ表示待测天线的中心工作频率的波长。优选地,所述近场扫描面通常选在离待测天线3至10个波长的距离范围内,所述波长是指待测天线的中心工作频率的波长。
条件三:所述近场扫描面的尺寸要大于截断的等效磁流面(equivalent magneticcurrent plane),以使得所述近场扫描面边缘位置的截断电平在-40dB以下。
所述幅度插值单元20用于通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值。
所述相位插值单元30用于通过双线性插值对待测天线近场扫描面上的电场相位进行插值。
所述远场计算单元一40用于基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场。
所述远场计算单元二50用于基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。
与现有技术相比,本申请具有如下有益的技术效果。
第一,本申请是对近场电场的幅度和相位分量进行插值处理,不需要像CN106291130A将探头天线直接采集到的幅度、相位数据转换为实部、虚部后再进行插值。因此本申请与CN106291130A相比省略了一次数学变换计算,减少了计算量。
第二,本申请针对电场幅度、相位分量各自的变化特点,分别采用三次样条插值法和双线性插值法对近场电场的幅度、相位数据进行插值处理。而CN106291130A不考虑插值数据自身的特性,只使用一种插值方法。因此本申请可以得到更接近理论值的结果,也就是计算准确度更高。
第三,本申请可在增大采样间隔后仍能保证计算的准确度,相邻扫描采样点的间隔可以大于二分之一波长,但小于或等于一个波长的范围内。这是因为将待测天线的近场区辐射场的数据通过近远场变换计算天线的远场方向图的计算过程本质上是方程组的求解。增大采样间隔后方程数减少甚至变为欠定方程组,影响计算质量。而通过插值的方法增加参与计算的近场数据量,方程数增多使方程组变为超定方程组,因此仍能保证计算的准确度。
第四,本申请通过插值的方法对近场电场数据进行补全,所以可以减少实际扫描采样点的数量。减少扫描采样点的数量意味着采样间隔的增大,当采样间隔很小时,需要测量设备定位非常准确,而增大采样间隔可以降低对设备精度的要求。减少扫描采样点的数量以及增大采样间隔,也意味着采集记录次数减少,这使得测量时间减少,从而提升了测量效率。
以上仅为本申请的优选实施例,并不用于限定本申请。对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (10)
1.一种利用插值算法的天线近场测量方法,其特征是,包括如下步骤;
步骤S10:采集并记录待测天线近场扫描面上的电场幅度、相位信息;所述近场扫描面为任意曲面形状,且需同时满足以下三个条件;
条件一:所述近场扫描面必须是位于电磁场没有置零的空间内;
条件三:所述近场扫描面的尺寸要大于截断的等效磁流面,以使所述近场扫描面边缘位置的截断电平在-40dB以下;
步骤S20:通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值;
步骤S30:通过双线性插值对待测天线近场扫描面上的电场相位进行插值;
步骤S40:基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场;
步骤S50:基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。
2.根据权利要求1所述的利用插值算法的天线近场测量方法,其特征是,所述条件二中,所述近场扫描面与待测天线的距离在3至10个λ的范围内。
3.根据权利要求1所述的利用插值算法的天线近场测量方法,其特征是,所述步骤S10中,将待测天线的位置固定,在近场扫描面上通过坐标定位扫描采样点,利用可移动的探头在待测天线前的近场扫描面上按设定的运动轨迹经过每一个扫描采样点,记录每一个扫描采样点位置所测量的电场的幅度、相位值,直至所有的扫描采样点的数据采集完毕。
4.根据权利要求3所述的利用插值算法的天线近场测量方法,其特征是,所述步骤S10中,所述近场扫描面为平面、弧面、半球面、半椭球面中的任一种。
5.根据权利要求3所述的利用插值算法的天线近场测量方法,其特征是,所述步骤S10中,相邻扫描采样点的间隔小于或等于λ。
7.根据权利要求1所述的利用插值算法的天线近场测量方法,其特征是,所述步骤S20中,基于三次样条插值法的电场幅度插值是:三次样条插值函数S(x)在[x0, xn]内存在n+1个数据节点,数据节点就是指待测天线的近场扫描面上的扫描采样点;在每个分段区间[xi, xi+1]内,S(x)=Si(x)都是一个三次多项式,i=0,1,2,…,n-1;
其中,x代表参与插值计算的数据节点,xi代表参与插值计算的第i个数据节点,满足S(xi)=yi,yi代表第i个数据节点采样得到的电场幅度值;ai、bi、ci、di是第i个分段区间的三次样条插值函数中的待求系数;
根据数据节点和边界条件计算出三次样条插值函数的二次微分值mi,进而计算出三次样条插值曲线的系数;
公式二:ai=yi;
公式四:ci=mi/2;
其中,mi代表第i个数据节点的二次微分值,mi+1代表第i+1个数据节点的二次微分值;
利用公式一至公式五完成水平方向上的电场幅度插值计算后,利用上述公式进行仰角方向上的电场幅度插值计算。
9.根据权利要求1所述的利用插值算法的天线近场测量方法,其特征是,所述步骤S40中,计算远场电场是:基于惠更斯等效原理,利用一个包含待测天线口径面的无限大电璧上的等效磁流源代替待测天线,利用近场电场数据计算等效磁流源,再利用等效磁流源计算远场电场;
其中,代表在面上的二重积分运算,s’表示截断的等效磁流面,表示梯度算子,用来区分等效磁流面和其他位置的梯度算子,代表观察点坐标,观察点表示不在等效磁流面上的点;代表等效源点坐标,等效源点表示在等效磁流面上离散的等效磁流源所在的点;代表等效磁流源,代表自由空间格林函数;
10.一种利用插值算法的天线近场测量装置,其特征是,包括近场测试单元、幅度插值单元、相位插值单元、远场计算单元一和远场计算单元二;
所述近场测试单元用于采集并记录待测天线近场扫描面上的电场幅度、相位信息;所述近场扫描面为任意曲面形状,且需同时满足以下三个条件;
条件一:所述近场扫描面必须是位于电磁场没有置零的空间内;
条件三:所述近场扫描面的尺寸要大于截断的等效磁流面,以使所述近场扫描面边缘位置的截断电平在-40dB以下;
所述幅度插值单元用于通过三次样条插值对待测天线近场扫描面上的电场幅度进行插值;
所述相位插值单元用于通过双线性插值对待测天线近场扫描面上的电场相位进行插值;
所述远场计算单元一用于基于插值得到的待测天线近场电场幅度、相位值,根据惠更斯等效原理计算待测天线的远场电场;
所述远场计算单元二用于基于计算出的待测天线的远场电场,计算待测天线的远场归一化方向图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110299350.2A CN112798876B (zh) | 2021-03-22 | 2021-03-22 | 一种利用插值算法的天线近场测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110299350.2A CN112798876B (zh) | 2021-03-22 | 2021-03-22 | 一种利用插值算法的天线近场测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112798876A CN112798876A (zh) | 2021-05-14 |
CN112798876B true CN112798876B (zh) | 2021-06-22 |
Family
ID=75817304
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110299350.2A Active CN112798876B (zh) | 2021-03-22 | 2021-03-22 | 一种利用插值算法的天线近场测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112798876B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113533867B (zh) * | 2021-07-14 | 2022-09-06 | 西安电子科技大学 | 基于Fourier插值的远场方向图快速测量方法 |
CN113932942B (zh) * | 2021-09-08 | 2024-04-26 | 西安理工大学 | 一种超声波检测电缆内部温度场的方法及装置 |
CN113917278A (zh) * | 2021-09-08 | 2022-01-11 | 西安理工大学 | 基于声波测温的电缆故障定位方法 |
CN114675088B (zh) * | 2022-05-27 | 2022-08-23 | 浙江大学 | 一种无监督学习的辐射源快速近场扫描方法 |
CN115455881B (zh) * | 2022-11-11 | 2023-03-24 | 上海霍莱沃电子***技术股份有限公司 | 一种金属介质混合结构的全空间散射场计算方法及装置 |
CN117665414B (zh) * | 2024-01-31 | 2024-04-26 | 深圳大学 | 一种近场测量方法、装置、终端及可读存储介质 |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2262783C1 (ru) * | 2004-05-05 | 2005-10-20 | Томский Государственный Университет Систем Управления И Радиоэлектроники (Тусур) | Высокочастотная приемопередающая антенная система |
CN101483273A (zh) * | 2009-02-24 | 2009-07-15 | 中国航天科技集团公司第五研究院第五○四研究所 | 一种幅度和相位可变的阵列天线的校准方法 |
CN102195695A (zh) * | 2010-03-02 | 2011-09-21 | 电信科学技术研究院 | 天线校准的方法及装置 |
CN102540183A (zh) * | 2010-12-14 | 2012-07-04 | 中国科学院电子学研究所 | 基于柱面几何的三维微波成像方法 |
CN103176172A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种机载干涉sar基于同步内定标信号的相位测量补偿方法 |
CN103616569A (zh) * | 2013-11-20 | 2014-03-05 | 中国电子科技集团公司第四十一研究所 | 一种毫米波平面近场测试相位修正方法 |
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
CN104931811A (zh) * | 2015-05-05 | 2015-09-23 | 惠州Tcl移动通信有限公司 | 一种移动终端天线测试方法及装置 |
CN106291130A (zh) * | 2016-07-29 | 2017-01-04 | 昆山瀚德通信科技有限公司 | 一种任意曲面扫描的近场天线测量方法 |
CN106486769A (zh) * | 2015-09-01 | 2017-03-08 | 株式会社万都 | 用于线性相控阵天线的空间插值方法和设备 |
CN206040960U (zh) * | 2016-08-31 | 2017-03-22 | 重庆大学 | 电小平面惠更斯源天线 |
DE102015115989A1 (de) * | 2015-09-22 | 2017-03-23 | Karsten Menzel | Verfahren zur Ermittlung der Sicherheitsabstände in der Nachbarschaft von Mobilfunkantennen |
CN106771672A (zh) * | 2017-02-21 | 2017-05-31 | 北京邮电大学 | 基于三次样条插值法进行天线远场测量的方法及装置 |
CN108268674A (zh) * | 2016-12-30 | 2018-07-10 | 深圳光启高等理工研究院 | 天线罩的电性能评估方法及装置 |
CN109155678A (zh) * | 2016-03-07 | 2019-01-04 | 萨迪斯飞英国有限公司 | 用于天线阵列的校准技术 |
CN109444561A (zh) * | 2018-10-26 | 2019-03-08 | 成都德杉科技有限公司 | 一种用于阵列天线校准的天线面测量方法 |
CN111553046A (zh) * | 2020-03-13 | 2020-08-18 | 电子科技大学 | 一种基于球面波展开和源重构的天线辐射计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8212716B2 (en) * | 2007-12-31 | 2012-07-03 | Elta Systems Ltd. | System and method for calibration of phased array antenna having integral calibration network in presence of an interfering body |
TWI603099B (zh) * | 2016-02-26 | 2017-10-21 | Near Field Antenna Measurement and Measurement System |
-
2021
- 2021-03-22 CN CN202110299350.2A patent/CN112798876B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2262783C1 (ru) * | 2004-05-05 | 2005-10-20 | Томский Государственный Университет Систем Управления И Радиоэлектроники (Тусур) | Высокочастотная приемопередающая антенная система |
CN101483273A (zh) * | 2009-02-24 | 2009-07-15 | 中国航天科技集团公司第五研究院第五○四研究所 | 一种幅度和相位可变的阵列天线的校准方法 |
CN102195695A (zh) * | 2010-03-02 | 2011-09-21 | 电信科学技术研究院 | 天线校准的方法及装置 |
CN102540183A (zh) * | 2010-12-14 | 2012-07-04 | 中国科学院电子学研究所 | 基于柱面几何的三维微波成像方法 |
CN103176172A (zh) * | 2013-02-06 | 2013-06-26 | 中国科学院电子学研究所 | 一种机载干涉sar基于同步内定标信号的相位测量补偿方法 |
CN103616569A (zh) * | 2013-11-20 | 2014-03-05 | 中国电子科技集团公司第四十一研究所 | 一种毫米波平面近场测试相位修正方法 |
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
CN104931811A (zh) * | 2015-05-05 | 2015-09-23 | 惠州Tcl移动通信有限公司 | 一种移动终端天线测试方法及装置 |
CN106486769A (zh) * | 2015-09-01 | 2017-03-08 | 株式会社万都 | 用于线性相控阵天线的空间插值方法和设备 |
DE102015115989A1 (de) * | 2015-09-22 | 2017-03-23 | Karsten Menzel | Verfahren zur Ermittlung der Sicherheitsabstände in der Nachbarschaft von Mobilfunkantennen |
CN109155678A (zh) * | 2016-03-07 | 2019-01-04 | 萨迪斯飞英国有限公司 | 用于天线阵列的校准技术 |
CN106291130A (zh) * | 2016-07-29 | 2017-01-04 | 昆山瀚德通信科技有限公司 | 一种任意曲面扫描的近场天线测量方法 |
CN206040960U (zh) * | 2016-08-31 | 2017-03-22 | 重庆大学 | 电小平面惠更斯源天线 |
CN108268674A (zh) * | 2016-12-30 | 2018-07-10 | 深圳光启高等理工研究院 | 天线罩的电性能评估方法及装置 |
CN106771672A (zh) * | 2017-02-21 | 2017-05-31 | 北京邮电大学 | 基于三次样条插值法进行天线远场测量的方法及装置 |
CN109444561A (zh) * | 2018-10-26 | 2019-03-08 | 成都德杉科技有限公司 | 一种用于阵列天线校准的天线面测量方法 |
CN111553046A (zh) * | 2020-03-13 | 2020-08-18 | 电子科技大学 | 一种基于球面波展开和源重构的天线辐射计算方法 |
Non-Patent Citations (5)
Title |
---|
A far-field transformation from phaseless near-field data on helicoidal geometry;S. Costanzo等;《IEEE Antennas and Propagation Society Symposium》;20040625;全文 * |
Prediction of far-field from under-sampled near-field data based on cubic spline interpolation in spherical near-field antenna measurements;Meng Li等;《 2017 IEEE 5th International Symposium on Electromagnetic Compatibility (EMC-Beijing)》;20171031;全文 * |
基于RPIM法的圆极化天线特性参数快速测试;张晓燕等;《微波学报》;20190831;第36卷(第4期);全文 * |
插值方法在天线平面近场测试中的应用;刘灵鸽;《空间电子技术》;20080325(第1期);全文 * |
球面多探头天线近场测试***校准方法研究;胡楚锋等;《仪器仪表学报》;20170515;第38卷(第5期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112798876A (zh) | 2021-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112798876B (zh) | 一种利用插值算法的天线近场测量方法及装置 | |
CN106291130B (zh) | 一种任意曲面扫描的近场天线测量方法 | |
JP5669337B2 (ja) | 比吸収率を測定するためのシステム及び方法 | |
CN106443587B (zh) | 一种高分辨率的快速反卷积声源成像算法 | |
US9335359B2 (en) | Far electromagnetic field estimation method and apparatus, and near electromagnetic field measurement apparatus | |
CN111366793B (zh) | 一种减小截断误差的平面天线近场测量方法及*** | |
JPH07294456A (ja) | X線コンピュータ断層写真装置 | |
CN110335301A (zh) | 一种基于激光雷达和结构光的三维扫描方法及装置 | |
CN113533867B (zh) | 基于Fourier插值的远场方向图快速测量方法 | |
CN107783086B (zh) | 用于诊断天线阵口径幅相场的畸变位置的方法 | |
US20090039285A1 (en) | Method and device for controlling and monitoring a position of a holding element | |
CN114966239B (zh) | 基于激励系数变量可分离的准远场测量方法 | |
Hatzis et al. | Implementation of a combination planar and spherical near-field antenna measurement system using an industrial 6-axis robot | |
CN116520035A (zh) | 一种采用平面近场快速反演方法的二维方向图测试方法 | |
CN115524541B (zh) | 基于准远场模式滤波的任意极化线天线阵方向图测量方法 | |
CN112710903A (zh) | 一种基于多探头球面场天线测试***的相位中心测试方法 | |
CN107239602B (zh) | 一种基于曲线拟合的探头天线模型快速计算方法 | |
EP4199034A1 (en) | Microscopy feedback for improved milling accuracy | |
CN115542024A (zh) | 一种天线方向图近场测量方法 | |
CN109374989B (zh) | 一种分时数据合成的共形阵列测试方法 | |
Kaslis et al. | Revising the Theory and Practice of Electrical Alignment Procedures for Spherical Antenna Measurement Facilities | |
TW202202858A (zh) | 準遠場量測系統、準遠場量測方法 | |
EP2333565B1 (en) | Three-dimensional base setting method for image data | |
CN111965602A (zh) | 一种相控阵雷达幅相一致性检测方法和*** | |
Yaccarino et al. | A comparison of conventional and phaseless planar near-field antenna measurements: the effect of probe position errors |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |