CN106127353A - 一种最优路线搜索*** - Google Patents

Abstract

本发明一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块。本发明预测精度较高且构造的预测模型更有针对性。

Description

一种最优路线搜索***

技术领域

本发明涉及智能交通领域，具体涉及一种最优路线搜索***。

背景技术

交通流量是指单位时间内通过道路某一断面的实际车辆数，是描述交通状态的重要特征参数。交通流量的变化又是一个实时、高维、非线性、非平稳的随机过程，相关因素的变化都可能影响下一时刻的交通流量。相关技术中，关于短时的预测装置局限性强，预测精度较低，实时预测未能取得令人满意的结果，未能对人们的实时道路选择提供有效建议，从而交通流量预测大部分停留在交通流量的中长期预测。

发明内容

针对上述问题，本发明提供一种最优路线搜索***。

本发明的目的采用以下技术方案来实现：

一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本发明的有益效果为：

1、设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；

2、设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，其中预测因子直接影响预测精度，相关系数是测量随机变量相关性的指标，能够帮助选取与预测点密切相关的变量作为预测模型的训练样本，选取多个相关系数作为预测因子，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；

3、相关系数计算模块中的空间相关系数反映了路网的可达性对预测模型的影响，时间相关系数能够表达流量序列的时间顺序，反映两序列时间上的因果关系，从而提高预测因子选取的效率；由于交通流量的周相似性，引入历史相关系数矩阵生成模块的历史相关系数，同时间相关系数和空间相关系数配合使用，为准确预测提供更多的数据支持。

附图说明

利用附图对本发明作进一步说明，但附图中的实施例不构成对本发明的任何限制，对于本领域的普通技术人员，在不付出创造性劳动的前提下，还可以根据以下附图获得其它的附图。

图1是本发明预测装置各模块的连接示意图。

图2是本发明搜索***的实现方法的流程图。

具体实施方式

结合以下实施例对本发明作进一步描述。

实施例1

参见图1，图2，本实施例一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本实施例设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；本实施例取值L＝8，M＝3，预测精度相对于相关技术提高了1.5％。

实施例2

参见图1，图2，本实施例一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本实施例设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；本实施例取值L＝9，M＝3，预测精度相对于相关技术提高了2％。

实施例3

参见图1，图2，本实施例一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本实施例设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；本实施例取值L＝10，M＝4，预测精度相对于相关技术提高了2.6％。

实施例4

参见图1，图2，本实施例一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本实施例设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；本实施例取值L＝11，M＝5，预测精度相对于相关技术提高了3.2％。

实施例5

参见图1，图2，本实施例一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

优选地，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

优选地，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

优选地，预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

$P\left(\mathrm{\τ}\right)=\frac{E\[(X_x-{\mathrm{\ν}}_i)(X_{x+\mathrm{\τ}}-{\mathrm{\ν}}_{x+\mathrm{\τ}})\]}{{\mathrm{\σ}}^2}$

其中，X_x表示待检验交通流量序列，ν_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，ν_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

${\mathrm{\ρ}}_{ij}\left(\mathrm{\τ}\right)=\frac{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right){\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})}{{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_i\left(t\right)-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_i\left(t\right)\]}}^2\×{\sqrt{{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}\[x_j(t+\mathrm{\τ})-\frac{1}{n-\mathrm{\τ}}{\mathrm{\Σ}}_{t=1}^{n-\mathrm{\τ}}{x}_j(t+\mathrm{\τ})\]}}^2}$

空间相关系数ρ_ij(w)的计算公式为：

优选地，预测装置还包括：

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

$\mathrm{\ρ}{\left(\mathrm{\τ}\right)}^{\′}=\left[\begin{array}{cccc}{\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& {\mathrm{\ρ}}_{1j}{\left(0\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(0\right)}^{\′}\\ {\mathrm{\ρ}}_{1j}{\left(1\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(1\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(1\right)}^{\′}\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {\mathrm{\ρ}}_{1j}{\left(L\right)}^{\′}& {\mathrm{\ρ}}_{2j}{\left(L\right)}^{\′}& \mathrm{...}& {\mathrm{\ρ}}_{Nj}{\left(L\right)}^{\′}\end{array}\right];$

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

${\mathrm{\ρ}}_{jm}\left(t\right)=\frac{{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right){\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)}{\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_j\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_j\left(t\right)\]}^2}\×\sqrt{{\mathrm{\Σ}}_{t=1}^n{\[x_{jm}\left(t\right)-\frac{1}{n}{\mathrm{\Σ}}_{t=1}^n{x}_{jm}\left(t\right)\]}^2}}$

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]|and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

$X^{\′}=\left[\begin{array}{cccc}{x_1}^{\′}\left(1\right)& {x_2}^{\′}\left(1\right)& \mathrm{...}& {x_p}^{\′}\left(1\right)\\ {x_1}^{\′}\left(2\right)& {x_2}^{\′}\left(2\right)& \mathrm{...}& {x_p}^{\′}\left(2\right)\\ \mathrm{...}& \mathrm{...}& \mathrm{...}& \mathrm{...}\\ {x_1}^{\′}(n-L_1)& {x_2}^{\′}(n-L_1)& \mathrm{...}& {x_p}^{\′}(n-L_1)\end{array}\right];$

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，

Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

其中，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

其中，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段的交通流量序列与预测路段的交通流量序列分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。

本实施例设置数据分类模块和平稳性检验模块，增加了数据的准确度，且使构造的预测模型更有针对性；设置相关系数计算模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块，消除了最初预测因子选取的主观性，能够增加预测精度，使预测模型构造模块更加稳定和准确；本实施例取值L＝12，M＝5，预测精度相对于相关技术提高了3.5％。

最后应当说明的是，以上实施例仅用以说明本发明的技术方案，而非对本发明保护范围的限制，尽管参照较佳实施例对本发明作了详细地说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和范围。

Claims (7)

1.一种最优路线搜索***，包括搜索***和与搜索***相连的预测装置，所述搜索***的实现方法包括：

a.接收用户的查询请求；

b.从常规索引库搜索与所述查询请求相关联的常规搜索结果，并从个性化索引库搜索与所述用户的识别信息以及所述查询请求相关联的个性化搜索结果；

c.对所述常规搜索结果以及所述个性化搜索结果进行合并，并发送合并后的所述常规搜索结果以及所述个性化搜索结果。

2.根据权利要求1所述的一种最优路线搜索***，其特征是，所述步骤b进一步包括：通过用户登录来获取所述识别信息。

3.根据权利要求2所述的一种最优路线搜索***，其特征是，所述步骤c进一步包括对所述常规搜索结果进行排序，以及对所述个性化搜索结果进行排序。

4.根据权利要求3所述的一种最优路线搜索***，其特征是，所述预测装置包括依次连接的采集模块、数据预处理模块、数据分类模块、平稳性检验模块、相关系数计算模块、阈值设定模块、时空相关系数矩阵生成模块、历史相关系数矩阵生成模块、预测因子选取模块和预测模型构造模块：

(1)采集模块，用于采集路网S内观测路段S_i、预测路段S_j对应各时间段的交通流量数据和通行情况；

(2)数据预处理模块，用于对所述交通流量数据进行数据预处理，并剔除不符合交通实际情况的数据；

(3)数据分类模块，用于对经过数据预处理的交通流量数据进行类型分类，所述类型包括节假日交通流量数据、周末交通流量数据和工作日交通流量数据；

(4)平稳性检验模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验，检验平稳性的自相关函数为：

其中，X_x表示待检验交通流量序列，v_i表示待检验交通流量序列的均值，X_x+τ表示X_x在时间延迟τ后的交通流量序列，v_x+τ为X_x+τ的均值，σ²为X_x与X_x+τ之间的方差；

当自相关函数P(τ)能快速衰减趋近于0或在0附近波动，则所述待检验交通流量序列通过平稳性检验；当自相关函数P(τ)不能快速衰减趋近于0或在0附近波动，则对所述待检验交通流量序列进行平稳处理后继续进行平稳性检验；

(5)相关系数计算模块，用于计算通过平稳性检验的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j在时间延迟τ下的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)，设路网S内有N个路段，交通流量序列X_i＝[x_i(1),x_i(2),...,x_i(n)]，交通流量序列x_i(t)表示观测路段S_i在t时刻的流量，x_j(t)表示预测路段S_j在t时刻的流量，t＝1,2,...n，时间相关系数ρ_ij(τ)的计算公式为：

空间相关系数ρ_ij(w)的计算公式为：

。

5.根据权利要求4所述的一种最优路线搜索***，其特征是，

(6)阈值设定模块，用于设定各路段之间的时间延迟最大值L、时空相关系数阈值T₁和历史相关系数阈值T₂；

(7)时空相关系数矩阵生成模块，用于根据各路段的时间相关系数ρ_ij(τ)和空间相关系数ρ_ij(w)构建各观测路段S_i与预测路段S_j在不同时间延迟τ下的时空相关系数矩阵ρ(τ)'，并计算各路段的时空相关系数ρ_ij(τ)'，其中i∈[1,N]且τ∈[0,L]，L的取值范围为[8,12]，时空相关系数矩阵ρ(τ)'的计算公式为：

时空相关系数ρ_ij(τ)'的计算公式为：

ρ_ij(τ)'＝ρ_ij(τ)ρ_ij(w)；

(8)历史相关系数矩阵生成模块，用于生成预测路段S_j的历史相关系数矩阵ρ(t)：

其中，选取近M周的同期且同一类型的历史流量作为交通流量序列X_j的历史相关序列，记为m＝1,2,...M，M的取值范围为[3,5]，所述历史相关系数ρ_jm(t)的计算公式为：

(9)预测因子选取模块，用于根据所述时空相关系数阈值T₁和历史相关系数阈值T₂选取与预测目标点相关的预测因子，并按照其所选空间位置j与时间延迟τ进行矩阵重构，选取原则为：

若ρ_ij(τ)'＞T₁，则将观测路段S_i的交通流量序列X_i中满足条件的交通流量组成新的序列并作为第一预测因子，记做X',X'＝(x₁',x₂',...,x_p')，其中p为所述满足条件的交通流量个数，设L₁为第一预测因子中时间延迟的最大值，L₁＝max{τ|τ∈[0,L]and ρ_ij(τ)'＞T₁}，则第一预测因子X'可表述成如下矩阵形式：

若ρ_jm(t)＞T₂，则将所有满足条件的历史相关序列X_jm(t)作为第二预测因子，记作Y'，Y'＝{y₁',y₂',...,y_q'}，其中q为满足条件的历史流量个数，第二预测因子Y'可表述成如下矩阵形式：

(10)预测模型构造模块，其通过将第一预测因子和第二预测因子作为训练样本来构造可预测路段在下一时刻的交通流量的预测模型。

6.根据权利要求5所述的一种最优路线搜索***，其特征是，所述数据预处理模块中，剔除所述不符合交通实际情况的数据的规则为：在一个数据更新周期内，分别设定各路段的总交通流量数据的阀值范围，若采集到的某路段的总交通流量数据落在对应的阈值范围内，则表明该组数据可靠，保留该组数据；若采集到的某路段的总交通流量数据落不在对应的阈值范围内，则表明该组数据不可靠，并将其剔除。

7.根据权利要求6所述的一种最优路线搜索***，其特征是，所述平稳性检验模块包括以下子模块：

(1)检验子模块，用于对处于同一类型的观测路段S_i的交通流量序列X_i与预测路段S_j的交通流量序列X_j分别进行平稳性检验；

(2)连续性检查子模块，与检验子模块连接，用于对不通过平稳性检验的待检验交通流量序列进行连续性检查，若不符合连续性，所述连续性检查子模块采用平均插值法对数据进行补齐；

(3)排错子模块，与连续性检查子模块连接，用于删除明显错误的数据，同时采用平均插值法对数据进行补齐；

(4)差分处理子模块，连接排错子模块和检验子模块，用于对补齐后的数据进行差分处理，并将差分处理后的数据传送到检验子模块。