(19)中华 人民共和国 国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202111570488.8 (22)申请日 2021.12.21 (71)申请人 中国人民解 放军国防科技大 学 地址 410073 湖南省长 沙市开福区德雅路 109号 (72)发明人 刘永祥　杨立儒　杨威　沈亲沐　 张文鹏　 (74)专利代理 机构 长沙国科天河知识产权代理 有限公司 432 25 代理人 唐品利 (51)Int.Cl. G06F 30/20(2020.01) G06F 17/18(2006.01) (54)发明名称 基于原点矩偏导的K分布杂波参数估计方法 和装置 (57)摘要 本申请涉及一种基于原 点矩偏导的K分布杂 波参数估计方法、 装置、 计算机设备和存储介质。 所述方法包括： 根据零记忆非线性变换法获取幅 度服从K分布的杂波序列数据， 通过核平滑直方 图对杂波序列数据进行处理， 得到统计概率密度 值； 利用原点矩偏导的K分布杂波参数估计方法 获取杂波序列数据的形状参数和尺度参数的估 计值； 根据形状参数和尺度参数的估计值进行概 率密度计算， 得到估计概率密度值； 对统计概率 密度值和估计概率密度值进行均方根误差计算， 得到最小均方根误差； 最小均方根误差对应的形 状参数和尺度参数的估计值为最优的K分布杂波 参数。 采用本方法能够提高参数估计精度。 权利要求书2页 说明书11页 附图6页 CN 114239294 A 2022.03.25 CN 114239294 A 1.一种基于原点矩偏导的K分布杂波参数估计方法， 其特 征在于， 所述方法包括： 根据零记 忆非线性变换法获取幅度服从K分布的杂波序列数据； 通过核平 滑直方图对所述杂波序列数据进行处 理， 得到统计概 率密度值； 利用原点矩偏导的K分布杂波参数估计方法获取所述杂波序列 数据的形状参数和尺度 参数的估计值； 根据所述形状参数和尺度参数的估计值进行概 率密度计算， 得到估计概 率密度值； 对所述统计概率密度值和估计概率密度值进行均方根误差计算， 得到最小均方根误 差； 所述最小均方根 误差对应的形状参数和尺度参数的估计值 为最优的K分布杂波参数。 2.根据权利要求1所述的方法， 其特征在于， 通过核平滑直方图对所述杂波序列数据进 行处理， 得到统计概 率密度值， 包括： 对所述核平 滑直方图的统计区间进行定义， 得到统计区间大小； 根据所述统计区间大小、 所述杂波序列数据的幅度最大值和幅度最小值进行计算， 得 到统计概 率密度值。 3.根据权利要求2所述的方法， 其特征在于， 根据所述统计区间大小、 所述杂波序列数 据的幅度最大值和幅度最小值进行计算， 得到统计概 率密度值， 包括： [fz0,ri]＝ksdensity(x,mi:stp:ma) 其中， fz0表 示统计概率密度值， r i表示统计区间的中心值， ksdensity()表 示核平滑密 度估计函数， mi表示杂波序列数据的幅度最小值， stp表示统计区间大小， ma表示杂波序列 数据的幅度最大值， x表示杂波幅度。 4.根据权利要求1至3任意一项所述的方法， 其特征在于， 利用原点矩偏导的K分布杂波 参数估计方法获取 所述杂波序列数据的形状参数和尺度参数的估计值， 包括： 利用原点矩偏导的K分布杂波参数估计方法获取所述杂波序列数据的k阶原点矩和k阶 原点矩的偏导； 根据所述k阶原点矩和k阶原点矩的偏导， 以及多个函数构建形状参数的估计值的求解 方程； 令所述求解方程为0进行求 解， 得到所述杂波序列数据的形状参数的估计值； 根据所述形状参数的估计值和k阶原点矩， 得到所述杂波序列数据的尺度参数的估计 值。 5.根据权利要求4所述的方法， 其特征在于， 根据所述k阶原点矩和k阶原点矩的偏导， 以及多个函数构建形状参数的估计值的求 解方程， 包括： fun＝@(v)0.5*k*psi(0.5*k+v) ‑log(gamma(0.5*k+v))+log(gamma(v)) ‑k*(dk/d0) ‑ log(d0)‑0.5*k*psi(0.5 *k+1)+log(gamma(0.5*k+1)) v＝fsolve(fun,0.1) 其中， @(v)表示求解v， v表示形状参数的估计值， k表示原点矩的阶数， d0表示k阶原点 矩， dk表示 k阶原点矩的偏导， psi()表示普西函数， g amma()表示伽玛函数。 6.根据权利要求5所述的方法， 其特征在于， 根据 所述形状参数的估计值和k阶原点矩， 得到所述杂波序列数据的尺度参数的估计值， 包括： σ ＝0.5*(d0*gamma(v)/(g amma(0.5*k+1)*gamma(0.5*k+v)))^(1/k)权　利　要　求　书 1/2 页 2 CN 114239294 A 2其中， σ 表示尺度参数的估计值。 7.根据权利要求6所述的方法， 其特征在于， 根据 所述形状参数和尺度参数的估计值进 行概率密度计算， 得到估计概 率密度值， 包括： fz＝(2/( σ *g amma(v))).*(ri/(2*σ ) ).^v.*bes selk(v‑1,ri/σ ) 其中， bes selk()表示第二类贝 塞尔函数。 8.一种基于原点矩偏导的K分布杂波参数估计装置， 其特 征在于， 所述装置包括： 获取杂波序列数据模块， 用于根据零记忆非线性变换法获取幅度服从K分布的杂波序 列数据； 杂波序列数据处理模块， 用于通过核平滑直方图对所述杂波序列数据进行处理， 得到 统计概率密度值； 杂波序列估计模块， 用于利用原点矩偏导的K分布杂波参数估计方法获取所述杂波序 列数据的形状参数和尺度参数的估计值； 概率密度计算模块， 用于根据所述形状参数和尺度参数的估计值进行概率密度计算， 得到估计概 率密度值； K分布杂波参数估计模块， 用于对所述统计概率密度值和估计概率密度值进行均方根 误差计算， 得到最小均方根误差； 所述最小均方根误差对应的形状参数和尺度参数 的估计 值为最优的K分布杂波参数。 9.一种计算机设备， 包括存储器和处理器， 所述存储器存储有计算机程序， 其特征在 于， 所述处 理器执行所述计算机程序时实现权利要求1至7中任一项所述方法的步骤。 10.一种计算机可读存储介质， 其上存储有计算机程序， 其特征在于， 所述计算机程序 被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。权　利　要　求　书 2/2 页 3 CN 114239294 A 3