CN103267961B - 一种移动终端的测向方法、系统及该移动终端 - Google Patents

一种移动终端的测向方法、系统及该移动终端 Download PDF

Info

Publication number
CN103267961B
CN103267961B CN201310143217.3A CN201310143217A CN103267961B CN 103267961 B CN103267961 B CN 103267961B CN 201310143217 A CN201310143217 A CN 201310143217A CN 103267961 B CN103267961 B CN 103267961B
Authority
CN
China
Prior art keywords
mobile terminal
lambda
sound wave
vector
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201310143217.3A
Other languages
English (en)
Other versions
CN103267961A (zh
Inventor
黄文超
熊焰
王行甫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN201310143217.3A priority Critical patent/CN103267961B/zh
Publication of CN103267961A publication Critical patent/CN103267961A/zh
Application granted granted Critical
Publication of CN103267961B publication Critical patent/CN103267961B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Telephone Function (AREA)

Abstract

本发明公开了一种移动终端的测向方法及、系统及该移动终端,其中,该方法包括:第一移动终端以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;所述第一移动终端根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;并根据所述相对速度与该第一移动终端在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得所述第一移动终端相对于所述第二移动终端的方位角通过采用本发明公开的方法及系统,实现了相邻手机用户之间的相互定位时不需要额外的装置来进行计算,且测向精度高。

Description

一种移动终端的测向方法、系统及该移动终端
技术领域
本发明涉及定位技术领域,尤其涉及一种移动终端的测向方法及、系统及该移动终端。
背景技术
现如今,智能手机用户之间一般采用全球定位系统来交互各自的地理位置信息,然而在室内场景中,全球定位系统失效,以致用户之间无法进行相互定位。目前,手机用户可以使用声波的传输时间差来得到他们之间的精确距离,但仍无任何方法来获取他们之间的精确方向。所以,智能手机之间的精确测向方法是手机用户之间相互定位的瓶颈问题。
发明内容
本发明的目的是提供一种移动终端的测向方法及、系统及该移动终端,从而实现相邻手机用户之间的相互定位。
本发明的目的是通过以下技术方案实现的:
一种移动终端的测向方法,该方法包括:
第一移动终端以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;
所述第一移动终端根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;
并根据所述相对速度与该第一移动终端在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得所述第一移动终端相对于所述第二移动终端的方位角。
一种测向的移动终端,包括:
多普勒频率偏移计算单元,用于以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;
相对速度计算单元,用于根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;
方位角获取单元,用于根据所述相对速度,与其自身在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得相对于所述第二移动终端的方位角。
一种移动终端的测向系统,该系统包括:前述的第一移动终端,以及第二移动终端;
其中,所述第二移动终端,用于通过扬声器发送的指定频率的声波信号。
由上述本发明提供的技术方案可以看出,通过利用多普勒效应并结合移动终端上的传感器来测定对方移动终端的方向,不需要额外的装置来进行计算,测向精度高,本发明的方案具有零成本、使用方便、和精确的优点。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他附图。
图1为本发明实施例一提供的一种移动终端的测向方法的流程图;
图2为本发明实施例二提供的一种测向的移动终端的示意图;
图3为本发明实施例三提供的一种移动终端的测向系统的示意图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明的保护范围。
实施例一
图1为本发明实施例一提供的一种移动终端的测向方法的流程图。如图1所示,该方法主要包括如下步骤:
步骤11、第一移动终端以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量。
本发明实施例中的移动终端可以是智能手机、平板电脑、掌上电脑等。
在进行定位之前可以约定扬声器声波信号的频率,以便准确的识别到所需侧向的移动终端。
本发明实施例基于声波信号的多普勒效应来进行侧向,因此,用于测向的移动终端需要以晃动的方法接收声波信号。
当手机接受到声波信号后,首先利用带通滤波器过滤多余的噪声,提取出每一扬声器发送的声波信号;然后,通过自动增益控制,使每一扬声器声波信号的振幅归一化;最后,利用锁相环计算振幅归一化后的每一扬声器声波信号的多普勒频率偏移量与相位偏移量。
具体的:多普勒频率偏移量可以使用如下方法进行计算:
首先,计算声波信号的相位偏移量,其公式为:
θ[k+1]=θ[k]-μLPF{rc[k]sin(2πfakTs+θ[k])};
其中,μ为预先设定的常量,LPF为低通滤波器,fa为声波的频率,Ts为声波的采样周期,rc[k]为经过带通滤波器和自动增益控制之后的声波,θ[k]为第k次声采样时相位的计算值。
再根据相位偏移量计算频率偏移量:
f ( t ) = 1 2 π dθ ( t ) dt ;
其中,θ(t)=θ(kTs)=θ[k]。
步骤12、所述第一移动终端根据所述多普勒频率偏移量,计算与第二移动终端的相对速度。
扬声器声波信号的多普勒频率偏移量与移动终端的相对速度v成正比,因此,可以用来计算相对速度,其公式为:
V = V a f a f ;
其中,Va为声波信号的传播速度;fa为对应扬声器发送声波信号的频率;f为多普勒偏移量。
步骤13、并根据所述相对速度与该第一移动终端在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得所述第一移动终端相对于所述第二移动终端的方位角。
通常,我们直接利用移动终端中自带的加速度传感器与陀螺仪计算其在世界坐标系下的速度矢量;具体的:可以通过加速度传感器及陀螺仪进行采样,求得世界坐标系下的加速度ax、ay与az从而可以估算出手机在世界系下的速度矢量vx、vy与vz
其公式为: v t [ k ] = v t [ 0 ] + Σ i = 0 k - 1 T [ i ] a t [ i ] + Σ i = 0 k - 1 T [ i ] σ t [ i ]
其中,vt[k]为第i+1个采样对应的速度矢量中的vx、vy或vz;vt[0]为对应的初始速度;T[i]表示第i个采样到第i+1个采样的时间间隔;at[i]表示第i+1个采样时所求的世界坐标系中对应的ax、ay或az;σt[i]表示at[i]的测量误差,并假定σt在摇晃手机的这段短时间内为恒定值ex
例如,计算X轴的速度矢量公式为: v x [ k ] = v x [ 0 ] + Σ i = 0 k - 1 T [ i ] a x [ i ] + Σ i = 0 k - 1 T [ i ] σ x [ i ] ,
其中,vx[0]为X轴的初始速度;T[i]表示第i个采样到第i+1个采样的时间间隔,ax[i]表示第i+1个采样时所求的世界坐标系X轴坐标的加速度,σx[i]表示ax[i]的测量误差,并假定σx在摇晃手机的这段短时间内为恒定值ex,vx[k]则为第i+1个采样所对应的速度值,同理vy与vz可使用类似的方法计算。
再将步骤12中计算获得的相对速度与及第一移动终端在世界坐标系下的速度矢量计算第二移动终端中扬声器声波信号的方向矢量λx、λy与λz,其公式为:
λ x v x [ k ] + λ y v y [ k ] + λ z v z [ k ] = v , ∀ k ;
其中,为任意值符号,vx[k]、vy[k]与vz[k]为第k次采样计算到的X轴、Y轴与Z轴的速度矢量;。#
进一步,为简化计算,可以仅计算二维方向的角度(不求λz),并去除vz误差带来的精度影响,则其公式为:
λ x v x [ k ] + λ y v y [ k ] = v = v a f a · f [ k ] , ∀ k ;
则解下列方程组,求出:λx与λy
w x [ 0 ] w y [ 0 ] 1 t [ 0 ] w x [ 1 ] w y [ 1 ] 1 t [ 1 ] · · · · · · · · · · · · w x [ n ] w y [ n ] 1 t [ n ] λ x λ y λ 0 λ 1 = v a f a · f [ 0 ] f [ 1 ] · · · f [ n ] ;
其中, w x [ k ] = Σ i = 0 k - 1 T [ i ] a x [ i ] ; w y [ k ] = Σ i = 0 k - 1 T [ i ] a y [ i ] ; λ0=λxvx[0]+λyvy[0];λ1=λxexyey,ey为测量ay所产生的误差。
根据λx与λy,求声源相对于手机的方位角α:
&alpha; = arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x &GreaterEqual; 0 &pi; + arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x < 0 .
本发明实施例通过利用多普勒效应并结合移动终端上的传感器来测定对方移动终端的方向,不需要额外的装置来进行计算,测向精度高,本发明的方案具有零成本、使用方便、和精确的优点。
实施例二
图2为本发明实施例二提供的一种测向的移动终端的示意图。如图2所示,该移动终端主要包括:
多普勒频率偏移计算单元21,用于以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;
相对速度计算单元22,用于根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;
方位角获取单元23,用于根据所述相对速度,与其自身在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得相对于所述第二移动终端的方位角。
所述方位角获取单元23包括:速度矢量计算单元231,用于计算世界坐标系下的速度矢量;且该单元包括:
加速度矢量获取单元2311,用于利用移动终端中的加速度传感器与陀螺仪进行采样,计算得到世界坐标系下的加速度矢量ax、ay与az
速度矢量估算单元2312,用于根据世界坐标系下的加速度矢量ax、ay与az估算世界坐标系下的速度矢量vx、vy与vz
其计算公式为: v t [ k ] = v t [ 0 ] + &Sigma; i = 0 k - T [ i ] a t [ i ] + &Sigma; i = 0 k - 1 T [ i ] &sigma; t [ i ] ;
其中,vt为速度矢量中的vx、vy或vz;T[i]表示第i个采样到第i+1个采样的时间间隔;at[i]表示第i+1个采样时所求的世界坐标系中对应的ax、ay或az;σt[i]表示at[i]的测量误差。
所述方位角获取单元23包括:声波信号的方向矢量计算单元232,用于计算所述第二终端中扬声器声波信号的方向矢量,其公式包括:
&lambda; x v x [ k ] + &lambda; y v y [ k ] = v , &ForAll; k ;
其中,λx与λy第二终端中扬声器的声波信号在X轴与Y轴的方向矢量;vx[k]与vy[k]为第k次采样计算到的X轴与Y轴的速度矢量;为任意值符号;
方位角计算单元233,用于根据计算到的第二终端中扬声器声波信号的方向矢量,计算所述第二移动终端相对于所述第一移动终端的方向,其公式包括:
&alpha; = arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x &GreaterEqual; 0 &pi; + arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x < 0 .
实施例三
图3为本发明实施例三提供的一种移动终端的测向系统的示意图。如图3所示,该系统主要包括:
实施例二所述的移动终端为第一移动终端31,以及第二移动终端32;其中,所述第二移动终端32,用于通过扬声器发送的指定频率的声波信号。
需要说明的是,上述设备、系统中包含的各个功能单元所实现的功能的具体实现方式在前面的各个实施例中已经有详细描述,故在这里不再赘述。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,仅以上述各功能单元的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能单元完成,即将设备、系统的内部结构划分成不同的功能单元,以完成以上描述的全部或者部分功能。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明披露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书的保护范围为准。

Claims (8)

1.一种移动终端的测向方法,其特征在于,该方法包括:
第一移动终端以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;
所述第一移动终端根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;
并根据所述相对速度与该第一移动终端在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得所述第一移动终端相对于所述第二移动终端的方位角;
其中,计算所述第一移动终端在世界坐标系下的速度矢量的步骤包括:
利用移动终端中的加速度传感器与陀螺仪进行采样,计算得到世界坐标系下的加速度矢量ax、ay与az;根据世界坐标系下的加速度矢量ax、ay与az估算世界坐标系下的速度矢量vx、vy与vz
其计算公式为: v t [ k ] = v t [ 0 ] + &Sigma; i = 0 k - 1 T [ i ] a t [ i ] + &Sigma; i = 0 k - 1 T [ i ] &sigma; t [ i ] ;
其中,vt[k]为第i+1个采样对应的速度矢量中的vx、vy或vz;vt[0]为对应的初始速度;T[i]表示第i个采样到第i+1个采样的时间间隔;at[i]表示第i+1个采样时所求的世界坐标系中对应的ax、ay或az;σt[i]表示at[i]的测量误差。
2.根据权利要求1所述的方法,其特征在于,计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量的步骤包括:
利用带通滤波器过滤噪声,提取出扬声器发送的正弦声波信号;
通过自动增益控制,使扬声器声波信号的振幅归一化;
利用锁相环计算振幅归一化后的扬声器声波信号的多普勒频率偏移量。
3.根据权利要求1所述的方法,其特征在于,计算第一移动终端与第二移动终端的相对速度的公式包括:
V = V a f a f ;
其中,Va为声波信号的传播速度;fa为对应扬声器发送声波信号的频率;f为多普勒偏移量。
4.根据权利要求1-3任一项所述的方法,其特征在于,计算所述第二终端中扬声器声波信号的方向矢量的公式包括:
&lambda; x v x [ k ] + &lambda; y v y [ k ] = v , &ForAll; k ;
其中,λx与λy第二终端中扬声器的声波信号在X轴与Y轴的方向矢量;vx[k]与vy[k]为第k次采样计算到的X轴与Y轴的速度矢量;为任意值符号。
5.根据权利要求4所述的方法,其特征在于,计算所述第一移动终端相对于所述第二移动终端的方位角的公式包括:
&alpha; = arcsin &lambda; &lambda; x 2 + &lambda; y 2 &lambda; x &GreaterEqual; 0 &pi; + arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x < 0 .
6.一种测向的移动终端,其特征在于,包括:
多普勒频率偏移计算单元,用于以晃动的方式接收第二移动终端中扬声器发送的指定频率的声波信号,并计算所述第二移动终端中扬声器声波信号的多普勒频率偏移量;
相对速度计算单元,用于根据所述多普勒频率偏移量,计算与第二移动终端的相对速度;
方位角获取单元,用于根据所述相对速度,与其自身在世界坐标系下的速度矢量计算所述第二终端中扬声器声波信号的方向矢量,获得相对于所述第二移动终端的方位角;
其中,所述方位角获取单元包括:速度矢量计算单元,用于计算世界坐标系下的速度矢量;且该单元包括:
加速度矢量获取单元,用于利用移动终端中的加速度传感器与陀螺仪进行采样,计算得到世界坐标系下的加速度矢量ax、ay与az
速度矢量估算单元,用于根据世界坐标系下的加速度矢量ax、ay与az估算世界坐标系下的速度矢量vx、vy与vz
其计算公式为: v t [ k ] = v t [ 0 ] + &Sigma; i = 0 k - 1 T [ i ] a t [ i ] + &Sigma; i = 0 k - 1 T [ i ] &sigma; t [ i ] ;
其中,vt[k]为第i+1个采样对应的速度矢量中的vx、vy或vz;T[i]表示第i个采样到第i+1个采样的时间间隔;at[i]表示第i+1个采样时所求的世界坐标系中对应的ax、ay或az;σt[i]表示at[i]的测量误差。
7.根据权利要求6所述的移动终端,其特征在于,所述方位角获取单元包括:声波信号的方向矢量计算单元,用于计算所述第二终端中扬声器声波信号的方向矢量,其公式包括:
&lambda; x v x [ k ] + &lambda; y v y [ k ] = v , &ForAll; k ;
其中,λx与λy第二终端中扬声器的声波信号在X轴与Y轴的方向矢量;vx[k]与vy[k]为第k次采样计算到的X轴与Y轴的速度矢量;为任意值符号;
方位角计算单元,用于根据计算到的第二终端中扬声器声波信号的方向矢量,计算所述第二移动终端相对于所述移动终端的方向,其公式包括:
&alpha; = arcsin &lambda; &lambda; x 2 + &lambda; y 2 &lambda; x &GreaterEqual; 0 &pi; + arcsin &lambda; y &lambda; x 2 + &lambda; y 2 &lambda; x < 0 .
8.一种移动终端的测向系统,其特征在于,该系统包括:权利要求6-7任一项所述的移动终端,以及第二移动终端;
其中,所述第二移动终端,用于通过扬声器发送的指定频率的声波信号。
CN201310143217.3A 2013-04-23 2013-04-23 一种移动终端的测向方法、系统及该移动终端 Expired - Fee Related CN103267961B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310143217.3A CN103267961B (zh) 2013-04-23 2013-04-23 一种移动终端的测向方法、系统及该移动终端

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310143217.3A CN103267961B (zh) 2013-04-23 2013-04-23 一种移动终端的测向方法、系统及该移动终端

Publications (2)

Publication Number Publication Date
CN103267961A CN103267961A (zh) 2013-08-28
CN103267961B true CN103267961B (zh) 2016-07-06

Family

ID=49011600

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310143217.3A Expired - Fee Related CN103267961B (zh) 2013-04-23 2013-04-23 一种移动终端的测向方法、系统及该移动终端

Country Status (1)

Country Link
CN (1) CN103267961B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107424393B (zh) * 2017-07-31 2020-11-13 江苏知链科技有限公司 移动终端寻找方法、移动终端及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1200618A (zh) * 1997-04-25 1998-12-02 罗宏基 寻呼通讯与测向定位系统
CN101137226A (zh) * 2007-09-30 2008-03-05 重庆邮电大学 单天线终端测向方法
CN101707805A (zh) * 2009-11-19 2010-05-12 淮阴工学院 移动终端多基站频移综合定位方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4450517B2 (ja) * 1999-04-20 2010-04-14 エヌエックスピー ビー ヴィ 携帯電話機におけるアンテナ方向の検出

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1200618A (zh) * 1997-04-25 1998-12-02 罗宏基 寻呼通讯与测向定位系统
CN101137226A (zh) * 2007-09-30 2008-03-05 重庆邮电大学 单天线终端测向方法
CN101707805A (zh) * 2009-11-19 2010-05-12 淮阴工学院 移动终端多基站频移综合定位方法

Also Published As

Publication number Publication date
CN103267961A (zh) 2013-08-28

Similar Documents

Publication Publication Date Title
US9817112B2 (en) Pairwise measurements for improved position determination
CN107003378B (zh) 便携式电子装置和确定便携式电子装置的地理位置的方法
US9681267B2 (en) Positioning techniques for narrowband wireless signals under dense multipath conditions
CN103439689B (zh) 一种分布式麦克风阵列中麦克风位置估计系统
Höflinger et al. Acoustic self-calibrating system for indoor smartphone tracking (assist)
CN106154222B (zh) 一种利用无线电射频信号检测人的行走方向的方法
CN105718064A (zh) 基于超声波的手势识别系统与方法
WO2021003757A1 (zh) 一种基于信号多径传播测量的室内定位方法及系统
KR20160087815A (ko) 운전자측 위치 검출
CN103257338B (zh) 一种室内定位的方法及系统
CN109725292B (zh) 水下作业多目标高精度超短基线定位方法及装置
CN113163486B (zh) 一种基于5g tdoa的相对定位方法及系统
CN110703191A (zh) 定位方法、系统、装置、定位设备及存储介质
CN106483495A (zh) 一种室内运动标签定位和测速方法
CN109683178B (zh) 卫星导航欺骗检测方法和装置
Granados et al. Redundant localization system for automatic vehicles
CN103267961B (zh) 一种移动终端的测向方法、系统及该移动终端
Suo et al. Single acoustic vector sensor DOA enhanced by unsaturated bistable stochastic resonance with linear amplitude response constrained
EP2912408B1 (en) Non-echo ultrasonic doppler for corrected inertial navigation
CN103591959A (zh) 室内人员的定位方法
CN115333904A (zh) 一种多普勒频移的校正方法、装置、设备及介质
WO2005017552A1 (ja) 情報処理装置およびgps測位方法
CN116106888A (zh) 无接触式运动目标速度测量方法、系统、设备及存储介质
CN109959921A (zh) 一种基于北斗授时的声信号距离估计方法
JP2012202896A (ja) 位置測定システム、位置測定装置および位置測定プログラム

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20160706

CF01 Termination of patent right due to non-payment of annual fee