CN104021289B - 一种非高斯非稳态噪声建模方法 - Google Patents
一种非高斯非稳态噪声建模方法 Download PDFInfo
- Publication number
- CN104021289B CN104021289B CN201410244453.9A CN201410244453A CN104021289B CN 104021289 B CN104021289 B CN 104021289B CN 201410244453 A CN201410244453 A CN 201410244453A CN 104021289 B CN104021289 B CN 104021289B
- Authority
- CN
- China
- Prior art keywords
- noise
- gauss
- gaussian
- sigma
- parameter
- 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
Landscapes
- Complex Calculations (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种非稳态非高斯观测噪声建模方法,用于目标跟踪系统的滤波器状态估计过程中对观测噪声的处理。本发明的主要特征在于:利用高斯混合模型对非稳态非高斯观测噪声进行近似,混合模型中高斯成员的分布参数融合在滤波器的迭代过程中进行计算和更新,能针对观测噪声统计特性的变化在线调节。在滤波器的每一次迭代处理中,高斯成员的分布参数分为先验参数和后验参数。先验参数由上一次迭代获得的参数估计结果进行计算,后验参数通过极大似然估计求取。利用后验分布参数构建高斯成员分布,并组成观测噪声的高斯混合近似模型,该模型能够保证滤波器精度维持在较高水平。本方法对于观测噪声的建模具有较高精度和较强鲁棒性。
Description
技术领域
本发明涉及目标跟踪技术领域,特别涉及一种非高斯非稳态噪声建模方法。
背景技术
目标跟踪是将传感器连续收到的目标运动信号数据对应于各种不确定信息源所产生的不同观测集合进行滤波处理,并相应地估计出机动目标的目标状态参数,如位置、速度和加速度等。由于传感器量测数据中含有大量的干扰成分,必须对量测信息进行处理,因此目标跟踪过程也是一个消除误差的处理过程。目标跟踪已经伴随着雷达、红外、激光、声纳等传感器系统的发展走过了几十年的历程,并由最初的“单目标跟踪”问题逐步发展为多个传感器跟踪多个目标的“多目标跟踪”问题。
贝叶斯滤波为目标跟踪提供了一个利用概率方法递归进行状态估计的框架,利用所有的可用信息,以一定信度确定运动目标处于某一确定状态。贝叶斯滤波通过时间更新和状态更新两个阶段完成目标状态的预测和传感器观测信息的融合,最终实现对目标状态的估计。在实际应用中,贝叶斯滤波框架有多种实现方式,从概率的表示方式上可将实现方法分成离散型和连续型。卡尔曼滤波和粒子滤波分别是连续性和离散型方法的典型代表。在过去的几十年中,研究人员在上述两类方法的基础上进行了深入研究,致力于滤波器估计性能的提高,并发现,如果能准确获取噪声的统计信息,就能够消除或减弱噪声对状态估计产生的干扰,提高估计性能。通常情况下,为了简化问题,时变噪声都被近似为零均值高斯白噪声。然而在很多实际问题中,噪声序列可能是非高斯的,如果近似为高斯噪声将会引入较大的估计误差。针对这一问题,产生了许多非高斯噪声的建模方法。其中一种典型的方法是将非高斯噪声的概率密度函数(probability density function,PDF)近似为混合高斯模型,也就是利用一组高斯概率密度的加权和形式表示非高斯噪声分布。混合高斯模型中高斯成员分布参数的选择决定了非高斯噪声的概率分布,通常的方法是使用期望最大化(expectation maximization,EM)、极大似然估计、非高斯自回归处理或相关函数法等技术确定每个高斯成员的分布参数(主要包括期望、方差、权值),利用这些高斯成员组成的混合模型对非高斯噪声进行近似,并将其用于滤波器的处理过程。
上述建模方法能够较好的描述噪声的非高斯特性,但它们有一个共同特点,这些方法都假设噪声的统计特性是固定不变的,用于近似噪声的高斯混合模型中的每个成员的分布参数一旦确定,在整个滤波过程中也不再改变,也就是将非高斯噪声作为平稳过程来处理。然而在实际问题中,传感器量测值常常受到非高斯非稳态噪声的影响。例如在雷达目标跟踪系统中,雷达量测容易受到闪烁噪声影响,闪烁噪声是一种具有重尾特性的非高斯非稳态噪声,它由目标表面不同部位产生的两个或以上反射间的干涉形成,使反射中心沿目标长度方向发生随机游走,导致波面阵传播过程(propogating)产生形状失真,引起测角误差进而干扰跟踪器的视线(line-of-sight)量测。对于这类非高斯非稳态噪声,以往的高斯混合模型由于使用了确定不变的分布参数,缺乏对噪声统计特性变化的有效跟踪,最终可能导致滤波器的估计性能下降,甚至估计发散。
有关文献:
Ibatoulline E A.Iterative method of maximum likelihood for estimationof non-Gaussian probability density of signals and interferences[C].InProc.IEEE Int.Symp.EMC,2003.
郝红霞,李红.非高斯非平稳背景的红外弱小目标检测[J].激光与红外,2006.
Liu W,Pokharel P P,Principe J C.Correntropy:Properties andapplications in non-Gaussian signal processing[J].IEEE Trans.Signal Process,2007.
公开号为CN103077303A的专利《电子倍增CCD噪声模型的参数估计方法》公开了一种电子倍增CCD噪声模型的参数估计方法,首先对噪声分布模型进行初始化设置,以拟作混合高斯分布模型进行处理;然后对噪声分布模型进行极大似然化迭代,将设定的初始值代入混合模型中,求样本值来自高斯源的后验概率密度,再将潜在数据代入不完全数据的对数函数中计算偏导数,求取极值,即可得到参数的迭代估计值;最后进行判断步,将迭代估计值与初始值比较,依据循环终止条件,判断是否满足终止条件,若满足则迭代停止,若不满足则把迭代值设为初始值,重新进行极大似然化迭代。该发明以简单的步骤实现了对电子倍增CCD图像噪声参数的极大似然估计,有效地降低了极大似然法的复杂度,能对电子倍增CCD图像噪声进行快速、准确的估计。公开号为CN102882820A的专利《认知无线电中非高斯噪声下数字调制信号识别方法》公开了一种在认知无线电中非高斯噪声背景下基于分数低阶循环谱相关系数的数字调制识别方法,解决在认知无线电中非高斯噪声背景下识别调制识别性能差、计算复杂度高的问题。其步骤为:对接收到的信号进行采样;计算采样后信号的分数低阶循环谱计算信号的分数低阶循环谱在循环频率ε=0处的截面和在频率f=0处的截面以及在循环频率ε面的投影和在频率f面的投影的相关系数ρ1、ρ2、ρ3、ρ4和ρ5;设置信号集的判决门限,通过采用基于判决树的分类器将不同调制方式的信号识别出来。在非高斯Alpha稳定分布噪声下,该发明不仅具有较高的识别率和良好的稳健性并且计算复杂度更低,更适合于认知无线电系统。
综上所述,现有的非高斯噪声模型只能应用于非高斯稳态噪声的近似,而无法表示具有非稳态特性的噪声,致使滤波器在处理具有非稳态噪声影响的跟踪问题时产生较大的估计误差。因此如何针对非高斯非稳态噪声设计一种建模方法,是提高滤波器估计性能亟待解决的问题。
发明内容
本发明的目的是,为具有非高斯非稳态特性的观测噪声提供一种具有自适应在线调节能力的建模方法。本发明中提供的方法区别于现有方法的显著特征在于:在将非高斯非稳态噪声建模为高斯混合模型后,各高斯成员的分布参数的取值并不是固定不变的,而是在滤波器每一次迭代过程中针对实际噪声的统计特性对分布参数进行计算和更新。在滤波器的每一次迭代处理中,将噪声模型每个高斯成员的分布参数分为先验参数和后验参数。先验参数在上一次迭代获得的参数估计结果基础上进行初始化,后验参数通过在整个高斯混合模型的观测似然函数取极大值处获得。利用这样获得的分布参数构建高斯分布成员,并组成对于真实噪声的高斯混合近似模型,该模型能够有效针对实际噪声的统计特性变化做出调节,保证滤波器的精度维持在较高水平。
技术方案:
本发明基于高斯混合模型,为具有非高斯非稳态特性的观测噪声提供了一种具有自适应在线调节能力的建模方法。模型建立的实质是对观测噪声的真实概率密度函数进行估计,其特征在于,噪声模型的参数确定是在滤波器估计过程的每一次迭代计算中完成的,以迭代过程中的任意k时刻为例,该方法包含以下步骤:
步骤1、对于任一k时刻的待建模观测噪声,将其近似为高斯混合模型形式,该模型由一系列服从高斯分布的噪声成员构成;
本发明所涉的混合模型中高斯成员的分布参数主要包括:期望、方差、分布权值,这些参数在滤波器的每一次迭代分为两种状态:先验参数和后验参数。先验参数由上一次迭代获得的后验分布参数初始化而来,后验参数根据当前观测噪声的实际分布情况由先验参数更新而来。以下步骤中,步骤2对高斯成员的先验分布参数进行设定,步骤7将先验分布参数更新到后验分布参数。
先验分布参数计算:
步骤2、对于滤波过程中的k时刻,若k=0,则需要进行先验分布参数的初始化,此时滤波器还未开始迭代;若k>0,则先验分布参数由上一次迭代(k-1时刻)获得的后验分布参数计算而来;
步骤3、对当前时刻观测噪声进行采样,形成噪声采样集合;
步骤4、利用先验分布参数计算每个高斯成员关于噪声采样的后验条件概率;
步骤5、建立高斯混合模型的似然函数及其对数形式;
步骤6、通过拉格朗日乘子构建对数似然函数的拉格朗日函数;
后验分布参数计算:
步骤7、通过极大似然估计分别求得混合模型中每个高斯成员的后验分布参数,包括以下步骤:
1)通过令对数似然函数的拉格朗日函数关于高斯成员期望的偏导数为0,求得高斯成员的期望;
2)通过令拉格朗日函数关于高斯成员标准差的偏导数为0,求得高斯成员的方差;
3)通过令拉格朗日函数关于高斯成员分布权值的偏导数为0,求得高斯成员的分布权值;
步骤8、利用每次迭代的后验高斯成员分布参数构建当前观测噪声的高斯混合近似模型,这一模型可用于各种滤波算法中对观测噪声的处理,提高滤波器的估计性能;
步骤9、如果滤波器迭代完成,则由滤波器输出状态估计结果,如果迭代未完成,则返回步骤2。
本发明利用高斯混合模型对非高斯非稳态观测噪声进行建模,模型中高斯成员的分布参数是在滤波器每次迭代过程中根据噪声的实际统计特性进行在线估计和更新的,最终获得的噪声模型用于各种滤波器,能够较为精确的反应实际噪声的分布情况,有效提高滤波器估计性能。
其中主要参数包括:混合高斯模型成员数、观测噪声采样数、拉格朗日乘子、高斯成员先验分布参数、高斯成员后验分布参数以及滤波器迭代次数。混合高斯模型成员数是指用于近似非高斯非稳态观测噪声的混合模型包含的高斯分布的数量,观测噪声采样数是指由观测噪声获得的采样集合包含的采样数量。拉格朗日乘子用于为高斯混合模型的对数似然函数构建拉格朗日函数。高斯成员先验分布参数是依据滤波器上一次迭代结果为噪声混合模型中每个高斯成员分配的先验分布参数,包含先验期望、先验方差、先验分布权值。高斯成员后验分布参数是滤波器每次迭代过程中依据当前观测噪声的实际分布情况对高斯成员先验分布参数进行更新获得的后验分布参数,包含后验期望、后验方差、后验分布权值。滤波器迭代次数用于判断滤波器是否终止迭代。
本发明在滤波器每次迭代过程中计算和更新噪声模型中高斯成员的分布参数,并将分布参数分为先验分布参数和后验分布参数两种状态。这与贝叶斯滤波框架中时间更新和状态更新两种阶段的划分是一致的,因此本发明提出的噪声模型能广泛适用于贝叶斯滤波的多种实现方式。
而且,本发明在利用极大似然估计求取高斯成员后验分布参数时,考虑到计算效率和难度,选取了在同一点处取极值的对数似然函数进行替换。同时通引入拉格朗日乘子,构建拉格朗日函数,在约束条件下通过求极值的方法获得高斯成员后验分布参数。
有益效果:本发明利用一系列加权高斯分布的混合模型对非高斯非稳态观测噪声进行建模,将噪声模型的分布参数估计融合在滤波器的迭代过程中进行,结合滤波器迭代中的时间更新和状态更新,依据实际噪声的统计特性实现噪声分布参数的计算和更新,最终获得的噪声模型能够更接近真实噪声的分布情况,将此噪声模型用于滤波器状态估计,能提高滤波器估计精度,加快收敛速度,是滤波器具有更强的鲁棒性。此外,本发明针对非高斯非稳态噪声进行建模,同时具有广泛的适用范围,可以推广到用于对任意概率分布的噪声进行建模。
附图说明
图1为本发明所述非高斯非稳态方法的系统结构图。
图2为本发明所述非高斯非稳态方法的流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行详细说明。
本发明公开了一种基于混合高斯模型的非高斯非稳态观测噪声建模方法,该方法的具体实施包括模型建立、观测噪声采样、拉格朗日函数构建以及高斯成员分布参数估计等关键内容。本发明所述的非高斯非稳态观测噪声建模方法通过计算机程序实施,图1所示是计算机实现的系统结构图。下面将按照流程详述本发明提出的技术方案的具体实施方式,流程如图2所示。本发明所述噪声模型是对观测噪声概率密度函数的近似,该实施方式主要包含以下关键内容:
对于式(1)描述的目标跟踪系统,为k时刻系统状态变量,服从分布p(xk|xk-1),表示过程噪声。观测变量是关于状态变量xk和观测噪声的函数,服从分布p(zk|xk)。f(·)和h(·)分别称为状态转移函数和观测函数。目标跟踪的任务是,在给定初始状态x0的情况下,通过不断获得和处理观测值zk,对每一时刻的状态值xk进行估计。本发明所述内容即对非高斯非稳态观测噪声vk的概率密度函数进行估计。
滤波器迭代开始前,执行步骤1建立观测噪声模型,并初始化:
步骤1、利用高斯混合模型对k时刻待建模观测噪声进行建模,将观测噪声的概率密度函数近似为一系列服从高斯分布的成员密度函数,如式(2)所示。
其中M为混合模型中高斯成员的数量,通常情况下M=10可以得到较好的建模效果,并保证滤波器估计效率。pi(vk)表示混合模型中第i个高斯成员的概率密度函数,μi,k、和αi,k分别为pi(vk)对应的期望、方差以及权值,αi,k满足式
本发明所涉的混合模型中高斯成员的分布参数在每一次滤波器迭代分为两种状态:先验参数和后验参数。以下步骤中,步骤2对高斯成员的先验分布参数进行设定,步骤7将先验分布参数更新到后验分布参数。
步骤2、若k=0,说明滤波器迭代还未开始,为这些高斯成员初始化先验分布参数:设定每个高斯成员的分布权值为αi,0=1/M,期望为μi,0=0,各高斯成员的方差可依据传感器的先验观测误差特性进行设定,无需考虑实际观测噪声分布情况,造成的偏差可由滤波过程中的参数确定和更新进行修正。若k>0,每次迭代中高斯成员的先验分布参数由上一次迭代获得的后验分布参数计算而来,具体方法如下:
先验权值为:
先验期望为:
先验方差为:
步骤3、对当前时刻噪声进行采样,形成包含N个噪声采样的采样集合,其中第j个噪声采样表示为其中j的取值范围是j=1,…,N。N的取值与噪声复杂程度有关,通常情况下可以取N=10
步骤4、计算每个高斯成员关于噪声采样的后验条件概率,方法如下:
步骤5、将式(3)、(4)、(5)代入式(6)可以得到高斯成员关于噪声采样的后验条件概率。
建立高斯混合模型的对数似然函数如式(7)所示:
步骤6、引入拉格朗日乘子λ,并由约束条件得到拉格朗日函数:
步骤7、通过极大似然估计分别求得混合模型中每个高斯成员的后验分布参数,包括以下步骤:
1)求解可以得到k时刻高斯成员的后验期望
2)求解可以得到高斯成员的后验方差
3)求解并由约束条件得到拉格朗日乘子取值和高斯成员对应的后验权值
λ=M (11)
将式(6)代入分别带入式(9)、(10)和(12),获得对混合模型中每个高斯成员的后验分布参数。
步骤8、将步骤7获得的后验高斯成员分布参数代入式(2),构建当前观测噪声的概率密度函数。这一噪声密度函数可用于各种滤波算法中对观测噪声的处理,提高滤波器的估计性能。
步骤9、判断滤波器的迭代状态,如果滤波器迭代完成,则由滤波器输出状态估计结果,如果迭代未完成,则返回步骤2。
Claims (1)
1.一种非高斯非稳态噪声建模方法,其特征在于,包括以下步骤:
步骤1、利用高斯混合模型对任一k时刻待建模观测噪声进行建模,将观测噪声的概率密度函数pi(vk)近似为一系列服从高斯分布的成员密度函数;
其中M为混合模型中高斯成员的数量,pi(vk)表示混合模型中第i个高斯成员的概率密度函数,μi,k、和αi,k分别为pi(vk)对应的期望、方差以及权值,αi,k满足
步骤2、若k=0,对高斯混合模型中的高斯成员密度函数进行先验分布参数初始化,设定每个高斯成员的分布权值为αi,0=1/M,期望为μi,0=0,设定为传感器的先验观测误差协方差;若k>0,则先验分布参数由k-1时刻迭代获得的后验分布参数计算,方法如下:
先验权值为:
先验期望为:
先验方差为:
步骤3、对当前时刻噪声进行采样,形成包含N个噪声采样的采样集合,其中第j个噪声采样表示为其中j的取值范围是j=1,…,N;
步骤4、计算每个高斯成员关于噪声采样的后验条件概率,方法如下:
步骤5、建立高斯混合模型的对数似然函数;
步骤6、引入拉格朗日乘子λ,并由得到拉格朗日函数:
步骤7、通过极大似然估计分别求得混合模型中每个高斯成员的后验分布参数,包括以下步骤:
1)通过令对数似然函数的拉格朗日函数关于高斯成员期望的偏导数为0,求得高斯成员的后验期望;
2)通过令拉格朗日函数关于高斯成员标准差的偏导数为0,求得高斯成员的后验方差;
3)通过令拉格朗日函数关于高斯成员分布权值的偏导数为0,求得高斯成员的后验权值;
步骤8、将步骤7获得的高斯成员后验权值、期望及方差代入式构建当前观测噪声的概率密度函数;
步骤9、如果滤波器迭代完成,则由滤波器输出状态估计结果,如果迭代未完成,则返回步骤2。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410244453.9A CN104021289B (zh) | 2014-06-04 | 2014-06-04 | 一种非高斯非稳态噪声建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410244453.9A CN104021289B (zh) | 2014-06-04 | 2014-06-04 | 一种非高斯非稳态噪声建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104021289A CN104021289A (zh) | 2014-09-03 |
CN104021289B true CN104021289B (zh) | 2017-04-26 |
Family
ID=51438041
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410244453.9A Expired - Fee Related CN104021289B (zh) | 2014-06-04 | 2014-06-04 | 一种非高斯非稳态噪声建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104021289B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106814608B (zh) * | 2016-12-08 | 2020-03-03 | 浙江中控软件技术有限公司 | 基于后验概率分布的预测控制自适应滤波算法 |
CN109117965B (zh) * | 2017-06-22 | 2022-03-01 | 毫末智行科技有限公司 | 基于卡尔曼滤波器的系统状态预测装置和方法 |
CN109240085B (zh) * | 2018-10-15 | 2021-07-27 | 太原理工大学 | 非高斯系统动态数据校正与系统控制性能优化方法 |
CN110414553B (zh) * | 2019-06-14 | 2021-07-16 | 中国人民解放军海军工程大学 | 一种融合多源信息的备件可靠性评估方法及系统 |
CN110568415B (zh) * | 2019-07-22 | 2022-04-08 | 广东工业大学 | 混合高斯模型下基于Arctan函数的信号检测方法 |
CN113381731B (zh) * | 2020-06-22 | 2023-10-20 | 沈阳工业大学 | 一种针对非高斯噪声的扩散式变步长自适应参数估计方法 |
CN112284518B (zh) * | 2020-09-30 | 2021-12-21 | 南京航空航天大学 | 非高斯多维随机振动波形生成方法 |
CN113625146B (zh) * | 2021-08-16 | 2022-09-30 | 长春理工大学 | 一种半导体器件1/f噪声SαS模型参数估计方法 |
CN114329342B (zh) * | 2021-12-20 | 2024-02-06 | 中交第二公路勘察设计研究院有限公司 | 基于极值-高斯混合分布模型的高速公路通行能力计算方法 |
CN117370737B (zh) * | 2023-12-08 | 2024-02-06 | 成都信息工程大学 | 一种基于自适应高斯滤波器的非稳态非高斯噪声去除方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101256715A (zh) * | 2008-03-05 | 2008-09-03 | 中科院嘉兴中心微系统所分中心 | 无线传感器网络中基于粒子滤波的多车辆声信号分离方法 |
CN102252672A (zh) * | 2011-04-28 | 2011-11-23 | 哈尔滨工程大学 | 一种用于水下导航的非线性滤波方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7668699B2 (en) * | 2005-10-20 | 2010-02-23 | Syracuse University | Optimized stochastic resonance method for signal detection and image processing |
-
2014
- 2014-06-04 CN CN201410244453.9A patent/CN104021289B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101256715A (zh) * | 2008-03-05 | 2008-09-03 | 中科院嘉兴中心微系统所分中心 | 无线传感器网络中基于粒子滤波的多车辆声信号分离方法 |
CN102252672A (zh) * | 2011-04-28 | 2011-11-23 | 哈尔滨工程大学 | 一种用于水下导航的非线性滤波方法 |
Non-Patent Citations (2)
Title |
---|
Weighted LS-SVM for function estimation applied to artifact removal in bio-signal processing;Caicedo A, Van Huffel S.;《Engineering in Medicine and Biology Society》;20101231;988-991 * |
自主式水下航行器同步定位与地图构建算法研究;杜航原;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20120115;63-67 * |
Also Published As
Publication number | Publication date |
---|---|
CN104021289A (zh) | 2014-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104021289B (zh) | 一种非高斯非稳态噪声建模方法 | |
CN105405151B (zh) | 基于粒子滤波和加权Surf的抗遮挡目标跟踪方法 | |
Maggio et al. | Efficient multitarget visual tracking using random finite sets | |
CN106683122A (zh) | 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法 | |
CN105182291B (zh) | 自适应目标新生强度的phd平滑器的多目标跟踪方法 | |
CN110632572A (zh) | 基于无意调相特性的雷达辐射源个体识别方法及装置 | |
CN101567084B (zh) | 基于水平集和分水岭相结合的图像轮廓检测方法 | |
CN103678949B (zh) | 基于密度分析和谱聚类的多扩展目标跟踪量测集划分方法 | |
CN106408591A (zh) | 一种抗遮挡的目标跟踪方法 | |
Fang et al. | Ensemble-based simultaneous input and state estimation for nonlinear dynamic systems with application to wildfire data assimilation | |
CN104732552B (zh) | 基于非平稳条件场的sar图像分割方法 | |
CN104680549B (zh) | 基于高阶邻域tmf模型的sar图像变化检测方法 | |
CN106199580A (zh) | 一种基于模糊推理系统的Singer模型改进算法 | |
CN111985566A (zh) | 一种基于Copula函数的分布式电源时序联合出力典型场景生成方法 | |
CN112213697B (zh) | 一种基于贝叶斯决策理论用于雷达欺骗干扰识别的特征融合方法 | |
CN108871365B (zh) | 一种航向约束下的状态估计方法及系统 | |
CN113780636A (zh) | 一种基于EMD-GRU-Attention的太阳能辐射预测方法 | |
CN104680536A (zh) | 利用改进的非局部均值算法对sar图像变化的检测方法 | |
Liu et al. | EM-based extended object tracking without a priori extension evolution model | |
CN108470189A (zh) | 基于对偶相似度模型的多域辐射源信息融合方法 | |
CN113406574B (zh) | 一种对多功能雷达工作模式序列的在线聚类方法 | |
CN107329131B (zh) | 一种利用粒子滤波的雷达微弱目标检测跟踪方法 | |
CN111262556B (zh) | 一种同时估计未知高斯测量噪声统计量的多目标跟踪方法 | |
Zhao et al. | Remarkable local resampling based on particle filter for visual tracking | |
CN111738339A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170426 Termination date: 20200604 |