基于声源阵列的空间麦克风定位方法技术

本发明专利技术基于声源阵列的空间麦克风定位方法，涉及应用声波用多条由路径差测量确定的位置线的配合确定信号源的位置，利用四个声源的声音信号、参考麦克风和各个声源之间的固定距离实现了三维空间中的麦克风的在线定位，步骤是：声源的四元三维阵列分布；设置参考麦克风；标记待定位麦克风；用麦克风定位单元定位待定位麦克风在空间中的位置。本发明专利技术克服了现有技术中存在的信号采集操作较繁琐，效率较低，只能离线计算，无法实现麦克风位置的在线估计以及只针对平面的麦克风阵列进行了位置估计，无法实现空间分布麦克风的位置估计的缺陷。

【技术实现步骤摘要】
基于声源阵列的空间麦克风定位方法
本专利技术的技术方案涉及应用声波用多条由路径差测量确定的位置线的配合确定信号源的位置，具体地说是基于声源阵列的空间麦克风定位方法。
技术介绍
伴随着科学技术日新月异的发展，传感器技术的应用领域逐步扩大，其所面对的工作环境也愈来愈复杂。对麦克风位置进行准确的定位，能够帮助分布式麦克风阵列信号处理，而麦克风阵列在语音增强、声源定位和声源分离的应用中起着重要作用。同时，麦克风定位可以帮助解决当前机器人领域众多问题。例如，当机器人处于未知的环境中，对自身位姿未知时，如何进行自定位和环境识别，即同时自定位和地图构建一直都是智能机器人技术研究的热点和难点，为了解决这一难题，从机器人听觉系统的角度出发，确定出待定位麦克风的位置坐标，将机器人技术与机器人所处环境中的分布的麦克风位置进行有机结合，解决机器人同时自定位与构建环境地图的问题，从而帮助实现机器人全自主移动等。因此，确定麦克风阵列中的麦克风位置是一个现今该
中必须解决的技术难题。成都电子科技大学肖华在2008年发表的《麦克风阵列的校正方法研究》中，采用三个不同位置的校正声源来校正麦克风阵列中麦克风位置，该方法将一个声源先后放置于三个位置，分三次进行信号的采集与计算，最后将三次计算结果融合得到麦克风位置，但该方法存在信号采集操作较繁琐，效率较低，只能离线计算，无法实现麦克风位置的在线估计的缺陷。CN103439689B公开了一种分布式麦克风阵列中麦克风位置估计系统，该方法通过选择三个不同频带的声音信号作为三个声源的发声信号，在麦克风位置估计过程中一次性采集并分离出三个声源，实现麦克风位置的在线估计，将基于能量和时延的测距方法结合使用，先后进行距离的粗、精估计，采用三角形质心算法获取阵列中各麦克风的位置，对各种无解情况均有相应的处理，但此方法存在只针对平面的麦克风阵列进行了位置估计，无法实现空间分布麦克风的位置估计的缺陷。
技术实现思路
本专利技术所要解决的技术问题是：提供基于声源阵列的空间麦克风定位方法，利用四个声源的声音信号、参考麦克风和各个声源之间的固定距离实现了三维空间中的麦克风的在线定位，克服了现有技术中存在的信号采集操作较繁琐，效率较低，只能离线计算，无法实现麦克风位置的在线估计以及只针对平面的麦克风阵列进行了位置估计，无法实现空间分布麦克风的位置估计的缺陷。本专利技术解决该技术问题所采用的技术方案是：基于声源阵列的空间麦克风定位方法，具体步骤如下：第一步，声源的四元三维阵列分布：将四个声源按照正四面体的形式空间分布，定义声源1空间坐标为S1(x1,y1,z1)，声源2空间坐标为S2(x2,y2,z2)，声源3空间坐标为S3(x3,y3,z3)，声源4空间坐标为S4(x4,y4,z4)，构建成为一个三维空间坐标系，由此完成声源的四元三维阵列分布，并且规定在此空间分布的麦克风位置的表示都将以此坐标系为基准，上述四个声源依次发声，共产生四个声音信号分别为si(t)(i＝1,2,3,4)；第二步，设置参考麦克风：设置参考麦克风，该参考麦克风到上述第一步中四个声源中的第j个声源的距离为Lj(j＝1,2,3,4)，即参考麦克风到上述第一步中四个声源的距离分别为L1、L2、L3和L4，这四个距离的数值确定为已知且固定不变，同时参考麦克风的位置在距离上述各个声源小于10米的范围之内，即10米≥Lj＞0，保证参考麦克风接收到的各个声源的声音信号强度均为强信号；第三步，标记待定位麦克风：在上述第一步的正四面体的形式空间中，在以长＝60米，宽＝60米，高＝20米的长方体空间范围内，随机分布M＝200个麦克风，上述第一步的四个声源依次产生声音信号，随机分布的麦克风中能够接收到四个声音信号的均被标记为待定位麦克风，待定位麦克风的数量记为N，M≥N≥1；第四步，用麦克风定位单元定位待定位麦克风在空间中的位置：所述麦克风定位单元的操作程序如下：(4.1)四个声源依次产生声音信号，参考麦克风与待定位麦克风接收到四个声音信号：上述第一步的四个声源依次产生声音信号，上述第二步中设置的参考麦克风与上述第三步中的待定位麦克风均会接收到四个声音信号；设定所有麦克风接收到的声音信号模型是理想模型，即这里环境噪声用高斯白噪声近似代替，设定四个声源依次产生的声音信号为si(t)(i＝1,2,3,4)，其中，t表示各个声源发声对应的时间，i＝1,2,3,4，表示四个声源，某一个待定位麦克风接收到的四个声音信号中的某一个声音信号x(t)如下公式(1)所示：x(t)＝αs(t-τ)+n(t)(1)，公式(1)中，x(t)表示声音信号x与时间t的函数关系，t为该待定位麦克风接收到声音信号对应的时间，α表示声音信号到麦克风的幅度衰减系数，τ表示声音信号到麦克风的时间延迟、n(t)表示环境噪声信号，s(t-τ)表示具有时间延迟τ的声源产生的声音信号，s(t)和n(t)互不相关，(4.2)用经典的时延估计法中的广义互相关函数法估计声达时间差：上述第一步的四个声源依次产生的声音信号中的某一个声音信号到达上述第二步中设置的参考麦克风与上述第三步中的某一个待定位麦克风时间差定义为声达时间差，用经典的时延估计法中的广义互相关函数法估计声达时间差，分别为TDOA1、TDOA2、TDOA3、TDOA4；(4.3)计算某一个声源与参考麦克风的距离和该声源与待定位麦克风的距离的差值：定义ΔTij(i＝1,2.....N；j＝1,2,3,4)为上述第二步中设置的参考麦克风与上述第三步的待定位麦克风中的第i个待定位麦克风接收到上述第一步中四个声源中的第j个声源发出的第j个声音信号的时间差，该时间差ΔTij由测定得到，定义声音在空气中传播速度为C，将前述时间差ΔTij乘以声速C，得到上述第一步中四个声源中的第j个声源与上述第二步中设置的参考麦克风的距离和与上述第三步的待定位麦克风中的第i个待定位麦克风的距离的差值ΔLij为如下公式(2)所示，ΔLij＝ΔTij·C(2)，(4.4)计算待定位麦克风与上述某一个声源的距离：根据上述第二步中设置的参考麦克风到上述第一步中四个声源中的第j个声源的距离Lj(j＝1,2,3,4)，以及由公式(2)计算得到的第j个声源到参考麦克风与上述第三步的待定位麦克风中的第i个待定位麦克风的距离差ΔLij，由公式(3)计算得到上述第三步的待定位麦克风中的第i个待定位麦克风到上述第一步中四个声源中的第j个声源的距离dij如下：dij＝Lj+ΔLij＝Lj+ΔTij·C(3)，根据几何关系，通过求解，分别以上述第一步中的四个声源位置为球心，以上述第三步的待定位麦克风到相对应的上述第一步中的每个声源的距离为半径分别绘制的四个球的交点，来确定上述第三步的待定位麦克风中的第i个待定位麦克风所在空间环境中的位置如下方程式：方程式(4)中，xi,yi,zi表示上述第三步中的待定位麦克风中的第i个待定位麦克风所在空间环境中的位置，x1、y1、z1；x2、y2、z2；x3、y3、z3；x4、y4、z4依次表示上述第一步中的四个声源的位置，di1、di2、di3、di4分别表示上述第三步中的待定位麦克风中的第i个待定位麦克风分别到上述第一步中的四个声源中的每个声源的距离，由于上述第一步中的本文档来自技高网...

【技术保护点】
1.基于声源阵列的空间麦克风定位方法，其特征在于具体步骤如下：第一步，声源的四元三维阵列分布：将四个声源按照正四面体的形式空间分布，定义声源1空间坐标为S1(x1,y1,z1)，声源2空间坐标为S2(x2,y2,z2)，声源3空间坐标为S3(x3,y3,z3)，声源4空间坐标为S4(x4,y4,z4)，构建成为一个三维空间坐标系，由此完成声源的四元三维阵列分布，并且规定在此空间分布的麦克风位置的表示都将以此坐标系为基准，上述四个声源依次发声，共产生四个声音信号分别为si(t)(i＝1，2，3，4)；第二步，设置参考麦克风：设置参考麦克风，该参考麦克风到上述第一步中四个声源中的第j个声源的距离为Lj(j＝1,2,3,4)，即参考麦克风到上述第一步中四个声源的距离分别为L1、L2、L3和L4，这四个距离的数值确定为已知且固定不变，同时参考麦克风的位置在距离上述各个声源小于10米的范围之内，即10米≥Lj＞0，保证参考麦克风接收到的各个声源的声音信号强度均为强信号；第三步，标记待定位麦克风：在上述第一步的正四面体的形式空间中，在以长＝60米，宽＝60米，高＝20米的长方体空间范围内，随机分布M＝200个麦克风，上述第一步的四个声源依次产生声音信号，随机分布的麦克风中能够接收到四个声音信号的均被标记为待定位麦克风，待定位麦克风的数量记为N，M≥N≥1；第四步，用麦克风定位单元定位待定位麦克风在空间中的位置：所述麦克风定位单元的操作程序如下：(4.1)四个声源依次产生声音信号，参考麦克风与待定位麦克风接收到四个声音信号：上述第一步的四个声源依次产生声音信号，上述第二步中设置的参考麦克风与上述第三步中的待定位麦克风均会接收到四个声音信号；设定所有麦克风接收到的声音信号模型是理想模型，即这里环境噪声用高斯白噪声近似代替，设定四个声源依次产生的声音信号为si(t)(i＝1，2，3，4)，其中，t表示各个声源发声对应的时间，i＝1,2,3,4，表示四个声源，某一个待定位麦克风接收到的四个声音信号中的某一个声音信号x(t)如下公式(1)所示：x(t)＝αs(t‑τ)+n(t)    (1)，公式(1)中，x(t)表示声音信号x与时间t的函数关系，t为该待定位麦克风接收到声音信号对应的时间，α表示声音信号到麦克风的幅度衰减系数，τ表示声音信号到麦克风的时间延迟、n(t)表示环境噪声信号，s(t‑τ)表示具有时间延迟τ的声源产生的声音信号，s(t)和n(t)互不相关，(4.2)用经典的时延估计法中的广义互相关函数法估计声达时间差：上述第一步的四个声源依次产生的声音信号中的某一个声音信号到达上述第二步中设置的参考麦克风与上述第三步中的某一个待定位麦克风时间差定义为声达时间差，用经典的时延估计法中的广义互相关函数法估计声达时间差，分别为TDOA1、TDOA2、TDOA3、TDOA4；(4.3)计算某一个声源与参考麦克风的距离和该声源与待定位麦克风的距离的差值：定义ΔTij(i＝1,2.....N；j＝1,2,3,4)为上述第二步中设置的参考麦克风与上述第三步的待定位麦克风中的第i个待定位麦克风接收到上述第一步中四个声源中的第j个声源发出的第j个声音信号的时间差，该时间差ΔTij由测定得到，定义声音在空气中传播速度为C，将前述时间差ΔTij乘以声速C，得到上述第一步中四个声源中的第j个声源与上述第二步中设置的参考麦克风的距离和与上述第三步的待定位麦克风中的第i个待定位麦克风的距离的差值ΔLij为如下公式(2)所示，ΔLij＝ΔTij·C    (2)，(4.4)计算待定位麦克风与上述某一个声源的距离：根据上述第二步中设置的参考麦克风到上述第一步中四个声源中的第j个声源的距离Lj(j＝1,2,3,4)，以及由公式(2)计算得到的第j个声源到参考麦克风与上述第三步的待定位麦克风中的第i个待定位麦克风的距离差ΔLij，由公式(3)计算得到上述第三步的待定位麦克风中的第i个待定位麦克风到上述第一步中四个声源中的第j个声源的距离dij如下：dij＝Lj+ΔLij＝Lj+ΔTij·C    (3)，根据几何关系，通过求解，分别以上述第一步中的四个声源位置为球心，以上述第三步的待定位麦克风到相对应的上述第一步中的每个声源的距离为半径分别绘制的四个球的交点，来确定上述第三步的待定位麦克风中的第i个待定位麦克风所在空间环境中的位置如下方程式：...

【技术特征摘要】
1.基于声源阵列的空间麦克风定位方法，其特征在于具体步骤如下：第一步，声源的四元三维阵列分布：将四个声源按照正四面体的形式空间分布，定义声源1空间坐标为S1(x1,y1,z1)，声源2空间坐标为S2(x2,y2,z2)，声源3空间坐标为S3(x3,y3,z3)，声源4空间坐标为S4(x4,y4,z4)，构建成为一个三维空间坐标系，由此完成声源的四元三维阵列分布，并且规定在此空间分布的麦克风位置的表示都将以此坐标系为基准，上述四个声源依次发声，共产生四个声音信号分别为si(t)(i＝1，2，3，4)；第二步，设置参考麦克风：设置参考麦克风，该参考麦克风到上述第一步中四个声源中的第j个声源的距离为Lj(j＝1,2,3,4)，即参考麦克风到上述第一步中四个声源的距离分别为L1、L2、L3和L4，这四个距离的数值确定为已知且固定不变，同时参考麦克风的位置在距离上述各个声源小于10米的范围之内，即10米≥Lj＞0，保证参考麦克风接收到的各个声源的声音信号强度均为强信号；第三步，标记待定位麦克风：在上述第一步的正四面体的形式空间中，在以长＝60米，宽＝60米，高＝20米的长方体空间范围内，随机分布M＝200个麦克风，上述第一步的四个声源依次产生声音信号，随机分布的麦克风中能够接收到四个声音信号的均被标记为待定位麦克风，待定位麦克风的数量记为N，M≥N≥1；第四步，用麦克风定位单元定位待定位麦克风在空间中的位置：所述麦克风定位单元的操作程序如下：(4.1)四个声源依次产生声音信号，参考麦克风与待定位麦克风接收到四个声音信号：上述第一步的四个声源依次产生声音信号，上述第二步中设置的参考麦克风与上述第三步中的待定位麦克风均会接收到四个声音信号；设定所有麦克风接收到的声音信号模型是理想模型，即这里环境噪声用高斯白噪声近似代替，设定四个声源依次产生的声音信号为si(t)(i＝1，2，3，4)，其中，t表示各个声源发声对应的时间，i＝1,2,3,4，表示四个声源，某一个待定位麦克风接收到的四个声音信号中的某一个声音信号x(t)如下公式(1)所示：x(t)＝αs(t-τ)+n(t)(1)，公式(1)中，x(t)表示声音信号x与时间t的函数关系，t为该待定位麦克风接收到声音信号对应的时间，α表示声音信号到麦克风的幅度衰减系数，τ表示声音信号到麦克风的时间延迟、n(t)表示环境噪声信号，s(t-τ)表示具有时间延迟τ的声源产生的声音信号，s(t)和n(t)互不相关，(4.2)用经典的时延估计法中的广义互相关函数法估计声达时间差：上述第一步的四个声源依次产生的声音信...

【专利技术属性】
技术研发人员：孙昊，陈晓辉，王硕朋，翟葆朔，张恒，
申请(专利权)人：河北工业大学，
类型：发明
国别省市：天津,12