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;}////////////////////////////////////////////////////////// |
