一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法技术

本发明专利技术提出一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，该方法为获取城市历史气象数据、历史臭氧数据及预测的未来气象数据，使用历史气象数据充当预报因子，臭氧浓度充当预报量，建立多元线性回归模型，并用最小二乘法估计回归系数，得到确定的多元线性回归方程并进行显著性检验，然后将预测的未来气象数据输入多元线性回归模型中，预测出未来城市近地面臭氧浓度，本发明专利技术通过对残差的控制，回归方程的显著性检验，保证了对城市未来近地面臭氧浓度预报的可靠性。近地面臭氧浓度预报的可靠性。近地面臭氧浓度预报的可靠性。

【技术实现步骤摘要】
一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法


[0001]本专利技术属于大气污染物浓度预测
，具体涉及一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法。

技术介绍

[0002]目前，全球臭氧浓度正处于持续上涨趋势，从我国的情况来看，近两年来，空气质量总体上有所改善，各种空气污染物的浓度降低了，但臭氧浓度不降反增。臭氧已成为影响AQI优良率的重要因素之一。
[0003]如中国专利公开了“一种近地面臭氧浓度估算方法”(专利号：CN202010090277.3)，该专利利用卫星遥感数据估算近地面臭氧浓度，通过将具有时间序列特点的长短期记忆网络模型和体现空间位置特征的胶囊模型结合，在时间和空间上均反映近地面臭氧浓度，提高了近地面臭氧估算精度，可为相关环保部门提供支撑服务。
[0004]但是，该专利无法解决不同地区因气象条件不同造成的臭氧浓度影响因素不同和应用于预测未来城市近地面层臭氧浓度的问题。

技术实现思路

[0005](一)解决的技术问题
[0006]针对现有技术的不足，本专利技术提供了一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，解决了不同地区因气象条件不同造成的臭氧浓度影响因素不同的问题，形成了一种简单、高效的的算法，能够广泛应用于不同地区。
[0007](二)技术方案
[0008]为实现以上目的，本专利技术通过以下技术方案予以实现：一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，包括以下步骤：
[0009]步骤1：获取城市一段时间序列内的历史气象数据、历史臭氧数据及预测的未来气象数据；
[0010]步骤2：对获取的城市一段时间序列内的历史气象数据、历史臭氧数据及预测的未来气象数据进行预处理，得到城市一段时间序列内的历史逐日气象数据、逐日最大臭氧八小时滑动平均值数据及预测的未来逐日气象数据，历史逐日气象数据充当预报因子，逐日最大臭氧八小时滑动平均值充当预报量；
[0011]步骤3：建立气象数据与逐日最大臭氧八小时滑动平均值的多元线性回归方程模型；
[0012]步骤4：用最小二乘法估计线性回归方程的回归系数，既残差平方和最小来估计多元线性回归方程模型中各气象预报因子的回归系数，得到确认的多元线性回归方程；
[0013]步骤5：根据确认的多元线性回归方程进行历史回代，进行回归效果分析，并计算复相关系数；
[0014]步骤6：对多元线性回归方程的回归效果进行显著性检验；
[0015]步骤7：预报值的置信区间；
[0016]步骤8：将城市预测的未来逐日气象数据输入确定的多元线性回归方程中，计算出城市未来的近地面逐日最大臭氧八小时滑动平均值。
[0017]优选的，所述步骤3中的建立的多元线性回归方程公式为：
[0018][0019]其中为预报量的估计值，x1，x2，
…
，x
p
为各气象预报因子，p为气象因子的个数，b0为回归常数，b1，b2，
…
，b
p
为各气象预报因子的回归系数。
[0020]优选的，所述步骤4中的用最小二乘法估计线性回归方程的回归系数，既残差平方和最小来估计多元线性回归方程模型中各气象预报因子的回归系数，其残差平方和公式为：
[0021][0022]其中y为预报量的实际观测值，为预报量的估计值，Q为预报量的实际观测值和估计值的残差平方和，实际观测值y既逐日最大臭氧八小时滑动平均值，i＝1，2，
……
，n，其中n为气象预报因子的时间序列长度，既历史气象数据的天数。
[0023]使预报量的实际观测值y与预报量的估计值的残差平方和Q达最小值，从而确定各气象预报因子的多元线性回归方程的回归系数b1,b2,
……
，b
p
和回归常数b0。
[0024]优选的，所述步骤5的根据确认的多元线性回归方程进行历史回代包括以下步骤：
[0025]A1:将城市历史逐日气象数据代入确定的多元线性回归方程中
[0026]A2:计算得出城市历史逐日最大臭氧八小时滑动平均值的估计值i＝1,2,
…
，n优选的，所述步骤5的回归效果分析，其公式为：
[0027]L
yy
＝Q+U
[0028]其中：
[0029]总离差平方和，表示样本本身的变化
[0030]回归平方和，表示各自变量与y的线性关系而引起的变化部分
[0031]其中：
[0032]残差平方和，其表示除去各自变量与y的线性关系以外的其他所有因素对y的贡献。
[0033]优选的，所述步骤5的复相关系数为衡量一个变量与多个变量之间的线性关系程度的量，即预报量与各预报因子之间的线性关系程度的量，因为变量之间的关系可归结为一个多元线性回归方程，所以复相关系数是衡量预报量y与估计量之间线性相关程度的
量，记为R，其公式为：
[0034][0035]优选的，所述步骤6回归效果的显著性检验包括以下步骤：
[0036]B1：假设变量x1，x2，
…
，x
p
预报量y之间无线性关系，则多元线性回归方程中的回归系数b1，b2，
…
，b
p
均应为0，即归结要检验假设：
[0037]H0:b1＝b2＝
…
＝b
p
＝0
[0038]是否成立
[0039]B2:进行F检验：
[0040][0041]在显著性水平
ɑ
＝0.05下，若F>F
ɑ
则否定假设，即认为回归方程是显著的，否则，接受原假设，认为回归方程不显著，则更换气象预报因子，直到回归方程显著。
[0042]优选的，所述步骤7预报值的置信区间为：
[0043][0044][0045]优选的，所述步骤8为将城市预测的未来逐日气象数据输入确定的多元线性回归方程中，计算出城市未来的近地面逐日最大臭氧八小时滑动平均值。
[0046](三)有益效果
[0047]本专利技术提供了一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法。具备以下有益效果：
[0048](1)该一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，本专利技术所采用的多元线性回归模型方法，通过控制残差，对回归方程进行显著性检验，保证了臭氧预测结果的精确性；可以选择与研究对象相关的气象因子作为预报因子参与多元回归，得到与研究对象有关的多元线性回归方程，从而为研究预报量与预报因子之间的关系提供依据；有效的建立了臭氧预测多元线性回归方程，能够预测未来城市近地面层臭氧浓度，并形成了固定的算法，能被应用于预测未来城市近地面层臭氧浓度，以便提前采取有效的大气污染防治相关措施。
附图说明
[0049]图1为本专利技术具体实施方式中基于多元线性回归模型的城市近地面本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，其特征在于，包括以下步骤：步骤1：获取城市一段时间序列内的历史气象数据、历史臭氧数据及预测的未来气象数据；步骤2：对获取的城市一段时间序列内的历史气象数据、历史臭氧数据及预测的未来气象数据进行预处理，得到城市一段时间序列内的历史逐日气象数据、逐日最大臭氧八小时滑动平均值数据及预测的未来逐日气象数据，历史逐日气象数据充当预报因子，逐日最大臭氧八小时滑动平均值充当预报量；步骤3：建立气象数据与逐日最大臭氧八小时滑动平均值的多元线性回归方程模型；步骤4：用最小二乘法估计线性回归方程的回归系数，既残差平方和最小来估计多元线性回归方程模型中各气象预报因子的回归系数，得到确认的多元线性回归方程；步骤5：根据确认的多元线性回归方程进行历史回代，进行回归效果分析，并计算复相关系数；步骤6：对多元线性回归方程的回归效果进行显著性检验；步骤7：预报值的置信区间；步骤8：将城市预测的未来逐日气象数据输入确定的多元线性回归方程中，预测出城市未来的近地面逐日最大臭氧八小时滑动平均值。2.根据权利要求1所述的一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，其特征在于，所述步骤3中的建立的多元线性回归方程公式为：其中为预报量的估计值，x1，x2，
…
，x
p
为各气象预报因子，p为气象因子的个数，b0为回归常数，b1，b2，
…
，b
p
为各气象预报因子的回归系数。3.根据权利要求1所述的一种基于多元线性回归模型的城市近地面臭氧浓度的预测方法，其特征在于，所述步骤4中的用最小二乘法估计线性回归方程的回归系数，既残差平方和最小来估计多元线性回归方程模型中各气象预报因子的回归系数，其残差平方和公式为：其中y为预报量的实际观测值，为预报量的估计值，Q为预报量的实际观测值和估计值的残差平方和，实际观测值y既逐日最大臭氧八小时滑动平均值，i＝1，2，
……
，n，其中n为气象预报因子的时间序列长度，既历史气象数据的天数。使预报量的实际观测值y与预报量的估计值的残差平方和Q达最小值，从而确定各气象预报因子的多元线性回归方程的回归系数b1,b2,
……
，b
p
和回归常数b0。4.根据权利要求1所述的一种基于多元线性回归模...

【专利技术属性】
技术研发人员：许益超，邵振华，胡文捷，李霖，
申请(专利权)人：碧空环境科技有限公司，
类型：发明
国别省市：