[发明专利]一种区域空气PM2.5浓度预测方法

摘要

本发明公开了一种区域空气PM2.5浓度预测方法，该方法首先通过历史数据来构造待训练的支持向量机回归模型的训练样本数据，然后通过训练样本数据得到训练的支持向量机回归模型，将该训练后的支持向量机回归模型作为PM2.5浓度预测模型；再将粒子群寻优算法与PM2.5浓度预测模型相结合，通过粒子群寻优算法不断寻优迭代，从而不断采用粒子的位置去重构PM2.5浓度预测模型的输入参数，直至迭代完成后得到粒子群最终的全局极性，采用粒子群最终的全局极值对应的粒子的位置去重构PM2.5浓度预测模型的输入参数，将该输入参数输入PM2.5浓度预测模型中得到的输出即为PM2.5浓度；优点是可降低PM2.5浓度预测模型输入参数的维数，提高PM2.5浓度预测准确率。

主权项

一种区域空气PM2.5浓度预测方法，其特征在于包括以下步骤：①构造待训练的支持向量机回归模型的训练样本数据：①‑1将当前时刻记为t时刻，t‑n时刻表示当前时刻的前n个时刻，n＝1,2,3,4,…；将t‑1时刻测量得到的区域空气的大气能见度记为xN(t‑1)，浊度记为xZ(t‑1)，气温记为xT(t‑1)，气压记为xp(t‑1)，风速记为xw(t‑1)；将t‑td时刻测量得到的区域空气中SO2的浓度记为xso2(t‑td)，NO2的浓度记为xNo2(t‑td)，PM10的浓度记为xpm10(t‑td)，PM2.5的浓度记为xpm25(t‑td)，O3的浓度记为xo3(t‑td)，CO的浓度记为xco(t‑td)；td为时滞因子且td＝1,2,…,d，d为大于等于2的整数；将t时刻预测的PM2.5浓度记为ypm25(t)；①‑2将待训练的支持向量机模型f(X)＝<ω,X>+b的输入参数记为Xin(T)，待训练的支持向量机回归模型的输出记为Yout(T)，Xin(T)和Yout(T)构成训练参数对，令Xin(T)=[xso2(T),xso2(T-1),...,xso2(T-d),xso2(T),xso2(T-1),...,xso2(T-d),xpm10(T),xpm10(T-1),...,xpm10(T-d),xpm25(T),xpm25(T-1),...,xpm25(T-d),xO3(T),xO3(T-1),...,xO3(T-d),xCO(T),xCO(T-1),...,xCO(T-d),xN(t-1),xZ(t-1),xT(t-1),xP(t-1),xw(t-1)]---(1)]]>Yout(T)＝xpm25(T+1)                                                     (2)Xin(t-1)=[xso2(t-1),xso2(t-2),...,xso2(t-1-d),xno2(t-1),xno2(t-2),...,xno2(t-1-d),xpm10(t-1),xpm10(t-2),...,xpm10(t-1-d),xpm25(t-1),xpm25(t-2),...,xpm25(t-1-d),xO3(t-1),xO3(t-2),...,xO3(t-1-d),xCO(t-1),xCO(t-2),...,xCO(t-1-d),xN(t-1),xZ(t-1),xT(t-1),xP(t-1),xw(t-1)]---(3)]]>其中T＝t‑2，t‑3，t‑4，…t‑501，[]为矩阵符号，Xin(T)为1×(6×d+5)的矩阵，符号×为乘运算符号；①‑3获取t时刻之前区域空气中SO2浓度、NO2浓度、PM10浓度、PM2.5浓度、O3浓度、CO浓度、大气能见度、浊度、气温记、气压和风速的历史测量数据，得到{Xin(t‑2)，Yout(t‑2)}，{Xin(t‑3)，Yout(t‑3)}，…，{Xin(t‑501)，Yout(t‑501)}共500组训练参数对，将该500组训练参数对作为待训练的支持向量机回归模型的训练样本数据；②采用训练样本数据对待训练的支持向量机回归模型f(X)＝<ω,X>+b进行训练，得到支持向量机回归模型的模型参数ω和b；将得到的支持向量机回归模型的模型参数ω和b代入支持向量机回归模型f(X)＝<ω,X>+b中得到训练后的支持向量机回归模型，将该训练后的支持向量机回归模型作为PM2.5浓度预测模型；③初始化粒子群参数：在[0,1]内随机产生种群规模为N的粒子群，N＝100，每个粒子均含有位置和速度属性，随机初始化粒子群中各粒子的速度和位置，将第i个粒子初始化后的当前位置记为Wi(0)，当前速度记为Vi(0)，i＝1,2,…,N，其中Wi(0)和Vi(0)均为1×m的矩阵，Wi(0)＝[w1i(0),w2i(0),…,wmi(0)],i＝1,2,…,N，w1i(0),w2i(0),…,wmi(0)∈[0,1]，Vi(0)＝[v1i(0),v2i(0),…,vmi(0)],i＝1,2,…,N，v1i(0),v2i(0),…,vmi(0)∈[0,1]，m＝6×d+5；将第i个粒子的在第k次迭代后的当前位置记为Wi(k)，Wi(k)＝[w1i(k),w2i(k),…,wmi(k)]，第k次迭代后的当前速度记为Vi(k)，Vi(k)＝[v1i(k),v2i(k),…,vmi(k)]，k＝1，2，…，kmax；kmax为粒子群的总迭代次数；④将粒子群中每个粒子的当前位置定义为PM2.5浓度预测模型输入参数的权系数矩阵；采用N个粒子的当前位置分别对PM2.5浓度预测模型的输入参数Xin(t‑1)进行重构，由此得到初始状态下N个重构后的PM2.5浓度预测模型的输入参数，将初始状态下第i个重构后的PM2.5浓度预测模型的输入参数记为Xin0(t-1,i)=[xso2(t-1)×w1i(0),xso2(t-2)×w2i(0),...,xso2(t-d)×wdi(0),]]>xno2(t‑1)×w(d+1)i(0),xno2(t‑2)×w(d+2)i(0),…,xno2(t‑d)×w(2d)i(0),xpm10(t‑1)×w(2d+1)i(0),xpm10(t‑2)×w(2d+2)i(0),…,xpm10(t‑d)×w(3d)i(0),xpm25(t‑1)×w(3d+1)i(0),xpm25(t‑2)×w(3d+2)i(0),…,xpm25(t‑d)×w(4d)i(0),xO3(t‑1)×w(4d+1)i(0),xO3(t‑2)×w(4d+2)i(0),…,xO3(t‑d)×w(5d)i(0),xCO(t‑1)×w(5d+1)i(0),xCO(t‑2)×w(5d+2)i(0),…,xCO(t‑d)×w(6d)i(0),xN(t‑1)w(6d+1)i(0),xZ(t‑1)w(6d+2)i(0),xT(t‑1)w(6d+3)i(0),xP(t‑1)w(6d+4)i(0),xw(t‑1)w(6d+5)i(0)]⑤将粒子群中第i个粒子的个体极值的初始值记为Fi(0)＝|xpm25(t‑1)‑youti(t,0)|，其中| |为取绝对值符号，youti(t,0)为作为输入参数输入PM2.5浓度预测模型中得到PM2.5浓度预测模型的输出；将粒子群的全局极值的初始值记为⑥采用粒子群每次迭代后更新后的位置对重构后的PM2.5浓度预测模型的输入参数进行更新，粒子群第k次迭代后重构后的PM2.5浓度预测模型的输入参数进行第k次更新，将第k次更新后的第i个重构后的PM2.5浓度预测模型的输入参数记为Xink(t-1,i)=[xso2(t-1)×w1i(k),xso2(t-2)×w2.1(k),...,xso2(t-d)×wdi(k),]]>xno2(t‑1)×w(d+1)i(k),xno2(t‑2)×w(d+2)i(k),…,xno2(t‑d)×w(2d)i(k),xpm10(t‑1)×w(2d+1)i(k),xpm10(t‑2)×w(2d+2)i(k),…,xpm10(t‑d)×w(3d)i(k),xpm25(t‑1)×w(3d+1)i(k),xpm25(t‑2)×w(3d+2)i(k),…,xpm25(t‑d)×w(4d)i(k),xO3(t‑1)×w(4d+1)i(k),xO3(t‑2)×w(4d+2)i(k),…,xO3(t‑d)×w(5d)i(k),xCO(t‑1)×w(5d+1)i(k),xCO(t‑2)×w(5d+2)i(k),…,xCO(t‑d)×w(6d)i(k),xN(t‑1)w(6d+1)i(k),xZ(t‑1)w(6d+2)i(k),xT(t‑1)w(6d+3)i(k),xP(t‑1)w(6d+4)i(k),xw(t‑1)w(6d+5)i(k)]将作为输入参数输入PM2.5浓度预测模型中，得到PM2.5浓度预测模型的输出，记为将代入公式Fi(k)＝|xpm25(t‑1)‑youti(t,k)|中，计算得到第k次更新后第i个粒子对应的适应度值Fi(k)；| |为取绝对值符号；第k次迭代后粒子群中第i个粒子的个体极值记为将第k次迭代后第i个粒子的个体极值对应的粒子位置记为Wi(a)，a为大于等于0且小于等于k的整数；将第k次迭代后粒子群的全局极值记为min表示取最小值；该全局极值对应的粒子的位置记为Wh(g)，h为大于等于1且小于等于N的整数，g为大于等于0且小于等于k的整数；⑦从k＝1开始对粒子群进行更新迭代，得到更新后的各粒子的速度和位置，按照步骤⑥更新粒子的个体极值和粒子群的全局极值，直至k＝kmax，结束更新迭代，得到最终的粒子群的全局极值；⑧采用最终的粒子群的全局极值对应的粒子的位置对重构后的PM2.5浓度预测模型的输入参数进行更新，将更新后的输入参数输入PM2.5浓度预测模型中，此时PM2.5浓度预测模型的输出即为预测到的PM2.5浓度。