CN102736057B - 一种基于距离和跟踪再生载波相位的二维角度解算方法 - Google Patents
一种基于距离和跟踪再生载波相位的二维角度解算方法 Download PDFInfo
- Publication number
- CN102736057B CN102736057B CN201210200067.0A CN201210200067A CN102736057B CN 102736057 B CN102736057 B CN 102736057B CN 201210200067 A CN201210200067 A CN 201210200067A CN 102736057 B CN102736057 B CN 102736057B
- Authority
- CN
- China
- Prior art keywords
- angle
- amin
- pitch
- value
- target
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于距离和跟踪再生载波相位的二维角度解算方法,包括如下步骤(1)将天线阵阵元的三维空间坐标、天线阵阵元的接收信号跟踪再生载波相位以及天线阵到目标的距离作为已知参数,将目标的方位角与俯仰角作为自变量,构造一个测角目标函数;(2)将二维测角问题转换为测角目标函数的最大值搜索问题,通过搜索最大值得到方位角与俯仰角度值;本发明通过搜索最优值即可得到方位角与俯仰角度值,实现了近距离目标的方位角与俯仰角的高精度测量,且有效抑制空间角度估计的模糊性。
Description
技术领域
本发明涉及一种基于距离和跟踪再生载波相位的二维角度解算方法,属于无线电测量技术领域。
背景技术
在公开刊物及公开渠道上了解到的方位、俯仰二维角度无线电测量中使用的方位角、俯仰角解算方法有2种:
(1)干涉法
干涉法的基本原理如图1所示:平面波以一定入射角到达2个天线阵元时,平面波必然先到达其中的一个阵元,后到达另一个阵元,到达时间不同导致两天线阵元的接收信号存在相位差,测量出相位差后,利用相位差与入射角之间的关系(式中,d为二天线阵元间距,θ为来波方向与天线法线方向间的夹角,为二天线接收目标信号之间的相位差,λ是入射波长)可以计算出入射角。
此方法假定入射波是平面波,应用于远距离目标测角且信号质量较好是具有较高的精度,但是应用于近距离目标测角时算法本身会引入较大的系统误差。
(2)空间谱估计法
空间谱估计法是一系列方法的总称,它的基本思想是利用协方差矩阵的特征结构,将任意阵列输出数据的协方差矩阵进行特征分解,从而得到与信号分量相对应的信号子空间和与信号分量相正交的噪声子空间,最后利用信号方向向量与噪声子空间的正交性来构成空间扫描谱,通过谱峰搜索,实现信号的入射方向估计。
此类方法假定入射波是平面波,在远距离目标测角时具有较高的精度和分辨率,应用于近距离目标测角时算法本身会引入较大的系统误差。同时此类方法计算量非常大,不适于处理器和存储器资源严重受限应用场景,例如深空探测中的飞行器间相对角度测量。
发明内容
本发明的目的在于克服现有技术的上述不足,提供一种基于距离和跟踪再生载波相位的二维角度解算方法,该方法使用测角天线阵多个阵元接收信号的跟踪再生载波相位和天线阵面几何中心与目标飞行器间的距离测量值作为输入,将二维测角问题转换为一个测角目标函数的最大值搜索问题,通过搜索最优值即可得到方位角与俯仰角度值,实现了近距离目标的方位角与俯仰角的高精度测量,且有效抑制空间角度估计的模糊性。
本发明的上述目的主要是通过如下技术方案予以实现的:
一种基于距离和跟踪再生载波相位的二维角度解算方法,包括如下步骤:
步骤(二)、计算复向量v,公式如下:
其中:py1、py2、p00、pz1、pz2分别为被测目标向测角天线阵发射单频信号时,阵元y1、y2、00、z1、z2接收的信号跟踪再生载波相位;j表示复数虚部;
步骤(三)、令方位角α=aminα、俯仰角β=aminβ,方位角α、俯仰角β的测角范围分别为[aminα,amaxα]、[aminβ,amaxβ],通过如下公式计算目标P到阵元y1、y2、00、z1、z2的距离Ry1、Ry2、R00、Rz1、Rz2:
其中:R为目标P到天线阵的距离;
(xy1,yy1,zy1)、(xy2,yy2,zy2)、(x00,y00,z00)、(xz1,yz1,zz1)、(xz2,yz2,zz2)分别为阵元y1、y2、00、z1、z2的坐标;
步骤(四)、计算载波相位差复向量函数s(α,β),公式如下:
j表示复数虚部;
步骤(五)、计算测角目标函数m(α,β),将计算出的m(α,β)记为m0×0,公式如下:
其中:H表示复向量的共轭转置操作;·表示复向量内积:
步骤(六)、令方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m,设立如下三个条件:(1)n、m是正整数;(2)aminα+stα×n≤amaxα;(3)aminβ+stβ×m≤amaxβ;重复步骤(三)~步骤(五),计算所有满足条件(1)、(2)、(3)的方位角α=aminα+stα×n、俯仰角β=aminβ+tβ×m组合上的测角目标函数值m(α,β),并将方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m组合上的测角目标函数值m(α,β)记为mn×m;
其中:stα为方位角的步长,stβ为俯仰角的步长;
步骤(七)、搜索步骤(五)计算出的m0×0和步骤(六)计算出的mn×m中的最大值,具体方法如下:
(1)设i、j、a、b、result为变量,且i、j、a、b、result的初始值置为0;
(2)如果mi×j大于等于result,
将mi×j的值赋给result,将i的值赋给a,将j的值赋给b;
否则,
result、a、b取值不变;
(3)如果aminα+stα×(i+1)≤amaxα;
i取值加1,执行第(2)步;
否则,
如果aminβ+stβ×(j+1)≤amaxβ:
i置为0,j取值加1,执行第(2)步;
否则,
执行第(4)步;
(4)输出a,b取值;
步骤(八)、计算目标方位角α和俯仰角β,目标方位角α=aminα+stα×a,俯仰角β=aminβ+stβ×b。
本发明与现有公开方法相比的优点在于:
(1)本发明使用测角天线阵多个阵元接收信号的跟踪再生载波相位和天线阵面几何中心与目标飞行器间的距离测量值作为输入,将二维测角问题转换为一个测角目标函数的最大值搜索问题,通过搜索最优值即可得到方位角与俯仰角度值,实现了近距离目标的方位角与俯仰角的高精度测量;
(2)本发明采用信号传输的球面波前特性,准确使用距离测量量,可以用于近距离目标的方位角与俯仰角的高精度测量;
(3)本发明使用相位信息作为复数的辐角构造测角目标函数,不会出现相位整周模糊问题;
(4)本发明同时运用面阵阵元接收信号的相位差和目标与面阵间的距离进行解算,能够有效抑制空间角度估计的模糊性,支持大角度范围的近距离高精度测角。
附图说明
图1为干涉法的基本原理示意图;
图2为本发明测量坐标系定义及阵元位置示意图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
测量坐标系定义及天线阵元位置如图2所示:OXYZ构成直角坐标系,符合右手系;P点为目标位置,R表示目标到坐标原点的距离;α和β分别为目标P在OXYZ直角坐标系中的方位角和俯仰角,即本发明需要解算输出的二维角度,在如图1的坐标系定义中:当α=0时,β定义为绕OY轴旋转;当β=0时,α定义为绕OZ轴旋转(图1所示的α<0、β<0);y1、y2、00、z1、z2表示5个测角天线阵元,位于YOZ平面,构成L型阵列,其中y1、y2、00在一条直线上,z1、z2、00在一条直线上。
当目标P点向天线阵元发射单频信号时,阵元00、y1、y2、z1、z2接收信号跟踪再生载波相位依次为p00、py1、py2、pz1、pz2。
使用上述5个跟踪再生载波相位构造复数向量v,j表示复数虚部:
记y1、y2、00、z1、z2阵元坐标记为(xy1,yy1,zy1)、(xy2,yy2,zy2)、(x00,y00,z00)、(xz1,yz1,zz1)、(xz2,yz2,zz2)。
记P点坐标为(xp,yp,zp),记目标P到阵元y1、y2、00、z1、z2的距离依次为Ry1、Ry2、R00、Rz1、Rz2。Ry1、Ry2、R00、Rz1、Rz2与α,β、R、(xy1,yy1,zy1)、(x2,yy2,zy2)、(x00,y00,z00)、(xz1,yz1,zz1)、(xz2,yz2,zz2)满足(2)式的关系。其中α,β是未知量,其余都是已知量。
在Ry1、Ry2、R00、Rz1、Rz2的基础上定义载波相位差复向量函数s(α,β)。其中Ry1、Ry2、R00、Rz1、Rz2通过(2)式计算,j表示复数虚部,k表示单频信号的波数(等于λ为波长),
在复数向量v和载波相位差复向量函数s(α,β)的基础上定义测角目标函数m(α,β),其中H表示复向量的共轭转置操作,·表示复向量内积,||表示复数求模操作。
记方位角α、俯仰角β测角范围分别是[aminα,amaxα]、[aminβ,amaxβ],其中aminα,amaxα是方位角α测角范围的下、上边界,aminβ,amaxβ是俯仰角β测角范围的下、上边界。
记方位角α、俯仰角β测角精度分别是δα、δβ。
实施例1
已知如下条件:
1)、被测目标向测角天线阵发射单频信号时,阵元y1、y2、00、z1、z2接收信号跟踪再生载波相位p00=2.0847、py1=1.1528、py2=7.1506e-005、pz1=-1.2997、pz2=-2.1579,跟踪再生载波相位单位是弧度;
2)、00、y1、y2、z1、z2阵元坐标依次是(x00,y00,z00)=(0,0.1249,-0.1249)、(xy1,yy1,zy1)=(0,0.0096,-0.1249)、(xy2,yy2,zy2)=(0,-0.1249,-0.1249)、(xz1,yz1,zz1)=(0,0.1249,-0.0096)、(xz2,yz2,zz2)=(0,0.1249,0.1249),坐标分量单位为米;
3)、目标P到天线阵距离R=5,距离的单位为米;
4)、方位角、俯仰角的测角范围[aminα,amaxα]=[-0.1745,0.1745]、[aminβ,amaxβ]=[-0.1745,0.1745],测角范围上下界的单位是弧度;
5)、方位角、俯仰角测角精度分别是δα=0.00175、δβ=0.00175,测角精度的单位是弧度;
6)、目标点发射到天线阵元的信号载波波长入射波λ=0.0096,波长的单位是米。
本实施例解算目标方位角α弧度,俯仰角β弧度的具体实施步骤如下:
步骤(二):计算复向量v。复向量v使用公式(1)计算,其中阵元y1、y2、00、z1、z2接收信号跟踪再生载波相位p00、py1、py2、pz1、pz2是已知量,j表示复数虚部:
步骤(三):令方位角α=aminα=-0.1745、俯仰角β=aminβ=-0.1745,计算目标P到阵元y1、y2、00、z1、z2的距离Ry1、Ry2、R00、Rz1、Rz2,这5个距离的单位是米。距离Ry1、Ry2、R00、Rz1、Rz2采用公式(2)进行计算,其中y1、y2、00、z1、z2阵元坐标(xy1,yy1,zy1)、(xy2,yy2,zy2)、(x00,y00,z00)、(xz1,yz1,zz1)、(xz2,yz2,zz2)是已知量:
步骤(四):计算载波相位差复向量函数s(α,β)。s(α,β)使用公式(3)计算,其中Ry1、Ry2、R00、Rz1、Rz2是步骤(三)计算的结果,j表示复数虚部,k表示单频信号的波数(等于π为圆周率,λ为波长):
步骤(五):计算测角目标函数m(α,β),将计算出的m(α,β)记为m0×0,m(α,β)使用公式(4)计算,其中v已由步骤(二)计算出,s(α,β)已由步骤(四)计算出,H表示复向量的共轭转置操作,·表示复向量内积,||表示复数求模操作:
步骤(六):令方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m,其中1)n、m是正整数;2)aminα+stα×n≤amaxα;3)aminβ+stβ×m≤amaxβ。计算所有满足条件1)、2)、3)的方位角、俯仰角组合(α=aminα+stα×n,β=aminβ+stβ×m)上的测角目标函数值m(α,β)。
每一个方位角、俯仰角组合(α,β)上的测角目标函数m(α,β)的计算过程与步骤(三)、步骤(四)和步骤(五)相同:首先采用步骤(三)中的公式(2)计算目标P到阵元y1、y2、00、z1、z2的距离Ry1、Ry2、R00、Rz1、Rz2,然后使用步骤(四)中的公式(3)计算载波相位差复向量函数s(α,β),最后使用步骤(五)中的公式(4)计算出m(α,β)。将方位角、俯仰角组合(α=aminα+stα×n,β=aminβ+stβ×m)上的测角目标函数值m(α,β)记为mn×m。
步骤(七):搜索步骤(五)计算出的m0×0和步骤(六)计算出的mn×m(n、m是正整数;aminα+stα×n≤amaxα;aminβ+stβ×m≤amaxβ)中的最大值。设j、j、a、b、result为变量,搜索方法如下:
第(1)步:i、j、a、b、result的初始值置为0;
第(2)步:
如果mi×j大于等于result:
将mi×j的值赋给result,将i的值赋给a,将j的值赋给b。
否则,
result、a、b取值不变。
第(3)步:
如果aminα+stα×(i+1)≤amaxα:
i取值加1,执行第(2)步。
否则,
如果aminβ+stβ×(j+1)≤amaxβ:
i置为0,j取值加1,执行第(2)步。
否则,
继续执行第(4)步。
第(4)步:输出a=19,b=182,搜索出最大值。
步骤(八):计算目标方位角α和俯仰角β:方位角α=aminα+stα×a=-0.1431,俯仰角β=aminβ+stβ×b=0.1414,方位角、俯仰角的单位是弧度。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。
Claims (1)
1.一种基于距离和跟踪再生载波相位的二维角度解算方法,其特征在于:包括如下步骤:
步骤(二)、计算复向量v,公式如下:
其中:py1、py2、p00、pz1、pz2分别为被测目标向测角天线阵发射单频信号时,阵元y1、y2、00、z1、z2接收的信号跟踪再生载波相位;j表示复数虚部;
步骤(三)、令方位角α=aminα、俯仰角β=aminβ,方位角α、俯仰角β的测角范围分别为[aminα,amaxα]、[aminβ,amaxβ],通过如下公式计算目标P到阵元y1、y2、00、z1、z2的距离Ry1、Ry2、R00、Rz1、Rz2:
其中:R为目标P到天线阵的距离;
(xy1,yy1,zy1)、(xy2,yy2,zy2)、(x00,y00,z00)、(xz1,yz1,zz1)、(xz2,yz2,zz2)分别为阵元y1、y2、00、z1、z2的坐标;
步骤(四)、计算载波相位差复向量函数s(α,β),公式如下:
j表示复数虚部;
步骤(五)、计算测角目标函数m(α,β),将计算出的m(α,β)记为m0×0,公式如下:
其中:H表示复向量的共轭转置操作;·表示复向量内积:
步骤(六)、令方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m,设立如下三个条件:(1)n、m是正整数;(2)aminα+stα×n≤amaxα;(3)aminβ+stβ×m≤amaxβ;重复步骤(三)~步骤(五),计算所有满足条件(1)、(2)、(3)的方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m组合上的测角目标函数值m(α,β),并将方位角α=aminα+stα×n、俯仰角β=aminβ+stβ×m组合上的测角目标函数值m(α,β)记为mn×m;
其中:stα为方位角的步长,stβ为俯仰角的步长;
步骤(七)、搜索步骤(五)计算出的m0×0和步骤(六)计算出的mn×m中的最大值,具体方法如下:
(1)设i、k’、a、b、result为变量,且i、k’、a、b、result的初始值置为0;
(2)如果mi×k′大于或等于result,
将mi×k′的值赋给result,将i的值赋给a,将k’的值赋给b;
否则,
result、a、b取值不变;
(3)如果aminα+stα×(i+1)≤amaxα;
i取值加1,执行第(2)步;
否则,
如果aminβ+stβ×(k′+1)≤amaxβ:
i置为0,k’取值加1,执行第(2)步;
否则,
执行第(4)步;
(4)输出a,b取值;
步骤(八)、计算目标方位角α和俯仰角β,目标方位角α=aminα+stα×a,俯仰角β=aminβ+stβ×b。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210200067.0A CN102736057B (zh) | 2012-06-18 | 2012-06-18 | 一种基于距离和跟踪再生载波相位的二维角度解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210200067.0A CN102736057B (zh) | 2012-06-18 | 2012-06-18 | 一种基于距离和跟踪再生载波相位的二维角度解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102736057A CN102736057A (zh) | 2012-10-17 |
CN102736057B true CN102736057B (zh) | 2014-05-28 |
Family
ID=46991891
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210200067.0A Active CN102736057B (zh) | 2012-06-18 | 2012-06-18 | 一种基于距离和跟踪再生载波相位的二维角度解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102736057B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108344969B (zh) * | 2018-01-08 | 2023-05-09 | 南京航空航天大学 | 一种球面波干涉测向方法 |
CN110873824A (zh) * | 2018-08-31 | 2020-03-10 | 南京捷希科技有限公司 | 一种Massive MIMO天线测试系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5457466A (en) * | 1993-09-07 | 1995-10-10 | Litton Systems Inc. | Emitter azimuth and elevation direction finding using only linear interferometer arrays |
CN101109799A (zh) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | 立体基线宽频带微波辐射源测向装置及测向方法 |
CN101963659A (zh) * | 2010-10-25 | 2011-02-02 | 哈尔滨工程大学 | 利用三维正交干涉仪实现宽频段测向的方法 |
CN102175988A (zh) * | 2011-01-21 | 2011-09-07 | 电子科技大学 | 一种基于维度拆分的相关干涉仪测向方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6784840B2 (en) * | 2002-12-23 | 2004-08-31 | Itt Manufacturing Enterprises, Inc. | Method for determining azimuth and elevation angles using a single axis direction finding system |
-
2012
- 2012-06-18 CN CN201210200067.0A patent/CN102736057B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5457466A (en) * | 1993-09-07 | 1995-10-10 | Litton Systems Inc. | Emitter azimuth and elevation direction finding using only linear interferometer arrays |
CN101109799A (zh) * | 2007-07-06 | 2008-01-23 | 哈尔滨工程大学 | 立体基线宽频带微波辐射源测向装置及测向方法 |
CN101963659A (zh) * | 2010-10-25 | 2011-02-02 | 哈尔滨工程大学 | 利用三维正交干涉仪实现宽频段测向的方法 |
CN102175988A (zh) * | 2011-01-21 | 2011-09-07 | 电子科技大学 | 一种基于维度拆分的相关干涉仪测向方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102736057A (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102411136B (zh) | 一种扩展基线解模糊的相位干涉仪测向方法 | |
CN105676171B (zh) | 单通道双基站超短波信号空间定位方法 | |
CN107015191B (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN103018730B (zh) | 分布式子阵波达方向估计方法 | |
CN102540138B (zh) | 一种多基线相位搜索式二维空间谱测向方法 | |
CN102735207B (zh) | 一种基于距离和载波相差的二维角度解算方法 | |
CN103267966B (zh) | 一种单个二维压差式矢量水听器方位估计的方法 | |
CN106802403B (zh) | 声矢量传感器二维阵列music解相干参数估计方法 | |
CN103576138B (zh) | 一种基于gnss-r信号几何关系的星载无源雷达定位方法 | |
WO2007016098A2 (en) | Methods, apparatuses and systems for locating non-cooperative objects | |
CN103823217A (zh) | 基于双频发射的双基地mimo雷达高速运动目标参数估计方法 | |
CN103969629A (zh) | 一种基于主瓣杂波配准的机载雷达杂波自适应抑制方法 | |
CN103605113B (zh) | 多发多收干涉合成孔径雷达空时二维信号波形设计方法 | |
CN103604944A (zh) | 一种基于单站船载高频地波雷达的表面流测量方法 | |
CN102749612A (zh) | 一种相关干涉仪的快速测向方法 | |
CN103760546A (zh) | 一种雷达用低空目标波达方向估计方法 | |
CN104793177B (zh) | 基于最小二乘法的麦克风阵列测向方法 | |
CN102478652B (zh) | 基于移动平台的多普勒频差测向法 | |
CN106772326A (zh) | 一种多子阵合成孔径声纳相位误差分析方法 | |
CN103293526A (zh) | 基于合成导向矢量的波束形成米波雷达测高方法 | |
CN103869298A (zh) | 一种分布式mimo天波超视距雷达海杂波仿真方法 | |
CN105929377A (zh) | 一种基于单极子交叉环天线的高频雷达船舶方位角估计方法 | |
CN104407318A (zh) | 一种48元均匀圆柱阵超短基线精确定向方法 | |
CN102736057B (zh) | 一种基于距离和跟踪再生载波相位的二维角度解算方法 | |
CN106872936A (zh) | 近场源l型声矢量传感器阵列解模糊多参数估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |