第12章 神经网络预测法

1、第第12章章 神经网络预测法神经网络预测法 12.1 神经网络基本理论 12.2 BP神经网络的MATLAB函数 12.3 案例分析 12.3.1 北京市房地产开发投资及销售分析练习与提高（练习与提高（12） 12.3.2 深证综合指数预测 12.1 神经网络基本理论12.1.1 人工神经网络人工神经网络人工神经网络(Artificial Neural Network,ANN)，是由大量处理单元(神经元Neurons)广泛互连而成的网络,由分布于若干层的节点组成。每个单节点都有自己的输入值、权重、求和与激活函数以及输出值，在处理之前，数据被分为训练数据集(Training Data set)和

2、测试数据集(Testing Data set)，然后将权重或输入，指派到第一层的每一个节点。每次重复时，系统处理输入，并与实际值相比较，得到度量后的误差，并反馈给系统，调整权重。大多数情形下，调整后的权重都能更好地预测实际值。当达到预定义的最小误差水平时，处理结束。12.1.2 BP神经网络的基本原理 BP(Back-Propagation Network) 是一种多层网络的“逆推”学习算法。其基本思想是: 学习过程由信号的正向传播与误差的反向传播组成。 正向传播时，输入样本从输入层传入，经隐层逐层处理后传向输出层。若输出层的实际输出与期望输出不符，则转向误差的反向传播阶段。 误差的反向传播是

3、将输出误差以某种形势通过隐层向输入层逐层反传，并将误差分摊给各层的所有单元，从而获得各层单元的误差信号，此误差信号即作为修正各单元权值的依据。 这种信号正向传播与误差反向传播的各层权值调整过程是周而复始地进行。权值不断调整的过程，也就是网络的学习训练过程。此过程一直进行到网络输出的误差减少到可以接受的程度，或进行到预订的设定的学习次数为止。12.1.3 BP神经网络过程1BP网络的结构BP网络是一种单向传播的多层前向网络，具有三层或三层以上的神经网络，包括输入层、中间层(隐层)和输入层。上下层之间实现全连接，每一层神经元之间无连接。 输入层 隐层 输出层2传递函数或激活函数一般输入层和隐层的传

4、递函数是S形函数（logsig）：正切S形函数（tansig）1( )1xf xe1( )1xxef xe输出层的是线性函数，用purelin表示3模拟过程首页首页网络通过对已知信息的反复学习训练，运用根据误差来逐步调整与改变神经元连接权重和神经元阈值的方法，使得相似的输入有相似的输出，从而达到处理信息、模拟输入输出关系的目的。（1）读入样本、设定初始权值和阈值；（2）设定参数；（3）计算隐含层输出；（4）计算输出层输出；（5）计算输出值与期望值的误差；（6）判断误差是否小于设定值，是则结束；（7）调整隐层到输出层的权值和阈值；（8）调整输入层到隐层的权值和阈值；（9）返回计算隐含层输出。4学

5、习过程首页首页（1）读入样本；（2）数据处理；（3）创建网络；（4）设定参数；（5）训练网络；（6）模拟输出；（7）调整参数：学习速率、动量系数、训练次数、误差精度等；（8）仿真预测：网络固定，输入新的样本集，模拟输出。12.1.4 BP神经网络预测神经网络预测首页首页BP神经网络通过对以往历史数据的学习，找出数据的变化趋势之间的非线性关系，并将其存储在网络具体的权值和阀值中，从而预测未来数据的走势。 用神经网络进行预测，即用神经网络通过一组数据点 来拟合函数f，得出未来 （k1）时刻数据的预测值。1(,)n m knnn mXf XXX1,nnn mXXXnmk1单步预测首页首页当k=1时，

6、且网络的所有输入数据都是时间序列的实际观测值时所做的预测就是单步预测。在进行预测时，把实际的时序观测值 ，这m个数据输入网络，输出是下一时刻的预测值 1,nnn mXXX1n mX若要继续对 的值进行预测，则用实际观测值 作为输入数据，得到预测值 2n mX121,nnn mXXX2n mX2多步预测首页首页当kl时，网络输入m个历史数据，输出 的预测值。多步预测用于股票价格预测误差较大。这是因为在网络运行调整权值和阀值时，每次迭代都要累加前一次k个预测值的误差，从而造成网络难以收敛的情况，甚至导致网络发生振荡。12,n mn mn m kXXX首页首页3滚动预测滚动预测，又可称为迭代一步预测

7、，是先进行单步预测，再将网络输出的预测值反馈给网络输入端作为输入的一部分，用于下一步的预测。若开始预测时输入数据是实际的时序观测值 输出是下一时刻的预测值 ，1,nnn mXXX1n mX将 与 一起作为输入数据对 项进行估计，得到输出的预测值 1n mX12,nnn mXXX2nmX2n mX如此反复迭代，就得到对未来一段时期的预测值。12.2 BP神经网络的MATLAB函数1数据的预处理和后处理数据的预处理和后处理是有效训练神经网络的关键步骤，直接影响到训练后神经网络的性能。常见的方法是将原始数据进行归一化处理，即通过一定的线性变换将输入和输出数据统一限制在0，1或-1，1区间内。利用pr

8、emnmx或prestd函数可以对输入和目标数据集进行归一化处理，使其落入-1，1区间。首页首页格式：格式： Pn，minp，maxp=premnmx(P)Pn，minp，maxp，Tn，mint，maxt=premnmx(P,T)说明：说明：premnmx函数用于对网络的输入数据或目标数据进行归一化，归一化后的数据将分布在-1，1区间内。归一化公式为：Pn=2*(P-minp)(maxp-minp)-l Tn=2*(T-mint)(maxt-mint)-l其中，P为原始输入数据，maxp和minp分别是P中的最大值和最小值，Pn为归一化后的输入数据。T是原始目标数据，maxt和mint分别是

9、T的最大值和最小值，Tn是归一化后的目标数据。格式：格式：P,T=postmnmx(Pn，minp，maxp，Tn，mint，maxt) 说明：说明：Postmnmx函数可将premnmx函数所归一化数据进行反归一化处理P=0.5*(Pn+1)*(maxp-minp)+minpT=0.5*(Tn+1)*(maxt-mint)+mint 将输入数据或目标数据转化为区间0,1的归一化处理公式为：Pn=(P-minp)(maxp-minp)Tn=(T-mint)(maxt-mint)其对应的反归一化处理公式 P=Pn*（maxp-minp)+ minp T=Tn*（maxt-mint)+ mint2

10、.创建网络（1）newff函数：用来建立一个前馈BP网络格式：格式：net=newff(PR，SN，TF，BTF，BLF，PF)；说明：说明：PR：表示由每组输入（共P组）元素的最大值和最小值 组成的P2维矩阵；或用函数minmax(P)表示；SN：表示网络隐含层和输出层神经元的个数；TF：表示网络隐含层和输出层的传递函数，tansig（默认）， logsig，purelin；BTF：表示网络的训练函数。普通训练traingdm：需设定 学习速率、动量系数，快速训练trainlm(默认)：BLF：表示网络权值学习函数，learngdf(默认)；PF：表示网络性能函数，mse(默认)，网络输出和

11、目标输出的均方误差。3设定参数net=init(net)；初始化网络权值和阈值（可不设定)net.trainparam.show=训练状态的显示幅度；（默认25）net.trainparam.lr=学习速率；（权值阈值的调整幅度）net.trainparam.mc=动量系数；（权阈值改变的重复度）net.trainparam.epochs=训练次数；（默认100）net.trainparam.goal=误差精度；（默认0）net.trainparam.time=训练秒数；（可不选）（4.训练网络：格式：格式：net，tr=train(net，P，T)说明：说明：P为输入样本矢量集；T为对应的目

12、标样本矢量集：等号左右两侧的net分别用于表示训练得到和训练以前的神经网络对象；tr存储训练过程中的步数信息和误差信息，并给出网络误差实时变化曲线。5.BP神经网络的仿真格式：格式：Y,Pf,Af,E,perf=sim(net，P，Pi，Ai，T)说明：说明：输入net为神经网络对象，P为网络输入，Pi为输入延迟的初始状态，Ai为层延迟的初始状态，T为目标矢量，Y为网络输出，Pf为训练终止时的输入延迟状态，Af为训练终止时的层延迟状态，E为输出和目标矢量之间的误差，perf为网络性能值。首页首页6.模拟输出图形输出：plot（横坐标，纵坐标，参数）查看参数：权值：net.IW（层序号）阈值：n

13、et.b（层序号）12.3 案例分析12.3.1 北京市房地产开发投资及销售分析北京市1999年至2010年房地产开发投资的资金来源中的国内贷款，与各种房屋类型销售额数据如表12-1所示，进而计算出2000年至2010年各项指标的增长率，数据如表12-2所示，试将住宅、办公楼、商业营业用房的销售额增长率作为输入元素，国内贷款额增长率作为目标函数，建立BP神经网络仿真模拟。日期国内贷款住宅销售额办公楼销售额商业营业用房销售额1999165.22232.0266.087.622000238.82409.3451.804.252001340.30531.7161.7312.462002382.777

14、16.5357.6825.952003586.86789.1640.5851.782004549.961085.1199.2150.592005676.921739.94249.62120.852006841.421626.30352.94162.5620071063.211845.97402.54237.072008889.371201.37230.72192.7620092367.772486.77431.14299.8620101439.082060.52487.32318.99 北京市房地产开发投资及销售原始数据 北京市房地产开发投资及销售额增长率日期国内贷款住宅销售额办公楼销售额商业

15、营业用房销售额20000.36840.5677-0.2435-0.583920010.35410.26160.17541.075620020.11760.2983-0.06790.733620030.42740.0965-0.35160.69082004-0.06490.31850.8940-0.023320050.20770.47220.92270.870820060.2175-0.06750.34640.296520070.23400.12670.13150.37732008-0.1785-0.4295-0.5566-0.206920090.97920.72750.62520.441920

16、10-0.4979-0.18800.12250.0618案例二：案例二： 中空保温玻璃的销售预测中空保温玻璃的销售预测MATLAB程序如下：clearp=0.5677 -0.2435 -0.58390.2616 0.1754 1.07560.2983 -0.0679 0.73360.0965 -0.3516 0.69080.3185 0.8940 -0.02330.4722 0.9227 0.8708-0.0675 0.3464 0.29650.1267 0.1315 0.3773-0.4295 -0.5566 -0.20690.7275 0.6252 0.4419-0.1880 0.1225

17、 0.0618;t=0.3684 0.3541 0.1176 0.4274 -0.0649 0.2077 0.2175 0.2340 -0.1785 0.9792 -0.4979;（ rand(state,0); %保证每次结果都相同net=newff(minmax(p),3,7,1,tansig,tansig,purelin,trainlm);net.trainparam.lr=0.1;net.trainparam.epochs=1000;net.trainparam.goal=0.001;net=train(net,p,t);y=sim(net,p)e=t-yres=norm(e) %整个

18、网络误差m=2000:2010;plot(m,t,-+,m,y,o)xlabel(日期)ylabel(国内贷款额增长率)00.511.522.533.5410-410-310-210-11001014 EpochsTraining-Blue Goal-BlackPerformance is 0.000176118, Goal is 0.00120002001200220032004200520062007200820092010-0.500.51日 期国内贷款额增长率确定最佳隐含层神经元个数网络的设计及训练代码如下：s=4:8res=1:5for i=1:5rand(state,0);net=newff(minmax(p),s(i),1,tansig,purelin);net.trainparam.lr=0.1;net.trainparam.epochs=1000;net.trainparam.goal=0.001;net=train(net,p,t);y=sim(net,p)e=t