基于多元线性回归分析的水位短期预测方法Word文档下载推荐.docx
- 文档编号:4640251
- 上传时间:2023-05-03
- 格式:DOCX
- 页数:4
- 大小:139.12KB
基于多元线性回归分析的水位短期预测方法Word文档下载推荐.docx
《基于多元线性回归分析的水位短期预测方法Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《基于多元线性回归分析的水位短期预测方法Word文档下载推荐.docx(4页珍藏版)》请在冰点文库上搜索。
(ChangjiangWuhanWaterwayBureau,Wuhan430014,China)
Abstract:
ToprovidemoreaccuratewaterdepthpredictiondatafortheterminalusersoftheYangtzeelectronicnavigationalchart,thispaperdiscussestheimplementationprocedureandpracticeofwaterlevelshort-termforecastingmethodbasedonthestatisticsandusingtheregressionanalysismethodbytakingtheYangtzeJianli
gaugingstationforexample,andhasachievedsatisfactoryresults.
Keywords:
regressionanalysis;
waterlevel;
short-termprediction
长江电子航道图系统是基于最新的河床地形
和实时水位数据向航行船舶提供可通航水域的智能化动态航行参考系统,能够为船舶提供准确的自身位置信息、实际水深数据、航标位置、周围航行船舶信息、最新航行通告、航道地物地貌信息、水深短期预测信息等[1]。
通过长江电子航道图
系统给营运船舶终端用户提供及时准确的航道水
深短期预测信息,能够为船舶合理配载、优化营运提供重要参考依据,具有一定的经济价值。
回归分析是研究两个变量或多个变量之间因果关系的统计方法,是一种“由因索果”的定量分析、预测技术。
其基本思想是,在相关分析基础上,对具有相关关系的两个或多个变量之间数量变化的一般关系进行测定,确立合适的数据模型,从而通过自变量的已知或设定值来估计和预测因变量均值。
本文以统计学为基础,应用回归
分析方法,以长江中游监利水位站为例,利用长
江中游宜昌至城陵矶各站点历史水位资料,建立回归分析预测模型,并运用此模型预测监利站短期内的水位变化情况。
1基于多元线性回归分析的水位短期预测方法实施步骤
某河段水位的涨落变化,主要受其上游来水、支流入汇和降雨等综合因素影响。
一方面,上游来水量的变化,会体现在上游水位的涨落变化上,而上游来水量的变化需经过一段流程时间,才能体现在其下游某河段的水位涨落上;
另一方面,某河段下游河段由于受支流入汇、降雨等多种原因产生的壅水作用,也会对其上游河段的水位变化有较大影响。
基于此思路,我们把水位涨落波动的差值作为回归分析的变量,再利用
收稿日期:
2013-04-02
作者简介:
范先友(1955—),男,高级工程师,主要从事航道管理工作。
•160• 水运工程 2014年
统计软件(SPSS软件)进行多元线性回归分析。
1.1 根据预测目标确定自变量和因变量根据需要预测监利水位的目标,把预测开始
时间作为当日,确定以预测未来1天(预测日)的监利水位与其前1日(当日)的水位差值作为因变量(y)。
通过前面分析和对长江水位统计资料的初步观察,找到可能影响监利水位涨落变化的因素,并经过初步筛选,选出主要的影响因素作为自变量。
确定以下因素作为自变量:
宜昌当日水位与其前1日水位差值(x1)、宜昌前1日水位
与前2日水位差值(x2)、沙市当日水位与其前
1日水位差值(x3)、沙市前1日水位与前2日水位差值(x4)、监利当日水位与其前1日水位差值(x5)、城陵矶当日水位与其前1日水位差值
(x6)。
1.2建立回归预测模型[2-3]
建立多元线性回归预测模型时,为了保证回归模型具有较好的解释能力和预测效果,在选择自变量时应遵循以下一些准则:
1)自变量对因变量必须有显著的影响,并呈密切的线性相关;
2)自变量与因变量之间的线性相关必须是真实的,而不是形式上的;
3)自变量之间应具有一定的互斥性,即自变量之间的相关程度不应高于自变量与因变量之间的相关程度;
4)自变量应具有完整的统计数据,其预测值容易确定。
遵照自变量的选择准则,对前文确定的自变量和因变量进行适应性分析:
1)长江监利水位站位于宜昌、沙市水位站下游,位于城陵矶水位站上游。
上游来水情况和下游壅水作用对监利水位有直接的、显著的影响,而且上下游的水位涨落关系有较密切的线性相关性。
2)以长江相邻各水位站点的水位涨落差值作为自变量和因变量,它们之间的依存关系和线性相关性是真实存在的。
3)由于长江沿程水位涨落相互影响、制约,所以确定的各自变量之间有一定的相关性,但自变量之间的相关程度低于自变量与因变量之间的相关程度。
4)长江中游重点港埠历史水位资料齐全、完
整,所有水位资料均逐日记录,其预测值也比较
容易确定,因此经计算处理可以作为自变量的统计数据。
如上所述,可以依据前文确定的自变量和因变量,建立经验回归方程,即多元线性回归分析预测模型为:
2 1 3 6 2
y=b0+b1x1+b2x2+b3x3+b4x4+b5x5+b6x6+e
(1)式中:
b0为常数项,e为随机误差,b1,b2,…,b6为回归系数。
b1为x2,x3,…,x6固定时,x1每增加一个单位对y的效应,即x1对y的偏回归系数;
同理b为x,x,…,x固定时,x每增加一个单位对y的效应,即x2对y的偏回归系数,等等。
1.3水位数据整理与计算所采用的样本水位数据为2008年9月1日—
2012年6月30日的长江宜昌、沙市、监利和城陵矶逐日水位数据资料。
按照确定的自变量、因变量和建立的回归模型,根据回归分析对水位数据的统计需要,对样本水位数据资料进行整理和计算,包括对水位资料的补缺和推算、对不可靠水位数据进行核实、对确定有误的水位值加以修正、对逐日水位涨落幅度进行差值计算、规范化
处理计算数据等。
1.4利用SPSS软件进行线性回归分析[4]
将经过规范化处理后的数据导入到SPSS软件中,选用进入法进行线性回归分析,并设置运行参数为:
1)按前文所述设置因变量和自变量;
2)选剔变量的准则为系统默认值,即使用回归方程的显著性统计量F所对应的概率P(Sig)值;
3)选入变量的显著性水平为0.05,即当
P≤0.05时,就将该变量选入回归方程;
4)剔除变量的显著性水平为0.10,即当
P≥0.10时,就该变量剔除回归方程。
通过设置相应参数,利用SPSS运行分析后,得到水位预测模型为:
y=(-4.121×
10-5)
+0.113x1-0.033x2+0.320x3+0.043x4-0.094x5+
0.469x6,具体输出结果见表1~3。
表1模型汇总
模型 R R2 调整R2 标准估计误差
1 0.907 0.822 0.822 0.075
第1期
范先友,等:
•161•
表1中,复相关系数R反映了模型自变量
(x1,…,x6)与因变量y之间线性回归关系的密切程度。
建立的模型越接近于实际,通过样本数据计算得出的构造统计量R2越接近于1。
由表中可知,模型调整R2系数为0.822,标准估计误差为
0.075,说明因变量与各自变量的相关度较高,回归模型比较接近实际情况。
表2
方差分析
模型
平方和
df
均方
F
Sig
回归
36.460
6
6.077
1070.318
0.000
残差
7.880
1388
0.006
总计
44.340
1394
表2列出了回归拟合的方差分析结果,其中包括回归平方和、残差平方和及自由度等。
通过方差分析构造统计量F可以检验回归方程的显著性,统计量F同样是通过样本数据计算得出。
从表2中可知,当预测因子为x1,…,x6时,检验值F=1070.318,显著性概率值Sig=0.000<0.05,说明x1,…,x6的回归系数不为零,模型具有统计学意义。
表3模型系数
出,自变量x1、x3、x6与因变量y的关系最密切,即
预测日监利水位的涨落幅度主要受宜昌当日水位涨落幅度、沙市当日水位涨落幅度和和城陵矶当日水位涨落幅度的影响。
1.5 计算并确定水位预测值利用前述分析得到的非标准化回归方程,通
过已知或设定自变量的值,计算得出水位涨落幅度的预测值,并对预测值进行综合分析,考虑局部流域降雨、水位大幅涨落等因素影响,估计随机误差,确定最终的预测值。
将最终确定的水位涨落幅度预测值,加上当日水位,得到未来一日的水位预测值。
1.6 模型实际预测检验回归预测模型能否用于实际预测,还必须对
模型的预测能力加以检验。
不难理解,回归方程对于样本期间来说是正确的,但是对用于实际预
测是否合适,还必须检验模型参数值的稳定性和相对于样本容量变化时的灵敏度,也必须研究确定模型是否可以用于样本数据以外的范围,具体做法是:
1)采用把增大样本容量以后模型分析的结果与原来的模型分析结果进行比较,并检验其差异
模型因变量:
y-监差0
非标准化系数 标准系数
B 标准误差试用版
t Sig
的显著性。
(常量)-4.121×
10-50.002-0.0200.984x1-宜差10.1130.0070.30317.1210.000x2-宜差2-0.0330.008-0.088-4.1140.000x3-沙差10.3200.0180.51818.2620.000x4-沙差20.0430.0130.0693.3820.001x5-监差1-0.0940.033-0.094-2.8520.004x6-城差10.4690.0200.48523.4360.000
从标准差分析,非标准化系数均具有较高可靠性;
从显著性分析,非标准化常量系数(b0)的显著性概率值Sig=0.984>0.10,不具有统计学
2)把预测模型用于样本以外某一段时间的实
际预测,并将预测值与实际观测值作一比较,然后检验其差异的显著性。
表4为增大样本容量(增加2012年7月1日—
9月30日水位数据作为样本)的模型分析系数,通过增大样本容量进行回归分析,得到的模型系数与原样本分析结果对比,差异不大,说明模型参数有很好的稳定性。
表4增大样本容量的模型系数
意义,剔除该常数项。
由此可以得到,非标准化
模型
因变量:
y-监差0
t Sig
方程为y=0.113x1-0.033x2+0.320x3+0.043x4-
0.094x5+0.469x6;
标准化方程为y=0.303x1-
1
0.088x2+0.518x3+0.069x4-0.094x5+0.485x6。
标
(常量)-0.0010.002-0.3950.693x-宜差10.1270.0060.34919.8620.000x2-宜差2-0.0230.008-0.063-2.9980.003
x3-
沙差1 0.282 0.017 0.462 16.3980.000
准系数表示为在方程其它变量不变的情况下,x1
改变一个单位,y改变0.303个单位;
x2改变一个单位,y改变-0.088个单位;
其它依此类推。
可以看
x4-沙差20.0500.0120.0824.2010.000x5-监差1-0.1160.032-0.116-3.6150.000x-城差10.4900.0200.49524.4130.000
•162• 水运工程 2014年
2 水位短期预测应用
长江电子航道图(2.0版)系统全线试运行后,根据电子航道图试运行工作需要,对长江中游水位站点开展了水位短期预测预报工作。
通过一段时间的水位预测预报,并将预测结果与实际水位观测值进行比较(图1)。
由图1可以看出,在水位涨落幅度不大的中枯水期和水位相对平稳的高洪水位期,水位预测准确性相对较高;
在水位经常有陡涨陡落的中洪水位期,预测误差较大。
总体来说,从水位实际预测应用情况看,短期内水位预测误差不大,说明模型有很好的适应性,能够用于对水位的短期预测。
13
12
3 结语
1)预测江河水位短期变化,在资料较充分的条件下,多元线性回归分析是一个有力的研究工具。
利用多元线性回归分析方法,对长江中游水位进行短期预测,具有一定的可操作性和可信度,其预报精度较高、误差较小。
2)回归预测方程是一种经验公式,其准确性决定于因变量与自变量之间的依存关系和依赖程度,选择适当的相关因子,对预测精度有很大的影响。
3)应用该方法进行水位预测,应综合考虑河道流域的自然特性和地理因素;
要进一步提高预测精度,还应结合实践经验和其它预测方法进行综合分析和验证。
ᄢѽපͮࠃᬄϘ1N
11
10
9
8
7
5
4
4 5 6 7
R2=0.998
8 9 10 11 12 13
参考文献:
[1]何明宪.长江电子航道图系统的研发与建设[EB/OL].
(2012-08-08)[2013-03-27]
main.do?
method=item&
id=43319.
[2]包凤达,翁心真.多元回归分析的软件求解和案例解读[J].数据统计和管理,2000(5):
56-61.
[3] 多元线性回归分析预测法[DB/OL].(2011-07-14)
ᄢѽᮔኃ1ܸපͮϘ1N
图1监利预测第1天水位与实际水位对比
需要说明的是,一个回归模型只是针对某一个水位站点的水位短期预测,要预测各个不同站点水位,需根据每个水位站的具体情况,分析影响因素,建立相应回归模型,进行回归分析和参数估计,进而得到各自的回归预测方程。
[2013-03-27]
NVSYwhNt4ksh9k4l6NCHUpzWDTmDyYNbp7nEKqvEkkNL3aF7zX01fjXgR5HC3c-sJQwjM3sXU1FEh-TiIF2dBw9mT-zveQqccbQ0PO3.
[4]顾伟欣,周红.基于逐步回归分析的孔隙度预测方法[J].
石油地质与工程,2008
(1):
37-39.
(本文编辑郭雪珍)
·
消 息·
中国交建2013年获国家授权专利512项再创新记录
中国交建2013年获国家授权专利总数再创新记录,全年共获国家授权专利512项,同比增长24%,其
中获发明专利70项。
公司获得授权专利数量已连续7年实现大幅增长,创新科技成果不断涌现,充分体现了公司及各单位对知识产权保护工作的高度重视,凝聚了广大科技人员的不懈努力。
摘编自《中国交通建设网》
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 多元 线性 回归 分析 水位 短期 预测 方法