CN104811867B - 基于阵列虚拟扩展的麦克风阵列空域滤波方法 - Google Patents
基于阵列虚拟扩展的麦克风阵列空域滤波方法 Download PDFInfo
- Publication number
- CN104811867B CN104811867B CN201510212768.XA CN201510212768A CN104811867B CN 104811867 B CN104811867 B CN 104811867B CN 201510212768 A CN201510212768 A CN 201510212768A CN 104811867 B CN104811867 B CN 104811867B
- Authority
- CN
- China
- Prior art keywords
- array
- virtual extended
- matrix
- data
- mrow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
Abstract
基于阵列虚拟扩展的麦克风阵列空域滤波方法,步骤为:获取频域分段麦克风阵列数据;计算语音信号处理的频率单元序列;组成与频率单元序列对应的频域麦克风阵列数据;构造目标语音信号的阵列导向矩阵及虚拟扩展阵列导向矩阵;估计阵列虚拟扩展矩阵;对频率单元序列对应的阵列数据进行虚拟扩展变换;对虚拟扩展阵列对应的阵列数据进行稳健自适应波束形成;将虚拟扩展阵列输出按照频域分段麦克风阵列数据格式进行组合,得到频域分段麦克风阵列输出数据;对频域分段麦克风阵列输出数据的每一子段数据进行逆傅里叶变换,得到时域分段阵列处理语音信号;进行输出语音信号合成,将相邻子段数据重叠部分进行平均,得到最终的麦克风阵列输出语音信号。
Description
技术领域
本发明属于信号处理技术领域,尤其涉及一种语音信号的空域滤波方法,可以降低空间噪声和其它话音干扰的影响,进而增强有用语音信号的信噪比,可用于通信、雷达和声呐等宽带信号的空域滤波等领域。
背景技术
在语音信号处理中,麦克风阵列作为提高语音质量的关键设备,主要用于噪声抑制、回声消除以及语音信号选择等语音信号增强方面,同时也用于声源的方位确定和空间位置估计。麦克风阵列由于具有空间选择性,因而可以通过空域处理来提高语音获取和处理的性能。然而由于语音信号的特殊性,也给基于麦克风阵列的语音信号处理带来了困难和挑战。对于麦克风阵列处理,由于语音信号基本占据了所有处理频带(0~8KHz),对于给定的麦克风阵列配置来说,阵元数和阵元间隔是确定的,对于高频段的语音分量,一般具有较好的处理效果,然而对于低频分量,由于信号波长趋于无穷大,因此相应的阵列处理性能将随着频率分量所处频带的降低而性能将急剧下降。
基于麦克风阵列的语音增强技术,决定其处理效果的核心是自适应波束形成器的实现,这是由语音信号的宽带特征所决定的。基于传统麦克风阵列处理的空域滤波很难获得理想的语音增强效果,这是因为低频段的主瓣波束很宽,几乎没有任何抑制空间噪声和回声效果。针对麦克风阵列语音增强处理的特殊需求,有必要提出可以改善空域滤波性能的方法,以改善基于麦克风阵列处理的噪声、干扰和回声的抑制性能。
发明内容
本发明的目的是提供一种用于麦克风阵列语音增强的空域滤波方法,解决麦克风阵列空域滤波时的全频段一致性,以提高麦克风阵列的空域滤波性能及语音增强性能。
为了实现上述目的,本发明采取如下的技术解决方案:
基于阵列虚拟扩展的麦克风阵列空域滤波方法,包括以下步骤:
步骤1、获取频域分段麦克风阵列数据;
对麦克风阵列的采样数据进行分段,并对每个通道的每一子段都分别进行离散傅里叶变换,得到频域分段麦克风阵列数据;
步骤2、计算语音信号处理的频率单元序列;
根据阵列数据的采样频率和子段数据长度计算麦克风阵列处理频段中对应的频率单元矢量,去掉零频率后得到语音信号处理的频率单元序列;
步骤3、组成与频率单元序列对应的频域麦克风阵列数据;
针对每一个频率单元,从频域分段麦克风阵列数据的每一子段数据中抽取对应的列数据,组成与频率单元序列对应的频域麦克风阵列数据
步骤4、构造目标语音信号的阵列导向矩阵A及虚拟扩展阵列导向矩阵
步骤5、根据目标语音信号的阵列导向矩阵A和虚拟扩展阵列导向矩阵估计阵列虚拟扩展矩阵P;
步骤6、根据阵列虚拟扩展矩阵P,对频率单元序列对应的频域麦克风阵列数据进行虚拟扩展变换,得到虚拟扩展阵列对应的阵列数据
步骤7、对虚拟扩展阵列对应的阵列数据进行稳健自适应波束形成,得到虚拟扩展阵列输出;
步骤7-1、根据虚拟扩展阵列对应的阵列数据估计样本协方差矩阵式中的Nseg为数据分段数;
步骤7-2、构造线性约束最小功率波束形成器的约束矩阵C和约束值矢量f;
步骤7-3、估计稳健LCMP波束形成器的加权矢量w;
其中,μ为加载电平,I为单位矩阵;
步骤7-4、根据加权矢量对虚拟扩展阵列对应的阵列数据进行加权处理,得虚拟扩展阵列输出
步骤8、将虚拟扩展阵列输出按照频域分段麦克风阵列的数据格式进行组合,得到频域分段麦克风阵列输出数据;
步骤9、对频域分段麦克风阵列输出数据的每一子段数据进行逆傅里叶变换,得到时域分段阵列处理语音信号;
步骤10、进行输出语音信号合成,将时域分段阵列处理语音信号中相邻子段数据重叠部分进行平均,得到最终的麦克风阵列输出语音信号。
更进一步的,所述步骤5中,当目标语音信号数量大于阵元数N时,根据下式估计阵列虚拟扩展矩阵P:
当目标语音信号数量小于阵元数N时,根据如下步骤估计阵列虚拟扩展矩阵P:
构造矩阵并对矩阵G进行奇异值分解:G=UΣVH,其中,U和V分别为以G的左右奇异矢量形成的N×N维和(N+ΔN)×(N+ΔN)维的酉矩阵,Σ为N×(N+ΔN)维的奇异特征值矩阵,ΔN为虚拟扩展的阵元个数;
计算阵列虚拟扩展矩阵其中,E=[IN 0N×ΔN]T,IN为N阶单位矩阵,0N×ΔN为N×ΔN维零矩阵,为阵列虚拟扩展增益。
更进一步的,所述步骤7-2中,约束矩阵C的构造方法是将目标语音信号中有用信号的导向矢量放在第一列,将其它干扰信号的导向矢量放在后续列;对应约束矩阵C的列将目标语音信号中有用信号的对应约束值设置为1,干扰信号的对应约束值设置为0,得到约束值矢量f。
本发明通过阵列虚拟扩展技术降低针对有用目标语音信号的主瓣宽度,同时降低旁瓣电平以减少反射回声和空间噪声的接收,进而提高有用目标信号的空间选择性,通过在低频段进行更大阵列孔径的虚拟扩展来保证低频段可以获得与高频端近似一致的空域滤波效果,进而保证全频段空域滤波的一致性,以提高基于麦克风阵列的语音增强效果。本发明实现简单,在进行空域滤波时直接进行阵列虚拟扩展处理,可用于麦克风阵列语音增强,也可用于其它宽带阵列自适应滤波处理领域。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中需要使用的附图做简单介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法的流程图。
具体实施方式
为了让本发明的上述和其它目的、特征及优点能更明显,下文特举本发明实施例,并配合所附图示,做详细说明如下。
本发明针对麦克风阵列语音增强处理的特殊需求,在阵列虚拟扩展基础上改善基于麦克风阵列处理的噪声、干扰和回声抑制性能,通过一种麦克风阵列虚拟扩展方法,以提高低频段的波束指向性能,进而获得较窄的主瓣波束宽度和较低的旁瓣电平,以降低空间噪声、干扰和回声信号的接收。同时实现麦克风阵列在语音处理全频段的空域滤波一致性。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其它不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制
参照图1,图1为本发明方法的流程图。本发明方法的步骤如下:由N个阵元组成的麦克风阵列接收空间远场声源信号,N为阵列的阵元数,
步骤1、获取频域分段麦克风阵列数据
对麦克风阵列的采样数据进行分段,并对每个通道的每一子段都分别进行离散傅里叶变换,得到频域分段麦克风阵列数据
其中,SubK为子段数据长度,子段数据长度为偶数,Nseg为数据分段数,Olap为相邻子段数据重叠长度,0≤Olap≤SubK-1,数据分段数TotalK为语音数据处理长度;
频域分段麦克分阵列数据中的每个通道的每一子段都由对应时域分段采样数据进行离散傅里叶变换得到,在频域分段麦克风阵列数据中,第l(l=1,…,Nseg)子段的表达式为:
步骤2、计算语音信号处理的频率单元序列;
根据阵列数据的采样频率fs和子段数据长度SubK计算麦克风阵列处理频段中的对应频率单元矢量由于阵列数据经过FFT变换后,而且SubK选择为偶数,故零频对应的阵列数据是不存在的,而且零频率对应的波长无穷大,也不满足阵列处理条件,因此,去掉零频率后得到语音信号处理的频率单元序列:
步骤3、组成与频率单元序列对应的频域麦克风阵列数据;
针对每一个频率单元,从频域分段麦克风阵列数据的每一子段数据中抽取对应的列数据,组成与频率单元序列对应的频域麦克风阵列数据
步骤4、构造目标语音信号的阵列导向矩阵A以及虚拟扩展阵列导向矩阵
根据目标语音信号的波达方向估计值、麦克风阵列配置参数以及虚拟扩展阵列参数列表构造阵列导向矩阵及虚拟扩展阵列导向矩阵,虚拟扩展阵列参数列表是根据阵列配置参数以及事先确定的空域滤波效果给出每个语音信号处理频率单元所对应的虚拟扩展阵元数量和阵元间距参数,该列表可以在麦克风阵列处理数据之前离线完成;
以均匀线阵为例,阵列导向矩阵
其中,为第p个语音信号的导向矢量,ωp为第p个语音信号的角频率,例如对于本发明,当进行频率单元fk所对应的阵列数据处理时,由于该频率单元的所有目标信号(包括有用信号和干扰信号)的中心频率都为fk,故语音信号的角频率ωp=2πfk,θp为第p个语音信号的方位角,d为阵元间距,c为语音信号的传播速度,p=1,…,P,n=1,…,N,P为语音信号的数量;
虚拟扩展阵列导向矩阵的构造方法同A,不同之处是阵列参数不同,如阵元间距、阵元数量等;
当阵列结构为其他类型时,按照相应阵列导向矢量公式进行构造;
步骤5、根据目标语音信号的阵列导向矩阵A和虚拟扩展阵列导向矩阵A~估计阵列虚拟扩展矩阵P;
步骤5-1、当目标语音信号数量大于阵元数N时,根据下式估计阵列虚拟扩展矩阵P:
步骤5-2、当目标语音信号数量小于阵元数N时,根据如下步骤估计阵列虚拟扩展矩阵P:
构造矩阵对矩阵G进行奇异值分解:G=UΣVH;
其中,U和V分别为以G的左右奇异矢量形成的N×N维和(N+ΔN)×(N+ΔN)维的酉矩阵,Σ为N×(N+ΔN)维的奇异特征值矩阵,由其前N列元素组成的分块矩阵ΣN为对角矩阵,而其余元素为零,ΔN为虚拟扩展的阵元个数;
阵列虚拟扩展矩阵
其中,E=[IN 0N×ΔN]T,IN为N阶单位矩阵,0N×ΔN为N×ΔN维零矩阵,为阵列虚拟扩展增益;
步骤6、根据阵列虚拟扩展矩阵P,对频率单元序列对应的阵列数据进行虚拟扩展变换,得到虚拟扩展阵列对应的阵列数据
步骤7、对虚拟扩展阵列对应的阵列数据进行稳健自适应波束形成;
步骤7-1、根据虚拟扩展阵列对应的阵列数据估计样本协方差矩阵
步骤7-2、构造线性约束最小功率(LCMP)波束形成器的约束矩阵C和约束值矢量f;
约束矩阵C的典型构造方法是将目标语音信号中有用信号的导向矢量放在第一列,将其它干扰信号的导向矢量放在后续列;对应约束矩阵C的列将目标语音信号中有用信号的对应约束值设置为1,干扰信号的对应约束值设置为0,得到约束值矢量f;根据场景信息和测向结果进行判断信号为有用信号或干扰信号;
步骤7-3、估计稳健LCMP波束形成器的加权矢量w;
其中,μ为加载电平,μ>0,加载电平可根据经验值进行选择,选择大于最大特征值的模值即可获得好的性能,I为单位矩阵;
步骤7-4、根据加权矢量对虚拟扩展阵列对应的阵列数据进行加权处理,得虚拟扩展阵列输出
步骤8、将虚拟扩展阵列输出按照频域分段麦克风阵列数据格式进行组合,得到频域分段麦克风阵列输出数据
步骤9、对频域分段麦克风阵列输出数据的每一子段数据进行逆傅里叶变换,得到时域分段阵列处理语音信号SegY1×(Nseg×SubK):
表示将第i(i=1,2,…,N)通道的第j(j=1,2,…Nseg)子段的频域数据变换到时域;
步骤10、进行输出语音信号合成,将相邻子段数据重叠部分进行平均,得到最终的麦克风阵列输出语音信号。
与现有技术相比,本发明具有以下优点:
1、利用奇异值分解法求解阵列虚拟扩展矩阵,使得阵列虚拟扩展的误差更小;
2、利用阵列虚拟扩展技术提高麦克风阵列语音处理低频段的空域滤波性能;
3、通过对不同频率单元利用不同的阵列虚拟扩展策略以保证麦克风阵列全频段空域滤波的一致性;
4、由于在阵列虚拟扩展变换中,干扰和噪声的协方差矩阵存在误差,通过稳健自适应波束形成算法抑制阵列虚拟扩展变换中的协方差矩阵失配,有效提高了麦克风阵列处理性能。
以上所述,仅是本发明的较佳实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (3)
1.基于阵列虚拟扩展的麦克风阵列空域滤波方法,其特征在于,包括以下步骤:
步骤1、获取频域分段麦克风阵列数据;
对麦克风阵列的采样数据进行分段,并对每个通道的每一子段都分别进行离散傅里叶变换,得到频域分段麦克风阵列数据;
步骤2、计算语音信号处理的频率单元序列;
根据阵列数据的采样频率和子段数据长度计算麦克风阵列处理频段中对应的频率单元矢量,去掉零频率后得到语音信号处理的频率单元序列;
步骤3、组成与频率单元序列对应的频域麦克风阵列数据;
针对每一个频率单元,从频域分段麦克风阵列数据的每一子段数据中抽取对应的列数据,组成与频率单元序列对应的频域麦克风阵列数据fk为频率单元,k=1,2,…,SubK,SubK为子段数据长度;
步骤4、构造目标语音信号的阵列导向矩阵A及虚拟扩展阵列导向矩阵
步骤5、根据目标语音信号的阵列导向矩阵A和虚拟扩展阵列导向矩阵估计阵列虚拟扩展矩阵P;
步骤6、根据阵列虚拟扩展矩阵P,对频率单元序列对应的频域麦克风阵列数据进行虚拟扩展变换,得到虚拟扩展阵列对应的阵列数据
步骤7、对虚拟扩展阵列对应的阵列数据进行稳健自适应波束形成,得到虚拟扩展阵列输出;
步骤7-1、根据虚拟扩展阵列对应的阵列数据估计样本协方差矩阵式中的Nseg为数据分段数;
步骤7-2、构造线性约束最小功率波束形成器的约束矩阵C和约束值矢量f;
步骤7-3、估计稳健LCMP波束形成器的加权矢量w;
<mrow>
<mi>w</mi>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>&mu;</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>C</mi>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>C</mi>
<mi>H</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mover>
<mover>
<mi>R</mi>
<mo>~</mo>
</mover>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>&mu;</mi>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>f</mi>
<mo>;</mo>
</mrow>
其中,μ为加载电平,I为单位矩阵;
步骤7-4、根据加权矢量对虚拟扩展阵列对应的阵列数据进行加权处理,得虚拟扩展阵列输出Olap为相邻子段数据重叠长度;
步骤8、将虚拟扩展阵列输出按照频域分段麦克风阵列的数据格式进行组合,得到频域分段麦克风阵列输出数据;
步骤9、对频域分段麦克风阵列输出数据的每一子段数据进行逆傅里叶变换,得到时域分段阵列处理语音信号;
步骤10、进行输出语音信号合成,将时域分段阵列处理语音信号中相邻子段数据重叠部分进行平均,得到最终的麦克风阵列输出语音信号。
2.根据权利要求1所述的基于阵列虚拟扩展的麦克风阵列空域滤波方法,其特征在于:所述步骤5中,当目标语音信号数量大于阵元数N时,根据下式估计阵列虚拟扩展矩阵P:
<mrow>
<mi>P</mi>
<mo>=</mo>
<mover>
<mi>A</mi>
<mo>~</mo>
</mover>
<msup>
<mi>A</mi>
<mi>H</mi>
</msup>
<msup>
<mrow>
<mo>(</mo>
<msup>
<mi>AA</mi>
<mi>H</mi>
</msup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>;</mo>
</mrow>
当目标语音信号数量小于阵元数N时,根据如下步骤估计阵列虚拟扩展矩阵P:
构造矩阵并对矩阵G进行奇异值分解:G=UΣVH,其中,U和V分别为以G的左右奇异矢量形成的N×N维和(N+ΔN)×(N+ΔN)维的酉矩阵,Σ为N×(N+ΔN)维的奇异特征值矩阵,ΔN为虚拟扩展的阵元个数;
计算阵列虚拟扩展矩阵其中,E=[IN 0N×ΔN]T,IN为N阶单位矩阵,0N×ΔN为N×ΔN维零矩阵,为阵列虚拟扩展增益。
3.根据权利要求1所述的基于阵列虚拟扩展的麦克风阵列空域滤波方法,其特征在于:所述步骤7-2中,约束矩阵C的构造方法是将目标语音信号中有用信号的导向矢量放在第一列,将其它干扰信号的导向矢量放在后续列;对应约束矩阵C的列将目标语音信号中有用信号的对应约束值设置为1,干扰信号的对应约束值设置为0,得到约束值矢量f。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510212768.XA CN104811867B (zh) | 2015-04-29 | 2015-04-29 | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510212768.XA CN104811867B (zh) | 2015-04-29 | 2015-04-29 | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104811867A CN104811867A (zh) | 2015-07-29 |
CN104811867B true CN104811867B (zh) | 2017-11-21 |
Family
ID=53696276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510212768.XA Expired - Fee Related CN104811867B (zh) | 2015-04-29 | 2015-04-29 | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104811867B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861114B (zh) * | 2017-11-07 | 2020-12-29 | 哈尔滨工程大学 | 一种基于水声阵列空域反转的噪声抑制方法 |
CN108169732B (zh) * | 2018-02-28 | 2021-08-20 | 哈尔滨工程大学 | 一种基于扩展孔径声纳的变换域波束形成方法 |
CN110310651B (zh) * | 2018-03-25 | 2021-11-19 | 深圳市麦吉通科技有限公司 | 波束形成的自适应语音处理方法、移动终端及存储介质 |
CN108828536B (zh) * | 2018-04-25 | 2022-03-22 | 西安电子科技大学 | 基于二阶锥规划的宽带发射数字波束形成干扰设计方法 |
CN109143190B (zh) * | 2018-07-11 | 2021-09-17 | 北京理工大学 | 一种零陷展宽的宽带稳健自适应波束形成方法 |
CN109493844A (zh) * | 2018-10-17 | 2019-03-19 | 南京信息工程大学 | 基于fir滤波器的恒定束宽波束形成方法 |
CN116017230A (zh) * | 2021-10-21 | 2023-04-25 | 广州视源电子科技股份有限公司 | 麦克风阵列及其信号处理方法、装置、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147458A (zh) * | 2010-12-17 | 2011-08-10 | 中国科学院声学研究所 | 一种针对宽带声源的波达方向估计方法及其装置 |
CN103454615A (zh) * | 2013-03-01 | 2013-12-18 | 西安电子科技大学 | 稳健盲信号参数估计方法 |
-
2015
- 2015-04-29 CN CN201510212768.XA patent/CN104811867B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102147458A (zh) * | 2010-12-17 | 2011-08-10 | 中国科学院声学研究所 | 一种针对宽带声源的波达方向估计方法及其装置 |
CN103454615A (zh) * | 2013-03-01 | 2013-12-18 | 西安电子科技大学 | 稳健盲信号参数估计方法 |
Non-Patent Citations (3)
Title |
---|
Adaptive Beamforming Applied to a Cylindrical Sensor Array Using an Interpolated Array Transformation;Blomberg A E A等;《IEEE Journal of Oceanic Engineering》;20111223;第37卷(第1期);25-34 * |
Constrained Sensor System Resolves Closely Spaced Obscured Low-SNR Sources Using Virtual Beamforming;Marino C S等;《IEEE Transactions on Instrumentation and Measurement》;20130826;第63卷(第1期);231-242 * |
Grating Lobe Suppression for Distributied Digital Subarrays Using Virtual Filling;Feng B K等;《IEEE Antennas and Wireless Propagation Letters》;20131009;第12卷;1323-1326 * |
Also Published As
Publication number | Publication date |
---|---|
CN104811867A (zh) | 2015-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104811867B (zh) | 基于阵列虚拟扩展的麦克风阵列空域滤波方法 | |
US8363850B2 (en) | Audio signal processing method and apparatus for the same | |
CN109490850B (zh) | 主瓣干扰下宽带阵列自适应波束形成方法 | |
CN102324237B (zh) | 麦克风阵列语音波束形成方法、语音信号处理装置及系统 | |
CN109407055B (zh) | 基于多径利用的波束形成方法 | |
CN108375763B (zh) | 一种应用于多声源环境的分频定位方法 | |
CN110148420A (zh) | 一种适用于噪声环境下的语音识别方法 | |
CN101795150B (zh) | 强弱信号的波达方向与信源数估计方法 | |
CN102830387B (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
CN103837861B (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
CN103885045B (zh) | 基于子阵划分的循环联合自适应波束形成方法 | |
CN102944870A (zh) | 协方差矩阵对角加载的稳健自适应波束形成方法 | |
CN106301498A (zh) | 子带处理方法及频空级联的宽带自适应波束获取方法 | |
CN111123192A (zh) | 一种基于圆形阵列和虚拟扩展的二维doa定位方法 | |
CN103728601B (zh) | 雷达信号运动干扰空域-极化域联合稳健滤波方法 | |
CN108538306B (zh) | 提高语音设备doa估计的方法及装置 | |
CN105182302A (zh) | 一种抗快速运动干扰稳健零陷展宽波束形成方法 | |
CN105204008A (zh) | 一种基于协方差矩阵扩展的自适应天线波束形成零陷展宽方法 | |
CN110888105A (zh) | 一种基于卷积神经网络和接收信号强度的doa估计方法 | |
CN104777450A (zh) | 一种两级music麦克风阵列测向方法 | |
CN113593596B (zh) | 一种基于子阵划分的鲁棒自适应波束形成定向拾音方法 | |
CN110082789B (zh) | 基于圆阵的空时域自适应宽线性降秩波束形成方法 | |
CN113050035B (zh) | 一种二维定向拾音方法及装置 | |
CN109541526A (zh) | 一种利用矩阵变换的圆环阵方位估计方法 | |
CN108680907A (zh) | 一种基于观测矩阵的压缩感知mimo雷达干扰抑制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171121 |