CN113064132A - 一种基于连续信任函数的鲁棒雷达目标检测方法 - Google Patents
一种基于连续信任函数的鲁棒雷达目标检测方法 Download PDFInfo
- Publication number
- CN113064132A CN113064132A CN202110331384.5A CN202110331384A CN113064132A CN 113064132 A CN113064132 A CN 113064132A CN 202110331384 A CN202110331384 A CN 202110331384A CN 113064132 A CN113064132 A CN 113064132A
- Authority
- CN
- China
- Prior art keywords
- function
- probability density
- density function
- probability
- equal
- 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
Links
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/04—Systems determining presence of a target
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computing Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及雷达目标检测技术领域,公开了一种基于连续信任函数的鲁棒雷达目标检测方法,包括如下步骤:根据应用需求,通过数值近似法计算并存储近似精度的标准概率密度函数等投注概率表,以获得的雷达回波幅度点观测为基础,计算点观测条件观测似然函数,然后基于广义贝叶斯定理使用条件观测似然函数求得基本信任分配函数,依据基本信任分配函数计算各假设的投注概率值,将结果直接用于决策判决或输出给下一级融合决策。本发明可以充分利用连续信任函数理论在不确定性信息表征和推理方面的优势,具有强大的信息表征能力和雷达目标检测鲁棒性。
Description
技术领域
本发明涉及雷达目标检测技术领域,更具体地说,涉及一种基于连续信任函数的鲁棒雷达目标检测方法。
背景技术
雷达目标检测是在含有噪声的雷达回波中,完成对“目标存在与否”假设的检验判断。鉴于噪声的随机特性以及雷达目标回波幅度随着姿态、频率和极化等多因素影响,雷达目标检测常在贝叶斯理论框架下建模为NP(Neyman-Pearson)问题,即将虚警概率限定在一个有限的水平,使检测概率最大化。
实际应用及理论研究中,大多假设雷达目标由多个独立散射体构成,每个散射体的雷达散射截面积独立且不变,但位置随机。根据电磁散射理论和中心极限定理,此类目标经线性检波器输出的雷达目标回波幅度服从瑞利分布。在贝叶斯理论框架下,NP检测的最终依据就是当线性检波器输出电压大于等于某门限,则判定有目标;当输出电压小于某门限,则判定无目标。其中,门限值由目标的平均幅度、噪声功率和虚警概率三个参数决定。一方面,针对上述NP检测原理,在雷达对抗领域提出了不少欺骗式干扰策略,也即提高干扰源的功率,渐进地诱偏跟踪锁定波束,从而实现战术逃逸;另一方面,由于目标回波起伏特性,极可能导致检测方法做出过于自信的决策,从而存在信息损失。因此,上述NP方法缺乏检测鲁棒性。
针对该问题,本发明提出了一种基于连续信任函数的鲁棒雷达目标检测方法,利用连续信任函数在不确定性信息表征方面的优势,设计了一种鲁棒检测器,使其同时对模型噪声和外来干扰都具备一定的判断力。
发明内容
本发明要解决的技术问题是雷达目标的鲁棒检测问题。针对传统NP方法缺乏鲁棒性,无法区分外来目标或干扰,在高度不确定的条件下依然会给出一个明确结论的问题,提出一种基于连续信任函数的鲁棒雷达目标检测方法,可以更好地区分模型噪声、目标和外来干扰等,具有更丰富的信息表征能力。
本发明解决上述问题的技术方案如下:
一种基于连续信任函数的鲁棒雷达目标检测方法,具体包括以下步骤:
步骤1:先验信息获取:选择雷达回波幅度起伏统计模型,即标准概率密度函数,所述标准概率密度函数为系数为1的瑞利概率密度函数,通过数值近似的方式获得标准概率密度函数的等投注概率密度表,以列表形式存储下来供数值求解应用,其中,系数为1的瑞利概率密度函数即且随机变量ξ表示雷达回波幅度取值,rls是系数为1的标准瑞利概率密度函数符号标志;等投注概率密度是指两个瑞利概率密度函数值相等,即对应的等投注概率对满足条件:u和为随机变量ξ的具体取值且
步骤3:依据点观测计算条件观测似然函数其中θc表示第c类假设,当c=0时表示模型噪声,当c=1时表示待检测目标;计算条件观测似然函数时,基于步骤1中提前获取的等投注概率密度表,然后采用数值近似的方式求解,以解决等投注概率对不存在解析解的问题;
式中:S表示辨别框架Θ={θ0,θ1}的任意子集,即且表示集合S的补集;当S={θ0}时,mΘ[z](S)表示点观测z来自于模型噪声的基本信任质量;当S={θ1}时,mΘ[z](S)表示表示点观测z来自于待检测目标的基本信任质量;当时,mΘ[z](S)表示点观测z来自于外来目标或非模型噪声的基本信任质量;当S={θ0,θ1}时,mΘ[z](S)表示点观测z来自于模型噪声或待检测目标的基本信任质量;
步骤5:根据基本信任分配函数mΘ[z](S)计算投注概率值BetP[z](θc),计算公式为:
式中|S|表示子集S的元素数目;
步骤6:将基本信任分配函数mΘ[z](S)或投注概率值BetP[z](θc)直接用于决策判决或输出给下一级融合决策。
优选的,所述步骤1中的雷达回波幅度起伏统计模型还可以是瑞利混合概率密度函数、或高斯概率密度函数、或高斯混合概率密度函数、或莱斯概率密度函数、或多自由度中心式卡方概率密度函数等概率密度函数。
优选的,所述步骤3中的条件观测似然函数数值近似法中累加求和下限利用观测值的线性变换或其对应的等投注概率对求得,累加求和上限根据等投注概率密度表的近似精度确定。
优选的,所述步骤4中的基本信任分配函数是在辨别框架的任意子集上赋值。
与现有技术相比,本发明的技术效果如下:
第一,鲁棒性强,相比贝叶斯理论框架中的NP准则,新检测器可以更好地区分模型噪声、目标和外来干扰等,具有更丰富的信息表征能力。
第二,扩展性好,该检测器需要求得瑞利概率密度函数的等投注概率对,针对其不存在解析解的难题,设计提出了一种数值近似的解算方法,可根据实际应用需求选择不同精度的等投注概率密度表;此外,只要标准概率密度函数具有单一峰值特性,或是多个单一峰值概率密度函数的加权混合,都可以采用本发明内容中的技术方案完成数值近似解算。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细说明。
图1为本发明实施步骤流程图;
图2为标准瑞利概率密度函数曲线图;
图3为雷达模型噪声和目标回波幅度起伏概率密度函数和对应条件观测似然函数曲线图;
图4为基本信任分配函数图;
图5为投注概率曲线图。
具体实施方式
为了更好地说明本发明的技术方案,以下结合实例对本发明的实施方式作进一步描述,以便对如何应用本发明技术手段来解决问题有更加深刻的理解,达到更好地解决实际问题的目的,所举实例只用于解释本发明,并非用于限定本发明的范围。
图1为本发明一种基于连续信任函数的鲁棒雷达目标检测方法,具体包括以下步骤:
步骤1:先验信息获取:通过数值近似的方式,获得标准瑞利概率密度函数的等投注概率对其中且γ(u)≤1≤u,以列表形式存储下来供数值求解应用。如图2所示,标准瑞利概率密度函数的单一峰值所处位置在1处,由此可将[0,1]区间等间隔分成N+1等份,在区间(0,1)内有N个点,分别对应于再利用数值法在区间(1,+∞)内寻找到对应的ui,使其建立如表1所示的等投注概率对。值得指出的是,该等投注概率表随着N的不同而变化,N越大表示近似精度越高。
表1.标准瑞利概率密度函数的等投注概率表
式中:θc表示第c类假设;累加求和上限COVrls表示近似的截断值,由等投注概率密度表确定,即表1中ui的最大值;累加求和下限ηi由观测值的线性变换或其对应的等投注概率对确定;Δη表示步进长度,即表1中相邻ui的间隔;表示第c类假设对应的条件瑞利概率密度函数系数;当c=0时,θ0表示模型噪声;当c=1时,θ1表示待检测目标;
式中:S表示辨别框架Θ={θ0,θ1}的任意子集,即当S={θ0}时,mΘ[z](S)表示点观测z来自于模型噪声的基本信任质量;当S={θ1}时,mΘ[z](S)表示表示点观测z来自于待检测目标的基本信任质量;当时,mΘ[z](S)表示点观测z来自于外来目标或非模型噪声的基本信任质量;当S={θ0,θ1}时,mΘ[z](S)表示点观测z来自于模型噪声或待检测目标的基本信任质量;
步骤5:根据基本信任分配函数mΘ[z](S)计算投注概率值BetP[z](θc),计算公式为:
式中|S|表示子集S的元素数目;
步骤6:将基本信任分配函数mΘ[z](S)或投注概率值BetP[z](θc)直接用于决策判决或输出给下一级融合决策。
本发明的效果通过以下实验进一步说明:
1.实验场景:
雷达回波幅度观测值为其模型噪声条件观测似然函数值为和目标条件观测似然函数值为辨别框架Θ={θ0,θ1}中的单一元素子集,其中θ0表示模型噪声假设,θ1表示待检测目标假设。设模型噪声的雷达回波幅度服从的标准瑞利概率密度函数,而待检测目标的雷达回波幅度服从的瑞利概率密度函数。
2.实验内容:
2.2)进一步按照本发明的步骤4,可以画出这些点观测所对应的基本信任分配函数,如图4所示。
2.3)最后按照本发明的步骤5,可以计算出其对应的投注概率值,如图5所示。
在贝叶斯理论框架下,由于模型噪声和待检测目标的先验概率未知,因此不同假设的概率与其条件概率密度函数成正比。
由图3可以看出,信任函数理论框架下的条件观测似然函数并不等于贝叶斯理论框架下的瑞利概率密度函数。
由图4可以看出,当点观测值z极小或者极大时,赋予空集的基本信任质量mΘ[z](S)都趋近于1,也即由于模型噪声和待检测目标的条件观测似然都很小,也即意味着模型噪声和待检测目标都不存在的可能性或外来干扰的可能性急速上升。当点观测z在模型噪声条件幅度起伏概率密度函数峰值位置附近时,赋予模型噪声假设S={θ0}的基本信任质量mΘ[z](S)相对最大;当点观测z在目标条件幅度起伏概率密度函数峰值位置附近时,赋予待检测目标假设S={θ1}的基本信任质量mΘ[z](S)相对最大;当点观测z在前述两个峰值之间时,模型噪声与待检测目标间的不可区分度又增加,因此赋予全集S={θ0,θ1}的基本信任质量mΘ[z](S)达到最大值。由此可以看出,基于连续信任函数的雷达目标检测器具有更强大的信息表征能力。
由图5曲线可以看出,两种不同的理论框架下得到的结果还是有区别的,其中P(θ0|z)和P(θ1|z)分别表示在贝叶斯理论框架中假设先验概率均匀分布,求解得到的后验概率。与贝叶斯理论相比,连续信任函数理论不需要对未知的先验概率进行配置,而且具有更丰富的信息表征能力,因此具有更高的鲁棒性
因此,通过瑞利目标检测实验验证了本发明方法的鲁棒性,本发明中的一种基于连续信任函数的鲁棒雷达目标检测方法可以分利用连续信任函数理论在不确定性信息表征和推理方面的优势实现鲁棒检测器。
上述说明示出并描述了发明应用的实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是其他实施例的排除,而可用于各种其他组合、修改和环境,并能在本文所述发明构想范围内,通过上述教导或相关领域的技术进行改动。而本领域人员所进行的改动和变化不脱离发明的精神和范围,则都应在发明所附权利要求的保护范围内。
Claims (4)
1.一种基于连续信任函数的鲁棒雷达目标检测方法,其特征在于,具体包括以下步骤:
步骤1:先验信息获取:选择雷达回波幅度起伏统计模型,即标准概率密度函数,所述标准概率密度函数为系数为1的瑞利概率密度函数,通过数值近似的方式获得标准概率密度函数的等投注概率密度表,以列表形式存储下来供数值求解应用,其中,系数为1的瑞利概率密度函数即且随机变量ξ表示雷达回波幅度取值,rls是系数为1的标准瑞利概率密度函数符号标志;等投注概率密度是指两个瑞利概率密度函数值相等,即对应的等投注概率对满足条件:u和为随机变量ξ的具体取值且
步骤3:依据点观测计算条件观测似然函数其中θc表示第c类假设,当c=0时表示模型噪声,当c=1时表示待检测目标;计算条件观测似然函数时,基于步骤1中提前获取的等投注概率密度表,然后采用数值近似的方式求解,以解决等投注概率对不存在解析解的问题;
式中:S表示辨别框架Θ={θ0,θ1}的任意子集,即且表示集合S的补集;当S={θ0}时,mΘ[z](S)表示点观测z来自于模型噪声的基本信任质量;当S={θ1}时,mΘ[z](S)表示表示点观测z来自于待检测目标的基本信任质量;当时,mΘ[z](S)表示点观测z来自于外来目标或非模型噪声的基本信任质量;当S={θ0,θ1}时,mΘ[z](S)表示点观测z来自于模型噪声或待检测目标的基本信任质量;
步骤5:根据基本信任分配函数mΘ[z](S)计算投注概率值BetP[z](θc),计算公式为:
式中|S|表示子集S的元素数目;
步骤6:将基本信任分配函数mΘ[z](S)或投注概率值BetP[z](θc)直接用于决策判决或输出给下一级融合决策。
2.根据权利要求1所述一种基于连续信任函数的鲁棒雷达目标检测方法,其特征在于,所述步骤1中的雷达回波幅度起伏统计模型还可以是瑞利混合概率密度函数、或高斯概率密度函数、或高斯混合概率密度函数、或莱斯概率密度函数、或多自由度中心式卡方概率密度函数等概率密度函数。
3.根据权利要求1所述一种基于连续信任函数的鲁棒雷达目标检测方法,其特征在于,所述步骤3中的条件观测似然函数数值近似法中累加求和下限利用观测值的线性变换或其对应的等投注概率对求得,累加求和上限根据等投注概率密度表的近似精度确定。
4.根据权利要求1所述一种基于连续信任函数的鲁棒雷达目标检测方法,其特征在于,所述步骤4中的基本信任分配函数是在辨别框架的任意子集上赋值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110331384.5A CN113064132B (zh) | 2021-03-29 | 2021-03-29 | 一种基于连续信任函数的鲁棒雷达目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110331384.5A CN113064132B (zh) | 2021-03-29 | 2021-03-29 | 一种基于连续信任函数的鲁棒雷达目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113064132A true CN113064132A (zh) | 2021-07-02 |
CN113064132B CN113064132B (zh) | 2022-08-16 |
Family
ID=76564179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110331384.5A Active CN113064132B (zh) | 2021-03-29 | 2021-03-29 | 一种基于连续信任函数的鲁棒雷达目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113064132B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113786176A (zh) * | 2021-08-17 | 2021-12-14 | 中国电子科技南湖研究院 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137424A (zh) * | 2015-07-17 | 2015-12-09 | 电子科技大学 | 一种杂波背景下实波束扫描雷达角超分辨方法 |
CN106324576A (zh) * | 2016-07-29 | 2017-01-11 | 西安电子科技大学 | 基于自回归模型的雷达目标自适应匹配滤波方法 |
CN112163373A (zh) * | 2020-09-23 | 2021-01-01 | 中国民航大学 | 基于贝叶斯机器学习的雷达系统性能指标动态评估方法 |
-
2021
- 2021-03-29 CN CN202110331384.5A patent/CN113064132B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105137424A (zh) * | 2015-07-17 | 2015-12-09 | 电子科技大学 | 一种杂波背景下实波束扫描雷达角超分辨方法 |
CN106324576A (zh) * | 2016-07-29 | 2017-01-11 | 西安电子科技大学 | 基于自回归模型的雷达目标自适应匹配滤波方法 |
CN112163373A (zh) * | 2020-09-23 | 2021-01-01 | 中国民航大学 | 基于贝叶斯机器学习的雷达系统性能指标动态评估方法 |
Non-Patent Citations (6)
Title |
---|
YUANYUAN LI等: "Distributed ISAR Subimage Fusion of Nonuniform Rotating Target Based on Matching Fourier Transform", 《SENSORS》 * |
宋莹等: "海杂波幅度分布模型的研究", 《现代电子技术》 * |
曾浩等: "复合高斯海杂波环境下雷达虚警概率分析", 《雷达科学与技术》 * |
李玮杰 等: "一种噪声环境下的雷达目标高分辨率距离像鲁棒识别方法", 《雷达学报》 * |
殷洁等: "基于反辐射导引头信息的近炸引信目标识别研究", 《电光与控制》 * |
龚树凤: "不确定条件下的机会阵雷达信号处理技术研究", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113786176A (zh) * | 2021-08-17 | 2021-12-14 | 中国电子科技南湖研究院 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
CN113786176B (zh) * | 2021-08-17 | 2024-03-15 | 中国电子科技南湖研究院 | 精确的毫米波雷达呼吸心跳测量方法、系统和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113064132B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Multi-target tracking algorithm based on PHD filter against multi-range-false-target jamming | |
Ayazgok et al. | Automatic antenna scan type classification for next‐generation electronic warfare receivers | |
CN113064132B (zh) | 一种基于连续信任函数的鲁棒雷达目标检测方法 | |
JP2018205175A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
CN103852747A (zh) | 增大雷达信号的目标探测概率的方法 | |
Ivković et al. | A new model of CFAR detector | |
Kim et al. | Radar scan pattern analysis for reduction of false identification in electronic warfare support systems | |
US8188910B2 (en) | Method for the multidimensional temporal kinematic filtering of radar blips, from antenna revolution to antenna revolution | |
CN107024680A (zh) | 一种基于双删除门限的目标检测方法 | |
Xie et al. | A novel method for deinterleaving radar signals: First‐order difference curve based on sorted TOA difference sequence | |
Sahin et al. | A multi-dimensional Hough transform algorithm based on unscented transform as a track-before-detect method | |
Wiley et al. | Performance comparison of PRF schedules for medium PRF radar | |
CN110031807A (zh) | 一种基于无模型强化学习的多阶段灵巧噪声干扰方法 | |
Musicki et al. | Linear multitarget finite resolution tracking in clutter | |
Tabet et al. | A generalized switching CFAR processor based on test cell statistics | |
Bae et al. | Joint estimation of multi‐target signal‐to‐noise ratio and dynamic states in cluttered environment | |
Ivković et al. | False alarm analysis of the CATM-CFAR in presence of clutter edge | |
Feng et al. | Detection performance analysis of recurrence quantification analysis measures for low observable target within sea clutter under different sea conditions | |
Tien et al. | An adaptive 2D-OS-CFAR thresholding in clutter environments: Test with real data | |
Yan-Juan et al. | A study on effectiveness modeling of multi-false-target jamming | |
CN113189557B (zh) | 对海雷达目标检测精细化处理方法及装置 | |
Kong et al. | A homogenous reference cells selector for CFAR detector in highly heterogeneous environment | |
Kim et al. | Spatial clutter measurement density estimation in nonhomogeneous measurement spaces | |
Cao et al. | An approach for high resolution radar target recognition based on BP neural network | |
Gilliam et al. | Covariance cost functions for scheduling multistatic sonobuoy fields |
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 |