CN107247961A - 一种应用模糊轨迹序列的轨迹预测方法 - Google Patents
一种应用模糊轨迹序列的轨迹预测方法 Download PDFInfo
- Publication number
- CN107247961A CN107247961A CN201710324975.3A CN201710324975A CN107247961A CN 107247961 A CN107247961 A CN 107247961A CN 201710324975 A CN201710324975 A CN 201710324975A CN 107247961 A CN107247961 A CN 107247961A
- Authority
- CN
- China
- Prior art keywords
- dimension
- sequence
- user
- sample
- trajectory
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23211—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with adaptive number of clusters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N7/00—Computing arrangements based on specific mathematical models
- G06N7/02—Computing arrangements based on specific mathematical models using fuzzy logic
Abstract
本发明提出一种应用模糊轨迹序列的轨迹预测方法,通过引入模糊轨迹序列模型,对用户的移动轨迹进行模糊化处理,并通过模糊时间序列的预测方法,对用户的运动轨迹进行预测。方法中对用户的轨迹进行了模糊化处理,设计了网格化的模糊化方法,由此实现了模糊轨迹序列模型。通过引入离群点处理的机制,本方法可以处理历史信息不足时的轨迹预测问题,通过引入离群点检测机制以及离群点预测模型,以降低离群点对算法整体性能的影响。
Description
技术领域
本发明涉及用户轨迹序列预测问题,特别涉及一种引入了模糊轨迹序列模型的轨迹预测方法。
背景技术
1、轨迹预测
随着定位技术的发展及广泛应用,基于位置的服务(Location-based Service,LBS)逐渐成为生活中不可获缺的部分。通过分析用户的轨迹信息,挖掘其中隐含的用户信息来完善用户的服务体验,成为数据挖掘的一个重要领域。通过用户的轨迹信息对用户的轨迹进行预测,对导航服务、交通管理以及基于位置的广告投放等应用至关重要,成为了轨迹挖掘的一个研究热点。
2、轨迹预测的常见方法
轨迹预测问题主要分为长期预测以及短期预测:长期预测适用于时间跨度长、采样间隔大的轨迹样本预测问题,用户的移动目标对用户轨迹起主要决定作用,主要通过历史用户轨迹的挖掘,发掘用户相似轨迹或规律用于样本的预测;短期预测则适用于时间跨度短、采样间隔小的轨迹样本预测问题,用户的移动习惯对用户轨迹起主要决定作用,主要通过对短期用户行为模式进行分析建模完成预测。
轨迹预测的方法可以分为两类,分别是统计类方法以及机器学习类方法。统计类方法主要通过数学领域的统计模型以及概率模型对用户轨迹进行建模,建立的模型简单,算法运行效率高,但因此难以有效的处理现实中复杂的数据集;机器学习类方法通过应用机器学习领域的预测技术,挖掘轨迹信息建立用户轨迹模型进行预测,能够建立起复杂的用户轨迹模型,但算法效率低,算法运行时间长。
现有的轨迹预测方法都立足于所得到的用户轨迹是大致精确的这一前提,不论是统计学的预测方法还是机器学习的方法,均对样本的数据不做过多的预处理而直接使用,而在实际生活中通过定位服务所得到的轨迹序列均包含误差,这种误差影响了现有方法的预测精度。
发明内容
本发明的目的在于针对不精确、存在偏移的用户轨迹序列的预测问题,提出一种应用模糊轨迹序列的轨迹预测方法,该方法通过引入模糊轨迹序列对用户的轨迹进行预测,使得预测方法能够处理不精确的、存在坐标偏移的轨迹序列,能够更广泛的应用于现实生活中的轨迹预测。
为了达到上述目的,本发明采用了以下技术方案:
1)针对任意一待预测用户,确定用户轨迹训练样本集的论域,所述用户轨迹训练样本为所述论域范围内若干个用户的轨迹序列,将用户轨迹训练样本集中的用户轨迹序列以截取的方式重构为若干个包含k个时刻用户位置坐标的高维轨迹序列,根据论域的网格化划分方案将所述高维轨迹序列模糊化为高维模糊轨迹序列,且满足使相同或相似的高维轨迹序列在模糊化后仍为相同或相似的高维模糊轨迹序列,得到用于高维模糊轨迹序列模型训练的集合;
2)在每一个高维模糊轨迹序列后加入下一刻对应用户的位置坐标,得到高维模糊轨迹序列模型的训练集,对高维模糊轨迹序列模型的训练集应用聚类方法获得簇,并由簇中心构建相应的模糊逻辑规则,依据模糊逻辑规则对由用户轨迹训练样本集中的用户轨迹序列重构得到的高维轨迹序列进行划分,对于任意一个高维轨迹序列:若该高维轨迹序列不满足模糊逻辑规则,则将该高维轨迹序列加入离群点样本训练集合,利用离群点样本训练得到离群点预测模型;
3)在所述论域范围内,将某一个待预测用户截止当前时刻的顺序采集的k个位置坐标转变为相应的高维轨迹序列,得到待预测样本,对待预测样本根据所述模糊逻辑规则进行划分,若待预测样本为满足模糊逻辑规则的一般样本,通过对簇中心最后S维分别求加权平均数获得用户轨迹预测结果,所述S的取值与用户位置坐标维数相同;若待预测样本为不满足模糊逻辑规则的离群点样本,则利用离群点预测模型进行轨迹预测。
所述步骤1)具体包括以下步骤:
1.1)计算所述轨迹训练样本集中用户坐标经、纬度的最值,则用户活动范围定义为经、纬度的最值所围成的矩形区域,根据该矩形区域确定用户轨迹训练样本集中用户轨迹序列的论域;
1.2)将所述论域划分为网格状区域,每个网格状区域对应不同的模糊集;将所述轨迹训练样本集中全部用户轨迹序列转化为高维轨迹序列;将所有高维轨迹序列中坐标模糊化为对应的模糊集,形成高维模糊轨迹序列集合;
1.3)对具有某一相同形式的高维模糊轨迹序列,计算这些高维模糊轨迹序列对应的模糊化前的高维轨迹序列中下一刻用户坐标的方差;将计算得到的所有不同形式的高维模糊轨迹序列对应的方差求和,作为对论域划分方案的评价指标;
1.4)以最小化所述评价指标为目标,利用启发式优化算法找出尽可能好的论域划分方案,作为最终的论域划分方案;
1.5)利用最终的论域划分方案,按照步骤1.2)生成用于高维模糊轨迹序列模型训练的集合。
所述步骤1.2)中,对论域进行网格状划分过程中,分别在经、纬度方向上选取划分点,并根据划分点对论域进行划分,在选择划分点时,对划分点的数目通过设置上限和下限进行限制。
选择等分方式对论域进行划分,在经、纬度方向上的划分点的数目均以指数形式增加,对由此形成的每个划分方案进行评价,选择使得评价指标最小的划分方案作为参考方案,计算上述每个划分方案的评价指标与参考方案评价指标之间差的绝对值作为对应划分方案的相对评价指标;在上述所有划分方案中,选择相对评价指标大于设定阈值且划分点比参考方案划分点数目少的方案中,划分点数目最多的方案的划分点数目作为划分点数目的下限,选择相对评价指标大于设定阈值且划分点比参考方案划分点数目多的方案中,划分点数目最少的方案的划分点数目作为划分点数目的上限。
所述步骤2)中,模糊逻辑规则的构建具体包括以下步骤:
2.1)对步骤1)中生成的用于高维模糊轨迹序列模型训练的集合中每个高维模糊轨迹序列,在序列后新增一维,值为该序列对应的用户轨迹序列中下一刻用户的2维平面位置坐标,生成用于聚类的高维模糊轨迹序列训练集;
2.2)对步骤2.1)中生成的用于聚类的高维模糊轨迹序列训练集,应用自适应高斯核的模糊C均值聚类算法生成高维轨迹序列簇,并计算簇中心作为基于一般样本的轨迹预测模型的模糊逻辑规则。
所述步骤2.2)中,自适应高斯核的模糊C均值聚类算法在完成聚类算法的每一次迭代后,依据现有样本相对于簇中心的分布情况,重新计算簇的核宽度并更新。
所述模糊逻辑规则定义为:假设x(t)是由F(t-k),F(t-k+1),…,F(t-1)决定的,且F(t)代表t时刻用户坐标对应的模糊集,x(t)为t时刻的用户坐标,那么就存在模糊逻辑规则F(t-k),F(t-k+1),...,F(t-1)→x(t)。
计算高维轨迹序列与所有簇间的隶属度,若所有隶属度均小于设定阈值η,则该样本不满足模糊逻辑规则;若存在大于等于设定阈值η的隶属度,则样本满足模糊逻辑规则。
所述离群点预测模型的建立方法包括以下步骤:将离群点样本对应的用户轨迹序列转化为间隔时间的位移序列以及位移方向的序列,然后采用Levenberg-Marquardt算法分别构建关于下一个间隔时间的位移以及位移方向的预测模型;Levenberg-Marquardt算法需要通过多次迭代来避免陷入局部极小,每次迭代均在设定范围内随机生成初始权向量,并选择最优的结果;Levenberg-Marquardt算法的初始权向量设定时,每一维度的值大于前一维度。
所述步骤3)中,对簇中心最后S维求加权平均数时,若待预测样本与簇间隶属度大于等于阈值η,则权值为该隶属度;否则,权值为0。
本发明具有以下有益效果:
本发明通过引入模糊轨迹序列模型,并通过模糊时间序列的预测方法,对用户的运动轨迹进行预测。方法中对用户的移动轨迹进行了模糊化处理,设计了网格化的模糊化方法,由此实现了模糊轨迹序列模型。通过引入离群点处理的机制,本方法可以处理历史信息不足时的轨迹预测问题,并引入离群点检测机制以及离群点预测模型,以降低离群点对算法整体性能的影响。本发明中通过引入模糊轨迹序列模型对用户轨迹序列进行建模,使得模型可以处理不精确的、存在坐标偏移的轨迹序列,为处理现实生活中定位服务存在的偏差的问题提供了新的解决思路。
进一步的,本发明采用离群点判定机制及离群点预测模型,可以识别训练样本不足导致的离群点这一特殊情况,同时通过对预测模型的调整,改善对离群点样本的预测精度,提升算法的鲁棒性以及预测的平均精度。
进一步的,本发明中采用的聚类方法区别于传统的核版本的模糊C均值聚类,采用动态调整的核宽度策略,同时引入核宽度估算机制,使得聚类方法能自适应于不同尺寸的簇,提升了预测算法的鲁棒性。
附图说明
图1为应用模糊轨迹序列的轨迹预测方法流程图;
图2为对用户轨迹论域的划分方案及模糊集示意图;
图3为论域划分方案搜索流程图;
图4为应用模糊轨迹序列的轨迹预测方法中模糊逻辑规则构建流程图;
图5为应用离群点预测模型的轨迹序列表征方式转换示意图;
图6为应用模糊轨迹序列的轨迹预测方法对南加州展会数据集的仿真结果。
具体实施方式
为了使本发明的内容、效果以及优点更加清楚明白,下面结合附图和实施例对本发明进行详细描述。
本发明是基于模糊逻辑的预测方法,通过对轨迹数据进行模糊化预处理,使得预测方法可以处理获得的不准确的、有偏差的轨迹预测问题。本发明设计了一种自适应的高斯核模糊C均值聚类方法,通过引入核宽度的调整算法,使聚类方法能够更适用于更为广泛的各种情况。同时本发明采用离群点模型应对训练样本稀疏引发的问题。图1中展示了本发明的策略流程,通过计算待预测样本与每个提取到的模糊逻辑规则对应的簇间隶属度,判定样本是否为离群点,对一般性样本采用加权平均数进行预测,而对离群点样本采用线性模型预测。该方法优势在于能够应对现实生活中轨迹预测问题,可以容忍数据采集中出现的偏差,同时也适用于训练样本数量不足、覆盖不全面的情况。
(一)定义用户轨迹论域划分方案及轨迹序列的模糊化处理
用户轨迹的论域定义了用户的活动范围以及用户可能到达的区域,为了使得定义的论域能够覆盖用户所有的轨迹,计算用户经度坐标的最值xmin、xmax以及纬度坐标的最值ymin、ymax,参见图2,用户的活动范围可以用计算得到的经纬度坐标最值围成的矩形区域[(xmin,ymin),(xmin,ymax),(xmax,ymin),(xmax,ymax)]来表示,考虑到用户可能到达的区域,将用户的活动范围扩大为[(x'min,y'min),(x'min,y'max),(x'max,y'min),(x'max,y'max)],其中:
Δx=0.1·(xmax-xmin)
Δy=0.1·(ymax-ymin)
x'min=xmin-Δx
x'max=xmax+Δx
y'min=ymin-Δy
y'max=ymax+Δy
该区域[(x'min,y'min),(x'min,y'max),(x'max,y'min),(x'max,y'max)]定义为用户轨迹的论域。
参见图3,对用户轨迹的论域进行划分包括以下步骤:
1)采用等时间间隔方式获取用户的位置坐标,构成用户轨迹序列。设定时间间隔t,对于用户的一次轨迹采样,每隔时间t收集一次位置坐标,并将坐标序列存储作为用户的一条轨迹。对于记录的用户的一条含有n个采样点的轨迹,第i个样本点对应的用户坐标为(xi,yi),其中,xi为经度坐标,yi为纬度坐标,则可得到:
xmin=min(x1,x2,...,xn)
xmax=max(x1,x2,...,xn)
ymin=min(y1,y2,...,yn)
ymax=max(y1,y2,...,yn)
由此可确定用户轨迹的论域[(x'min,y'min),(x'min,y'max),(x'max,y'min),(x'max,y'max)]。
2)参见图3,对于用户轨迹的论域[(x'min,y'min),(x'min,y'max),(x'max,y'min),(x'max,y'max)],分别在经度和纬度范围内选取各自的划分点和由这些划分点确定论域划分的网格。每个网格对应一个模糊集Si(i=1,2,...,(u+1)·(v+1)),这些模糊集构成一种用户轨迹论域划分方案。
3)对用户的一条轨迹序列T=[(x1,y1),(x2,y2),...,(xn,yn)],对每个时刻,将该时刻用户位置坐标与前k-1个连续时刻的位置坐标合成一个k维轨迹序列,将一条轨迹序列转化为n-k+1条k维轨迹序列:
T1=[(x1,y1),(x2,y2),...,(xk,yk)]
T2=[(x2,y2),(x3,y3),...,(xk+1,yk+1)]
...
Tn-k+1=[(xn-k+1,yn-k+1),(xn-k+2,yn-k+2),...,(xn,yn)]
4)对得到的每条k维轨迹序列Ti=[(xi,yi),(xi+1,yi+1),...,(xi+k-1,yi+k-1)],将每个轨迹中每个坐标点替换为所在网格区域对应的模糊集,得到k维模糊轨迹序列:
其中,f(t)代表t时刻用户坐标所在的模糊集下标,t=i,i+1,…,i+k-1。
5)对所得到的划分方案进行评估,并通过利用启发式优化算法寻找到一个可行划分方案。
6)根据得到的可行的论域划分方案,执行步骤3)和步骤4)生成k维模糊轨迹序列,这些序列构成了用于高维模糊轨迹序列模型训练的集合。
在步骤2)中,划分点的数目需要满足设置的上、下限,划分点数目上、下限的确定方法为:利用等分的方法,依次在论域的经度和纬度范围内选择1、10、…、10b个等分划分点;并分别根据步骤2)中描述生成等分用户轨迹论域划分方案,其中b取值为使得等分后的每个区间长度小于数据集的分辨率的最小值;针对每个等分用户轨迹论域划分方案,执行步骤3)和步骤4),并对这些等分用户轨迹论域划分方案进行评价,并给出评价指标值;假设包含10l个等份划分点的论域划分方案的评价指标值为el,l=0,1,2,…,b,且所有方案中评价指标的最小值为ek:则对于l<k,若满足且对于任意整数p(l<p≤k),则10l为划分点数目的下限;而对于l>k,若满足且对于任意整数p(k≤p<l),则10l为划分点数目的上限;的取值范围为(0,1],优选为0.2。
(二)用户轨迹论域划分方案的评价标准及可行方案的搜索方法
良好的用户论域划分方案,应满足使相同或相似的k维用户轨迹序列,模糊化后仍为相同或相似的模糊轨迹序列。因而对一种用户轨迹论域的划分方案评价过程包括以下步骤:
1)针对(一)中步骤6)得到的用于高维模糊轨迹序列模型训练的集合,将其中形式相同的k维模糊轨迹序列归入同一个集合;
2)对其中所有具有和相同形式的k维模糊轨迹序列所构成的集合计算集合中所有序列对应的模糊化前的高维(k维)轨迹序列的下一刻用户位置的方差,即计算集合的方差,作为形式对应的方差;
3)针对用于高维模糊轨迹序列模型训练的集合中所有不同的k维模糊轨迹序列形式,利用步骤2)计算该形式对应的方差,并求得所有的方差之和;当得到的方差之和越小,则证明划分方案性能越好。
参见图3,采用启发式优化算法搜索得到最佳划分方案,为了建立问题的解空间,将用户轨迹的经纬度范围分别等分为由u'、v'个划分点隔开的小区域,其中,u'>>u且v'>>v。则问题转化为分别在u'个经度划分点和v'个纬度划分点中分别选取u、v个划分点的问题,通过启发式算法搜索该问题的最优解,并作为最终采用的划分方案。
(三)应用用户k维模糊轨迹序列构建k维模糊逻辑规则
k维模糊逻辑规则定义为:假设x(t)是由F(t-k),F(t-k+1),…,F(t-1)决定的,且F(t)代表t时刻用户坐标对应的模糊集,x(t)为t时刻的用户坐标,那么就存在模糊逻辑规则F(t-k),F(t-k+1),...,F(t-1)→x(t)。
k维模糊逻辑规则构建方法包括以下步骤:
1)对所有的k维用户轨迹序列{Ti|i=1,2,...,n-k+1},应用所选的论域划分方式,进行模糊化处理,得到k维模糊轨迹序列集
2)对任意一条k维模糊轨迹序列将该k维模糊轨迹序列后加入下一刻用户位置(xj+k,yj+k),构成k+2维模糊轨迹序列:
T'j=[Sf(j),Sf(j+1),...,Sf(j+k-1),xj+k,yj+k]
将所有k维模糊轨迹序列转化为这样的k+2维模糊轨迹序列并构成集合。
3)对所有的k+2维模糊轨迹序列,应用自适应高斯核的模糊C均值聚类算法,获得聚类结果,得到的每个簇中心即为一条k+2维模糊逻辑规则。
由于模糊C均值聚类算法容易陷入局部极小,所以最终得到的结果受初始给定的簇中心影响,为了消除这种影响而得到性能较好的k+2维模糊逻辑规则,需要迭代该过程,每次均随机给予初始聚类中心,并最终选择性能最好的结果,并由此最终构建得到k+2维模糊逻辑规则;对得到的聚类结果进行性能评价的方法为,利用聚类结果构建的k+2维模糊逻辑规则对训练集中的样本进行预测,并依据预测结果的均方根误差来判定性能,均方根误差越小,代表所得到的模型性能越好。
(四)自适应高斯核模糊C均值聚类算法流程
参见图4,自适应高斯核模糊C均值聚类算法通过引入高斯核函数作为隶属度函数,一方面减小了离群点样本对簇本身的影响,另一方面提升了簇的鲁棒性;通过自适应的核宽度设置,使得到的簇能够更好的适应不同的样本分布。算法包括以下步骤:
1)迭代k+2次,第h(0<h≤k)次迭代中,在用户轨迹论域范围内随机生成rh个值,后两次迭代中,分别在用户轨迹论域的经度和纬度范围内生成rk+1和rk+2个值,并组合成个簇中心;对第i个簇中心,为其核宽度赋一个初始值
2)通过模糊C均值聚类方法迭代过程迭代一次,迭代过程中,采用簇对应核宽度下的高斯核函数作为该簇的隶属度函数;在第h次迭代过程中,k+2维模糊轨迹序列T'j对簇中心Ci对应的簇的隶属度可由下式计算:
其中||T'j-Ci||表示T'j与Ci间的欧式距离;
3)利用Parzen窗估计法,对每个得到的簇,根据样本关于簇中心的分布情况,计算新的核宽度并更新;
4)迭代直至聚类结果收敛,并给出包含k+2维模糊轨迹序列的簇中心的聚类结果。
由于模糊轨迹序列中包含模糊集,在计算模糊化后的样本对簇的隶属度以及样本间距离时,定义在计算含有模糊轨迹序列的欧式距离时,将模糊集用模糊集对应的矩形区域的中心坐标值来代替。
(五)离群点样本判定及离群点预测模型构建
k维模糊逻辑规则构建完成后,对训练集中样本Ti,计算样本对得到的簇中心Cj对应的簇的隶属度,隶属度函数采用高斯核函数memk(Ti,Cj),核宽度为聚类结果给出的簇中心Cj的核宽度,由于Ti为k维向量,因而memk(Ti,Cj)的计算过程中仅考虑Cj的前k维数据;对任意样本Ti,若存在簇中心Cj,使得memk(Ti,Cj)≥η,则表明样本可以被所构建的k维模糊逻辑规则所表征;若对于任意簇中心Cj,memk(Ti,Cj)<η,则样本不被所构建的k维模糊逻辑规则所表征,样本不适合于一般样本预测模型,将样本加入离群点样本集合。
对由训练集中得到的离群点样本集合,囿于样本数目不足,因此采用简单的线性模型进行学习训练,使用Levenberg-Marquardt算法得到离群点的预测模型;利用Levenberg-Marquardt算法构建模型时,包括以下步骤:
1)参见图5,将用户轨迹序列T=[(x1,y1),(x2,y2),...,(xn,yn)]转化为位移与位移方向表征的序列:
Tp=[(ρ2,θ2),(ρ3,θ3),...,(ρn,θn)]
其中,(ρi,θi)为(xi,yi)关于以(xi-1,yi-1)为极点以x轴正方向为极轴的极坐标系表示,i=2,3,…,n;
2)对离群点样本集合中任意样本T'i,生成离群点位移预测模型训练样本[ρi+1,ρi+2,...,ρi+k]以及位移方向预测模型训练样本[θi+1,θi+2,...,θi+k],并分别构成位移预测模型训练集和位移方向预测模型训练集;
3)采用Levenberg-Marquardt算法构建预测模型,即分别在离群点位移预测模型训练集以及位移方向预测模型训练集上进行训练,得到移预测模型和位移方向预测模型;两个预测模型得到的预测结果可以运算出下一刻用户位置的预测值。
由于Levenberg-Marquardt算法会陷入局部极小,因而需要随机给出符合条件的多个初始参数,随机生成的初始参数需要符合的条件为:对于生成的初始参数,变量的每一维度的值必须大于前一维度的值。迭代找到使得Levenberg-Marquardt算法自身目标函数取得最小值的初始参数,用来构建离群点预测模型。
(六)利用模糊轨迹序列模型预测用户下一刻位置
对用户轨迹序列进行预测,即是根据用户前k时刻用户位置坐标,预测下一刻用户的位置。则可得到待预测k维用户轨迹序列对得到的k维模糊逻辑规则对应的簇中心Ci,计算待预测样本得到的k维轨迹序列与每个簇间隶属度memk(Tp,Ci),若存在i,使得memk(Tp,Ci)≥η,则为一般样本,则对各个簇中心最后2维分别求加权平均数,得到用户下一刻位置,其中,若待预测样本与簇间隶属度大于等于阈值η,则权值为该隶属度;否则,权值为0:
其中,r表示簇中心的个数,Ci(k+1)、Ci(k+2)表示簇中心Ci的第k+1、k+2维值;若不存在i,使得memk(Tp,Ci)≥η,待预测样本视为离群点样本,Tp转化为k-1维的用户位移序列以及位移角度序列,根据Levenberg-Marquardt算法得到的参数分别求得用户下一刻的位移ρp以及位移角度θp,用户下一刻位置计算为:
为第k时刻用户位置坐标;
(七)仿真实验
参见图6,对于用户的轨迹数据,采用奥兰多迪士尼、纽约以及南加州展会上采集到的用户轨迹数据集,数据采样时间间隔为30秒一次,3个数据集分别包含41条、39条、19条用户轨迹,分别包含44386个、39554个、5861个用户坐标,用户坐标以某固定点为原点的直角坐标系表示,单位为英尺。
算法分别在奥兰多、纽约以及南加州展会上的3个用户轨迹数据集上运行,对于每个数据集上,选取其中每条轨迹前90%的用户轨迹序列作为学习算法的训练集,选取剩余的10%的用户轨迹作为算法验证的测试集,算法在每个数据集间独立运行。
仿真参数设计如表1所示:
表1.仿真参数
本发明应用模糊轨迹序列模型的轨迹预测方法,例如,在奥兰多用户数据集上,对测试集样本的预测均方根误差为13.2746英尺、13.4603英尺、13.9036英尺。
经过多次仿真实验,用户轨迹序列维度k的取值范围为3≤k≤10,k<3时,预测误差较大,当k>10时,运算复杂度较高,且会引起过拟合问题;离群点判定阈值η的取值范围为0.1≤η≤0.6,阈值过小则会无法有效识别离群点,阈值过大会降低模型的适用范围。
本发明通过模糊时间序列的预测方法,对用户的运动轨迹进行预测,方法中对用户的轨迹进行了模糊化处理,设计了网格化的模糊化方法,由此实现了模糊轨迹序列模型,使得预测方法能够处理现实中获取的不精确的、有偏差的用户轨迹信息。通过引入离群点处理的机制,本方法可以处理历史信息不足时的轨迹预测问题,并引入离群点检测机制以及离群点预测模型,以降低离群点对算法整体性能的影响,提升了算法的总体性能,改善了算法的鲁棒性。
Claims (10)
1.一种应用模糊轨迹序列的轨迹预测方法,其特征在于:该预测方法包括以下步骤:
1)针对任意一待预测用户,确定用户轨迹训练样本集的论域,所述用户轨迹训练样本为所述论域范围内若干个用户的轨迹序列,将用户轨迹训练样本集中的用户轨迹序列以截取的方式重构为若干个包含k个时刻用户位置坐标的高维轨迹序列,根据论域的网格化划分方案将所述高维轨迹序列模糊化为高维模糊轨迹序列,且满足使相同或相似的高维轨迹序列在模糊化后仍为相同或相似的高维模糊轨迹序列,得到用于高维模糊轨迹序列模型训练的集合;
2)在每一个高维模糊轨迹序列后加入下一刻对应用户的位置坐标,得到高维模糊轨迹序列模型的训练集,对高维模糊轨迹序列模型的训练集应用聚类方法获得簇,并由簇中心构建相应的模糊逻辑规则,依据模糊逻辑规则对由用户轨迹训练样本集中的用户轨迹序列重构得到的高维轨迹序列进行划分:对于任意一个高维轨迹序列,若该高维轨迹序列不满足模糊逻辑规则,则将该高维轨迹序列加入离群点样本训练集合,利用离群点样本训练得到离群点预测模型;
3)在所述论域范围内,将某一个待预测用户截止当前时刻的顺序采集的k个位置坐标转变为相应的高维轨迹序列,得到待预测样本,对待预测样本根据所述模糊逻辑规则进行划分,若待预测样本为满足模糊逻辑规则的一般样本,通过对簇中心最后S维分别求加权平均数获得用户轨迹预测结果,所述S的取值与用户位置坐标维数相同;若待预测样本为不满足模糊逻辑规则的离群点样本,则利用离群点预测模型进行轨迹预测。
2.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述步骤1)具体包括以下步骤:
1.1)计算所述轨迹训练样本集中用户坐标经、纬度的最值,则用户活动范围定义为经、纬度的最值所围成的矩形区域,根据该矩形区域确定用户轨迹训练样本集中用户轨迹序列的论域;
1.2)将所述论域划分为网格状区域,每个网格状区域对应不同的模糊集;将所述轨迹训练样本集中全部用户轨迹序列转化为高维轨迹序列;将所有高维轨迹序列中坐标模糊化为对应的模糊集,形成高维模糊轨迹序列集合;
1.3)对具有某一相同形式的高维模糊轨迹序列,计算这些高维模糊轨迹序列对应的模糊化前的高维轨迹序列中下一刻用户坐标的方差;将计算得到的所有不同形式的高维模糊轨迹序列对应的方差求和,作为对论域划分方案的评价指标;
1.4)以最小化所述评价指标为目标,利用启发式优化算法找出尽可能好的论域划分方案,作为最终的论域划分方案;
1.5)利用最终的论域划分方案,按照步骤1.2)生成用于高维模糊轨迹序列模型训练的集合。
3.根据权利要求2所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述步骤1.2)中,对论域进行网格状划分过程中,分别在经、纬度方向上选取划分点,并根据划分点对论域进行划分,在选择划分点时,对划分点的数目通过设置上限和下限进行限制。
4.根据权利要求3所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:选择等分方式对论域进行划分,在经、纬度方向上的划分点的数目均以指数形式增加,对由此形成的每个划分方案进行评价,选择使得评价指标最小的划分方案作为参考方案,计算上述每个划分方案的评价指标与参考方案评价指标之间差的绝对值作为对应划分方案的相对评价指标;在上述所有划分方案中,选择相对评价指标大于设定阈值且划分点比参考方案划分点数目少的方案中,划分点数目最多的方案的划分点数目作为划分点数目的下限,选择相对评价指标大于设定阈值且划分点比参考方案划分点数目多的方案中,划分点数目最少的方案的划分点数目作为划分点数目的上限。
5.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述步骤2)中,模糊逻辑规则的构建具体包括以下步骤:
2.1)对步骤1)中生成的用于高维模糊轨迹序列模型训练的集合中每个高维模糊轨迹序列,在序列后新增一维,值为该序列对应的用户轨迹序列中下一刻用户的2维平面位置坐标,生成用于聚类的高维模糊轨迹序列训练集;
2.2)对步骤2.1)中生成的用于聚类的高维模糊轨迹序列训练集,应用自适应高斯核的模糊C均值聚类算法生成高维轨迹序列簇,并计算簇中心作为基于一般样本的轨迹预测模型的模糊逻辑规则。
6.根据权利要求5所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述步骤2.2)中,自适应高斯核的模糊C均值聚类算法在完成聚类算法的每一次迭代后,依据现有样本相对于簇中心的分布情况,重新计算簇的核宽度并更新。
7.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述模糊逻辑规则定义为:假设x(t)是由F(t-k),F(t-k+1),…,F(t-1)决定的,且F(t)代表t时刻用户坐标对应的模糊集,x(t)为t时刻的用户坐标,那么就存在模糊逻辑规则F(t-k),F(t-k+1),...,F(t-1)→x(t)。
8.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:计算高维轨迹序列与所有簇间的隶属度,若所有隶属度均小于设定阈值η,则不满足模糊逻辑规则;若存在大于等于设定阈值η的隶属度,则满足模糊逻辑规则。
9.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述离群点预测模型的建立方法包括以下步骤:将离群点样本对应的用户轨迹序列转化为间隔时间的位移序列以及位移方向的序列,然后采用Levenberg-Marquardt算法分别构建关于下一个间隔时间的位移以及位移方向的预测模型;Levenberg-Marquardt算法需要通过多次迭代来避免陷入局部极小,每次迭代均在设定范围内随机生成初始权向量,并选择最优的结果;Levenberg-Marquardt算法的初始权向量设定时,每一维度的值大于前一维度。
10.根据权利要求1所述一种应用模糊轨迹序列的轨迹预测方法,其特征在于:所述步骤3)中,对簇中心最后S维求加权平均数时,若待预测样本与簇间隶属度大于等于阈值η,则权值为该隶属度;否则,权值为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710324975.3A CN107247961B (zh) | 2017-05-10 | 2017-05-10 | 一种应用模糊轨迹序列的轨迹预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710324975.3A CN107247961B (zh) | 2017-05-10 | 2017-05-10 | 一种应用模糊轨迹序列的轨迹预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107247961A true CN107247961A (zh) | 2017-10-13 |
CN107247961B CN107247961B (zh) | 2019-12-24 |
Family
ID=60016474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710324975.3A Active CN107247961B (zh) | 2017-05-10 | 2017-05-10 | 一种应用模糊轨迹序列的轨迹预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107247961B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268597A (zh) * | 2017-12-18 | 2018-07-10 | 中国电子科技集团公司第二十八研究所 | 一种动目标活动概率图谱构建及行为意图识别方法 |
CN108492355A (zh) * | 2018-03-05 | 2018-09-04 | 天津工业大学 | 基于线段传递的三维物体表面可变性度量方法 |
CN108509981A (zh) * | 2018-03-05 | 2018-09-07 | 天津工业大学 | 基于序列顶点特征的三维物体内部部件自动划分方法 |
CN108520044A (zh) * | 2018-04-03 | 2018-09-11 | 西安交通大学 | 一种基于公度性思想的最小二乘法路径预测算法 |
CN110163460A (zh) * | 2018-03-30 | 2019-08-23 | 腾讯科技(深圳)有限公司 | 一种确定应用分值的方法及设备 |
CN111260950A (zh) * | 2020-01-17 | 2020-06-09 | 清华大学 | 一种基于轨迹预测的轨迹跟踪方法、介质和车载设备 |
CN112019320A (zh) * | 2019-05-30 | 2020-12-01 | 中国科学技术大学 | 侧信道分析中的能量轨迹提取方法及系统 |
CN112380400A (zh) * | 2020-11-20 | 2021-02-19 | 浙江大学计算机创新技术研究院 | 一种基于后缀树误差有界的轨迹预测方法 |
CN112667763A (zh) * | 2020-12-29 | 2021-04-16 | 电子科技大学 | 一种基于自适应时间戳与多尺度特征提取的轨迹预测方法 |
CN113365298A (zh) * | 2021-04-22 | 2021-09-07 | 山东师范大学 | 基于高阶模糊认知图的mr网络信号强度预测方法及系统 |
CN113553516A (zh) * | 2021-09-18 | 2021-10-26 | 南京森根科技股份有限公司 | 一种基于模糊路径的频繁轨迹挖掘方法 |
CN114091630A (zh) * | 2022-01-24 | 2022-02-25 | 阿里巴巴达摩院(杭州)科技有限公司 | 获取目标对象移动轨迹信息的方法、装置及系统 |
CN114793460A (zh) * | 2019-12-10 | 2022-07-26 | 标致雪铁龙汽车股份有限公司 | 用于创建用于对交通参与者进行计算机模拟的算法的方法、用于训练至少一个用于机动车的控制设备的算法的方法、计算机程序产品及机动车 |
CN115297441A (zh) * | 2022-09-30 | 2022-11-04 | 上海世脉信息科技有限公司 | 一种大数据环境下个体时空活动鲁棒性的计算方法 |
CN115688682A (zh) * | 2022-12-29 | 2023-02-03 | 浙江工业大学 | 一种基于模糊预测的车辆轨迹数据压缩方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101256673A (zh) * | 2008-03-18 | 2008-09-03 | 中国计量学院 | 用于在实时视频跟踪系统中跟踪手臂运动的方法 |
US20090018918A1 (en) * | 2004-11-04 | 2009-01-15 | Manyworlds Inc. | Influence-based Social Network Advertising |
CN102567497A (zh) * | 2011-12-23 | 2012-07-11 | 浙江大学 | 一种最匹配模糊轨迹问题的查询方法 |
-
2017
- 2017-05-10 CN CN201710324975.3A patent/CN107247961B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090018918A1 (en) * | 2004-11-04 | 2009-01-15 | Manyworlds Inc. | Influence-based Social Network Advertising |
CN101256673A (zh) * | 2008-03-18 | 2008-09-03 | 中国计量学院 | 用于在实时视频跟踪系统中跟踪手臂运动的方法 |
CN102567497A (zh) * | 2011-12-23 | 2012-07-11 | 浙江大学 | 一种最匹配模糊轨迹问题的查询方法 |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108268597A (zh) * | 2017-12-18 | 2018-07-10 | 中国电子科技集团公司第二十八研究所 | 一种动目标活动概率图谱构建及行为意图识别方法 |
CN108492355A (zh) * | 2018-03-05 | 2018-09-04 | 天津工业大学 | 基于线段传递的三维物体表面可变性度量方法 |
CN108509981A (zh) * | 2018-03-05 | 2018-09-07 | 天津工业大学 | 基于序列顶点特征的三维物体内部部件自动划分方法 |
CN110163460A (zh) * | 2018-03-30 | 2019-08-23 | 腾讯科技(深圳)有限公司 | 一种确定应用分值的方法及设备 |
CN110163460B (zh) * | 2018-03-30 | 2023-09-19 | 腾讯科技(深圳)有限公司 | 一种确定应用分值的方法及设备 |
CN108520044A (zh) * | 2018-04-03 | 2018-09-11 | 西安交通大学 | 一种基于公度性思想的最小二乘法路径预测算法 |
CN112019320B (zh) * | 2019-05-30 | 2022-04-19 | 中国科学技术大学 | 侧信道分析中的能量轨迹提取方法及系统 |
CN112019320A (zh) * | 2019-05-30 | 2020-12-01 | 中国科学技术大学 | 侧信道分析中的能量轨迹提取方法及系统 |
CN114793460A (zh) * | 2019-12-10 | 2022-07-26 | 标致雪铁龙汽车股份有限公司 | 用于创建用于对交通参与者进行计算机模拟的算法的方法、用于训练至少一个用于机动车的控制设备的算法的方法、计算机程序产品及机动车 |
CN111260950A (zh) * | 2020-01-17 | 2020-06-09 | 清华大学 | 一种基于轨迹预测的轨迹跟踪方法、介质和车载设备 |
CN111260950B (zh) * | 2020-01-17 | 2021-03-26 | 清华大学 | 一种基于轨迹预测的轨迹跟踪方法、介质和车载设备 |
CN112380400A (zh) * | 2020-11-20 | 2021-02-19 | 浙江大学计算机创新技术研究院 | 一种基于后缀树误差有界的轨迹预测方法 |
CN112667763B (zh) * | 2020-12-29 | 2022-09-13 | 电子科技大学 | 一种基于自适应时间戳与多尺度特征提取的轨迹预测方法 |
CN112667763A (zh) * | 2020-12-29 | 2021-04-16 | 电子科技大学 | 一种基于自适应时间戳与多尺度特征提取的轨迹预测方法 |
CN113365298B (zh) * | 2021-04-22 | 2022-08-02 | 山东师范大学 | 基于高阶模糊认知图的mr网络信号强度预测方法及系统 |
CN113365298A (zh) * | 2021-04-22 | 2021-09-07 | 山东师范大学 | 基于高阶模糊认知图的mr网络信号强度预测方法及系统 |
CN113553516A (zh) * | 2021-09-18 | 2021-10-26 | 南京森根科技股份有限公司 | 一种基于模糊路径的频繁轨迹挖掘方法 |
CN114091630B (zh) * | 2022-01-24 | 2022-05-06 | 阿里巴巴达摩院(杭州)科技有限公司 | 获取目标对象移动轨迹信息的方法、装置及系统 |
CN114091630A (zh) * | 2022-01-24 | 2022-02-25 | 阿里巴巴达摩院(杭州)科技有限公司 | 获取目标对象移动轨迹信息的方法、装置及系统 |
CN115297441A (zh) * | 2022-09-30 | 2022-11-04 | 上海世脉信息科技有限公司 | 一种大数据环境下个体时空活动鲁棒性的计算方法 |
CN115297441B (zh) * | 2022-09-30 | 2023-01-17 | 上海世脉信息科技有限公司 | 一种大数据环境下个体时空活动鲁棒性的计算方法 |
CN115688682A (zh) * | 2022-12-29 | 2023-02-03 | 浙江工业大学 | 一种基于模糊预测的车辆轨迹数据压缩方法及装置 |
CN115688682B (zh) * | 2022-12-29 | 2023-11-14 | 浙江工业大学 | 一种基于模糊预测的车辆轨迹数据压缩方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN107247961B (zh) | 2019-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107247961A (zh) | 一种应用模糊轨迹序列的轨迹预测方法 | |
CN110928993B (zh) | 基于深度循环神经网络的用户位置预测方法及系统 | |
Li et al. | Block-wisely supervised neural architecture search with knowledge distillation | |
CN107396322A (zh) | 基于路径匹配与编码译码循环神经网络的室内定位方法 | |
CN109409499A (zh) | 一种基于深度学习和卡尔曼滤波修正的轨迹恢复方法 | |
CN112969144B (zh) | 一种面向移动边缘计算的微服务预部署方法及系统 | |
CN102901500A (zh) | 基于概率a星与智能体混合的飞行器最优路径确定方法 | |
Xu et al. | HighAir: A hierarchical graph neural network-based air quality forecasting method | |
CN107293115A (zh) | 一种用于微观仿真的交通流量预测方法 | |
CN110334589A (zh) | 一种基于空洞卷积的高时序3d神经网络的动作识别方法 | |
CN101196905A (zh) | 一种智能图形检索方法 | |
Tastan et al. | Learning to intercept opponents in first person shooter games | |
CN106022471A (zh) | 基于粒子群优化算法的小波神经网络模型船舶横摇实时预测方法 | |
CN112699596B (zh) | 一种基于学习的广域电磁法激电信息非线性提取方法 | |
CN110443285A (zh) | 相似轨迹的确定方法、装置及计算机存储介质 | |
CN114047770A (zh) | 一种多内心搜寻改进灰狼算法的移动机器人路径规划方法 | |
CN110149595A (zh) | 一种基于hmm的异构网络用户行为预测方法 | |
CN109961129A (zh) | 一种基于改进粒子群的海上静止目标搜寻方案生成方法 | |
CN110956684A (zh) | 基于残差网络的人群运动疏散仿真方法及系统 | |
CN114461931A (zh) | 一种基于多关系融合分析的用户轨迹预测方法和系统 | |
CN107104747A (zh) | 无线时变信道中的多径分量的分簇方法 | |
CN107562837A (zh) | 一种基于道路网的机动目标跟踪算法 | |
CN113139644B (zh) | 一种基于深度蒙特卡洛树搜索的信源导航方法及装置 | |
Wang et al. | Indoor sound source localization algorithm based on BP neural network | |
CN114646313A (zh) | 一种用户轨迹定位方法、电子设备及计算机存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |