[发明专利]一种声源数量的估计方法及系统在审
| 申请号: | 202210445311.3 | 申请日: | 2022-04-26 |
| 公开(公告)号: | CN114724585A | 公开(公告)日: | 2022-07-08 |
| 发明(设计)人: | 潘超;赵旭东;陈景东 | 申请(专利权)人: | 西北工业大学 |
| 主分类号: | G10L25/18 | 分类号: | G10L25/18;G10L25/21;G10L21/0216;G06F17/18 |
| 代理公司: | 北京众达德权知识产权代理有限公司 11570 | 代理人: | 吴莹 |
| 地址: | 710072 陕西*** | 国省代码: | 陕西;61 |
| 权利要求书: | 查看更多 | 说明书: | 查看更多 |
| 摘要: | |||
| 搜索关键词: | 一种 声源 数量 估计 方法 系统 | ||
1.一种声源数量的估计方法,其特征在于,所述方法包括:
通过麦克风阵列接收获得第一信号ym(t);
对所述第一信号进行短时傅里叶变换,获得变换域信号Ym(k,l);
在所述变换域信号内的多个频点上分别进行单声源DOA估计,获得估计结果;
在所述估计结果内,统计获得预设时间周期内不同声源方位被估计获得的频率,获得统计频率结果Q(θ);
查找获得所述统计频率结果Q(θ)内的峰值点,获取被估计获得的频率大于门限值τ1的多个声源方位,获得声源数量估计结果;
基于所述声源数量估计结果,更新当前活跃的声源个数。
2.根据权利要求1所述的方法,其特征在于,所述短时傅里叶变换,通过下式进行计算:
其中,m为所述麦克风阵列内第m个麦克风,k为频率索引,l为帧索引,w(n)为一窗函数,N为窗长,L为跳跃长度。
3.根据权利要求1所述的方法,其特征在于,所述单声源DOA估计包括:
计算:Pk,l(θ)=dH(θ)y(k,l)y(k,l)Hd(θ),其中,d(θ)为所述麦克风阵列的流型矢量,θ∈{θ1,θ2,…,θD},θ为声源可能会出现的方位;
计算:y(k,l)=[Y1(k,l) Y2(k,l)…YM(k,l)]T,其中,M为所述麦克风阵列中麦克风的个数;
寻找出现最大峰值点的角度,如下式:
即为在第l帧、第k个频点上的单声源波达方位估计结果;
计算获得全部时间内,多个所述频点上的全部单声源波达方位估计结果,获得所述估计结果。
4.根据权利要求3所述的方法,其特征在于,所述统计获得预设时间周期内不同声源方位被估计获得的频率,包括:
根据所述估计结果,获得全部时间内,多个所述频点上的全部所述单声源波达方位估计结果;
基于全部所述单声源波达方位估计结果,统计获得全部单声源波达方位被估计获得的频率Q(θ),其中,θ∈{θ1,θ2,…,θD}。
5.根据权利要求1所述的方法,其特征在于,所述查找获得所述统计频率结果Q(θ)内的峰值点,获取被估计获得的频率大于门限值τ1的多个声源方位,包括:
查找获取所述统计频率结果Q(θ)内所有的峰值点;
选取大于所述门限值τ1的全部所述峰值点,记作θi1,θi2,…,θiI;
构建获得第一向量θcand,θcand长度为D,如下:
将所述第一向量作为所述声源数量估计结果。
6.根据权利要求5所述的方法,其特征在于,所述更新当前活跃的声源个数,包括:
采用激活、保持和中止策略进行更新,更新当前活跃的声源个数。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西北工业大学,未经西北工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202210445311.3/1.html,转载请声明来源钻瓜专利网。
- 上一篇:一种曲面铝板幕墙系统
- 下一篇:半导体结构及其制备方法





