CN108828508B - 一种超视距辐射源直接定位模型的性能分析方法 - Google Patents
一种超视距辐射源直接定位模型的性能分析方法 Download PDFInfo
- Publication number
- CN108828508B CN108828508B CN201810631219.XA CN201810631219A CN108828508B CN 108828508 B CN108828508 B CN 108828508B CN 201810631219 A CN201810631219 A CN 201810631219A CN 108828508 B CN108828508 B CN 108828508B
- Authority
- CN
- China
- Prior art keywords
- radiation source
- over
- base station
- horizon
- subscript
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/0009—Transmission of position information to remote stations
- G01S5/0018—Transmission from mobile station to base station
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提出一种超视距辐射源直接定位模型的性能分析方法,其方法包括:超视距辐射源定位系统的数学建模;入射角参量的费舍尔矩阵计算;辐射源位置参量的克拉美罗界计算。本发明所述方法在传统的视距辐射源定位性能分析中引入了超视距模型,给出了超视距条件下的辐射源定位的理论误差下限,突破了传统分析方法的局限性,能够有效地指导超视距定位系统设计。
Description
技术领域
本发明属于辐射源定位技术、雷达信号处理、阵列信号处理以及参数估计理论技术领域,特别是涉及一种超视距辐射源直接定位模型的性能分析方法。
背景技术
辐射源定位技术在雷达、声呐和无线通信领域都是一个重要的研究课题。在辐射源定位领域,多数研究将重点放在了辐射源定位方法的设计上,力求获得更有效的定位结果。而值得注意的是,同具体的定位方法一样,定位模型的理论性能分析也是辐射源定位技术的重要组成部分。尤其是在系统设计、阵型优化等工作中,模型的理论性能分析更是能起到至关重要的指导作用。
现有的辐射源直接定位模型的性能分析方法均只适用于视距范围内的目标,在实际应用中存在局限性。本发明中给出了一种适用于超视距辐射源直接定位模型的性能分析方法,该方法从克拉美罗界的角度给出了超视距定位系统的估计误差下限,能够对超视距定位系统设计提供有效的理论支撑。
发明内容
本发明的目的在于改善现有技术的局限性,提供了一种超视距辐射源直接定位模型的性能分析方法,该方法能够有效地给出超视距辐射源直接定位系统的理论误差下限。
本发明的目的是通过以下技术方案实现:一种超视距辐射源直接定位模型的性能分析方法,包括如下步骤:
步骤1、超视距辐射源定位系统的数学建模;
步骤2、入射角参量的费舍尔矩阵计算;
步骤3、辐射源位置参量的克拉美罗界计算。
进一步地,所述步骤1具体为:
假设地球为一半径为r的理想球体,球面上存在单一目标辐射源;定位系统共配置有L个基站,每个基站包含一个具有M个天线单元的均匀线阵;令代表球坐标系下辐射源的位置,其中上标s代表球坐标系,下标t代表目标,θt代表目标的俯仰角,代表目标的方位角;nl(k)代表功率为的加性高斯白噪声矢量,其中k代表第k个快拍,下标n代表噪声;sl(k)代表入射信号,则第l个基站的接收数据表示为:
其中
表示阵列的导向矢量,d为阵元间距,λ为入射信号波长;对于各基站阵列,规定指向北极点方向为0°方向,阵面相对于0°方向的顺时针偏角由αl表示;令表示球坐标系下各基站位置,其中上标s代表球坐标系,下标l代表第l个基站,θl代表基站l的俯仰角,代表基站l的方位角;入射信号与阵面夹角记为ψl,其与的关系可由下式给出:
其中
u′0,l=PUu0,l u′t,l=PUut,l
PU=U(UHU)-1UH
为北极点坐标;和分别为直角坐标系下辐射源位置和基站位置坐标,其中上标c代表直角坐标系;和为过的球面切平面xxl+yyl+zzl=r2上不与共线的任意两点。
进一步地,所述步骤2具体为:
入射角参量ψl的费舍尔矩阵具有如下形式:
F=diag[f1(ψ1) f2(ψ2) … fL(ψL)]
其中
下标s和下标n用以区分信号和噪声,K代表快拍总数,al(ψl)代表以入射角参量为自变量的阵列导向矢量。
进一步地,所述步骤3具体为:
入射角参量ψl的费舍尔矩阵F与辐射源位置参量的克拉美罗界矩阵C有如下关系:
C=(GFGH)-1
其中G代表角度参量与位置参量的转换矩阵,其表达式为:
取C的对角线元素即得到辐射源位置参量的估计误差下限。
本发明所述方法在传统的视距辐射源定位性能分析中引入了超视距模型,给出了超视距条件下的辐射源定位的理论误差下限,突破了传统分析方法的局限性,能够有效地指导超视距定位系统设计。
附图说明
图1为本发明所述方法流程图;
图2为超视距传播条件下的辐射源定位系统示意图;
图3为超视距辐射源定位理论误差下限与最大似然定位算法均方误差曲线的对比图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
结合图1,本发明提出一种超视距辐射源直接定位模型的性能分析方法,包括如下步骤:
步骤1、超视距辐射源定位系统的数学建模;
步骤2、入射角参量的费舍尔矩阵计算;
步骤3、辐射源位置参量的克拉美罗界计算。
结合图2,所述步骤1具体为:
为方便表述与理解,现将所用符号统一说明如下:矩阵与矢量以粗斜体符号表示;上标(·)T、(·)H和(·)-1分别表示转置、共轭转置和取逆操作符;符号||·||,和diag(·)分别表示2范数和对角化操作;I代表单位阵。
考虑地球为一半径为r的理想球体,球面上存在单一目标辐射源。定位系统共配置有L个基站,每个基站包含一个具有M个天线单元的均匀线阵。令代表球坐标系下辐射源的位置,其中上标s代表球坐标系,下标t代表目标,θt代表目标的俯仰角,代表目标的方位角;nl(k)代表功率为的加性高斯白噪声矢量,其中k代表第k个快拍,下标n代表噪声;sl(k)代表入射信号,则第l个基站的接收数据可表示为:
其中
表示阵列的导向矢量,d为阵元间距,λ为入射信号波长。对于各基站阵列,规定指向北极点方向为0°方向,阵面相对于0°方向的顺时针偏角由αl表示,令表示球坐标系下各基站位置,其中上标s代表球坐标系,下标l代表第l个基站,θl代表基站l的俯仰角,代表基站l的方位角。入射信号与阵面夹角记为ψl,其与的关系可由下式给出:
其中
u′0,l=PUu0,l u′t,l=PUut,l
PU=U(UHU)-1UH
为北极点坐标;和分别为直角坐标系下辐射源位置和基站位置坐标,其中上标c代表直角坐标系;和为过的球面切平面xxl+yyl+zzl=r2上不与共线的任意两点。
所述步骤2具体为:
入射角参量ψl的费舍尔矩阵具有如下形式:
F=diag[f1(ψ1) f2(ψ2) … fL(ψL)]
其中
下标s和下标n用以区分信号和噪声,K代表快拍总数,al(ψl)代表以入射角参量为自变量的阵列导向矢量。
所述步骤3具体为:
入射角参量ψl的费舍尔矩阵F与辐射源位置参量的克拉美罗界矩阵C有如下关系:
C=(GFGH)-1
其中G代表角度参量与位置参量的转换矩阵,其表达式为:
取C的对角线元素即得到辐射源位置参量的估计误差下限。
图3给出了超视距辐射源定位理论误差下限与最大似然定位算法均方误差曲线的对比,其中ML代表最大似然算法,CRB代表克拉美罗界,SNR代表信噪比。可以看出,按本发明所述方法推得的理论误差下限与定位算法的性能曲线贴合得很好,因此所述方法的有效性得到了验证。
以上对本发明所提供的一种超视距辐射源直接定位模型的性能分析方法,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (3)
1.一种超视距辐射源直接定位模型的性能分析方法,其特征在于,包括如下步骤:
步骤1、超视距辐射源定位系统的数学建模;
步骤2、入射角参量的费舍尔矩阵计算;
步骤3、辐射源位置参量的克拉美罗界计算;
所述步骤1具体为:
假设地球为一半径为r的理想球体,球面上存在单一目标辐射源;定位系统共配置有L个基站,每个基站包含一个具有M个天线单元的均匀线阵;令代表球坐标系下辐射源的位置,其中上标s代表球坐标系,下标t代表目标,θt代表目标的俯仰角,代表目标的方位角;nl(k)代表功率为的加性高斯白噪声矢量,其中k代表第k个快拍,下标n代表噪声;sl(k)代表入射信号,则第l个基站的接收数据表示为:
其中
表示阵列的导向矢量,d为阵元间距,λ为入射信号波长;对于各基站阵列,规定指向北极点方向为0°方向,阵面相对于0°方向的顺时针偏角由αl表示;令表示球坐标系下各基站位置,其中上标s代表球坐标系,下标l代表第l个基站,θl代表基站l的俯仰角,代表基站l的方位角;入射信号与阵面夹角记为ψl,其与的关系由下式给出:
其中
u′0,l=PUu0,l u′t,l=PUut,l
PU=U(UHU)-1UH
为北极点坐标;和分别为直角坐标系下辐射源位置和基站位置坐标,其中上标c代表直角坐标系;和为过的球面切平面xxl+yyl+zzl=r2上不与共线的任意两点。
2.根据权利要求1所述的一种超视距辐射源直接定位模型的性能分析方法,其特征在于,所述步骤2具体为:
入射信号与阵面夹角ψl的费舍尔矩阵具有如下形式:
F=diag[f1(ψ1) f2(ψ2) … fL(ψL)]
其中
下标s和下标n用以区分信号和噪声,表示信号功率,I表示单位矩阵,K代表快拍总数,al(ψl)代表以入射角参量为自变量的阵列导向矢量。
3.根据权利要求2所述的一种超视距辐射源直接定位模型的性能分析方法,其特征在于,所述步骤3具体为:
入射信号与阵面夹角ψl的费舍尔矩阵F与辐射源位置参量的克拉美罗界矩阵C有如下关系:
C=(GFGH)-1
其中G代表角度参量与位置参量的转换矩阵,其表达式为:
取C的对角线元素即得到辐射源位置参量的估计误差下限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810631219.XA CN108828508B (zh) | 2018-06-19 | 2018-06-19 | 一种超视距辐射源直接定位模型的性能分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810631219.XA CN108828508B (zh) | 2018-06-19 | 2018-06-19 | 一种超视距辐射源直接定位模型的性能分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108828508A CN108828508A (zh) | 2018-11-16 |
CN108828508B true CN108828508B (zh) | 2019-06-21 |
Family
ID=64142674
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810631219.XA Active CN108828508B (zh) | 2018-06-19 | 2018-06-19 | 一种超视距辐射源直接定位模型的性能分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108828508B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117630993B (zh) * | 2024-01-15 | 2024-04-12 | 华中科技大学 | 一种基于sair多快拍的rfi源地理定位方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE69936354T2 (de) * | 1999-01-08 | 2008-06-19 | Trueposition, Inc. | Synthese der bandbreite für ein drahtloses positionsbestimmungssystem |
US7342536B2 (en) * | 2005-09-12 | 2008-03-11 | Lockheed Martin Corporation | System and method for determining the location of emitters having well-behaved timing features |
CN104375121A (zh) * | 2014-01-28 | 2015-02-25 | 大连大学 | 基于先验信息的mimo雷达波形与有偏估计器的联合优化方法 |
CN104808179A (zh) * | 2015-04-09 | 2015-07-29 | 大连大学 | 杂波环境下基于克拉美罗界的mimo雷达波形优化方法 |
CN105445718B (zh) * | 2015-11-19 | 2018-01-26 | 哈尔滨工业大学 | 一种基于阵列重构的分布式多载舰超视距雷达的doa估计方法 |
CN105676217B (zh) * | 2016-03-29 | 2017-11-14 | 电子科技大学 | 一种改进的ml天波雷达机动目标参数估计方法 |
-
2018
- 2018-06-19 CN CN201810631219.XA patent/CN108828508B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108828508A (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10571544B2 (en) | Direction finding using signal power | |
CN107015191B (zh) | 一种在多径干扰环境下单偶极子极化敏感阵列降维doa估计方法 | |
CN106526530B (zh) | 基于传播算子的2-l型阵列二维doa估计算法 | |
CN1664611A (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN105204008B (zh) | 一种基于协方差矩阵扩展的自适应天线波束形成零陷展宽方法 | |
CN105137430B (zh) | 一种前视阵列sar的回波稀疏获取及其三维成像方法 | |
CN103728601B (zh) | 雷达信号运动干扰空域-极化域联合稳健滤波方法 | |
US8633850B2 (en) | Identifying a location of a target object using a monopulse radar system and space-time adaptive processing (STAP) | |
CN103837867A (zh) | 一种利用ais信息进行高频雷达天线通道校正的方法 | |
CN106526531A (zh) | 基于三维天线阵列的改进传播算子二维doa估计算法 | |
CN107124216A (zh) | 一种针对阵列误差的Capon稳健自适应波束形成方法及系统 | |
CN106872934B (zh) | L型电磁矢量传感器阵列解相干esprit参数估计方法 | |
Wan et al. | High Accuracy 2D‐DOA Estimation for Conformal Array Using PARAFAC | |
CN109031197B (zh) | 一种基于超视距传播模型的辐射源直接定位方法 | |
CN107085202A (zh) | 一种针对互耦误差的Capon稳健自适应波束形成方法及系统 | |
CN109116295A (zh) | 基于相控阵选取基线的无源测向算法 | |
CN108872930A (zh) | 扩展孔径二维联合对角化doa估计方法 | |
CN108828508B (zh) | 一种超视距辐射源直接定位模型的性能分析方法 | |
CN117459176A (zh) | 数字相控阵天线多方向的噪声调制方法 | |
CN108089161A (zh) | 一种基于功率反馈的天线阵列合成波束空域选择方法 | |
CN100585429C (zh) | 一种基于非直线天线阵列的无源通道校正方法 | |
CN105182332B (zh) | 一种二维波束锐化方法 | |
CN109932690A (zh) | 一种基于接收信号强度的多目标水声定位方法 | |
CN112995888B (zh) | 一种基于阵列天线的定位方法、系统、电子设备及存储介质 | |
CN112013833B (zh) | 一种基于深度神经网络的victs天线对星跟踪方法 |
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 |