CN106209269A - 一种射频仿真***中球面复合阵列近场效应的校准方法 - Google Patents
一种射频仿真***中球面复合阵列近场效应的校准方法 Download PDFInfo
- Publication number
- CN106209269A CN106209269A CN201610601087.7A CN201610601087A CN106209269A CN 106209269 A CN106209269 A CN 106209269A CN 201610601087 A CN201610601087 A CN 201610601087A CN 106209269 A CN106209269 A CN 106209269A
- Authority
- CN
- China
- Prior art keywords
- tlv triple
- antenna
- composite array
- field effect
- center
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/10—Monitoring; Testing of transmitters
- H04B17/11—Monitoring; Testing of transmitters for calibration
- H04B17/12—Monitoring; Testing of transmitters for calibration of transmit antennas, e.g. of the amplitude or phase
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种射频仿真***中球面复合阵列近场效应的校准方法,属于仪器设备校准技术领域。球面复合阵列是由微波与毫米波天线阵列组成,从而产生微波与毫米波频段的电磁波信号的装置。校准方法的步骤:联合全波算法和高频算法精确计算复合阵列中给定子阵(三元组)天线在周围同频、异频天线单元存在的环境中发射时,接收天线位置处(通常位于发射三元组天线的辐射近场)的电场分布;基于此电场相位分布,应用相位梯度法求出三元组辐射场的等效相位中心,并计算与基于重心公式所得等效相位中心估计值的偏差;根据偏差值校准三元组各单元输入功率比例。本发明可制成校准表,便于在工程中对复合阵列近场效应进行快速、实时校准。
Description
技术领域
本发明属于仪器设备校准技术领域,涉及一种射频仿真***中球面复合阵列近场效应的校准方法。
背景技术
射频半实物仿真***在导弹制导***性能评估、***参数调试与验证、飞行试验结果预测与故障复现等方面有着不可替代的作用,不但能够有效的缩短导弹研制周期,还可大大节约科研经费。半实物仿真能够在雷达***仿真中,解决数学仿真中某些部分不易建模和物理仿真中费用较高,结构难以进行修改的问题。半实物射频仿真***一般包括以下六个部分:电磁波暗室、射频目标仿真器、计算机及其接口、目标和干扰环境模型的数据库及相应的软件、监控操作台及显示设备、校准***。射频仿真中,电磁波信号由暗室中的球面天线阵列中的一个相应子阵(由三个辐射天线组成,以下称三元组)辐射,通过馈电装置给予目标位置三元组一定比例的功率,射频目标仿真器便能产生目标信号,提供试验所需要的外界电磁环境。由于电磁波暗室和仿真测试***造价昂贵,通常半实物仿真***并不能满足测试所需的“远场条件”,待测***一般是工作于三元组的“辐射近场”区,为保证测试精度,必须对仿真***的近场效应进行有效的校正;随着应用发展需要,试验要求射频半实物仿真***能够提供更加复杂的电磁环境,在同一个发射天线阵列上实现多个频段的目标信号的球面复合阵列也已提出,单阵列近场效应的校正方法不能满足球面复合阵列近场效应的校正精度。
发明内容
本发明的目的是提供一种射频仿真***中球面复合阵列近场效应的校准方法。
本发明要解决的技术问题是:
在被测单元处实现较高测角精度,校准了因为球面复合阵列近场效应所带来的误差,通过误差校准表,实现半实物仿真***中各发射“三元组”的快速、准确控制。
本发明所采用的技术方案是:
一种射频仿真***中球面复合阵列近场效应的校准方法,该方法包括以下具体步骤:
步骤1:以发射天线球面阵列球心即被测目标的接收天线位置为坐标原点,根据所需电磁波信号的方位角与俯仰角,用重心公式计算阵列中对应三元组(球面复合阵列上任意相邻的三个同频天线组成的子阵)各单元的输入功率初值,以及在被测目标处来波等效相位中心的估计值;
步骤2:使用重心公式所得三元组各单元输入功率的初值,用全波方法和一致性几何绕射理论(UTD)相结合的混合算法,计算工作三元组及近旁同频和异频单元存在时的复杂电磁环境下接收天线附近的三元组辐射近场的幅度和相位分布;
步骤3:根据计算所得的三元组辐射场的相位分布,
用相位梯度法计算三元组电磁波信号的等效相位中心的准确值与估计值的方位偏差lx和俯仰偏差ly,lx与ly分别除以阵面半径便得到方位方向和俯仰方向的测角偏差,即所谓“测角误差”;
步骤4:如测角误差小于所给定的误差范围,则以重心公式所得的功率作为三元组各单元实际输入功率;否则按修正公式对三元组各单元的输入功率值进行修正,得到更精确的输入功率值,重复上述步骤1至步骤3直到获得满足给定误差范围的测角误差;其中,修正公式如下所示:
三元组中三个天线分别命名为Mic1、Mic2以及Mic3,其相应的输入幅度分别为M1、M2以及M3,以Mic1的输入幅度为基准,即M1=1,L为三元组中任意两天线连线长度,M3和M2为指定等效相位中心位置处Q(τQ,υQ)并考虑了由于近场效应引起的方位方向和俯仰方向的偏差(lx,ly)的修正后的输入幅度;根据修正后的输入幅度能够计算三元组各单元相应的修正后输入功率。
与现有技术相比,本发明的有益效果是:克服了射频半实物仿真***中球面复合阵列近场效应所带来的误差。
附图说明
图1球面复合阵列局部示意图;
图2本发明流程图;
图3收发天线(以接收天线中点为坐标原点)关系图;
图4收发天线(以三元组等效相位中心点为坐标原点)关系图;
图5 Q点位置修正计算原理图;
图6球面复合阵列中微波三元组近场效应修正前后比较图;
图7球面复合阵列中毫米波三元组近场效应修正前后比较图。
具体实施方式
本发明所述的球面复合阵列由微波阵列与毫米波阵列组成,如图1所示,显示了球面复合阵列局部;图中,●为微波天线,■为毫米波天线,a为微波三元组,b为毫米波三元组。
参阅图2,本发明的流程如下:
(1)以接收天线中心为坐标原点O建立笛卡尔坐标系,天线Mic1、Mic2、Mic3组成三元组并且分布在半径为R的阵面上,接收天线轴线指向Mic1,如图3所示,描述了接收天线与发射天线之间的关系,图中,◇为接收天线,●为微波天线,■为毫米波天线,微波天线与毫米波天线均表示发射天线。选择三元组所在面上任意一点Q(如图5所示)作为等效相位中心点,即目标信号的发射位置。根据以下重心公式即可获得归一化的三元组各单元的激励幅度比例,根据归一化的幅度比例计算相应的功率,作为三元组各单元输入功率初值。
式中的φ与θ分别表示该坐标系下等效相位中心的方位角和俯仰角,φi与θi分别表示三元组各单元的方位角与俯仰角,等效相位中心方位角和俯仰角的变化范围为固定值分别为Δφ和Δθ,Ei表示三元组各单元的幅度。
(2)以上一步等效相位中心Q点为坐标原点O建立笛卡尔坐标系,建立三元组仿真模型,接收天线位于辐射近场区,如图4所示,其中O’为仿真计算得到的等效相位中心的位置。给予三元组上一步计算的输入功率初值,接着仿真得到接收天线处的幅度和相位分布,根据相位梯度法计算三元组仿真得到的新的等效相位中心点。
上式的含义是仿真计算得到的新的等效相位中心点的方位偏差和俯仰偏差。该坐标系中的ψ和η分别表示接收天线中心处的方位角和俯仰角,α和β分别表示接收天线位置处相位沿着η和ψ方向的梯度除以传播常数k。
上式中dx和dy分别表示方位方向和俯仰方向的测角误差,其中R表示球面复合阵列的阵面半径。
(3)判断测角误差是否超过给定的误差范围,如果超过了给定的误差范围则需要求取该等效相位中心点的修正的幅度比例。当等效相位中心点位于点Q(τQ,υQ)时,方位偏差和俯仰偏差分别为lx和ly,σ值为Q点位置的角度,L为三元组天线中任意两天线之间的间距。具体如图5所示,横坐标τ表示天线1与天线3连线上不同位置距离天线1的长度所占总长度的比例值,纵坐标υ表示天线1天线3连线中点与天线2连线上不同位置与连线中点的长度所占总长度的比例值。
以Mic1的输入幅度作为基准,Mic1幅度M1=1;先修正等效相位中心点在方位方向的偏差,再修正因为Mic2的引入而存在的沿着过Q点的斜线方向的偏差。方位方向的偏差值为lx与ly的方位方向分量的和为lx+ly.tanσ,Mic2与Q连线方向存在的偏差值为俯仰方向偏差值沿着该连线方向的一个分量为ly/cosσ。经过化简整理可得修正后的天线2,天线3的辅助修正公式:
求修正后归一化幅度比例对应的功率值并且作为三元组实际的输入功率,返回到(2)计算测角偏差,循环直到获得满足给定测角偏差范围的三元组输入功率。
(4)由以上三个步骤可以求取三个天线连线区域内任意位置作为等效相位中心点的修正后的三个天线幅度比例,接着归一化幅度比例,给出相对应的功率比例并且制成修正表以供试用。
根据上述球面复合阵列近场效应校准的流程分析的实例见实施例1。
本发明结合以下实施例1作具体说明,但不用来限制本发明的保护范围。
实施例1
球面复合阵列中微波/毫米波三元组近场效应校准,包括以下内容:
首先建立仿真模型,选择球面复合阵列中任意一组相邻的微波/毫米波天线组成三元组并且给予其输入功率初值,其它位置的微波天线和毫米波天线不给予输入功率。接收天线轴线指向三元组天线Mic1。
将等效相位中心设置为三元组所在平面两天线连线中垂线上,中垂线分为十等份,九个间隔点便是等效相位中心的位置。根据重心公式获得目标信号来源于中垂线不同位置时三元组各单元应给予的输入功率初值,仿真、计算获得目标信号来源于中垂线上不同位置时的测角误差值,接着根据球面复合阵列近场效应修正原理得到三元组输入功率的修正值,仿真、计算获得修正后的目标信号来源于中垂线上不同位置时的测角误差值。因为目标信号来源于中垂线上不同位置时天线1和天线2的对称性,所以测角误差主要发生在俯仰方向,方位方向测角误差可以忽略不计。如图6所示,显示了微波三元组修正前后中垂线上俯仰测角误差的变化情况,其中▲线条表示修正前俯仰方向测角误差,◆线条表示修正后俯仰方向测角误差,通过近场效应修正使得最大等效相位中心俯仰测角误差降低了97.50%;如图7所示,显示了毫米波三元组修正前后中垂线上俯仰测角误差的变化情况,其中▲线条表示修正前俯仰方向测角误差,◆线条表示修正后俯仰方向测角误差,通过近场效应修正使得最大等效相位中心俯仰测角误差降低了99.76%。
Claims (1)
1.一种射频仿真***中球面复合阵列近场效应的校准方法,其特征在于,该方法包括以下步骤:
步骤1:以发射天线球面阵列球心即被测目标的接收天线位置为坐标原点,根据所需电磁波信号的方位角与俯仰角,用重心公式计算阵列中对应三元组各单元的输入功率初值,以及在被测目标处来波等效相位中心的估计值;其中,所述三元组为球面复合阵列上任意相邻的三个同频天线组成的子阵;
步骤2:使用重心公式所得三元组各单元输入功率的初值,用全波方法和一致性几何绕射理论(UTD)相结合的混合算法,计算工作三元组及近旁同频和异频单元存在时的复杂电磁环境下接收天线附近的三元组辐射近场的幅度和相位分布;
步骤3:根据计算所得的三元组辐射场的相位分布,用相位梯度法计算三元组电磁波信号的等效相位中心的准确值与估计值的方位偏差lx和俯仰偏差ly;lx与ly分别除以阵面半径便得到方位方向和俯仰方向的测角偏差,即所谓“测角误差”;
步骤4:若测角误差小于所给定的误差范围,则以重心公式所得的功率作为三元组各单元实际输入功率;否则按修正公式对三元组各单元的输入功率值进行修正,得到更精确的输入功率值,重复上述步骤1至步骤3直至获得满足给定误差范围的测角误差;其中,修正公式如下所示:
三元组中三个天线分别命名为Mic1、Mic2以及Mic3,其相应的输入幅度分别为M1、M2以及M3,以Mic1的输入幅度为基准,即M1=1,L为三元组中任意两天线连线长度,M3和M2为指定等效相位中心位置处Q(τQ,υQ)并考虑了由于近场效应引起的方位方向和俯仰方向的偏差(lx,ly)的修正后的输入幅度;根据修正后的输入幅度能够计算三元组各单元相应的修正后输入功率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610601087.7A CN106209269B (zh) | 2016-07-27 | 2016-07-27 | 一种射频仿真***中球面复合阵列近场效应的校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610601087.7A CN106209269B (zh) | 2016-07-27 | 2016-07-27 | 一种射频仿真***中球面复合阵列近场效应的校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106209269A true CN106209269A (zh) | 2016-12-07 |
CN106209269B CN106209269B (zh) | 2018-09-28 |
Family
ID=57495482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610601087.7A Active CN106209269B (zh) | 2016-07-27 | 2016-07-27 | 一种射频仿真***中球面复合阵列近场效应的校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106209269B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107239044A (zh) * | 2017-06-05 | 2017-10-10 | 北京仿真中心 | 一种限定性求解矢量控制近场效应修正初值的方法 |
CN107942145A (zh) * | 2017-11-22 | 2018-04-20 | 北京仿真中心 | 一种角变化球面三元组天线辐射电磁场测量***和方法 |
CN108021037A (zh) * | 2017-11-24 | 2018-05-11 | 上海机电工程研究所 | 三元组目标阵列近场效应修正方法 |
CN108051801A (zh) * | 2017-12-05 | 2018-05-18 | 南京长峰航天电子科技有限公司 | 一种基于微波与毫米波复合的高精度射频仿真*** |
CN108062036A (zh) * | 2017-12-21 | 2018-05-22 | 南京长峰航天电子科技有限公司 | 一种用于Ku/Ka/W波段仿真试验的复合阵列***及实现方法 |
CN110632399A (zh) * | 2019-09-19 | 2019-12-31 | 电子科技大学 | 用于球面近场测量数据的修正方法及天线方向图测量方法 |
CN110717234A (zh) * | 2019-10-17 | 2020-01-21 | 上海机电工程研究所 | 非规则布局三元组角位置模拟方法、***及介质 |
CN112710903A (zh) * | 2019-10-25 | 2021-04-27 | 上海益麦电磁技术有限公司 | 一种基于多探头球面场天线测试***的相位中心测试方法 |
CN112947119A (zh) * | 2021-03-08 | 2021-06-11 | 中国人民解放军63892部队 | 一种射频半实物仿真数字化阵列实施***及方法 |
CN114280957A (zh) * | 2021-12-22 | 2022-04-05 | 南京长峰航天电子科技有限公司 | 一种射频阵列的数字化建模方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682165A (en) * | 1996-05-02 | 1997-10-28 | Hughes Electronics | Active array self calibration |
CN1740812A (zh) * | 2005-06-17 | 2006-03-01 | 武汉大学 | 高频表面波雷达均匀直线阵接收通道的近场校正方法 |
US20110133982A1 (en) * | 2007-12-31 | 2011-06-09 | Elta Systems Ltd. | System and Method for Calibration of Phased Array Antenna Having Integral Calibration Network in Presence of an Interfering Body |
CN103926474A (zh) * | 2014-03-18 | 2014-07-16 | 中国电子科技集团公司第十研究所 | 相控阵天线单元特性近场测量方法 |
CN104063544A (zh) * | 2014-06-25 | 2014-09-24 | 电子科技大学 | 近场效应误差分析方法 |
CN105162536A (zh) * | 2015-08-21 | 2015-12-16 | 西安空间无线电技术研究所 | 一种相控阵天线在轨幅相校正***及方法 |
CN105353229A (zh) * | 2015-10-20 | 2016-02-24 | 上海无线电设备研究所 | 一种基于一维旋转的相控阵幅相误差近场校准方法 |
-
2016
- 2016-07-27 CN CN201610601087.7A patent/CN106209269B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5682165A (en) * | 1996-05-02 | 1997-10-28 | Hughes Electronics | Active array self calibration |
CN1740812A (zh) * | 2005-06-17 | 2006-03-01 | 武汉大学 | 高频表面波雷达均匀直线阵接收通道的近场校正方法 |
US20110133982A1 (en) * | 2007-12-31 | 2011-06-09 | Elta Systems Ltd. | System and Method for Calibration of Phased Array Antenna Having Integral Calibration Network in Presence of an Interfering Body |
CN103926474A (zh) * | 2014-03-18 | 2014-07-16 | 中国电子科技集团公司第十研究所 | 相控阵天线单元特性近场测量方法 |
CN104063544A (zh) * | 2014-06-25 | 2014-09-24 | 电子科技大学 | 近场效应误差分析方法 |
CN105162536A (zh) * | 2015-08-21 | 2015-12-16 | 西安空间无线电技术研究所 | 一种相控阵天线在轨幅相校正***及方法 |
CN105353229A (zh) * | 2015-10-20 | 2016-02-24 | 上海无线电设备研究所 | 一种基于一维旋转的相控阵幅相误差近场校准方法 |
Non-Patent Citations (2)
Title |
---|
S.XIAO,ETC: ""Focusing pulsed electromagnetic radiation in the near field"", 《PULSED POWER CONFERENCE,2009.PPC 09.IEEE》 * |
马静等: ""面向引用的射频仿真***近场效应误差修正"", 《***工程与电子技术》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107239044A (zh) * | 2017-06-05 | 2017-10-10 | 北京仿真中心 | 一种限定性求解矢量控制近场效应修正初值的方法 |
CN107239044B (zh) * | 2017-06-05 | 2020-06-02 | 北京仿真中心 | 一种限定性求解矢量控制近场效应修正初值的方法 |
CN107942145A (zh) * | 2017-11-22 | 2018-04-20 | 北京仿真中心 | 一种角变化球面三元组天线辐射电磁场测量***和方法 |
CN108021037B (zh) * | 2017-11-24 | 2021-01-22 | 上海机电工程研究所 | 三元组目标阵列近场效应修正方法 |
CN108021037A (zh) * | 2017-11-24 | 2018-05-11 | 上海机电工程研究所 | 三元组目标阵列近场效应修正方法 |
CN108051801A (zh) * | 2017-12-05 | 2018-05-18 | 南京长峰航天电子科技有限公司 | 一种基于微波与毫米波复合的高精度射频仿真*** |
CN108062036A (zh) * | 2017-12-21 | 2018-05-22 | 南京长峰航天电子科技有限公司 | 一种用于Ku/Ka/W波段仿真试验的复合阵列***及实现方法 |
CN110632399A (zh) * | 2019-09-19 | 2019-12-31 | 电子科技大学 | 用于球面近场测量数据的修正方法及天线方向图测量方法 |
CN110717234A (zh) * | 2019-10-17 | 2020-01-21 | 上海机电工程研究所 | 非规则布局三元组角位置模拟方法、***及介质 |
CN112710903A (zh) * | 2019-10-25 | 2021-04-27 | 上海益麦电磁技术有限公司 | 一种基于多探头球面场天线测试***的相位中心测试方法 |
CN112947119A (zh) * | 2021-03-08 | 2021-06-11 | 中国人民解放军63892部队 | 一种射频半实物仿真数字化阵列实施***及方法 |
CN112947119B (zh) * | 2021-03-08 | 2022-10-25 | 中国人民解放军63892部队 | 一种射频半实物仿真数字化阵列实施***及方法 |
CN114280957A (zh) * | 2021-12-22 | 2022-04-05 | 南京长峰航天电子科技有限公司 | 一种射频阵列的数字化建模方法及*** |
CN114280957B (zh) * | 2021-12-22 | 2024-02-13 | 南京长峰航天电子科技有限公司 | 一种射频阵列的数字化建模方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN106209269B (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106209269A (zh) | 一种射频仿真***中球面复合阵列近场效应的校准方法 | |
US11309977B2 (en) | Calibration method and communications device | |
CN107329003B (zh) | 优化sar天线方向图测试的方法 | |
CN109495189A (zh) | 一种阵列天线校准方法及装置 | |
US10302735B2 (en) | Wireless terminal antenna directional characteristic measurement system and wireless terminal antenna directional characteristic measurement method | |
CN108037374B (zh) | 一种阵列天线近场标定方法 | |
CN103217589B (zh) | 一种相控阵列天线等效隔离度测试方法 | |
WO2020001661A1 (zh) | 一种测量方法及设备 | |
US11519949B2 (en) | Measurement system and method of performing an over-the-air test | |
CN107339991A (zh) | 一种飞行器航向角的检测方法及装置 | |
CN106526562A (zh) | 一种基于天线方向图修正计算目标rcs的方法 | |
CN104063544A (zh) | 近场效应误差分析方法 | |
CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
CN110045340A (zh) | 平面相控阵天线中信标天线的坐标测算方法 | |
CN106680789B (zh) | 一种单脉冲测角设备的测角功能验证方法和验证*** | |
CN104375129A (zh) | 一种分布式阵列相参合成雷达发射相参参数校准方法 | |
CN102608591A (zh) | 一种基于相位修正的rcs外推方法 | |
CN102621560B (zh) | 一种基于微波暗室的大场景卫星星座模拟方法 | |
Kim et al. | Phase effects analysis of patch antenna CRPAs for JPALS | |
CN103544376A (zh) | 短波固定监测站测向数据修正方法 | |
CN111965602B (zh) | 一种相控阵雷达幅相一致性检测方法和*** | |
CN105676226A (zh) | 一种射频阵列天线校准装置 | |
CN112711040B (zh) | 卫星导航天线性能评估***及方法 | |
CN110717234A (zh) | 非规则布局三元组角位置模拟方法、***及介质 | |
Gu et al. | Analysis of amplitude-phase error of phased array calibration in mid-field |
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 |