[发明专利]一种基于多源数据融合的公共自行车租借预测方法

权利要求书

1.一种基于多源数据融合的公共自行车租借预测方法，其特征在于包括如下步骤：

步骤(1)数据采集与处理：

读取公共自行车系统中存储的用户租车记录，计算每个站点每隔时间段time剩余公共自行车数量，针对每个站点建立向量N_date＝[n₁,n₂,...,n_j,...,n_M]；从互联网上获取历史天气状况、气温和节假日情况，建立环境属性向量S＝{S₁,S₂,...,S_i,...,S_L}，S_i＝[d,w,t],L为数据样本集中样本的个数；M、L、i和j均为正整数，其中j≤M,i≤L；

所述的用户租车记录包括租借时间、租借站点、租借车桩号、还车时间、还车站点、还车车桩号；

n_j表示第j个时间段time该站点内剩余自行车数量；

S_i表示第i天的环境属性；

d表示当天的节假日属性，d的取值为

w表示天气状况，w的取值为

t表示气温，t的取值为

步骤(2)K-means聚类：

N_date＝[n₁,n₂,...,n_j,...,n_M]作为待分类的数据样本集中的一个数据元组；从数据样本集中随机选取k′个点作为初始类中心，分别用μ₁,μ₂,...,μ_k′表示类c₁,c₂,...,c_k′的类中心；其中，k′为正整数；

K-means聚类算法包括如下步骤：

2-1.计算其余的数据元组N_i′与类中心μ_k′的距离，将其余的数据元组N_i′分到与其距离最短的类中；

2-2.重新计算各类的新的类中心μ_k′^*；

2-3.重复步骤2-1和步骤2-2，直到各数据元组所属类别都不再变化；

步骤(3)构建贝叶斯分类器：使用k-means聚类算法对数据样本集进行分类后，得到k个类，数据样本集中的每个数据元组都将得到一个类标号；贝叶斯分类法基于贝叶斯定理：

P(H|X)=P(X|H)P(H)P(X)]]>

计算P(X)，P(H)和P(X|H)；

步骤(4)自适应粒子群神经网络模型训练

使用分类后的数据样本训练各类别所属的自适应粒子群神经网络模型，设定神经网络的输入层单元个数为n，隐含层单元个数为p，输出层单元个数为q；则神经网络中需要调整的权值个数N＝n×p+p×q，用需要调整的权值个数N表示问题空间的维度；用向量X_i＝[x_i1,x_i2,...,x_iN]表示神经网络所有权值的一个可能取值，也表示粒子群优化算法中粒子的位置；粒子群神经网络模型训练步骤如下：

4-1.初始化粒子群种群中各粒子的位置和速度，粒子的位置和速度位于[0,1]之间，即随机生成G个向量X_i＝[x_i1,x_i2,...,x_iN],1≤i≤G作为粒子的位置，x_iN介于[0,1]之间；再随机生成G个向量V_i＝[v_i1,v_i2,...,v_iN],G+1≤i≤2×G作为粒子的速度，v_iN介于[0,1]之间；设定最大迭代次数itermax，当前迭代次数k＝1；

4-2.计算各粒子的适应度值，存储各粒子当前的位置为自身最优值P_bi，存储所有粒子中适应度值最小的位置为全局最优值P_bg；

4-3.判断k是否达到最大迭代次数itermax，若是，转到步骤4-7，否则转到4-4；

4-4.更新各粒子速度和位置，计算各粒子新的适应度值；若第i个粒子的当前适应度值优于P_bi，则更新P_bi；若当前粒子群的最优位置优于P_bg，则更新P_bg；设置k＝k+1；

4-5.根据惯性权重计算公式更新惯性权重；

4-6.若全局最优位置P_bg在10次迭代中没有变化，则转到步骤4-7，否则转到步骤4-4；

4-7.输出P_bg作为训练好的粒子群神经网络的权值；

步骤(4)中使用的计算公式及变量含义如下：

适应度值fitness(X_g)的计算公式为：

fitness(Xg)=12M′Σl=1M′Σk=1q((Σj=1pwjk(1/(1+exp(Σi=1nwijxli-θhj)))-θok)-dlk)]]>

M′为数据样本集中数据元组个数；

x_li为第l个数据样本的第i个输入；

θ_hj为隐含层第j个单元的阈值；

θ_ok为输出层第k个单元的阈值；

w_ij为第i个输入层与第j个隐含层的连接权值；

w_jk为第j个隐含层与第k个输出层的连接权值；

d_lk为第l个数据样本的第k个输出；

速度计算公式为：

Vik+1=w×Vik+c1×r1×(Pik-Xik)+c2×r2×(Pgk-Xik)]]>

V_i^k为粒子X_i在第k次迭代的速度；

r₁,r₂是取值在[0,1]之间的随机数；

c₁，c₂为一组随机数，取c₁＝c₂＝2.5；

ω为惯性权重；

惯性权重ω_k计算公式为：

ωk=ωmax-ωmax-ωminitmax×k]]> 2

ωmax＝0.9，ωmin＝0.2；

粒子的位置计算公式为：

Xik+1=Xik+Vik+1;]]>

X_i^k为第k次迭代时第i个粒子的位置；

步骤(5)租借模式预测：根据输入的待预测环境属性向量S_pre和计算得到的P(H)，P(X|H)的值，计算出该数据元组属于各类别的概率，输出后验概率最大的类别；

步骤(6)公共自行车租借预测

根据步骤(5)租借模式预测中输出的类别，选择相应的自适应粒子群神经网络模型，将前n′个时间段的剩余公共自行车数量输入该模型，输出即为预测的公共自行车数量。