CN107167781B - 海杂波幅度对数正态分布参数的分位点估计方法 - Google Patents
海杂波幅度对数正态分布参数的分位点估计方法 Download PDFInfo
- Publication number
- CN107167781B CN107167781B CN201710206425.1A CN201710206425A CN107167781B CN 107167781 B CN107167781 B CN 107167781B CN 201710206425 A CN201710206425 A CN 201710206425A CN 107167781 B CN107167781 B CN 107167781B
- Authority
- CN
- China
- Prior art keywords
- sea clutter
- amplitude
- quantile
- clutter amplitude
- 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.)
- Active
Links
Images
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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/414—Discriminating targets with respect to background clutter
Abstract
Description
技术领域
本发明属于雷达目标检测技术领域,具体涉及到一种海杂波幅度对数正态分布参数的分位点估计方法,用于海杂波背景下目标检测中确定海杂波幅度分布模型位置参数和尺度参数。
背景技术
海杂波是影响海面目标检测和目标运动参数估计的主要因素。海杂波背景下的最优目标检测方法依赖于海杂波幅度分布模型的模型参数,如何从复杂海杂波场景中有效精确估计出海杂波幅度分布模型的模型参数是海面目标检测的关键。海杂波的幅度分布模型随着雷达分辨率和海况而改变。在低分辨率条件下,复的海杂波一般服从复高斯分布,幅度分布模型为单参数的瑞利分布模型。瑞利分布与每个散射体的振幅分布无关,只要求散射体的数目足够多,并且所有散射体中没有一个起主导作用。对于低分辨力雷达,在高仰角和平稳环境时,瑞利分布的杂波模型可以得到较为精确的结果。但是,随着对雷达杂波幅度分布特性分析的逐步深入,人们发现,对于海浪杂波和地物杂波,瑞利分布模型并不能给出令人满意的结果。特别是随着距离分辨力的提高,杂波分布出现了比瑞利分布更长的拖尾,即出现高振幅的概率相当大。因而,如果继续采用瑞利分布模型,将出现较高的虚警概率。海杂波幅度的分布不仅是脉冲宽度的函数,而且也与雷达极化方式、工作频率、天线视角以及海情、风向和风速等因素有关。对于高分辨力雷达,在低仰角或恶劣海情下,海杂波不再服从瑞利分布,而通常能用对数正态分布来描述。
对于海杂波幅度对数正态分布模型,矩估计是常用的参数估计方法。幅度分布的矩能够表示为模型参数的非线性函数,因此联立方程可以用样本矩表示位置和尺度参数。在文献“对数正态分布参数的最大似然估计,”参见“于洋,孙月静.对数正态分布参数的最大似然估计[J].九江学院学报,2007,26(6):55-57.”中,提出了对数正态分布参数的最大似然的估计方法,它具有比矩估计更高的精度。矩估计和最大似然估计方法都是基于基本假定:所使用的是纯杂波数据。然而,实际应用中这一基本假定是不成立的,获得的海杂波数据样本常常包含了少量具有很大幅度的岛礁回波、目标回波构成的异常散射单元等。这些异常样本使得矩估计和最大似然估计方法的估计精度急剧下降,导致最优检测无法有效实现。
发明内容
本发明的目的在于针对上述现有技术的不足,提出一种海杂波幅度对数正态分布参数的分位点估计方法,以提高在岛礁、目标回波组成的异常散射单元或样本存在条件下,海杂波幅度对数正态分布位置和尺度参数的估计精确和稳健性。
为实现上述技术目的,本发明的技术方案包括如下:
(1)对海杂波幅度对数正态分布模型的概率密度函数fR(r)进行积分,得到海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r):
r>0,μ∈R,σ>0,
其中,Prob(·)表示事件的概率值,r′为样本幅度,0≤r′≤r,erfc(·)表示互补误差函数,μ和σ分别表示海杂波幅度对数正态分布的位置参数和尺度参数;
(2)根据(1)中得到的海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r),得到参数估计所需的方程:
其中erfc-1()为互补误差函数的逆函数。
本发明与现有技术相比具有如下优点:
本发明由于在估计参数时只选取了样本的中位数和一个分位数,排除了异常散射单元或样本对参数估计的影响,使得本发明的方法具有抗异常散射单元或样本的能力,在海杂波样本包含少量大幅度的岛礁、目标回波情况下,仍然能够获得位置参数和尺度参数的高精度稳健估计。
附图说明
图1为本发明的实现流程图;
图2为实测海杂波幅度在对数正态分布模型下,采用本发明和现有方法得到的幅度概率密度函数曲线图;
图3为实测海杂波幅度在对数正态分布模型下,采用本发明和现有方法得到的幅度累积分布函数曲线图。
具体实施方式
下面结合附图对本发明作进一步说明:
参照图1,本发明的实现步骤如下:
步骤1,对海杂波幅度对数正态分布模型的概率密度函数fR(r)进行积分,得到海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r)的表达形式。
(1.1)将海杂波幅度对数正态分布模型的概率密度函数fR(r),表示如下:
r>0,μ∈R,σ>0,
其中r表示海杂波幅度,μ是位置参数,σ是尺度参数;
(1.2)对式<1>进行积分,得到海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r):
其中,Prob(·)表示事件的概率值,r′为样本幅度,0≤r′≤r,erfc(·)表示互补误差函数。
步骤2,根据步骤1中得到的海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r),得到参数估计所需的方程。
rm=eμ <3>
(2.2)利用累积分布函数Fμ,σ(r),将满足α=Fμ,σ(rα),0<α<1的rα定义为分位点α对应的分位数,根据该定义,得到方程:
本发明实例中取α=0.75。
(3.1)利用雷达发射机发射连续的脉冲信号,经过海面散射形成回波,利用雷达接收机接收回波数据矩阵X,其中X是一个K×N的矩阵,K表示回波数据矩阵的距离单元数,N表示回波数据矩阵的积累脉冲数;
(3.2)在回波数据矩阵X中,选取只包含杂波数据的距离单元,该距离单元的数据为N个杂波幅度数据,分别表示为r1,r2,....,rN,对这N个杂波幅度数据进行升序排列,得到一个递增的序列r(1),r(2),....,r(N);
其中,round(·)表示最接近·的整数。
步骤4,根据步骤3中得到的中位数rm和分位数rα的估计值,计算位置参数和尺度参数的估计值。
联立<3>和<4>,得到如下方程组:
其中erfc-1(·)表示互补误差函数的逆函数。
基于步骤1到步骤4,完成了对海杂波幅度对数正态分布模型的参数估计。
下面结合仿真实验对本发明的效果做进一步说明。
1、仿真参数
仿真实验使用了IPIX雷达录取的实测海杂波数据19931109_191449.mat中第1个距离单元HH极化数据,该数据为纯杂波数据,不含目标。雷达工作在驻留模式,脉冲重复频率为1000Hz,距离分辨率为30m;数据中包含了14个距离单元和131072个连续的相参脉冲串数据。
2、仿真实验内容
从纯杂波数据中取N个幅度数据r1,r2,....,rN,其中N=131072。仿真实验中,分别使用矩估计法和本发明的分位点估计法对实测海杂波数据进行参数估计,并拟合概率密度函数曲线和累积分布函数曲线。
实验一:利用矩估计法和本发明的分位点估计法进行参数估计,并拟合实测海杂波数据的概率密度函数曲线,实验结果如图2,其中横轴表示海杂波幅度,纵轴表示海杂波幅度对数正态分布模型概率密度的值。图2中“·”标记的曲线表示实测杂波数据的经验概率密度函数曲线,“---”标记的曲线表示采用本发明得到概率密度函数曲线,“—”标记的曲线表示采用矩估计法得到的概率密度函数曲线,
从图2可见,利用矩估计法得到位置参数估计值为-0.7106,尺度参数估计值为1.0282,利用本发明方法得到的位置参数估计值为-0.8049,尺度参数估计值为1.1273。本发明得到的概率密度函数曲线更加接近实测海杂波数据概率密度函数曲线。
实验二:利用矩估计法和本发明的分位点估计法拟合实测海杂波数据的累积分布函数曲线,实验结果如图3,其中横轴表示海杂波幅度,纵轴表示海杂波幅度对数正态分布模型累积分布函数的值。图3中“·”标记的曲线表示实测杂波数据样本的经验累积分布函数,“---”标记的曲线表示采用本发明得到的累积分布函数曲线图,“—”标记的曲线表示采用矩估计法得到累积分布函数曲线图。
从图3可见,两种方法得到的累积分布函数均接近于实测数据的经验累积分布函数,于是对于两种估计方法,引入KS检验法来证明哪一种估计法更接近实测海杂波数据经验累积分布函数,KS距离越小表示对应的海杂波幅度对数正态分布模型越接近于实测海杂波幅度模型。计算得知,本发明对应的累积分布函数与经验累积分布函数之间的KS距离为0.0293,而矩估计法对应的分布函数与经验累积分布函数之间的KS距离为0.0758,明显高于本发明。
综上,本发明提出了海杂波幅度对数正态分布模型参数的分位点估计方法具有抗异常散射单元或样本的能力,在海杂波样本包含少量大幅度的岛礁、目标回波情况下,仍然能够获得形状参数和尺度参数的高精度稳健估计。
Claims (3)
1.一种海杂波幅度对数正态分布参数的分位点估计方法,其特征在于,包括以下步骤:
(1)对海杂波幅度对数正态分布模型的概率密度函数fR(r)进行积分,得到海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r):
其中,Prob(·)表示事件的概率值,r′为样本幅度,r为海杂波幅度,0≤r′≤r,erfc(·)表示互补误差函数,μ和σ分别表示海杂波幅度对数正态分布的位置参数和尺度参数;
(2)根据(1)中得到的海杂波幅度对数正态分布模型累积分布函数Fμ,σ(r),得到参数估计所需的方程:
其中erfc-1()为互补误差函数的逆函数。
(3a)利用雷达发射机发射连续的脉冲信号,经过海面散射形成回波,利用雷达接收机接收回波数据矩阵X,其中X是一个K×N的矩阵,K表示回波数据矩阵的距离单元数,N表示回波数据矩阵的积累脉冲数;
(3b)在回波数据矩阵X中,选取只包含杂波数据的距离单元,该距离单元的数据为N个杂波幅度数据,分别表示为r1,r2,....,rN,对这N个杂波幅度数据进行升序排列,得到一个递增的序列r(1),r(2),....,r(N);
其中,round(·)表示最接近·的整数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710206425.1A CN107167781B (zh) | 2017-03-31 | 2017-03-31 | 海杂波幅度对数正态分布参数的分位点估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710206425.1A CN107167781B (zh) | 2017-03-31 | 2017-03-31 | 海杂波幅度对数正态分布参数的分位点估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107167781A CN107167781A (zh) | 2017-09-15 |
CN107167781B true CN107167781B (zh) | 2020-04-14 |
Family
ID=59849089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710206425.1A Active CN107167781B (zh) | 2017-03-31 | 2017-03-31 | 海杂波幅度对数正态分布参数的分位点估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107167781B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107741581B (zh) * | 2017-09-22 | 2020-10-09 | 西安电子科技大学 | 基于截断矩的广义帕累托分布参数估计方法 |
CN109446481B (zh) * | 2018-09-17 | 2023-05-12 | 中国人民解放军海军工程大学 | 一种对数正态型单元寿命分布参数估计方法 |
CN109143196B (zh) * | 2018-09-25 | 2023-01-06 | 西安电子科技大学 | 基于k分布海杂波幅度模型的三分位点参数估计方法 |
CN110275148B (zh) * | 2019-06-20 | 2021-08-27 | 中国人民解放军海军航空大学 | 一种海杂波幅度分布参数估计方法及系统 |
CN110796381A (zh) * | 2019-10-31 | 2020-02-14 | 深圳前海微众银行股份有限公司 | 建模数据评价指标的处理方法、装置、终端设备及介质 |
CN110879927B (zh) * | 2019-11-20 | 2023-08-08 | 中国海洋大学 | 一种用于海目标检测的海杂波幅度统计分布现场建模方法 |
CN111381216B (zh) * | 2020-03-10 | 2022-03-11 | 中国科学院微电子研究所 | 混合分布的雷达海杂波分析方法及装置 |
CN111929671B (zh) * | 2020-06-30 | 2022-02-18 | 湖南华诺星空电子技术有限公司 | 基于正态匹配的无线电侦测的接收幅度的估计方法及系统 |
CN111830479B (zh) * | 2020-07-09 | 2023-03-07 | 中国人民解放军海军航空大学 | 雷达回波双分量混合幅度分布模型参数估计方法及装置 |
CN112881988A (zh) * | 2021-01-11 | 2021-06-01 | 西北工业大学 | 导航雷达模拟训练系统中杂波模拟显示方法 |
CN113640763B (zh) * | 2021-07-29 | 2023-08-04 | 西安邮电大学 | 基于分数阶矩的对数正态纹理海杂波幅度分布形状参数的估计方法 |
CN116736781B (zh) * | 2023-08-15 | 2023-11-03 | 国网浙江省电力有限公司杭州供电公司 | 一种工业自动化控制设备安全状态监测方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535981A (zh) * | 2015-01-21 | 2015-04-22 | 西安电子科技大学 | 海杂波Pareto幅度分布参数的双分位点估计方法 |
CN104749564A (zh) * | 2015-04-10 | 2015-07-01 | 西安电子科技大学 | 海杂波Weibull幅度分布参数的多分位点估计方法 |
CN105699952A (zh) * | 2016-01-25 | 2016-06-22 | 西安电子科技大学 | 海杂波k分布形状参数的双分位点估计方法 |
CN106199537A (zh) * | 2016-06-29 | 2016-12-07 | 西安电子科技大学 | 基于逆高斯纹理海杂波幅度分布参数的分位点估计方法 |
-
2017
- 2017-03-31 CN CN201710206425.1A patent/CN107167781B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104535981A (zh) * | 2015-01-21 | 2015-04-22 | 西安电子科技大学 | 海杂波Pareto幅度分布参数的双分位点估计方法 |
CN104749564A (zh) * | 2015-04-10 | 2015-07-01 | 西安电子科技大学 | 海杂波Weibull幅度分布参数的多分位点估计方法 |
CN105699952A (zh) * | 2016-01-25 | 2016-06-22 | 西安电子科技大学 | 海杂波k分布形状参数的双分位点估计方法 |
CN106199537A (zh) * | 2016-06-29 | 2016-12-07 | 西安电子科技大学 | 基于逆高斯纹理海杂波幅度分布参数的分位点估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107167781A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107167781B (zh) | 海杂波幅度对数正态分布参数的分位点估计方法 | |
CN109116311B (zh) | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 | |
CN111913158B (zh) | 复杂杂波背景下检测低慢小目标的雷达信号处理方法 | |
CN105699952B (zh) | 海杂波k分布形状参数的双分位点估计方法 | |
CN102955150B (zh) | 基于天线主瓣占优强度约束的sar方位模糊抑制方法 | |
CN109324315B (zh) | 基于双层次块稀疏性的空时自适应处理雷达杂波抑制方法 | |
Martorella et al. | On the fractal dimension of sea surface backscattered signal at low grazing angle | |
CN104749564A (zh) | 海杂波Weibull幅度分布参数的多分位点估计方法 | |
CN108387884B (zh) | 基于知识辅助稀疏渐进最小方差的机载雷达杂波抑制方法 | |
CN106569193B (zh) | 基于前-后向收益参考粒子滤波的海面小目标检测方法 | |
CN110673108B (zh) | 一种基于迭代Klett的机载海洋激光雷达信号处理方法 | |
CN109143224A (zh) | 一种多目标关联方法和装置 | |
CN108469608A (zh) | 一种运动平台雷达多普勒质心精确估计方法 | |
CN107741581B (zh) | 基于截断矩的广义帕累托分布参数估计方法 | |
CN109709526B (zh) | 一种知识辅助分组广义似然比检测方法 | |
CN107064893B (zh) | 基于对数矩的广义帕累托分布参数估计方法 | |
CN108761417B (zh) | 基于知识辅助最大似然的机载雷达杂波抑制方法 | |
CN106199552A (zh) | 一种局部均匀海杂波背景下的分组广义似然比检测方法 | |
CN111796288B (zh) | 一种基于杂波频谱补偿技术的三坐标雷达动目标处理方法 | |
CN106093927B (zh) | 基于雷达大脉宽信号的目标测速测距方法 | |
CN110244289A (zh) | 一种自适应粒子滤波地波雷达目标一体化探测方法 | |
CN108008374B (zh) | 基于能量中值的海面大型目标检测方法 | |
CN114296046B (zh) | 基于人工神经网络的hfswr多海况有效波高提取方法与装置 | |
CN114252878A (zh) | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 | |
CN114518564A (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 |