matlab 广义帕累托分布,基于对数矩的广义帕累托分布参数估计方法与流程

本发明属于信号处理技术领域，具体涉及一种广义帕累托分布参数估计方法，可用于海杂波背景下的目标检测。

背景技术：

海杂波背景下的目标检测技术是雷达应用技术中一个至关重要的研究方向，在军事和民用领域已经得到广泛应用。而对于海杂波统计特性的准确分析是海杂波背景下目标检测技术能否取得良好效果的重要因素。因此，给出合适的模型并对于其模型参数进行准确估计成为我们需要解决的重要问题。

随着现代雷达系统距离分辨力的提高，雷达回波出现以往低分辨力雷达系统所没有的统计特性，通常表现为其回波包络的拖尾变长，异常值变多的特点。而广义帕累托分布作为复合高斯模型的一种，在对于高分辨低擦地角海杂波的功率分布拟合上取得了很好的效果。因此在海杂波统计特性的研究中占据重要地位。而在海杂波背景下的目标检测中，杂波模型参数的估计质量又对于目标检测效果有很大影响，因此在重拖尾的杂波数据下给出广义帕累托的分布参数具有重要的研究意义。

近年来，很多研究者对广义帕累托分布的参数估计方法，提出了一些基于特定条件下的广义帕累托分布参数估计理论。

文献“Castillo,E.,Hadi,A.S.,1997.Fitting the generalized Pareto distribution to data.J.Amer.Statist.Assoc.92,1609–1620.”中给出广义帕累托分布的矩估计以及最大似然估计方法，分别根据样本矩以及似然函数对于参数进行估计，但是由于矩估计本身容易受到样本数量和异常数据的影响，其估计精度难以保证。而最大似然估计的估计精度虽然能够满足要求，但是算法时间复杂度高，因此工程实现较为困难。

文献“Arnold,B.C.,Press,S.J.,1989.Bayesian estimation and prediction for Pareto data.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的一阶对数矩估计量κ1、二阶对数矩估计量κ2和均值κ3：

(4)利用功率归一化后海杂波数据Y的一阶对数矩估计量κ1、二阶对数矩估计量κ2和均值κ3计算形状参数的估计值和尺度参数的估计值

本发明通过利用对数域内的样本信息，实现广义帕累托分布参数的估计，与现有技术相比具有以下优点：

1)相比于矩估计方法，降低了矩估计的阶数，提高了参数估计的精度；

2)相比于最大似然估计方法，本发明具有解析表达式，无需通过搜索的方式获得其最优解，运算速度快，能够适应雷达系统信号实时处理的要求；

附图说明

图1为本发明的实现流程图；

图2为采用本发明和现有两种估计方法在不同参数取值下的估计效果对比；

图3为采用本发明和现有两种估计方法在不同样本数量下的估计效果对比。

具体实施方式

下面结合附图对本发明作进一步说明：

参照图1，本发明的实现步骤如下：

步骤1，利用雷达发射机发射脉冲信号，利用雷达接收机接收经过海面散射形成的回波数据。

回波数据是一个包括脉冲维，距离维和波位维的三维矩阵，每个距离维和波位维构成一个分辨单元，每个分辨单元中的回波序列为X：

X＝[x1,x2,...,xi,...,xN]

其中xi表示第i个回波数据，N表示脉冲数。

步骤2，获取当前杂波数据的功率信息，并将其按功率进行归一化，得到功率归一化后的杂波数据样本Y。

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的平均功率

其中pXi表示PX的第i个数据；

2c)根据平均功率和当前杂波样本数据X的功率PX，得到功率归一化后的杂波数据样本Y：

其中，表示功率归一化后的杂波数据样本Y的第i个数据。

步骤3，计算功率归一化后的海杂波数据Y的一阶对数矩估计量κ1、二阶对数矩估计量κ2和均值κ3：

步骤4，利用功率归一化后的海杂波数据Y的一阶对数矩估计量κ1、二阶对数矩估计量κ2和均值κ3，计算形状参数的估计值和尺度参数的估计值

本发明的效果可以通过以下仿真实验进一步说明：

1.仿真参数

仿真实验采用仿真产生的广义帕累托数据。

2.仿真实验内容

仿真实验中分别采用本发明、矩估计以及最大似然估计方法对于仿真产生的帕累托分布数据进行参数的估计，通过相对误差和均方根误差比较三种不同方法的估计效果。

实验1，使用matlab软件中的gprnd函数分别产生不同形状参数和尺度参数下的帕累托分布数据，测试样本数量为1000，分别使用本发明、矩估计以及最大似然估计对于仿真产生的帕累托分布数据的形状参数和尺度参数进行估计，通过比较参数估计的相对误差和均方根误差RMSE比较不同估计方法的效果，每个参数取值下的实验重复2000次，最终给出2000次实验相对误差和均方根误差RMSE的平均值，结果如图2，其中，

图2(a)为用三种方法对于尺度参数估计的相对误差随尺度参数取值的变化曲线，其中横坐标表示尺度参数取值，纵坐标表示相对误差

图2(b)为用三种方法对于形状参数估计的相对误差随形状参数取值的变化曲线，其中横坐标表示形状参数取值，纵坐标表示相对误差。

图2(c)为用三种方法对于尺度参数估计的均方根误差RMSE随尺度参数取值的变化曲线，其中横坐标表示尺度参数取值，纵坐标表示均方根误差RMSE。

图2(d)为用三种方法对于形状参数估计的均方根误差RMSE随形状参数取值的变化曲线，其中横坐标表示形状参数取值，纵坐标表示均方根误差RMSE。

实验2，使用matlab软件中的gprnd函数分别产生不同样本数量下的帕累托分布数据，形状参数取值为0.4，尺度参数取值为0.6，分别使用本发明、矩估计以及最大似然估计对于仿真产生的帕累托分布数据的形状参数和尺度参数进行估计，通过比较参数估计的相对误差和均方根误差RMSE比较不同估计方法的效果，每个样本数量下的实验重复2000次，最终给出2000次实验相对误差和均方根误差RMSE的平均值。结果如图3，其中，

图3(a)为用三种方法对于尺度参数估计的相对误差随样本数量的变化曲线，其中横坐标表示样本数量，纵坐标表示相对误差。

图3(b)为用三种方法对于形状参数估计的相对误差随样本数量的变化曲线，其中横坐标表示样本数量，纵坐标表示相对误差。

图3(c)为用三种方法对于尺度参数估计的均方根误差RMSE随样本数量的变化曲线，其中横坐标表示样本数量，纵坐标表示均方根误差RMSE。

图3(d)为用三种方法对于形状参数估计的均方根误差RMSE随样本数量的变化曲线，其中横坐标表示样本数量，纵坐标表示均方根误差RMSE。

从图2和图3中可以看出，本发明得到的参数估计精度高于矩估计且接近于最大似然估计。表明本发明提出的基于对数矩的广义帕累托分布参数估计方法，可以通过降低估计的阶数，提高帕累托分布参数的估计精度，并且计算速度快，能够满足雷达系统的实时处理要求，有利于后续海杂波背景下目标检测性能的提高。