CN104063544B - 近场效应误差分析方法 - Google Patents
近场效应误差分析方法 Download PDFInfo
- Publication number
- CN104063544B CN104063544B CN201410293567.2A CN201410293567A CN104063544B CN 104063544 B CN104063544 B CN 104063544B CN 201410293567 A CN201410293567 A CN 201410293567A CN 104063544 B CN104063544 B CN 104063544B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- array elements
- radiation
- 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
- 230000000694 effects Effects 0.000 title claims abstract description 11
- 238000004458 analytical method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims abstract description 21
- 238000012937 correction Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000005855 radiation Effects 0.000 claims description 37
- 230000005684 electric field Effects 0.000 claims description 18
- 230000005669 field effect Effects 0.000 claims description 17
- 238000009826 distribution Methods 0.000 claims description 16
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 16
- 230000001808 coupling effect Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000005672 electromagnetic field Effects 0.000 description 2
- 238000005290 field theory Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明涉及近场效应误差分析方法,属于射频仿真领域。国内现有的三元阵天线近场效应误差分析等工作均是基于解析方法,尽管简单,但计算结果不精确、所适用的频段范围有限等,因此为了满足射频仿真中的高频率、高精度要求。本发明从三阵元的三个阵元出发,针对实际应用,发展计算电磁学数值方法同时结合电磁仿真软件模拟三元阵天线的近场效应,生成方位角和俯仰角的近场效应误差修正表格,通过近场效应误差修正得到正确的目标复现位置。本发明方法简单易行,满足了射频仿真的高频率、高精度要求。
Description
技术领域
本发明涉及射频仿真系统,主要涉及三阵元近场效应误差导致的模拟目标位置误差分析方法。
背景技术
阵列式射频目标仿真具有便于模拟复杂目标和环境的优点,是目前国内外都相继采用的仿真方法。国外对此技术的应用已经比较成熟。美国射频仿真系统先后建成了微波波段和两个毫米波波段的半实物仿真系统,我国也建设了各射频仿真中心,但是总体上和美、俄等国还存在差距。国内通过多年的研究对这类系统的工作原理、基本设计要求、方法等都已掌握。目前研究方向集中在解决和完善系工程设计所面临的关键技术,如近场效应误差及修正技术。
发明内容
本发明的目的在于克服现有技术中所存在的上述不足,为了满足射频仿真中的高频率、高精度要求,还需要根据严格的电磁场理论进行精确计算等。本项目针对实际应用,拟发展计算电磁学数值方法同时结合电磁仿真软件模拟三元阵天线的近场效应,实时直观地研究电磁波传播的演进过程,通过近场效应误差修正得到正确的目标复现位置,深入研究其近场效应机理,提供一种近场效应误差分析方法。
为了实现上述发明目的,本发明提供了以下技术方案:
本发明近场效应误差分析方法,分步骤地逐步推进工作,首先研究基于三阵元幅度控制方法近场效应误差修正理论,继而基于计算电磁学方法结合误差修正理论研究射频仿真系统中的误差修正方案,探索复杂目标条件下的近场效应误差修正理论及修正方法,主要包括以下步骤:
步骤一:
(1)、建立坐标系:三个辐射阵元呈等边三角形分布在二维直角坐标系x0y平面内,上述三个辐射阵元的中心位于坐标原点O;三个辐射阵元在直角坐标系x0y中的坐标为(xi,yi),i=1、2、3;固定观测点P在以O为坐标原点的三维坐标系中的位置为(0,0,R);任意观测点P′在以O点为球心的球坐标系2中的位置为(α,β);任意观测点P′在以O为坐标原点的三维坐标系中的位置为:
(2)、将观测点从固定观测点P运动到任意观测点P′,得到观测点到三个辐射阵元的距离变化量ri为:
(3)、定义三个阵元辐射的电磁波角频率为ω,初始相位为φi,在任意观测点P′处接收到三个辐射阵元的辐射信号为:
Ei=Aiexp[j(ωt-kR-kri+φi)],
其中,Ai(i=1,2,3)为三个辐射阵元的馈电幅度;j为虚数单位,j的平方为-1;t是时间变量;k为相位常数,且k=2π/λ,λ表示三阵元辐射电磁波的波长;φi(i=1,2,3)为三个辐射阵元馈电的初始相位;
(4)、将三个阵元等效为阵面上的一个辐射源,任意观测点P′处接收的辐射信号E为:
E=Aexp[j(ωt-kR+δ)],
其中,A为三阵元等效辐射源的馈电幅度;k为相位常数,且k=2π/λ;R为三阵元中心O点到固定观测点P的距离;δ为在任意观测点P′处三阵元辐射信号合成后的等相位面相对于在三阵元中心O处的参考源辐射的等相位面的倾斜;
(5)、改变三阵元的馈电幅度和相位中的一种或者两种,从而控制在任意观察点P′处测得的等相位面的位置,所述的等相位面的位置与馈电幅度和相位的关系如下式:
其中,n=3;
(6)、改变三阵元的馈电幅度和相位中的一种或者两种,从而控制目标的位置,所述的目标的位置与馈电幅度和相位的关系如下式:
其中,(i=1,2,3)是三个辐射阵元在以P点为球心的球坐标系1中的坐标,是三个辐射阵元的等效辐射源在以P点为球心的球坐标系1中的坐标,n=3;
当φi=0(i=1,2,3)时,上式可简化为:
步骤二:
三阵元辐射阵元i的馈电幅度Ai(i=1,2,3)和相位φi(i=1,2,3),得到辐射阵元i在接收天线上接收到的平均电场m表示近场效应修正算法流程中的迭代序列,其中是第m次迭代中天线阵元i在接收天线口径面上产生的场强,Es(ρ,ζ)是接收天线的口径面场分布,S口径表示天线口径面的面积,∫∫口径面表示对口径面上的积分计算;
步骤三:
令m=m+1,导引头的口径面场分布为:
利用电磁仿真可以得到导引头口径面上若干个离散点处的电场,将计算平均电场分布的积分离散化,避免推导繁复的解析公式,N值越大,计算结果越精确。通过分析阵元间残留互耦,具有一定口径的接收天线,辐射单元的间距、工作波长、接收天线口径的尺寸以及天线的形式、方向图因素,得到误误差产生的原因和怎样消去或者减小误差。
三阵元近场条件下,三阵元天线辐射信号由于光程差会在接收天线口面处形成相位差,等效辐射相位波前不是理想球面,而是近似球面。要通过三阵元幅度控制获得最好的相位波前球面近似,就是要寻找一个最佳幅度控制样本,使得三阵元合成信号在接收天线口面处形成的最大相位差最小。为了提高精度,需要根据严格的电磁场理论进行精确计算,然后制成表格对数据进行实时修正。从三阵元的三个阵元出发,求出接收天线阵面上四个象限接收到的合成场强,利用接收天线的和信号方向图、方位差信号和俯仰差信号方向图,得到接收天线的方位和俯仰误差角,通过调整三元天线组振幅幅值,使到达角误差趋于零等生成近场效应误差修正表格。同时分析阵元间残留互耦的影响,综合考虑具有一定口径的接收天线,多个辐射单元的合成场在接收天线口径面上各个点的相位波前畸变不同且不对称,产生近场效应误差的辐射单元的间距、工作波长、接收天线口径的尺寸以及天线的形式、方向图等因素
与现有技术相比,本发明的有益效果:满足了射频仿真的高频率、高精度要求。本项目针对实际应用,拟发展计算电磁学数值方法同时结合电磁仿真软件模拟三元阵天线的近场效应,实时直观地研究电磁波传播的演进过程,通过近场效应误差修正得到正确的目标复现位置,深入研究其近场效应机理。
附图说明:
图1三阵元模拟目标的方位角和俯仰角;
图2初步算法流程图;
图3阵元在接收天线口径面上产生的场;
图4天线阵元i在导引头口径平面的场分布求解示意图;
图5有源阵元方向图原理示意图;
图6子阵列示意图;
图7改进后的算法流程图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现的技术均属于本发明的范围。
如图1,左边为接收天线的口径平面,O′点为转台中心。右边为目标阵列中的一个三阵元,其中心为O。所模拟的目标为D,用O′D和O′O的夹角来表示目标的位置。本算法中用到两个坐标系,两个坐标系之间可以通过一定的对应关系相互转换,下面介绍两个坐标系的定义。
球坐标系,以转台中心O′为原点,球面阵所模拟的目标位置D由表示,三阵元重心O的对应的坐标为θ=0。
极坐标系,以导引头口径面中心P为原点,口径平面上任一点g的坐标为(ρ,ζ),ρ和ζ分别表示极坐标系下的极径和极角。
假设三阵元的三个天线阵元编号为i=1,2,3。它们的馈电信号幅度和相位分别为Ai、φi,三个天线阵元在球坐标系中的位置表示为已知三阵元的馈电相位与幅度时,可以利用如下公式(1)得出所模拟的目标在球坐标系中的坐标。
已知三阵元辐射阵元i的馈电幅度Ai和相位φi时,由以下公式(2)可以得到第m+1次迭代辐射阵元i在接收天线上接收到的平均电场其中是第m次迭代天线阵元i在接收天线口径面上(ρ,ζ)点产生的场强,Es(ρ,ζ)是接收天线的口径场分布。
本方案所采用的算法流程图如图2所示,下面将详细介绍算法的各个步骤所需量的定义及计算方法。
从流程图中可知,在第m+1次迭代中,输入三阵元的馈电方案可以由第m次迭代中求出的阵元i在接收天线上口径面上接收到的平均电场得到,又可以通过天线i在接收天线上的电场分布得到。
算法迭代过程中,每一步都需更新的参数由导引头接收到来自天线阵元i的场强Ei,m和接收天线的口径面场分布Es(ρ,ζ)。在仿真频率较低的时候,以上两个量都是通过解析公式来近似得到的。然而在毫米波频段,如果再采用传统的解析公式方法,会造成非常大的误差。而且不同的导引头上所采用的接收天线不相同,也就会有不同的口径面场分布,对于结构复杂的天线,推导Es(ρ,ζ)的解析表达式非常困难,推导过程中往往会采取很多近似手段,从而造成解析公式的不准确性。而且,根据解析公式推导出来的Es(ρ,ζ)不具有普遍适用性,从而导致开发的近场效应误差修正软件不具有普遍适用性。这样就大大增加了导引头校准的复杂度,增加了测试成本。因此考虑采用数值仿真方法来获取Ei,m、Es(ρ,ζ),从而提高算法精度。采用数值仿真的方法还有一个好处是,针对不同的接收天线类型,可以通用同一套程序来修正近场效应误差。下面将分别介绍Ei,m、Es(ρ,ζ)的获取方法。
阵元在接收天线上的场分布:
如上图3所示,三阵元天线阵元1、2、3分别在接收天线口径面上产生一个电场分布。由于辐射阵元不是各向同性辐射,它在导引头口径面上所产生的电场也并不均匀。
通常在射频仿真系统中,目标阵列和接收天线分别放在暗室的两端,而暗室的尺寸少则几十个波长,多则上百个波长。因此想直接用仿真软件或计算电磁学的方法直接计算阵元i在口径面上产生的电场分布是不大现实的,但可以间接的计算出口径面上产生的电场分布,如图4所示。
Es(ρ,ζ):根据天线的互易原理,Es(ρ,ζ)可以通过仿真接收天线的远场方向图得到。
Ei,m:以天线i在口径面上g点产生的电场为例,天线i到点g的距离为d。在距离天线i为R处的电场可由仿真软件或者计算电磁学的方法得到,设为Fi,那么阵元i在g点产生的电场为Fi(1/r)e-jβr。
实际应用中,辐射阵元i处于目标阵列中。相邻的辐射阵元对其具有耦合效应,此耦合不可忽视。因此需要修正上述所求得的Fi(1/r)e-jβr。
计入阵列耦合效应的阵元方向图:
在电磁仿真中,可以有两种方法可以计入阵列耦合效应对阵元方向图的影响:其一是采用有源阵元方向图,其二是利用周期性边界条件,模拟无限大阵列中间的一个阵元。本方案中采取有源阵元方向图的方法。下面介绍有源阵元方向图的基本原理。
图5中给出了天线阵有源阵元方向图,如何计入阵列其他阵元对天线1的耦合效应。以均匀直线阵为例,天线1馈电,而天线2,3,…,M接匹配负载时,研究阵列的方向图,所得到的方向图即为天线阵元1的有源阵元方向图。由于此时所测得的阵列辐射,已经包含了耦合电流IM1的作用,所以认为利用有源阵元方向图来计入阵列耦合的方法可行。
要进行三阵元近场误差效应修正,就需要求得目标阵列上每个阵元的有源阵元方向图。目标阵列是一个电大尺寸的阵列,对整个阵列进行一次仿真已经非常耗时,因此需要想法解决此问题。如图6,阵列中,阵元M与馈电阵元1的距离越远,其感应电流IM1就越小。因此,在仿真有源阵元方向图时,不必仿真整个阵列,只需仿真由阵元1附近的辐射阵元组成的子阵列即可.根据所需精度大小,选取子阵列的大小。
修正算法细节:
导引头的口径面场分布式:
利用电磁仿真可以得到导引头口径面上若干个离散点处的电场,于是可以如公式(3)将计算平均电场分布的积分离散化,避免推导繁复的解析公式。N值越大,计算结果越精确。
改进后的算法流程图如图7所示。
Claims (2)
1.近场效应误差分析方法,其包括以下步骤:
步骤一:
(1)、建立坐标系:三个辐射阵元呈等边三角形分布在二维直角坐标系x0y平面内,上述三个辐射阵元的中心位于坐标原点O;三个辐射阵元在直角坐标系x0y中的坐标为(xi,yi),i=1、2、3;固定观测点P在以O为坐标原点的三维坐标系中的位置为(0,0,R);任意观测点P′在以O点为球心的球坐标系2中的位置为(α,β);任意观测点P′在以O为坐标原点的三维坐标系中的位置为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
<mo>=</mo>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>=</mo>
<mi>R</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&beta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>z</mi>
<mo>=</mo>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
<mi>cos</mi>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
(2)、将观测点从固定观测点P运动到任意观测点P′,得到观测点到三个辐射阵元的距离变化量ri为:
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&alpha;</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mi> </mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&beta;</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mi> </mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&beta;</mi>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<msqrt>
<mrow>
<msup>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>;</mo>
</mrow>
(3)、定义三个阵元辐射的电磁波角频率为ω,初始相位为φi,在任意观测点P′处接收到三个辐射阵元的辐射信号为:
Ei=Aiexp[j(ωt-kR-kri+φi)],
其中,Ai(i=1,2,3)为三个辐射阵元的馈电幅度;j为虚数单位,j的平方后的值为-1;t是时间变量;k为相位常数,且k=2π/λ,λ表示三阵元辐射电磁波的波长;φi(i=1,2,3)为三个辐射阵元馈电的初始相位;
(4)、将三个阵元等效为阵面上的一个辐射源,任意观测点P′处接收的辐射信号E为:
E=Aexp[j(ωt-kR+δ)],
其中,A为三阵元等效辐射源的馈电幅度;k为相位常数,且k=2π/λ;R为三阵元中心O点到固定观测点P的距离;δ为在任意观测点P′处三阵元辐射信号合成后的等相位面相对于在三阵元中心O处的参考源辐射的等相位面的倾斜;
(5)、改变三阵元的馈电幅度和相位中的一种或者两种,从而控制在任意观察点P′处测得的等相位面的位置,所述的等相位面的位置与馈电幅度和相位的关系如下式:
其中,n=3;
(6)、改变三阵元的馈电幅度和相位中的一种或者两种,从而控制目标的位置,所述的目标的位置与馈电幅度和相位的关系如下式:
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>cos&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>cos&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>+</mo>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>sin&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>&theta;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>sin&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>sin&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msub>
<mi>cos&phi;</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,(i=1,2,3)是三个辐射阵元在以P点为球心的球坐标系1中的坐标,是三个辐射阵元的等效辐射源在以P点为球心的球坐标系1中的坐标,n=3;
当φi=0(i=1,2,3)时,上式可简化为:
<mrow>
<mi>&theta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>E</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
</mrow>
<mrow>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>E</mi>
<mn>3</mn>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
步骤二:
三阵元辐射阵元i的馈电幅度Ai(i=1,2,3)和相位φi(i=1,2,3),得到辐射阵元i在接收天线上接收到的平均电场m表示近场效应修正算法流程中的迭代序列,其中是第m次迭代中天线阵元i在接收天线口径面上产生的场强,Es(ρ,ζ)是接收天线的口径面场分布,S口径表示天线口径面的面积,∫∫口径面表示对口径面上的积分计算;
步骤三:
令m=m+1,导引头的口径面场分布为:
利用电磁仿真可以得到导引头口径面上若干个离散点处的电场,将计算平均电场分布的积分离散化,避免推导繁复的解析公式,N值越大,计算结果越精确。
2.根据权利要求1中所述的近场效应误差分析方法,其特征在于:分析阵元间残留互耦,具有一定口径的接收天线,辐射单元的间距、工作波长、接收天线口径的尺寸以及天线的形式、方向图因素,得到误差产生的原因和怎样消去或者减小误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410293567.2A CN104063544B (zh) | 2014-06-25 | 2014-06-25 | 近场效应误差分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410293567.2A CN104063544B (zh) | 2014-06-25 | 2014-06-25 | 近场效应误差分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104063544A CN104063544A (zh) | 2014-09-24 |
CN104063544B true CN104063544B (zh) | 2017-11-24 |
Family
ID=51551257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410293567.2A Active CN104063544B (zh) | 2014-06-25 | 2014-06-25 | 近场效应误差分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104063544B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536291B (zh) * | 2014-11-25 | 2017-03-29 | 北京理工大学 | 基于射频系统模拟弹性振动对导引头测量信号影响的方法 |
CN106123889B (zh) * | 2016-06-14 | 2019-03-01 | 中国科学院上海天文台 | 适用于近场目标成像的空间坐标系建立方法 |
CN106209269B (zh) * | 2016-07-27 | 2018-09-28 | 华东师范大学 | 一种射频仿真系统中球面复合阵列近场效应的校准方法 |
CN107239044B (zh) * | 2017-06-05 | 2020-06-02 | 北京仿真中心 | 一种限定性求解矢量控制近场效应修正初值的方法 |
CN107783086B (zh) * | 2017-10-17 | 2021-01-05 | 中国电子科技集团公司第三十八研究所 | 用于诊断天线阵口径幅相场的畸变位置的方法 |
CN108959806B (zh) * | 2018-07-23 | 2022-03-15 | 电子科技大学 | 一种基于球面近场测量和球模式源的等效辐射建模方法 |
CN111562445B (zh) * | 2020-05-06 | 2022-08-19 | 中国人民解放军63892部队 | 一种射频仿真试验系统角模拟精度实时监测的方法 |
CN111664816B (zh) * | 2020-07-13 | 2021-03-23 | 上海交通大学 | Qtt天线主反射面变形检测方法、校正方法及检测装置 |
CN113486550B (zh) * | 2021-06-23 | 2023-07-14 | 北京科技大学 | 一种半实物射频仿真中的三元组馈电系数确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103022699A (zh) * | 2012-12-27 | 2013-04-03 | 北京理工大学 | 星载反射面稀疏相控阵多波束天线 |
CN103682677A (zh) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | 一种艇载雷达共形稀疏阵列天线及其信号处理方法 |
-
2014
- 2014-06-25 CN CN201410293567.2A patent/CN104063544B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103022699A (zh) * | 2012-12-27 | 2013-04-03 | 北京理工大学 | 星载反射面稀疏相控阵多波束天线 |
CN103682677A (zh) * | 2013-11-14 | 2014-03-26 | 中国科学院电子学研究所 | 一种艇载雷达共形稀疏阵列天线及其信号处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104063544A (zh) | 2014-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104063544B (zh) | 近场效应误差分析方法 | |
US7119739B1 (en) | Near field to far field DF antenna array calibration technique | |
CN106209269B (zh) | 一种射频仿真系统中球面复合阵列近场效应的校准方法 | |
CN108445303B (zh) | 一种近场电磁散射特性模拟方法 | |
CN104391183A (zh) | 近场测量天线远区场特性快速计算方法 | |
CN103777186B (zh) | 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 | |
CN105606906B (zh) | 一种毫米波相控阵测试标定方法 | |
CN103364645A (zh) | 虚拟馈电网络的天线阵列近场测量方法 | |
US11131701B1 (en) | Multi-probe anechoic chamber for beam performance testing of an active electronically steered array antenna | |
CN104375129B (zh) | 一种分布式阵列相参合成雷达发射相参参数校准方法 | |
CN108959806A (zh) | 一种基于球面近场测量和球模式源的等效辐射建模方法 | |
Zeng et al. | Low side lobe pattern synthesis using projection method with genetic algorithm for truncated cone conformal phased arrays | |
CN107942145A (zh) | 一种角变化球面三元组天线辐射电磁场测量系统和方法 | |
Meng et al. | Robotic radiation pattern measurement system for 6–110 GHz based in both near field and far field | |
CN116593981A (zh) | 一种单发多收柱面波圆周稀疏阵局域扫描的小角度rcs精准测量方法 | |
CN107015065B (zh) | 窄波束天线电轴、相位中心和时延的远场联合标定方法 | |
CN111965602B (zh) | 一种相控阵雷达幅相一致性检测方法和系统 | |
Geise et al. | On the Imaging of Large Antenna Array Navigation Systems | |
CN104133376B (zh) | 一种用于射频仿真中目标的矢量控制方法 | |
CN113252999A (zh) | 天线平面近场测试方法 | |
Lu et al. | Modeling and simulation of a tri-reflector CATR system at 300GHz | |
Li et al. | Investigation on method of improving precision for radio frequency simulation system | |
Guo et al. | Simulation of Electromagnetic Target Based on Movement of Quite Zone | |
CN113486550B (zh) | 一种半实物射频仿真中的三元组馈电系数确定方法 | |
Fu et al. | An imaging algorithm for burden surface with T-shaped MIMO radar in the blast furnace |
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 |