CN105301593B - 一种基于双频rcs信息融合的飞行目标高度估计方法 - Google Patents
一种基于双频rcs信息融合的飞行目标高度估计方法 Download PDFInfo
- Publication number
- CN105301593B CN105301593B CN201510325679.6A CN201510325679A CN105301593B CN 105301593 B CN105301593 B CN 105301593B CN 201510325679 A CN201510325679 A CN 201510325679A CN 105301593 B CN105301593 B CN 105301593B
- Authority
- CN
- China
- Prior art keywords
- rcs
- frequency
- information
- height
- 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.)
- Expired - Fee Related
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
- G01S13/00—Systems 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/933—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft
- G01S13/935—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft for terrain-avoidance
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于双频RCS信息融合的飞行目标高度估计方法,包括:通过融合高频地波雷达双频信号回波强度信息,构造单基地体制高频地波雷达目标飞行高度和RCS信息完全可观测融合估计模型。在高频地波雷达双频工作模式下,对两个工作频率下低空飞行目标RCS信息进行组合和转换,实现双频信号回波信息的融合,提高了高度估计模型的可观测性,实现飞行目标高度和RCS的精确估计。
Description
技术领域
本发明涉及雷达技术领域,尤其涉及一种基于双频RCS信息融合的飞行目标高度估计方法。
背景技术
高频地波雷达因其独特的点波传播衰减特性,可以超视距探测到地平线以外的舰船和飞机目标。
高频地波雷达因波长较长,天线阵地较大,难以在俯仰方向上构建垂直阵列,不能在俯仰方向上获得目标角度信息,因而不能区别两个处于同一距离但不同高度上的飞行目标。在工程应用中,还没有实现高频地波雷达飞行目标的高度估计。
现有目标飞行高度估计技术中,信号回波强度是目标高度和RCS的函数,存在多解的问题。现有高度估计技术只适用于低空飞行目标,不适用于高空飞行目标的高度估计。另外,现有技术都是在假定目标的高低空属性是已知的情况下进行高度估计的,但是真实情况下,目标的高低空属性信息却是未知的,因此得到的高度估计结果的准确性较差。
发明内容
针对现有技术中存在的问题,本发明的目的为提供一种基于双频RCS信息融合的飞行目标高度估计方法,以解决现有技术中高度估计模型中多解以及估计结果准确性较差的技术问题。
为实现上述目的,本发明提供一种基于双频RCS信息融合的飞行目标高度估计方法,包括:
在高频地波雷达双频工作模式下,对两个工作频率下低空飞行目标RCS信息进行组合和转换,并对双频信号回波信息进行融合;
利用双频模式下的两个信号回波强度构造目标飞行高度和RCS信息完全可观测融合估计模型;
通过状态滤波方法对目标飞行高度和RCS信息进行精确估计。
进一步的,假定高频地波雷达双频工作频率为f1和f2,那么在两个频率下的RCS比值定义为
其中,θ=1°~180°。
进一步的,当得到某一频率下的目标RCS估计值,通过缓慢变化的RCS比值信息可以得到另一频率下目标的RCS信息。即假设得到估计值那么可以近似为
对于不同的飞机目标,可以通过数值建模得到随姿态角度变化的RCS比值信息库,RCS比值信息为两个不同频率下的RCS信息转换提供个了可能。
进一步的,所述高度和RCS融合估计模型的目标状态向量为
其中,h(k),vh(k)分别为k时刻目标飞行高度和高度方向上的变化速度,σt(k),为k时刻目标相对于T/R站的RCS和RCS变化速度。
进一步的,在高频波段,目标RCS相对于姿态角是一个慢变化信息,定义状态方程为
X(k)=FX(k-1)+v(k-1)
其中,v(k-1)为零均值高斯白噪声,过程协方差为Q。状态转移矩阵表示为
进一步的,以地波雷达两个频率下的信号回波强度为量测值,定义量测方程为
其中,是工作频率为f1下的信号回波强度,是工作频率为f2下的信号回波强度,具体形式分别为
w(k)为均值为零,方差为R(k)的高斯白噪声。σ(f1)与σ(f2)的转换关系如下
σ(f2,θ)=σc(f1,f2)·σ(f1,θ);
转换系数σc(f1,f2)可通过对不同飞机用矩量法建模分析得到。那么,可以表示为
进一步的,以和做为量测值的量测向量是目标状态X(k)的非线性函数,利用扩展卡尔曼滤波算法对飞机目标飞行高度和RCS进行估计,其迭代形式为
状态预测和量测预测:
其中,H(k)为线性化量测矩阵
协方差预测:
P(k/k-1)=FP(k-1/k-1)F'+Q(k)
量测预测协方差矩阵:
S(k)=HkP(k/k-1)H'k+R(k)
增益:
K(k)=P(k/k-1)H'kS-1(k)
状态更新方程:
更新协方差方程:
P(k/k)=[I-K(k)Hk]P(k/k-1)[I+K(k)Hk]'-K(k)R(k)K(k)'
本发明的有益效果在于,对两个工作频率下低空飞行目标RCS信息进行组合和转换,实现双频信号回波信息的融合,提高了高度估计模型的可观测性,实现飞行目标高度和RCS的精确估计。
附图说明
图1为目标在低空区域不同高度上的传播衰减曲线图。
图2为本发明实施提供的一种双频RCS信息融合的飞行目标高度估计方法的步骤流程图。
图3为本发明实施例中提供的目标RCS随频率变化的特性图。
具体实施方式
体现本发明特征与优点的典型实施例将在以下的说明中详细叙述。应理解的是,本发明能够在不同的实施例上具有各种的变化,其皆不脱离本发明的范围,且其中的说明及附图在本质上是当作说明之用,而非用以限制本发明。
本发明实施例一提供了一种估计飞行目标高度方法,步骤流程如图2所示,包括以下步骤:
步骤S1、在高频地波雷达双频工作模式下,对两个工作频率下低空飞行目标RCS信息进行组合和转换,并对双频信号回波信息进行融合;
步骤S2、利用双频模式下的两个信号回波强度构造目标飞行高度和RCS信息完全可观测融合估计模型;
步骤S3、通过状态滤波方法对目标飞行高度和RCS信息进行精确估计。
步骤S1中假定高频地波雷达双频工作频率为f1和f2,那么在两个频率下的RCS比值定义为
其中,θ=1°~180°。
进一步的,目当得到某一频率下的目标RCS估计值,通过缓慢变化的RCS比值信息可以得到另一频率下目标的RCS信息。即假设得到估计值那么可以近似为
对于不同的飞机目标,可以通过数值建模得到随姿态角度变化的RCS比值信息库,RCS比值信息为两个不同频率下的RCS信息转换提供个了可能。
进一步的,所述高度和RCS融合估计模型的目标状态向量为
其中,h(k),vh(k)分别为k时刻目标飞行高度和高度方向上的变化速度,σt(k),为k时刻目标相对于T/R站的RCS和RCS变化速度。
进一步的,在高频波段,目标RCS相对于姿态角是一个慢变化信息,定义状态方程为
X(k)=FX(k-1)+v(k-1)
其中,v(k-1)为零均值高斯白噪声,过程协方差为Q。状态转移矩阵表示为
进一步的,以地波雷达两个频率下的信号回波强度为量测值,定义量测方程为
其中,是工作频率为f1下的信号回波强度,是工作频率为f2下的信号回波强度,具体形式分别为
w(k)为均值为零,方差为R(k)的高斯白噪声。σ(f1)与σ(f2)的转换关系如下
σ(f2,θ)=σc(f1,f2)·σ(f1,θ)
转换系数σc(f1,f2)可通过对不同飞机用矩量法建模分析得到。那么,可以表示为
进一步的,以和做为量测值的量测向量是目标状态X(k)的非线性函数,利用扩展卡尔曼滤波算法对飞机目标飞行高度和RCS进行估计,其迭代形式为
状态预测和量测预测:
其中,H(k)为线性化量测矩阵
协方差预测:
P(k/k-1)=FP(k-1/k-1)F'+Q(k)
量测预测协方差矩阵:
S(k)=HkP(k/k-1)H'k+R(k)
增益:
K(k)=P(k/k-1)H'kS-1(k)
状态更新方程:
更新协方差方程:
P(k/k)=[I-K(k)Hk]P(k/k-1)[I+K(k)Hk]′-K(k)R(k)K(k)'
在高频地波雷达双频工作模式下,对两个工作频率下低空飞行目标RCS信息进行组合和转换,实现双频信号回波信息的融合,提高了高度估计模型的可观测性,实现飞行目标高度和RCS的精确估计。
本领域技术人员应当意识到在不脱离本发明所附的权利要求所揭示的本发明的范围和精神的情况下所作的改动与润饰,均属本发明的权利要求的保护范围之内。
Claims (2)
1.一种基于双频RCS信息融合的飞行目标高度估计方法,其特征在于,包括:
在高频地波雷达双频工作模式下,对两个工作频率下低空飞行目标RCS信息进行组合和转换,并对双频信号回波信息进行融合;
利用双频模式下的两个信号回波强度构造目标飞行高度和RCS融合估计模型;
通过状态滤波方法对目标飞行高度和RCS信息进行精确估计;
假定高频地波雷达双频工作频率为f1和f2,那么在两个频率下的RCS比值定义为
其中,θ=1°~180°;
当得到某一频率下的目标RCS估计值,通过缓慢变化的RCS比值信息可以得到另一频率下目标的RCS信息;即假设得到估计值那么近似为
对于不同的飞机目标,通过数值建模得到随姿态角度变化的RCS比值信息库,RCS比值信息为两个不同频率下的RCS信息转换提供了可能。
2.如权利要求1所述的方法,其特征在于,所述高度和RCS融合估计模型的目标状态向量为
其中,h(k),vh(k)分别为k时刻目标飞行高度和高度方向上的变化速度,σt(k),为k时刻目标相对于T/R站的RCS和RCS变化速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510325679.6A CN105301593B (zh) | 2015-06-12 | 2015-06-12 | 一种基于双频rcs信息融合的飞行目标高度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510325679.6A CN105301593B (zh) | 2015-06-12 | 2015-06-12 | 一种基于双频rcs信息融合的飞行目标高度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105301593A CN105301593A (zh) | 2016-02-03 |
CN105301593B true CN105301593B (zh) | 2017-10-10 |
Family
ID=55199072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510325679.6A Expired - Fee Related CN105301593B (zh) | 2015-06-12 | 2015-06-12 | 一种基于双频rcs信息融合的飞行目标高度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105301593B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107271992A (zh) * | 2017-06-01 | 2017-10-20 | 无锡孚嘉航海科技有限公司 | 一种高频雷达空中目标高度估测系统及其估测方法 |
CN111339493B (zh) * | 2020-03-06 | 2022-05-24 | 电子科技大学 | 一种非均匀量化状态转移特征矩阵提取方法 |
CN112162277B (zh) * | 2020-09-29 | 2022-05-24 | 哈尔滨工业大学 | 一种基于混合传播高频地波雷达的飞机目标高度估计方法 |
CN113391287B (zh) * | 2021-06-10 | 2023-09-01 | 哈尔滨工业大学 | 基于时间序列的高频地波雷达海态数据融合方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941296A (zh) * | 2014-04-04 | 2014-07-23 | 中国人民解放军理工大学 | 一种双频复合天线探地雷达 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7460052B2 (en) * | 2004-01-20 | 2008-12-02 | Bae Systems Information And Electronic Systems Integration Inc. | Multiple frequency through-the-wall motion detection and ranging using a difference-based estimation technique |
-
2015
- 2015-06-12 CN CN201510325679.6A patent/CN105301593B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941296A (zh) * | 2014-04-04 | 2014-07-23 | 中国人民解放军理工大学 | 一种双频复合天线探地雷达 |
Non-Patent Citations (3)
Title |
---|
The advantage of dual-frequency operation in ship tracking by HF surface wave radar;Hank Leong;《Proc. of the International Radar Conference》;20041231;1-6 * |
双基地高频地波雷达飞行目标高度估计研究;赵孔瑞;《中国博士学位论文全文数据库 信息科技辑》;20141215;I136-96 * |
双频率高频地波雷达船只目标点迹关联与融合处理;纪永刚 等;《系统工程与电子技术》;20140228;266-271 * |
Also Published As
Publication number | Publication date |
---|---|
CN105301593A (zh) | 2016-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105301593B (zh) | 一种基于双频rcs信息融合的飞行目标高度估计方法 | |
JP6328789B2 (ja) | レーダー警報受信機における到来角(aoa)を判定する方法及び装置 | |
CN108536171A (zh) | 一种多约束下多无人机协同跟踪的路径规划方法 | |
Liu et al. | An adaptive cubature Kalman filter algorithm for inertial and land-based navigation system | |
CN104035095B (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
EP2656101B1 (en) | Target altitude estimation based on measurements obtained by means of a passive radar | |
US9864380B1 (en) | Formation flight system extremum-seeking-control using blended performance parameters | |
CN110209180B (zh) | 一种基于HuberM-Cubature卡尔曼滤波的无人水下航行器目标跟踪方法 | |
CN110231620B (zh) | 一种噪声相关系统跟踪滤波方法 | |
CN105182326B (zh) | 一种利用方位信息的目标跟踪快速方法及装置 | |
CN107818324B (zh) | 将功能性网格单元添加到随机稀疏树形网格以便空间滤波的系统和方法 | |
CN107390166B (zh) | 一种自适应干扰源定位飞行校验方法 | |
RU130410U1 (ru) | Радиолокационное устройство идентификации воздушных объектов | |
CN105116393A (zh) | 一种基于位置指纹的高空目标飞行高度和雷达截面积估计方法 | |
CN107621632A (zh) | 用于nshv跟踪滤波的自适应滤波方法及系统 | |
CN104914421A (zh) | 基于和差波束的低空风切变风速估计方法 | |
CN108107416B (zh) | 一种二次雷达低空探测杂波抑制方法 | |
CN116106824A (zh) | 一种基于认知学习的无人机多阶段信号源定位方法和系统 | |
CN112835034B (zh) | 一种双通道雷达对地测高系统及方法 | |
CN104931932A (zh) | 一种改进的去偏坐标转换卡尔曼滤波方法 | |
CN116225053A (zh) | 一种面向任务分级目标运动特性的跟踪轨迹优化方法 | |
CN111679270B (zh) | 一种反射点不确定场景下多路径融合目标检测算法 | |
Kohlleppel | Ground target tracking with signal adaptive measurement error covariance matrix | |
CN113093171A (zh) | 基于目标跟踪的机载雷达路径与辐射资源联合优化方法 | |
CN108363055B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20190226 Address after: 266000 No. 98 Xiangjiang Road, Huangdao District, Qingdao City, Shandong Province Patentee after: China Electronics Technology Instrument and Meter Co., Ltd. Address before: 266000 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: The 41st Institute of CETC |
|
TR01 | Transfer of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171010 Termination date: 20200612 |
|
CF01 | Termination of patent right due to non-payment of annual fee |