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