[发明专利]一种基于声音位置指纹的室内声源定位方法

摘要

本发明一种基于声音位置指纹的室内声源定位方法，涉及应用声波确定信号源的位置的技术，操作步骤是离线采样，构建一个关于声达时间差和声音采样点位置坐标间关系的数据库，即声音位置指纹数据库；在线定位，确定待定位的声源的坐标位置。本发明方法利用定位算法比较待定位测点的声源信号与声音位置指纹数据库中每个参考点的邻近性，从而实现对声源待测点位置坐标的确定，克服了现有的声音定向方法只能定出声源的方向，而不能定出声源的具体位置，不能实现对室内声源进行定位的缺陷。

主权项

一种基于声音位置指纹的室内声源定位方法，其特征在于：构建一个关于声达时间差和声音采样点位置坐标间关系的声音位置指纹数据库，引入了位置指纹的定位算法，声音位置指纹数据库中存储的是参考点的位置坐标和对应的声达时间差，具体步骤如下：第一步，离线采样，构建一个关于声达时间差和声音采样点位置坐标间关系的数据库，即声音位置指纹数据库：(1.1)确定采样布局图：测定所指定的室内的面积和确定定位精度，根据所指定的室内的面积和定位精度，将所指定的室内划分成0.5m×0.5m的网格，把每个网格的节点作为参考点，得到n个参考节点，由此确定采样布局图；(1.2)设定声音采集点：在上述室内设定四个声音采集点，在每个采集点处分别布置麦克风0、麦克风1、麦克风2和麦克风3，以在采集点处采集声源信号；(1.3)采集声源信号：在上述(1.1)步中的每一个参考点处分别使用同一声源发出声音，用上述(1.2)步设置的麦克风0、麦克风1、麦克风2和麦克风3采集声音信号，声音信号经过声音放大模块的放大，用数据采集模块对放大的数据进行采集并实时的传输到上位机；(1.4)求麦克风1、麦克风2和麦克风3与基准麦克风0的声达时间差：用广义互相关函数法计算声源的声音到达不同麦克风的声达时间差，先求两个电压信号间的互功率谱函数并在频域内进行加权运算，再将频域的函数利用反傅里叶变换的方法转换到时域，得到两组信号的互相关函数，这个互相关函数的峰值所对应的时间为麦克风1、麦克风2和麦克风3与基准麦克风0的的声达时间差，具体计算方法如下：在第i个参考节点处，设s(t)为目标声源，x0(t)、x1(t)、x2(t)和x3(t)分别为麦克风0、麦克风1、麦克风2和麦克风3接收到的信号，n0(t)、n1(t)、n2(t)和n3(t)分别为麦克风0、麦克风1、麦克风2和麦克风3接收到的附加噪声信号，麦克风0、麦克风1、麦克风2和麦克风3得到的信号分别表示为：x0(t)＝α0s(t)+n0(t)x1(t)＝α1s(t‑τ)+n1(t)x2(t)＝α2s(t‑τ)+n2(t)x3(t)＝α3s(t‑τ)+n3(t)其中α0、α1、α2和α3分别是声音信号到达麦克风后的衰减系数，τ为两个麦克风拾取信号的时间延迟，即声达时间差TDOA,此处声达时间差TDOA采用传统的互相关方法进行估计，其互相关函数为：Rx0x1(τ)=∫0πψ01(ω)Gx0x1(ω)ejωτdω]]>Rx0x2(τ)=∫0πψ02(ω)Gx0x2(ω)ejωτdω]]>Rx0x3(τ)=∫0πψ03(ω)Gx0x3(ω)ejωτdω]]>其中，是麦克风0和麦克风1获取信号的互相关功率谱，而X0(ω)是x0(t)由时域经过傅里叶变换转到频域的信号，是频域内信号1的共轭，信号ψ01(ω)是权函数，在此权函数取1，是广义互相关功率谱，最后再经过傅立叶反变换到时域得到时延：的峰值处即为麦克风1和基准麦克风0的声达时间差TDOAi1；是麦克风0和麦克风2获取信号的互相关功率谱，而X0(ω)是x0(t)由时域经过傅里叶变换转到频域的信号，是频域内信号2的共轭，信号ψ02(ω)是权函数，在此权函数取1，是广义互相关功率谱，最后再经过傅立叶反变换到时域得到时延：的峰值处即为麦克风2和基准麦克风0的声达时间差TDOAi2；是麦克风0和麦克风3获取信号的互相关功率谱，而X0(ω)是x0(t)由时域经过傅里叶变换转到频域的信号，是频域内信号3的共轭，信号ψ03(ω)是权函数，在此权函数取1，是广义互相关功率谱，最后再经过傅立叶反变换到时域得到时延：的峰值处即为麦克风3和基准麦克风0的声达时间差TDOAi3；由此，对每一个参考点进行数据采集，把参考点的位置坐标和对应的声达时间差TDOA存入到声音位置指纹库中，构成声音位置指纹库；该声音位置指纹库，用M表示如下：M=L1TDOA11TDOA12TDOA13............LiTDOAi1TDOAi2TDOAi3............LnTDOAn1TDOAn2TDOAn3]]>其中Li表示第i个参考点的位置信息，位置信息表示如下：Li＝(xi,yi)，TDOAi1、TDOAi2、TDOAi3分别表示在第i个参考点处麦克风1、麦克风2、麦克风3到基准麦克风0的声达时间差，n表示参考点的个数，n个参考点的指纹组成了声音位置指纹库；第二步，在线定位，确定待定位的声源的坐标位置：(2.1)测量并采集待定位点声源信号到麦克风1、麦克风2和麦克风3与基准麦克风0的声达时间差：在待定位点用同样的声源发出声音，用广义互相关函数法计算声源的声音到达不同麦克风的声达时间差，先求两个电压信号间的互功率谱函数并在频域内进行加权运算，再将频域的函数利用反傅里叶变换的方法转换到时域，得到两组信号的互相关函数，这个互相关函数的峰值所对应的时间为麦克风1、麦克风2和麦克风3与基准麦克风0的的声达时间差(TDOAj1、TDOAj2、TDOAj3)，具体计算过程同上述(1.4)步；(2.2)利用定位算法将采集的待定位声源声音的声达时间差与第一步构建的声音位置指纹数据库中的信息进行匹配定位，计算出该声源的位置坐标：用K最近邻法进行信息匹配定位，具体过程如下：①计算待测点与所有参考点的欧氏距离：di=Σm=13(TDOAjm-TDOAim)2]]>其中，di表示待测点到第i个参考点的欧氏距离，m表示测的的声达时间差的个数，TDOAjm表示待测点测得的第m个麦克风与基准麦克风的声达时间差，TDOAim表示第i个参考点的第m个声达时间差；②对求出的上述(1.1)步中的n个参考节点的di进行比较，取出3个最小的欧氏距离da,db,dc，提取其坐标分别为A(xa,ya)、B(xb,yb)、C(xc,yc)；③求取A(xa,ya)、B(xb,yb)、C(xc,yc)这三个点的质心，质心的坐标就是待定位的声源的坐标，求取质心的方法是：(xj,yj)=db+dcda+db+dc(xa,ya)+da+dcda+db+dc(xb,yb)+da+dbda+db+dc(xc,yc).]]>其中(xj,yj)是待测点的坐标，da,db,dc是上述步骤②中测得的欧氏距离；由此，完成待定位的室内声源定位。