一种基于可变窗的局部拟合优度检验方法技术

技术编号:17344738 阅读:31 留言:0更新日期:2018-02-25 09:36
本发明专利技术公开了一种基于可变窗的局部拟合优度检验方法。拟合优度检验是对海杂波幅度分布模型的定量判断,克服了均方误差检验(MSD)和修正均方检验(MMSD)方法对幅度观测区间较为单一的问题,通过选取一个或多个雷达感兴趣的局部幅度观测区间,计算出所选局部区间拟合优度检验的加权结果,该方法增强了统计模型拟合优度检验结果的全面性,为海杂波幅度序列统计模型的判断和海面目标恒虚警检测(CFAR)方法的设计提供了更为可靠的依据。

【技术实现步骤摘要】
一种基于可变窗的局部拟合优度检验方法
本专利技术涉及一种基于可变窗的海杂波幅度统计分布的局部拟合优度检验方法,属于雷达

技术介绍
海上目标探测是雷达的一项重要任务,然而海杂波是影响雷达目标探测的一个重要因素,为了能够有效地实现海杂波背景下的目标检测,必须要分析海杂波的统计分布模型。对于实测海杂波数据而言,仅通过肉眼观测从多种分布模型来选择拟合最好的海杂波幅度模型是不准确的,因此需要一种定量描述海杂波幅度模型与实测数据拟合程度的统计量,即拟合优度检验。传统的经典拟合优度检验方法主要包括:MSD法和MMSD法。然而MSD法仅仅关注了海杂波幅度统计的全局统计拟合,忽略了幅度局部区间的拟合效果;MMSD法仅仅关注了海杂波大幅度序列的拟合效果,即“拖尾”拟合效果。以上两种典型拟合优度检验方法仅单一地考虑了海杂波幅度的全局或尾部的拟合结果,忽略了不同海况条件下海杂波局部幅度区间内的拟合情况。雷达在实际工作时,会面临不同海况条件下的复杂海杂波回波信号,仅单一地从全局或尾部去分析海杂波的模型拟合结果是不全面的,也是不准确的。
技术实现思路
本专利技术的技术解决问题是:克服现有技术的不足,提供一种基于可变窗的局部拟合优度检验方法,能够对不同海况条件下海杂波局部幅度区间内的拟合情况进行检验,增强了拟合优度检验的可靠性和全面性,为目标恒虚警检测提供更加可靠的依据。本专利技术的技术解决方案是:一种基于可变窗的局部拟合优度检验方法,包括如下步骤:(1)从雷达实测的海杂波回波数据中提取海杂波幅度序列X={xk,k=1,2,3,...,N},将海杂波幅度序列等分为L个区间,得到海杂波幅度统计的实际概率密度分布pe(xl),其中N为海杂波幅度序列的长度,l=1,2,3,…,L;(2)根据海杂波幅度序列,采用矩估计方法计算出序列的幅度统计分布模型参数,利用常用拟合方法拟合出统计模型的海杂波理论概率分布函数pt(xl);(3)确定局部观测窗的个数、位置和长度;(4)针对每种拟合方法拟合出的海杂波理论概率分布函数pt(xl),计算每个观测窗的局部拟合优度检验结果,并通过加权处理得到每种拟合方法拟合出的海杂波理论概率分布函数pt(xl)与实际概率密度分布pe(xl)的最终拟合结果;(5)最终拟合结果的最小值对应的海杂波理论概率分布函数pt(xl)即为拟合最好的海杂波幅度序列统计分布模型。所述步骤(1)中,将海杂波幅度序列等分为L个区间,得到海杂波幅度统计的实际概率密度分布pe(xl)的方法如下:按照海杂波幅度序列中的最大值和最小值将海杂波幅度序列等分为L个区间,统计落入每个区间的海杂波幅度个数,得到海杂波幅度统计的实际概率密度分布函数pe(xl),pe(xl)=落入第l个区间的海杂波幅度个数/N。其中,区间L>20。所述步骤(2)中,局部观测窗的个数等于雷达感兴趣的局部幅度序列区域个数,第i个局部观测窗的位置[Ni1,Ni2]与雷达感兴趣的第i个局部幅度序列区域重合,Ni1,Ni2分别表示第i个局部观测窗的起始位置和终止位置,且第i个局部观测窗的长度leni=Ni2-Ni1≥4,Ni2、Ni1均为正整数。所述步骤(4)中,利用如下公式计算第i个局部观测窗的拟合优度检验结果:所述步骤(4)中,利用如下公式计算第j种拟合方法的最终拟合结果Dlmsd:Dlmsd=α1Dlmsd,1+α2Dlmsd,2+…+αiDlmsd,i+…+αMDlmsd,M0<αi<1,且α1+α2+…+αi+…+αM<1,M为局部观测窗的总个数。αi的取值与海况等级成正比。与现有技术相比,本专利技术具有如下有益效果:本专利技术克服了传统方法幅度观测区间单一的问题,提出了加窗处理的思想,可根据雷达实际观测需求来选择局部观测窗的个数、位置和长度,通过计算加窗后的拟合优度检验结果,匹配出最适合的分布模型,能够对不同海况条件下海杂波局部幅度观测区间内的拟合情况进行检验。当雷达需兼顾多组不同区间的拟合结果(多个局部区间或多个局部区间和整个幅度区间),可对多组区间的联合拟合结果进行加权评估,增强了统计模型拟合优度检验结果的全面性,为海面目标恒虚警检测算法的设计提供有力的支撑。附图说明图1为基于可变窗的局部统计拟合优度方法流程图;图2为实测海杂波幅度序列的统计拟合结果;图3为实测海杂波数据的拟合优度检验结果,(a)为四种常见分布的MSD检验结果,(b)为四种常见分布的MMSD检验结果,(c)为四种常见分布的局部拟合优度检验结果。具体实施方式考虑到不同海况条件下的模型拟合情况,本专利技术提出一种基于可变窗的局部拟合优度检验(LMSD)方法,克服了现有方法对幅度观测区域单一的问题,可以根据雷达的观测需要,通过加窗处理来计算局部幅度观测区间的拟合优度检验结果;当对多个不同幅度观测区间感兴趣时,可以通过对多个观测区间加权处理,计算多个区间共同作用下的幅度统计拟合优度检验结果,增强了拟合优度检验的可靠性和全面性,为海杂波幅度统计模型的判断和海面目标恒虚警检测方法的设计提供了更加可靠的依据。如图1所示,本专利技术主要实现步骤如下:步骤1从雷达实测的海杂波回波数据中提取海杂波幅度序列X={xk,k=1,2,3,...,N},N是海杂波幅度序列的长度。按照海杂波幅度序列中的最大值和最小值将海杂波幅度序列等分为L个区间,统计海杂波幅度在每个区间的出现频次,得到海杂波幅度统计的实际概率密度分布函数pe(xl)。其中,L选择不能太小,L较小会导致幅度分布函数的拟合结果不准确,因此,一般选择L>20。l=1,2,3,…,L。相邻两个区间的间隔为(幅度序列中的最大值-最小值)/L,pe(xl)=落入第l个区间的海杂波幅度个数除以海杂波幅度总数N。步骤2海杂波幅度统计拟合:根据海杂波幅度序列,采用矩估计方法计算出序列的幅度统计分布模型参数,利用常用拟合方法拟合出统计模型的海杂波理论概率分布函数pt(xl)。步骤3确定局部观测窗的个数、位置和长度:3.1选择观测窗个数M。若雷达仅对海杂波幅度统计的单一局部区域感兴趣,则选择单个矩形窗,即M=1。若雷达对海杂波幅度统计的多个局部区域感兴趣,则为每个区域选择一个矩形窗,即M>1,然后对多个局部区域拟合结果进行加权处理。本专利技术主要以一组或两组观测区间的拟合结果为例进行分析,若对更多的观测区间感兴趣,M还可以选更大的值,分析方法类似于M=2的情况。3.2根据3.1确定观测窗个数后,可根据雷达感兴趣区域选择观测窗位置。当雷达对海杂波某一局部幅度统计分布感兴趣时,可将矩形窗加在任意观测区间上;当雷达对海杂波的全局拟合统计分布感兴趣时可将矩形窗加在整个观测区间上。3.3局部观测窗长度的选择。设第i个局部观测窗的位置为[Ni1,Ni2],Ni1,Ni2分别表示第i个局部观测窗的起始位置和终止位置,Ni2、Ni1均为正整数,且第i个局部观测窗的长度leni=Ni2-Ni1≥4,即加窗后序列的幅值范围应覆盖整个序列幅值的4N/L以上,这样能保证足够的幅度样本参与运算,降低噪声及样本数太少导致的随机误差影响。步骤4计算每个观测窗的局部拟合优度检验(LMSD)结果对于第j种拟合方法,其第i个局部观测窗的拟合优度检本文档来自技高网...
一种基于可变窗的局部拟合优度检验方法

【技术保护点】
一种基于可变窗的局部拟合优度检验方法,其特征在于包括如下步骤:(1)从雷达实测的海杂波回波数据中提取海杂波幅度序列X={xk,k=1,2,3,...,N},将海杂波幅度序列等分为L个区间,得到海杂波幅度统计的实际概率密度分布pe(xl),其中N为海杂波幅度序列的长度,l=1,2,3,…,L;(2)根据海杂波幅度序列,采用矩估计方法计算出序列的幅度统计分布模型参数,利用常用拟合方法拟合出统计模型的海杂波理论概率分布函数pt(xl);(3)确定局部观测窗的个数、位置和长度;(4)针对每种拟合方法拟合出的海杂波理论概率分布函数pt(xl),计算每个观测窗的局部拟合优度检验结果,并通过加权处理得到每种拟合方法拟合出的海杂波理论概率分布函数pt(xl)与实际概率密度分布pe(xl)的最终拟合结果;(5)最终拟合结果的最小值对应的海杂波理论概率分布函数pt(xl)即为拟合最好的海杂波幅度序列统计分布模型。

【技术特征摘要】
1.一种基于可变窗的局部拟合优度检验方法,其特征在于包括如下步骤:(1)从雷达实测的海杂波回波数据中提取海杂波幅度序列X={xk,k=1,2,3,...,N},将海杂波幅度序列等分为L个区间,得到海杂波幅度统计的实际概率密度分布pe(xl),其中N为海杂波幅度序列的长度,l=1,2,3,…,L;(2)根据海杂波幅度序列,采用矩估计方法计算出序列的幅度统计分布模型参数,利用常用拟合方法拟合出统计模型的海杂波理论概率分布函数pt(xl);(3)确定局部观测窗的个数、位置和长度;(4)针对每种拟合方法拟合出的海杂波理论概率分布函数pt(xl),计算每个观测窗的局部拟合优度检验结果,并通过加权处理得到每种拟合方法拟合出的海杂波理论概率分布函数pt(xl)与实际概率密度分布pe(xl)的最终拟合结果;(5)最终拟合结果的最小值对应的海杂波理论概率分布函数pt(xl)即为拟合最好的海杂波幅度序列统计分布模型。2.根据权利要求1所述的一种基于可变窗的局部拟合优度检验方法,其特征在于:所述步骤(1)中,将海杂波幅度序列等分为L个区间,得到海杂波幅度统计的实际概率密度分布pe(xl)的方法如下:按照海杂波幅度序列中的最大值和最小值将海杂波幅度序列等分为L个区间,统计落入每个区间的海杂波幅度个数,得到海杂波幅度统计的实际概率密度分布函数pe(xl),pe(xl)=落入第l...

【专利技术属性】
技术研发人员:范一飞王伟伟杨晓超李财品黎薇萍李奇李渝
申请(专利权)人:西安空间无线电技术研究所
类型:发明
国别省市:陕西,61

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

1