基于改进Sinc插值的小范围高精度定位方法
【专利摘要】本发明提供了一种基于改进Sinc插值的小范围高精度定位方法，包括：L个固定的参考节点接收目标节点发出锯齿调频波v(t)信号；对接收的信号进行幅度均衡并对锯齿调频波以时间间隔T进行采样,得到样值函数xi(n)；对样值函数进行重构,得到重构函数y(k)；对发送的锯齿调频波v(t)以时间间隔T/(M+1)采样得到样值函数u(k),把y(k)和u(k)进行互相关运算得相关峰值位置；利用相关峰值位置的差分别计算信号到达各参考节点和到达第1个参考节点的时间差，得到TDOA测量值；将TDOA测量值和参考节点的坐标代入Chan氏算法计算目标节点的位置；求目标节点的定位计算结果和真实位置的均方根误差得到定位精度RMSE。本发明通过降低采样频率，保证了定位的高精度。
【专利说明】基于改进Sinc插值的小范围高精度定位方法
【技术领域】
[0001]本发明涉及无线通信领域，尤其涉及一种基于改进Sinc插值的小范围高精度定位方法。
【背景技术】
[0002]全球定位系统(GPS)的使用大幅度的提高了无线定位精度，使定位精度可达到米级，但是在小范围环境内GPS定位效果并不好，米级的定位精度远远不足。在一些精度要求比较高的应用系统中，现有定位系统往往不能满足用户的要求。如机场大厅、超市、图书馆、停车场等复杂的小范围环境，由于受到室内环境等条件的影响，目前的定位技术在精度上往往达不到需求。如何在小范围定位系统中提高定位精度是各国学者研究的热点。在小范围条件下对于复杂环境而言，室内布置，材料结构，建筑物尺度的不同导致了信号的路径损耗很大，与此同时，小范围复杂的环境会引起信号的反射，绕射，折射和散射，形成多径现象，使得接收信号的幅度，相位和到达时间发生变化，造成信号的损失，定位的难度大。相比于GPS系统，小范围的定位系统具有比较鲜明的特点与优势。小范围的定位系统其精度可达到厘米级，而且对于复杂的环境具有较强的适应性，应用范围更加的广阔。而且随着嵌入式技术、通信技术、信息处理技术的发展和成熟，大量适用于复杂环境的小型无线设备的出现也给小范围无线定位的研究和应用带来了方便。因此，小范围定位技术在近几年逐渐成为人们关注的一个研究方向，在将来也会继续成为研究的热点。
[0003]小范围无线定位技术通过对接收到的无线电波的一些参数进行测量，根据待定的算法判断出被测物体的位置，其主要机理是根据接收到的定位信号的强度(SOA =Strengthof Arrival)、信号到达的角度(AOA:Angle of Arrival)、信号到达的时间(Τ0Α:Time ofArrival)和信号到达的时间差(TDOA:Time Difference of Arrival)以及这些参数的不同组合来进行移动终端的定位过程。到达时间差的定位技术TDOA是目前无线定位系统中应用最广泛的一项技术。与TOA定位技术相比，它不要求参考节点和目标节点之间的精确同步，与AOA定位技术相比，它不需要特殊的天线阵列，因而更加经济。由于SOA技术受周围环境影响很大，因此定位精度不及TODA高。在TODA技术中，主要有基于最小二乘算法(LS算法)，Chan氏算法,泰勒级数展开法等定位算法。当对定位精度要求不是很高的时候，需要对目标节点位置做初始估计时，采用TDOA最小二乘算法。Chan氏算法是对LS算法的一种改进，先采用加权最小二乘算法(WLS)得到一个初始解，通过利用第一次得到的估计位置坐标及附加变量等已知的约束条件来进行第二次WLS估计能得到改进的估计位置。
[0004]TDOA作为一种无线定位技术，是基于时间测量的定位方法。它实现的原理可分为两个基本步骤:1、利用时延估计技术将不同参考节点所接收到的目标节点发出的信号到达时间差进行精确估计，同时利用距离和时间的关系构造非线性双曲线方程；2、将估计出的参数带入方程，通过适当的算法求解方程，得到移动台位置参数。TDOA技术主要是用于构建定位的基础数学模型。TDOA是通过检测信号到达两个基站的时间差，而不是到达的绝对时间来确定移动台的位置，降低了时间同步要求。采用三个不同的基站可以测到两个TD0A，移动站位于两个TDOA决定的双曲线的交点上。
[0005]传统的基于最小二乘算法的TDOA定位技术实现流程如附图1所示。其最小二乘算法为:假设待定位的目标节点坐标为(X，y)，(Xi, Yi)为第i个参考节点的已知位置，共L个参考节点，如果要实现对移动台的二维坐标的定位，至少需要三个基站参与TDOA的测量，也即获得两个TDOA测量值，从而构成非线性方程组:
[0006]
【权利要求】
1.一种基于改进Sinc插值的小范围高精度定位方法，其特征在于，包括: 步骤一:L个固定的参考节点接收目标节点发送的锯齿调频波v(t)的信号，其中L为正整数； 步骤二:对接收的信号进行幅度均衡，并对M个连续的锯齿调频波以时间间隔T进行采样，得到样值函数Xi (η)，其中i = 1，2，*..,Μ ；n = O, I,..., N_l, N为每个锯齿调频波的样值点个数，其中M为正整数； 步骤三:利用Sinc插值算法对样值函数进行重构，得到重构函数y (k)，其中k = O，1，…，(N-1) (M+l)+M。 步骤四:对发送的锯齿调频波v(t)直接以间隔时间T/(M+1)采样，得到样值函数u(k),其中k = 0，I,..., (N-1) (M+l)+M，把重构函数y(k)和u(k)进行互相关运算得到相关峰值位置Ai，其中i取1，2，3...?; 步骤五:利用相关峰值位置的差分别计算信号到达第2，3，4，5，…，L个参考节点和到达第I个参考节点的时间差t2，1； t3；1, -, ti； i，…\ i，得到TDOA测量值，其中ti; i =(A1-A1)*!1, T为y(k)的样值点时间间隔； 步骤六:将TDOA测量值,和参考节点的坐标(X1J1), (x2, y2),…，(XiJi),…，(xL, yL)带入Chan氏算法计算目标节点的位置； 步骤七:求目标节点的定位计算结果和真实位置(x，y)的均方根误差得到定位精度RMSE。
2.如权利要求1所述的基于改进Sinc插值的小范围高精度定位方法，其特征在于，所述时间间隔T的选择满足:1/T>2fv，其中fv为v(t)的信号带宽，fv的计算公式如下:
fv = 2(mv+l)Fm 式中mv为v(t)的调制指数，Fm为的锯齿调制信号的最大频率分量，Fm = 7F，其中F为锯齿波峰出现的频率。
3.如权利要求2所述的基于改进Sinc插值的小范围高精度定位方法，其特征在于，步骤三具体包括: A:用Sinc插值算法对步骤二得到的样值函数Xi (η)进行重构，得到重构函数(h)，i=1,2,…，M;h = 0，l，…，(N-1) (M+l)+M，其中 N 为正整数； Xi (n), (i = 1,2,..., M ；n = 0,2,..., N_l)为以T为间隔经过采样得到的已知序列，在两个样值点之间插入M个点，构成长度为M的序列Ci (m, n), (i，m = 1,2,..., M ；n = O,I,..., N_l),则可得: f 、 X-' /nsin[(/i + w/(M +1)-/)λγ], ^、 c, (m, η) = 2.^仍^-//1# ? n^( 5 )
(rt + m /(M + 1)- 1}π 把样值函数Xi (η)和样值点之间的M个插值点结合，得到重构函数7101)；
Ii (O) = Xi (O)，Ji (M+1) = Xi (I)，…，Ji ((N-1) (M+l)) = Xi (N-1)，
Ii ⑴=C1 (i, O)，Ji (M+2) = C1 (i, I)，…，Ji ((N-1) (M+1) +1) = C1 (i, N-1)，
4.如权利要求3所述的基于改进Sinc插值的小范围高精度定位方法,其特征在于,步骤六具体包括: A:TD0A构造构成非线性方程组: 假设待定位的目标节点坐标为(X，y)，(Xi, Yi)为第i个参考节点的已知位置，共L个参考节点，获得两个TDOA测量值，构成非线性方程组:
进行线性化处理:
(6)式中 xirl = X1-X1, yul = Y1-Y1, Ki =xf +yf+J12 ,且 i 取 2，3---L ;x，y, I11 是未知数，(7)式为线性方程组； 利用Chan氏算法求解目标节点的坐标； B:Chan氏算法求解第一次加权平均:令23为未知矢量，其中ζρ = [X，y]T，假定Za的元素间相互独立，则Za的ML估
进一步的近似式为: VkGlQMj'GlQ、 利用所述近似式得到初始解，利用所得到的初始解计算B矩阵，再利用2fl =(Gy..4Ga)...' Gy-1Z1求得第一次的WLS计算结果； C:Chan氏算法求解第二次加权平均求出定位结果: 进行第二次WLS计算，首先计算估计位置Za的协方差矩阵；其协方差矩阵为:cov(ru) = (Gy 1G11) 1 ;矢量Za是均值为实际值、协方差矩阵由aw(2,) = (Gj Ψ iG, 1确定的随机矢量； 求得方程组:11^ =h, -G' V Λ,
定义Ψ, Sza的误差矢量；Ψ,的协方差矩阵为Ψ, =[Ψ' Ψ, T]=4B' cov(za)B'， 其中，
B'= O /-V1 O ，B'中的x°，y°，^，即为第一次的WLS计算结果;
ψ'为高斯分布，Z' a的最大似然估计为:4=(0,【ψΗ0(:),求得目标节定位计算结果为
求得目标节点和第i参考节点与第I个参考节点的距离差ri;1,其公式为= ti；1*c,其中c为电磁波的传播速度和(Xi, Yi)代入到Chan氏算法中，求解非线性方程组式(6)，获得目标节点的坐标(x，y)，其中(Xi，yi)为第i参考节点的坐标。
【文档编号】G01S5/06GK104198991SQ201410390828
【公开日】2014年12月10日 申请日期:2014年8月10日 优先权日:2014年8月10日
【发明者】李争平, 王景中, 曾勇波, 秦朝亮 申请人:北方工业大学