CN113723253A - 基于杂波预分类的运动目标检测方法 - Google Patents
基于杂波预分类的运动目标检测方法 Download PDFInfo
- Publication number
- CN113723253A CN113723253A CN202110969902.6A CN202110969902A CN113723253A CN 113723253 A CN113723253 A CN 113723253A CN 202110969902 A CN202110969902 A CN 202110969902A CN 113723253 A CN113723253 A CN 113723253A
- Authority
- CN
- China
- Prior art keywords
- distribution
- vector
- clutter
- detected
- matrix
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 73
- 238000009826 distribution Methods 0.000 claims abstract description 138
- 239000013598 vector Substances 0.000 claims abstract description 138
- 239000011159 matrix material Substances 0.000 claims abstract description 100
- 238000000034 method Methods 0.000 claims abstract description 47
- 238000013528 artificial neural network Methods 0.000 claims abstract description 37
- 238000012549 training Methods 0.000 claims abstract description 34
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 14
- 238000011176 pooling Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000002372 labelling Methods 0.000 claims description 3
- 210000002569 neuron Anatomy 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 2
- 230000004913 activation Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 230000001629 suppression Effects 0.000 abstract description 4
- 238000007781 pre-processing Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 36
- 238000012545 processing Methods 0.000 description 8
- 238000012360 testing method Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 229910003460 diamond Inorganic materials 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000001537 neural effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
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
- 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/50—Systems of measurement based on relative movement of 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Databases & Information Systems (AREA)
- Computational Linguistics (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
一种基于杂波预分类的运动目标检测方法,其步骤是:对建立的训练集进行预处理;设置残差神经网络各层的初始参数;生成待检向量集和辅助向量集;将雷达杂波矩阵中的待检向量进行预分类,采用均值估计法、韦‑对协方差估计法和更新协方差估计法,分别计算预分类后的瑞利分布、韦布尔分布、对数正态分布和K分布的协方差矩阵,应用所估计的协方差矩阵计算待检向量的自适应检测统计量;根据检测统计量的大小确定是否存在运动目标。本发明提高了非匀质杂波下的运动目标检测性能,可应用于在实际复杂杂波目标检测场景发生变化时的杂波抑制和自适应运动目标检测。
Description
技术领域
本发明属于通信技术领域,更进一步涉及信号处理技术领域中的一种基于杂波预分类的运动目标检测方法。本发明可用于雷达对实际复杂杂波场景中运动目标的自适应检测。
背景技术
实际复杂杂波场景中的运动目标检测是对雷达接收的回波数据进行信号建模,采用信号处理技术对雷达工作场景中感兴趣的目标进行检测。实际复杂杂波场景中的运动目标检测是雷达应用技术的重要组成部分,可以在卫星定位和民航控制中,识别目标,为目标定位、追踪提供精准信息,在军用和民用领域至关重要。目前,实际复杂杂波场景中的运动目标检测方法主要有广义似然比、Rao、自适应匹配滤波方法等,但上述方法中用来估计杂波协方差矩阵的雷达接收回波数据均未经杂波幅度模型预分类处理,导致杂波信号建模不够准确,给杂波抑制和运动目标检测带来困难。
何友,简涛,苏峰等人在其发表的论文“非高斯杂波协方差矩阵估计方法及CFAR特性分析”(《中国科学(信息科学)》)中公开了一种运动目标自适应检测方法。该方法利用雷达回波数据对未知的杂波协方差矩阵进行估计,以求得具有恒虚警率CFAR(ConstantFalse Alarm Rate)特性的目标检测器,并以上述检测器完成对实际杂波场景中运动目标的自适应检测。该方法采用的两种不同杂波协方差矩阵估计方法分别为:采样协方差矩阵SCM(Sample Covariance Matrix)和归一化采样协方差矩阵NSCM(Normalized SampleCovariance Matrix)。但是,该方法存在的不足之处是,两种协方差矩阵估计方法都仅针对一种杂波幅度分布进行杂波信号建模,当实际场景中存在服从其他幅度分布的杂波信号时,所建立的模型无法根据杂波特性进行灵活调节,导致杂波建模不够准确,杂波抑制效果不理想,该目标检测器并不能实现完全意义上的自适应目标检测。
中国人民解放军海军航空工程学院在其申请的专利文献“基于变参数广义结构的距离扩展目标自适应检测方法”(申请号:201710284894.5,申请公布号:CN 106932766 A)中提出了一种运动目标自适应检测方法。该方法首先利用杂波信号的非高斯特征信息构建杂波协方差矩阵估计模型,然后利用特定杂波环境中现有运动目标检测统计量的共性特征和上述矩阵估计模型共同构建运动目标自适应检测器,实现雷达在实际杂波场景中对运动目标的自适应检测。但是,该方法仍然存在的不足之处是,在建立杂波协方差矩阵估计模型时,仅根据当前运动目标检测场景中杂波的非高斯特征信息进行信号建模,当运动目标检测场景发生变化时,所建立的杂波协方差矩阵估计模型难以快速与发生变化的运动目标检测场景相匹配,导致雷达信号处理实时性差,运动目标检测性能下降。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种基于杂波预分类的运动目标检测方法,旨在解决非匀质杂波环境下和目标检测场景发生变化时杂波抑制效果不理想和运动目标检测性能差的问题。
实现本发明目的的思路是,利用训练集中已知幅度特性的雷达杂波矩阵,在训练网络的过程中逐步更新残差神经网络的权重参数,使得网络能模仿人脑神经元的连接以及神经信号的传导,从而能充分考虑雷达杂波数据中可能存在的待检向量类型,并将其准确的分为瑞利分布、韦布尔分布、对数正态分布和K分布四类,再针对不同分布的待检向量分别采用均值估计法、韦-对协方差估计法和更新协方差估计法准确的估计其协方差矩阵;最后根据所估计的准确协方差矩阵计算每个待检向量的自适应检测统计量,实现对运动目标的检测。
其实现步骤包括如下:
步骤1,生成训练集:
(1a)选取至少10000个已知幅度特性的雷达杂波矩阵组成数据集;对矩阵集中的每个雷达杂波矩阵进行标注,将所有的标签文件组成标签集;将数据集和标签集组成训练样本集;
(1b)在区间[0,1]内,为训练样本集中的每个矩阵随机生成一个实数;
(1c)将对应实数值在区间[0.5,1]内的每个矩阵,分别以该矩阵的中心列和中心行为轴,依次将该矩阵在水平和垂直方向翻转后得到训练集;
步骤2,训练残差神经网络:
设置残差神经网络中各层的参数;将训练集输入到残差神经网络中,利用交叉熵损失函数,计算网络的输出与输入矩阵的标签之间的交叉熵损失值,通过随机梯度下降算法SGD迭代更新网络参数,直到交叉熵损失函数收敛为止,得到训练好的残差神经网络;
步骤3,生成待检向量和辅助向量集:
(3a)将实测雷达接收的杂波回波数据生成N×(K+L)维的回波矩阵,其中,N表示雷达在时间维发射的采样脉冲的总数,K表示辅助距离门的总数,L表示待检距离门的总数;
(3b)将杂波回波矩阵中对应待检距离门的每一列元素组成一个待检向量,其维度为N×1,将所有待检向量组成待检向量集;
(3c)将杂波回波矩阵中除待检向量所在列外的每一列元素组成一个辅助向量,其维度为N×1,将所有辅助向量组成辅助向量集;
步骤4,对雷达杂波矩阵进行预分类:
将待检向量集和辅助向量集分别输入到训练好的残差神经网络中,输出待检向量集和辅助向量集中每个向量的瑞利分布、对数正态分布、韦布尔分布和K分布四个概率值,将每个向量中最大概率值对应的分布类型作为该向量的分布类型;
步骤5,根据预分类的类型估计协方差矩阵:
(5a)从辅助向量集中选出所有服从瑞利分布的向量,利用均值估计公式,估计服从瑞利分布的协方差矩阵;
(5b)从辅助向量集中选出所有服从韦布尔分布或对数正态分布的向量,利用归一化估计公式,估计服从韦布尔分布或对数正态分布的协方差矩阵;
(5c)从辅助向量集中选出所有服从K分布的向量,采用迭代估计法,估计服从K分布的协方差矩阵;
步骤6,根据预分类的类型计算待检向量的自适应检测统计量:
(6a)按照下式,计算每个服从瑞利分布的待检向量的自适应检测统计量:
其中,表示第a个服从瑞利分布向量的自适应检测统计量,|·|表示取模操作,ya表示第a个服从瑞利分布的向量,a表示服从瑞利分布的向量的序号,(·)H表示取共轭转置操作,s表示导向矢量,s=[1,…,e-j·0.2·π·(N-1)],e(·)表示以e为底的指数操作,表示服从瑞利分布的协方差矩阵,(·)-1表示取逆操作;
(6b)按照下式,计算每个服从韦布尔分布或对数正态分布的待检向量的自适应检测统计量:
其中,表示第b个服从韦布尔分布或对数正态分布向量的自适应检测统计量,yb表示第b个服从韦布尔分布或对数正态分布的向量,b表示服从韦布尔分布或对数正态分布向量的序号,表示服从韦布尔分布或对数正态分布的协方差矩阵;
(6c)按照下式,计算每个服从K分布的待检向量的自适应检测统计量:
步骤7,对运动目标进行检测:
判定雷达杂波数据中的待检向量的自适应统计量的最高值是否大于或等于次高值的两倍,若是,则雷达杂波数据中存在运动目标,否则,雷达杂波数据中不存在运动目标。
本发明与现有的技术相比,具有以下优点:
第一,本发明利用残差神经网络对雷达杂波数据中的待检向量进行预分类,根据预分类结果准确的估计实际目标检测场景中的待检向量的协方差矩阵,克服了现有技术中所估计的协方差矩阵难以与发生变化的实际运动目标检测场景相匹配的问题,使得本发明可以根据目标检测场景灵活调节估计协方差矩阵的方式,提高运动目标检测性能。
第二,本发明将雷达杂波数据中的待检向量分为瑞利分布、K分布、韦布尔分布和对数正态分布四类,针对不同分布类型的待检向量,采用不同的协方差矩阵估计方法分别估计其协方差矩阵,克服了现有技术中只针对一种待检向量分布特性估计协方差矩阵,导致协方差矩阵估计不够准确,运动目标检测效果不理想的问题,使得本发明能准确的估计协方差矩阵,提高运动目标检测性能。
附图说明
图1为本发明的流程图;
图2为本发明仿真实验1的结果对比图;
图3为本发明仿真实验2的准确率的变化图。
具体实施方式
下面结合附图和实施例对本发明做进一步的描述。
结合图1,对本发明的实现步骤做进一步的描述。
步骤1,建立训练集。
选取至少10000个已知幅度特性的雷达杂波矩阵组成数据集。
对矩阵集中的每个雷达杂波矩阵进行标注,将所有的标签文件组成标签集。
将数据集和标签集组成训练集。
步骤2,对训练集进行预处理。
在区间[0,1]内,为训练集中的每个矩阵随机生成一个实数。
将对应实数值在区间[0.5,1]内的每个矩阵,分别以该矩阵的中心列和中心行为轴,依次将该矩阵在水平和垂直方向翻转。
步骤3,设置残差神经网络初始参数。
应用一个20层的残差神经网络并设置该网络各层的参数;该残差神经网络的结构依次为:输入层,卷积层,最大池化层,卷积模块,平均池化层,全连接层,Softmax分类器;所述卷积模块是由16个串联的卷积层组成。
设置残差神经网络各层的参数如下:将最大池化层和平均池化层的池化窗口均设置为3×3;卷积层的卷积核大小为7×7×64,步长为2;所述卷积模块中第一至第四卷积层的卷积核均设置为3×3×64,步长均为1;第五至第八卷积层的卷积核均设置为3×3×128,步长依次为2,1,1,1;第九至第十三卷积层的卷积核均设置为3×3×256,步长依次为2,1,1,1;第十四至第十七卷积层的卷积核均设置为3×3×512,步长依次为2,1,1,1;卷积层和卷积模块中的所有激活函数均为Relu函数;将全连接层的神经元个数设置为4。
步骤4,训练残差神经网络。
将训练集输入到残差神经网络中,利用交叉熵损失函数,计算网络的输出与输入矩阵的标签之间的交叉熵损失值,通过随机梯度下降算法SGD迭代更新网络参数,直到交叉熵损失函数收敛为止,得到训练好的残差神经网络。
所述交叉熵损失函数如下:
其中,Loss表示交叉熵损失函数,R表示训练集中矩阵的总数,∑表示求和操作,r表示训练集中矩阵的序号,zr表示训练集中第r个矩阵的类型标注,log表示以2为底的对数操作,sr表示将训练集中第r个矩阵经残差神经网络后的输出矩阵。
步骤5,生成待检向量和辅助向量集。
将实测雷达接收的杂波回波数据生成N×(K+L)维的回波矩阵,其中,N表示雷达在时间维发射的采样脉冲的总数,K表示辅助距离门的总数,L表示待检距离门的总数。
将杂波回波矩阵中对应待检距离门的每一列元素组成一个待检向量,其维度为N×1,将所有待检向量组成待检向量集。
将杂波回波矩阵中除待检向量所在列外的每一列元素组成一个辅助向量,其维度为N×1,将所有辅助向量组成辅助向量集。
步骤6,对雷达杂波矩阵进行分类。
将待检向量集和辅助向量集分别输入到已经训练好的残差神经网络中,再分别输出待检向量集和辅助向量集中每一个向量的瑞利分布、对数正态分布、韦布尔分布和的K分布四个概率值,将最大概率值对应的分布类型作为该向量的分布类型。
步骤7,根据分布类型估计协方差矩阵。
从辅助向量集中选出所有服从瑞利分布的向量,利用下述均值估计公式,估计服从瑞利分布的协方差矩阵:
其中,A表示服从瑞利分布的向量的总数,Σ表示求和操作。
从辅助向量集中选出所有服从韦布尔分布或对数正态分布的向量,利用下述韦-对协方差估计公式,估计服从韦布尔分布或对数正态分布的协方差矩阵:
其中,B表示服从韦布尔分布或对数正态分布的向量的总数,max{·}表示取最大值操作,min{·}表示取最小值操作。
从辅助向量集中选出所有服从K分布的向量,采用更新协方差估计法,估计服从K分布的协方差矩阵。
所述更新协方差估计法的步骤如下。
其中,D表示服从K分布向量的总数。
第四步,判断协方差矩阵迭代更新的次数是否大于或等于3,若是,则将当前迭代更新后的协方差矩阵作为服从K分布的协方差矩阵后执行步骤8,否则,执行本步骤的第三步。
步骤8,根据分布类型计算待检向量的自适应检测统计量。
按照下式,计算每个服从瑞利分布的待检向量的自适应检测统计量:
其中,表示第a个服从瑞利分布向量的自适应检测统计量,|·|表示取模操作,ya表示第a个服从瑞利分布的向量,a表示服从瑞利分布的向量的序号,(·)H表示取共轭转置操作,s表示导向矢量,s=[1,…,e-j·0.2·π·(N-1)],e(·)表示以e为底的指数操作,表示服从瑞利分布的协方差矩阵,(·)-1表示取逆操作。
按照下式,计算每个服从韦布尔分布或对数正态分布的待检向量的自适应检测统计量:
其中,表示第b个服从韦布尔分布或对数正态分布向量的自适应检测统计量,yb表示第b个服从韦布尔分布或对数正态分布的向量,b表示服从韦布尔分布或对数正态分布向量的序号,表示服从韦布尔分布或对数正态分布的协方差矩阵。
按照下式,计算每个服从K分布的待检向量的自适应检测统计量:
步骤9,对运动目标进行检测。
判定雷达杂波数据中的待检向量的自适应统计量的最高值是否大于或等于次高值的两倍,若是,则雷达杂波数据中存在运动目标,否则,雷达杂波数据中不存在运动目标。
下面结合仿真实验对本发明的效果做进一步的说明:
本发明的仿真实验有两个。
仿真实验1为对运动目标进行检测的仿真实验。
仿真实验1的仿真条件:
本发明仿真实验1的运动目标检测仿真实验的硬件平台为:处理器为Intel Corei5-8500 CPU,主频为3.00GHz,内存8GB。
本发明仿真实验1的运动目标检测仿真试验的软件平台为:Windows 10操作***、MATLAB R2016a。
本发明仿真实验1的运动目标检测仿真实验的参数设置为:雷达在时间维发射的采样脉冲的总数N=4,辅助距离门的总数K=11,待检距离门的总数L=34。
本发明仿真实验1的杂波环境为:申请人利用MATLAB仿真软件构造了一个服从正态高斯分布的噪声环境,信杂比为SNR=20dB。
仿真实验1的仿真内容及其结果分析:
本发明的运动目标检测仿真实验是采用本发明和两个现有技术(采样协方差方法和归一化采样协方差方法)分别对申请人构造的杂波环境下的运动目标自适应检测性能进行了仿真实验,得到了本发明方法和两个现有技术在各距离单元的检测统计量大小,并绘制各距离单元与检测统计量大小的关系曲线如图2所示。
在运动目标检测仿真实验中采用的两个现有技术是指:
现有技术采样协方差法是指,何友,简涛等人在“非高斯杂波协方差矩阵估计方法及CFAR特性分析”(中国科学,2011,41(1))中提出的采样协方差法。
现有技术归一化采样协方差法是指,刘明,水鹏朗等人在“基于功率中值和归一化采样协方差矩阵的自适应匹配滤波检测器”(电子与信息学报,2015,37(6))中提出的归一化采样协方差法。
下面结合图2的仿真图对本发明的效果做进一步的描述。
图2为利用本发明方法和两个现有技术分别得到的每个待检单元的检测统计量的结果对比图。图2中的横轴表示距离单元。图2中的纵轴表示检测统计量的大小。图2中以菱形符号标示的曲线表示采用本发明方法对各待检单元的检测统计量计算结果的曲线。图2中以三角形形符号标示的曲线表示采用采样协方差方法对各待检单元的检测统计量计算结果的曲线。图2中以圆形符号标示的曲线表示采用归一化采样协方差方法对各待检单元的检测统计量计算结果的曲线。
由图2的仿真结果图可见,在第17个待检单元处,本发明方法在该待检单元处的检测统计量明显高于两个现有技术在该待检单元处的检测统计量,由此可见,本发明能适应实际复杂杂波场景的变化,在复杂杂波背景下进行运动目标的检测,具有更优的运动目标检测性能。
本发明的仿真实验2为对杂波数据进行预分类的仿真实验。
本发明仿真实验2的仿真实验条件:
本发明仿真实验2的杂波预分类处理仿真试验的硬件平台为:处理器为IntelXeon 3070 GPU,主频为2.70GHz,内存8GB。
本发明仿真实验2的杂波预分类处理仿真试验的软件平台为:PYTORCH1.7.0。
本发明仿真实验2的杂波预分类处理仿真实验的参数设置为:batchsize=128,epoches=100。
仿真实验2的仿真内容及其结果分析:
本发明的预分类处理仿真实验是将雷达杂波矩阵输入到残差神经网络中,采用随机梯度下降算法迭代更新残差神经网络的参数100次。在每一次迭代完成后,利用混淆矩阵准确率计算方法,计算当前迭代后的残差神经网络的预分类准确率。将100次迭代后计算得到的残差神经网络的预分类准确率,绘制成残差神经网络的预分类准确率与迭代次数的关系曲线,如图3所示。
下面结合图3对本发明的效果做进一步的描述。
图3为残差神经网络在更新的过程中,对待检单元进行分类的准确率的变化图。图3中的横轴表示训练的迭代次数,图3中的纵轴表示残差神经网络预测分类准确率的大小,单位为百分率(%)。图3中以五角星符号标示的曲线表示神经网络在训练过程中准确率的变换曲线。
由图3的仿真结果图可见,残差神经网络的预分类准确率随着网络训练次数的增大而增大,残差神经网络的分类准确率由65%上升至94%。由此可见,当达到最大迭代次数时,残差神经网络能准确的进行待检向量的预分类,以便在实际场景中根据杂波特性的不同,采用相应的方法估计杂波的协方差矩阵,实现对杂波的充分抑制。
Claims (6)
1.一种基于杂波预分类的运动目标检测方法,其特征在于,利用残差神经网络对实际复杂杂波场景下的雷达杂波数据进行预分类,基于预分类结果分别对服从不同幅度分布的杂波数据建立不同的杂波协方差矩阵;该方法的步骤包括如下:
步骤1,生成训练集:
(1a)选取至少10000个已知幅度特性的雷达杂波矩阵组成数据集;对矩阵集中的每个雷达杂波矩阵进行标注,将所有的标签文件组成标签集;将数据集和标签集组成训练样本集;
(1b)在区间[0,1]内,为训练样本集中的每个矩阵随机生成一个实数;
(1c)将对应实数值在区间[0.5,1]内的每个矩阵,分别以该矩阵的中心列和中心行为轴,依次将该矩阵在水平和垂直方向翻转后得到训练集;
步骤2,训练残差神经网络:
设置残差神经网络中各层的参数;将训练集输入到残差神经网络中,利用交叉熵损失函数,计算网络的输出与输入矩阵的标签之间的交叉熵损失值,通过随机梯度下降算法SGD迭代更新网络参数,直到交叉熵损失函数收敛为止,得到训练好的残差神经网络;
步骤3,生成待检向量和辅助向量集:
(3a)将实测雷达接收的杂波回波数据生成N×(K+L)维的回波矩阵,其中,N表示雷达在时间维发射的采样脉冲的总数,K表示辅助距离门的总数,L表示待检距离门的总数;
(3b)将杂波回波矩阵中对应待检距离门的每一列元素组成一个待检向量,其维度为N×1,将所有待检向量组成待检向量集;
(3c)将杂波回波矩阵中除待检向量所在列外的每一列元素组成一个辅助向量,其维度为N×1,将所有辅助向量组成辅助向量集;
步骤4,对雷达杂波矩阵进行预分类:
将待检向量集和辅助向量集分别输入到训练好的残差神经网络中,输出待检向量集和辅助向量集中每个向量的瑞利分布、对数正态分布、韦布尔分布和K分布四个概率值,将每个向量中最大概率值对应的分布类型作为该向量的分布类型;
步骤5,根据预分类的类型估计协方差矩阵:
(5a)从辅助向量集中选出所有服从瑞利分布的向量,利用均值估计公式,估计服从瑞利分布的协方差矩阵;
(5b)从辅助向量集中选出所有服从韦布尔分布或对数正态分布的向量,利用归一化估计公式,估计服从韦布尔分布或对数正态分布的协方差矩阵;
(5c)从辅助向量集中选出所有服从K分布的向量,采用迭代估计法,估计服从K分布的协方差矩阵;
步骤6,根据预分类的类型计算待检向量的自适应检测统计量:
(6a)按照下式,计算每个服从瑞利分布的待检向量的自适应检测统计量:
其中,表示第a个服从瑞利分布向量的自适应检测统计量,|·|表示取模操作,ya表示第a个服从瑞利分布的向量,a表示服从瑞利分布的向量的序号,(·)H表示取共轭转置操作,s表示导向矢量,s=[1,…,e-j·0.2·π·(N-1)],e(·)表示以e为底的指数操作,表示服从瑞利分布的协方差矩阵,(·)-1表示取逆操作;
(6b)按照下式,计算每个服从韦布尔分布或对数正态分布的待检向量的自适应检测统计量:
其中,表示第b个服从韦布尔分布或对数正态分布向量的自适应检测统计量,yb表示第b个服从韦布尔分布或对数正态分布的向量,b表示服从韦布尔分布或对数正态分布向量的序号,表示服从韦布尔分布或对数正态分布的协方差矩阵;
(6c)按照下式,计算每个服从K分布的待检向量的自适应检测统计量:
步骤7,对运动目标进行检测:
判定雷达杂波数据中的待检向量的自适应统计量的最高值是否大于或等于次高值的两倍,若是,则雷达杂波数据中存在运动目标,否则,雷达杂波数据中不存在运动目标。
2.根据权利要求1所述的基于杂波预分类的运动目标检测方法,其特征在于,步骤(2)中所述残差神经网络的结构依次为:输入层,卷积层,最大池化层,卷积模块,平均池化层,全连接层,Softmax分类器;所述卷积模块是由16个串联的卷积层组成;
设置残差神经网络各层的参数如下:将最大池化层和平均池化层的池化窗口均设置为3×3;卷积层的卷积核大小为7×7×64,步长为2;所述卷积模块中第一至第四卷积层的卷积核均设置为3×3×64,步长均为1;第五至第八卷积层的卷积核均设置为3×3×128,步长依次为2,1,1,1;第九至第十三卷积层的卷积核均设置为3×3×256,步长依次为2,1,1,1;第十四至第十七卷积层的卷积核均设置为3×3×512,步长依次为2,1,1,1;卷积层和卷积模块中的所有激活函数均为Relu函数;将全连接层的神经元个数设置为4。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110969902.6A CN113723253B (zh) | 2021-08-23 | 2021-08-23 | 基于杂波预分类的运动目标检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110969902.6A CN113723253B (zh) | 2021-08-23 | 2021-08-23 | 基于杂波预分类的运动目标检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113723253A true CN113723253A (zh) | 2021-11-30 |
CN113723253B CN113723253B (zh) | 2023-06-20 |
Family
ID=78677358
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110969902.6A Active CN113723253B (zh) | 2021-08-23 | 2021-08-23 | 基于杂波预分类的运动目标检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113723253B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2949882A1 (fr) * | 2009-09-08 | 2011-03-11 | Thales Sa | Methode de traitement du signal comportant une etape d'estimation de la matrice mediane d'un echantillon de matrices |
US20120249361A1 (en) * | 2011-04-04 | 2012-10-04 | Zafer Sahinoglu | Method for Detecting Targets Using Space-Time Adaptive Processing |
CN104215944A (zh) * | 2014-09-22 | 2014-12-17 | 南京信息工程大学 | 一种基于参数选择的雷达海杂波自适应抑制处理方法 |
CN106468770A (zh) * | 2016-09-23 | 2017-03-01 | 西安电子科技大学 | K分布杂波加噪声下的近最优雷达目标检测方法 |
CN107102302A (zh) * | 2017-04-27 | 2017-08-29 | 中国人民解放军海军航空工程学院 | 基于非高斯参数的雷达目标融合检测方法 |
CN109116311A (zh) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 |
CN110763187A (zh) * | 2019-09-30 | 2020-02-07 | 中国科学院测量与地球物理研究所 | 一种稳健的基于雷达分布式目标的地面沉降监测方法 |
CN110879927A (zh) * | 2019-11-20 | 2020-03-13 | 中国海洋大学 | 一种用于海目标检测的海杂波幅度统计分布现场建模方法 |
CN111999714A (zh) * | 2020-09-02 | 2020-11-27 | 中国人民解放军海军航空大学 | 基于多散射点估计和杂波知识辅助的自适应融合检测方法 |
CN112147584A (zh) * | 2020-09-18 | 2020-12-29 | 长沙驰芯半导体科技有限公司 | 一种基于非均匀杂波的mimo雷达扩展目标检测方法 |
CN113256667A (zh) * | 2021-04-14 | 2021-08-13 | 大连海事大学 | 一种基于Otsu与联合分布的SAR图像舰船目标检测方法 |
-
2021
- 2021-08-23 CN CN202110969902.6A patent/CN113723253B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2949882A1 (fr) * | 2009-09-08 | 2011-03-11 | Thales Sa | Methode de traitement du signal comportant une etape d'estimation de la matrice mediane d'un echantillon de matrices |
US20120249361A1 (en) * | 2011-04-04 | 2012-10-04 | Zafer Sahinoglu | Method for Detecting Targets Using Space-Time Adaptive Processing |
CN104215944A (zh) * | 2014-09-22 | 2014-12-17 | 南京信息工程大学 | 一种基于参数选择的雷达海杂波自适应抑制处理方法 |
CN106468770A (zh) * | 2016-09-23 | 2017-03-01 | 西安电子科技大学 | K分布杂波加噪声下的近最优雷达目标检测方法 |
CN107102302A (zh) * | 2017-04-27 | 2017-08-29 | 中国人民解放军海军航空工程学院 | 基于非高斯参数的雷达目标融合检测方法 |
CN109116311A (zh) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 |
CN110763187A (zh) * | 2019-09-30 | 2020-02-07 | 中国科学院测量与地球物理研究所 | 一种稳健的基于雷达分布式目标的地面沉降监测方法 |
CN110879927A (zh) * | 2019-11-20 | 2020-03-13 | 中国海洋大学 | 一种用于海目标检测的海杂波幅度统计分布现场建模方法 |
CN111999714A (zh) * | 2020-09-02 | 2020-11-27 | 中国人民解放军海军航空大学 | 基于多散射点估计和杂波知识辅助的自适应融合检测方法 |
CN112147584A (zh) * | 2020-09-18 | 2020-12-29 | 长沙驰芯半导体科技有限公司 | 一种基于非均匀杂波的mimo雷达扩展目标检测方法 |
CN113256667A (zh) * | 2021-04-14 | 2021-08-13 | 大连海事大学 | 一种基于Otsu与联合分布的SAR图像舰船目标检测方法 |
Non-Patent Citations (5)
Title |
---|
NINGYUAN SU 等: "A deep learning method of moving target classification in clutter background", pages 1 - 9 * |
XIAOQIAN MOU 等: "Motion classification for radar moving target via STFT and convolution neural network", vol. 2019, no. 2019, pages 6287 - 6290, XP006085844, DOI: 10.1049/joe.2019.0179 * |
张彦飞;孙文杰;孙玉梅;孟祥伟;陈祥光;: "非高斯杂波中机载多进多出雷达动目标检测", vol. 47, no. 03, pages 23 - 31 * |
方明 等: "基于联合稀疏恢复的宽带雷达动目标检测方法", vol. 37, no. 37, pages 2977 - 2983 * |
薛健: "复合高斯海杂波背景雷达目标检测算法", vol. 2021, no. 2021, pages 136 - 50 * |
Also Published As
Publication number | Publication date |
---|---|
CN113723253B (zh) | 2023-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112001270B (zh) | 基于一维卷积神经网络的地面雷达自动目标分类识别方法 | |
CN111160176B (zh) | 基于融合特征的一维卷积神经网络地面雷达目标分类方法 | |
CN111965632B (zh) | 一种基于黎曼流形降维的雷达目标检测方法 | |
CN108447057B (zh) | 基于显著性和深度卷积网络的sar图像变化检测方法 | |
CN111639746A (zh) | 一种基于cnn神经网络的gnss-r海面风速反演方法及*** | |
CN112986950A (zh) | 基于深度学习的单脉冲激光雷达回波特征提取方法 | |
CN111880158A (zh) | 一种基于卷积神经网络序列分类的雷达目标检测方法及*** | |
CN111381216B (zh) | 混合分布的雷达海杂波分析方法及装置 | |
CN115685096B (zh) | 一种基于逻辑回归的二次雷达副瓣抑制方法 | |
CN112597820A (zh) | 一种基于雷达信号分选的目标聚类方法 | |
CN110033043B (zh) | 基于条件生成式对抗网络的雷达一维距离像拒判方法 | |
CN112882010B (zh) | 基于信噪比领域知识网络的高分辨距离像目标识别方法 | |
CN112835008B (zh) | 基于姿态自适应卷积网络的高分辨距离像目标识别方法 | |
CN113887583A (zh) | 基于深度学习的低信噪比下雷达rd图像目标检测方法 | |
CN117451055A (zh) | 一种基于基追踪降噪的水下传感器定位方法和*** | |
CN113723253A (zh) | 基于杂波预分类的运动目标检测方法 | |
CN105738882B (zh) | 对实测杂波协方差矩阵估计性能的白化度评价方法 | |
CN113126086B (zh) | 一种基于状态预测积累的生命探测雷达弱目标检测方法 | |
CN115345216A (zh) | 一种融合先验信息的fmcw雷达干扰消除方法 | |
CN111144316B (zh) | 基于堆栈自编码器的海面目标检测方法 | |
CN117849751B (zh) | 一种针对星凸形不规则形状扩展目标的跟踪方法及*** | |
CN111665500A (zh) | 一种基于单比特压缩感知的脉冲穿墙雷达成像方法 | |
CN114943286B (zh) | 一种基于时域特征与空域特征融合的未知目标判别方法 | |
CN117289262B (zh) | 穿墙雷达目标检测方法及*** | |
CN111537976B (zh) | 一种机载雷达的运动目标检测方法及装置 |
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 |