一种移动源信号混叠的盲分离方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN110111806A(43)申请公布日2019.08.09(21)申请号201910234231.1(22)申请日2019.03.26(71)申请人广东工业大学地址510062广东省广州市大学城外环西路100号(72)发明人解元谢胜利谢侃吴宗泽(74)专利代理机构广东广信君达律师事务所44329代理人杜鹏飞杨晓松(51)Int.Cl.G10L21/0216(2013.01)G10L21/0232(2013.01)G10L21/0272(2013.01)G10L25/69(2013.01)权利要求书3页说明书7页附图4页(54)发明名称一种移动源信号混叠的盲分离方法(57)摘要本发明公开了一种移动源信号混叠的盲分离方法，由于在现实生活中，采集的混叠信号具有一定的可变性，它的混叠滤波器是时变的，与传统的时不变混叠滤波器相比，源信号的分离过程将变得更加困难，本发明针对这种时变的混叠情况，提出一种鲁棒的盲分离算法；首先利用时间差估计算法定位源信号的位置，给出了严格的数学理论推导；然后对混叠滤波器进行重构，利用期望最大化算法更新模型参数，再利用维纳滤波法分离源信号；最后，通过仿真实验验证所提算法的有效性，同时对比一种鲁棒的移动源信号混叠盲分离算法(Full-rank算法)，证明本发明所提算法对移动的源信号混叠盲分离具有更好的鲁棒性。CN110111806ACN110111806A权利要求书1/3页1.一种移动源信号混叠的盲分离方法，其特征在于，包括下述步骤：步骤一，对移动的声音源信号混叠进行数学建模；假设利用J个固定的麦克风去接收I个移动的声音源信号，产生的混叠信号的数学表达式为：其中，xj(t)是第j个麦克风接收到的混叠信号(j＝1,...,J)，hijt(τ)是时变的空间脉冲响应，τ是时间延迟；为了在频域上进行源信号的分离，对上述公式(1)进行短时傅里叶变换，则时变的空间脉冲响应hijt(τ)的混叠性质随着时间的变化是缓慢的，因此可以认为在一个很短的时间窗内混叠过程是时不变；所以，时变的混叠模型公式(1)可以近似于频域上的混叠模型如下述公式(2)所示：T其中，xfn＝[xfn1,...,xfnJ]表示混叠信号的短时傅里叶变换，sfn,i表示第i个源信号的T单通道短时傅里叶变换，hfn,i＝[hfn1,...,hfnJ]是频域上的空间脉冲响应；步骤二，时间延迟τ的估计；给定一组麦克风m和n，它们在笛卡尔坐标系的位置分别为m33和n，即m∈R,n∈R，k0指代两个麦克风的中间位置p到声音源信号方向上的单位向量，即k03T∈R,||k0||＝1，d是两个麦克风之间的距离；同时定义p＝[0,0,0]是笛卡尔坐标系上的原点，θ指代仰角，且-90°≤θ≤90°；利用向量乘积，可得：其中，||n||表示向量n的范数，τn(k0)表示声音传输到麦克风n和位置p时产生的时间延迟，v是声音在空气中的传播速度，取v＝340m/s；根据上式可得：因此，麦克风m和n之间的时间延迟为：其中，m＝1,...,J,n＝1,...,J；步骤三，混叠滤波器的重构；基于步骤二估计的每组麦克风之间的延迟τ(m,n)，对应于频域上的相位差是：2CN110111806A权利要求书2/3页hmn,i＝exp(-jωfτ(m,n))，其中，ωf＝2π(f-1)Fs/N，Fs是样本频率，N是短时傅里叶变换窗长度；步骤四，源信号的分离；首先，定义源信号的功率谱密度的非负矩阵分解如下：其中，是正实值幅度谱图，wfk表示单个非负矩阵分解分量k的振幅谱，hkn表示每一帧n的分量增益；然后，将上述定义的非负矩阵分解源信号模型代入混叠模型公式(2)中，可得：选用平方Frobenius范数作为代价函数，其表达式如下：然后，利用期望最大化算法即EM算法迭代更新模型参数Hfn,wfk,hkn；具体细节如下：E-step：自然统计的条件期望：Rc,fn＝diag([wfkhkn]k)M-step：更新模型参数的公式：其中，而且cfn＝TK[c1,fn,…,cK,fn]∈C表示成分分量系数的向量，每个分量ck,fn遵循适当的多元复高斯分布，即，ck,fn～Nc(0,wfkhkn)；3CN110111806A权利要求书3/3页同时，利用维纳滤波法从混叠信号中分离源信号，得到频域上的估计源信号：最后，利用傅里叶变换的逆运算对估计的源信号进行重构，得到时域上的源信号。4CN110111806A说明书1/7页一种移动源信号混叠的盲分离方法技术领域[0001]本发明涉及盲信号分离技术领域，具体涉及一种移动源信号混叠的盲分离方法。背景技术[0002]在经典的鸡尾酒会问题中，由于聚会上的人处于移动状态，采集的混叠声音比较复杂。如何仅仅从接收