发明名称 一种麦克风阵列及实现语音信号增强的方法
摘要 本发明公开了一种麦克风阵列及实现语音信号增强的方法,其特征在于,先将N个具有相同单元方向因子的阵元均匀径向地向外排列在一个半径为R的圆周上;形成2π弧度范围的信号接收处理的圆形阵列;或者将所述阵元均匀径向地向外排列在一个半径为R、圆心角为φ的一段圆弧上构成φ弧度范围信号接收处理的的部分圆形阵列;然后根据期望的信号方位,按照阵列的需要处理的音频信号的中心频率计算出各个阵元上的加权值以及某一2π/N弧度范围内的入射方向对应的加权向量,随之可以利用循环移位操作得到整个圆形阵列或部分圆形阵列其他入射方向的加权向量;最后对各个阵元上接收到的信号按计算出来的加权值进行加权合并。
申请公布号 CN101351058B 申请公布日期 2012.01.04
申请号 CN200810150881.X 申请日期 2008.09.09
申请人 西安交通大学 发明人 殷勤业;陈捷;邓科;王晨;张渭乐
分类号 H04R1/20(2006.01)I;H04R3/00(2006.01)I;G10L21/02(2006.01)I 主分类号 H04R1/20(2006.01)I
代理机构 西安通大专利代理有限责任公司 61200 代理人 朱海临
主权项 1.一种用麦克风阵列实现语音信号增强的方法,其特征在于,包括下述步骤:a)将多个具有相同单元方向因子的阵元均匀径向地向外排列在一个半径为R的圆周上以形成2π弧度范围的信号接收处理的圆形阵列或者将所述阵元均匀径向地向外排列在一个半径为R、圆心角为φ的一段圆弧上构成φ弧度范围信号接收处理的部分圆形阵列;所述圆形阵列或部分圆形阵列的相邻两个阵元与圆心连线的夹角为2π/N弧度,N是阵元的总数目,对圆形阵列而言,N取值范围为3至500;对部分圆形阵列而言,N取值范围为2至500;b)根据期望的信号入射角度,按照阵列的设计工作频率计算出各个阵元上的加权值,包括相位加权或/和幅度加权,所述阵列的设计工作频率为需要处理的音频信号的中心频率;对于加权值的计算,只需计算出某一2π/N弧度范围内的入射方向对应的加权向量,随之利用循环移位操作得到整个圆形阵列或部分圆形阵列其他入射方向的加权向量;所述某一2π/N弧度范围内的入射方向对应的加权向量为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mo>=</mo><mi>A</mi><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>,</mo><msup><mi>e</mi><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac><mi>r</mi><mo>[</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>d</mi></msub><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mo>-</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></msup><mo>,</mo><msup><mi>e</mi><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac><mi>r</mi><mo>[</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>d</mi></msub><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mn>2</mn><mo>-</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></msup><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msup><mi>e</mi><mrow><mi>j</mi><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><msub><mi>&lambda;</mi><mn>0</mn></msub></mfrac><mi>r</mi><mo>[</mo><msub><mrow><mi>cos</mi><mi>&theta;</mi></mrow><mi>d</mi></msub><mo>-</mo><mi>cos</mi><mrow><mo>(</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi></mrow><mi>N</mi></mfrac><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></msup><mo>)</mo></mrow><mi>T</mi></msup></mrow></math>]]></maths>式中:θ<sub>d</sub>表示为期望的信号入射角度;f<sub>0</sub>表示阵列的设计工作频率;<img file="FSB00000601053300012.GIF" wi="153" he="108" />表示该频率对应的波长;T表示矩阵转置;c)对各个阵元上接收到的信号按计算出来的加权值进行加权合并;步骤b)所述循环移位操作是指在期望的信号入射角度为θ<sub>d</sub>时计算出的加权向量为<img file="FSB00000601053300013.GIF" wi="131" he="58" />如果期望的信号入射角度改变<img file="FSB00000601053300014.GIF" wi="123" he="109" />k是任意一整数,即对应入射方向为<img file="FSB00000601053300015.GIF" wi="190" he="108" />时,加权向量为<img file="FSB00000601053300016.GIF" wi="407" he="67" />这里<img file="FSB00000601053300017.GIF" wi="380" he="69" />定义为<img file="FSB00000601053300018.GIF" wi="99" he="63" />的循环移位,即:k>0时<maths num="0002"><![CDATA[<math><mrow><mi>SHIFT</mi><mrow><mo>(</mo><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></msub></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow></math>]]></maths>k<0时<maths num="0003"><![CDATA[<math><mrow><mi>SHIFT</mi><mrow><mo>(</mo><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msup><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow></mrow></msub></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mrow><mo>(</mo><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mi>T</mi></msup></mrow></math>]]></maths>k=0时<maths num="0004"><![CDATA[<math><mrow><mi>SHIFT</mi><mrow><mo>(</mo><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>W</mi><mrow><msub><mi>f</mi><mn>0</mn></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow></msub></mrow></math>]]></maths>其中<img file="FSB00000601053300024.GIF" wi="156" he="58" />是加权向量<img file="FSB00000601053300025.GIF" wi="100" he="57" />的第i个元素;步骤c)所述加权合并的方法为将各阵元所接收到的信号按步骤b)得到的加权向量加权相加,加权合并后得到的信号为:<maths num="0005"><![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Real</mi><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>W</mi><mrow><msub><mi>f</mi><mi>c</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow><mi>H</mi></msubsup><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>或<maths num="0006"><![CDATA[<math><mrow><mi>y</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mi>Real</mi><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>A</mi><mi>k</mi></msub><msubsup><mi>W</mi><mrow><msub><mi>f</mi><mi>c</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>d</mi></msub></mrow><mi>H</mi></msubsup><mi>s</mi><mrow><mo>(</mo><mi>t</mi><mo>-</mo><msub><mi>t</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,(6)式为相位加权表达式;(7)式为幅度加权表达式;式中:t<sub>k</sub>表示信号到第k个阵元的延时;A<sub>k</sub>表示第k个阵元上的幅度加权值;Real[.]表示取实部;H表示矩阵的共轭转置。
地址 710049 陕西省西安市碑林区咸宁路28号