发明名称 近场径向干扰源抑制算法及其球麦克风阵列语音采集装置
摘要 一种近场径向干扰源抑制算法及其球麦克风阵列语音采集装置,当期望源位于近场且期望源和干扰源的来波方向和频率相同时利用空间径向抵消滤波算法有效抑制干扰,该球麦克风阵列语音采集装置采用拓扑结构为正十二面体的球麦克风阵列,阵列冗余度小,设计、计算简单,同时可精确确定麦克风在球面上的位置,不会产生位置误差。
申请公布号 CN102855409B 申请公布日期 2014.02.26
申请号 CN201210350168.6 申请日期 2012.09.20
申请人 辽宁工业大学 发明人 王冬霞;程宗明;郭庆龙
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 锦州辽西专利事务所 21225 代理人 李辉
主权项 1.一种近场径向干扰源抑制算法,其特征在于,包括如下具体步骤:第一步,利用球麦克风阵列语音采集装置采集期望源的声压信号,所述球麦克风阵列语音采集装置的麦克风阵列的拓扑结构采用正十二面体,所述正十二面体的每个顶点处均设有一路麦克风,利用基于球阵列的MUSIC算法和基于低频不同模式差近场期望源距离估计算法估计出期望源的位置<img file="2012103501686100001DEST_PATH_IMAGE001.GIF" wi="65" he="25" />,<i>r</i><sub><i>s</i></sub>表示期望源距球心距离,<i>θ</i><sub><i>s</i></sub>表示入射俯仰角,<i>φ</i><sub><i>s</i></sub>表示入射方位角;第二步,将球麦克风阵列接收到的声压信号进行球傅里叶变换,声压的球傅里叶变换公式如下:<img file="2012103501686100001DEST_PATH_IMAGE002.GIF" wi="220" he="48" />上式中<img file="DEST_PATH_IMAGE003.GIF" wi="46" he="25" />表示声压的球傅里叶变换系数,<img file="2012103501686100001DEST_PATH_IMAGE004.GIF" wi="21" he="26" />为与采样方案相关的权值参数,<i> J</i>表示球阵列所含的阵元数,<img file="DEST_PATH_IMAGE005.GIF" wi="74" he="26" />表示第<i>j</i>个麦克风所采集的声压,<img file="2012103501686100001DEST_PATH_IMAGE006.GIF" wi="86" he="26" />为第<i>j</i>个麦克风在球阵表面所处的位置,<img file="DEST_PATH_IMAGE007.GIF" wi="18" he="26" />和<img file="2012103501686100001DEST_PATH_IMAGE008.GIF" wi="18" he="26" />分别表示球坐标的俯仰角和方位角,<img file="DEST_PATH_IMAGE009.GIF" wi="49" he="26" />表示n阶m度的球谐函数;第三步,根据第一步所估计出的期望源位置<img file="204495DEST_PATH_IMAGE001.GIF" wi="65" he="25" />,对径向陷波滤波器和多尔夫-切比雪夫径向滤波器进行设计,得到两个波束形成系数<i>d</i><sub><i>n</i>1</sub>(<i>k</i>)和<i>d</i><sub><i>n</i>2</sub>(<i>k</i>),相应滤波器的阵列权值分别为<img file="DEST_PATH_IMAGE011.GIF" wi="52" he="28" />和<img file="2012103501686100001DEST_PATH_IMAGE012.GIF" wi="53" he="28" />,阵列权值按照以下公式计算为:<img file="DEST_PATH_IMAGE013.GIF" wi="302" he="49" />其中<img file="2012103501686100001DEST_PATH_IMAGE014.GIF" wi="46" he="26" />(g = 1, 2),控制空间滤波器的径向行为和波束图,<img file="DEST_PATH_IMAGE015.GIF" wi="80" he="25" />表示阵列的指向方向,<img file="2012103501686100001DEST_PATH_IMAGE016.GIF" wi="45" he="25" />由阵列边界决定,即<img file="DEST_PATH_IMAGE017.GIF" wi="204" he="45" /><img file="2012103501686100001DEST_PATH_IMAGE018.GIF" wi="41" he="22" />为球汉克尔函数,可表示球面波前的径向衰减,<img file="DEST_PATH_IMAGE019.GIF" wi="41" he="22" />为球贝塞尔函数;第四步,空间径向抵消滤波器的主支路利用多尔夫-切比雪夫径向滤波器对信号初步处理,抑制一部分噪声和干扰,辅支路用径向陷波滤波器对信号进行处理,完全抑制掉期望信号保留噪声和干扰,两支路输出结果分别为<img file="2012103501686100001DEST_PATH_IMAGE020.GIF" wi="40" he="25" />和<img file="DEST_PATH_IMAGE021.GIF" wi="41" he="25" />,分别为:<img file="2012103501686100001DEST_PATH_IMAGE022.GIF" wi="184" he="46" /><img file="DEST_PATH_IMAGE023.GIF" wi="186" he="46" />第五步,将两支路的输出结果进行加权求和,得到空间径向抵消滤波器的输出为<img file="2012103501686100001DEST_PATH_IMAGE024.GIF" wi="34" he="22" />:<img file="DEST_PATH_IMAGE025.GIF" wi="156" he="25" />式中的<img file="2012103501686100001DEST_PATH_IMAGE026.GIF" wi="15" he="14" />为利用自适应算法得出的最优权值。
地址 121000 辽宁省锦州市古塔区士英街169号