CN106338655B - 一种平面近场探头安装精度引起的测试误差修正方法 - Google Patents
一种平面近场探头安装精度引起的测试误差修正方法 Download PDFInfo
- Publication number
- CN106338655B CN106338655B CN201610710111.0A CN201610710111A CN106338655B CN 106338655 B CN106338655 B CN 106338655B CN 201610710111 A CN201610710111 A CN 201610710111A CN 106338655 B CN106338655 B CN 106338655B
- Authority
- CN
- China
- Prior art keywords
- probe
- data
- field
- fft
- installation accuracy
- 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
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
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
一种平面近场探头安装精度引起的测试误差修正方法,包括步骤如下:一、按照平面近场测试方法将待测天线进行安装,场地校准;二、在探头初始位置,设定扫描范围、测试距离,进行第一个极化分量数据的采集,得到近场数据D0,探头旋转90度得到近场数据D90,探头继续旋转90度得到近场数据D180,探头继续旋转90度得到近场数据D270;三、进行近远场变换与数据修正,得到最终修正后的待测天线方向图数据。本发明避免了探头安装的机械调整工作,通过数学算法对探头安装精度引起的测量误差进行修正,得到修正探头安装精度后的方向图测试数据,极大的降低平面近场对于探头安装的机械精度要求。
Description
技术领域
本发明涉及一种平面近场误差修正方法,特别是一种平面近场探头安装精度引起的测试误差修正方法,属于天线测量技术领域。
背景技术
平面近场是天线测试方法中的一种。近场测量是在被测目标的近场测量电磁场的一种技术。利用已知辐射特性的天线作为探头在待测目标近场区的一个表面上作扫描。再通过数学变换,得到待测目标的远场特性。平面近场即扫描区域为一个平面。平面近场起源于上世纪50年代,其优良的特性很快受到了广泛的好评,近几年更是得到了很大的推广。
随着近场测试技术的广泛使用,2012年美国IEEE标准局发布了关于近场测试基本误差源,总共划分为18项。其中探头安装误差就是18项误差中作为探头部分相关重要误差源中的一项。
该部分的误差主要由两部分引起,第一部分由于探头自身直线度以及机械组装精度导致的探头倾斜,另一部分由于探头安装误差无法保证探头口面中心与安装轴的延伸线不相交,导致偏心误差。随着天线测试频率越来越高,该误差相对于波长就会越来对大,对于天线测试影响就越来越大,会成为平面近场高频天线测试的障碍。
该部分造成的误差现阶段均采用机械方式进行调整。探头垂直度一般在探头加工完成后进行垂直度标定,由于探头为金属矩形腔体结构,对于矩形内表面平整度也有一定要求,加工完成后很难进行高精度的垂直度调整,只能在探头安装时根据倾斜方向使用非常薄的垫片进行整体倾斜度调整,调整难度很大。对于探头安装轴与探头口面中心不相交导致的偏心误差调整难度更大,为了保证探头在运动时的稳定度,探头安装在安装轴上时通道都采用销钉、螺母等硬连接方式。探头安装时调整此误差的可动范围极小,调整起来非常困难。
发明内容
本发明的技术解决问题是:克服现有技术的不足,本发明提供了一种平面近场探头安装精度引起的测试误差修正方法,避免了探头安装的机械调整工作,采用数据后处理的方式,通过数学算法对探头安装精度引起的测量误差进行修正,得到修正探头安装精度后的方向图测试数据,降低了平面近场对于探头安装的机械精度要求。
本发明的技术解决方案是:一种平面近场探头安装精度引起的测试误差修正方法,包括步骤如下:
步骤一、按照平面近场测试方法将待测天线进行安装,场地校准;
步骤二、在探头初始位置,设定扫描范围、测试距离,进行第一个极化分量数据的采集,得到近场数据D0;探头旋转90度,进行第二极化分量数据的测量,得到近场数据D90;探头继续旋转90度,进行第三极化分量数据的测量,得到近场数据D180;探头继续旋转90度,进行第四极化分量数据的测量,得到近场数据D270;
步骤三、进行近远场变换与数据修正,得到最终修正后的待测天线方向图数据,计算公式如下:
其中,am_h=(am_fft(D0)+am_fft(D180))/2,
ph_h=(ph_fft(D0)+ph_fft(D180))/2-π/2,
am_v=(am_fft(D90)+am_fft(D270))/2,
ph_v=(ph_fft(D90)+ph_fft(D270))/2-π/2,
am_h表示修正后的水平极化的幅度,ph_h表示修正后的水平极化的相位,am_v表示修正后的垂直极化的幅度,ph_v表示修正后的垂直极化的相位;
am_fft表示傅里叶变换后取幅度部分,ph_fft表示傅里叶变换后取相位部分;e为自然底数;为待测天线方向图,为探头0度位置方向图,为探头90度位置方向图,i为复数的虚部单位量,i2=-1。
所述步骤一中扫描范围采用截断电平小于-35dB。
所述步骤一中测试距离为距离测试天线口面3~10个波长。
本发明与现有技术相比的优点在于:
(1)本发明采用数据后处理的方式通过数学算法对探头安装精度引起的测量误差进行修正,得到修正探头安装精度后的方向图测试数据,对探头安装精度引起的测量误差进行修正,得到修正探头安装精度后的方向图测试数据,实现难度远远低于探头安装的机械调整。
(2)本发明的修正算法在工程应用中简单方便,仅仅利用平面近场测试系统本身即可完成修正所需数据的采集工作,无需借助额外的测量工具与测量手段;规避了传统探头安装精度通过机械调整的难点,如操作难度大,调整精度要求高,调整时间长等。
(3)本发明采用数据后处理修正的方式,精度远高于机械调整精度,机械校准精度如何调整都存在一定的偏差,数据后处理修正方面在理论上经过了严密的数学推导,对探头安装精度所引起的误差进行了修正。
附图说明
图1探头安装误差模拟示意图。
图2探头旋转四个方向运动示意图。
具体实施方式
如图1、图2所示,建立坐标系OXYZ,Z轴方向为探头的旋转轴,Z轴正方向为探头开口向外方向,一种平面近场探头安装精度引起的测试误差修正方法,包括步骤如下:
(1)按照常规平面近场测试方法将待测天线进行安装,场地校准;
(2)在探头初始位置,选择合适扫描范围,一般采用截断电平小于-35dB的范围,以及合适的测试距离,一般选择距离测试天线口面3~10个波长,进行第一个极化分量数据的采集,得到近场数据D0;此时探头位于0度位置;
(3)探头按照顺时针方向(沿Z轴正方向看)旋转90度,进行第二极化分量数据的测量,得到近场数据D90;此时探头位于90度位置;
(4)探头继续旋转90度,进行第三极化分量数据的测量,得到近场数据D180;此时探头位于180度位置;
(5)探头继续旋转90度,进行第四极化分量数据的测量,得到近场数据D270;此时探头位于270度位置;
(6)利用修正算法结合(2)到(5)所得数据,进行近远场变换与数据修正。得到最终修正后的方向图数据。
二、算法简介
探头安装误差或者探头自身倾斜会引起探头两次扫描数据网格原点不重合以及探头倾斜所导致的探头补偿方向图变化。这两种情况基本涵盖了由于探头安装误差对于圆极化测试所带来的误差源。如图1所示,建立坐标系OXYZ,Z轴方向为探头的旋转轴,探头既有位置偏心,又有微小倾斜角度δ。这两种情况是平时普遍存在的,由于其量非常小,通常在测试时被忽略。
1探头旋转运动轨迹方程
如图2所示,探头由于偏心和倾斜,探头口面中心运动轨迹为圆C,偏移量L1为探头根部与旋转轴偏心所致,L2为探头自身倾斜导致。此假设式具有一般性。
运动轨迹曲线可以用下式表示:
其中,x为探头口面中心横坐标,y为探头口面中心纵坐标,r为运动半径即圆C的半径,θ0为探头口面中心初始位置(x0,y0)决定的初始角即初始位置(x0,y0)与X轴的夹角。
则探头在步骤(3)、(4)、(5)、(6)测试时四个位置的口面中心坐标可以表示为:
x90=rcos(θ0+90)=-rsinθ0
y90=rsin(θ0+90)=rcosθ0 (式4)
x180=rcos(θ0+180)=-rcosθ0
y180=rsin(θ0+180)=-rsinθ0 (式5)
x270=rcos(θ0+270)=rsinθ0
y270=rsin(θ0+270)=-rcosθ0 (式6)
其中,探头按照顺时针方向旋转90度的口面中心坐标位置为(x90,y90),探头继续旋转90度的口面中心坐标位置为(x180,y180),探头继续旋转90度的口面中心坐标位置为(x270,y270)。
可以看出相隔180度探头位置总是关于轨迹圆心中心对称。
2实际探头方向图计算
探头在初始位置的倾斜角度可以利用测量工具,如关节臂,垂直千分尺,激光跟踪仪等测量出探头的倾斜角度,此处不再赘述测量方法与测量过程。不失一般性的假设倾斜角度在AzoverEl坐标系下为(Az,El),因为在OXYZ坐标系中,探头沿着Z轴旋转,所以与Z轴夹角θ不变,在XOY面投影与X轴夹角线性增加90度。通过式7,求得旋转90度后的AzoverEl坐标系下的倾斜角度(Az90,El90):
求得El90=arcsin(cosElsinAz),同理可以求出其余两个位置时探头的倾斜角度。
El180=arcsin(-sinAz),Az180=arcsin(-cosElsinAz/cos(arcsin(-sinAz)));
El270=arcsin(-cosElsinAz),
再根据方向图公式,分别求出方向图旋转(Az,El)、(Az90,El90)、(Az180,El180)、(Az270,El270)的方向图数据。
3数据修正算法
以图2所示的坐标原点为方向图输出原点,根据平面近场近远场变换理论,步骤(3)、(4)、(5)、(6)所采集的数据与天线方向图关系如式8所示。
其中,为待测天线方向图,因为是同一天线所以远场方向图都为 分别为探头在四个位置时该坐标系下的方向图,其方向图是通过“2实际探头方向图计算”部分计算所得;D0、D90、D180、D270分别为探头在0度、90度、180度、270度极化采集所得到的近场数据,fft表示近远场变化中的傅里叶变换。为相位补偿因子,以第一个位置为例,实际网格中心为(x0,y0),在图2中所示的坐标系输出时相位需要原点平移。e为自然底数,k为波数,u为坐标系中单位矢量与X轴之间夹角的余弦,v分别为坐标系中单位矢量与Y轴之间夹角的余弦。i为复数的虚部单位量,i2=-1。
联立式8中的所有方程以及上面已经求出的探头方向图,即可求出修正后的待测天线两个极化分量方向图结果。
4工程应用算法简化
通常采用开口波导作为平面近场探头方向图,其方向图波束非常宽,在正前方正负10度范围内,角度变化0.1度,幅度最大变化0.0071dB。对于关心角度较窄的天线,在探头倾斜很小的情况下,导致的幅度影响可以忽略不计。
式8就可化简为
对式9-1与式9-3左侧幅度、相位分别作平均运算,等号左边为同一个物理量所以平均后不变。实际中D0与D180为同样的极化的近场不同网格原点的数据,根据乃奎斯特采样定理,其经过傅里叶变换后幅度理论上应该完全一样,只是由于网格原点不同导致输出相位不同。整理得到式10。
同理对式9-2与式9-4进行同样的操作,得到式11
其中,am_h、ph_h表示修正后的水平极化的幅度和相位,am_v、ph_v表示修正后的垂直极化的幅度和相位;am_fft表示傅里叶变换后取幅度部分,ph_fft表示傅里叶变换后取相位部分;e为自然底数;为待测天线方向图;为探头0度位置方向图,为探头90度位置方向图。
式10与式11即为修正探头偏心距离后的结果,联立即可求得修正后的待测天线方向图。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知常识。
Claims (3)
1.一种平面近场探头安装精度引起的测试误差修正方法,其特征在于,包括步骤如下:
步骤一、按照平面近场测试方法将待测天线进行安装,场地校准;
步骤二、在探头初始位置,设定扫描范围、测试距离,进行第一个极化分量数据的采集,得到近场数据D0;探头旋转90度,进行第二极化分量数据的测量,得到近场数据D90;探头继续旋转90度,进行第三极化分量数据的测量,得到近场数据D180;探头继续旋转90度,进行第四极化分量数据的测量,得到近场数据D270;
步骤三、进行近远场变换与数据修正,得到最终修正后的待测天线方向图数据,计算公式如下:
其中,am_h=(am_fft(D0)+am_fft(D180))/2,
ph_h=(ph_fft(D0)+ph_fft(D180))/2-π/2,
am_v=(am_fft(D90)+am_fft(D270))/2,
ph_v=(ph_fft(D90)+ph_fft(D270))/2-π/2,
am_h表示修正后的水平极化的幅度,ph_h表示修正后的水平极化的相位,am_v表示修正后的垂直极化的幅度,ph_v表示修正后的垂直极化的相位;am_fft表示傅里叶变换后取幅度部分,ph_fft表示傅里叶变换后取相位部分;e为自然底数;为待测天线方向图,为探头0度位置方向图,为探头90度位置方向图,i为复数的虚部单位量,i2=-1。
2.根据权利要求1所述的一种平面近场探头安装精度引起的测试误差修正方法,其特征在于:所述步骤二中扫描范围采用截断电平小于-35dB。
3.根据权利要求1或2所述的一种平面近场探头安装精度引起的测试误差修正方法,其特征在于:所述步骤二中测试距离为距离测试天线口面3~10个波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610710111.0A CN106338655B (zh) | 2016-08-23 | 2016-08-23 | 一种平面近场探头安装精度引起的测试误差修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610710111.0A CN106338655B (zh) | 2016-08-23 | 2016-08-23 | 一种平面近场探头安装精度引起的测试误差修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106338655A CN106338655A (zh) | 2017-01-18 |
CN106338655B true CN106338655B (zh) | 2018-10-09 |
Family
ID=57824709
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610710111.0A Active CN106338655B (zh) | 2016-08-23 | 2016-08-23 | 一种平面近场探头安装精度引起的测试误差修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106338655B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107219410B (zh) * | 2017-06-21 | 2019-09-06 | 西安空间无线电技术研究所 | 一种基于探头扫频位移偏移量的平面近场测量修正方法 |
CN108008331B (zh) * | 2017-11-17 | 2020-05-22 | 中国信息通信研究院 | 一种针对多探头球面近场的探头一致性检测方法 |
CN108663577B (zh) * | 2018-03-29 | 2020-06-09 | 北京空间飞行器总体设计部 | 一种多探头球面近场交叉极化校正方法 |
CN110631523B (zh) * | 2019-10-16 | 2021-04-20 | 大连理工大学 | 一种精密二维转台轴间位置误差的测量装置及方法 |
CN112269070B (zh) * | 2020-08-31 | 2023-12-12 | 西安空间无线电技术研究所 | 一种天线测试双极化探头的通道平衡方法和系统 |
CN113640589B (zh) * | 2021-07-21 | 2023-08-25 | 上海机电工程研究所 | 基于辐射信号监测的偏心测量补偿系统、方法及介质 |
CN116559745A (zh) * | 2023-04-12 | 2023-08-08 | 成都飞机工业(集团)有限责任公司 | 一种平面近场法天线测量中的扫描探头修正方法 |
CN117233683B (zh) * | 2023-11-14 | 2024-01-26 | 陕西拾贝通讯技术有限公司 | 基于距离点云数据对天线平面近场采集数据的校准方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323685A (zh) * | 2013-07-16 | 2013-09-25 | 南京顺仕祥电子有限公司 | 天线平面近场测试扫描仪 |
CN103411528A (zh) * | 2013-08-26 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | 利用圆极化天线轴比方向图计算电场探头旋转偏移的方法 |
CN103424066A (zh) * | 2013-08-26 | 2013-12-04 | 中国科学院空间科学与应用研究中心 | 利用圆极化天线远场相位差计算电场探头旋转偏移的方法 |
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
CN105353229A (zh) * | 2015-10-20 | 2016-02-24 | 上海无线电设备研究所 | 一种基于一维旋转的相控阵幅相误差近场校准方法 |
CN105445566A (zh) * | 2015-11-13 | 2016-03-30 | 西北工业大学 | 球面多探头天线测试数据处理方法 |
JP6028964B2 (ja) * | 2012-06-19 | 2016-11-24 | 国立研究開発法人情報通信研究機構 | ニアフィールド測定装置 |
-
2016
- 2016-08-23 CN CN201610710111.0A patent/CN106338655B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6028964B2 (ja) * | 2012-06-19 | 2016-11-24 | 国立研究開発法人情報通信研究機構 | ニアフィールド測定装置 |
CN103323685A (zh) * | 2013-07-16 | 2013-09-25 | 南京顺仕祥电子有限公司 | 天线平面近场测试扫描仪 |
CN103411528A (zh) * | 2013-08-26 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | 利用圆极化天线轴比方向图计算电场探头旋转偏移的方法 |
CN103424066A (zh) * | 2013-08-26 | 2013-12-04 | 中国科学院空间科学与应用研究中心 | 利用圆极化天线远场相位差计算电场探头旋转偏移的方法 |
CN104391183A (zh) * | 2014-10-11 | 2015-03-04 | 中国船舶重工集团公司第七二四研究所 | 近场测量天线远区场特性快速计算方法 |
CN105353229A (zh) * | 2015-10-20 | 2016-02-24 | 上海无线电设备研究所 | 一种基于一维旋转的相控阵幅相误差近场校准方法 |
CN105445566A (zh) * | 2015-11-13 | 2016-03-30 | 西北工业大学 | 球面多探头天线测试数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106338655A (zh) | 2017-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106338655B (zh) | 一种平面近场探头安装精度引起的测试误差修正方法 | |
CN105404238B (zh) | 一种在机激光测量中测头位姿的线性化标定方法 | |
CN110618408B (zh) | 一种精密测距系统天线相位中心的系统标定方法 | |
CN101551240B (zh) | 基于激光跟踪技术的大型齿轮测量方法 | |
CN108037374B (zh) | 一种阵列天线近场标定方法 | |
CN106052556B (zh) | 一种三坐标测量机空域坐标修正方法 | |
CN110926364B (zh) | 基于线结构光的叶片检测方法 | |
CN108168465A (zh) | 一种共光路激光外差干涉法滚转角高精度测量装置及方法 | |
CN110926365B (zh) | 一种基于线结构光检测物标定方法 | |
CN112526486B (zh) | 一种基于轴系误差模型的三维激光雷达空间坐标校准方法 | |
CN109884659B (zh) | 基于激光追踪仪多站位测量系统的大型精密转台标定方法 | |
CN109884658B (zh) | 基于激光追踪仪多站位测量系统的激光追踪仪站位方法 | |
CN111044077B (zh) | 星敏感器测量坐标系与星敏感器立方镜坐标系间的标定方法 | |
CN104535974A (zh) | 一种飞机雷达系统校靶装置及其使用方法 | |
CN114994420A (zh) | 相控阵天线阵元通道相位偏转特性测试方法 | |
CN109141868B (zh) | 精密轴系误差运动的测量装置及测量方法 | |
CN110640546B (zh) | 用于大型齿轮在机旁置测量的被测齿轮回转轴线测定方法 | |
CN105937885B (zh) | 一种自由曲面子孔径拼接干涉检测中被测面定位匹配方法 | |
CN116256705A (zh) | 一种机相扫高精度跟踪雷达轴系匹配校准方法 | |
CN116359615A (zh) | 一种平面近场测量中的准直装置及方法 | |
CN107544554B (zh) | 一种复合天线电轴重合简易调整方法 | |
CN110631523B (zh) | 一种精密二维转台轴间位置误差的测量装置及方法 | |
CN102445162B (zh) | 一种旋转激光束呈锥面扫描的测量装置的参数标定方法 | |
CN115014257B (zh) | 一种二维线激光传感器标定方法 | |
CN110879382A (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 |