一种用于噪声环境下声源分离和定位方法

专利类型：发明专利 

语 言：中文 

申 请 号：CN201410249428.X 

申 请 日：20140606 

发 明 人：柴毅屈剑锋王坤朋杨志敏董勇 

申 请 人：重庆大学 

申请人地址：400044 重庆市沙坪坝区沙正街174号重庆大学自动化学院 

公 开 日：20180605 

公 开 号：CN104053107B 

代 理 人：赵荣之 

代理机构：北京同恒源知识产权代理有限公司 11275 

摘　　要：本发明涉及助听装置及方法，具体指一种用于噪声环境下声源分离和定位的助听装置及方法，包括麦克风阵列、声源位置显示器、声源分离与定位模块、信号采集器、模拟信号放大器和声源选择键盘，声源分离与定位模块采用互相关法对采集的麦克风信号进行处理，得到8个声源相对于麦克风阵列的初始方位，采用基于空间搜索的盲源分离方法对采集的麦克风信号进行处理，得到包含交谈对象和其他人说话声音的声源，用户通过声源选择键盘选择出属于交谈对象的声源，根据选择的声源计算出交谈对象相对于麦克风阵列的准确方位，并在声源位置显示器中显示该声源的方位；实现噪声环境下交谈双方语音信号自动分离，同时增加了语音信号源的定位功能，使用方便。 

主 权 项：一种用于噪声环境下声源分离和定位方法，其特征在于：包括如下步骤：a.采用助听装置对交谈过程中的语音信号进行采集；所述助听装置包括麦克风阵列(1006)、声源位置显示器(1001)、声源分离与定位模块(1002)、信号采集器(1003)、模拟信号放大器(1004)和声源选择键盘(1005)；所述麦克风阵列(1006)、模拟信号放大器(1004)、信号采集器(1003)、声源分离与定位模块(1002)和声源位置显示器(1001)依次连接；所述声源选择键盘(1005)与声源分离与定位模块(1002)连接；所述声源分离与定位模块(1002)采用互相关法对采集的麦克风信号进行处理，得到8个声源相对于麦克风阵列的初始方位；采用基于空间搜索的盲源分离方法对采集的麦克风信号进行处理，得到包含交谈对象(1009)和其他人(1010)说话声音的声源，用户通过声源选择键盘(1005)选择出属于交谈对象(1009)的声源，根据选择的声源计算出交谈对象(1009)相对于麦克风阵列(1006)的准确方位，并在声源位置显示器(1001)中显示该声源的方位；所述麦克风阵列(1006)由微型麦克风(1007)构成：所述麦克风阵列(1006)由8个独立的微型麦克风(1007)均匀排列在正方形边上，正方形边长为10厘米，任意两个微型麦克风(1007)中心点之间的距离(1008)为5厘米；所述微型麦克风(1007)的指向性为全向型，无指向，直径小于1厘米；a1.所述助听装置采用模拟信号放大器(1004)对麦克风阵列(1006)中8个微型麦克风(1007)的语音信号进行放大，放大后的信号传输至信号采集器(1003)；a2.所述信号采集器(1003)对模拟信号放大器(1004)放大后的信号进行滤波、模/数转换和采集，采样频率为16KHz，采样精度为16位，并对采集的信号进行储存；b.所述声源分离与定位模块(1002)采用基于空间搜索的盲源分离方法对采集的麦克风信号进行处理，具体步骤如下：b1.信号采集器得到8路观测信号x(t)＝[x₁(t),x₂(t),…,x_s(t)]^T，t为采样时刻；b2.观测信号x(t)被建模为8个声源s(t)＝[s₁(t),s₂(t),…,s_s(t)]^T的线性混合：$\begin{array}{c}{x}_1\left(t\right)=\frac{1}{a_{11}}·s_1\left(t-\frac{a_{11}}{v}\right)+\frac{1}{a_{12}}·s_2\left(t-\frac{a_{12}}{v}\right)+\mathrm{...}+\frac{1}{a_{18}}·s_8\left(t-\frac{a_{18}}{v}\right)+{\eta }_1\left(t\right)\\ x_2\left(t\right)=\frac{1}{a_{21}}·s_1\left(t-\frac{a_{21}}{v}\right)+\frac{1}{a_{22}}·s_2\left(t-\frac{a_{22}}{v}\right)+\mathrm{...}+\frac{1}{a_{28}}·s_8\left(t-\frac{a_{28}}{v}\right)+{\eta }_2\left(t\right)\\ \begin{array}{c}.\\ .\\ .\end{array}\\ x_8\left(t\right)=\frac{1}{a_{81}}·s_1\left(t-\frac{a_{81}}{v}\right)+\frac{1}{a_{82}}·s_2\left(t-\frac{a_{82}}{v}\right)+\mathrm{...}+\frac{1}{a_{88}}·s_8\left(t-\frac{a_{88}}{v}\right)+{\eta }_8\left(t\right)\end{array}---\left(1\right)$式(1)中v为空气中声音传播的速度，约为340米/秒；a_mn，m,n＝1,2,…,8为第n个声源到达第m个微型麦克风的距离，声源信号的幅值衰减与距离a_mn成反比关系，声源s_n，n＝1,2,…,8到达各微型麦克风的时间延迟与距离a_mn成正比关系：η(t)＝(η₁(t),η₂(t),…,η₈(t))^T为微型麦克风的噪声；公式(1)简写为x(t)＝As(t)+η(t)(2)式(2)中A为8×8声源混合矩阵；b3.采用互相关法对采集的麦克风信号进行处理，计算得到声源相对于麦克风阵列的初始方位，具体如下：b301：为了降低声源定位和分离算法的复杂度，构建一个栅格化的正方体，并存储每个栅格顶点至微型麦克风(2003)的距离；将麦克风阵列正前方的正方体区域进行栅格化处理，正方体的边长为10厘米，栅格的边长为1厘米，最小栅格的边长为1厘米，记栅格顶点(2001)为g_ijk，i,j,k＝1,2,…,10；b302：栅格顶点g_ijk距离第m个微型麦克风(2003)的距离(2002)为d_m，将g_ijk至8个微型麦克风(2003)的距离写作列向量形式为d，将所有栅格顶点至微型麦克风的距离写作矩阵形式为D，矩阵的第l列表示第(k?1)×100+(j?1)×10+i个栅格点；b303：由公式(1)知，任意两路观测信号x_m(t)与x_m'(t)，m,m'∈{1,2,…,8}进行延迟相关运算C_τ＝x_m(t)·x_m'(t?τ)，τ∈[?0.414/v,0.414/v]，τ为延迟时间，必然存在8个极大值n＝1,2,…,8；由于τ非采样时间间隔的整倍数，在计算时需要首先对观测信号x(t)进行自然三次样条插值，对于采样时刻t₀,t₁,…,t_p，对应的采样点x_m＝[y₀,y₁,…,y_P]^T，P为采样点数目，插值计算后得到的观测信号为插值函数为：$z_p\left(t\right)=\frac{b_{p+1}·{(t-t_p)}^3+b_p·{(t_{p+1}-t)}^3}{6h_p}+(\frac{y_{p+1}}{h_p}-\frac{h_p}{6}{b}_{p+1})(t-t_p)+(\frac{y_p}{h_p}-\frac{h_p}{6}{b}_p)(t_{p+1}-t)---\left(3\right)$h_p＝t_p+1?t，0≤p≤P系数b₀,b₁,…,b_P通过求解式(4)得到；$\left\{\begin{array}{c}{b}_0=0\\ h_{p-1}{b}_{p-1}+2(h_{p-1}+h_p)+h_p{b}_{p+1}=6\left(\right(y_{p+1}-y_p)/h_p-(y_p-y_{p-1})/h_{p-1})\\ b_P=0\end{array}\right.---\left(4\right)$b304：分别计算微型麦克风(3001)与(3002)、(3003)、(3004)采集到观测信号之间的延迟相关系数C_τ的极大值，并利用极大值处对应的延迟时间和计算出声源相对于微型麦克风的方位；b305：利用延迟时间τ与声速v相乘分别计算8个声源与4个微型麦克风(3001)、(3002)、(3003)、(3004)的距离d₁,d₃,d₆,d₈，并通过计算求取与声源距离最近的栅格点g_ijk，进而得到8个声源距离8个微型麦克风的距离，即声源混合矩阵A；b4.采用空间搜索的方法以最小化声源之间互信息I为目标函数，同时分离出混合在观测信号x(t)中的所有声源n＝1,2,…,8，由于语音信号为短时平稳信号，通常认为在20毫秒的采样时长内等效为平稳信号，因此，以下过程需要在有限时间T内完成，T≤20ms，具体如下：b401：初始化空间搜索次数r←0，声源混合矩阵A⁽⁰⁾←A，互信息最小值I^*←+∞，时间tc←0；b402：对声源混合矩阵A^(r)进行求逆运算得到矩阵W^(r)，利用计算得到真实声源信号的估计信号并对其进行标准化处理，使得各声源估计信号均值为0方差为1；b403：采用互信息度量声源信号之间的相互独立性，声源之间的互信息是一个非负指标，当且仅当各信号之间相互独立时互信息为零；分别记x_m＝(x_m(1),…,x_m(P))为第m个观测信号，为真实声源信号的估计，P是采样点数目；由此，估计声源信号之间的互信息写为：$I^{\left(r\right)}=\underset{n=1}{\overset{N}{\Sigma }}H\left({\hat{s}}_n\right)-H\left(\hat{s}\right)---\left(5\right)$式中，为估计声源信号的边缘熵，为：$H\left(\hat{s}\right)=H\left(x\right)+log\left|\det A^{\left(r\right)}\right|---\left(7\right)$若I^(r)＜I^*，则置I^*←I^(r)，否则，进入下一步；b404：随机搜索8个声源附近的栅格点，得到一组新的栅格点及其对应的到达微型麦克风的距离，通过查询距离矩阵D，得到新的声源混合矩阵A；b405：若计算时间tc＜T，置迭代次数r←r+1，置下一次迭代的混合矩阵为A^(r)←A，返回b402进行迭代计算，否则，结束迭代过程。 

关 键 词：声源;麦克风阵列;交谈对象;噪声环境;麦克风信号;定位模块;声源位置;助听装置;显示器;模拟信号放大器;采集;信号采集器;语音信号源;定位功能;键盘选择;空间搜索;盲源分离;说话声音;语音信号;准确方位;自动分离;互相关;源选择;和声;键盘 

法律状态：授权 

IPC专利分类号：H04R25/00(2006.01)I,H04R27/02(2006.01)I,G10L21/0272(2013.01)I