CN110492917B - 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法 - Google Patents

一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法 Download PDF

Info

Publication number
CN110492917B
CN110492917B CN201910709956.1A CN201910709956A CN110492917B CN 110492917 B CN110492917 B CN 110492917B CN 201910709956 A CN201910709956 A CN 201910709956A CN 110492917 B CN110492917 B CN 110492917B
Authority
CN
China
Prior art keywords
channel
calibration
channels
transmitting
amplitude
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
Application number
CN201910709956.1A
Other languages
English (en)
Other versions
CN110492917A (zh
Inventor
洪伟
郭翀
田玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201910709956.1A priority Critical patent/CN110492917B/zh
Publication of CN110492917A publication Critical patent/CN110492917A/zh
Application granted granted Critical
Publication of CN110492917B publication Critical patent/CN110492917B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/10Monitoring; Testing of transmitters
    • H04B17/11Monitoring; Testing of transmitters for calibration
    • H04B17/12Monitoring; Testing of transmitters for calibration of transmit antennas, e.g. of the amplitude or phase
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/20Monitoring; Testing of receivers
    • H04B17/21Monitoring; Testing of receivers for calibration; for correcting measurements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/084Equal gain combining, only phase adjustments
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0868Hybrid systems, i.e. switching and combining

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明公开了一种全数字波束成形阵列以及实现收发通道幅度、相位自动校准方法,在天线阵列口径面上嵌入一条传输线,传输线两端接有校准收发通道;发射校准时,依次激励每个发射通道,同时记录下两校准接收通道采样得到的信号幅度和相位值,计算可得所有发射通道相对于参考通道的幅度、相位差;接收校准时,系统所有接收通道处于工作状态,相继激励两侧的校准发射通道,记录所有接收通道采样得到的校准信号幅度、相位值,可推得所有接收通道相对于参考通道的幅度、相位差。本发明提供的全数字波束成形阵列收发通道自动校准的方法,没有场地和人力的限制,在硬件上减少了大量定向耦合器的引入,适合用于大规模全数字波束成形阵列的自校准。

Description

一种全数字波束成形阵列及实现收发通道幅度、相位自动校 准方法
技术领域
本发明属于微波毫米波通信系统的校准技术领域,尤其涉及一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法。
背景技术
移动通信技术经过2G,3G,4G时代的演变,如今将逐步进入5G时代。相较于4G技术,5G希望在网络容量上实现1000倍提升,同时实现小于1ms的数据传输延时。为了实现上述目标,将采用超密集组网(Ultra-densification)、毫米波(mm-Wave)、大规模天线阵(MassiveMIMO)等关键技术。其中,全数波束成形系统因为能在数字域内方便、快速地实现波束成形、多波束和下行数据的预编码,减少了模拟移相器的使用,所以是实现大规模MIMO最直接有效的系统方案。
数字波束成形系统是幅度、相位相参系统。为了实现精确的波束成形,除了要对整个系统进行频率相参和时钟同步考虑外,还需要对整个系统进行幅度、相位校准,简称通道的幅/相校准。由于通道数众多,如何实现系统快速、准确的幅/相校准值得研究。现有的校准技术大部分针对有源相控系统。由于全数字多波束系统有其固有的架构特点,在校准方式上会有所不同。数字波束成形系统中存在大量的有源器件。受温度、时间等外在因素的变化,会造成器件电特性的改变,从而改变系统收发通道的幅度和相位,影响波束赋形的精度。另一方面,对于架设于高处的基站,对系统再次进行暗室校准或地面校准不切实际。所以如何实现系统通道的幅/相自动校准,保证精确的波束成形,值得研究。
目前,常用的数字波束成形系统的自校准方式多针对接收通道的校准,对于发射通道或者收发通道同时校准的方式鲜有提及。同时在硬件上,通常在每个射频通道末级放置定向耦合器,耦合部分能量用来实现系统校准,增加了系统硬件的复杂性。
在天线阵列中嵌入传输线进行阵列通道的检测的方法早期被应用于大型相控阵校准中。由于全数字波束成形系统不需要通道合路,所以其在校准方法上与相控阵系统存在些许差异。利用阵列中嵌入传输线耦合通道能量的方式相较于通道末级加耦合器的方式在硬件上更为简单。同时,利用阵列路径的对称性和单元与传输线耦合的相似性,可以更好地完成全数字波束成形阵列通道的自动校准。
发明内容
发明目的:针对以上问题,本发明提出一种能实现全数字波束成形系统收发通道幅度、相位自动校准的方法,该方法针对收发通道均可以实现自动校准,并尽可能简化了校准硬件电路。
技术方案:本发明所述的一种全数字波束成形系统收发通道幅度、相位自动校准方案在硬件上需要在天线阵列口径面上加一条微带传输线。天线单元采用矩形贴片形式,每个单元结构相同且等间距排布成阵,整个阵列排布在一张印刷电路板上。耦合传输线采用微带传输线形式,同样印刷在该印刷电路板上,且与天线阵列有一定距离间隔,其到每个天线单元位置一定。天线结构对称,整个阵列结构是单元结构的重复排列。在耦合传输线两端接有相应的校准通道,校准通道具有收发模式,用开关进行切换。
本发明的天线阵列结构如下:一种全数字波束成形阵列,该阵列天线包括天线单元(1)、耦合传输线(2)、具有收发模式的校准通道(3)和(4)、阵列收发通道(5),阵列规模为1×N,N为通道总数;所述天线单元(1)与收发通道(5)连接,且一一对应,并且天线单元(1)和收发通道(5)周期性重复排列;耦合传输线(2)位于阵列天线口径面上,耦合传输线(2)到每个天线单元(1)位置一定,且整个天线结构对称,在耦合传输线(2)两端均接有具有收发模式的校准通道(3)和(4)。
并且,在上述天线阵列结构的基础上实现的收发通道幅度、相位自动校准方法,该方法包括如下步骤:
(1)记校准通道(3)和(4)为P1、P2,天线系统处于发射状态发射校准信号,两校准通道P1、P2处于接收状态,依次激励各发射通道,并记录耦合传输线两端的校准通道接收到的信号幅度大小和相位;
(2)利用路径对称和各单元到耦合传输线的耦合量约相等的特性,计算相邻两个发射通道的通道幅度和相位差,用矢量βt n,n+1表示;
(3)计算出所有发射通道相对于参考通道的幅度和相位差,用矢量αt L,n表示,完成发射通道自校准;
(4)天线系统处于接收状态,首先只激励一侧校准通道P1,另一侧校准通道P2处于接收状态用于校准信号匹配,记录阵列中每个接收通道采样得到的耦合信号幅度大小和相位值;
(5)只激励校准通道P2,校准通道P1处于接收状态,同样记录每个接收通道采样得到的耦合信号幅度大小和相位,再次利用路径对称和各单元到耦合传输线的耦合量约相等的特性,计算出相邻接收通道幅度和相位差,用矢量βr n,n+1表示;
(6)计算出所有接收通道相对于参考通道的幅度和相位差,用矢量αr L,n表示,完成接收通道的自校准。
进一步,步骤(2)中,计算出相邻两个发射通道的通道幅度和相位差,方法如下:
(2.1)将全数字波束成形阵列收发通道从左至右依次编号1,2,…,N,N表示通道总数,随机选择一个编号为L的通道作为参考通道,根据此参考通道进行校准得到的所有发射、接收通道的幅度、相位值;
(2.2)假设Tn表示阵列中第n个发射通道的通道状态信息,n∈N,其相位和幅度是随机量,需校准求解,Cn表示第n个天线单元到传输线的耦合量,Ljn表示第n个发射通道的发射信号耦合到传输线上之后,经过传输线的传输,直到被第j个校准通道采样的路径,j=1,2,也即P1、P2
(2.3)在t1时刻,只有第n个发射通道被激励,其余通道处于未被激励状态或通道衰减值最大状态,天线与耦合传输线间的耦合信号将被校准通道P1和P2采集,该采样信号分别用an,bn表示为:
Figure BDA0002153384660000031
Figure BDA0002153384660000032
其中,an表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn+1,n表示第n个通道相对于第n+1个通道的发射校准信号由于传输路径不同造成的校准通道P2采样信号值之间的差异,下标有方向性,其与天线S散射参数关系可表示为:
Figure BDA0002153384660000033
其中,
Figure BDA0002153384660000034
表示第n个发射通道发射的信号传输到P2校准通道的传输量,
Figure BDA0002153384660000035
表示第n+1个发射通道的发射信号传输到P2校准通道的传输量;
(2.4)在t2时刻,只有第n+1个通道被激励,其余通道处于未被激励状态或通道衰减值最大状态,天线与传输线间的耦合信号将被校准通道P1和P2采集,an+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn,n+1表示第n+1个通道相对于第n个通道的发射信号由于传输路径不同造成的校准通道P1采样信号值之间的差异,下标有方向性,an+1,bn+1的值可以表示为:
Figure BDA0002153384660000041
Figure BDA0002153384660000042
在阵列中天线单元是相同的,且整个阵列结构对称,有以下假设成立:
Cn≈Cn+1,δn+1,n≈δn,n+1 (8)
结合假设(8),用(6)/(3),(7)/(4),得到:
Figure BDA0002153384660000043
Figure BDA0002153384660000044
进一步,得到:
Figure BDA0002153384660000045
将κn,n+1的幅度用dB表示为(12),角度用rad表示为(13):
Figure BDA0002153384660000046
Figure BDA0002153384660000047
其中,mag_dB表示矢量信号取幅度,且表示为dB形式,arg表示矢量信号取角度;相邻通道的校准因子βt n,n+1表示为:
Figure BDA0002153384660000048
Figure BDA00021533846600000411
Figure BDA0002153384660000049
βt n,n+1是矢量也是相对量,其表达形式为
Figure BDA00021533846600000410
其幅度代表第n+1个发射通道相对于第n个发射通道的幅度差,相位代表第n+1个发射通道相对于第n个发射通道的相位差;
假设βt n,n+1=Tn+1/Tn=κn,n+1,将其代入公式(9)中,得到真正要求的δn,n+1的推测值δ’n,n+1,该推测值δ’n,n+1与真实值δn,n+1的关系表示为:
Figure BDA0002153384660000051
即δn,n+1的推测值δ’n,n+1与其真实值比较,幅度相同,相位上存在π的模糊;
(2.5)通过仿真得到天线阵列的S参数,将仿真结果作为判断依据,排除掉推测值δ’n,n+1中相位存在π模糊的值,从而得到真正的δn,n+1与βt n,n+1
(2.6)依次类推,得到所有发射通道的相邻通道校准因子βt i,i+1,其中,i=1,2,..,N-1,N代表通道总数。
进一步的,步骤(3)中完成发射通道自校准,方法如下:阵列中第n个发射通道相对于参考通道的通道校准系数αt L,n可以表示为:
Figure BDA0002153384660000052
至此,可以完成所有发射通道相对于参考通道L的幅度、相位校准,也即完成了整个阵列发射通道的幅度、相位自校准过程。
进一步的,步骤(4)的具体方法如下:
设置阵列处于接收状态,所有通道均处于工作状态,在t1时刻,校准通道P1处于发射状态,通过耦合线向阵列注入校准信号,校准通道P2处于接收状态,只作信号匹配用途,所有接收通道同时采样来自传输线耦合得到的信号,该过程中第n个和第n+1个接收通道采样得到的信号pn,pn+1可以表示为:
Figure BDA0002153384660000053
Figure BDA0002153384660000054
其中,Rn表示阵列中第n个接收通道的通道状态信息,其相位和幅度为未知的,需要校准得到,Cn表示传输线到第n个天线单元的耦合量,Ln,1表示第1个校准通道发射的校准信号经过传输线的传输,直到第n个接收通道采样所经过的传输路径,δn+1,n表示第n个通道与第n+1个通道采样得到的校准信号,由于不同路径传输引起的信号幅度、相位差,下标有方向性,n=1,2,…,N-1。
进一步的,步骤(5)的具体方法如下:
设置在t2时刻,校准通道P2处于发射状态,通过耦合线向阵列注入校准信号,校准通道P1处于接收状态,只作信号匹配用途,所有接收通道同时采样来自传输线耦合得到的信号,这样,第n个和第n+1个接收通道采样得到的信号qn,qn+1可以表示为:
Figure BDA0002153384660000061
Figure BDA0002153384660000062
根据公式(17)-(20),并结合假设(8),可以得到:
Figure BDA0002153384660000063
Figure BDA0002153384660000064
同发射通道校准过程(11)-(15),可以得到相邻接收通道校准因子
Figure BDA0002153384660000066
同发射校准因子βt n,n+1一样,其也是矢量且为相对值,其幅度代表第n+1个接收通道相对于第n个接收通道的幅度差,相位代表第n+1个接收通道相对于第n个接收通道的相位差。
进一步的,步骤(6)的具体方法如下:阵列中第n个接收通道相对于参考通道(L通道)的校准系数
Figure BDA0002153384660000067
可以表示为:
Figure BDA0002153384660000065
至此,可以完成所有接收通道幅度、相位自校准过程。
有益效果:与现有技术相比,本发明的技术方案具有以下有益技术效果:
本发明公开了一种全数字波束成形阵列收发通道幅度、相位自动校准方法,可以同时实现对阵列收、发通道的通道幅度和相位的自动校准。同时,与其它自校准方案相比,简化了校准硬件电路,校准时间较短,这就使得在全数字大规模天线阵列中采用该自动校准方案硬件成本较低,且几乎不影响系统的正常工作状态。
附图说明
图1为本发明校准硬件实施方案图;
图2为本发明具体实施方式中发射通道校准过程图;
图3为本发明具体实施方式中接收通道校准过程图;
图4为本发明工作流程图;
图5为本发明实测得到的一个1×8全数字波束阵列经过发射通道自校准后得到的波束合成方向图;
图6为本发明实测得到的一个1×8全数字波束阵列经过接收通道自校准后得到的波束合成方向图。
具体实施方式
下面结合具体实施方式和附图对本发明的技术方案作进一步的介绍。
本具体实施方式公开了一种全数字波束成形阵列收发通道幅度、相位自动校准方法,如图1所示。1表示天线单元,单元采用矩形贴片形式,每个单元结构相同且等间距排布成阵,整个阵列排布在一张印刷电路板上。天线结构对称,整个阵列结构是单元结构的重复排列。2表示耦合传输线,耦合传输线采用微带传输线形式,同样印刷在该印刷电路板上,且与天线阵列有一定距离间隔。通过空间电磁耦合,每个天线单元辐射的能量均可以耦合到传输线上。同理,传输线上的能量也可以耦合到每个天线单元上。该间距决定了天线单元到传输线的耦合量大小,并且耦合传输线到每个天线单元间隔相同。3、4表示两个校准通道,它们相应地接在耦合传输线两端,采样由系统发射通道耦合到传输线上的校准信号(发射通道校准)或给系统接收通道提供校准信号源(接收通道校准)。校准通道具有收发模式,用开关进行切换。5表示系统收发通道,具有接收和发射两种模式,用开关作模式切换。它们和天线单元一一对应,通过天线单元将通道能量辐射出去。系统通道从左至右编号1,2,…,N,参考通道为第L通道。全数字波束成形系统需要知道各通道相对于参考通道的幅度、相位差量,该通道幅度、相位差量是未知量,可利用该自动校准方法计算得到。
下面结合附图对阵列收发通道自校准流程作更进一步说明:
(1)如图2所示,先对系统发射通道进行自校准。假设Tn表示阵列中第n个发射通道的通道状态信息(n∈N),其相位和幅度为未知的,Cn表示第n个天线单元到传输线的耦合量,P1、P2分别表示图1中的校准通道3和校准通道4,Ljn表示第n个发射通道的发射信号耦合到传输线上之后,经过传输线的传输,直到被第j个校准通道采样的路径,j=1,2。
如图2(a)(b),在t1时刻,只有第n个发射通道被激励,其余通道处于未被激励状态或通道衰减值最大状态,黑色实框表示被激励状态,斜线框表示未被激励状态,以下相同。这样,天线与传输线间的耦合信号将被校准通道P1和P2采集,该采样信号分别用an,bn表示为:
Figure BDA0002153384660000081
Figure BDA0002153384660000082
其中,an表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn+1,n表示第n个通道相对于第n+1个通道的发射校准信号由于传输路径不同造成的校准通道P2采样信号值之间的差异,下标有方向性,其与天线S散射参数关系可表示为:
Figure BDA0002153384660000083
其中,
Figure BDA0002153384660000084
表示第n个发射通道发射的信号传输到P2校准通道的传输量(当其他端口匹配)。同理,
Figure BDA0002153384660000085
表示第n+1个发射通道的发射信号传输到P2校准通道的传输量(当其他端口匹配)。
(2)如图2(c)(d)所示,在t2时刻,只有第n+1个通道被激励,其余通道处于未被激励状态或通道衰减值最大状态。这样,天线与传输线间的耦合信号将被校准通道P1和P2采集,an+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn,n+1表示第n+1个通道相对于第n个通道的发射信号由于传输路径不同造成的校准通道P1采样信号值之间的差异,下标有方向性,an+1,bn+1的值可以表示为:
Figure BDA0002153384660000086
Figure BDA0002153384660000087
因为在阵列中天线单元是相同的,且整个阵列结构对称,有以下假设成立:
Cn≈Cn+1,δn+1,n≈δn,n+1 (8)
结合假设(8),用(6)/(3),(7)/(4),得到:
Figure BDA0002153384660000088
Figure BDA0002153384660000091
进一步,得到:
Figure BDA0002153384660000092
将κn,n+1的幅度用dB表示为(12),角度用rad表示为(13):
Figure BDA0002153384660000093
Figure BDA0002153384660000094
其中,mag_dB表示矢量信号取幅度,且表示为dB形式,arg表示矢量信号取角度。
相邻通道的校准因子βt n,n+1就可以表示为:
Figure BDA0002153384660000095
Figure BDA0002153384660000099
Figure BDA0002153384660000096
βt n,n+1是矢量也是相对量,有
Figure BDA0002153384660000098
的表达形式。其幅度代表第n+1个发射通道相对于第n个发射通道的幅度差,相位代表第n+1个发射通道相对于第n个发射通道的相位差。
可以看到,计算得到的校准因子的相位存在π模糊。为了解决这个问题,加一步判断,即假设βt n,n+1=Tn+1/Tn=κn,n+1,将其代入公式(9)中,得到真正要求的δn,n+1的推测值δ’n,n+1。该推测值δ’n,n+1与真实值δn,n+1的关系可以表示为:
Figure BDA0002153384660000097
即δn,n+1的推测值δ’n,n+1与其真实值比较,幅度相同,相位上存在π的模糊。
由(5)可知δn,n+1与天线S散射参数存在关系,而天线阵列的S参数可以事先通过全波仿真软件得到。因此,将仿真结果作为判断依据,可以排除掉推测值δ’n,n+1中相位存在π模糊的值,从而得到真正的δn,n+1与βt n,n+1
依次类推,得到所有发射通道的相邻通道校准因子βt i,i+1,其中,i=1,2,..,N-1,N代表通道总数。
(4)阵列中第n个发射通道(n=1,2,..,N)相对于参考通道(第L通道)的通道校准系数αt L,n可以表示为:
Figure BDA0002153384660000101
至此,可以完成所有发射通道相对于L通道的幅度、相位校准,也即完成了整个阵列发射通道的幅度、相位自校准过程。
(5)完成发射通道校准后,下一步进行阵列的接收通道校准。阵列处于接收状态,所有通道均处于工作状态,如图3所示。如图3(a)所示,在t1时刻,校准通道1处于发射状态,通过耦合线向阵列注入校准信号。校准通道2处于接收状态,只作信号匹配用途。所有接收通道同时采样来自传输线耦合得到的信号。这样,该过程中第n个和第n+1个接收通道采样得到的信号pn,pn+1可以表示为:
Figure BDA0002153384660000102
Figure BDA0002153384660000103
其中,Rn表示阵列中第n个接收通道的通道状态信息(n∈N),其相位和幅度为未知的,Cn表示第传输线到第n个天线单元的耦合量,Ln1表示第1个校准通道发射的校准信号经过传输线的传输,直到第n个接收通道(n=1,2,..,N)采样所经过的传输路径。δn+1,n表示第n个通道与第n+1个通道采样得到的校准信号,由于不同路径传输引起的信号幅度、相位差,下标有方向性,
(6)在t2时刻,校准通道2处于发射状态,通过耦合线向阵列注入校准信号。校准通道1处于接收状态,只作信号匹配用途。所有接收通道同时采样来自传输线耦合得到的信号。这样,第n个和第n+1个接收通道采样得到的信号qn,qn+1可以表示为:
Figure BDA0002153384660000104
Figure BDA0002153384660000105
根据公式(17)-(20)并结合假设(8),可以得到:
Figure BDA0002153384660000111
Figure BDA0002153384660000112
同发射通道校准过程(11)-(15),可以得到相邻接收通道校准因子
Figure BDA0002153384660000113
(7)阵列中第n个接收通道(n=1,2,..,N)相对于参考通道(L通道)的校准系数
Figure BDA0002153384660000114
可以表示为:
Figure BDA0002153384660000115
至此,可以完成所有接收通道幅度、相位自校准过程。
上述自校准流程图如图4所示。
将上述自校准过程在一个1×8全数字波束合成阵列中做验证,其经过发射通道自校准过程后的波束合成方向图如图5所示,经过接收通道自校准过程后的波束合成方向图如图6所示。通过未校准与校准后的方向图对比结果,证明了上述自校准过程的可行性。
以上所述仅是本发明的实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (4)

1.根据一种全数字波束成形阵列实现的收发通道幅度、相位自动校准方法,其特征在于,该方法包括如下步骤:
(1)记校准通道(3)和(4)为P1、P2,天线系统处于发射状态发射校准信号,两校准通道P1、P2处于接收状态,依次激励各发射通道,并记录耦合传输线两端的校准通道接收到的信号幅度大小和相位;
(2)利用路径对称和各单元到耦合传输线的耦合量约相等的特性,计算相邻两个发射通道的通道幅度和相位差,用矢量βt n,n+1表示;
(3)计算出所有发射通道相对于参考通道的幅度和相位差,用矢量αt L,n表示,完成发射通道自校准;
(4)天线系统处于接收状态,首先只激励一侧校准通道P1,另一侧校准通道P2处于接收状态用于校准信号匹配,记录阵列中每个接收通道采样得到的耦合信号幅度大小和相位值;
(5)只激励校准通道P2,校准通道P1处于接收状态,同样记录每个接收通道采样得到的耦合信号幅度大小和相位,再次利用路径对称和各单元到耦合传输线的耦合量约相等的特性,计算出相邻接收通道幅度和相位差,用矢量βr n,n+1表示;
(6)计算出所有接收通道相对于参考通道的幅度和相位差,用矢量αr L,n表示,完成接收通道的自校准;
所述一种全数字波束成形阵列包括天线单元(1)、耦合传输线(2)、具有收发模式的校准通道(3)和(4)、阵列收发通道(5),阵列规模为1×N,N为通道总数;所述天线单元(1)与收发通道(5)一一对应连接,并且天线单元(1)和收发通道(5)作为整体等间距周期性重复排列成阵;耦合传输线(2)位于阵列天线口径面上,耦合传输线(2)到每个天线单元(1)位置一定,在耦合传输线(2)两端均接有具有收发模式的校准通道(3)和(4),且整个天线结构对称;
步骤(2)中,计算出相邻两个发射通道的通道幅度和相位差,方法如下:
(2.1)将全数字波束成形阵列收发通道从左至右依次编号1,2,…,N,N表示通道总数,随机选择一个编号为L的通道作为参考通道,根据此参考通道进行校准得到的所有发射、接收通道的幅度、相位值;
(2.2)假设Tn表示阵列中第n个发射通道的通道状态信息,n∈N,其相位和幅度是随机量,需校准求解,Cn表示第n个天线单元到传输线的耦合量,Ljn表示第n个发射通道的发射信号耦合到传输线上之后,经过传输线的传输,直到被第j个校准通道采样的路径,j=1,2,也即P1、P2
(2.3)在t1时刻,只有第n个发射通道被激励,其余通道处于未被激励状态或通道衰减值最大状态,天线与耦合传输线间的耦合信号将被校准通道P1和P2采集,采样信号分别用an,bn表示为:
Figure FDA0003657641600000021
Figure FDA0003657641600000022
其中,an表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn表示第n个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn+1,n表示第n个通道相对于第n+1个通道的发射校准信号由于传输路径不同造成的校准通道P2采样信号值之间的差异,下标有方向性,其与天线S散射参数关系可表示为:
Figure FDA0003657641600000023
其中,
Figure FDA0003657641600000024
表示第n个发射通道发射的信号传输到P2校准通道的传输量,
Figure FDA0003657641600000025
表示第n+1个发射通道的发射信号传输到P2校准通道的传输量;
(2.4)在t2时刻,只有第n+1个通道被激励,其余通道处于未被激励状态或通道衰减值最大状态,天线与传输线间的耦合信号将被校准通道P1和P2采集,an+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P1采集得到的信号;bn+1表示第n+1个发射通道的信号耦合到传输线后,经过传输被校准通道P2采集得到的信号;δn,n+1表示第n+1个通道相对于第n个通道的发射信号由于传输路径不同造成的校准通道P1采样信号值之间的差异,下标有方向性,an+1,bn+1的值可以表示为:
Figure FDA0003657641600000026
Figure FDA0003657641600000027
在阵列中天线单元是相同的,且整个阵列结构对称,有以下假设成立:
Cn≈Cn+1,δn+1,n≈δn,n+1(8)
结合假设(8),用(6)/(3),(7)/(4),得到:
Figure FDA0003657641600000028
Figure FDA0003657641600000031
进一步,得到:
Figure FDA0003657641600000032
将κn,n+1的幅度用dB表示为(12),角度用rad表示为(13):
Figure FDA0003657641600000033
Figure FDA0003657641600000034
其中,mag_dB表示矢量信号取幅度,且表示为dB形式,arg表示矢量信号取角度;相邻通道的校准因子βt n,n+1表示为:
Figure FDA0003657641600000035
Figure FDA0003657641600000036
Figure FDA0003657641600000037
βt n,n+1是矢量也是相对量,其表达形式为
Figure FDA0003657641600000038
其幅度代表第n+1个发射通道相对于第n个发射通道的幅度差,相位代表第n+1个发射通道相对于第n个发射通道的相位差;
假设βt n,n+1=Tn+1/Tn=κn,n+1,将其代入公式(9)中,得到真正要求的δn,n+1的推测值δ’n,n+1,该推测值δ’n,n+1与真实值δn,n+1的关系表示为:
Figure FDA0003657641600000039
即δn,n+1的推测值δ’n,n+1与其真实值比较,幅度相同,相位上存在π的模糊;
(2.5)通过仿真得到天线阵列的S参数,将仿真结果作为判断依据,排除掉推测值δ’n,n+1中相位存在π模糊的值,从而得到真正的δn,n+1与βt n,n+1
(2.6)依次类推,得到所有发射通道的相邻通道校准因子βt i,i+1,其中,i=1,2,..,N-1,N代表通道总数;
步骤(3)中完成发射通道自校准,方法如下:阵列中第n个发射通道相对于参考通道的通道校准系数αt L,n可以表示为:
Figure FDA0003657641600000041
至此,可以完成所有发射通道相对于参考通道L的幅度、相位校准,也即完成了整个阵列发射通道的幅度、相位自校准过程。
2.根据权利要求1所述的一种全数字波束成形阵列收发通道幅度、相位自动校准方法,其特征在于,步骤(4),具体方法如下:
设置阵列处于接收状态,所有通道均处于工作状态,在t1时刻,校准通道P1处于发射状态,通过耦合线向阵列注入校准信号,校准通道P2处于接收状态,只作信号匹配用途,所有接收通道同时采样来自传输线耦合得到的信号,该过程中第n个和第n+1个接收通道采样得到的信号pn,pn+1可以表示为:
Figure FDA0003657641600000042
Figure FDA0003657641600000043
其中,Rn表示阵列中第n个接收通道的通道状态信息,其相位和幅度为未知的,需要校准得到,Cn表示传输线到第n个天线单元的耦合量,Ln,1表示第1个校准通道发射的校准信号经过传输线的传输,直到第n个接收通道采样所经过的传输路径,δn+1,n表示第n个通道与第n+1个通道采样得到的校准信号,由于不同路径传输引起的信号幅度、相位差,下标有方向性,n=1,2,…,N-1。
3.根据权利要求1所述的一种全数字波束成形阵列收发通道幅度、相位自动校准方法,其特征在于,步骤(5),具体方法如下:
设置在t2时刻,校准通道P2处于发射状态,通过耦合线向阵列注入校准信号,校准通道P1处于接收状态,只作信号匹配用途,所有接收通道同时采样来自传输线耦合得到的信号,这样,第n个和第n+1个接收通道采样得到的信号qn,qn+1可以表示为:
Figure FDA0003657641600000044
Figure FDA0003657641600000045
根据公式(17)-(20),并结合假设(8),可以得到:
Figure FDA0003657641600000051
Figure FDA0003657641600000052
同发射通道校准过程(11)-(15),可以得到相邻接收通道校准因子
Figure FDA0003657641600000053
同发射校准因子βt n,n+1一样,其也是矢量且为相对值,其幅度代表第n+1个接收通道相对于第n个接收通道的幅度差,相位代表第n+1个接收通道相对于第n个接收通道的相位差。
4.根据权利要求3所述的一种全数字波束成形阵列收发通道幅度、相位自动校准方法,其特征在于,步骤(6)的具体方法如下:阵列中第n个接收通道相对于参考通道的校准系数
Figure FDA0003657641600000054
可以表示为:
Figure FDA0003657641600000055
至此,可以完成所有接收通道幅度、相位自校准过程。
CN201910709956.1A 2019-07-29 2019-07-29 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法 Active CN110492917B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910709956.1A CN110492917B (zh) 2019-07-29 2019-07-29 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910709956.1A CN110492917B (zh) 2019-07-29 2019-07-29 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法

Publications (2)

Publication Number Publication Date
CN110492917A CN110492917A (zh) 2019-11-22
CN110492917B true CN110492917B (zh) 2022-07-29

Family

ID=68549123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910709956.1A Active CN110492917B (zh) 2019-07-29 2019-07-29 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法

Country Status (1)

Country Link
CN (1) CN110492917B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111123220B (zh) * 2019-11-29 2023-11-07 瞬联软件科技(北京)有限公司 一种毫米波雷达的多通道幅相校准方法及系统
CN111585669B (zh) * 2020-04-29 2022-04-22 电子科技大学 一种用于天线阵列的接收发射通道同时校准方法
CN113296059B (zh) * 2021-05-25 2022-07-19 四川九洲空管科技有限责任公司 一种基于tcas定向天线发射波束控制合成全向方向图的方法
CN114512805B (zh) * 2022-01-19 2023-08-25 中国电子科技集团公司第十研究所 埋藏式夹心宽带耦合校准网络
CN115694683B (zh) * 2023-01-03 2023-03-21 成都实时技术股份有限公司 一种基于Lasso优化的数模异频多通道发射校准方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102426300B (zh) * 2011-08-31 2013-08-28 西安空间无线电技术研究所 一种星载波束形成接收通道幅相误差校准系统及其方法
CN107315183B (zh) * 2017-06-01 2020-06-26 西南电子技术研究所(中国电子科技集团公司第十研究所) 导航卫星阵列天线接收系统的校准方法
CN108155958A (zh) * 2017-11-22 2018-06-12 西南电子技术研究所(中国电子科技集团公司第十研究所) 大规模mimo天线阵列远场校准系统
CN108549058B (zh) * 2018-02-28 2022-05-13 四川九洲电器集团有限责任公司 一种二次雷达宽带有源相控阵系统及其动态校准方法

Also Published As

Publication number Publication date
CN110492917A (zh) 2019-11-22

Similar Documents

Publication Publication Date Title
CN110492917B (zh) 一种全数字波束成形阵列及实现收发通道幅度、相位自动校准方法
US10663563B2 (en) On-site calibration of array antenna systems
US6507315B2 (en) System and method for efficiently characterizing the elements in an array antenna
EP1468300B1 (en) Antenna measurement system
US8754811B1 (en) Digital beamforming phased array
EP2449626B1 (en) Self calibrating conformal phased array
CN105785328B (zh) 基于子阵划分的fda距离-角度解耦合波束形成方法
Guo et al. A hybrid adaptive antenna array for long-range mm-wave communications [antenna applications corner]
JP2019132836A (ja) Esa計測学の方法とシステム
CN109067439B (zh) 一种数字多波束阵发射装置采用的测试方法
US5926135A (en) Steerable nulling of wideband interference signals
JP2011014980A (ja) 電波伝搬パラメータ推定装置、電波伝搬パラメータ推定方法
CN109599677A (zh) 多极化数字波束形成天线
US20180375550A1 (en) Simultaneous communication through multiple beams of a wireless node
Medina-Sanchez Beam steering control system for low-cost phased array weather radars: Design and calibration techniques
CN115021833A (zh) 相控阵天线阵元通道一致性多模并行处理标校方法
GB2289798A (en) Improvements relating to radar antenna systems
CN111953429B (zh) 相控阵天线测试系统及测试方法
CN117155486A (zh) 一种桌面与外场相结合的相控阵阵面校准方法
GB2346013A (en) Calibration method for a phased array
CN114512805A (zh) 埋藏式夹心宽带耦合校准网络
Sharma et al. Mutual coupling compensation in circular dipole array for FM based passive coherent location radar
RU2393597C1 (ru) Антенна
Yin et al. Performance evaluation of a hybrid-beamforming sounder for 26 GHz channel measurements
Zhang et al. Parameter estimation for reconfigurable holographic surfaces enabled radars

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