人工神经网络1..docVIP

人工神经网络——BP网络 BP神经网络逼近函数 步骤1：假设频率参数k=1，绘制要逼近的非线性函数的曲线。函数的曲线如图1所示 k=1; p=[-1:.05:8]; t=1+sin(k*pi/4*p); plot(p,t) title(要逼近的非线性函数); xlabel(时间); ylabel(非线性函数); 图1要逼近的非线性函数曲线 步骤2：网络的建立 应用newff()函数建立BP网络结构。隐层神经元数目n可以改变，暂设为n=3，输出层有一个神经元。选择隐层和输出层神经元传递函数分别为tansig函数和purelin函数，网络训练的算法采用Levenberg – Marquardt算法trainlm。 n=3; net = newff(minmax(p),[n,1],{tansig purelin},trainlm); 对于初始网络，可以应用sim（）函数观察网络输出 y1=sim(net,p); figure; plot(p,t,-,p,y1,:) title(未训练网络的输出结果); xlabel(时间); ylabel(仿真输出--原函数-); 同时绘制网络输出曲线，并与原函数相比较，结果如图2所示。 图2未训练网络的输出结果 其中 “ ” 代表要逼近的非线性函数曲线； “‥‥‥” 代表未经训练的函数曲线； 因为使用newff( )函数建立函数网络时，权值和阈值的初始化是随机的，所以网络输出结构很差，根本达不到函数逼近的目的，每次运行的结果也有时不同。 步骤3：网络训练 应用train（）函数对网络进行训练之前，需要预先设置网络训练参数。将训练时间设置为50，训练精度设置为0.01，其余参数使用缺省值。训练后得到的误差变化过程如图3所示。 图3 net.trainParam.epochs=50; %网络训练步骤设置为50 net.trainParam.goal=0.01;%网络训练精度设置为0.01 net=train(net,p,t);%开始训练网络 TRAINLM, Epoch 0/50, MSE 2.26463/0.01, Gradient 364.666/1e-010 TRAINLM, Epoch 9/50, MSE 0.0096751/0.01, Gradient 16.4769/1e-010 TRAINLM, Performance goal met. 从以上结果可以看出，网络训练速度很快，经过一次循环跌送过程就达到了要求的精度0.01。 步骤4： 网络测试 对于训练好的网络进行仿真： y2=sim(net,p); figure; plot(p,t,-,p,y1,:,p,y2, --) title(训练后网络的输出结果); xlabel(时间); ylabel(仿真输出); 绘制网络输出曲线，并与原始非线性函数曲线以及未训练网络的输出结果曲线相比较， 比较出来的结果如图4所示。 图4训练后网络的输出结果 其中 “ ” 代表要逼近的非线性函数曲线； “‥‥‥” 代表未经训练的函数曲线； “―――” 代表经过训练的函数曲线； 从图中可以看出，得到的曲线和原始的非线性函数曲线很接近。这说明经过训练后，BP网络对非线性函数的逼近效果比较好。 不同频率下的逼近效果 改变非线性函数的频率和BP函数隐层神经元的数目，对于函数逼近的效果有一定的影响。网络非线性程度越高，对于BP网络的要求越高，则相同的网络逼近效果要差一些；隐层神经元的数目对于网络逼近效果也有一定影响，一般来说隐层神经元数目越多，则BP网络逼近非线性函数的能力越强。 下面通过改变频率参数和非线性函数的隐层神经元数目来加以比较证明。 频率参数设为k=1，当隐层神经元数目分别取n=3、n=6、n=8，时得到了训练后的网络输出结果如图5、6、7所示。 图5 当n=3时训练后网络的输出结果 图6 当n=6时训练后网络的输出结果 图7当n=8时训练后网络的输出结果 其中 “ ” 代表要逼近的非线性函数曲线； “‥‥‥” 代表未经训练的函数曲线； “―――” 代表经过训练的函数曲线； （2）频率参数设为k=2，当隐层神经元数目分别取n=3、n=6、n=8时，得到了训练后的网络输出结果如图8、9、10所示。 图8当n=3时训练后网络的输出结果 图9当n=6时训练后网络的输出结果 图10当n=8时训练后网络的输出结果 其中 “ ” 代表要逼近的非线性函数曲线；