CN105445709A - 一种稀布阵列近场无源定位幅相误差校正方法 - Google Patents

一种稀布阵列近场无源定位幅相误差校正方法 Download PDF

Info

Publication number
CN105445709A
CN105445709A CN201510770276.2A CN201510770276A CN105445709A CN 105445709 A CN105445709 A CN 105445709A CN 201510770276 A CN201510770276 A CN 201510770276A CN 105445709 A CN105445709 A CN 105445709A
Authority
CN
China
Prior art keywords
theta
phase error
amplitude
array
iteration
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.)
Granted
Application number
CN201510770276.2A
Other languages
English (en)
Other versions
CN105445709B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201510770276.2A priority Critical patent/CN105445709B/zh
Publication of CN105445709A publication Critical patent/CN105445709A/zh
Application granted granted Critical
Publication of CN105445709B publication Critical patent/CN105445709B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

该发明公开了一种稀布阵列近场无源定位幅相误差校正方法,涉及相控阵雷达技术领域。在目标位于近场的情况下,固定俯仰角将目标距离和方位角进行交替循环估计,直至两个参数都收敛至真实值,作为误差校正的初始位置;然后根据,对目标方位角进行一维搜索并用上述幅相误差自校正的方法对阵列幅相误差进行校正,估计出阵列幅相误差矩阵Γθ以及方位角的一次迭代值θ′;估计出阵列幅相误差矩阵Γr以及目标距离的一次迭代值;以此进行循环迭代,直至||Гθr||<ε2,即参数估计值收敛于真实值。本发明在一定的布阵范围内,阵元位置随机稀疏摆放,只需要较少的阵元就可达到相同的角度分辨率,在实际工程中有利于节约成本。

Description

一种稀布阵列近场无源定位幅相误差校正方法
技术领域
本发明涉及相控阵雷达技术领域,具体涉及无源定位的阵列误差校正方法。
背景技术
与常规均匀线阵相比,稀布阵列用相同的阵元使得阵列孔径增大,以较少的阵元达到较高的分辨率,从而简化阵列结构,降低造价。目前,稀布阵列天线在抗环境干扰的卫星接收天线,高频地面雷达天线和射电天文中的干涉阵列等领域中正在得到越来越广泛的应用。另外,波达方向(DirectionofArrival,DOA)估计在雷达等领域广泛应用,通常高分辨DOA估计阵列处理均基于理想阵列信号模型,但是具体在实际工程应用中,各通道间硬件参数的不同,使得目标回波信号在不同的阵元位置有不同的幅度和相位加权,即阵列中存在多通道幅相误差问题。
WeissA.J.和FriedlanderB.等人在“Eigenstructuremethodsfordirectionfindingwithsensorgainandphaseuncertainties”(Acoustics,Speech,andSignalProcessing,1988.ICASSP-88.,1988InternationalConferenceon.IEEE,1988:2681-2684)中提出了将阵列通道幅相误差参数与信号源方位交替联合迭代的阵列误差自校正算法(称为WF算法),以上文献只针对了在远场情况下对于MUSIC算法进行的幅相误差估计,但是在信源位于近场的情况下,信号的波前不能近似为平面波而要看作球面波,实现定位需要在距离和角度上进行联合二维搜索。本发明提供了一种基于迭代的稀布阵列近场无源定位幅相误差校正方法,在距离维和角度维上进行轮换迭代,即使在幅相误差较大时,该方法也能精确校正稀布阵列通道幅相误差。
发明内容
为了克服上述现有技术的不足,本发明提供了一种稀布阵列近场无源定位幅相误差校正方法,从而达到定位精度高,成本低廉的目的。
本发明提出的基于迭代的稀布阵列近场无源定位福相误差校正方法基本思想如下:在目标位于近场的情况下,固定俯仰角将目标距离r和方位角θ进行交替循环估计,直至两个参数都收敛至真实值,即:首先在存在误差的情况下粗略估计出目标距离r0和方位角θ0,作为误差校正的初始位置;然后根据r0,对目标方位角θ进行一维搜索并用上述幅相误差自校正的方法对阵列幅相误差进行校正,估计出阵列幅相误差矩阵Γθ以及方位角θ的一次迭代值θ′;最后根据θ′,对目标距离r进行一维搜索并用上述幅相误差自校正的方法对阵列幅相误差进行校正,估计出阵列幅相误差矩阵Γr以及目标距离r的一次迭代值r′;以此进行循环迭代,直至||Γθr||<ε2,即参数估计值收敛于真实值,迭代过程如附图2所示。因而本发明是一种稀布阵列近场无源定位幅相误差校正方法,该方法具体包括:
步骤1:获得定位场中信号接收机即定位阵元的个数及其位置信息;
步骤2:以第一个阵元为参考阵元,各阵元接收P个位于定位场中的目标信号源发送的信号;
步骤3:用xm(t)表示第m个阵元在第t个时刻接收信号的采样值,则天线阵列接收信号
x ( t ) = &lsqb; x 1 ( t ) , x 2 ( t ) , ... , x M ( t ) &rsqb; T = &Sigma; p = 1 P a ( r p , &theta; p ) s p ( t ) + n ( t ) , t = 1 , 2 , ... , L
其中,(·)T表示转置操作符,sp(t)表示第p个近场信号源,n(t)表示天线阵列接收的噪声矢量,L表示采样总数,a(rpp)表示第p个近场信号源导向矢量,rp表示第p个近场信号源到参考阵元的距离;
步骤4:通过接收天线阵列数据x(t)估算协方差矩阵(·)H表示对矩阵进行共轭转置计算;再估算协方差矩阵的逆矩阵R-1,(·)-1表示对矩阵进行求逆运算;
步骤5:采用最小方差波束形成算法估计出目标角度和距离的初始值 &theta; ^ ( 0 ) = ( &theta; ^ 1 ( 0 ) , &theta; ^ 2 ( 0 ) , ... , &theta; ^ P ( 0 ) ) , 令二维迭代次数N=0;
步骤6:首先在角度维上进行迭代,固定目标距离为则可以在角度维上进行一维幅相误差校正,设定角度维上迭代次数变量为k,从角度维上估计出的稀布阵列雷达幅相误差对角矩阵为初始值为
定义代价函数
J k = &Sigma; p = 1 P a H ( r p ( k ) , &theta; p ( k ) ) ( &Gamma; &theta; ( k ) ) H R - 1 &Gamma; &theta; ( k ) a ( r p ( k ) , &theta; p ( k ) )
其中表示第p个近场信号源的导向矢量,(·)H表示对矩阵进行共轭转置计算;增加迭代次数,直到Jk-1-Jk≤ε1结束迭代,获得角度估计值和幅相误差估计值为
步骤7:在距离维上进行迭代,固定目标角度为则可以在距离维上进行一维幅相误差校正;设定从距离维上估计出的稀布阵列雷达幅相误差对角矩阵初始值为进行与步骤7相同的一维迭代过程,可求得距离估计值和幅相误差估计值
步骤8:比较两个维度上估计出的幅相误差对角矩阵,若满足则停止迭代,其中ε2为根据实际情况预设的门限,即得到幅相误差的估计值否则,令N=N+1,转到步骤6,继续迭代,其中N为角度维和距离维共同迭代次数的变量。
进一步的,所述步骤6的具体步骤为:
步骤6.1:定义空间谱公式为:
p ( &theta; ^ ( k ) | &Gamma; ^ &theta; ( k ) , r ^ ( k ) ) = 1 a ^ H ( r ^ ( k ) , &theta; ^ ( k ) ) ( &Gamma; ^ &theta; ( k ) ) H R - 1 &Gamma; ^ &theta; ( k ) a ^ ( r ^ ( k ) , &theta; ^ ( k ) )
根据利用上述空间谱公式求取对应的P个最大值时的角度便可得到第k次迭代的P个波达方向的估计值;
步骤6.2:令w=[1,0,...,0]T,则满足约束条件则用拉格朗日乘子法,可求得幅相误差矩阵如下:
&Gamma; ^ &theta; ( k + 1 ) = d i a g { &delta; ^ ( k + 1 ) } &delta; ^ ( k + 1 ) = ( Q k ) - 1 w w T ( Q k ) - 1 w
其中, Q k = &Sigma; p = 1 P a ~ H ( r ^ p ( k ) , &theta; ^ p ( k ) ) R - 1 a ~ ( r p ( k ) , &theta; ^ p ( k ) ) , a ~ ( r ~ p ( k ) , &theta; ^ p ( k ) ) = d i a g { a ( r ^ p ( k ) , &theta; ^ p ( k ) ) } ;
步骤6.3判断是否收敛
J k = ( &delta; ^ ( k + 1 ) ) H Q k &delta; ^ ( k + 1 )
在满足Jk-1-Jk1,时,继续迭代,ε1为根据实际情况预设的门限值;否则,结束迭代;最后一次迭代得到角度估计值和幅相误差估计值为
与现有技术相比,本发明方法能够对近场目标进行阵列幅相误差的二维校正,并利用近场无源定位的二维特点,对估计出的误差矩阵进行修正并能实现准确估计。此外,本发明基于一种随机稀布阵列,即在一定的布阵范围内,阵元位置随机稀疏摆放,只需要较少的阵元就可达到相同的角度分辨率,在实际工程中有利于节约成本。
附图说明
图1为近场窄带信号无源定位模型;
图2为一种基于迭代的稀布阵列近场无源定位幅相误差校正方法流程图;
图3为存在阵列误差情况下,稀布阵列近场无源定位的方向图;
图4为存在阵列误差情况下,采用本发明的幅相误差校正方法,对阵列接收信号进行补偿后,稀布阵列近场无源定位的方向图;
图5是图3和图4中当方位角为-22.5°时,误差校正前后的切面对比图;
图6是图3和图4中当距离为1700m时,误差校正前后的切面对比图。
具体实施方式
下面结合附图对本发明进一步说明。
步骤1:在笛卡尔坐标系中的x轴上的[0,D]米范围内随机布置M个天线,形成一个稀布阵,其位置产生方式如下:在位置0以及位置D处各布置1个阵元且位置分别设为d1与dM,即d1=0,dM=D,固定阵列孔径。然后在(0,D)范围以均匀分布方式随机产生M-2个数,并从小到大进行排序,将剩下的M-2个阵元分别布置在M-2个随机数表示的位置上,并设阵元位置分别为d2,…,dm,…,dM-1
步骤2:以d1所在的位置的阵元为参考阵元,考虑P个近场不相关目标信号源入射到步骤1产生的阵列,对任意的p=1,2,…,P,用(rpp)描述第p个信号源的位置,其中,如图1所示,rp表示第p个信号源到参考阵元的距离,θp表示第p个信号源相对于参考阵元与z轴的夹角;此外,P个近场不相关目标信号源要满足如下条件;
近场条件:0<rp≤2D2p,对任意的p=1,2,…,P,其中,λp表示第p个信号源的波长;
不相关条件:
r = = E { &lsqb; s i - E { s i } &rsqb; &lsqb; s q - E { s q } &rsqb; } E { &lsqb; s i - E { s i } &rsqb; 2 } E { &lsqb; s q - E { s q } &rsqb; 2 } = 0 i = 1 , 2 , ... , P q = 1 , 2 , ... , P i &NotEqual; q
其中,E[·]表示求期望,si表示第i个近场源信号,sq表示第q个近场源信号,r表示相关系数,对任意的i=1,2,…,Pq=1,2,…,P且i≠q,r=0表示si与sq不相关。
步骤3:则由几何知识可得第p个信号源到第m个阵元的距离rmp可表示为
r m p = ( r p sin&theta; p - d m ) 2 + ( r p cos&theta; p ) 2
则第p个信号源到第m个阵元与第p个信号源到参考阵元的距离差△rmp可以表示为如下公式
&Delta;r m p = ( r p sin&theta; p - d m ) 2 + ( r p cos&theta; p ) 2 - r p = r p ( 1 + d m 2 r p 2 - 2 d m sin&theta; p r p ) 1 2 - r p
若当距离参数满足0.62(D3p)1/2<rp<2D2p,即第p个近场信源位于菲涅耳区,则距离差△rmp可由泰勒级数展开近似表示为
&Delta;r m p = - d m sin&theta; p + d m 2 cos 2 &theta; p 2 r p
则第m个阵元在第t个采样时刻接收的数据表示为
x m ( t ) = &Sigma; p = 1 P r p r m p s p ( t ) e - j 2 &pi; &lambda; p ( - d m sin&theta; p + d m 2 cos 2 &theta; p 2 r p ) + n ( t ) , t = 1 , 2 , ... , L
其中,rmp表示第p个信号源到第m个阵元的距离,n(t)表示t时刻噪声信号,L表示总的采样数。
将M个阵元在t时刻采样的数据排成一个列矢量x(t),可以表示为
x ( t ) = &lsqb; x 1 ( t ) , x 2 ( t ) , ... , x M ( t ) &rsqb; T = &Sigma; p = 1 P a ( r p , &theta; p ) s p ( t ) + n ( t ) ,
其中,n(t)是独立的高斯白噪声信号矢量,a(rpp)表示第p个近场信号源导向矢量,其具体表达式为
a ( r p , &theta; p ) = &lsqb; 1 , r p r 2 p e - j 2 &pi; &lambda; p ( - d 2 sin&theta; p + d 2 2 cos 2 &theta; p 2 r p ) , ... , r p r M p e - j 2 &pi; &lambda; p ( - d M sin&theta; p + d M 2 cos 2 &theta; p 2 r p ) &rsqb; T
步骤4:估计接收数据的协方差矩阵可以表示为
R ^ = 1 L &Sigma; t = 1 L x ( t ) x H ( t )
计算协方差矩阵的逆其中(·)-1表示对矩阵求逆。
步骤5:设定迭代次数变量为k,第k次迭代中从距离维上估计出的稀布阵列雷达幅相误差对角矩阵为从角度维上估计出的稀布阵列雷达幅相误差对角矩阵为令初始值k=0,取其中Γ0表示雷达幅相误差对角矩阵初始值。
步骤6:进行迭代之前先结合最小方差波束形成算法估计出目标角度和距离的初始值
r ^ ( 0 ) = ( r ^ 1 ( 0 ) , r ^ 2 ( 0 ) , ... , r ^ p ( 0 ) ) , &theta; ^ ( 0 ) = ( &theta; ^ 1 ( 0 ) , &theta; ^ 2 ( 0 ) , ... , &theta; ^ P ( 0 ) ) .
最小方差波束形成算法的基本思想是:
通过对接收信号进行加权求和使得期望方向的功率最强,从而实现目标方向信号最大化接收并抑制其他方向干扰,其期望方向的加权矢量跟方向矢量相同。最小均方误差准则可以总结
为:
从数学上可以理解,最小方差算法就是在满足的所以权向量w条件下,找到满足使得阵列输出的功率wHRxw最小的权向量w。从物理意义上理解,上式保证了θ0方向上的信号正确接收,抑制了其它方向上的干扰信号。
对上式采用拉格朗日乘子法可得最优的权矢量为
则阵列的谱函数为
对谱函数进行搜索,得到的谱峰值就是来波方向以及位置。
步骤7:首先在角度维上进行迭代,固定目标距离为则可以在角度维上进行一维幅相误差校正。
定义代价函数
J k = &Sigma; p = 1 P a H ( &theta; p ( k ) , r ^ ( k ) ) ( &Gamma; &theta; ( k ) ) H R - 1 &Gamma; &theta; ( k ) a ( &theta; p ( k ) , r ^ ( k ) )
式中:“||·||”表示Frobenius范数。
7.1定义空间谱公式为:
p ( &theta; ^ ( k ) | &Gamma; ^ &theta; ( k ) , r ^ ( k ) ) = 1 a ^ H ( &theta; ^ ( k ) , r ^ ( k ) ) ( &Gamma; ^ &theta; ( k ) ) H R - 1 &Gamma; ^ &theta; ( k ) a ^ ( &theta; ^ ( k ) , r ^ ( k ) )
根据利用上述空间谱公式求取对应的P个最大值时的角度便可得到第k次迭代的P个波达方向的估计值。
7.2令w=[1,0,...,0]T,则满足约束条件则利用拉格朗日乘子法,可求得幅相误差矩阵如下:
&Gamma; ^ &theta; ( k + 1 ) = d i a g { &delta; ^ ( k + 1 ) } &delta; ^ ( k + 1 ) = ( Q k ) - 1 w w T ( Q k ) - 1 w
其中, Q k = &Sigma; p = 1 P a ~ H ( &theta; ^ p ( k ) ) R - 1 a ~ ( &theta; ^ p ( k ) ) .
7.3判断是否收敛
J k = ( &delta; ^ ( k + 1 ) ) H Q k &delta; ^ ( k + 1 )
在满足Jk-1-Jk>ε(ε为预设的门限值)时,继续迭代;否则,结束迭代。最后一次迭代得到的角度估计值和幅相误差估计值为
步骤8:在距离维上进行迭代,固定目标角度为即为步骤7中估计所得到的角度估计值,则可以在距离维上进行一维幅相误差校正。进行与步骤7相同的迭代过程,可求得距离估计值和幅相误差估计值
同样定义代价函数
J k = &Sigma; p = 1 P a H ( &theta; ^ ( k ) , r ^ p ( k ) ) ( &Gamma; R ( k ) ) H R - 1 &Gamma; r ( k ) a ( &theta; ^ ( k ) , r ^ p ( k ) )
式中:“||·||”表示Frobenius范数。
8.1定义空间谱公式为:
p ( r ^ ( k ) | &Gamma; ^ R ( k ) , &theta; ^ ( k ) ) = 1 a ^ H ( &theta; ^ ( k ) , r ^ ( k ) ) ( &Gamma; ^ R ( k ) ) H R - 1 &Gamma; ^ R ( k ) a ^ ( &theta; ^ ( k ) , r ^ ( k ) )
根据利用上述空间谱公式求取对应的P个最大值时的距离便可得到第k次迭代的P个目标距离的估计值。
8.2令w=[1,0,...,0]T,则满足约束条件则利用拉格朗日乘子法,可求得幅相误差矩阵如下:
&Gamma; ^ &theta; ( k + 1 ) = d i a g { &delta; ^ ( k + 1 ) } &delta; ^ ( k + 1 ) = ( Q k ) - 1 w w T ( Q k ) - 1 w
其中, Q k = &Sigma; p = 1 P a ~ H ( R ^ p ( k ) ) R - 1 a ~ ( R ^ p ( k ) ) .
8.3判断是否收敛
J k = ( &delta; ^ ( k + 1 ) ) H Q k &delta; ^ ( k + 1 )
在满足Jk-1-Jk>ε(ε为预设的门限值)时,继续迭代;否则,结束迭代。最后一次迭代得到的角度估计值和幅相误差估计值为
步骤9
由于阵列幅相误差来自阵列通道的不一致性,比较两个维度上估计出的幅相误差对角矩阵,若满足则停止迭代,即得到幅相误差的估计值否则,令k=k+1,转到步骤7,继续迭代。
本发明的效果可以通过下列仿真实验加以说明:
仿真条件
本例中的无源定位阵列采用步骤2中所述方法产生的稀布阵列,如图1所示,沿着x轴方向,在0~100m范围内稀疏布置8个阵元,近场目标窄带不相关信号源个数P=1,位置分别(-22.5°,1700m),波长λ=0.3m,信噪比均为SNR=20dB,快拍数L=400,角度维在-45°~45°范围内以0.1°为步长进行搜索,距离维在1200m~2000m范围内以1m为步长进行搜索。设相位均方根误差δp=10°,幅度均方根误差δa=1dB,即
&Delta; &phi; ~ U ( - 12 &sigma; &phi; 2 / 2 , 12 &sigma; &phi; 2 / 2 ) &sigma; &phi; = &delta; &phi; &pi; 180
&Delta; a ~ U ( - 12 &sigma; a 2 / 2 , 12 &sigma; a 2 / 2 ) &delta; a = 10 lg&sigma; a
计算可得相位误差△φ~U(-0.3023,0.3023),单位rad,或者△φ~U(-17.3°,17.3°)幅度误差△a~U(-2.1085,2.1085),即接收信号幅度为a,则加入幅度误差后的幅度a'=a+△a。
仿真结果
图3给出了存在幅相误差的情况下,近场信号无源定位的二维方向图,旁瓣峰值达到-5dB左右。采用本发明的方法进行幅相误差校正之后,获得的二维方向图如4所示,旁瓣峰值降低到-20dB左右,有非常明显的校正效果。图5和图6还分别给出了,校正误差前后在二维方向图的峰值位置上距离维和角度维的切面图对比,显然,经过本发明的幅相误差校正后,即使在存在幅相误差的情况下仍然有效保证了近场无源定位二维方向图的性能。

Claims (2)

1.一种稀布阵列近场无源定位幅相误差校正方法,该方法具体包括:
步骤1:获得定位场中信号接收机即定位阵元的个数及其位置信息;
步骤2:以第一个阵元为参考阵元,各阵元接收P个位于定位场中的目标信号源发送的信号;
步骤3:用xm(t)表示第m个阵元在第t个时刻接收信号的采样值,则天线阵列接收信号
x ( t ) = &lsqb; x 1 ( t ) , x 2 ( t ) , ... , x M ( t ) &rsqb; T = &Sigma; p = 1 P a ( r p , &theta; p ) s p ( t ) + n ( t ) , t = 1 , 2 , ... , L
其中,(·)T表示转置操作符,sp(t)表示第p个近场信号源,n(t)表示天线阵列接收的噪声矢量,L表示采样总数,a(rpp)表示第p个近场信号源导向矢量,rp表示第p个近场信号源到参考阵元的距离;
步骤4:通过接收天线阵列数据x(t)估算协方差矩阵(·)H表示对矩阵进行共轭转置计算;再估算协方差矩阵的逆矩阵R-1,(·)-1表示对矩阵进行求逆运算;
步骤5:采用最小方差波束形成算法估计出目标角度和距离的初始值 &theta; ^ ( 0 ) = ( &theta; ^ 1 ( 0 ) , &theta; ^ 2 ( 0 ) , . . . , &theta; ^ P ( 0 ) ) , 令二维迭代次数N=0;
步骤6:首先在角度维上进行迭代,固定目标距离为则可以在角度维上进行一维幅相误差校正,设定角度维上迭代次数变量为k,从角度维上估计出的稀布阵列雷达幅相误差对角矩阵为初始值为
定义代价函数
J k = &Sigma; p = 1 P a H ( r p ( k ) , &theta; p ( k ) ) ( &Gamma; &theta; ( k ) ) H R - 1 &Gamma; &theta; ( k ) a ( r p ( k ) , &theta; p ( k ) )
其中a(rp (k)p (k))表示第p个近场信号源的导向矢量,(·)H表示对矩阵进行共轭转置计算;增加迭代次数,直到Jk-1-Jk≤ε1结束迭代,获得角度估计值和幅相误差估计值为
步骤7:在距离维上进行迭代,固定目标角度为则可以在距离维上进行一维幅相误差校正;设定从距离维上估计出的稀布阵列雷达幅相误差对角矩阵初始值为进行与步骤7相同的一维迭代过程,可求得距离估计值和幅相误差估计值
步骤8:比较两个维度上估计出的幅相误差对角矩阵,若满足则停止迭代,其中ε2为根据实际情况预设的门限,即得到幅相误差的估计值否则,令N=N+1,转到步骤6,继续迭代,其中N为角度维和距离维共同迭代次数的变量。
2.如权利要求1所述的一种稀布阵列近场无源定位幅相误差校正方法,其特征在于所述步骤6的具体步骤为:
步骤6.1:定义空间谱公式为:
p ( &theta; ^ ( k ) | &Gamma; ^ &theta; ( k ) , r ^ ( k ) ) = 1 a ^ H ( &theta; ^ ( k ) , r ^ ( k ) ) ( &Gamma; ^ &theta; ( k ) ) H R - 1 &Gamma; ^ &theta; ( k ) a ^ ( &theta; ^ ( k ) , r ^ ( k ) )
根据利用上述空间谱公式求取对应的P个最大值时的角度便可得到第k次迭代的P个波达方向的估计值;
步骤6.2:令w=[1,0,...,0]T,则满足约束条件则用拉格朗日乘子法,可求得幅相误差矩阵如下:
&Gamma; ^ &theta; ( k + 1 ) = d i a g { &delta; ^ ( k + 1 ) } &delta; ^ ( k + 1 ) = ( Q k ) - 1 w w T ( Q k ) - 1 w
其中, Q k = &Sigma; p = 1 P a ~ H ( r ^ p ( k ) , &theta; ^ p ( k ) ) R - 1 a ~ ( r p ( k ) , &theta; ^ p ( k ) ) , a ~ ( r ^ p ( k ) , &theta; ^ p ( k ) ) = d i a g { a ( r ^ p ( k ) , &theta; ^ p ( k ) ) } ;
步骤6.3判断是否收敛
J k = ( &delta; ^ ( k + 1 ) ) H Q k &delta; ^ ( k + 1 )
在满足Jk-1-Jk1,时,继续迭代,ε1为根据实际情况预设的门限值;否则,结束迭代;最后一次迭代得到角度估计值和幅相误差估计值为
CN201510770276.2A 2015-11-12 2015-11-12 一种稀布阵列近场无源定位幅相误差校正方法 Active CN105445709B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510770276.2A CN105445709B (zh) 2015-11-12 2015-11-12 一种稀布阵列近场无源定位幅相误差校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510770276.2A CN105445709B (zh) 2015-11-12 2015-11-12 一种稀布阵列近场无源定位幅相误差校正方法

Publications (2)

Publication Number Publication Date
CN105445709A true CN105445709A (zh) 2016-03-30
CN105445709B CN105445709B (zh) 2017-11-17

Family

ID=55556128

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510770276.2A Active CN105445709B (zh) 2015-11-12 2015-11-12 一种稀布阵列近场无源定位幅相误差校正方法

Country Status (1)

Country Link
CN (1) CN105445709B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093898A (zh) * 2016-08-23 2016-11-09 中国电子科技集团公司第四十研究所 一种分区域式的mimo阵列校准方法
CN107783087A (zh) * 2017-10-27 2018-03-09 西南电子技术研究所(中国电子科技集团公司第十研究所) 球面相控阵天线近场通道标校链路的自校正方法
CN108107276A (zh) * 2017-10-27 2018-06-01 西南电子技术研究所(中国电子科技集团公司第十研究所) 球面阵天线通道相对幅相特性的近场标校方法
CN108828667A (zh) * 2018-03-23 2018-11-16 中国矿业大学(北京) 一种微地震复杂地表高程校正方法
CN109613474A (zh) * 2018-12-17 2019-04-12 东南大学 一种适用于短距离车载雷达的测角补偿方法
CN109752705A (zh) * 2017-11-03 2019-05-14 中电科海洋信息技术研究院有限公司 高频水声阵列性能参数测量方法及系统、设备及存储介质
CN110516286A (zh) * 2019-07-09 2019-11-29 西安电子科技大学 基于稀疏重构算法的阵列幅相误差估计方法
CN113362856A (zh) * 2021-06-21 2021-09-07 国网上海市电力公司 一种应用于电力物联网的声音故障检测方法以及装置
CN113625225A (zh) * 2021-08-12 2021-11-09 电子科技大学 一种基于多个均匀线阵的近场三维空间定位方法
CN115833972A (zh) * 2023-02-17 2023-03-21 华清瑞达(天津)科技有限公司 对阵列馈电系统输出的目标角位置进行校准的系统及方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644765A (zh) * 2009-04-23 2010-02-10 中国科学院声学研究所 一种用于水声换能器线阵的幅相误差校正方法
US20120268325A1 (en) * 2011-04-21 2012-10-25 Hansen Charles T Maximum likelihood angle estimation of wideband signals using phased array antennas
CN103353588A (zh) * 2013-06-13 2013-10-16 西安电子科技大学 基于天线均匀平面阵的二维波达方向角估计方法
CN103885048A (zh) * 2014-03-20 2014-06-25 西安电子科技大学 双基地mimo雷达收发阵列幅相误差的校正方法
CN104535987A (zh) * 2014-12-15 2015-04-22 南京航空航天大学 适用于均匀圆阵列声纳系统的幅相误差自校正方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101644765A (zh) * 2009-04-23 2010-02-10 中国科学院声学研究所 一种用于水声换能器线阵的幅相误差校正方法
US20120268325A1 (en) * 2011-04-21 2012-10-25 Hansen Charles T Maximum likelihood angle estimation of wideband signals using phased array antennas
CN103353588A (zh) * 2013-06-13 2013-10-16 西安电子科技大学 基于天线均匀平面阵的二维波达方向角估计方法
CN103885048A (zh) * 2014-03-20 2014-06-25 西安电子科技大学 双基地mimo雷达收发阵列幅相误差的校正方法
CN104535987A (zh) * 2014-12-15 2015-04-22 南京航空航天大学 适用于均匀圆阵列声纳系统的幅相误差自校正方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
鲁祖坤 等: ""基于子空间的幅相误差自校正算法"", 《现代导航》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106093898A (zh) * 2016-08-23 2016-11-09 中国电子科技集团公司第四十研究所 一种分区域式的mimo阵列校准方法
CN106093898B (zh) * 2016-08-23 2018-05-25 中国电子科技集团公司第四十一研究所 一种分区域式的mimo阵列校准方法
CN107783087A (zh) * 2017-10-27 2018-03-09 西南电子技术研究所(中国电子科技集团公司第十研究所) 球面相控阵天线近场通道标校链路的自校正方法
CN108107276A (zh) * 2017-10-27 2018-06-01 西南电子技术研究所(中国电子科技集团公司第十研究所) 球面阵天线通道相对幅相特性的近场标校方法
CN109752705A (zh) * 2017-11-03 2019-05-14 中电科海洋信息技术研究院有限公司 高频水声阵列性能参数测量方法及系统、设备及存储介质
CN108828667A (zh) * 2018-03-23 2018-11-16 中国矿业大学(北京) 一种微地震复杂地表高程校正方法
CN108828667B (zh) * 2018-03-23 2019-09-03 中国矿业大学(北京) 一种微地震复杂地表高程校正方法
CN109613474A (zh) * 2018-12-17 2019-04-12 东南大学 一种适用于短距离车载雷达的测角补偿方法
CN109613474B (zh) * 2018-12-17 2022-09-30 东南大学 一种适用于短距离车载雷达的测角补偿方法
CN110516286A (zh) * 2019-07-09 2019-11-29 西安电子科技大学 基于稀疏重构算法的阵列幅相误差估计方法
CN113362856A (zh) * 2021-06-21 2021-09-07 国网上海市电力公司 一种应用于电力物联网的声音故障检测方法以及装置
CN113625225A (zh) * 2021-08-12 2021-11-09 电子科技大学 一种基于多个均匀线阵的近场三维空间定位方法
CN115833972A (zh) * 2023-02-17 2023-03-21 华清瑞达(天津)科技有限公司 对阵列馈电系统输出的目标角位置进行校准的系统及方法
CN115833972B (zh) * 2023-02-17 2023-04-28 华清瑞达(天津)科技有限公司 对阵列馈电系统输出的目标角位置进行校准的系统及方法

Also Published As

Publication number Publication date
CN105445709B (zh) 2017-11-17

Similar Documents

Publication Publication Date Title
CN105445709B (zh) 一种稀布阵列近场无源定位幅相误差校正方法
CN110113085B (zh) 一种基于协方差矩阵重构的波束形成方法及系统
CN104977558B (zh) 一种基于贝叶斯压缩感知的分布源中心波达方向估计方法
CN108459307B (zh) 基于杂波的mimo雷达收发阵列幅相误差校正方法
CN104407335B (zh) 一种3轴交叉阵列的doa估计方法
CN107390197B (zh) 基于特征空间的雷达自适应和差波束测角方法
CN107092007A (zh) 一种虚拟二阶阵列扩展的波达方向估计方法
CN111487599B (zh) 增益-相位误差背景下双基地mimo雷达角度估计方法及装置
CN113835063B (zh) 一种无人机阵列幅相误差与信号doa联合估计方法
CN109212466B (zh) 一种基于量子蜻蜓演化机制的宽带测向方法
CN107315162A (zh) 基于内插变换和波束形成的远场相干信号doa估计方法
CN112379327A (zh) 一种基于秩损估计的二维doa估计与互耦校正方法
CN113189592B (zh) 考虑幅相互耦误差的车载毫米波mimo雷达测角方法
CN111046591A (zh) 传感器幅相误差与目标到达角度的联合估计方法
CN110138430A (zh) 期望信号doa误差下基于导向矢量估计的稳健宽带波束形成方法
CN111352063A (zh) 一种均匀面阵中基于多项式求根的二维测向估计方法
CN107966676B (zh) 复杂噪声环境下阵列天线角度和信源个数联合估计方法
CN113671439A (zh) 基于非均匀智能超表面阵列的无人机集群测向系统及方法
CN112255629A (zh) 基于联合uca阵列的序贯esprit二维不相干分布源参数估计方法
CN106125039B (zh) 基于局域联合处理的改进空时自适应单脉冲测角方法
CN107167803A (zh) 基于导向矢量失配估计的鲁棒波束域自适应波束形成方法
Abdullah et al. Comparative Study of Super-Performance DOA Algorithms based for RF Source Direction Finding and Tracking
CN115980721A (zh) 一种无误差协方差矩阵分离的阵列自校正方法
Mills et al. Fast iterative interpolated beamforming for interference DOA estimation in GNSS receivers using fully augmentable arrays
CN112881971A (zh) 电磁方向性互耦效应下相干干扰源测向方法

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