CN106486769B - 用于线性相控阵天线的空间插值方法和设备 - Google Patents

用于线性相控阵天线的空间插值方法和设备 Download PDF

Info

Publication number
CN106486769B
CN106486769B CN201610796509.0A CN201610796509A CN106486769B CN 106486769 B CN106486769 B CN 106486769B CN 201610796509 A CN201610796509 A CN 201610796509A CN 106486769 B CN106486769 B CN 106486769B
Authority
CN
China
Prior art keywords
antenna element
bad
received signal
antenna elements
spatial interpolation
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
CN201610796509.0A
Other languages
English (en)
Other versions
CN106486769A (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.)
Halla Creedong Electronics Co ltd
SNU R&DB Foundation
Original Assignee
Seoul University School-Industry Cooperation Group
Co Wandu
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 Seoul University School-Industry Cooperation Group, Co Wandu filed Critical Seoul University School-Industry Cooperation Group
Publication of CN106486769A publication Critical patent/CN106486769A/zh
Application granted granted Critical
Publication of CN106486769B publication Critical patent/CN106486769B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • G01S7/4004Means for monitoring or calibrating of parts of a radar system
    • G01S7/4021Means for monitoring or calibrating of parts of a radar system of receivers
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q3/00Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
    • H01Q3/26Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06Systems determining position data of a target
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00Antenna arrays or systems
    • H01Q21/06Arrays of individually energised antenna units similarly polarised and spaced apart
    • H01Q21/08Arrays of individually energised antenna units similarly polarised and spaced apart the units being spaced along or adjacent to a rectilinear path
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q3/00Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
    • H01Q3/26Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture
    • H01Q3/30Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture varying the relative phase between the radiating elements of an array
    • H01Q3/34Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture varying the relative phase between the radiating elements of an array by electrical means

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

用于线性相控阵天线的空间插值方法和设备。根据本公开的实施方式,一种用于线性相控阵天线的空间插值方法涉及一种用于线性相控阵天线的空间插值方法,该线性相控阵天线包括多个发射天线单元和多个接收天线元件,并且该空间插值方法包括以下步骤:步骤1,在该步骤中,多个接收天线元件接收从目标反射的反射波;步骤2,在该步骤中,使用角估计算法来估计反射波入射到多个接收天线元件上的入射角;第3步,在该步骤中,从多个接收天线元件当中选择不良状况天线元件;以及步骤4,在该步骤中,对不良状况天线元件的接收信号进行补偿并且使用角估计算法来重新估计反射波入射到多个接收天线元件上的入射角。

Description

用于线性相控阵天线的空间插值方法和设备
技术领域
本公开涉及一种用于天线的空间插值方法和设备,并且更具体地,涉及用于线性相控阵天线的空间插值方法。
背景技术
在过去的几十年期间,已经对利用雷达系统检测目标进行了许多研究。知道与目标的距离、速度和角度有关的信息对于确定目标的精确位置是至关重要的。使天线机械旋转的方法已经被常规用于识别从目标反射的信号的到达方向(在下文中,DOA)。然而,利用相控阵天线,可以通过电导引相控阵天线的波束来估计入射信号的DOA。
在相控阵天线的情况下,用于处理大量阵列信号的概念和技术已经被建议用于同时估计多个入射信号的DOA。具体地,已经尝试使用高分辨率角估计算法利用均匀线性相控阵天线来识别准确的DOA。然而,由于到目前为止提出的高分辨率角估计算法对接收的信号的信噪比(下文中,SNR)非常敏感,因此需要对所接收的信号进行校准以估计准确的DOA。
[现有技术文献]
[专利文献]
(专利文献1)韩国未审查专利申请公开No.10-2011-0080218
发明内容
为了解决上述问题,根据本公开的实施方式的用于线性相控阵天线的空间插值方法和设备具有以下方面。
本公开的一方面在于提供一种能够通过利用线性相控阵天线保证精确的DOA的估计来确定目标的精确位置的用于线性相控阵天线的空间插值方法和装置。
本公开的各方面不限于上述所提到的,并且本领域普通技术人员应当从下面的描述清楚地理解其它未提及的方面。
根据本公开的实施方式,一种用于线性相控阵天线的空间插值方法涉及一种用于包括多个发射天线单元和多个接收天线元件的线性相控阵天线的空间插值方法,并且包括以下步骤:步骤1,在该步骤中,所述多个接收天线元件接收从目标反射的反射波;步骤2,在该步骤中,使用角估计算法来估计所述反射波入射到所述多个接收天线元件上的入射角;第3步,在该步骤中,从所述多个接收天线元件当中选择不良状况(bad-conditioned)天线元件;以及步骤4,在该步骤中,对不良状况天线元件的接收信号进行补偿并且使用角估计算法来重新估计所述反射波入射到所述多个接收天线元件上的入射角。
优选地,步骤3包括以下步骤:步骤3-1,在该步骤中,计算由所述多个接收天线元件中的每一个接收的反射波的接收信号功率;步骤3-2,在该步骤中,将由所述多个接收天线元件中的每一个接收的反射波的接收信号功率与阈值进行比较;以及步骤3-3,在该步骤中,基于在步骤3-2中的比较的结果从所述多个接收天线元件中选择不良状况天线单元。
优选地,步骤4包括以下步骤:步骤4-1,在该步骤中,基于另一接收天线元件的接收信号来对不良状况天线元件的接收信号进行补偿;以及步骤4-2,在该步骤中,通过将经补偿的不良状况天线元件的接收信号反映到估计中,使用角估计算法来重新估计反射波入射到所述多个接收天线元件上的入射角。
优选地,所述角估计算法是多信号分类(MUSIC)算法或者通过旋转不变性技术(ESPRIT)的参数估计。
根据本公开的实施方式,一种用于线性相控阵天线的空间插值装置包括:多个接收天线元件,其接收从目标反射的反射波;入射角估计器,其估计所述反射波入射到所述多个接收天线元件上的入射角;检测器,其基于所述反射波的入射角来从所述多个接收天线元件当中选择不良状况天线元件;以及补偿器,其用于基于由所述检测器选择的结果来对所述不良状况天线元件的接收信号进行补偿,其中,所述入射角估计器基于所述不良状况天线元件的由所述补偿器补偿后的接收信号来重新估计所述反射波入射到所述多个接收天线元件的入射角。
优选地,所述检测器基于由所述多个接收天线元件中的每一个接收的反射波的接收信号功率来选择不良状况天线元件。
优选地,所述检测器包括:功率计算器,其用于计算由所述多个接收天线元件中的每一个接收的反射波的接收信号功率;比较器,其将所计算的接收信号功率与阈值进行比较;以及确定器,其用于基于由所述比较器比较的结果来确定所述多个接收天线元件当中的不良状况天线元件。
优选地,所述补偿器基于另一接收天线元件的接收信号来对所述不良状况天线元件的接收信号进行补偿。
优选地,所述入射角估计器使用角估计算法来估计和重新估计反射波入射到所述多个接收天线元件上的入射角。
优选地,所述角估计算法是MUSIC算法或ESPRIT。
附图说明
通过结合附图对本发明的示例性实施方式进行详细描述,本发明的以上和其它目的、特征和优点对于本领域普通技术人员而言将是更显而易见的,其中:
图1是按时间序列例示了根据本公开的实施方式的用于线性相控阵天线的空间插值方法的流程图;
图2和图3是图1的方法的一些步骤被进一步细分并详细说明的流程图;以及
图4是简要地例示了根据本公开的实施方式的用于线性相控阵天线的空间插值装置的框图。
具体实施方式
将参照附图详细地描述本公开的示例性实施方式。相同或相似的元件将被赋予相同的附图标记,并且将省略重复的描述。
此外,当公知的相关技术的详细描述被认为在描述本公开中使本发明的要旨模糊不清时,其详细描述将被省略。另外,附图仅用于帮助容易地理解本公开的精神,并且本公开的精神不应被解释为受附图限制。
在下文中,将参照图1至图3对根据实施方式的用于线性相控阵天线的空间插值方法进行说明。图1是按时间序列例示了根据本公开的实施方式的用于线性相控阵天线的空间插值方法的流程图,并且图2和图3是图1的方法的一些步骤被进一步细分并详细说明的流程图。
根据本公开的实施方式,一种用于线性相控阵天线的空间插值方法涉及一种用于包括多个发射天线单元110和多个接收天线元件120在内的线性相控阵天线100的空间插值方法,并且如图1中所例示的,该空间插值方法包括以下步骤:步骤1(S100),在该步骤中,多个接收天线元件120接收反射波;步骤2(S200),在该步骤中,估计该反射波的入射角;步骤3(S300),在该步骤中,选择不良状况天线元件;以及步骤4(S400),在该步骤中,重新估计该反射波的入射角。
具体地,在步骤1(S100)中,多个发射天线元件110向目标10发射信号,所发射的信号被目标10反射,并且接收天线元件120接收从目标10反射的反射波。此外,关于步骤1(S100)至步骤4(S400)的详细描述,将首先在下文描述相控阵天线的接收信号的模型。
假定L个平面波利用已知的载波频率fc从L个方向θ12,...,θL入射到以相等间隔d线性间隔开的N个接收天线元件120上。θi(i=1,2,...,L)从面对阵列天线的前轴线开始被限定。此外,假定存在窄带长距离源。接下来,阵列天线在时刻k所接收的信号矢量被表示为下面的式1。
[式1]
x(k)=As(k)+n(k)=[x1(k),x2(k),…,xN(k)]T
在上面的式1中,(·)T是转置运算符,并且A=[a(θ1),a(θ2),...,a(θL)]是由导引矢量形成的导引矩阵,a(θi)(i=1,2,...,L)在下面的式2中被示出。
[式2]
Figure BDA0001106885890000041
λ是与载波频率fc对应的波长。此外,s(k)=[s1(k),s2(k),...,sN(k)]T是入射信号矢量,并且这里,si(K)(i=1,2,...,L)是在时刻k来自第i个源的入射信号的复数振幅。假定振幅是零均值复高斯的并且针对各个振幅样本是独立的。假定n(k)是由[n1(k),n2(k),...,nN(k)]T形成的噪声矢量,并且该噪声矢量也是零均值复高斯的并独立于每个样本。相关矩阵被提供为σi 2IN(i=1,2,…,N),并且这里,IN是N*N单位矩阵。此外,假定来自入射信号矢量和噪声矢量的样本彼此独立。
在步骤1(S100)之后,执行步骤2(S200),在该步骤中,对多个接收的反射波的入射角进行估计。角估计算法被用于入射角的估计。角估计算法包括多信号分类(MUSIC)算法、经由旋转不变性技术(ESPRIT)的参数估计等,并且下面将主要描述使用MUSIC算法来估计反射波的入射角的实施方式。
为了使用MUSIC算法,首先,应该对上面的式1中的接收信号矢量的相关矩阵进行运算。相关矩阵由下面的式3限定。
[式3]
Rxx=E[x(k)x(k)H]
=AE[s(k)s(k)H]AHN
=ARssAHN
在式3中,E[·]是随机过程的总体平均,并且(·)H是复共轭转置运算符。
RXX的秩遵循ARSSAH的秩,并且该秩为L。RXX的特征值可以按照在下面的式4中表达的递减顺序排列。
[式4]
λ1>λ2>…>λL>λL+1>…>λN
前L个特征值与信号子空间相关,而剩余的特征值与噪声子空间相关。当各个特征值的归一化特征矢量被给定为vi(i=1,2,…,N)时,MUSIC算法的伪谱P(θ)被限定为下面的式5。
[式5]
Figure BDA0001106885890000051
这里,分母中的
Figure BDA0001106885890000052
形成噪声子空间,并且与入射信号正交。入射信号的估计的DOA由θ的使P(θ)的分母接近零的值来确定。然而,实际上不可能知道对于信号和噪声的准确统计,因此难以运算式3的总体平均值。在这种情况下,当假定过程是遍历方法时,时间平均相关矩阵在存在噪音的同时通过有限的K次测量来计算,并且被限定如下。
[式6]
Figure BDA0001106885890000053
应用MUSIC算法时的另一个显著问题是如何确定入射信号的数目。然而,由于不可能精确地知道入射平面波的数目,因此通过对相关矩阵的特征值的大小进行比较来估计入射平面波的数目,或者根据经验设置的设置值被用作入射平面波的数目。这里,设置值可以根据视场(FOV)的大小而改变,FOV是DOA被期望知道的范围。
以上所述的与式6相关的假定被应用到步骤2(S200),并且MUSIC算法被应用到最初接收的入射到多个接收天线元件120上的信号。
此外,如上所述,常规的MUSIC算法的性能对阵列的接收信号的SNR敏感。当所有元件的接收信号的SNR高并且几乎彼此相等时,可以充分地操作MUSIC算法。当接收信号的SNR不一致时,存在MUSIC算法的分解性能劣化的问题。因此,需要空间插值以提高分辨率。因此,在根据本公开的实施方式的用于线性相控阵天线的空间插值方法中,在步骤2(S200)之后执行步骤3(S300)以及步骤4(S400),在步骤3(S300)中,从多个接收天线元件120中选择不良状况天线元件,在步骤4(S400)中,对该不良状况天线元件的接收信号进行补偿并且使用角估计算法重新估计反射波的入射角。此后,将参照图2和图3详细地描述步骤3(S300)和步骤4(S400)。
如图2中所例示,步骤3(S300)包括:步骤3-1(S310),在该步骤中,计算反射波的接收信号功率;步骤3-2(S320),在该步骤中,将反射波的接收信号功率与阈值进行比较;以及步骤3-3(S330),在该步骤中,选择不良状况天线元件。
在步骤2(S200)中,首先测量作为反射波的接收信号的SNR。考虑到难以从所接收的信号中将期望信号和噪声分离的事实,利用步骤3-1(S310)中的接收信号功率来替换SNR值。
因为假定入射信号来自接收等量功率的不相关的源,所以接收信号功率主要受噪声功率水平影响。因此,SNR值在这种情况下与接收信号的功率成反比。例如,当预定天线元件具有高的噪底(noise floor)时,测量的接收信号功率高,而测得的SNR低。由于这个原因,接收信号功率被用作替代SNR值的替代参数。当第i(i=1,2,...,N)个天线元件关于K个时间样本的平均接收信号功率被表达为Pi时,接收信号功率矢量可以被限定为下面的式7。
[式7]
当Pi之间几乎不存在差异时,算法可以针对两个相邻的入射信号估计两个单独的DOA(即,
Figure BDA0001106885890000071
Figure BDA0001106885890000072
)。当Pi之间存在差异时,两个信号被整合成单个信号,并且该算法仅估计一个DOA
Figure BDA0001106885890000073
在步骤后3-1(S310)之后,执行步骤3-2(S320)以及步骤3-3(S330),在步骤3-2(S320)中,针对仅一个估计的DOA将由多个接收天线元件120中的每一个接收的反射波的接收信号功率与阈值进行比较,在步骤3-3(S330)中,基于步骤3-2(S320)中的比较的结果从多个接收天线元件120中选择不良状况天线元件。
步骤3-2(S320)中的阈值可以是预定值,或者可以由与最大值的比率来确定。在下文中,将描述使用对分别由多个接收天线元件120接收的反射波的接收信号功率求平均来执行步骤3-2(S320)和步骤3-3(S330)的示例。
首先,针对对所有元件的Pi求平均,
Figure BDA0001106885890000075
是按照以下的式8的一个估计的DOA。
[式8]
Figure BDA0001106885890000076
然后,从多个接收天线元件120当中选择不良状况天线元件。可以通过选择满足以下的式9中所示的条件的元件来选择不良状况天线元件。
[式9]
Figure BDA0001106885890000077
在上面的式9中,i*被限定为不良状况天线元件,并且Pthr是指应该被设置的阈值功率值。此外,Pthr可以是随机设置的值或比率,或者还可以使用与最大功率值的比率来设置。
此外,虽然上面已经描述了使用式8和式9来选择不良状况天线元件的实施方式,但是也可以使用诸如中值和平均最大值等这样的其它各种统计方法来选择不良状况天线元件。
此外,时域中的接收信号可以利用快速傅里叶变换(FFT)被变换成频域的接收信号,并且还可以通过对接收信号功率的峰值进行比较来选择不良状况天线元件。
具体地,可以选择由多个天线元件接收的具有接收信号功率的峰值的频带,并且,通过对该频带中的接收信号功率的峰值彼此进行比较,可以选择接收到具有这些峰值当中的相对最小峰值的接收信号功率的天线元件作为不良状况天线元件。
与此不同,时域的接收信号可以利用FFT变换成频域的接收信号,可以针对各个天线元件对峰值附近的周边值的最低水平与接收信号功率的峰值之间的差进行比较,并且可以选择具有最小差的天线元件作为不良状况天线元件,或者可以预置预定的基准值并且可以选择周边值的最低水平和峰值之间的差为该基准值或更小的天线元件作为不良状况天线元件。这里,可以利用对除峰值以外的周边值求平均来计算周边值的最低水平。
当如上面那样选择不良状况天线元件时,执行重新估计反射波的入射角的步骤4(S400)。如图3中所例示,步骤4(S400)可以被划分成:步骤4-1(S410),其中对不良状况天线元件的接收信号进行补偿;以及步骤4-2(S420),其中基于经补偿的不良状况天线元件的接收信号来重新估计反射波的入射角。
具体地,步骤4-1(S410)基于另一接收天线元件的接收信号对不良状况天线元件的接收信号进行补偿。在下文中,将描述基于与不良状况天线元件相邻的两个天线元件的接收信号来对不良状况天线元件的接收信号进行补偿的示例。
来自第i*元件的K个时间采样的接收信号被表达如下。
[式10]
Figure BDA0001106885890000081
这里,如参照式1描述的,
Figure BDA0001106885890000082
是由第i*元件在时刻k接收到的接收信号。信号矢量
Figure BDA0001106885890000083
具有比剩余的信号矢量低的SNR。因此,当
Figure BDA0001106885890000084
被修改时,MUSIC算法的分辨率性能提高。通过使用从第(i*-1)元件和第(i*+1)元件接收到的信号矢量,
Figure BDA0001106885890000085
可以被修改为如下面的式11所示的
[式11]
Figure BDA0001106885890000087
这里,
Figure BDA0001106885890000088
是在步骤2(S200)中估计的DOA。当两个信号s1(k)和s2(k)从方向θ1和θ2入射到阵列上时,
Figure BDA0001106885890000089
Figure BDA00011068858900000810
被表达为下面的式12。
[式12]
Figure BDA0001106885890000091
Figure BDA0001106885890000092
Figure BDA0001106885890000093
因此,当两个信号源彼此非常接近时,建立根据式13的关系。
[式13]
Figure BDA0001106885890000094
依据式13,式11适合于对接收信号进行插值,并且因此,
Figure BDA0001106885890000095
被变换为
Figure BDA0001106885890000096
如图3中例示,在如上所述对不良状况天线元件的接收信号进行补偿之后,执行步骤4-2(S420),其中通过将经补偿的不良状况天线元件的接收信号反映到估计中,使用角估计算法来重新估计反射波入射到多个接收天线元件120的入射角。在下文中,将描述如在步骤2(S200)那样应用MUSIC算法来重新估计当步骤4-2(S420)被执行时反射波的入射角的示例。
首先,在步骤4-1(S410)中新生成的信号矢量被用来生成如同下面的式14表达的新的相关矩阵
Figure BDA0001106885890000098
[式14]
Figure BDA0001106885890000099
通过使用新生成的相关矩阵
Figure BDA00011068858900000910
以应用MUSIC算法,来对特征值分解进行处理。然后,使用式5中的MUSIC算法的伪谱来估计新的DOA。另一方面,当所选择的第i*元件是阵列的第一个元件或最后一个元件时,优选地使用一个相邻的接收信号矢量。
在下文中,将参照图4描述根据本公开的实施方式的用于线性相控阵天线的空间插值装置,同时将省略与根据本公开的实施方式的用于线性相控阵天线的空间插值方法的以上描述重复的内容。
如图4中例示,根据本公开的实施方式的用于线性相控阵天线的空间插值装置包括多个接收天线元件120、入射角估计器200、检测器300和补偿器400。
多个接收天线元件120是被如下配置的:该配置与发射天线单元110中的多个线性相控阵天线100一起被设置,并且执行接收由从多个发射天线元件110发射的信号产生的正被目标10反射的反射波的功能。
入射角估计器200是被如下配置的:该配置估计反射波入射到接收天线元件120上的入射角,并且通过使用角估计算法来执行估计反射波入射到接收天线元件120上的入射角的功能。MUSIC算法可以被用作角度估计算法。
检测器300是被如下配置的:该配置执行基于由入射角估计器200估计的反射波的入射角在多个接收天线元件120当中选择不良状况天线元件的功能。具体地,检测器300可以基于由多个接收天线元件120中的每一个接收的反射波的接收信号功率来选择不良状况天线元件。为了执行该功能,检测器300还可以被配置有功率计算器310、比较器320和确定器330。功率计算器310计算由多个接收天线元件120中的每一个接收到的反射波的接收信号功率,比较器320将所计算的接收信号功率与阈值进行比较,并且确定器330基于比较器320的比较的结果来确定多个接收天线元件120当中的不良状况天线元件。
补偿器400对所选择的不良状况天线元件的接收信号进行补偿,并且因为已经给出了根据本公开的实施方式的用于线性相控阵天线的空间插值方法的描述,所以将省略对其的详细描述。
此外,在根据本公开的实施方式的用于线性相控阵天线的空间插值设备中,入射角估计器200通过将不良状况天线元件的由补偿器400补偿的接收信号反映到估计中来重新估计反射波入射到接收天线元件120上的入射角。如上所述,可以应用诸如MUSIC算法这样的高分辨率角估计算法来估计反射波的入射角。
因此,当入射角估计器200第一次估计反射波的入射角时,检测器300使用所估计的入射角来选择不良状况天线元件,补偿器400对不良状况天线元件的接收信号进行补偿,并且入射角估计器200通过将经补偿的不良状况天线元件的接收信号反映到估计中来第二次估计反射波的入射角,入射信号的入射角可以利用精细且高的分辨率来计算。
根据本公开的实施方式,用于线性相控阵天线的空间插值方法和设备通过使用角估计算法来估计反射波的入射角,当正在估计反射波的入射角时选择处于不良状况的天线元件,并且对所选择的天线元件的接收信号进行插值以重新估计反射波的入射角,因此提高DOA估计算法的分辨率性能。
本公开的有利效果不限于以上提及的效果,并且本领域普通技术人员应该从上面的描述清楚地理解其它未提及的有利效果。
本文中所描述的实施方式和附图仅仅示意性地描述了属于本公开的技术精神的一部分。因此,由于本文中所公开的实施方式是用于描述而不是限制本公开的技术精神,应该不言而喻的是,本公开的技术精神的范围不受所述实施方式限制。可以由本领域普通技术人员容易地在包含在本公开的说明书和附图中的技术精神的范围内推断出的改进的示例和具体实施方式都应该被解释为属于本公开的范围。
[标号说明]
10:目标 100:线性相控阵天线
110:发射天线元件 120:接收天线单元
200:入射角估计器 300:检测器
310:功率计算器 320:比较器
330:确定器 400:补偿器
相关申请的交叉引用
本申请要求于2015年9月1日申请的韩国专利申请No.10-2015-0123802的优先权和利益,该韩国专利申请的公开内容通过引用方式被完整地并入到本文中。

Claims (16)

1.一种用于线性相控阵天线(100)的空间插值方法,该线性相控阵天线(100)包括多个发射天线元件(110)和多个接收天线元件(120),该空间插值方法包括以下步骤:
由所述多个接收天线元件(120)接收从目标(10)反射的反射波;
估计所述反射波入射到所述多个接收天线元件(120)上的入射角;
通过使用由所述多个接收天线元件(120)中的每一个接收的反射波的接收信号功率来从所述多个接收天线元件(120)当中选择不良状况天线元件;以及
通过使用另一接收天线元件的接收信号来对所述不良状况天线元件的接收信号进行补偿,并且通过将经补偿的所述不良状况天线元件的接收信号反映到估计中来重新估计所述反射波入射到所述多个接收天线元件(120)上的入射角,
其中,重新估计所述反射波的入射角的步骤包括:对所述不良状况天线元件的接收信号进行插值以重新估计所述反射波的入射角。
2.根据权利要求1所述的空间插值方法,其中,在选择所述不良状况天线元件的步骤中,将由所述多个接收天线元件(120)接收的反射波的接收信号功率的平均值与阈值进行比较,并且选择所述反射波的接收信号功率的平均值小于所述阈值的接收天线元件作为所述不良状况天线元件。
3.根据权利要求1所述的空间插值方法,其中,在选择所述不良状况天线元件的步骤中,将由所述多个接收天线元件(120)接收的反射波的接收信号功率的峰值彼此进行比较,并且选择具有所述接收信号功率的最小峰值的接收天线元件作为所述不良状况天线元件。
4.根据权利要求1所述的空间插值方法,其中,在选择所述不良状况天线元件的步骤中,将由所述多个接收天线元件(120)接收的接收信号功率的峰值附近的周边值的最低水平与所述峰值之间的差彼此进行比较,并且选择具有最小差的天线元件作为所述不良状况天线元件。
5.根据权利要求1所述的空间插值方法,其中,在选择所述不良状况天线元件的步骤中,预置预定的基准值,选择由所述多个接收天线元件(120)接收的接收信号功率的峰值附近的周边值的最低水平与所述峰值之间的差为所述基准值或小于所述基准值的天线元件作为所述不良状况天线元件。
6.根据权利要求1所述的空间插值方法,其中,在重新估计所述反射波的入射角的步骤中,通过使用与所述不良状况天线元件相邻的接收天线元件的接收信号来对所述不良状况天线元件的接收信号进行补偿。
7.根据权利要求1所述的空间插值方法,其中,在估计所述反射波的入射角的步骤和重新估计所述反射波的入射角的步骤中,使用角估计算法。
8.根据权利要求7所述的空间插值方法,其中,所述角估计算法是多信号分类MUSIC算法或者通过旋转不变性技术ESPRIT的参数估计。
9.一种用于线性相控阵天线的空间插值装置,该空间插值装置包括:
多个接收天线元件(120),所述多个接收天线元件(120)被配置为接收从目标(10)反射的反射波;
入射角估计器(200),所述入射角估计器(200)被配置为估计所述反射波入射到所述多个接收天线元件(120)上的入射角;
检测器(300),所述检测器(300)被配置为通过使用由所述多个接收天线元件(120)中的每一个接收的反射波的接收信号功率来在所述多个接收天线元件(120)当中选择不良状况天线元件;以及
补偿器(400),所述补偿器(400)被配置为使用另一接收天线元件的接收信号来对所述不良状况天线元件的接收信号进行补偿,
其中,所述入射角估计器(200)通过将经补偿的所述不良状况天线元件的接收信号反映到估计中来重新估计所述反射波入射到所述多个接收天线元件(120)上的入射角,并且
其中,所述入射角估计器(200)对所述不良状况天线元件的接收信号进行插值以重新估计所述反射波的入射角。
10.根据权利要求9所述的空间插值装置,其中,所述检测器(300)将由所述多个接收天线元件(120)接收的反射波的接收信号功率的平均值与阈值进行比较,并且选择所述反射波的接收信号功率的平均值小于所述阈值的接收天线元件作为所述不良状况天线元件。
11.根据权利要求9所述的空间插值装置,其中,所述检测器(300)将由所述多个接收天线元件(120)接收的反射波的接收信号功率的峰值彼此进行比较,并且选择具有所述接收信号功率的最小峰值的接收天线元件作为所述不良状况天线元件。
12.根据权利要求9所述的空间插值装置,其中,所述检测器(300)将由所述多个接收天线元件(120)接收的接收信号功率的峰值附近的周边值的最低水平与所述峰值之间的差彼此进行比较,并且选择具有最小差的天线元件作为所述不良状况天线元件。
13.根据权利要求9所述的空间插值装置,其中,所述检测器(300)预置预定的基准值,选择由所述多个接收天线元件(120)接收的接收信号功率的峰值附近的周边值的最低水平与所述峰值之间的差为所述基准值或小于所述基准值的天线元件作为所述不良状况天线元件。
14.根据权利要求9所述的空间插值装置,其中,所述补偿器(400)通过使用与所述不良状况天线元件相邻的接收天线元件的接收信号来对所述不良状况天线元件的接收信号进行补偿。
15.根据权利要求9所述的空间插值装置,其中,所述入射角估计器(200)通过使用角估计算法来估计或重新估计所述反射波入射到所述多个接收天线元件(120)上的入射角。
16.根据权利要求15所述的空间插值装置,其中,所述角估计算法是多信号分类MUSIC算法或者通过旋转不变性技术ESPRIT的参数估计。
CN201610796509.0A 2015-09-01 2016-08-31 用于线性相控阵天线的空间插值方法和设备 Active CN106486769B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR10-2015-0123802 2015-09-01
KR1020150123802A KR102422396B1 (ko) 2015-09-01 2015-09-01 선형 위상 어레이 안테나의 공간 보간 방법 및 보간 장치

Publications (2)

Publication Number Publication Date
CN106486769A CN106486769A (zh) 2017-03-08
CN106486769B true CN106486769B (zh) 2020-01-07

Family

ID=58010856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610796509.0A Active CN106486769B (zh) 2015-09-01 2016-08-31 用于线性相控阵天线的空间插值方法和设备

Country Status (4)

Country Link
US (1) US10215843B2 (zh)
KR (1) KR102422396B1 (zh)
CN (1) CN106486769B (zh)
DE (1) DE102016010554A1 (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107329110B (zh) * 2017-08-24 2019-08-30 浙江大学 基于稀疏阵列直接内插的波达方向估计方法
CN107958312B (zh) * 2017-12-12 2021-08-31 国网湖南省电力有限公司 基于反演算法的输电线路舞动预测方法、系统及存储介质
KR102323655B1 (ko) * 2018-01-08 2021-11-09 현대모비스 주식회사 가상 수신신호 생성을 이용한 도래각 추정 방법 및 장치
KR102182293B1 (ko) * 2019-03-11 2020-11-24 현대모비스 주식회사 Mimo 시스템에서 도래각 추정 장치 및 방법
CN110954877B (zh) * 2019-11-05 2022-04-26 中国电波传播研究所(中国电子科技集团公司第二十二研究所) 一种车载雷达散射特性测量区域地形生成方法
DE102020204759A1 (de) 2020-04-15 2021-10-21 Robert Bosch Gesellschaft mit beschränkter Haftung Radarempfängereinrichtung, Radarempfängervorrichtung und Radarempfangsverfahren
CN112798876B (zh) * 2021-03-22 2021-06-22 上海莱天通信技术有限公司 一种利用插值算法的天线近场测量方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1913402A (zh) * 2005-08-11 2007-02-14 中兴通讯股份有限公司 一种智能天线故障检测的方法
CN101904051A (zh) * 2007-12-17 2010-12-01 波音公司 相控阵列天线的准确自校准
CN103760519A (zh) * 2014-01-24 2014-04-30 深圳大学 高分辨率doa估计方法及系统
CN104303072A (zh) * 2012-05-16 2015-01-21 株式会社电装 雷达装置以及雷达装置的入射波处理方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5657023A (en) * 1996-05-02 1997-08-12 Hughes Electronics Self-phase up of array antennas with non-uniform element mutual coupling and arbitrary lattice orientation
JP2003218621A (ja) * 2002-01-21 2003-07-31 Nec Corp アレーアンテナの校正装置及び校正方法
JP2006003097A (ja) * 2004-06-15 2006-01-05 Fujitsu Ten Ltd レーダ装置
US6992622B1 (en) * 2004-10-15 2006-01-31 Interdigital Technology Corporation Wireless communication method and antenna system for determining direction of arrival information to form a three-dimensional beam used by a transceiver
US7576686B2 (en) * 2006-08-07 2009-08-18 Garmin International, Inc. Method and system for calibrating an antenna array for an aircraft surveillance system
US7477192B1 (en) * 2007-02-22 2009-01-13 L-3 Communications Titan Corporation Direction finding system and method
JP4484892B2 (ja) * 2007-03-14 2010-06-16 三菱電機株式会社 車載用レーダ装置
US7522096B2 (en) * 2007-04-09 2009-04-21 Honeywell International Inc Method for phase calibrating antennas in a radar system
JP4978323B2 (ja) 2007-06-13 2012-07-18 株式会社日立製作所 Ofdmセルラ無線通信基地局
US8154452B2 (en) * 2009-07-08 2012-04-10 Raytheon Company Method and apparatus for phased array antenna field recalibration
KR101137038B1 (ko) 2010-01-05 2012-04-19 주식회사 만도 레이더 장치, 안테나 장치 및 데이터 획득 방법
WO2011112083A1 (en) * 2010-03-08 2011-09-15 Nederlandse Organisatie Voor Toegepast- Natuurwetenschappelijk Onderzoek Tno A method of compensating sub-array or element failure in a phased array radar system, a phased array radar system and a computer program product
DE102011055693A1 (de) * 2011-11-24 2013-05-29 Hella Kgaa Hueck & Co. Radargerät und Verfahren zur Erkennung eines Ausfalls eines Empfangskanals eines Radargerätes
CN105026963B (zh) 2013-03-07 2018-01-23 东丽株式会社 黑矩阵基板
US9291659B2 (en) * 2013-12-19 2016-03-22 Ford Global Technologies, Llc Antenna blockage detection

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1913402A (zh) * 2005-08-11 2007-02-14 中兴通讯股份有限公司 一种智能天线故障检测的方法
CN101904051A (zh) * 2007-12-17 2010-12-01 波音公司 相控阵列天线的准确自校准
CN104303072A (zh) * 2012-05-16 2015-01-21 株式会社电装 雷达装置以及雷达装置的入射波处理方法
CN103760519A (zh) * 2014-01-24 2014-04-30 深圳大学 高分辨率doa估计方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"均匀圆阵部分阵元失效情况下的DOA估计方法";张宏谋;《电子设计工程》;20131031;全文 *
"基于线性预测和最小二乘估计的缺陷阵波束形成";崔琳;《系统仿真学报》;20091231;第7782-7785页 *

Also Published As

Publication number Publication date
KR20170027375A (ko) 2017-03-10
KR102422396B1 (ko) 2022-07-20
US20170059694A1 (en) 2017-03-02
DE102016010554A1 (de) 2017-03-02
CN106486769A (zh) 2017-03-08
US10215843B2 (en) 2019-02-26

Similar Documents

Publication Publication Date Title
CN106486769B (zh) 用于线性相控阵天线的空间插值方法和设备
CN107064861B (zh) 用于估计到达角的设备以及用于波束成形的设备
CN111373282B (zh) 用于fmcw雷达系统的雷达处理链
CN110018452B (zh) 使用虚拟接收信号的生成来估计到达方向的方法和设备
US7847733B2 (en) Direction-of-arrival estimating device and program
JP5600866B2 (ja) 探知測距装置および探知測距方法
US7379020B2 (en) Apparatus and method for estimating direction of arrival of radio wave
JP7117064B2 (ja) レーダ装置、レーダシステム
CN109407093B (zh) 解决宽孔径雷达到达角模糊度的多普勒测量
US9470782B2 (en) Method and apparatus for increasing angular resolution in an automotive radar system
EP3324205B1 (en) Decentralised radar system
JP4972852B2 (ja) レーダ装置
US10732273B2 (en) Radar device for vehicle and method for estimating angle of target using same
JP6362262B2 (ja) 角度推定装置および角度推定方法
KR101468548B1 (ko) 레이더 장치 및 레이더 장치에서의 방위각 추정 방법
US10663574B2 (en) Angle estimating method and radar system
EP4014061A1 (en) Radio station for client localization in multipath indoor environment
JP2010223895A (ja) レーダ装置
Lee et al. Enhanced performance of MUSIC algorithm using spatial interpolation in automotive FMCW radar systems
US11520004B2 (en) System and method for generating point cloud data in a radar based object detection
JP5025170B2 (ja) 到来波数検出装置
JP4275425B2 (ja) レーダ装置
JP3928551B2 (ja) 到来波推定装置及び到来波推定方法、並びに波源位置推定装置
WO2020148802A1 (ja) ビーム形成装置、レーダ装置及びビーム形成方法
TWI808874B (zh) 用於交通工具的雷達系統及偵測方法

Legal Events

Date Code Title Description
C06 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
TR01 Transfer of patent right

Effective date of registration: 20220106

Address after: Gyeonggi Do, South Korea

Patentee after: Wandu Mobile System Co.,Ltd.

Patentee after: SNU R&DB FOUNDATION

Address before: Gyeonggi Do, South Korea

Patentee before: MANDO Corp.

Patentee before: SNU R&DB FOUNDATION

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220907

Address after: Inchon, Korea

Patentee after: Halla Creedong Electronics Co.,Ltd.

Patentee after: SNU R&DB FOUNDATION

Address before: Gyeonggi Do, South Korea

Patentee before: Wandu Mobile System Co.,Ltd.

Patentee before: SNU R&DB FOUNDATION

TR01 Transfer of patent right