一种基于散射矩阵原理的球面近远场方法与流程

1.本发明涉及天线测量技术领域，具体涉及一种基于散射矩阵原理的球面近远场方法。


背景技术：

2.在当下，想要提高天线设计理念和改善天线分析的理论方法就必须要提高天线辐射的测量精确度，这个必要性已经体现在各种各样的领域里，例如在特殊的太空应用里，对天线的各项性能和规格都有严格的要求，这就需要对天线进行深入的分析，而天线测量正好为此提供了一种强有力的手段。因此，在近几年里天线测量已经引起了人们极大的重视，世界各国的学者们都做了大量的研究，提出了新的发展方向，研究了不同的测量方法，构建了不同种类的测量系统。
3.随着对天线的要求越来越高，人们急需找到一种测量精度高、能适应各种环境并且能提高测量效率的测量方法，近场测量法由此诞生。近场测量是一种间接方法，它是在微波暗室内，把待测天线放在转台上，使用一个小探头在待测天线的近场测量面的规则网格上进行扫描采样通过对扫描得到的近场数据进行合适的处理，就可获得待测天线的远场特性。


技术实现要素：

4.发明目的：本发明旨在提出一种由近场测量数据外推待测天线的远场特性的方法，实现球面近场远场的变换以及得到近场实测数据外推远区方向图。
5.技术方案：为实现上述技术效果，本发明提出的技术方案为：
6.由近场测量数据外推待测天线的远场特性的方法，包括依次执行的以下步骤：
7.步骤1，利用待测天线的散射矩阵推导探头的接收信号与待测天线的激励信号之间的传输方程。
8.步骤1.1，利用赫姆霍兹方程推导出球面波函数。
9.步骤1.2，推导出待测天线的散射矩阵。
10.步骤1.3，将待测天线远场的电场和磁场采用球面波函数进行表示，得到基于球面波函数的待测天线远场表达式。
11.步骤1.4，根据天线的散射矩阵以及基于球面波函数的待测天线远场表达式，推导出探头的接收信号与待测天线的激励信号之间的传输方程。
12.步骤2，将探头的接收信号与待测天线的激励信号之间的传输方程代入到天线传输系数方程，通过解天线传输系数方程得到待测天线传输系数vt
1mn
和vt
2mn
。
13.步骤3，把传输系数vt
1mn
和vt
2mn
代入远场的传输公式，求出待测天线的远场表达式并绘制近远场方向图。
14.进一步的，还包括探头补偿，具体过程包括推导出探针响应常数公式并且进行探头响应补偿修正，可消除探头本身给测量数据带来的误差。
15.有益效果：本发明对球面近场测量的近远场变换算法进行了研究，与传统的基于模式展开的近远场变换算法不同，本发明采用的是散射矩阵来推导近远场变换算法。该方法是基于微波网络的思想，把待测天线和探头分别看成一个开放的二端口网络，通过探头的接收信号与待测天线的激励信号之间的传输方程，把待测天线的近场数据代入传输公式，求出待测天线的传输系数。由于近场和远场的传输系数是相同的，所以把求得的传输系数代入远场的传输公式，就可以求出待测天线的远场表达式。因此，本算法的关键是求出待测天线的传输系数。
附图说明
16.图1是本发明的一种球面近场测量外推远场的步骤流程图；
17.图2是球面近场测量外推远场计算流程框图；
18.图3是待测天线场分布图；
19.图4是待测天线的波导与球面波模式图；
20.图5是原始近场数据图；
21.图6是远场外推结果图；
22.图7是反推远场数据与实测近场数据对比图。
具体实施方式
23.目前为止，我们无法直接测量待测天线的远场特性，故而我们采取了一种间接测量的方式：近场测量外推远场数据。近场测量是一种间接方法，它是在微波暗室内，把待测天线放在转台上，使用一个小探头在待测天线的近场测量面的规则网格上进行扫描采样通过对扫描得到的近场数据进行处理，就可获得待测天线的远场特性。
24.图1所示为本发明的一种球面近场测量外推远场的整体流程图，包括以下步骤：
25.步骤1，利用待测天线的散射矩阵推导探头的接收信号与待测天线的激励信号之间的传输方程；
26.步骤1.1，利用无源区的赫姆霍兹方程推导出球面波函数
27.如图4，r＝r0是包围待测天线的一个最小球面半径。则在r≥r0的区域，待测天线产生的场是满足无源区的赫姆霍兹方程的:
[0028][0029][0030]
其中是电场矢量，是磁场矢量。
[0031]
若f(r,θ,φ)是标量赫姆赫兹方程在球坐标系中的解，也就是:
[0032]
则可由
[0033]
所组成的两组矢量函数构成赫姆霍兹方程的解：
[0034][0035]
其中表示矢量波函数。
[0036]
式的基本解为：
[0037][0038]
其中是归一化连带勒让德函数,其中
[0039][0040][0040]
表示径向函数，其具体表达式如下：
[0041]
z
n(1)
＝j
n
(kr)(球贝塞尔函数)
[0042]
z
n(2)
＝n
n
(kr)(球诺依曼函数)
[0043]
z
n(3)
＝h
n(1)
(kr)＝j
n
(kr) i*n
n
(kr)(第一类球汉克尔函数)
[0044]
z
n(4)
＝h
n(2)
(kr)＝j
n
(kr)
‑
i*n
n
(kr)(第二类球汉克尔函数)
[0045]
c＝1和c＝2，j
n
(kr)及n
n
(kr)代表驻波，而c＝4,h
n(1)
(kr)代表出射行波，c＝3,h
n(2)
(kr)代表入射行波。在r＝0处n
n
(kr)具有奇异性，因此，对于0≤r≤r0的有限区域的内的场选取j
n
(kr)为波函数，对于r0≤r≤∞的外部区域选取h
n
(kr)为波函数。此时，矢量波函数则变为:
[0046][0047]
矢量波函数变为：
[0048]
[0049]
为球面波函数，是归一化连带勒让德函数,
[0050]
其中
[0051][0051]
表示径向函数。
[0052]
步骤1.2，推导出待测天线的散射矩阵
[0053]
在图4中待测天线最小球面外的总电场为：
[0054][0055]
其中，j＝2{n(n 1) m
‑
1} s，1} s，n、m表示了不同主模模式的矩形波导，n＝kr0 10为最高模式阶数,r0为包围待测天线的最小球面半径。
[0056]
由于阶数n的限制，球面波的模式数有限，传输线模型数也将有限，则出射波和入射波的线性关系可以写成j 1次矩阵:
[0057][0058]
式中γ表示天线反射系数，是复数，r表示天线接收系数，是一个元素为jr的1
×
j的行矩阵，t表示天线传输系数，是一个元素为jt的j
×
1的列矩阵，j＝1,2,
……
,j，s表示天线散射系数，是元素为s
ij
的方阵。列矩阵a和b分别为接收端口处的入射球面波和出射球面波的波系数，j
×
1的列矩阵a和b元素分别为a1、a2、
……
a
j
和b1、b2……
b
j
。v和w分别为发射天线端口的输入波和输出波的波系数，为待测天线的散射矩阵。
[0059]
步骤1.3，将待测天线远场的电场和磁场采用如下球面波函数进行表示，待测天线场分布如图3所示。
[0060][0061][0062]
其中，以及为引入符号，s＝1时表示函数，s＝2时表示函数；k表示波数；η表示波阻抗。
[0063]
引入定义远场方向图函数
[0064]
得到基于球面波函数的待测天线远场的电场和磁场表达式变为：
[0065]
[0066][0067]
v表示激励的幅度，是绝对远场图，定义为:
[0068][0069]
t
smn
表示天线传输系数，表示方向图函数。
[0070]
步骤1.4，根据天线的散射矩阵以及待测天线远场表达式，推导出探头的接收信号与待测天线的激励信号之间的传输方程
[0071]
定义a为探针坐标系原点到原坐标系原点的距离；
[0072]
θ为探针坐标系原点到原坐标系原点连线与z轴相交的角度；
[0073]
φ为探针坐标系原点到原坐标系原点连线在xoy平面投影与x轴相交的角度。
[0074]
如图2所示，输入探头采样数据待测天线的激励信号w1(a，θ，φ)和探头的接收信号w
‑1(a，θ，φ)，根据两次数据变换计算出w
1m
(a，θ)和w
‑
1m
(a，θ)，所述两次数据变换的过程就是对激励信号w1(a，θ，φ)和探头的接收信号w
‑1(a，θ，φ)降维的过程。
[0075]
具体为使用到方程组：
[0076]
一次变换数据：
[0077][0078]
二次变数数据：
[0079][0080]
三次数据变换推导出和
[0081][0082]
其中为旋转系数：
[0083][0084]
步骤2、计算待测天线的传输系数vt
1mn
和vt
2mn
[0085]
根据待测天线的激励信号和探头的接收信号之间的传输方程，利用探头在待测天线近场的采样数据，计算出待测天线的传输系数vt
1mn
和vt
2mn
，其中v表示激励的幅度，进而
求出待测天线的远场方向图。
[0086]
近场测量不同于远场测量，远场测量要求待测天线与测量探头之间的距离必须大于2d2/λ，其中d表示探头的口径，λ表示电磁波波长；而球面近场测量则要求测量探头距离包围待测天线的最小球面3～5个波长即可。
[0087]
测量时，待测天线固定在(x,y,z)坐标系的原点上，而测量探头位于它自身的(x’,y’,z’)坐标系里，探头将在以(x,y,z)坐标系的原点为中心、a为半径的球面上移动。因此，在待测天线坐标系(x,y,z)里，探头的坐标就是(a,θ,ф)。当探头在测量球面上移动时，它会一直指向(x,y,z)坐标系的原点，并且它还需要围绕本身的轴旋转一个角度χ。
[0088]
通过解解天线传输系数方程可得到待测天线传输系数vt
1mn
和vt
2mn
。
[0089]
通过解天线传输系数方程：
[0090][0091]
可以得到待测天线传输系数，如下式所示。将步骤1.4求得待测天线的激励信号和探头的接收信号之间的传输方程，代入天线传输系数方程，可得到天线传输系数。
[0092][0093][0094]
步骤3、探头补偿
[0095]
在平面近场测量中，用探头在扫描面上进行数据采集，作为采样探头的开口矩形波导探头对测量结果有直接的影响，由于探头具有一定的大小和方向性，它的存在会对待测场有扰动，从而导致近远场变换后的数据不仅包含了待测天线的参数信息，也包含了采样探头的参数信息。要得到准确的待测天线参数信息，必须将探头信息从近远场变换后的采样数据中剔除掉。
[0096]
步骤3.1，推导探针响应常数公式
[0097]
对于μ＝
±
1的线极化探针，利用平移系数和球面波展开系数、探针响应常数的关系可得到球面te波(s＝1)和tm波(s＝2)的探针响应常数，具体的：
[0098]
平移系数：
[0099]
其中表示s＝1、1次方的平移系数，表示s＝2、2次方的平移系数，表示s＝2、1次方的平移系数，表示s＝1、2次方的平移系数.
[0100]
球面波展开系数：q
1，
‑
1，v
＝
‑
q
11v
，q
2，
‑
1，v
＝
‑
q
21v
；
[0101]
其中，q
1，
‑
1，v
表示s＝1、μ＝
‑
1时球面波展开系数，q
11v
表示s＝1、μ＝1时的球面波展开系数，q
2，
‑
1，v
表示s＝2、μ＝
‑
1时的球面波展开系数，q
21v
表示s＝1、μ＝1时的球面波展开系数。
[0102]
探针响应常数：p
s，
‑
1，n
(ka)＝(
‑
1)
s 1
p
s1n
(ka)。
[0103]
其中，p
s，
‑
1，n
表示μ＝
‑
1时的探针响应常数，p
s1n
表示μ＝1时的探针响应常数。
[0104]
可得到球面te波(s＝1)和tm波(s＝2)的探针响应常数公式如下：
[0105][0106]
若探针为圆锥喇叭、开口圆波导和矩形波导，只需测得探针远区e面和h面方向图，则上式q
11v
、q
21v
为探针的球面波展开系数，可利用非探针校准方法解：
[0107][0108]
其中表示探针为圆锥喇叭、开口圆波导和矩形波导时的探针响应常数，是归一化连带勒让德函数，其中是归一化连带勒让德函数，其中
[0109]
步骤3.2，探头响应补偿修正
[0110]
测得探头辐射方向图带入到探头响应常数公式中，得到修正后的探头平移系数：
[0111][0112][0113]
其中，相比于矩形开口波导天线，vivaldi天线的辐射特性复杂，故在计算其探头响应系数时需要更多的辐射模式，一般取i＝k r0 10，r0为包围探针口径的最小球面半径。
[0114]
探头信息会影响实际测量数据，即探头发射、接受的数据。对探头响应补偿修正本发明中是指对探头平移系数修正。
[0115]
步骤4、把传输系数代入远场的传输公式，求出待测天线的远场表达式。
[0116]
将探头在待测天线近场的采样数据(现场测量数据)带入远场传输公式：
[0117][0118]
其中，t
smn
表示天线未知的传输系数，表示旋转系数，p
sμn
表示探头接收数据，v表示激励信号。
[0119]
可由vt
1mn
和vt
2mn
计算出待测天线远场在φ方向上的表达式w
φ
(θ，φ)和待测天线远场在θ方向上的表达式w
θ
(θ，φ)，并绘出近远场方向图，如图5和图6所示，并通过数据绘制反推远场数据与实测近场数据对比图，如图7所示。
[0120]
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。