广义帕累托分布参数的迭代最大似然估计方法技术

技术编号:15327711 阅读:623 留言:0更新日期:2017-05-16 11:46
本发明专利技术公开了一种广义帕累托分布参数的迭代最大似然估计方法,主要解决现有参数估计方法的估计精度差和执行效率低的问题。其技术方案是:1通过海面监视雷达获取海杂波数据样本;2将获取到的杂波数据样本按杂波功率归一化;3迭代计算数据归一化后的杂波数据样本的似然估计量;4利用似然估计量计算广义帕累托分布参数。本发明专利技术降低了广义帕累托分布最大似然估计方法的时间复杂度,且估计精度高,计算速度快,能够适应雷达系统信号实时处理的要求,可以用于海杂波背景下的目标检测。

Iterative maximum likelihood estimation method for parameters of generalized Pareto distribution

The invention discloses an iterative maximum likelihood estimation method for generalized Pareto distribution parameters, mainly solving the problems of poor estimation accuracy and low execution efficiency of the existing parameter estimation methods. The technical scheme is: 1 by sea surveillance radar sea clutter data acquisition; clutter data sample 2 will get to the clutter power according to the normalized likelihood calculation; clutter data sample data normalized 3 iterative estimators; calculation of generalized Pareto distribution parameters using 4 likelihood estimator. The invention reduces the generalized Pareto distribution of maximum likelihood estimation method of time complexity and high estimation accuracy, fast calculation speed and can adapt to the real-time radar signal processing system, can be used for target detection in sea clutter background.

【技术实现步骤摘要】
广义帕累托分布参数的迭代最大似然估计方法
本专利技术属于信号处理
,具体涉及一种广义帕累托分布参数估计方法,可用于海杂波背景下的目标检测。
技术介绍
海杂波背景下的目标检测技术是雷达应用技术中一个至关重要的研究方向,在军事和民用领域已经得到广泛应用。而对于海杂波统计特性的准确分析是海杂波背景下目标检测技术能否取得良好效果的重要因素。因此,给出合适的模型并对于其模型参数进行准确估计成为我们需要解决的重要问题。随着现代雷达系统距离分辨力的提高,雷达回波出现以往低分辨力雷达系统所没有的统计特性,通常表现为其回波包络的拖尾变长,异常值变多的特点。而广义帕累托分布作为复合高斯模型的一种,在对于高分辨低擦地角海杂波的功率分布拟合上取得了很好的效果。因此在海杂波统计特性的研究中占据重要地位。而在实验数据下给出概率分布模型的模型参数又作为模型研究的关键部分,因此在重拖尾的杂波数据下给出广义怕累托的分布参数具有重要的研究意义。近年来,很多研究者对广义帕累托分布的参数估计方法,提出了一些基于特定条件下的广义帕累托分布参数估计理论。文献“Castillo,E.,Hadi,A.S.,1997.FittingthegeneralizedParetodistributiontodata.J.Amer.Statist.Assoc.92,1609–1620.”中给出广义帕累托分布的矩估计以及最大似然估计方法,分别根据样本矩以及似然函数对于参数进行估计,但是由于矩估计本身容易受到样本数量和异常数据的影响,其估计精度难以保证。而最大似然估计的估计精度虽然能够满足要求,但是算法时间复杂度高,因此工程实现较为困难。文献“Arnold,B.C.,Press,S.J.,1989.BayesianestimationandpredictionforParetodata.J.Amer.Statist.Assoc.84,1079–1084.”给出了基于先验信息的广义帕累托分布参数估计方法,但是其计算相对复杂,并且估计效果受到先验信息准确程度的影响,应用较为不便。
技术实现思路
本专利技术的目的在于提出一种广义帕累托分布参数的迭代最大似然估计方法,以提高最大似然估计的执行效率,提升后续海杂波背景下目标检测的性能。实现本专利技术目的的技术方案是:通过将杂波样本功率归一化,获取其形状参数和尺度参数之间的确定关系,并利用样本的迭代似然函数进行广义帕累托分布参数的估计,其步骤包括如下:(1)利用雷达发射机发射脉冲信号,利用雷达接收机接收经过海面散射形成的回波数据,该回波数据的每个分辨单元中的回波序列为X=[x1,x2,…xi,…xN],其中xi表示第i个回波数据,i=1,2,...,N,N表示脉冲数;(2)获取当前杂波数据的功率信息,并将其按功率进行归一化,得到功率归一化后的海杂波数据:Y=[y1,y2,…yi,…yN],其中yi是Y的第i个数据,是杂波样本功率PX的第i个数据,是杂波样本功率PX的平均值。(3)根据功率归一化后的杂波数据Y,利用矩估计方法计算尺度参数的迭代似然估计初值形状参数的迭代似然估计初值(4)根据尺度参数的迭代似然估计初值形状参数的迭代似然估计初值迭代计算尺度参数的迭代似然估计值和形状参数的迭代似然估计值4a)选取合适的参数估计精度ε;4b)计算尺度参数的第m次迭代似然估计值形状参数的第m次迭代似然估计值其中,m表示迭代次数。4c)判断是否满足收敛性条件:若这两式同时成立,终止迭代过程,得到尺度参数的迭代似然估计值和形状参数的迭代似然估计值否则,令m=m+1,返回步骤(4b)进行第m+1次迭代。本专利技术通过迭代广义帕累托分布的似然函数,实现了广义帕累托分布参数的估计,与现有技术相比具有以下优点:1)相比于矩估计方法,提高了矩估计的精度,扩展了参数估计方法的适用范围;2)相比于最大似然估计,该方法具有很快的收敛速度,并且无需通过搜索的方式获得其最优解,运算速度快,能够适应雷达系统信号实时处理的要求;附图说明图1为本专利技术的实现流程图;图2为采用本专利技术和现有两种估计方法在不同参数取值下的估计效果对比;图3为采用本专利技术和现有两种估计方法在不同样本数量下的估计效果对比。具体实施方式下面结合附图对本专利技术作进一步说明:参照图1,本专利技术的实现步骤如下:步骤1,利用雷达发射机发射脉冲信号,利用雷达接收机接收经过海面散射形成的回波数据。回波数据是一个包括脉冲维,距离维和波位维的三维矩阵,每个距离维和波位维构成一个分辨单元,每个分辨单元中的回波序列为X:X=[x1,x2,...,xi,...,xN],其中xi表示第i个回波数据,N表示脉冲数。步骤2,获取当前杂波数据的功率信息,并将其按功率进行归一化。2a)计算当前样本数据X的功率PX:PX=|X|2=[|x1|2,|x2|2,…|xi|2,…|xn|2]其中xi表示第i个回波数据,i=1,2,...,N,N表示脉冲数,PX服从广义帕累托分布,广义帕累托分布定义式如下:其中,σ表示尺度参数,k表示形状参数;2b)计算当前杂波数据样本功率PX的平均功率其中表示PX的第i个数据;2c)根据平均功率和当前杂波样本数据X的功率PX,得到功率归一化后的杂波数据样本Y:其中,表示功率归一化后的杂波数据样本Y的第i个数据。步骤3,计算尺度参数的迭代似然估计初值形状参数的迭代似然估计值计算尺度参数的迭代似然估计初值形状参数的迭代似然估计初值可采用现有的参数估计方法,包括但不限于矩估计、最大似然估计,双分位点估计。本实例采用矩估计方法计算尺度参数的迭代似然估计初值形状参数的迭代似然估计初值其计算公式如下:其中表示Y的均值,表示Y的方差。步骤4,根据尺度参数的迭代似然估计初值形状参数的迭代似然估计初值迭代计算尺度参数的迭代似然估计值和形状参数的迭代似然估计值4a)设定参数估计精度ε;4b)计算尺度参数的第m次迭代似然估计值形状参数的第m次迭代似然估计值其中,m表示迭代次数。4c)判断是否满足收敛性条件:若这两式同时成立,则终止迭代过程,得到尺度参数的迭代似然估计值和形状参数的迭代似然估计值否则,令m=m+1,返回步骤(4b)进行第m+1次迭代。本专利技术的效果可以通过以下仿真实验进一步说明:1.仿真参数仿真实验采用仿真产生的广义帕累托数据。2.仿真实验内容仿真实验中分别采用本专利技术、矩估计以及最大似然估计方法对于仿真产生的帕累托分布数据进行参数的估计,通过相对误差和均方根误差比较三种不同方法的估计效果。实验1,使用matlab软件中的gprnd函数分别产生不同形状参数和尺度参数下的帕累托分布数据,测试样本数量为1000,分别使用本专利技术、矩估计以及最大似然估计对于仿真产生的帕累托分布数据的形状参数和尺度参数进行估计,通过比较参数估计的相对误差和均方根误差RMSE比较不同估计方法的效果,每个参数取值下的实验重复2000次,最终给出2000次实验相对误差和均方根误差RMSE的平均值,结果如图2,其中,图2(a)为用三种方法对于形状参数估计的相对误差随形状参数取值的变化曲线,其中横坐标表示形状参数取值,纵坐标表示相对误差图2(b)为用三种方法对于尺度参数估计的相对误差随尺度参数取值的变化曲线,其中横坐标表示尺度参数取值,纵坐标表示相对误差。本文档来自技高网
...
广义帕累托分布参数的迭代最大似然估计方法

【技术保护点】
一种广义帕累托分布参数的迭代最大似然估计方法,包括如下步骤:(1)利用雷达发射机发射脉冲信号,利用雷达接收机接收经过海面散射形成的回波数据,该回波数据的每个分辨单元中的回波序列为X=[x

【技术特征摘要】
1.一种广义帕累托分布参数的迭代最大似然估计方法,包括如下步骤:(1)利用雷达发射机发射脉冲信号,利用雷达接收机接收经过海面散射形成的回波数据,该回波数据的每个分辨单元中的回波序列为X=[x1,x2,…xi,…xN],其中xi表示第i个回波数据,i=1,2,...,N,N表示脉冲数;(2)获取当前杂波数据的功率信息,并将其按功率进行归一化,得到功率归一化后的海杂波数据:Y=[y1,y2,…yi,…yN],其中yi是Y的第i个数据,是杂波样本功率PX的第i个数据,是杂波样本功率PX的平均值。(3)根据功率归一化后的杂波数据Y,利用矩估计方法计算尺度参数的迭代似然估计初值形状参数的迭代似然估计初值(4)根据尺度参数的迭代似然估计初值形状参数的迭代似然估计初值迭代计算尺度参数的迭代似然估计值和形状参数的迭代似然估计值4a)设定参数估计精度ε;4b)计算尺度参数的第m次迭代似然估计值和形状参数的第m次迭代似然估计值其中,m表示迭代次数;4c)判断是否满足收敛性条件:若这两式同时成立,则终止迭代过程,得到尺度参数的迭代似然估计...

【专利技术属性】
技术研发人员:许述文王乐水鹏朗黎鑫
申请(专利权)人:西安电子科技大学
类型:发明
国别省市:陕西,61

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1