BP神经网络预测的matlab代码.doc

附录5：BP神经网络预测的matlab代码：P= 0 0.1386 0.2197 0.2773 0.3219 0.3584 0.3892 0.4159 0.4394 0.4605 0.4796 0.4970 0.5278 0.5545 0.5991 0.6089 0.6182 0.6271 0.6356 0.6438 0.6516 0.6592 0.6664 0.6735 0.7222 0.7275 0.7327 0.7378 0.7427 0.7475 0.7522 0.7568 0.7613 0.76570.7700T=0.44550.3230.41160.32550.44860.29990.49260.22490.48930.23570.48660.22490.48190.22170.49970.22690.50270.2170.51550.19180.50580.23950.45410.24080.40540.27010.39420.33160.21970.29630.55760.10610.49560.2670.51260.22380.53140.20830.51910.2080.51330.18480.50890.2420.48120.21290.49270.2870.48320.27420.59690.24030.50560.21730.53640.19940.52780.20150.51640.22390.44890.24040.48690.29630.48980.19870.50750.29170.49430.2902threshold=0 1net=newff(threshold,11,2,tansig,logsig,trainlm);net.trainParam.epochs=6000net.trainParam.goal=0.01LP.lr=0.1;net=train(net,P,T)P_test= 0.7742 0.77840.78240.78640.79020.7941 out=sim(net,P_test)友情提示：以上面0.7742为例0.7742=ln(47+1)/5因为网络输入有一个元素，对应的是测试时间，所以P只有一列，Pi=log(t+1)/10，这样做的目的是使得这些数据的范围处在0 1区间之内，但是事实上对于logsin命令而言输入参数是正负区间的任意值，而将输出值限定于0到1之间。可能只有当输入值在0到1之间，输出的函数才可以通过返归依化得到正确的对应的答案，因此在后面有threshold来限制输入值的区间n=-5:0.1:5;a=logsig(n);plot(n,logsig(n)图像如右图所示T是指CD4Count*0.02，与HIV分别对应的值，是输出值，有两组数据，T的第一列数据= 表格当中（CD4Count*0.02）0.02，再取对数，再10T的第一列数据= 表格当中（HIV）的对数，再5这样就完成了最开始的数据导入工作，即为了训练神经而进行的输入与输出。threshold=0 1- 变量threshold用于规定输入向量的最大值和最小值，最大值为1，最小值为0net=newff(threshold,11,2,tansig,logsig,trainlm);网络中间层神经元函数采用S型正切函数tansig，设定网络的训练函数为trainlm，它采用Levenberg-Marquardt算法进行网络学习。输出层神经元函数采用S型对数函数logsignet.trainParam.epochs=6000net.trainParam.goal=0.01LP.lr=0.1;经多次尝试，当模型的训