B-P算法对函数y=0.4sin(2*PI*x)+0.5的拟合 #include #include #include #define M 8 //隐节点个数 #define N 20 //学习样本个数 #define PI 3.14159 ////////////////////////////////////////////////////////// static float x[N]; //x输入 static float y[N]; //y输出 static float t[N]; //t教师 static float h[M][N]; // 记录隐单元的输出 static float b[N]; //输出误差 static float b1[M][N]; //记录隐单元误差 float w1[M]={0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8}; //输入与第一隐层的权 float w2[M]={0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1}; //第一隐层与输出层的权 float q[M]={0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8}; //第一隐层的阈值 float p=0.3; //输出层的阈值 float r=0.2; //步长 int n=0; //记录训练次数 ////////////////////////////////////////////////////////////////////////// float f(float u) { return 1/(1+exp(-u)); } /////////////////////////////////////////////////////////////////////////// //赋值函数 void fuzhi(void) { int k; for(k=0;k x[k]=(float)k/20; // 输入值 for(k=0;k { t[k]=0.4*sin(2*PI*x[k])+0.5; //计算教师的值 cout<<"x["< } } //////////////////////////////////////////////////////////// //计算隐单元输出，输出单元 void jisuan(void) { int i,j; static float a[N]; //存放中间结果 for( i=0;i for( j=0;j h[i][j]=f(w1[i]*x[j]-q[i]); //计算隐单元输出 for(j=0;j { a[j]=0; //注意初值 for(i=0;i a[j]+=w2[i]*h[i][j]; y[j]=f(a[j]-p); } /////////////////////////////////////////////////////////// //计算误差 void wucha(void) { int i,j; for( j=0;j b[j]=y[j]*(1-y[j])*(t[j]-y[j]); //输出 for( i=0;i for(j=0;j b1[i][j]=b[j]*w2[i]*h[i][j]*(1-h[i][j]);//隐单元 } ////////////////////////////////////////////////////////////// //权的修改 void xiugai(void) { float c,d; //临时变量 int i,j; for( i=0;i { c=0; for( j=0;j c+=b[j]*h[i][j]; w2[i]+=r*c; } c=0; for(j=0;j c+=b[j]*(-1); p+=c*r; for(i=0;i { c=0; d=0; for(j=0;j { c+=b1[i][j]*x[j]; d+=b1[i][j]*(-1); } w1[i]+=r*c; q[i]+=r*d; } } //////////////////////////////////////////////////////////// //计算输出的总误差 float shuchu(void) { float E=0; for(int j=0;j E+=(t[j]-y[j])*(t[j]-y[j]); return E; } /////////////////////////////////////////////////////////// //主函数 void main() { char s; float Ez; fuzhi(); do{ ++n; jisuan(); wucha(); xiugai(); Ez=shuchu();cout< } while(Ez/2>0.01); // for(int j=0;j // cout<<"y["< cin>>s; } //////////////////////////////////////////////////////////
本文地址：http://com.8s8s.com/it/it1170.htm