模拟.cpp

用BP算法对某个表达式进行模拟

#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; 

} 

//////////////////////////////////////////////////////////