《基于多因子小波神经网络模型的水文时间序列预测方法.pdf》由会员分享,可在线阅读,更多相关《基于多因子小波神经网络模型的水文时间序列预测方法.pdf(16页珍藏版)》请在专利查询网上搜索。
1、10申请公布号CN102323970A43申请公布日20120118CN102323970ACN102323970A21申请号201110130040422申请日20110519G06F19/00201101G06N3/0220060171申请人河海大学地址210098江苏省南京市鼓楼区西康路1号72发明人朱跃龙李士进王继民范青松冯钧万定生74专利代理机构南京经纬专利商标代理有限公司32200代理人许方54发明名称基于多因子小波神经网络模型的水文时间序列预测方法57摘要本发明公开了一种基于多因子小波神经网络模型的水文时间序列预测方法。本发明提出了一种多因子小波神经网络预测模型用于水文时间序列预。
2、测,该模型以多时间序列信息作为输入,不仅包括预测目标时间序列的当前小波系数,还包括与此时间序列相关的其它时间序列的当前小波系数,并根据其与预测目标时间序列之间互信息作为评判两者相关性的度量,选择相关性强的其它时间序列;进一步利用基于加权相关系数的小波函数选择标准,为该模型选择最优的小波函数。相比现有技术,本发明方法具有更高的预测准确性,以及更好的可扩展性和实用价值。51INTCL19中华人民共和国国家知识产权局12发明专利申请权利要求书2页说明书11页附图2页CN102323981A1/2页21一种基于多因子小波神经网络模型的水文时间序列预测方法,该方法首先根据待预测的水文时间序列建立小波神经。
3、网络水文时间序列预测模型,然后按照建立的预测模型进行水文时间序列预测,其特征在于,所述小波神经网络水文时间序列预测模型为多因子小波神经网络水文时间序列预测模型,以多时间序列信息作为输入,不仅包括预测目标时间序列的当前小波系数,还包括与此时间序列相关的其它时间序列的当前小波系数。2如权利要求1所述基于多因子小波神经网络模型的水文时间序列预测方法,其特征在于,所述与预测目标时间序列相关的其它时间序列,是根据其与预测目标时间序列之间互信息作为评判两者相关性的度量,互信息值较大,则相关性强;具体按照以下方法确定步骤1、提供与待预测序列O可能有关的若干原始输入序列I1,I2,IN,N是输入序列数目;步骤。
4、2、对待预测序列O和原始输入序列I1,I2,IN进行离散化,得到离散序列DO,D1,D2,DN;步骤3、分别计算DO与D1,D2,DN之间的互信息,结果记为M1,M2,MN步骤4、根据M1,M2,MN,选择MITH的原始输入序列II作为待预测序列小波神经网络建模时的相关序列,I为1到N之间的整数;TH为预先设定的阈值。3如权利要求2所述基于多因子小波神经网络模型的水文时间序列预测方法,其特征在于,采用等概率划分法对待预测序列O和原始输入序列I1,I2,IN进行离散化。4如权利要求2所述基于多因子小波神经网络模型的水文时间序列预测方法,其特征在于,所述阈值TH的值为015。5如权利要求14任一项。
5、所述基于多因子小波神经网络模型的水文时间序列预测方法,其特征在于,所述多因子小波神经网络水文时间序列预测模型中所使用的小波函数按照以下方法确定步骤1、分别使用待选的小波函数构建小波神经网络水文时间序列预测模型并进行预测;步骤2、对于每一种待选的小波函数,分别按照以下方法得到相关系数向量对不同水平上的小波分解序列,分别统计预测目标时间序列的小波序列自相关系数,以及预测目标时间序列和与预测目标时间序列相关的其它时间序列的小波系数序列的相关系数;最后采用加权的方法综合各个水平上的相关系数,得到相关系数向量;步骤3、根据采用各小波函数预测所得到的相关系数向量,确定最终使用的小波函数。6如权利要求5所述。
6、基于多因子小波神经网络模型的水文时间序列预测方法,其特征在于,两个时间序列的相关系数,按照以下方法得到假设两个时间序列分别为、,时间序列的长度为N;则这两个时间序列之间的相关系数按照下式计算,权利要求书CN102323970ACN102323981A2/2页3式中,是两个序列时滞为T的相关系数,是第一个序列的均值,是第二个序列的均值;单个时间序列的时滞为T的自相关系数,按照以下方法得到将该时间序列向右移动T个位置生成一个新的时间序列,然后按照上式计算该时间序列与新生成的时间序列的相关系数,得到的相关系数即为该时间序列时滞为T的自相关系数。权利要求书CN102323970ACN102323981。
7、A1/11页4基于多因子小波神经网络模型的水文时间序列预测方法技术领域0001本发明涉及一种复杂时间序列预测方法,尤其涉及一种基于多因子小波神经网络模型的水文时间序列预测方法,属于水文预报技术领域。背景技术0002时间序列数据挖掘研究主要包括预测、分类、相似性搜索及序列模式挖掘,而复杂时间序列预测是数据挖掘领域具有挑战性的问题之一。解决复杂时间序列预测问题的一个较好的方法是基于小波神经网络模型的时间序列预测方法。0003小波分析是FOURIER分析发展史上的一个里程碑式的进展,具有时、频同时局部化的优点,因此被誉为数学“显微镜”。与傅立叶分析仅能够提供频域表示相比较,小波变换能够同时提供时频局。
8、部化特征,且克服了短时傅立叶变换的分辨率有限的缺点。离散小波分解可以对时间序列进行多尺度分解,能够提取不同频率区间的成分序列,实现对时间序列的分频研究。借助于小波的多分辨分解功能,可以从原始序列得到分辨率“从粗到细”的多尺度特征,通过对这些小波特征的学习可以对时间序列的潜在的多因素变化过程有更加显性的描述,所以网络更加容易“抓住”输入和输出数据之间的内在的联系和规律。相反,基于传统的单分辨率学习的网络通常很难描述复杂的时间序列,例如混沌序列,其收敛速度慢,泛化能力差。小波变换将复杂的时间序列分解成若干不同频率的细节信号和背景信号,细节信号是高频部分,背景信号是低频部分,它们是原始时间序列在不同。
9、频率区间上的表现,相对原始序列的比重不同,因而对原始时间序列的预测起到的作用不同。小波网络模型结合了小波分析的在信号处理上的时间域频率域的多分辨性,和神经网络的强大的非线性逼近特性,兼有二者的优点。以上述方式建立的模型称为小波网络模型WAVELETNETWORKMODEL,简记为WNN。0004很多学者对小波网络模型对不同领域的时间序列预测进行了研究,例如CHENYUEHUI,YANGBO,DONGJIWENTIMESERIESPREDICTIONUSINGALOCALLINEARWAVELETNEURALNETWORKNEUROCOMPUTING,2006,696449465、DASHPK,。
10、NAYAKMAYA,SENAPATIMR,ETAL,MININGFORSIMILARITIESINTIMESERIESDATAUSINGWAVELETBASEDFEATUREVECTORSANDNEURALNETWORKS,ENGINEERINGAPPLICATIONSOFARTIFICIALINTELLIGENCE,2007,202185201、GANXIAOBING,LIUYING,ANDAUSTINFRANCISRAPREDICTIONMETHODFORTIMESERIESBASEDONWAVELETNEURALNETWORKSPROCEEDINGSOFCIS2005,PP902908、。
11、BENAOUDAD,MURTAGHFIONN,ANDSTARCKJEANLUC,ETALWAVELETBASEDNONLINEARMULTISCALEDECOMPOSITIONMODELFORELECTRICITYLOADFORECASTINGNEUROCOMPUTING,2006,7013139154、CHENDI;FENGHAIHANG;LINQINGJIA;ETAL,MULTISCALEINTERNETTRAFFICPREDICTIONUSINGWAVELETNEURALNETWORKCOMBINEDMODEL,PROCEEDINGSOFFIRSTINTERNATIONALCONFERE。
12、NCEONCOMMUNICATIONSANDNETWORKINGINCHINA,2006,PP15、李爱国,覃征,滑动窗口二次自回归模型预测非线性时间序列,计算机学报,2004,270710041008说明书CN102323970ACN102323981A2/11页5等。但这些文献仅仅是利用小波分解对单道时间序列进行分解。仅仅依靠预测时间序列自身提供的信息,对于预测精度提高的作用具有一定局限性。而在物理意义上,一个系统内的多时间序列往往具有相关性。0005水文预报例如径流预报对水资源的调配管理和防洪减灾的决策有着重要的意义。然而由于气候和下垫面等因素的综合影响,径流表现出复杂的非线性和非平稳性。
13、特征,一般包括确定性成分和随机成分。神经网络具有很强的处理大规模复杂非线性动力学系统的能力,被广泛用于河道水情预报,可以识别水流运动变化过程与其影响因子之间的复杂非线性关系,为河道水情预报提供了一条新的途径。径流序列自身的组成比较复杂,在一个径流中含有多种频率成份,所以有必要对其进行分频率研究,而小波变换的方法提供了一种便利的时频分析技术。小波和神经网络结合使用的方法近年来得到广泛关注,利用此组合模型进行水文预测成为研究热点。文献王文圣,熊华康,丁晶,日流量预测的小波网络模型初探,水科学进展,2004,153382386中提出了一种小波网络模型用于对日流量进行预测。该小波网络模型是在三层神经网。
14、络的基础上,以T时刻原始信号分解的小波系数作为输入,以TT时刻的原始信号FTT为输出而构建的模型T是预见期。然而该小波网络模型也仅利用单时间序列分解出来的小波系数,对未来进行预测,模型的预测能力是有限的。由于从原始时间序列分解出来的小波系数,已经包含对时间序列前期变化中长时间记忆和短时间记忆,加入更多的前期小波系数往往不能够对预测目标起到更好的辅助作用。而且前期小波系数的加入,使得模式输入个数增加,相应的训练模式集的复杂度也增加,对神经网络结构的确定和训练时间提出更高要求。对于更加复杂时间序列,仅仅根据本时间序列的前期小波系数,不能对时间序列未来的演变进行更好的刻画。径流实际上跟系统内的其他水。
15、文要素具有相关关系,如上下游关系站或相关站的水文要素包括流量、水位、气温和降水等对径流产生具有明显的影响关系。0006不同领域的时间序列数据具有不同特性,仅利用同一种小波函数进行分解是不合适的。文献桑燕芳,王栋,水文序列小波分析中小波函数选择方法,水利学报,2008,393295300从序列重构角度,应用蒙特卡罗方法,通过模拟数据以及黄河利津站和浙江白溪流域的实测水文序列,探讨小波分析影响小波函数选择的因素,进而建立合理选择小波函数的依据。研究结果表明序列自身的变化特性是小波函数选择的重要影响因素。文献刘素一,权先璋,张勇传,不同小波函数对径流分析结果的影响,水电能源科学,2003,21129。
16、31通过训练多尺度预测模型的实验方法,对比了几种小波函数对多尺度框架下径流预测的影响。对于不同数据,特别是大量数据情况下,使用这种类似特征选择中的WRAPPER方法确定最优的小波函数比较耗时。发明内容0007本发明所要解决的技术问题在于克服现有技术以单一时间序列作为小波神经网络模型的输入进行水文时间序列预测的不足,提供一种基于多因子小波神经网络模型的水文时间序列预测方法。0008本发明的基于多因子小波神经网络模型的水文时间序列预测方法,首先根据待预测的水文时间序列建立小波神经网络水文时间序列预测模型,然后按照建立的预测模型进行水文时间序列预测,其特征在于,所述小波神经网络水文时间序列预测模型为。
17、多因子小说明书CN102323970ACN102323981A3/11页6波神经网络水文时间序列预测模型,以多时间序列信息作为输入,不仅包括预测目标时间序列的当前小波系数,还包括与此时间序列相关的其它时间序列的当前小波系数。0009进一步地,所述与预测目标时间序列相关的其它时间序列,是根据其与预测目标时间序列之间互信息MUTUALINFORMATION作为评判两者相关性的度量,互信息值较大,则相关性强;具体按照以下方法确定0010步骤1、提供与待预测序列O可能有关的若干原始输入序列I1,I2,IN,N是输入序列数目;0011步骤2、对待预测序列O和原始输入序列I1,I2,IN进行离散化,得到离。
18、散序列DO,D1,D2,DN;0012步骤3、分别计算DO与D1,D2,DN之间的互信息,结果记为M1,M2,MN;0013步骤4、根据M1,M2,MN,选择MITH的原始输入序列II作为待预测序列小波神经网络建模时的相关序列,I为1到N之间的整数;TH为预先设定的阈值。0014更进一步地,所述多因子小波神经网络水文时间序列预测模型中所使用的小波函数按照以下方法确定0015步骤1、分别使用待选的小波函数构建小波神经网络水文时间序列预测模型并进行预测;0016步骤2、对于每一种待选的小波函数,分别按照以下方法得到相关系数向量0017对不同水平上的小波分解序列,分别统计预测目标时间序列的小波序列自。
19、相关系数,以及预测目标时间序列和与预测目标时间序列相关的其它时间序列的小波系数序列的相关系数;最后采用加权的方法综合各个水平上的相关系数,得到相关系数向量;0018步骤3、根据采用各小波函数预测所得到的相关系数向量,确定最终使用的小波函数。0019其中,两个时间序列的相关系数,按照以下方法得到0020假设两个时间序列分别为、,时间序列的长度为N;则这两个时间序列之间的相关系数按照下式计算,00210022式中,是两个序列时滞为T的相关系数,是第一个序列的均值,是第二个序列的均值;0023单个时间序列的时滞为T的自相关系数,按照以下方法得到0024将该时间序列向右移动T个位置生成一个新的时间序列。
20、,然后按照上式计算该时间序列与新生成的时间序列的相关系数,得到的相关系数即为该时间序列时滞为T的自相关系数。0025本发明提出了一种多因子小波神经网络预测模型用于水文时间序列预测,该模型以多时间序列信息作为输入,不仅包括预测目标时间序列的当前小波系数,还包括与此时间序列相关的其它时间序列的当前小波系数,并根据其与预测目标时间序列之间互信息作说明书CN102323970ACN102323981A4/11页7为评判两者相关性的度量,选择相关性强的其它时间序列;进一步利用基于加权相关系数的小波函数选择标准,为该模型选择最优的小波函数。相比现有技术,本发明方法具有更高的预测准确性,以及更好的可扩展性和。
21、实用价值。附图说明0026图1为本发明的多因子小波神经网络预测模型的结构示意图;0027图2为本发明确定相关的其它时间序列的流程示意图;0028图3为本发明选择最优小波函数的流程示意图。具体实施方式0029下面结合附图对本发明的技术方案进行详细说明0030本发明的多因子小波神经网络预测模型以多时间序列信息作为输入,不仅包括预测目标时间序列的当前小波系数,还包括与此时间序列相关的其它时间序列的当前小波系数,其结构如图1所示,其中要预测的站点为A,与A相关的站点为B至I。也可将此多因子小波神经网络预测模型称为河道小波网络模型,简记为RWNN。0031本发明的基于多因子小波神经网络模型的水文时间序列。
22、预测方法,按照以下步骤进行0032步骤10、确定相关输入时间序列;0033在信息论中,熵HX被用来度量随机变量X的不确定性。假设X是一离散随机变量,熵HX被定义为00340035与熵只能衡量单一变量的不确定性不同,互信息可用来衡量两个以上随机变量之间共有的信息量。假设Y是另一离散随机变量,则X和Y之间的互信息IX;Y被定义为00360037可以证明互信息和熵有如下关系0038IX;YHXHX|YHYHY|X,0039IX;YHXHYHX,Y0040显然,两个随机变量之间互信息越大,它们之间相关性越强。关于互信息的详细内容可参见文献COVERANDTHOMAS,1991TMCOVERANDJAT。
23、HOMAS,ELEMENTSOFINFORMATIONTHEORY,WILEY,NEWYORK1991。从时间序列预测角度看,两个序列之间互信息越大,则通过其中一个序列预测另一个序列越准确。并且,由于小波变换之前的序列之间可能存在非线性相关,而互信息恰好能够反映不同变量之间的非线性相关性。因此,本发明将时间序列之间互信息作为评判两者相关性的度量,提出基于互信息的输入序列确定方法,如图2所示,具体按照以下方法确定相关的输入时间序列0041步骤11、提供与待预测序列O可能有关的若干原始输入序列I1,I2,IN,N是输入序列数目;0042步骤12、对待预测序列O和原始输入序列I1,I2,IN进行离散。
24、化,得到离散说明书CN102323970ACN102323981A5/11页8序列DO,D1,D2,DN;常用的连续变量离散化方法主要有两种等距离划分法和等概率划分法,一般后者效果比较好,本发明优选文献FRASERANDSWINNEY,1986AMFRASERANDHLSWINNEY,INDEPENDENTCOORDINATESFORSTRANGEATTRACTORSFROMMUTUALINFORMATION,PHYSREVA3321986,PP11341140中提出的等概率划分法对时间序列进行离散化;0043步骤13、分别计算DO与D1,D2,DN之间的互信息,结果记为M1,M2,MN;00。
25、44步骤14、根据M1,M2,MN,选择MITH的原始输入序列II作为待预测序列小波神经网络建模时的相关序列,I为1到N之间的整数;TH为预先设定的阈值,该阈值优选015。0045步骤20、对确定的相关输入序列和目标序列进行小波分解,根据加权相关系数进行小波选择;0046小波分解是一种有效的问题分解和特征提取工具,下表1列举了目前常见的几种小波函数及其特性。0047表1常用小波函数0048小波正交性滤波系数长度支撑宽度对称性HAARYES21NODB2YES43FARFROMDB3YES65FARFROMBIOR44NO98YESB3SPLINENO54YES0049但是在具体领域中,存在不同。
26、性质的时间序列数据,用不同的小波函数对这些数据进行分解,所产生的成分序列对建立的预测模型性能会有不同程度的影响。由于小波滤波器具有不同的支撑长度和自身的滤波性质,所以应该选择一个合适的小波对具体的时间序列进行分解,使得分解后的小波系数序列可以提高预测模型的性能。本发明使用T时刻的不同时间序列的小波系数作为模型的输入,而输出为预测时间序列的TT原始数据。对于多因子预测模型,设T时刻的模型输入包括预测时间序列的小波系数和扩展的输入输出为原始时间序列的未来值XATT。预测目标XATT也含有各个频率区间上的对应成分,即不同水平上的小波系数0050根据时间序列分析中关于平稳时间序列预测理论具体可参见文献。
27、PBROCKWELL,ADAVIS,TIMESERIESTHEORYANDMETHODS,2NDEDITION,NEWYORK,SPRINGER,SPRINGERSERIESINSTATISTICS,2002,设有预测目标时间序列XT,T1,N的T时间间隔的自相关系数又存在与XT,T1,N相关的时间序列YT,T1,N且两者之间的T时间间隔的相关系数如果和越大,那么预说明书CN102323970ACN102323981A6/11页9测目标时间序列XT,T1,N的预测方差就越小。0051根据上述思想,本发明针对所提出的多因子小波神经网络模型,提出了一种基于加权相关系数的小波函数选择标准,对不同水平。
28、上的小波分解序列,统计预测目标时间序列的小波序列自相关系数和预测目标时间序列和多因子时间序列的小波系数序列相关系数,最后采用加权的方法综合各个水平上的相关系数,如图3所示,具体按照以下方法选择最优的小波函数0052步骤21、分别使用待选的小波函数构建小波神经网络水文时间序列预测模型并进行预测;0053步骤22、对于每一种待选的小波函数,分别按照以下方法得到相关系数向量0054对不同水平上的小波分解序列,分别统计预测目标时间序列的小波序列自相关系数,以及预测目标时间序列和与预测目标时间序列相关的其它时间序列的小波系数序列的相关系数;最后采用加权的方法综合各个水平上的相关系数,得到相关系数向量;0。
29、055步骤23、根据采用各小波函数预测所得到的相关系数向量,确定最终使用的小波函数。0056其中,两个时间序列的相关系数,按照以下方法得到0057假设两个时间序列分别为、,时间序列的长度为N;则这两个时间序列之间的相关系数按照下式计算,00580059式中,是两个序列时滞为T的相关系数,是第一个序列的均值,是第二个序列的均值;0060单个时间序列的时滞为T的自相关系数,按照以下方法得到0061将该时间序列向右移动T个位置生成一个新的时间序列,然后按照上式计算该时间序列与新生成的时间序列的相关系数,得到的相关系数即为该时间序列时滞为T的自相关系数。0062例如,对于输入为和输出为XAT1的预测模。
30、型,计算方法如表2,利用最后一行的加权向量来比较不同小波函数的分解结果。0063表2二因子小波网络输入的小波函数选择指标计算方法说明书CN102323970ACN102323981A7/11页1000640065假设一个时间序列为,另一相关时间序列为,时间序列的长度为N。则有相关系数定义如公式1,具体计算方法如公式2。006600670068计算单道序列的时滞为T自相关系数,只需要令向右移动T个位置生成新的序列,然后以新序列代替按照公式2计算即可。0069步骤30、采用选择的小波函数,对输入时间序列进行小波分解,建立多因子小波神经网络模型;0070步骤40、训练多因子小波神经网络模型;0071。
31、步骤50、根据给定输入时刻的多序列小波变换系数,作为多因子小波神经网络模型的输入,获得最终的预测值。0072下面以王家坝站的汛期日流量预测为例进一步说明本发明的技术方案及有益效果。0073首先用本发明的基于加权相关系数的小波函数选择标准进行小波函数的选择。小波函数的选择应当和具体研究对象有关,对比不同年份王家坝站的汛期日流量过程,发现该站2003年日流量序列较为典型,选择2003年汛期日流量序列作为小波函数适应性选择的数据资料,实验资料还包括王家坝站的三个主要上游站班台、息县和潢川的2003年同期的汛期日流量数据。本实验选用如下小波对上述四个站点的2003年日流量序列数据进行三水平小波分解,参。
32、与比较的小波与其相应的低通滤波器如表3。0074表3小波函数及其低通滤波器0075说明书CN102323970ACN102323981A8/11页110076因为王家坝站点的日流量序列跟自身的成分相关性较大,所以令其权重系数WJB06,又因为班台和潢川站均在小支流上,息县站在大支流上,所以令XIX02,BT01,HC01。即权重分配为BT,XIX,HC,WJB01,02,01,06。令T1时,各个不同的小波函数对时间序列分解结果的相关性统计结果见表4由于篇幅所限,省略了加权之前的具体数据,文中只给出加权后的相关系数。0077表4各小波分解后加权相关系数0078D1D2D3C3B3SPLINE0。
33、4933076810853809627HAAR03813070750829009572DB202286062990800909492DB301652062170803009427COIFL01835061720798209470BIOR1302875068430821009518BIOR2203268049970762709245BIOR4401669056110787309419RBIO440173006472080800944000790080从表4可知,B3SPLINE小波分解后序列与目标序列各对应子序列相关系数最大,说明书CN102323970ACN102323981A9/11页12H。
34、AAR小波比B3SPLINE的结果略差。我们还可以发现,DB2比HAAR小波的统计结果要差,这说明同一家族中支撑宽度大的小波分解出来的小波系数,并不会使预测更加容易。如果选择一个不适合的小波对王家坝等四站的日流量数据进行分解,将导致很差的统计结果,例如BIOR22小波,在各个水平上的统计结果都较差,说明其对此三站日流量数据进行分解未能有效提取出对预测有利的特征。0081以王家坝站日流量作为预测对象,选取与之主要相关的三个观测站班台、息县、潢川,此三站分别为上游三大支流上的站点,数据时间跨度为2000年2007年,仅包括每年汛期的日流量数据,具体为20002005年数据区间均为61至101日,2。
35、006年数据区间51至101日,2007年数据区间515至101。本文将2000年至2005年的数据作为训练集,2006和2007年的数据作为测试集。0082一河道小波网络模型的建立0083根据上面的小波函数选择结果,HAAR小波和B3SPLINE小波分解的数据对王家坝站的河道小波网络模型的适应度在相关系数的意义上是较好的,所以首先采用这两种小波对王家坝、班台、息县和潢川四站的日流量序列的分别做3水平分解。然后以此四站的当日小波系数作为神经网络的输入,以王家坝站的预测目标1D、2D、3D预测期的日流量作为输出,对两种小波分解的预测结果进行比较。0084对神经网络模型训练之前,对输入作MIDRA。
36、NGE标准化,对输出流量做自然对数变换后也将其规范化到1,1,然后利用MATLAB的神经网络工具箱来实现神经网络模型的训练。在标准BP算法的基础上,使用LM数值优化方法来加快神经网络训练速度,并使用贝叶斯正则化方法保证训练网络的泛化能力。采用试错法确定最优的神经网络隐层节点数目。对于应用HAAR小波的神经网络模型,1D、2D和3D预测期最优隐层节点分别为4、5和5;而对于应用B3SPLINE小波的神经网络模型,1D、2D和3D预测期最优隐层节点分别为5、4和6。然后对不同的初始权值,采用多次训练的方法选取产生最优模型。最后用训练得到的网络在测试集上进行模拟,产生预测输出。网络的原始输出需要进行。
37、反MIDRANGE标准化,得到数据Y,由于原始流量作为输出做了自然对数变换,所以需要对Y进行反自然对数变换,最终得到原始数据的预测结果。0085二结果评价指标0086对于时间序列预测,常用的评价标准有平均绝对误差MAE和平均相对误差MRE,其计算公式如公式3和4。008700880089其中PI,AI分别是T个观测中第I个观测的预测值和实际值。0090另外,按照水文情报预报规范合格率定义,计算相对误差在20之内的合格预测数目,得到合格率,记为ER_ZP20。又由于是洪水预报,我们还考虑了预测方向正确率以及发生洪水时的高流量方向正确率。方向正确率,又叫方向变差对称值DVS,如式5,其中P为预测期。
38、。为了比较高流量方向判断的准确度,对流量大于500的预测统计其方向正确率,记为DVS_HF。说明书CN102323970ACN102323981A10/11页1300910092三HAAR小波和B3SPLINE小波分解建模结果比较0093表5给出了基于HAAR小波分解和B3SPLINE小波分解的小波网络模型在测试集上的预测结果。0094表5HAAR和B3SPLINE小波网络模型在测试集上不同预测期结果的比较00950096从表5可知,在1D预测期上,B3SPLINE的ER_ZP20指标要好于HAAR的指标;在MAE和MRE指标上基于B3SPLINE小波分解的小波网络模型要比基于HAAR小波分解。
39、的模型稍差;B3SPLINE的DVS和DVS_HF更好些;在2D预测期上,B3SPLINE的MRE和合格率指标均要好于HAAR的指标,但是在合格率仅约为60的情况下,从DVS和DVS_HF指标上看,HAAR小波的预测结果稍好于B3SPLINE小波的结果;在3D预测期上,在合格率很低的情况下,DVS_HF显示B3SPLINE小波结果稍优于HAAR小波的结果。0097总体上B3SPLINE小波预测结果要好于HAAR小波,这也验证了本发明的小波函数选择方法的正确性。0098四与单序列的小波网络模型结果比较0099表6显示了本发明的河道小波网络模型RWNN和小波网络模型WNN在测试集合上在不同评价指标。
40、上的统计结果。0100表6测试集上WNN和RWNN在不同预测期上的结果比较0101说明书CN102323970ACN102323981A11/11页140102从表6可知,RWNN模型在不同预测期上的所有指标均都优于WNN模型,这充分体现了扩展模型预测性能的提高。按照水文情报预报规范评价标准,对于1D预测期,河道小波网络模型满足甲级预报的标准,小波网络模型仅满足乙级预报的标准;对于2D的预测期,河道小波网络模型满足丙级预报标准,可用于参考性预报,而小波网络模型精度低于丙级;对于3D预测期情况,两个模型精度均低于丙级,但是从DVS和DVS_HF可以看出,河道小波网络模型在方向预测性能上表现更好,特别是在高流量的方向预报正确率DVS_HF上,河道小波网络模型高于小波网络模型15个百分点,适合于做3D预测期的参考性估报。同时,在DVS_HF指标上,1D预测期的RWNN比WNN高出12个百分点,达到86,能更加可靠地回答在洪峰过程中关注的“次日流量是升高还是降低”的问题。说明书CN102323970ACN102323981A1/2页15图1图2说明书附图CN102323970ACN102323981A2/2页16图3说明书附图CN102323970A。