[发明专利]一种确定信号信源个数的估计方法在审

专利信息
申请号: 202111590451.1 申请日: 2021-12-23
公开(公告)号: CN114282167A 公开(公告)日: 2022-04-05
发明(设计)人: 李海鹏;张居成;韩云峰;郑翠娥 申请(专利权)人: 哈尔滨工程大学青岛船舶科技有限公司
主分类号: G06F17/16 分类号: G06F17/16;G06F17/18;G06N7/00;H04B17/00
代理公司: 哈尔滨市松花江专利商标事务所 23109 代理人: 张换男
地址: 266000 山东*** 国省代码: 山东;37
权利要求书: 查看更多 说明书: 查看更多
摘要:
搜索关键词: 一种 确定 信号 信源 个数 估计 方法
【说明书】:

一种确定信号信源个数的估计方法,它属于阵列信号处理技术领域。本发明解决了采用现有方法对信源个数进行估计时,获得的信源个数估计结果的准确率低的问题。本发明根据观测数据的概率密度和观测数据的特征结构的概率密度,来构造广义贝叶斯信息量准则表达式,计算所有可能的信源个数值对应的信息量值,得到最大的信息量值所对应的信源个数值。与现有方法相比,本发明方法可以获得更高的信源个数估计的准确率,为基于子空间的算法提供先验知识。本发明方法可以应用于对确定信号进行信源个数的估计。

技术领域

本发明属于阵列信号处理技术领域,具体涉及一种确定信号信源个数的估计方法。

背景技术

子空间估计是阵列信号处理中的重要研究内容,信源个数估计算法是子空间类算法的重要组成。基于信息量准则的算法是目前最常用的信源个数估计算法,常用的估计方法包括赤池信息量准则、经典贝叶斯信息量准则等,这些方法将信源个数的估计问题转化为模型参数的选择问题,在参数估计问题中,通常会采用似然函数作为目标函数,通过增加观测数据量提高模型精度,同时加入惩罚项来避免过拟合问题。

但是现有方法在构造目标函数时,没有充分考虑观测数据的概率密度和观测数据的特征结构的概率密度,使采用现有方法进行信源个数估计时获得的信源个数估计结果的准确率较低。

发明内容

本发明的目的是为解决采用现有方法对信源个数进行估计时,获得的信源个数估计结果的准确率低的问题,而提出了一种确定信号信源个数的估计方法。

本发明为解决上述技术问题所采取的技术方案是:

一种确定信号信源个数的估计方法,所述方法具体包括以下步骤:

步骤一、采用传感器阵列获得观测数据;

步骤二、对观测数据的协方差矩阵进行特征值分解,获得观测数据的特征值;

步骤三、根据特征值计算观测数据的概率密度以及观测数据的特征结构的概率密度,再根据概率密度构造广义贝叶斯信息量准则表达式;

所述观测数据的概率密度f1为:

其中,r表示信源个数,N表示传感器阵列中阵元的个数,L表示传感器阵列的每个阵元接收的数据的采样点数,λi表示观测数据的协方差矩阵的第i个特征值;

观测数据的特征结构的概率密度f2为:

其中,表示观测数据的协方差矩阵的第i个特征值的真值,λj表示观测数据的协方差矩阵的第j个特征值,σ2表示高斯白噪声的方差,Γ表示一种数学运算;

步骤四、将最大的广义贝叶斯信息量值所对应的r值作为信源个数的估计结果。

进一步地,所述观测数据中包含高斯白噪声和多个窄带信号,窄带信号的个数即信源个数。

进一步地,所述广义贝叶斯信息量准则表达式为:

B=-2lg(f1+f2)+Klg(L)

其中:B表示广义贝叶斯信息量值,K表示子空间的自由度。

进一步地,所述子空间的自由度K为:

更进一步地,所述数学运算Γ的定义为:

对于任意正整数n,Γ(n)=(n-1)·(n-2)·2·1。

本发明的有益效果是:

下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学青岛船舶科技有限公司,未经哈尔滨工程大学青岛船舶科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/pat/books/202111590451.1/2.html,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code back_top