本发明提供一种基于声音信号利用双曲线交叉定位的室内定位方法，在室内不同位置安装n个发声设备；n≥3；n个发声设备的位置坐标均为已知值；n个发声设备按照预设时间间隔T1依次发出特定音频信号段；手持移动设备按采样频率f采集原始音频信号段，并发送到主机；主机对原始音频信号段和特定音频信号段进行相关性匹配，然后利用声音差和两个扬声器的位置建立单支双曲线，通过计算两个单支双曲线的交点得到手持移动设备的坐标，间接定位人的位置；具有定位精度高、无需其他专业设备等优点。
1.一种基于声音信号利用双曲线交叉定位的室内定位方法，其特征在于，包括以下步骤：S1，在室内不同位置安装n个发声设备；即：发声设备1、发声设备2…发声设备n；其中，n≥3；
在该室内任选坐标原点建立直角坐标系，则n个发声设备的位置坐标均为已知值，分别为：(x1，y1)、(x2，y2)...(xn，yn)；
S2，n个发声设备按照预设时间间隔T1依次发出特定音频信号段，其中，所述特定音频信号段的持续时间长度为T2，并且，T1＞T2；
手持移动设备为需要被进行室内定位的设备，所述手持移动设备按采样频率f采集时间长度为T3的原始音频信号段，并且，T3＞T2；
S3，所述手持移动设备将所述原始音频信号段发送到主机；所述主机读取预存储的所述特定音频信号段，然后，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，计算得到各个发声设备所发出的特定音频信号段从被所述发声设备发出到被所述手持移动设备接收所经过的时间间隔，分别为：t1、t2...tn；其中，t1与发声设备1对应、t2与发声设备2对应...tn与发声设备n对应；
S4，对于任意的i、j、k，设：ti与发声设备i对应、tj与发声设备j对应...tk与发声设备k对应；其中，i、j、k∈(1、2…n)；
则：计算时间差：δtij＝ti-tj；δtjk＝tj-tk；
计算距离差：δsij＝V×δtij；δsjk＝V×δtjk；其中，V为声音在空气中传播速度；
以发声设备i和发声设备j的位置坐标分别作为双曲线焦点，以δsij作为双曲线的半实轴长，根据δtij的正负值，建立双曲线分支ω1；
同理，以发声设备j和发声设备k的位置坐标分别作为双曲线焦点，以δsjk作为双曲线的半实轴长，根据δtjk的正负值，建立双曲线分支ω2；
S5，判断双曲线分支ω1和双曲线分支ω2是否合理，如果合理，则求解双曲线分支ω1和双曲线分支ω2在室内范围内的交点坐标，该交点坐标即为所述手持移动设备的位置坐标；
如果不合理，则改变i、j、k的值，重新计算获得两条双曲线分支，然后重复本步骤；
所述特定音频信号段通过以下方式产生：首先产生时间长度为5ms的2kHz正弦音频信号；再紧接着产生时间长度为20ms的频率随时间线性增加的音频信号，其中，音频信号随时间变化的函数为： 其中，t为事件，φ0为初始相位，f0为初始频率，k为频率变化率；此时，特定音频信号段的持续时间长度T2＝25ms；
其中，S5中，通过以下方法判断双曲线分支ω1和双曲线分支ω2是否合理：对于双曲线分支ω1，设发声设备i和发声设备j之间的直线连接距离为Lij，如果δsij＞Lij，则证明双曲线分支ω1不合理；
对于双曲线分支ω2，设发声设备j和发声设备k之间的直线连接距离为Ljk，如果δsjk＞Ljk，则证明双曲线分支ω2不合理。
2.根据权利要求1所述的基于声音信号利用双曲线交叉定位的室内定位方法，其特征在于，S2中，采样频率f＝44.1kHz；T3＝1s。
3.根据权利要求1所述的基于声音信号利用双曲线交叉定位的室内定位方法，其特征在于，S3中，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配具体为：所述主机对所述原始音频信号段进行滤波处理，然后将滤波处理后的音频信号段与所述特定音频信号段进行相关性匹配；其中，滤波处理具体为：对所述原始音频信号段进行快速傅利叶变换，去除位于2kHz以下和4kHz以上的音频分量，之后进行快速傅利叶逆变换，通过傅利叶逆变换得到的波形即为滤波处理后的音频信号段。
4.根据权利要求1所述的基于声音信号利用双曲线交叉定位的室内定位方法，其特征在于，S3中，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，计算得到各个发声设备所发出的特定音频信号段从被所述发声设备发出到被所述手持移动设备接收所经过的时间间隔，分别为：t1、t2...tn；具体为：S31，将所述原始音频信号段按手持移动设备接收时间先后顺序存入数组，得到数组B；
设数组B的长度为LB；数组B中各元素接收时间分别为：TB1、TB2…TBLB；
将特定音频信号段按发声时间先生顺序存入数组，得到数组A；其中，数组A的长度短于数组B的长度；设数组A的长度为LA；
S32，创建相关度数组C，设AvEvC1分别为数组A、B、C的第i个元素，则：数组C的长度Lc＝LB-LA+1；
即：数组C的Lc个元素分别为：C1、C1...CLc；各元素分别对应时刻：TB1、TB2…TBLc；
S33，将(TB1、C1)、(TB2、C2)…(TBLc、CLc)共Lc个点表示到直角坐标系上，其中，横坐标为TBi，纵坐标为Ci；
S34，将Lc个点连成曲线，将该曲线峰值最大的n个点所对应的TBi按时间从小到大顺序分别记为P1、P2…Pn；
则：S3中，t1＝P1-P0；t2＝P2-P0-T1；tn＝Pn-P0-T1(n-1)；其中，T1为发声设备发声的时间间隔；P0为手持移动设备开始采样的时间。
5.根据权利要求1所述的基于声音信号利用双曲线交叉定位的室内定位方法，其特征在于，S5之后，还包括：S6，创建一个窗口，显示室内地图；将得到的手持移动设备的位置坐标用有颜色的点显示到室内地图的对应位置上。
基于声音信号利用双曲线交叉定位的室内定位方法
技术领域
[0001] 本发明属于室内定位技术领域，具体涉及一种基于声音信号利用双曲线交叉定位的室内定位方法。
背景技术
[0002] GPS作为一种成熟的定位技术，已被各行各业广泛使用。但是，由于GPS信号存在室内信号强度较弱、易受遮挡等等特点，因此，难以应用于室内定位。
[0003] 目前，关于室内定位的研究成果主要包括以下几点：(1)谷歌应用GPS、WIFI信号、手机基站信号进行室内定位，发布了一款名为“Google Map Floor Plan Marker”的手机应用，但存在定位精度较低的不足。(2)诺基亚使用蓝牙模块进行室内定位，但是，需要在室内安装一个蓝牙装置作为发射台，通过手机和蓝牙发射台的通信完成定位，精度在30cm-100cm之间。由于需要在室内安装蓝牙装置，从而增加了室内定位成本。(3)博通公司使用WIFI、蓝牙、NFC技术以及手机内置的加速度传感器、陀螺仪完成室内定位，但是，需要手机内置博通的芯片，不方便广泛普及。
发明内容
[0004] 针对现有技术存在的缺陷，本发明提供一种基于声音信号利用双曲线交叉定位的室内定位方法，不需要使用特殊专业设备即能完成室内定位，并且，还具有定位精度高的优点。
[0005] 本发明采用的技术方案如下：
[0006] 本发明提供一种基于声音信号利用双曲线交叉定位的室内定位方法，包括以下步骤：
[0007] S1，在室内不同位置安装n个发声设备；即：发声设备1、发声设备2...发声设备n；其中，n≥3；
[0008] 在该室内任选坐标原点建立直角坐标系，则n个发声设备的位置坐标均为已知值，分别为：(x1，y1)、(x2，y2)...(xn，yn)；
[0009] S2，n个发声设备按照预设时间间隔T1依次发出特定音频信号段，其中，所述特定音频信号段的持续时间长度为T2，并且，T1>T2；
[0010] 手持移动设备为需要被进行室内定位的设备，所述手持移动设备按采样频率f采集时间长度为T3的原始音频信号段，并且，T3>T2；
[0011] S3，所述手持移动设备将所述原始音频信号段发送到主机；所述主机读取预存储的所述特定音频信号段，然后，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，计算得到各个发声设备所发出的特定音频信号段从被所述发声设备发出到被所述手持移动设备接收所经过的时间间隔，分别为：t1、t2...tn；其中，t1与发声设备1对应、t2与发声设备2对应...tn与发声设备n对应；
[0012] S4，对于任意的i、j、k，设：ti与发声设备i对应、tj与发声设备j对应...tk与发声设备k对应；其中，i、j、k∈(1、2...n)；
[0013] 则：计算时间差：δtij=ti-tj；δtjk=tj-tk；
[0014] 计算距离差：δsij=V×δtij；δsjk=V×δtjk；其中，V为声音在空气中传播速度；
[0015] 以发声设备i和发声设备j的位置坐标分别作为双曲线焦点，以δsij作为双曲线的半实轴长，根据δtij的正负值，建立双曲线分支ω1；
[0016] 同理，以发声设备j和发声设备k的位置坐标分别作为双曲线焦点，以δsjk作为双曲线的半实轴长，根据δtjk的正负值，建立双曲线分支ω2；
[0017] S5，判断双曲线分支ω1和双曲线分支ω2是否合理，如果合理，则求解双曲线分支ω1和双曲线分支ω2在室内范围内的交点坐标，该交点坐标即为所述手持移动设备的位置坐标；
[0018] 如果不合理，则改变i、j、k的值，重新计算获得两条双曲线分支，然后重复本步骤。
[0019] 优选的，S2中，时间间隔T1为100ms；
[0020] 所述特定音频信号段通过以下方式产生：首先产生时间长度为5ms的2kHz正弦音频信号；再紧接着产生时间长度为20ms的频率随时间线性增加的音频信号，其中，音频信号随时间变化的函数为： 其中，t为事件，φ0为初始相位，f0为初始频率，k为频率变化率；此时，特定音频信号段的持续时间长度T2=25ms。
[0021] 优选的，S2中，采样频率f=44.1kHz；T3=1s。
[0022] 优选的，S3中，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配具体为：
[0023] 所述主机对所述原始音频信号段进行滤波处理，然后将滤波处理后的音频信号段与所述特定音频信号段进行相关性匹配；其中，滤波处理具体为：
[0024] 对所述原始音频信号段进行快速傅利叶变换，去除位于2kHz以下和4kHz以上的音频分量，之后进行快速傅利叶逆变换，通过傅利叶逆变换得到的波形即为滤波处理后的音频信号段。
[0025] 优选的，S3中，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，计算得到各个发声设备所发出的特定音频信号段从被所述发声设备发出到被所述手持移动设备接收所经过的时间间隔，分别为：t1、t2...tn；具体为：
[0026] S31，将所述原始音频信号段按手持移动设备接收时间先后顺序存入数组，得到数组B；设数组B的长度为LB；数组B中各元素接收时间分别为：TB1、TB2...TBLB；
[0027] 将特定音频信号段按发声时间先生顺序存入数组，得到数组A；其中，数组A的长度短于数组B的长度；设数组A的长度为LA；
[0028] S32，创建相关度数组C，设Ai，Bi，Ci分别为数组A、B、C的第i个元素，则：数组C的长度Lc=LB-LA+1；
[0029] 即：数组C的Lc个元素分别为：C1、C1...CLc；各元素分别对应时刻：TB1、TB2...TBLc；
[0030] S33，将(TB1、C1)、(TB2、C2)...(TBLc、CLc)共Lc个点表示到直角坐标系上，其中，横坐标为TBi，纵坐标为Ci；
[0031] S34，将Lc个点连成曲线，将该曲线峰值最大的n个点所对应的TBi按时间从小到大顺序分别记为P1、P2...Pn；
[0032] 则：S3中，t1=P■P0；t2=P2-P0-T1；tn=Pn-P0-T1(n-1)；其中，T1为发声设备发声的时间间隔；P0为手持移动设备开始采样的时间。
[0033] 优选的，S5中，通过以下方法判断双曲线分支ω1和双曲线分支ω2是否合理：
[0034] 对于双曲线分支ω1，设发声设备i和发声设备j之间的直线连接距离为Lij，如果δsij>Lij，则证明双曲线分支ω1不合理；
[0035] 对于双曲线分支ω2，设发声设备j和发声设备k之间的直线连接距离为Ljk，如果δsjk>Ljk，则证明双曲线分支ω2不合理。
[0036] 优选的，S5之后，还包括：
[0037] S6，创建一个窗口，显示室内地图；将得到的手持移动设备的位置坐标用有颜色的点显示到室内地图的对应位置上。
[0038] 本发明的有益效果如下：
[0039] 本发明提供的基于声音信号利用双曲线交叉定位的室内定位方法，利用声音差和两个扬声器的位置建立单支双曲线，通过计算两个单支双曲线的交点得到手持移动设备的坐标，间接定位人的位置；具有定位精度高、无需其他专业设备等优点。
附图说明
[0040] 图1为本发明提供的基于声音信号利用双曲线交叉定位的室内定位方法的流程示意图；
[0041] 图2为本发明提供的基于声音信号利用双曲线交叉定位的室内定位方法的演示示意图。
具体实施方式
[0042] 以下结合附图对本发明进行详细说明：
[0043] 如图1所示，本发明提供一种基于声音信号利用双曲线交叉定位的室内定位方法，包括以下步骤：
[0044] S1，在室内不同位置安装n个发声设备；即：发声设备1、发声设备2...发声设备n；其中，n≥3；
[0045] 在该室内任选坐标原点建立直角坐标系，则n个发声设备的位置坐标均为已知值，分别为：(x1，y1)、(x2，y2)...(xn，yn)；
[0046] 本发明中，需要在室内至少布置3台发声设备，才能进行室内定位。其中，发声设备可以为能够发出特定频率声音的任何设备，例如，扬声器等，本发明对此并不进行限制。
[0047] S2，n个发声设备按照预设时间间隔T1依次发出特定音频信号段，其中，所述特定音频信号段的持续时间长度为T2，并且，T1>T2；
[0048] 手持移动设备为需要被进行室内定位的设备，所述手持移动设备按采样频率f采集时间长度为T3的原始音频信号段，并且，T3>T2；
[0049] 例如，实际应用中，特定音频信号段可以通过以下方式产生：首先产生时间长度为5ms的2kHz正弦音频信号；再紧接着产生时间长度为20ms的频率随时间线性增加的音频信号，其中，音频信号随时间变化的函数为： 其中，t为
事件，φ0为初始相位，f0为初始频率，k为频率变化率；此时，特定音频信号段的持续时间长度T2=25ms。采用特定音频信号段的声音，即：5ms的正弦信号和20ms的频率线性增加信号，解决了普遍手持移动设备中初始采集声音信号失真的问题。
[0050] 此时，可以采用以下参数值：T1=100ms；T2=25ms，T3=1s；采样频率f=44.1kHz。使用100ms作为发声间隔，大于T2的25ms，可有效防止两个扬声器的声音混叠，使分析定位更简单。
[0051] 以室内布置三台发声设备为例，即：发声设备1、发声设备2和发声设备3，则当T1=100ms、T2=25ms、T3=1s时，这三台发声设备按照以下方式发声：今初始时刻t=0，则，从t=0时刻开始，发声设备1持续发出长度为25ms的特定音频信号段，然后发声设备1停止发声；当t=100ms时，发声设备2持续发出长度为25ms的特定音频信号段，然后发声设备2停止发声；当t=200ms时，发声设备3持续发出长度为25ms的特定音频信号段，然后发声设备3停止发声；然后当t=300ms时，又轮到发声设备1发声。因此，发声设备1、发声设备2和发声设备3按上述顺序循环发声。从t=0时刻开始，手持移动设备按照44.1kHz的频率采集1s长度的原始音频信号段。需要说明的是，由于各发声设备所发出的声音经过反射或折射等才有可能到达手持移动设备，因此，手持移动设备所采集到的原始音频信号段即为环境声音。
[0052] 另外，本发明中，手持移动设备指目前市场上常见的可以进行音频采集和网络通信的手持移动设备，如智能手机、平板电脑等。
[0053] S3，所述手持移动设备将所述原始音频信号段发送到主机；所述主机读取预存储的所述特定音频信号段，然后，所述主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，计算得到各个发声设备所发出的特定音频信号段从被所述发声设备发出到被所述手持移动设备接收所经过的时间间隔，分别为：t1、t2...tn；其中，t1与发声设备1对应、t2与发声设备2对应...tn与发声设备n对应；
[0054] 本步骤中，主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配具体为：
[0055] 所述主机对所述原始音频信号段进行滤波处理，然后将滤波处理后的音频信号段与所述特定音频信号段进行相关性匹配；其中，滤波处理具体为：对所述原始音频信号段进行快速傅利叶变换，去除位于2kHz以下和4kHz以上的音频分量，之后进行快速傅利叶逆变换，通过傅利叶逆变换得到的波形即为滤波处理后的音频信号段。
[0056] 另外，主机对所述原始音频信号段和所述特定音频信号段进行相关性匹配，还需要计算得到各个发声设备所发出的特定音频信号段从被发声设备发出到被手持移动设备接收所经过的时间间隔，今其分别为：t1、t2...tn；则具体步骤为：
[0057] S31，将所述原始音频信号段按手持移动设备接收时间先后顺序存入数组，得到数组B；设数组B的长度为LB；数组B中各元素接收时间分别为：TB1、TB2...TBLB；
[0058] 将特定音频信号段按发声时间先生顺序存入数组，得到数组A；其中，数组A的长度短于数组B的长度；设数组A的长度为LA；
[0059] S32，创建相关度数组C，设Ai，Bi，Ci分别为数组A、B、C的第i个元素，则：数组C的长度Lc=LB-LA+1；
[0060] 即：数组C的Lc个元素分别为：C1、C1...CLc；各元素分别对应时刻：TB1、TB2...TBLc；
[0061] 例如，当采样频率为44.1kHz，即在一秒内采集44100个数据，存入数组，得到数组B；数组B的长度LB=44100，数组B中各元素接收时间分别为：TB1、TB2...TB44100；
[0062] 当采样频率为44.1kHz时，由于T2=25ms，则可以计算出，数组A的长度LA=1102。
[0063] 相关度数组C的长度Lc=44100-1102+1=42999；今Ai，Bi，Ci分别为数组A、B、C的第i个元素，则：
[0064] S33，将(TB1、C1)、(TB2、C2)...(TBLc、CLc)共Lc个点表示到直角坐标系上，其中，横坐标为TBi，纵坐标为Ci；
[0065] S34，将Lc个点连成曲线，将该曲线峰值最大的n个点所对应的TBi按时间从小到大顺序分别记为P1、P2...Pn；
[0066] 则：S3中，t1=P1-P0；t2=P2-P0-T1；tn=Pn-P0-T1(n-1)；其中，T1为发声设备发声的时间间隔；P0为手持移动设备开始采样的时间。
[0067] S4，对于任意的i、j、k，设：ti与发声设备i对应、tj与发声设备j对应...tk与发声设备k对应；其中，i、j、k∈(1、2...n)；
[0068] 则：计算时间差：δtij=ti-tj；δtjk=tj-tk；
[0069] 计算距离差：δsij=V×δtij；δsjk=V×δtjk；其中，V为声音在空气中传播速度；在标准情况下，V为340m／s。δsij指：发声设备i到手持移动设备之间的距离与发声设备j到手持移动设备之间的距离之差。
[0070] 根据双曲线定义，以发声设备i和发声设备j的位置坐标分别作为双曲线焦点，以δsij作为双曲线的半实轴长，根据δtij的正负值，建立双曲线分支ω1；需要说明的是，由于δtij非正即负，因此，最终只能得到双曲线的一个分支。
[0071] 同理，以发声设备j和发声设备k的位置坐标分别作为双曲线焦点，以δsjk作为双曲线的半实轴长，根据δtjk的正负值，建立双曲线分支ω2；
[0072] S5，判断双曲线分支ω1和双曲线分支ω2是否合理，如果合理，则求解双曲线分支ω1和双曲线分支ω2在室内范围内的交点坐标，该交点坐标即为所述手持移动设备的位置坐标；
[0073] 如果不合理，则改变i、j、k的值，重新计算获得两条双曲线分支，然后重复本步骤。
[0074] 可以采用以下方法求取两个双曲线分支的交点：将matlab的solve函数进行封装，将两个双曲线分支的相关参数传入，计算后即可得到输出结果。
[0075] 如图2所示，两条双曲线分支共有三个交点，分别为：交点A、交点B和交点C；由于交点A和交点B位于房间外，为不可行解，舍弃；而交点C位于房间内，为可行解，因此，最终计算得到的交点C坐标即为定位得到的手持移动设备的位置坐标。
[0076] 可以采用以下方法判断双曲线分支ω1和双曲线分支ω2是否合理：
[0077] 对于双曲线分支ω1，设发声设备i和发声设备j之间的直线连接距离为Lij，如果δsij>Lij，则证明双曲线分支ω1不合理；
[0078] 对于双曲线分支ω2，设发声设备j和发声设备k之间的直线连接距离为Ljk，如果δsjk>Ljk，则证明双曲线分支ω2不合理。
[0079] 为了使得至少存在两条合理的双曲线分支，可以布置多于三个扬声器，如果前两组扬声器中有一组不合理，则可以从后面的扬声器组中选择一组合理的进行替代。
[0080] S5之后，还包括：
[0081] S6，创建一个窗口，显示室内地图；将得到的手持移动设备的位置坐标用有颜色的点显示到室内地图的对应位置上。
[0082] 综上所述，本发明提供的基于声音信号利用双曲线交叉定位的室内定位方法，利用声音差和两个扬声器的位置建立单支双曲线，通过计算两个单支双曲线的交点得到手持移动设备的坐标，间接定位人的位置；具有定位精度高、无需其他专业设备等优点。
[0083] 以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视本发明的保护范围。
