Matlab 版 (精华区)

发信人: zjliu (秋天的萝卜), 信区: Matlab
标  题: bp网络算法对y=0.4sin(2*PI*x)+0.5的拟合
发信站: 哈工大紫丁香 (Wed May  5 20:25:09 2004), 站内信件

标  题: bp网络算法对y=0.4sin(2*PI*x)+0.5的拟合
发信站: 兵马俑BBS (Wed May  5 13:09:04 2004), 本站(bbs.xjtu.edu.cn)


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;

}

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

--
╔═══════════════════╗
║★★★★★友谊第一  比赛第二★★★★★║
╚═══════════════════╝

※ 来源:．哈工大紫丁香 bbs.hit.edu.cn [FROM: 202.118.229.162]