//?BPnetwork.cpp?:?定義控制臺應用程序的入口點。
//

#include?“stdafx.h“
#include?
#include?
#include?
#include?
#include?
#define?Data?75?//樣本容量
#define?In?4?//樣本參數個數
#define?Out?3??//輸出參數個數
#define?Neuron?10?//神經元個數
#define?TrainC?2000?//訓練次數
#define?StudyRate_A?0.2//輸出學習率?基于經驗
#define?StudyRate_B?0.4//輸出學習率?基于事實
#define?StudyRate_a?0.2//輸入學習率?基于經驗
#define?StudyRate_b?0.3//輸入學習率?基于事實
double?d_in[Data][In];?//樣本輸入
double?d_out[Out][Data];?//樣本輸出
double?w[Neuron][In];?//輸入對神經元的權
double?v[Out][Neuron];?//神經元對輸出的權
double?dw[Neuron][In];??//w的修正
double?dv[Out][Neuron];?//v的修正
double?o[Neuron];?//神經元通過激活函數的輸出
double?OutputData[Out];?//BP網絡的輸出
double?Maxin[In]Minin[In]Maxout[Out]Minout[Out];

//void?ReadData（）
//{
//?for?（int?i?=?0;?i?//?{
//??d_in[i][0]=rand（）%100;
//??d_in[i][1]=rand（）%100;
//??d_out[0][i]=d_in[i][0]+d_in[i][1];
//??//printf（“data%d:?%lf?+?%lf?=?%lf\n“i+1d_in[i][0]d_in[i][1]d_out[0][i]）;
//?}
//}

void?ReadData（）
{
	int?type=0;
	FILE*??file;
	errno_t?err=fopen_s（&file“Iris-train.txt““r“）;
	if（err==0）
		for（int?i=0;i		{
			if（feof（file））
				break;
			else
			{
				fscanf_s（file“%lf?%lf?%lf?%lf?%d?“&d_in[i][0]&d_in[i][1]&d_in[i][2]&d_in[i][3]&type）;
				d_out[0][i]=0;
				d_out[1][i]=0;
				d_out[2][i]=0;
				d_out[type][i]=1;
			}
		}
	else
		printf_s（“%d“err）;
}

void?InitBPNework（）
{
?int?ij;
?
?for?（i?=?0;?i??{
??Minin[i]=Maxin[i]=d_in[0][i];
??for?（j?=?0;?j???{
???Maxin[i]=Maxin[i]>d_in[j][i]?Maxin[i]:d_in[j][i];
???Minin[i]=Minin[i]??}
?}
?for?（i?=?0;?i??{
??Minout[i]=Maxout[i]=d_out[0][i];
??for?（j?=?0;?j???{
???Maxout[i]=Maxout[i]>d_out[i][j]?Maxout[i]:d_out[i][j];
???Minout[i]=Minout[i]??}
?}
?//printf（“max?of?in:%lf%lf\nmin?of?in:%lf%lf\n“Maxin[0]Maxin[1]Minin[0]Minin[1]）;
?//printf（“max?of?out:%lf\nmin?of?out:%lf\n“Maxout[0]Minout[0]）;
?
?for?（i?=?0;?i???for（j?=?0;?j????d_in[j][i]=（d_in[j][i]-Minin[i]+1）/（Maxin[i]-Minin[i]+1）;
?for?（i?=?0;?i???for（j?=?0;?j????d_out[i][j]=（d_out[i][j]-Minout[i]+1）/（Maxout[i]-Minout[i]+1）;
?
?for?（i?=?0;?i???for?（j?=?0;?j???{???
???w[i][j]=（rand（）%2-1）;
???dw[i][j]=0;
??}
?for?（i?=?0;?i???for?（j?=?0;?j???{
???v[j][i]=（rand（）%2-1）;
???dv[j][i]=0;
??}
}

void?ComputO（int?var）
{?
?int?ij;
?double?sumy;
????
?for?（i?=?0;?i??{
??sum=0;
??for?（j?=?0;?j???{
???sum+=w[i][j]*d_in[var][j];
???o[i]=1/（1+exp（-1*sum））;
??}
?}
?
????for?（i?=?0;?i??{
??sum=0;
??for?（j?=?0;?j????sum+=v[i][j]*o[j];
??OutputData[i]=sum;
??//printf（“data%d:output?%lf?ans?%lf\n“varOutputData[i]d_out[i][var]）;
?}???
}

void?BackUpdate（int?va

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1424??2015-10-23?14:29??Iris-test.txt

?????文件???????1425??2015-10-23?14:29??Iris-train.txt

?????文件???????6092??2015-10-31?17:37??BPnetwork.cpp

-----------?---------??----------?-----??----

?????????????????8941????????????????????3