CN109818689B - 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 - Google Patents
一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 Download PDFInfo
- Publication number
- CN109818689B CN109818689B CN201711168211.6A CN201711168211A CN109818689B CN 109818689 B CN109818689 B CN 109818689B CN 201711168211 A CN201711168211 A CN 201711168211A CN 109818689 B CN109818689 B CN 109818689B
- Authority
- CN
- China
- Prior art keywords
- array
- antenna
- port
- array element
- calibration
- 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
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000003860 storage Methods 0.000 title claims description 6
- 238000010586 diagram Methods 0.000 claims abstract description 71
- 239000011159 matrix material Substances 0.000 claims abstract description 69
- 238000005259 measurement Methods 0.000 claims abstract description 64
- 230000005284 excitation Effects 0.000 claims abstract description 59
- 238000010168 coupling process Methods 0.000 claims description 74
- 238000005859 coupling reaction Methods 0.000 claims description 74
- 230000008878 coupling Effects 0.000 claims description 69
- 238000004088 simulation Methods 0.000 claims description 16
- 230000009466 transformation Effects 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 238000003491 array Methods 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 239000000523 sample Substances 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
Images
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种阵列天线的校准方法,包括如下步骤:S1.获得阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;S2.馈入T组线性无关的端口激励I1,I2…IT,获得相应的T组口径场激励I′1,I′2…I′T;S3.根据T组线性无关的端口激励I1,I2…IT和T组口径场激励I′1,I′2…I′T计算得到校准矩阵C,其中对于任意一组对应的端口激励It和口径场激励I′t都有I′t=C×It;S4.根据校准矩阵C对阵列天线各阵元进行校准。本发明通过较少的测量数据,结合阵列天线已知的先验知识,能够快速、高效地对阵列天线进行校准。与现有的校准方法相比,本发明的方法测量数据少,效率高,且系统组成简单,易于工程实现,成本低。
Description
技术领域
本发明涉及天线技术领域,具体而言,涉及一种阵列天线的校准方法、设备、系统以及计算机可读存储介质。
背景技术
天线在通信、广播、电视、雷达和导航等无线电系统中被广泛的应用,起到了传播无线电波的作用,是有效地辐射和接受无线电波必不可少的装置。
阵列天线是一类由不少于两个天线阵元规则或随机排列并通过适当激励获得预定辐射特性的天线。近年来,阵列天线作为民用和军用天线技术的一个重要发展方向备受关注。
阵列天线由很多天线阵元组成,每个阵元馈以一定幅度和相位的信号以形成特定的波束并实现波束扫描,各阵元的信号叠加形成阵列天线的信号。一般情况下,通过控制连接阵元的衰减器对阵元的信号幅度进行调节并形成所需要的波束,通过改变连接阵元的移相器的相位对阵元信号的相位进行控制以实现波束扫描。
阵列天线在实际加工制造过程中,由于加工精度等造成的结构不对称,以及器件本身的不一致性,另外天线本身的波动、天线阵元之间的互耦等使得每一个天线阵元的幅度和相位可能与预期值不同,从而会使阵列性能达不到最佳工作状态,甚至达不到所要求的技术条件,影响通信质量。因此,需要对阵列天线进行校准,目的就是尽量消除误差,使阵列性能达到所要求的技术条件或最佳工作状态。
目前阵列天线的校准方法主要包括以下几种:
1.近场测量法:该方法通过探头对围绕阵列闭合面上电场的采样,经过数值计算获得远场方向图和阵列幅相分布。该方法测量精度高,可用于各种体制的阵列天线,但对测量仪器的同步性要求高、扫描时间长,数据量大,测量效率低。
2.旋转矢量法,该方法只要测得合信号功率随天线阵元相位变化的正弦曲线,就可计算出每个阵元的幅相值,而无需矢量测量仪器。但在实际应用中,若阵列阵元数较多,则单个阵元造成的曲线变化不明显,且需要校准的时间也较长。
3.互相校准法,该方法是基于大型阵列的阵中相邻阵元的互相系数是相同的这一基本原理,通过对阵列中相邻阵元进行收发测试,由测试数据计算出各有源通道的幅相信息,再根据理想分布进行阵列校准。互相校准法虽然无需外场测量装置,测试过程简单,但其仅适用于收发共口面的相控阵天线,且辐射阵元之间的隔离度不能太大。
4.换相测量法,该方法的基本思想是在相控阵天线和测量探头均保持固定的情况下,测量相控阵天线不同配相状态下探头的接收信号幅相,然后用数学算法对实验数据进行处理即可确定任意配相状态下各通道激励的幅相,进而根据一次试验结果复原所有的方向图。换相测量法的不足之处在于,一般情况下系数矩阵不满秩,在解测量方程时,必须引入其它的先验知识,这样也会导致其测量时间较长,测试效率不高。
发明内容
本发明的主要目的在于克服现有技术的不足,提供一种阵列天线的校准方法,通过较少的测量,快速、精确地对阵列天线进行校准。
为实现上述目的,本发明一方面提出了一种阵列天线的校准方法,其特征在于,所述阵列天线包括N个阵元,所述校准方法包括如下步骤:
S1.获得所述阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;
S2.馈入T组线性无关的端口激励I1,I2…IT,获得相应的T组口径场激励I′1,I′2…I′T,具体包括如下步骤:
S201.馈入端口激励It,t=1,2,…,T,
S202.获得Mt个测量点的位置以及阵列天线在Mt个测量点处的电/磁场的测量数据Et,所述测量数据Et包含幅度和相位信息,Mt≥N/3;
S203.根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据Et获得口径场激励I′t;
S3.根据T组线性无关的端口激励I1,I2…IT和T组口径场激励I′1,I′2…I′T计算得到校准矩阵C,其中对于任意一组对应的端口激励It和口径场激励I′t都有I′t=C×It;
S4.根据所述校准矩阵C对阵列天线各阵元进行校准。
作为本发明的进一步限定,所述步骤S2中,T≥N,构建矩阵II,II=[I1 I2 … IT],构建矩阵II′,II′=[I′1 I′2 … I′T],所述步骤S3中根据II′=C×II计算得到校准矩阵C。
作为本发明的进一步限定,校准矩阵C表达为
其中,Cn=[C1n C2n … CNn]T为校准矩阵C中的第n列,称为第n个阵元的端口耦合系数,表示第n个阵元端口对各阵元端口的耦合,Cdn称为第n个阵元端口对第d个阵元端口的耦合元素,n=1,2,…,N,d=1,2,…,N,()T表示转置;
定义以单个阵元为圆心,xλ为半径的圆为该单个阵元的耦合区,其中x为不小于1的实数,λ为阵列天线工作频率的波长,对任意两个阵元而言,如果其耦合区内的阵元数量和分布相同,则所述两个阵元互为等耦合阵元,令所述阵列天线包括U组互不相同的等耦合阵元,在所述U组互不相同的等耦合阵元的每组中分别任意选取1个阵元组成包括U个阵元的子阵列,所述步骤S2中,向所述子阵列馈入T组线性无关的端口激励I1,I2…IT,T≥U,获得相应的所述子阵列的T组口径场激励I′1,I′2…I′T,构建矩阵III,III=[I1 I2 … IT],构建矩阵III′,III′=[I′1 … I′T],所述步骤S3包括如下两个步骤:
S301.根据III′=C′×III计算得到所述子阵列的校准矩阵C′,
S302.根据校准矩阵C′中的端口耦合系数以及各阵元的相对位置,计算得到校准矩阵C。
作为本发明的进一步限定,所述步骤S302包括如下步骤:
S3021.在所述包括U个阵元的子阵列中任意指定一个阵元A,阵元A的耦合区为A′,提取阵元A的端口耦合系数CA中阵元A端口对耦合区A′中各阵元端口的耦合元素;
S3022.在所述阵元A的等耦合阵元中任意指定一个阵元B,阵元B的耦合区为B′,阵元B端口对耦合区B′中任意一个阵元端口的耦合元素与阵元A端口对耦合区A′中同样相对位置的阵元端口的耦合元素相等,据此可确定阵元B的端口耦合系数CB中阵元B端口对耦合区B′中的阵元端口的耦合元素,端口耦合系数CB中阵元B端口对耦合区B′之外的阵元端口的耦合元素为0,其中,所述“同样相对位置”指耦合区中的阵元相对于该耦合区圆心的位置。
作为本发明的进一步限定,所述阵中方向图通过测量获得,或者基于阵列天线的物理参数或/和机械模型或/和仿真模型通过仿真获得,所述阵列天线的物理参数包括天线形式和阵列结构。
作为本发明的进一步限定,所述任意两个等耦合阵元的阵中方向图相同。
作为本发明的进一步限定,所述步骤S203中阵中方向图、阵中方向图中心的位置、测量点的位置、测量数据Et以及口径场激励I′t满足关系式:Et=YI′t,其中Et是所述M个测量点测得的电/磁场,为M×1的矩阵,Y是阵元到测量点的幅相变换矩阵,Y根据阵中方向图、阵中方向图中心的位置以及测量点的位置获得。
作为本发明的进一步限定,以任意参考点为原点建立球坐标系,第n个阵元的阵中方向图中心的位置的坐标为(Rn,θn,φn),n=1,2,…,N,第n个阵元的阵中方向图表示为fn(θ,φ),第m个测量点的位置的坐标为(R′m,θ′m,φ′m),m=1,2,…,M,所述阵元到测量点的幅相变换矩阵Y为
是第n个阵元的阵中方向图在第m个测量点的位置的幅相变换因子,其中,(θ′mn,φ′mn)是第m个测量点的位置相对于第n个阵元的阵中方向图中心的位置所在的角度,fn(θ′mn,φ′mn)是第n个阵元在(θ′mn,φ′mn)角度的阵中方向图信息,包含幅度和相位信息,是对第n个阵元的阵中方向图在第m个测量点的位置进行的相位修正,是第m个测量点的位置指向第n个阵元的阵中方向图中心的位置的矢量的模长,k是电磁波传播常数。
作为本发明的进一步限定,所述阵列天线各阵元的阵中方向图相同,f1(θ,φ)=f2(θ,φ)=…=fN(θ,φ)=f(θ,φ),所述阵元到测量点的幅相变换矩阵Y为
作为本发明的进一步限定,所述测量点位于阵元的辐射远场。
作为本发明的进一步限定,当Mt>N/3时,通过最小二乘法计算得到所述口径场激励I′t;当T>N时,通过最小二乘法计算得到所述校准矩阵C;当T>U时,通过最小二乘法计算得到所述校准矩阵C′。
本发明另一方面提出了一种阵列天线的校准设备,其特征在于,所述校准设备包括:
阵中方向图获取模块,用于获取所述阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;
馈电模块,用于向所述阵列天线馈入端口激励It;
信号收发模块,与所述测量天线连接,用于获得Mt个测量点的位置、通过测量天线在Mt个测量点处向所述阵列天线发射测量信号并获得所述阵列天线电/磁场的测量数据,所述测量数据包含幅度和相位信息,Mt≥N/3,N为所述阵列天线的阵元数;
口径场激励获取模块,用于根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据获得口径场激励I′t;
校准矩阵获取模块,用于根据端口激励It与口径场激励I′t的关系I′t=C×It计算得到校准矩阵C;
校准模块,用于根据所述校准矩阵C对阵列天线各阵元进行校准。
作为本发明的进一步限定,所述阵中方向图获取模块包括:
测量单元,用于测量获得所述阵列天线的阵中方向图;或/和
仿真单元,用于基于阵列天线的物理参数或/和机械模型或/和仿真模型通过仿真获得所述阵列天线的阵中方向图,所述物理参数包括天线形式和阵列结构。
本发明另一方面提出了一种阵列天线的校准设备,其特征在于,所述校准设备包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现前述方法的步骤。
本发明另一方面提出了一种阵列天线的校准系统,包括电波暗室以及测量天线,其特征在于,在所述校准系统中集成前述的校准设备。
本发明另一方面提出了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现前述方法的步骤。
本发明通过较少的电/磁场测量数据,结合测量点的位置信息、阵列天线的阵元的阵中方向图及阵中方向图中心的位置信息等阵列天线的先验知识,反演获得口径场激励I′,再计算得到校准矩阵,对阵列天线各阵元进行校准。与现有的校准方法相比,本发明的方法测量数据少,测试效率高,且系统组成简单,易于工程实现,成本低。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明,其中:
图1为本发明第一实施方式的阵列天线的校准方法的流程图。
图2为本发明第二实施方式的阵列天线的校准设备的框图。
具体实施方式
以下对本发明的实施方式作详细说明。应该强调的是,下述说明仅仅是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
此外,在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
由N个阵元构成的阵列天线,其阵列天线的方向图可认为是在口径场激励下,所有阵元的方向图的迭加。本发明通过阵列天线的电/磁场的测量数据、测量点的位置信息、阵列天线的阵元的阵中方向图以及阵中方向图中心的位置,反演得到口径场激励,再根据端口激励和口径场激励的数学关系获得校准矩阵,据此对阵列天线各阵元进行校准。图1示意了本发明第一实施方式的阵列天线的校准方法流程,所述方法包括如下步骤:
S1.获得阵列天线的阵元的阵中方向图以及阵中方向图中心的位置。所述阵中方向图通过测量获得,或者基于阵列天线的物理参数(包括天线形式和阵列结构)或/和机械模型或/和仿真模型通过仿真获得。
S2.馈入T组线性无关的端口激励I1,I2…IT,获得相应的T组口径场激励I′1,I′2…I′T,具体包括如下步骤:
S201.馈入端口激励It,t=1,2,…,T;端口激励It是已知的,定义端口激励It如下:
S202.在Mt个测量点处通过测量天线对阵列天线进行辐射测量,获得Mt个测量点的位置以及Mt个测量点处的电/磁场的测量数据Et,所述测量数据E包含幅度和相位信息,Mt≥N/3。这里可采用多种测量方式,如常见的球面扫描、平面扫描、柱面扫描等,或其他测量方式;测量点位于阵元的辐射远场。
S203.根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据Et获得口径场激励I′t;定义口径场激励I′t如下:
本步骤中阵中方向图、阵中方向图中心的位置、测量点的位置、测量数据Et以及口径场激励I′t满足关系式:Et=YI′t,其中Et是所述M个测量点测得的电/磁场,为M×1的矩阵,Y是阵元到测量点的幅相变换矩阵,Y根据阵中方向图、阵中方向图中心的位置以及测量点的位置获得。
S3.根据T组线性无关的端口激励I1,I2…IT和T组口径场激励I′1,I′2…I′T计算得到校准矩阵C,其中对于任意一组对应的端口激励It和口径场激励I′t都有I′t=C×It;
S4.根据所述校准矩阵C对阵列天线各阵元进行校准。
下面具体阐述本实施方式中获得校准矩阵C的两种方式。
第一种方法中,T≥N,构建矩阵II,II=[I1 I2 … IT],构建矩阵II′,II′=[I′1I′2 … I′T],所述步骤S3中根据II′=C×II计算得到校准矩阵C。本方法中,至少需要进行N组测量,测量数据较多,测量时间较长。
第二种方法测量数据少,测量时间短,具体实施方式为:
校准矩阵C表达为
其中,Cn=[C1n C2n … CNn]T为校准矩阵C中的第n列,称为第n个阵元的端口耦合系数,表示第n个阵元端口对各阵元端口的耦合,Cdn称为第n个阵元端口对第d个阵元端口的耦合元素,n=1,2,…,N,d=1,2,…,N,()T表示转置;
定义以单个阵元为圆心,xλ为半径的圆为该单个阵元的耦合区,其中x为不小于1的实数,λ为阵列天线工作频率的波长,对任意两个阵元而言,如果其耦合区内的阵元数量和分布相同,则所述两个阵元互为等耦合阵元,令所述阵列天线包括U组互不相同的等耦合阵元,在所述U组互不相同的等耦合阵元的每组中分别任意选取1个阵元组成包括U个阵元的子阵列,所述步骤S2中,向所述子阵列馈入T组线性无关的端口激励I1,I2…IT,T≥U,获得相应的所述子阵列的T组口径场激励I′1,I′2…I′T,构建矩阵III,III=[I1 I2 … IT],构建矩阵III′,III′=[I′1 … I′T],根据III′=C′×III计算得到所述子阵列的校准矩阵C′,在所述包括U个阵元的子阵列中任意指定一个阵元A,阵元A的耦合区为A′,提取阵元A的端口耦合系数CA中阵元A端口对耦合区A′中各阵元端口的耦合元素;在所述阵元A的等耦合阵元中任意指定一个阵元B,阵元B的耦合区为B′,阵元B端口对耦合区B′中任意一个阵元端口的耦合元素与阵元A端口对耦合区A′中同样相对位置的阵元端口的耦合元素相等,据此可确定阵元B的端口耦合系数CB中阵元B端口对耦合区B′中的阵元端口的耦合元素,端口耦合系数CB中阵元B端口对耦合区B′之外的阵元端口的耦合元素为0,其中,所述“同样相对位置”指耦合区中的阵元相对于该耦合区圆心的位置。由此计算得到校准矩阵C。
需要说明的是,以上所述实施顺序并非不可改变的。
下面具体阐述本实施方式中阵元到测量点的幅相变换矩阵Y的计算方法。
以任意参考点为原点建立球坐标系,第n个阵元的阵中方向图中心的位置的坐标为(Rn,θn,φn),n=1,2,…,N,第n个阵元的阵中方向图表示为fn(θ,φ),第m个测量点的位置的坐标为(R′m,θ′m,φ′m),m=1,2,…,M,所述阵元到测量点的幅相变换矩阵Y为
是第n个阵元的阵中方向图在第m个测量点的位置的幅相变换因子,其中,(θ′mn,φ′mn)是第m个测量点的位置相对于第n个阵元的阵中方向图中心的位置所在的角度,fn(θ′mn,φ′mn)是第n个阵元在(θ′mn,φ′mn)角度的阵中方向图信息,包含幅度和相位信息,是对第n个阵元的阵中方向图在第m个测量点的位置进行的相位修正,是第m个测量点的位置指向第n个阵元的阵中方向图中心的位置的矢量的模长,k是电磁波传播常数;
所述阵列天线各阵元的阵中方向图相同,f1(θ,φ)=f2(θ,φ)=…=fN(θ,φ)=f(θ,φ),所述阵元到测量点的幅相变换矩阵Y为
在此,本实施方式中有3点需要说明:
(1)本实施方式中,当任意两个阵元互为等耦合阵元时,认为所述两个阵元的阵中方向图相同。据此可以将耦合区内阵元数量和位置分布相同的阵元等效处理,对于阵元数目较多的阵列天线而言,这种等效处理可以大大降低阵元的测量次数或仿真计算量,大大提高测量速度。
(2)当Mt=N/3时,在计算口径场激励I′t时,方程数等于待求解的未知变量数,可以通过求解线性方程组得到口径场激励I′t;当Mt>N/3时,在计算口径场激励I′t时,方程数大于待求解的未知变量数,可以通过最小二乘法计算得到口径场激励I′t。同理,当T>N时,通过最小二乘法计算得到所述校准矩阵C;当T>U时,通过最小二乘法计算得到所述校准矩阵C′。
(3)本实施方式中所采用的球坐标系只是为了本发明的描述方便,本领域普通技术人员应当理解,也可采用其他坐标系进行描述,例如可按照公知的标准球坐标-直角坐标变换规则转换为直角坐标系进行描述,不影响本发明的实质内容,也应落入本发明的保护范围。
参照图2,本发明的第二实施方式为阵列天线的校准设备200,本实施方式中,校准设备200包括存储器201和处理器202,存储器201连接到处理器202,以便存储操作系统、应用、计算机程序代码、数据等,特别说明的是,存储器201存储有可在处理器202上运行的计算机程序,处理器202执行所述计算机程序时实现如前述第一实施方式所述方法的步骤,并且处理器202与以下模块相连:
阵中方向图获取模块203,用于获取所述阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;本模块包括:测量单元,用于测量获得所述阵列天线的阵中方向图;或/和仿真单元,用于基于阵列天线的物理参数(包括天线形式和阵列结构)或/和机械模型或/和仿真模型通过仿真获得所述阵列天线的阵中方向图;
馈电模块204,用于向所述阵列天线馈入端口激励It;
信号收发模块205,与所述测量天线连接,用于获得Mt个测量点的位置、通过测量天线在Mt个测量点处向所述阵列天线发射测量信号并获得所述阵列天线电/磁场的测量数据,所述测量数据包含幅度和相位信息,Mt≥N/3,N为所述阵列天线的阵元数;
口径场激励获取模块206,用于根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据获得口径场激励I′t;
校准矩阵获取模块207,用于根据端口激励It与口径场激励I′t的关系I′t=C×It计算得到校准矩阵C;
校准模块208,用于根据所述校准矩阵C对阵列天线各阵元进行校准。
需要说明的是,校准设备200是为了方便描述而示出的,校准设备200还可以包括其他必需的模块。此外,还可以将校准设备200中的至少部分模块进行合并或细分。
本发明的第三实施方式为阵列天线的校准系统,包括电波暗室以及测量天线,在所述校准系统中集成了如第二实施方式中所述的校准设备。
本发明的第四实施方式为一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如前述第一实施方式所述方法的步骤。
需要说明的是,可通过多种方式来实现本发明的实施方式,例如,硬件、固件、软件或其组合。
以上内容是结合具体/优选的实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,其还可以对这些已描述的实施方式做出若干替代或变型,而这些替代或变型方式都应当视为属于本发明的保护范围。
Claims (16)
1.一种阵列天线的校准方法,其特征在于,所述阵列天线包括N个阵元,所述校准方法包括如下步骤:
S1.获得所述阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;
S2.馈入T组线性无关的端口激励I1,I2...IT,获得相应的T组口径场激励I′1,I′2...I′T,具体包括如下步骤:
S201.馈入端口激励It,t=1,2,...,T,
S202.获得Mt个测量点的位置以及阵列天线在Mt个测量点处的电/磁场的测量数据Et,所述测量数据Et包含幅度和相位信息,Mt≥N/3;
S203.根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据Et获得口径场激励I′t;
S3.根据T组线性无关的端口激励I1,I2...IT和T组口径场激励I′1,I′2...I′T计算得到校准矩阵C,其中对于任意一组对应的端口激励It和口径场激励I′t都有I′t=C×It;
S4.根据所述校准矩阵C对阵列天线各阵元进行校准。
2.根据权利要求1所述的阵列天线的校准方法,其特征在于,所述步骤S2中,T≥N,构建矩阵II,II=[I1 I2 ... IT],构建矩阵II′,II′=[I′1 I′2 ... I′T],所述步骤S3中根据II′=C×II计算得到校准矩阵C。
3.根据权利要求1所述的阵列天线的校准方法,其特征在于,校准矩阵C表达为
其中,Cn=[C1n C2n ... CNn]T为校准矩阵C中的第n列,称为第n个阵元的端口耦合系数,表示第n个阵元端口对各阵元端口的耦合,Cdn称为第n个阵元端口对第d个阵元端口的耦合元素,n=1,2,...,N,d=1,2,...,N,()T表示转置;
定义以单个阵元为圆心,xλ为半径的圆为该单个阵元的耦合区,其中x为不小于1的实数,λ为阵列天线工作频率的波长,对任意两个阵元而言,如果其耦合区内的阵元数量和分布相同,则所述两个阵元互为等耦合阵元,令所述阵列天线包括U组互不相同的等耦合阵元,在所述U组互不相同的等耦合阵元的每组中分别任意选取1个阵元组成包括U个阵元的子阵列,所述步骤S2中,向所述子阵列馈入T组线性无关的端口激励I1,I2...IT,T≥U,获得相应的所述子阵列的T组口径场激励I′1,I′2...I′T,构建矩阵III,III=[I1 I2 ... IT],构建矩阵III′,III′=[I′1 ... I′T],所述步骤S3包括如下两个步骤:
S301.根据III′=C′×III计算得到所述子阵列的校准矩阵C′,
S302.根据校准矩阵C′中的端口耦合系数以及各阵元的相对位置,计算得到校准矩阵C。
4.根据权利要求3所述的阵列天线的校准方法,其特征在于,所述步骤S302包括如下步骤:
S3021.在所述包括U个阵元的子阵列中任意指定一个阵元A,阵元A的耦合区为A′,提取阵元A的端口耦合系数CA中阵元A端口对耦合区A′中各阵元端口的耦合元素;
S3022.在所述阵元A的等耦合阵元中任意指定一个阵元B,阵元B的耦合区为B′,阵元B端口对耦合区B′中任意一个阵元端口的耦合元素与阵元A端口对耦合区A′中同样相对位置的阵元端口的耦合元素相等,据此可确定阵元B的端口耦合系数CB中阵元B端口对耦合区B′中的阵元端口的耦合元素,端口耦合系数CB中阵元B端口对耦合区B′之外的阵元端口的耦合元素为0,其中,所述“同样相对位置”指耦合区中的阵元相对于该耦合区圆心的位置。
5.根据权利要求1-4中任一所述的阵列天线的校准方法,其特征在于,所述阵中方向图通过测量获得,或者基于阵列天线的物理参数或/和机械模型或/和仿真模型通过仿真获得,所述阵列天线的物理参数包括天线形式和阵列结构。
6.根据权利要求5所述的阵列天线的校准方法,其特征在于,所述任意两个等耦合阵元的阵中方向图相同。
7.根据权利要求1,2,3,4,6中任一所述的阵列天线的校准方法,其特征在于,所述步骤S203中阵中方向图、阵中方向图中心的位置、测量点的位置、测量数据Et以及口径场激励I′t满足关系式:Et=YI′t,其中Et是所述Mt个测量点测得的电/磁场,为Mt×1的矩阵,Y是阵元到测量点的幅相变换矩阵,Y根据阵中方向图、阵中方向图中心的位置以及测量点的位置获得。
8.根据权利要求7所述的阵列天线的校准方法,其特征在于,以任意参考点为原点建立球坐标系,第n个阵元的阵中方向图中心的位置的坐标为(Rn,θn,φn),n=1,2,...,N,第n个阵元的阵中方向图表示为fn(θ,φ),第m个测量点的位置的坐标为(R′m,θ′m,φ′m),m=1,2,...,M,所述阵元到测量点的幅相变换矩阵Y为
10.根据权利要求1,2,3,4,6,8,9中任一所述的阵列天线的校准方法,其特征在于,所述测量点位于阵元的辐射远场。
11.根据权利要求1,2,3,4,6,8,9中任一所述的阵列天线的校准方法,其特征在于,当Mt>N/3时,通过最小二乘法计算得到所述口径场激励I′t;当T>N时,通过最小二乘法计算得到所述校准矩阵C;当T>U时,通过最小二乘法计算得到所述校准矩阵C′。
12.一种阵列天线的校准设备,其特征在于,所述校准设备包括:
阵中方向图获取模块,用于获取所述阵列天线的阵元的阵中方向图以及阵中方向图中心的位置;
馈电模块,用于向所述阵列天线馈入端口激励It;
信号收发模块,与测量天线连接,用于获得Mt个测量点的位置、通过测量天线在Mt个测量点处向所述阵列天线发射测量信号并获得所述阵列天线电/磁场的测量数据,所述测量数据包含幅度和相位信息,Mt≥N/3,N为所述阵列天线的阵元数;
口径场激励获取模块,用于根据所述阵中方向图、阵中方向图中心的位置、测量点的位置以及测量数据获得口径场激励I′t;
校准矩阵获取模块,用于根据端口激励It与口径场激励I′t的关系I′t=C×It计算得到校准矩阵C;
校准模块,用于根据所述校准矩阵C对阵列天线各阵元进行校准。
13.根据权利要求12所述的阵列天线的校准设备,其特征在于,所述阵中方向图获取模块包括:
测量单元,用于测量获得所述阵列天线的阵中方向图;或/和
仿真单元,用于基于阵列天线的物理参数或/和机械模型或/和仿真模型通过仿真获得所述阵列天线的阵中方向图,所述物理参数包括天线形式和阵列结构。
14.一种阵列天线的校准设备,其特征在于,所述校准设备包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如权利要求1,2,3,4,6,8,9中任一所述方法的步骤。
15.一种阵列天线的校准系统,包括电波暗室以及测量天线,其特征在于,在所述校准系统中集成如权利要求12-14中任一所述的校准设备。
16.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1,2,3,4,6,8,9中任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711168211.6A CN109818689B (zh) | 2017-11-21 | 2017-11-21 | 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711168211.6A CN109818689B (zh) | 2017-11-21 | 2017-11-21 | 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109818689A CN109818689A (zh) | 2019-05-28 |
CN109818689B true CN109818689B (zh) | 2021-06-25 |
Family
ID=66600629
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711168211.6A Active CN109818689B (zh) | 2017-11-21 | 2017-11-21 | 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109818689B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112583498B (zh) * | 2020-11-25 | 2023-06-30 | 上海毫微太科技有限公司 | 阵列天线系统校准方法、装置、设备及存储介质 |
CN112904094B (zh) * | 2021-02-04 | 2022-08-23 | 中国人民解放军国防科技大学 | 一种基于无人机的口面天线外场测试方法 |
CN113824512B (zh) * | 2021-09-13 | 2023-10-10 | 中信科移动通信技术股份有限公司 | 大规模天线调测方法、测试设备及计算机设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102386983A (zh) * | 2011-10-17 | 2012-03-21 | 中国舰船研究设计中心 | 船舶大型阵列天线间电磁耦合反衍预测方法 |
CN103856272A (zh) * | 2012-12-03 | 2014-06-11 | 深圳市通用测试系统有限公司 | Mimo无线终端的无线性能测试方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3224644A4 (en) * | 2014-11-12 | 2018-11-07 | Emscan Corporation | Reactive near-field antenna measurement |
-
2017
- 2017-11-21 CN CN201711168211.6A patent/CN109818689B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102386983A (zh) * | 2011-10-17 | 2012-03-21 | 中国舰船研究设计中心 | 船舶大型阵列天线间电磁耦合反衍预测方法 |
CN103856272A (zh) * | 2012-12-03 | 2014-06-11 | 深圳市通用测试系统有限公司 | Mimo无线终端的无线性能测试方法 |
Non-Patent Citations (1)
Title |
---|
全极化相控阵雷达精密测量技术;李棉全;《CNKI中国优秀博士学位论文全文数据库 信息科技辑》;20151115;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109818689A (zh) | 2019-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108417999B (zh) | 多模相控阵天线和用于加宽其波束的方法 | |
CN109813967B (zh) | 一种阵列天线方向图的测量方法、设备、系统 | |
CN109818689B (zh) | 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 | |
CN109813969B (zh) | 一种阵列天线的诊断方法、设备、系统 | |
CN114531182A (zh) | 阵列天线的校准方法、装置及存储介质 | |
US11811461B2 (en) | Calibration method for a phased array of antennas | |
CN105606906B (zh) | 一种毫米波相控阵测试标定方法 | |
CN112083234A (zh) | 阵列天线总辐射功率测量方法、装置以及计算机存储介质 | |
WO2016089738A1 (en) | Satellite-based phased array calibration | |
CN104375129B (zh) | 一种分布式阵列相参合成雷达发射相参参数校准方法 | |
JP6678554B2 (ja) | アンテナ測定装置 | |
CN109818688B (zh) | 一种阵列天线的校准方法、设备、系统以及计算机可读存储介质 | |
CN115360518A (zh) | 一种相控阵天线的通道幅相标校方法 | |
Grundmann et al. | Investigation of direction of arrival estimation using characteristic modes | |
CN104852775A (zh) | 一种大型相控阵天线波束性能的快速评估方法 | |
CN109813968B (zh) | 一种阵列天线方向图的测量方法、设备、系统 | |
CN109815509B (zh) | 一种阵列天线的诊断方法、设备、系统以及计算机可读存储介质 | |
CN114915356B (zh) | 相位阵列天线校正方法 | |
CN117031418A (zh) | Sar卫星方位敏捷观测模式扫描指向检测方法及系统 | |
RU2392634C1 (ru) | Способ определения направлений на источники излучения и углового разрешения источников | |
CN111965602B (zh) | 一种相控阵雷达幅相一致性检测方法和系统 | |
CN118074826B (zh) | 一种基于cnn-lstm采用功率测量的相控阵校准方法 | |
Zhang et al. | Structural mutual coupling calibration experiments of full polarization conformal phased array | |
Wei et al. | Research on fast measurement method for active element pattern of phased array antenna | |
Ronggui et al. | Electrical compensation for deformed antenna arrays and its experimental verification |
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 |