发明名称 一种基于声矢量传感器的机动声源方位估计方法
摘要 本发明涉及一种基于声矢量传感器的机动声源方位估计方法,属于信号处理技术领域。本发明首先估计空域噪声协方差矩阵,通过噪声白化方法固化最大能量定位方法中的加权参数,以避免最优加权参数选择的一维搜索过程,提高了最大能量定向算法的估计精度,然后结合最大能量定向估计子输出和声源匀速运动的先验信息,在极坐标系下采用卡尔曼滤波技术进一步提高机动声源的方位估计精度。通过理论分析和仿真研究,本发明的基于声矢量传感器的机动声源方位估计与跟踪方法的估计精度优于原最大能量定向方法,并且由于采用了卡尔曼滤波,本发明的均方角度误差低于静态声源定位情况下的克拉美-罗下界。
申请公布号 CN104330768A 申请公布日期 2015.02.04
申请号 CN201410315852.X 申请日期 2014.07.03
申请人 河南科技大学 发明人 侯云山;金勇;翟红村;徐向艺;李佩佩;汤艳红;翟普杰
分类号 G01S3/802(2006.01)I 主分类号 G01S3/802(2006.01)I
代理机构 郑州睿信知识产权代理有限公司 41119 代理人 胡泳棋
主权项 一种基于声矢量传感器的机动声源方位估计方法,其特征在于,该方法的步骤如下:1)利用改进最大能量梯度法获取初始时刻声源方向矢量<img file="FDA0000532421500000011.GIF" wi="52" he="66" />和k时刻声源方向矢量<img file="FDA0000532421500000012.GIF" wi="73" he="66" />2)求出<img file="FDA0000532421500000013.GIF" wi="55" he="66" />对应极坐标下的俯仰角<img file="FDA0000532421500000014.GIF" wi="51" he="76" />与方位角<img file="FDA0000532421500000015.GIF" wi="92" he="66" />初始状态向量<img file="FDA0000532421500000016.GIF" wi="367" he="101" />同时计算k时刻声源方向矢量<img file="FDA0000532421500000017.GIF" wi="53" he="66" />对应极坐标下的俯仰角<img file="FDA0000532421500000018.GIF" wi="52" he="76" />与方位角<img file="FDA0000532421500000019.GIF" wi="79" he="66" />3)根据卡尔曼滤波方法,获取k时刻的状态预测<img file="FDA00005324215000000110.GIF" wi="201" he="73" />及其误差协方差矩阵P(k,k‑1),<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>&Phi;</mi><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>X</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005324215000000111.GIF" wi="592" he="79" /></maths>P(k,k‑1)=Φ(k,k‑1)P(k‑1,k‑1)Φ(k,k‑1)<sup>Τ</sup>+Γ(k‑1)Q(k‑1)Γ(k‑1)<sup>Τ</sup>其中,<img file="FDA00005324215000000112.GIF" wi="436" he="96" />为k时刻声源目标的状态向量,θ为声源信号入射向量的俯仰角,<img file="FDA00005324215000000113.GIF" wi="44" he="60" />为声源的方位角,Φ(k,k‑1)为k‑1时刻至k时刻的状态转移矩阵,Γ(k‑1)为系统噪声驱动阵,W(k‑1)为系统激励噪声序列;4)利用k时刻先验信息对以上预测信息进行修正,从而得到k时刻波达方向估计<img file="FDA00005324215000000114.GIF" wi="142" he="73" />及其误差协方差矩阵P(k,k):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mi>K</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>Z</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mi>H</mi><mover><mi>X</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005324215000000115.GIF" wi="903" he="79" /></maths>P(k,k)=P(k,k‑1)‑K(k)HP(k,k‑1)其中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>H</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA00005324215000000116.GIF" wi="378" he="153" /></maths>Z(k)为k时刻声源方向矢量<img file="FDA00005324215000000117.GIF" wi="51" he="66" />对应的极坐标下的量测方程,K(k)为卡尔曼滤波器的增益。
地址 471003 河南省洛阳市涧西区西苑路48号