CN109212519A - 基于bf-dlstm的窄带雷达目标跟踪方法 - Google Patents
基于bf-dlstm的窄带雷达目标跟踪方法 Download PDFInfo
- Publication number
- CN109212519A CN109212519A CN201810977596.9A CN201810977596A CN109212519A CN 109212519 A CN109212519 A CN 109212519A CN 201810977596 A CN201810977596 A CN 201810977596A CN 109212519 A CN109212519 A CN 109212519A
- Authority
- CN
- China
- Prior art keywords
- narrow
- target
- band radar
- long term
- dlstm
- 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
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/66—Radar-tracking systems; Analogous systems
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于贝叶斯深度长短期记忆网络BF‑DLSTM的窄带雷达目标跟踪方法,主要解决类带雷达目标跟踪中,待跟踪目标运动类型与预先假设的目标运动类型不匹配造成的难以有效跟踪的问题。本发明实现的步骤如下:(1)构建多种运动类型组成的训练数据集;(2)构建贝叶斯深度长短期记忆网络BF‑DLSTM;(3)对训练集进行归一化处理;(4)训练贝叶斯深度长短期记忆网络BF‑DLSTM;(5)窄带雷达目标跟踪。本发明通过基于贝叶斯深度长短期记忆网络BF‑DLSTM的窄带雷达目标跟踪方法,能够对多种运动类型的目标及机动目标实现精度更高地有效跟踪。
Description
技术领域
本发明属于通信技术领域,更进一步涉及窄带雷达通信技术领域中的一种基于贝叶斯深度长短期记忆网络BF-DLSTM(Bayesian Filter based Deep Long Short TermMemory Network)的窄带雷达目标跟踪方法。本发明可用于对窄带雷达观测下的目标数据序列中的单个运动目标进行跟踪,并且可以实现对多种运动类型目标的有效跟踪。
背景技术
窄带雷达观测下的运动目标跟踪的主要任务是,从窄带雷达对目标连续的观测中提取出目标在每一帧回波中的真实状态。随着人们对窄带雷达通信领域的不断深入了解,运动目标跟踪在该领域得到广泛的应用和发展,目前已经存在大量跟踪算法来实现运动目标跟踪。但是由于目标运动类型的不确定性容易造成跟踪滤波器失配,使得准确地进行目标跟踪目标仍面临极大挑战。
电子科技大学在其申请的专利文献“一种用于目标跟踪的强跟踪容积卡尔曼滤波方法”(专利申请号201410654687.0,申请公开号104408744A)中公安开了一种用于目标跟踪的强跟踪容积卡尔曼滤波方法。该方法实现的具体步骤是,(1)建立离散非线性动态***模型;(2)进行***初始化;(3)进行时间更新,引入时变渐消因子λk;(4)进行量测更新;(5)进行滤波更新。该方法的不足之处是,建立完成的离散非线性动态***模型只有一种,在实际应用中若目标运动与其不符合,则容易产生滤波器失配、丢失目标的问题。
东南大学在其申请的专利文献“一种基于交互多模型的目标跟踪方法”(专利申请号201410778057.4,申请公开号104573190A)中公开了一种基于交互式多模型的目标跟踪方法。该方法实现的具体步骤是,(1)根据目标动态条件,设定五组辛格尔Singer模型参数,构建五个Singer模型;(2)根据观测到的目标位置信息,将五个Singer模型进行交互式多模型非线性滤波,估计目标的运动速度、加速度和位置信息;(3)利用目标的运动速度和加速度,计算目标转弯角速率;(4)将目标转弯角速率与设定的阈值进行比较,判断是否发生转弯运动,如果未发生转弯运动,则将步骤(2)中得到的位置信息作为目标跟踪结果;(5)如果发生转弯运动,选取三个相邻时刻的转弯角速率构建三个转弯模型来近似目标的实际运动,进行交互式多模型非线性滤波得到目标位置作为目标跟踪结果。该方法存在的不足之处在于,观测误差会影响目标转弯角速率计算的准确性,进而影响关于目标是否发生转弯运动的判断,并且在目标运动特性发生变化时,用于近似目标实际运动的多个模型近似误差会增大,使得跟踪滤波误差增大。
发明内容
本发明的目的是针对上述现有技术的不足,提出一种基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,以实现更加准确地对多种不同运动类型的目标的跟踪。
为实现本发明目的思路是,首先,根据待跟踪目标的运动类型构建训练数据集,利用窄带雷达***的观测范围归一化训练数据集。其次,构建出贝叶斯深度长短期记忆网络BF-DLSTM,并用归一化后的训练数据集来训练贝叶斯深度长短期记忆网络BF-DLSTM。最后,利用窄带雷达***的观测范围归一化窄带雷达对待跟踪的目标的观测,输入到训练好的贝叶斯深度长短期记忆网络BF-DLSTM,再分别利用窄带雷达的观测范围将贝叶斯深度长短期记忆网络BF-DLSTM关于滤波值和预测值的输出逆归一化,得到待跟踪目标滤波和预测的状态,完成目标跟踪。
本发明的具体步骤如下:
(1)构建多种运动类型组成的训练数据集:
(1a)利用距离单元公式,计算窄带雷达观测场景中距离单元的总数Nr;
(1b)利用目标运动轨迹的构建方法,构建10Nr条训练目标的运动轨迹,将每一个训练目标每次被窄带雷达扫描时的坐标值,组成二维向量,作为当前时刻该训练目标的真实运动状态;
(1c)利用窄带雷达观测方程,生成对应于每一个训练目标的真实运动状态的窄带雷达观测值;
(1d)将所有训练目标的真实运动状态与相应的窄带雷达观测值组成训练数据集;
(2)构建贝叶斯深度长短期记忆网络BF-DLSTM:
(2a)搭建一个10层的贝叶斯深度长短期记忆网络BF-DLSTM,其结构依次为:观测输入编码层→含有3个隐藏层的滤波长短期记忆网络→滤波输出解码层→滤波输入编码层→含有3个隐藏层的预测长短期记忆网络→预测输出解码层;
(2b)设置每一个网络各层参数如下:
将观测输入编码层设置为2×16的全连接层;
将滤波长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256;
将滤波输出解码层设置为256×2的全连接层;
将滤波输入编码层设置为2×16的全连接层;
将预测长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256;
将预测输出解码层设置为256×2的全连接层;
(3)对训练集进行归一化处理:
(3a)用训练目标的真实运动状态除以窄带雷达观测的最远距离,得到归一化后的目标真实运动状态;
(3b)用窄带雷达观测训练目标的径向距离除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测训练目标的径向距离;用窄带雷达观测训练目标的方位角除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测训练目标的方位角;
(4)训练贝叶斯深度长短期记忆网络BF-DLSTM:
(4a)将归一化后的训练集中径向距离观测值与方位角观测值,分别输入到贝叶斯深度长短期记忆网络BF-DLSTM中得到网络的实际输出值,将归一化后的训练集中实际运动状态作为贝叶斯深度长短期记忆网络BF-DLSTM的期望输出值;
(4b)将深度长短期记忆网络的实际输出值与期望输出值,组成贝叶斯深度长短期记忆网络BF-DLSTM的代价函数;
(4c)使用最小批梯度下降法,训练贝叶斯深度长短期记忆网络BF-DLSTM,调整梯度下降法中梯度更新的步长,直到贝叶斯深度长短期记忆网络BF-DLSTM的代价函数趋于不变,得到训练好的贝叶斯深度长短期记忆网络BF-DLSTM;
(5)窄带雷达目标跟踪:
(5a)用窄带雷达对待跟踪目标径向距离的观测值,除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测待跟踪目标的径向距离;用窄带雷达对待跟踪目标方位角的观测值,除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测待跟踪目标的方位角;
(5b)将归一化后待跟踪目标位置的径向距离和方位角,输入训练好的贝叶斯深度长短期记忆网络BF-DLSTM,得到相应的归一化后待跟踪目标位置的滤波值和预测值;
(5c)用窄带雷达观测的最远距离,分别乘以归一化后待跟踪目标位置的滤波值和预测值,得到待跟踪目标位置的滤波值和预测值,完成窄带雷达目标跟踪。
发明与现有的技术相比具有以下优点:
第一,由于本发明构建多种运动类型组成的训练数据集,克服了现有技术中只能跟踪预先假设的单个运动类型,目标实际运动类型与其不同则容易发生滤波器失配、丢失目标的问题,使得本发明能够更有效地适用于不同运动类型目标的跟踪。
第二,由于本发明构建贝叶斯深度长短期记忆网络BF-DLSTM,克服了现有技术中采用多模型方法跟踪的目标运动特性发生变化时,用于近似目标实际运动的多个模型近似误差会增大,使得跟踪滤波误差增大的问题,使得本发明能够适用于机动目标的跟踪,提高了滤波精度。
附图说明
图1为本发明的流程图;
图2为本发明贝叶斯深度长短期记忆网络BF-DLSTM的结构示意图;
图3为本发明的仿真图。
具体实施方式
下面结合附图对本发明做进一步的描述。
参照图1,对本发明的具体步骤做进一步的描述。
步骤1,构建多种运动类型组成的训练数据集。
利用距离单元公式,计算窄带雷达观测场景中距离单元的总数Nr。
所述距离单元公式如下:
其中,Nr表示窄带雷达观测场景中距离单元的总数,c表示光速,B表示窄带雷达发射信号的带宽,Δρ表示窄带雷达观测的距离范围,Δθ表示窄带雷达观测的方位角范围,θ3dB表示窄带雷达发射信号的3dB波瓣宽度。
利用目标运动轨迹的构建方法,构建10Nr条训练目标的运动轨迹,将每一个训练目标每次被窄带雷达扫描时的坐标值,组成二维向量,作为当前时刻该训练目标的真实运动状态。
所述目标运动轨迹的构建方法的步骤如下:
第1步,在待跟踪的目标运动类型中随机选择训练目标的运动类型。
第2步,在窄带雷达扫描范围内,随机设定训练目标运动的初始位置。
第3步,根据待跟踪目标的机动特性,设定训练目标运动的过程噪声的方差。
第4步,若训练目标运动为匀速运动,在待跟踪目标运动速度范围内随机选择其运动速度,若训练目标运动为匀变速运动,在待跟踪目标运动加速度范围内随机选择其加速度,若训练目标运动为协同转弯运动,在待跟踪目标运动转弯率范围内随机选择其转弯率。
第5步,根据不同运动类型的状态转移方程,仿真训练目标的真实运动状态。
利用窄带雷达观测方程,生成对应于每一个训练目标的真实运动状态的窄带雷达观测值。
所述窄带雷达观测方程如下:
ρm=ρt+ερ
θm=θt+εθ
其中,ρm表示窄带雷达观测训练目标的径向距离,ρt表示训练目标真实的径向距离,ερ表示窄带雷达观测训练目标的径向距离的观测噪声,θm表示窄带雷达观测训练目标的方位角,θt表示训练目标真实的方位角,εθ表示窄带雷达观测训练目标方位角的观测噪声。
将所有训练目标的真实运动状态与相应的窄带雷达观测值组成训练数据集。
步骤2,构建贝叶斯深度长短期记忆网络BF-DLSTM。
搭建一个10层的贝叶斯深度长短期记忆网络BF-DLSTM,其结构依次为:观测输入编码层→含有3个隐藏层的滤波长短期记忆网络→滤波输出解码层→滤波输入编码层→含有3个隐藏层的预测长短期记忆网络→预测输出解码层。
设置每一个网络各层参数如下:
将观测输入编码层设置为2×16的全连接层。
将滤波长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256。
将滤波输出解码层设置为256×2的全连接层。
将滤波输入编码层设置为2×16的全连接层。
将预测长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256。
将预测输出解码层设置为256×2的全连接层。
参照图2,对贝叶斯深度长短期记忆网络BF-DLSTM的结构做进一步描述。
图2中的hk-1和ck-1分别表示k-1时刻滤波长短期记忆网络的隐状态和单元状态,zk表示第k次窄带雷达对目标的观测,经过滤波长短期记忆网络的隐状态和单元状态再次输入到预测长短期记忆网络中,表示输入训练集中窄带雷达对目标在第k次的观测值时贝叶斯深度长短期记忆网络BF-DLSTM关于滤波值的输出,表示输入训练集中窄带雷达对目标在第k次的观测值时贝叶斯深度长短期记忆网络BF-DLSTM关于预测值的输出。
步骤3,对训练集进行归一化处理。
用训练目标的真实运动状态除以窄带雷达观测的最远距离,得到归一化后的目标真实运动状态。
用窄带雷达观测训练目标的径向距离除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测训练目标的径向距离;用窄带雷达观测训练目标的方位角除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测训练目标的方位角。
步骤4,训练贝叶斯深度长短期记忆网络BF-DLSTM。
将归一化后的训练集中径向距离观测值与方位角观测值,分别输入到贝叶斯深度长短期记忆网络BF-DLSTM中得到网络的实际输出值,将归一化后的训练集中实际运动状态作为贝叶斯深度长短期记忆网络BF-DLSTM的期望输出值。
将贝叶斯深度长短期记忆网络BF-DLSTM的实际输出值与期望输出值,组成贝叶斯深度长短期记忆网络BF-DLSTM的代价函数。
所述贝叶斯深度长短期记忆网络BF-DLSTM的代价函数如下:
其中,J(θ)表示贝叶斯深度长短期记忆网络BF-DLSTM的代价函数,θ表示贝叶斯深度长短期记忆网络BF-DLSTM中连接每一层网络的参数,其在训练贝叶斯深度长短期记忆网络BF-DLSTM的过程中随着代价函数趋于不变而趋于最优,
n表示训练集中目标的序号,k表示每一个目标被窄带雷达观测的序号,N表示训练集中目标的总数,Σ表示求和操作,Ln表示训练集中第n个目标被窄带雷达观测的次数,表示训练集中第n个目标在第k次被窄带雷达观测时的真实状态,表示输入训练集中窄带雷达对第n个目标在第k次的观测值时滤波输出解码层的输出,∈表示贝叶斯深度长短期记忆网络BF-DLSTM预测精度相对于滤波精度在代价函数中所占的权重,表示训练集中第n个目标在第k+1次被窄带雷达观测时的真实状态,表示输入训练集中窄带雷达对第n个目标在第k次的观测值时预测输出解码层的输出。
使用最小批梯度下降法,训练贝叶斯深度长短期记忆网络BF-DLSTM,调整梯度下降法中梯度更新的步长,直到贝叶斯深度长短期记忆网络BF-DLSTM的代价函数趋于不变,得到训练好的贝叶斯深度长短期记忆网络BF-DLSTM。
步骤5,窄带雷达目标跟踪。
用窄带雷达对待跟踪目标径向距离的观测值,除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测待跟踪目标的径向距离;用窄带雷达对待跟踪目标方位角的观测值,除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测待跟踪目标的方位角。
将归一化后待跟踪目标位置的径向距离和方位角,输入训练好的贝叶斯深度长短期记忆网络BF-DLSTM,得到相应的归一化后待跟踪目标位置的滤波值和预测值。
用窄带雷达观测的最远距离,分别乘以归一化后待跟踪目标位置的滤波值和预测值,得到待跟踪目标位置的滤波值和预测值,完成窄带雷达目标跟踪。
下面结合仿真实验对本发明的效果做进一步说明。
1.仿真实验条件:
本发明仿真实验的硬件测试平台是:处理器为CPU Xeon E5-2643,主频为3.4GHz,内存64GB;软件平台为:Ubuntu 16.04LTS,64位操作***,Python 2.7。
2.仿真内容及仿真结果分析:
为证明贝叶斯深度长短期记忆网络BF-DLSTM方法可以同时实现对不同运动类型目标以及机动目标的跟踪,且能够获得更高的滤波精度,本发明的仿真实验采用本发明的方法贝叶斯深度长短期记忆网络BF-DLSTM和现有技术的交互多模型方法,分别针对做匀速运动和机动的目标进行了跟踪的仿真实验。
本发明的仿真实验训练目标为民航飞机,其所做的运动主要为匀速运动、匀加速运动及协同转弯运动,窄带雷达的采样时间为10s,径向距离观测误差标准差为30m,方位角观测误差标准差为0.5°,窄带雷达观测径向距离范围为15km~150km,方位角范围为0~90°。在构建多种运动类型组成的训练数据集时,按照初速度为100m/s~300m/s,最大加速度30m/s2,最大转弯率4.5°/s产生训练目标的真实运动状态与窄带雷达观测。
在窄带雷达目标跟踪阶段,窄带雷达对匀速运动的目标和机动目标分别进行持续50步地跟踪。其中,匀速运动目标的起始位置为[30km,80km],运动速度为机动目标的起始位置为[30km,10k]m,初始速度为[0,150m/s]。在100s的匀速运动之后,目标按照[1.5m/s2,-1.5m/s2]的加速度做持续100s的匀加速运动,目标继续按照[150m/s,]0的速度匀速运动100s,然后按照1°/s的速度朝正北方向转弯,持续90s之后,最后以[0,150m/s]的速度继续做110s的匀速直线运动。
针对上述测试场景,分别使用传统的包含匀速运动类型和匀加速运动类型的交互多模型方法(Interactive Multiple Model,IMM)和本发明方法,通过5000次蒙特卡洛结果,比较滤波后的结果。
图3(a)为本发明仿真实验中使用贝叶斯深度长短期记忆网络BF-DLSTM与交互多模型方法跟踪一个匀速运动目标时,滤波位置的均方误差的仿真结果对比图。图3(a)中x轴表示跟踪的步数,y轴表示两个方法关于目标位置滤波的均方误差。其中以实线标示的为采用现有技术交互多模型方法进行跟踪的滤波位置均方误差曲线,点划线标示的为采用本发明方法贝叶斯深度长短期记忆网络BF-DLSTM进行跟踪的滤波位置均方误差曲线。图3(a)中整个跟踪阶段本发明方法贝叶斯深度长短期记忆网络BF-DLSTM均可以获得比现有技术交互多模型方法更低的滤波位置均方误差。可见针对匀速运动目标,本发明方法贝叶斯深度长短期记忆网络BF-DLSTM可以获得更高的滤波精度。
图3(b)为本发明仿真实验中使用贝叶斯深度长短期记忆网络BF-DLSTM与交互多模型方法跟踪一个机动目标时,滤波位置的均方误差的仿真结果对比图。图3(b)中x轴表示跟踪的步数,y轴表示两个方法关于目标位置滤波的均方误差。其中以实线标示的为采用现有技术交互多模型方法进行跟踪的滤波位置均方误差曲线,以点划线标示的为采用本发明方法贝叶斯深度长短期记忆网络BF-DLSTM进行跟踪的滤波位置均方误差曲线。图3(b)中整个跟踪阶段本发明方法贝叶斯深度长短期记忆网络BF-DLSTM均可以获得比现有技术交互多模型方法更低的滤波位置均方误差,且本发明方法的滤波位置均方误差曲线起伏更小。可见针对做机动的目标,本发明方法在整个跟踪阶段滤波都可以获得更高的滤波精度,且本发明方法对目标的机动不敏感,在目标运动类型发生变化时仍可以保持有效的跟踪。
综合图3(a)和图3(b),本发明方法贝叶斯深度长短期记忆网络BF-DLSTM针对匀速运动目标和机动目标均能够实现相比于现有技术交互多模型方法更有效的跟踪,可以看出:本发明方法,能够适用于不同运动类型目标的跟踪,并且利用本发明方法贝叶斯深度长短期记忆网络BF-DLSTM能够获得的滤波精度更高,可以实现更有效地跟踪。
Claims (5)
1.一种基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,其特征在于,构建多种运动类型组成的训练数据集,构建贝叶斯深度长短期记忆网络BF-DLSTM,该方法的具体步骤包括如下:
(1)构建多种运动类型组成的训练数据集:
(1a)利用距离单元公式,计算窄带雷达观测场景中距离单元的总数Nr;
(1b)利用目标运动轨迹的构建方法,构建10Nr条训练目标的运动轨迹,将每一个训练目标每次被窄带雷达扫描时的坐标值,组成二维向量,作为当前时刻该训练目标的真实运动状态;
(1c)利用窄带雷达观测方程,生成对应于每一个训练目标的真实运动状态的窄带雷达观测值;
(1d)将所有训练目标的真实运动状态与相应的窄带雷达观测值组成训练数据集;
(2)构建贝叶斯深度长短期记忆网络BF-DLSTM:
(2a)搭建一个10层的贝叶斯深度长短期记忆网络BF-DLSTM,其结构依次为:观测输入编码层→含有3个隐藏层的滤波长短期记忆网络→滤波输出解码层→滤波输入编码层→含有3个隐藏层的预测长短期记忆网络→预测输出解码层;
(2b)设置每一个网络各层参数如下:
将观测输入编码层设置为2×16的全连接层;
将滤波长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256;
将滤波输出解码层设置为256×2的全连接层;
将滤波输入编码层设置为2×16的全连接层;
将预测长短期记忆网络的隐藏层数设置为3,输入单元的个数设置为16,隐单元的个数设置为256;
将预测输出解码层设置为256×2的全连接层;
(3)对训练集进行归一化处理:
(3a)用训练目标的真实运动状态除以窄带雷达观测的最远距离,得到归一化后的目标真实运动状态;
(3b)用窄带雷达观测训练目标的径向距离除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测训练目标的径向距离;用窄带雷达观测训练目标的方位角除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测训练目标的方位角;
(4)训练贝叶斯深度长短期记忆网络BF-DLSTM:
(4a)将归一化后的训练集中径向距离观测值与方位角观测值,分别输入到贝叶斯深度长短期记忆网络BF-DLSTM中得到网络的实际输出值,将归一化后的训练集中实际运动状态作为贝叶斯深度长短期记忆网络BF-DLSTM的期望输出值;
(4b)将深度长短期记忆网络的实际输出值与期望输出值,组成贝叶斯深度长短期记忆网络BF-DLSTM的代价函数;
(4c)使用最小批梯度下降法,训练贝叶斯深度长短期记忆网络BF-DLSTM,调整梯度下降法中梯度更新的步长,直到贝叶斯深度长短期记忆网络BF-DLSTM的代价函数趋于不变,得到训练好的贝叶斯深度长短期记忆网络BF-DLSTM;
(5)窄带雷达目标跟踪:
(5a)用窄带雷达对待跟踪目标径向距离的观测值,除以窄带雷达观测的最远距离,得到归一化后的窄带雷达观测待跟踪目标的径向距离;用窄带雷达对待跟踪目标方位角的观测值,除以窄带雷达观测的最大方位角,得到归一化后的窄带雷达观测待跟踪目标的方位角;
(5b)将归一化后待跟踪目标位置的径向距离和方位角,输入训练好的贝叶斯深度长短期记忆网络BF-DLSTM,得到相应的归一化后待跟踪目标位置的滤波值和预测值;
(5c)用窄带雷达观测的最远距离,分别乘以归一化后待跟踪目标位置的滤波值和预测值,得到待跟踪目标位置的滤波值和预测值,完成窄带雷达目标跟踪。
2.根据权利要求1所述的基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,其特征在于:步骤(1a)中所述的距离单元公式如下:
其中,Nr表示窄带雷达观测场景中距离单元的总数,c表示光速,B表示窄带雷达发射信号的带宽,Δρ表示窄带雷达观测的距离范围,Δθ表示窄带雷达观测的方位角范围,θ3dB表示窄带雷达发射信号的3dB波瓣宽度。
3.根据权利要求1所述的基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,其特征在于,步骤(1b)中所述的目标运动轨迹的构建方法的步骤如下:
第一步,在待跟踪的目标运动类型中随机选择训练目标的运动类型;
第二步,在窄带雷达扫描范围内,随机设定训练目标运动的初始位置;
第三步,根据待跟踪目标的机动特性,设定训练目标运动的过程噪声的方差;
第四步,若训练目标运动为匀速运动,在待跟踪目标运动速度范围内随机选择其运动速度,若训练目标运动为匀变速运动,在待跟踪目标运动加速度范围内随机选择其加速度,若训练目标运动为协同转弯运动,在待跟踪目标运动转弯率范围内随机选择其转弯率;
第五步,根据不同运动类型的状态转移方程,仿真训练目标的真实运动状态。
4.根据权利要求1所述的基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,其特征在于:步骤(1c)中所述的窄带雷达观测方程如下:
ρm=ρt+ερ
θm=θt+εθ
其中,ρm表示窄带雷达观测训练目标的径向距离,ρt表示训练目标真实的径向距离,ερ表示窄带雷达观测训练目标的径向距离的观测噪声,θm表示窄带雷达观测训练目标的方位角,θt表示训练目标真实的方位角,εθ表示窄带雷达观测训练目标方位角的观测噪声。
5.根据权利要求1所述的基于贝叶斯深度长短期记忆网络BF-DLSTM的窄带雷达目标跟踪方法,其特征在于:步骤(4b)中所述的贝叶斯深度长短期记忆网络BF-DLSTM的代价函数如下:
其中,J(θ)表示贝叶斯深度长短期记忆网络BF-DLSTM的代价函数,θ表示贝叶斯深度长短期记忆网络BF-DLSTM中连接每一层网络的参数,其在训练贝叶斯深度长短期记忆网络BF-DLSTM的过程中随着代价函数趋于不变而趋于最优,n表示训练集中目标的序号,k表示每一个目标被窄带雷达观测的序号,N表示训练集中目标的总数,Σ表示求和操作,Ln表示训练集中第n个目标被窄带雷达观测的次数,表示训练集中第n个目标在第k次被窄带雷达观测时的真实状态,表示输入训练集中窄带雷达对第n个目标在第k次的观测值时滤波输出解码层的输出,∈表示贝叶斯深度长短期记忆网络BF-DLSTM预测精度相对于滤波精度在代价函数中所占的权重,表示训练集中第n个目标在第k+1次被窄带雷达观测时的真实状态,表示输入训练集中窄带雷达对第n个目标在第k次的观测值时预测输出解码层的输出。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810977596.9A CN109212519B (zh) | 2018-08-27 | 2018-08-27 | 基于bf-dlstm的窄带雷达目标跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810977596.9A CN109212519B (zh) | 2018-08-27 | 2018-08-27 | 基于bf-dlstm的窄带雷达目标跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109212519A true CN109212519A (zh) | 2019-01-15 |
CN109212519B CN109212519B (zh) | 2023-04-07 |
Family
ID=64989222
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810977596.9A Active CN109212519B (zh) | 2018-08-27 | 2018-08-27 | 基于bf-dlstm的窄带雷达目标跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109212519B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110212891A (zh) * | 2019-06-02 | 2019-09-06 | 西安电子科技大学 | 基于nrp-lstm的非线性滤波方法 |
CN110780290A (zh) * | 2019-11-01 | 2020-02-11 | 西安电子科技大学 | 基于lstm网络的多机动目标跟踪方法 |
CN113534128A (zh) * | 2020-10-21 | 2021-10-22 | 中国人民解放军空军预警学院 | 一种机载预警雷达海面机动舰船目标自适应跟踪方法 |
CN114266925A (zh) * | 2021-12-30 | 2022-04-01 | 华北电力大学 | 一种基于dlstm-rf的用户窃电检测方法及*** |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09113610A (ja) * | 1995-10-23 | 1997-05-02 | Japan Radio Co Ltd | 目標追尾レーダ装置 |
CN102254176A (zh) * | 2011-04-11 | 2011-11-23 | 西安电子科技大学 | 基于状态空间模型的雷达高分辨距离像目标识别方法 |
CN102608586A (zh) * | 2012-02-11 | 2012-07-25 | 中国人民解放军陆军军官学院 | 一种基于神经网络的合成孔径雷达射频干扰抑制方法 |
CN104020466A (zh) * | 2014-06-17 | 2014-09-03 | 西安电子科技大学 | 基于变结构多模型的机动目标跟踪方法 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
US9297886B1 (en) * | 2013-03-12 | 2016-03-29 | Lockheed Martin Corporation | Space time adaptive technique for suppression of spaceborne clutter |
CN105717505A (zh) * | 2016-02-17 | 2016-06-29 | 国家电网公司 | 利用传感网进行多目标跟踪的数据关联方法 |
CN106257302A (zh) * | 2016-07-01 | 2016-12-28 | 电子科技大学 | 一种杂波环境下的认知雷达目标跟踪方法 |
CN108226889A (zh) * | 2018-01-19 | 2018-06-29 | 中国人民解放军陆军装甲兵学院 | 一种雷达目标识别的分类器模型训练方法 |
JP2018115974A (ja) * | 2017-01-19 | 2018-07-26 | 沖電気工業株式会社 | 情報処理装置、情報処理方法、およびプログラム |
-
2018
- 2018-08-27 CN CN201810977596.9A patent/CN109212519B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09113610A (ja) * | 1995-10-23 | 1997-05-02 | Japan Radio Co Ltd | 目標追尾レーダ装置 |
CN102254176A (zh) * | 2011-04-11 | 2011-11-23 | 西安电子科技大学 | 基于状态空间模型的雷达高分辨距离像目标识别方法 |
CN102608586A (zh) * | 2012-02-11 | 2012-07-25 | 中国人民解放军陆军军官学院 | 一种基于神经网络的合成孔径雷达射频干扰抑制方法 |
US9297886B1 (en) * | 2013-03-12 | 2016-03-29 | Lockheed Martin Corporation | Space time adaptive technique for suppression of spaceborne clutter |
CN104020466A (zh) * | 2014-06-17 | 2014-09-03 | 西安电子科技大学 | 基于变结构多模型的机动目标跟踪方法 |
CN104408744A (zh) * | 2014-11-17 | 2015-03-11 | 电子科技大学 | 一种用于目标跟踪的强跟踪容积卡尔曼滤波方法 |
CN105717505A (zh) * | 2016-02-17 | 2016-06-29 | 国家电网公司 | 利用传感网进行多目标跟踪的数据关联方法 |
CN106257302A (zh) * | 2016-07-01 | 2016-12-28 | 电子科技大学 | 一种杂波环境下的认知雷达目标跟踪方法 |
JP2018115974A (ja) * | 2017-01-19 | 2018-07-26 | 沖電気工業株式会社 | 情報処理装置、情報処理方法、およびプログラム |
CN108226889A (zh) * | 2018-01-19 | 2018-06-29 | 中国人民解放军陆军装甲兵学院 | 一种雷达目标识别的分类器模型训练方法 |
Non-Patent Citations (7)
Title |
---|
HAMID IZADINIA: "Multi-pose multi-target tracking for activity understanding", 《2013 IEEE WORKSHOP ON APPLICATIONS OF COMPUTER VISION (WACV)》 * |
Y. BAR-SHALOM: "Tracking a maneuvering target using input estimation versus the interacting multiple model algorithm", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
YUN ZHANG: "A novel multi-target track initiation method based on convolution neural network", 《2017 INTERNATIONAL WORKSHOP ON REMOTE SENSING WITH INTELLIGENT PROCESSING (RSIP)》 * |
刘婧逸等: "基于JEM的雷达窄带目标识别方法研究", 《信息化研究》 * |
夏双志等: "基于贝叶斯理论的一种接力跟踪方法", 《电子与信息学报》 * |
祝依龙等: "雷达目标特征信号与运动模式的相关性研究", 《中国科学:信息科学》 * |
韦文等: "基于递推贝叶斯的检测跟踪一体化方法", 《宇航学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110212891A (zh) * | 2019-06-02 | 2019-09-06 | 西安电子科技大学 | 基于nrp-lstm的非线性滤波方法 |
CN110780290A (zh) * | 2019-11-01 | 2020-02-11 | 西安电子科技大学 | 基于lstm网络的多机动目标跟踪方法 |
CN110780290B (zh) * | 2019-11-01 | 2022-10-21 | 西安电子科技大学 | 基于lstm网络的多机动目标跟踪方法 |
CN113534128A (zh) * | 2020-10-21 | 2021-10-22 | 中国人民解放军空军预警学院 | 一种机载预警雷达海面机动舰船目标自适应跟踪方法 |
CN114266925A (zh) * | 2021-12-30 | 2022-04-01 | 华北电力大学 | 一种基于dlstm-rf的用户窃电检测方法及*** |
CN114266925B (zh) * | 2021-12-30 | 2022-09-30 | 华北电力大学 | 一种基于dlstm-rf的用户窃电检测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109212519B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109212519A (zh) | 基于bf-dlstm的窄带雷达目标跟踪方法 | |
CN107396322A (zh) | 基于路径匹配与编码译码循环神经网络的室内定位方法 | |
CN106372646B (zh) | 基于srck-gmcphd滤波的多目标跟踪方法 | |
CN104237879B (zh) | 一种雷达***中的多目标跟踪方法 | |
CN105043380A (zh) | 基于微机电传感器、WiFi定位、磁场匹配的室内导航方法 | |
CN106646395B (zh) | 一种飞行目标的雷达回波推演方法 | |
Tepecik et al. | A novel hybrid model for inversion problem of atmospheric refractivity estimation | |
Yang et al. | UGHF for acoustic tracking with state-dependent propagation delay | |
CN109188420A (zh) | 基于深度长短期记忆网络的窄带雷达目标跟踪方法 | |
Potapov | Fractal and topological sustainable methods of overcoming expected uncertainty in the radiolocation of low-contrast targets and in the processing of weak multi-dimensional signals on the background of high-intensity noise: A new direction in the statistical decision theory | |
CN103954964A (zh) | 多角度合成孔径雷达数据获取的方法 | |
CN109655827A (zh) | 一种针对单建筑墙角的隐蔽目标定位方法 | |
CN108549062A (zh) | 一种用于海面搜索雷达的***平台及多模型目标跟踪方法 | |
CN109507654A (zh) | 一种基于ls的复杂环境下相位信息计算方法 | |
CN102830394A (zh) | 基于多谱线积累的弱目标探测方法 | |
Lai et al. | Bearings-only tracking and Doppler-bearing tracking with inequality constraint | |
CN114488104B (zh) | 基于交互一致性的天波超视距雷达目标跟踪方法 | |
CN106546982B (zh) | 一种基于模拟退火算法的层状介质的波形反演方法 | |
Garapati Vaishnavi et al. | Underwater bearings-only tracking using particle filter | |
CN114660587A (zh) | 基于Jerk模型的跳跃滑翔弹道目标跟踪方法及*** | |
Potapov | Physical bases and principles of construction of fractal radars and fractal sensors: a new direction-fractal analysis and its application in the theory of statistical solutions and in statistical radio engineering | |
Ermakov et al. | The use of mathematical modeling to solve the problem of detecting a low-flying target | |
KR20150055260A (ko) | 분산형 수색 레이더 환경에서 획득한 데이터 간의 시간 동기화 방법 | |
CN110212891A (zh) | 基于nrp-lstm的非线性滤波方法 | |
CN113514809B (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 |