#include“bp_rbf.h“
double?randnumber（）//-0.1~0.1的隨機函數?
{
	return?（sin（rand（））/10）?;
}
//求某個節點的輸出值
double?f（int?nodeint?layer）{
	double?sum=0;
	for?（int?node2=0;node2yer-1];node2++）
	{
		sum=sum+w[layer][node][node2]*u[layer-1][node2];
	}
	sum=1.0/（1+exp（0-sum））;
	return?sum;
}
//求誤差函數對u的導數
double?gf（int?layerint?node）{
	double?sum=0;
	for?（int?m=0;myer+1];m++）
		sum=sum+e[layer+1][m]*u[layer+1][m]*（1-u[layer+1][m]）*w[layer+1][m][node];
	return?sum;
}

//求前饋
void?feed_forward（）//前饋輸出值
{
	for?（int?layer=1;layer<=2;layer++）
		for（int?node=0;nodeyer];node++）
			u[layer][node]=f（nodelayer）;
}
//求梯度
void?compute_gradient（）
{
	for?（int?layer=2;layer>0;layer--）
	{
		for（int?node=0;nodeyer];node++）
		{
			if?（layer==2）?e[2][node]=u[2][node]-output[node];
			else?e[layer][node]=gf（layernode）;
		}
		
		for?（int?j=0;jyer];j++）
			for（int?i=0;iyer-1];i++）
			{
????????????????g[layer][j][i]=e[layer][j]*u[layer][j]*（1.0-u[layer][j]）*u[layer-1][i];??
				
			}
	}
}
//修改權值
void?update_weights（）
{
	for?（int?layer=1;layer<=2;layer++）
		for?（int?j=0;jyer];j++）
			for（int?i=0;iyer-1];i++）
			{
				temp_w_1[layer][j][i]=w[layer][j][i];
				w[layer][j][i]=w[layer][j][i]-STUDY_SPEED*g[layer][j][i]+ALPHA*delta_w[layer][j][i];
				temp_w_2[layer][j][i]=w[layer][j][i];
				delta_w[layer][j][i]=temp_w_2[layer][j][i]-temp_w_1[layer][j][i];
			}
}

double?erro（）//誤差
{
	double?sum=0;
	sum=pow（fabs（u[2][0]-output[0]）2）+pow（fabs（u[2][1]-output[1]）2）+pow（fabs（u[2][2]-output[2]）2）;
	return?sum;
}

//////////初始化，裝載樣例
void?initial（FILE?*fp）
{
	
	if（feof（fp））?
	{
		rewind（fp）;
		k_k=（++k_k）%7;
		if（!k_k）??fprintf（fpwucha“\n“）;
		fprintf（fpwucha“%16f“wucha）;
		wucha=0;
	}
	input[0]=1.0;
	for（int?i=1;i		fscanf（fp“%f“input+i）;
	
	for（i=0;i		fscanf（fp“%d“output+i）;
	//	cout<	//////////////////////至此讀出一個樣本的數據（inputoutput）;
	for?（int?m=0;m		u[0][m]=input[m];
}

void?test（）
{
	
	FILE?*fpt=NULL;
	fpt=fopen（“test.txt““r“）;
	for?（int?i=0;i<=7;i++）
	{
		initial（fpt）;
		feed_forward（）;
		u[2][0]=floor（u[2][0]+0.5）;
		u[2][1]=floor（u[2][1]+0.5）;
		u[2][2]=floor（u[2][2]+0.5）;
		cout<		if（!（（output[0]==u[2][0]）&&（output[1]==u[2][1]）&&（output[2]==u[2][2]）））;
		else?
		{
			count++;
			//	cout<		}
	}
	fclose（fpt）;
	
}
//////////////RBF的東東///////////////////////////////////////////////////////////////////////
/////////////////////文件定位////////////////////////////////
void??fseeking（FILE*?fpint?row_no）
{??
	rewind（fp）;
	if?（!row_no）?return;
	l:while（（!feof（fp））&&（fgetc（fp）!=‘\n‘））
	{
	??file_time++;
	};
	//if（feof（fp））?return;
	row_no--;
	if?（row_no）?goto?l;?
	else	return;
}

////徑向距離/////
void?min_juli（int?*set_number）//*j返回輸入向量所在的聚類
{???
	
	//	size_t?size

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

?????文件???????4383??2007-12-28?20:51??水果智能識別系統\bisheji.dsp

?????文件????????537??2005-05-17?08:52??水果智能識別系統\bisheji.dsw

?????文件??????66560??2008-01-15?09:53??水果智能識別系統\bisheji.ncb

?????文件??????53760??2008-01-15?09:53??水果智能識別系統\bisheji.opt

?????文件????????248??2008-01-14?16:23??水果智能識別系統\bisheji.plg

?????文件???????1367??2008-01-14?16:10??水果智能識別系統\bp_rbf.h

?????文件??????20758??2008-01-14?16:20??水果智能識別系統\bp結果\權值.txt

?????文件????????734??2008-01-14?16:20??水果智能識別系統\bp結果\誤差變化序列.txt

?????文件?????278604??2008-01-14?16:10??水果智能識別系統\Debug\bisheji.exe

?????文件?????623616??2008-01-14?16:10??水果智能識別系統\Debug\bisheji.pdb

?????文件??????48618??2008-01-14?16:10??水果智能識別系統\Debug\mainApp.obj

?????文件??????????0??2008-01-14?16:10??水果智能識別系統\Debug\mainApp.sbr

?????文件???????4927??2007-12-28?12:57??水果智能識別系統\Debug\rbf.obj

?????文件?????110592??2008-01-14?16:10??水果智能識別系統\Debug\vc60.pdb

?????文件??????13200??2007-12-28?20:54??水果智能識別系統\mainApp.cpp

?????文件???????3559??2005-05-18?09:14??水果智能識別系統\matrix_op.h

?????文件?????????78??2005-05-23?09:40??水果智能識別系統\rbf.cpp

?????文件??????20319??2008-01-14?16:21??水果智能識別系統\rbf結果\權值.txt

?????文件?????306074??2008-01-14?16:21??水果智能識別系統\rbf結果\誤差.txt

?????文件??????21550??2008-01-14?16:09??水果智能識別系統\study.txt

?????文件???????7182??2008-01-14?16:10??水果智能識別系統\test.txt

?????目錄??????????0??2009-02-16?22:57??水果智能識別系統\bp結果

?????目錄??????????0??2009-02-16?22:57??水果智能識別系統\Debug

?????目錄??????????0??2009-02-16?22:57??水果智能識別系統\rbf結果

?????目錄??????????0??2009-02-16?22:57??水果智能識別系統

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

??????????????1586666????????????????????25