BP网络练习

利用BP算法及Sigmoid函数，研究函数的逼近问题：

f(x)=1/x,  1<=x<=100;

BP神经网络：

clc;close all;clear;
N = 7;
P = rand(1,N);
P = P*100
T = P.^(-1)
w1 = rand(1,2);
w2 = rand(2,1);
E = 0;
Eta = 0.5;
epoch = 0;
while (epoch==0||0.5*(1/N)*E>0.00001)  
    E = 0;
    for n=1:N
            v1 = (w1')*P(n);                 
            y1 = (exp(-1*v1)+1).^(-1);
            v2 = (w2')*y1;
            O = (exp(-1*v2)+1).^(-1);
            e=T(n)-O;
            E=E+0.5*(e.^2);
            Delta2 = e*O*(1-O);
            w2 = w2+Eta*Delta2*y1;
            Delta1(1,1) = y1(1,1)*(1-y1(1,1))*([1 1]*(Delta2*w2));
            Delta1(2,1) = y1(2,1)*(1-y1(2,1))*([1 1]*(Delta2*w2));
            w1(1,1) = w1(1,1)+Eta*Delta1(1,1)*P(n);
            w1(1,2) = w1(1,2)+Eta*Delta1(2,1)*P(n);