[发明专利]处理声音信号的方法及系统有效
申请号: | 201811645765.5 | 申请日: | 2018-12-29 |
公开(公告)号: | CN109727605B | 公开(公告)日: | 2020-06-12 |
发明(设计)人: | 袁斌 | 申请(专利权)人: | 苏州思必驰信息科技有限公司 |
主分类号: | G10L21/0232 | 分类号: | G10L21/0232;G10L21/0208 |
代理公司: | 北京商专永信知识产权代理事务所(普通合伙) 11400 | 代理人: | 方挺;车江华 |
地址: | 215123 江苏省苏州市苏*** | 国省代码: | 江苏;32 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 处理 声音 信号 方法 系统 | ||
1.一种处理声音信号的方法,其特征在于,包括:
获取待处理声音信号,所述待处理声音信号包括目标声音信号和干扰声音信号,所述干扰声音信号包括噪声信号和回声信号;
确定所述干扰声音信号的功率谱密度;
将所述待处理声音信号转换为频域信号E(Ω);
根据以下公式确定后验信噪比PostSNR(Ω):
PostSNR(Ω)=|E(Ω)|2/(Rbb(Ω)+Rnn(Ω)),
其中,Rbb(Ω)为所述回声信号的功率谱密度,Rnn(Ω)为所述噪声信号的功率谱密度;
根据以下公式推导出先验信噪比PrioriSNR(Ω):
PrioriSNR(Ωi)=(1-alpha)*P(PostSNR(Ωi)-1)+alpha*|S’(Ωi-1)|2/Rbb(Ω);
其中,alpha为平滑因子,P(x)=(|x|+x)/2,S’(Ωi-1)为上一帧声音信号的频谱估计;
进一步计算加权系数HLSA(Ω),并得到所述目标声音信号的频谱估计S’(Ω):
S’(Ω)=E(Ω)*HLSA(Ω),
其中,theta=PostSNR(Ω)*PrioriSNR(Ω)/(PrioriSNR(Ω)+1);
根据所述频谱估计确定掩蔽阈值;
确定所述待处理声音信号中干扰声音信号的频谱成分大于所述掩蔽阈值的情况下,对所述待处理声音信号进行滤波处理。
2.根据权利要求1所述的方法,其特征在于,确定所述待处理声音信号中干扰声音信号的频谱成分大于所述掩蔽阈值的情况下,对所述待处理声音信号进行滤波处理的步骤包括:
根据回声信号的功率谱密度和噪声信号的功率谱密度确定出滤波处理的加权系数H(Ω):
H(Ω)=min(1,sqrt(RTT(Ω)/(Rbb(Ω)+Rnn(Ω)))+(zeta_b*Rbb(Ω)+zeta_n*Rnn(Ω))/(Rbb(Ω)+Rnn(Ω))),
其中,Rbb(Ω)为所述回声信号的功率谱密度,Rnn(Ω)为所述噪声信号的功率谱密度,zeta_b为回声衰减系数,zeta_n为噪声衰减系数。
3.根据权利要求1所述的方法,其特征在于,根据所述频谱估计确定掩蔽阈值的步骤包括:
根据频谱估计,确定所述待处理声音信号的临界频带的功率谱密度B(k)和扩展临界频带频谱C(k):
C(k)=B(k)*SF(k),
其中,SF(k)=15.81+7.5*k+0.474-17.5*sqrt(1+(k+0.474)2),bh,bl分别为各临界频带的上下限频率;
根据扩展临界频带频谱C(k)和偏移函数O(k),确定初步掩蔽阈值T(k):
T(k)=10lg(C(k))-(O(k)/10),
其中,偏移函数O(k)=belta*(14.5+k)+(1-belta)*5.5;belta为音调系数;
根据初步掩蔽阈值T(k)和绝对听阈Tabs(k),确定掩蔽阈值RTT(Ω):
RTT(Ω)=min(T(k),Tabs(k)),
其中,Tabs(k)=3.64f-0.8-6.5exp(f-3.3)2+10-3f4。
4.根据权利要求1所述的方法,其特征在于,获取待处理声音信号的步骤包括:
接收初始声音信号;
对所述初始声音信号进行回声消除,以得到所述待处理声音信号。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于苏州思必驰信息科技有限公司,未经苏州思必驰信息科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201811645765.5/1.html,转载请声明来源钻瓜专利网。