?? somnet.cpp
字號:
// Somnet.cpp: implementation of the CSomnet class.
//
//////////////////////////////////////////////////////////////////////
#include "stdafx.h"
#include "Somnet.h"
#include <fstream.h>
#include <iomanip.h>
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////
CSomnet::CSomnet()
{
current_iteration_num=0;
min_distance=1000.0;
study_efficiency=efficiency;
neighbor_width=OutputLayerColumn;
neighbor_height=OutputLayerRow;
winner_row=OutputLayerRow/2;
winner_column=OutputLayerColumn/2;
}
CSomnet::~CSomnet()
{
}
int CSomnet::getcurrentnum()
{
return current_iteration_num;
}
long double CSomnet::getmindistance()
{
return min_distance;
}
/****************************************************************/
//該函數(shù)初始化輸入數(shù)組
/****************************************************************/
CSomnet::CSomnet(int layer,int nd)
{
int i;
layernumber=layer;
/*midianumber=midia;*/
ndnumber=nd;
inputMode=new double*[layernumber];
for(i=0;i<layernumber;i++)
{
*(inputMode+i)=new double[ndnumber];
}
// for(i=0;i<layernumber;i++)
// for(j=0;j<midianumber;j++)
// {
// *(*(inputMode+i)+j)=new int[ndnumber];
// }
InputLayerNum=ndnumber;
weight=new double*[OutputLayerColumn*OutputLayerRow];
for(i=0;i<OutputLayerRow*OutputLayerColumn;i++)
{
weight[i]=new double[InputLayerNum];
}
}
/****************************************************************/
//該函數(shù)初始化距離變量為0,初始化保存勝出節(jié)點的位置的變量
/****************************************************************/
void CSomnet::init_distance()
{
for(int i=0;i<OutputLayerRow;i++)
for(int j=0;j<OutputLayerColumn;j++)
distance[i*OutputLayerColumn+j]=0.0;
}
/****************************************************************/
//該函數(shù)用于計算歐氏距離,并找到獲勝神經(jīng)元
/****************************************************************/
void CSomnet::eula_distance(int z)
{
int i,j,n;
int ttLow,ttUp,ppLow,ppUp;
ttLow=winner_column-neighbor_width/2;
ttUp=winner_column+neighbor_width/2;
ppLow=winner_row-neighbor_height/2;
ppUp=winner_row+neighbor_height/2;
if(ttLow<0)
ttLow=0;
if(ttUp>=OutputLayerColumn)
ttUp=OutputLayerColumn-1;
if(ppLow<0)
ppLow=0;
if(ppUp>=OutputLayerRow)
ppUp=OutputLayerRow-1;
for( i=ppLow;i<=ppUp;i++)
for( j=ttLow;j<=ttUp;j++)
{
if(!(flag[i][j]==100)) //用一個變化的外部變量控制循環(huán)的方法
{
/*for(m=0;m<midianumber;m++)*/
for(n=0;n<ndnumber;n++)
distance[i*OutputLayerColumn+j]+=
pow((inputMode[z][n]-weight[i*OutputLayerColumn+j][n]),2);
if(distance[i*OutputLayerColumn+j]<min_distance)
{
min_distance=distance[i*OutputLayerColumn+j];
temp_row=i;
temp_column=j;
}
}
}
if(current_iteration_num>0)
{
if(min_distance<=error_limit)
{
row[temp_row]=temp_row;
column[temp_column]=temp_column;
flag[temp_row][temp_column]=100;
}
}
}
/****************************************************************/
//調(diào)整權(quán)值
/****************************************************************/
void CSomnet::weight_change(int z)
{
int i,j,n;
int ttLow,ttUp,ppLow,ppUp;
winner_row=temp_row;
winner_column=temp_column;
ttLow=winner_column-neighbor_width/2;
ttUp=winner_column+neighbor_width/2;
ppLow=winner_row-neighbor_height/2;
ppUp=winner_row+neighbor_height/2;
if(ttLow<0)
ttLow=0;
if(ttUp>=OutputLayerColumn)
ttUp=OutputLayerColumn-1;
if(ppLow<0)
ppLow=0;
if(ppUp>=OutputLayerRow)
ppUp=OutputLayerRow-1;
for(i=ppLow;i<=ppUp;i++)
for(j=ttLow;j<=ttUp;j++)
{
if(!(flag[i][j]==100))
{
/* for(m=0;m<midianumber;m++)*/
for(n=0;n<ndnumber;n++)
weight[i*OutputLayerColumn+j][n]=
weight[i*OutputLayerColumn+j][n]+
study_efficiency*(inputMode[z][n]-weight[i*OutputLayerColumn+j][n]);
}
}
}
/****************************************************************/
//調(diào)整學習效率以及獲勝節(jié)點的鄰域大小
/****************************************************************/
void CSomnet::paraChange()
{
study_efficiency=
study_efficiency*(1.0-((double)current_iteration_num)/total_iteration_Num);
neighbor_width=
int(neighbor_width*(1.0-((double)current_iteration_num)/total_iteration_Num));
neighbor_height=
int(neighbor_height*(1.0-((double)current_iteration_num)/total_iteration_Num));
}
/*****************************************************************/
//該函數(shù)用于將所有輸入模式從文件中讀入,并存放到數(shù)組inputMode中;
//同時進行權(quán)值的初始化,采用隨機賦值的方法;
/*****************************************************************/
void CSomnet::initialize()
{
int i,j,k;
for(i=0;i<OutputLayerRow;i++)
row[i]=100;
for(j=0;j<OutputLayerColumn;j++)
column[j]=100;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
flag[i][j]=0;
//從文件中將所有輸入模式讀入,并存放到數(shù)組inputMode中
FILE *pf=fopen("相關(guān)數(shù)據(jù)\\輸入數(shù)據(jù)\\input.txt","a+");
if(pf==NULL)
{
cout<<"Can not open input file!\n";
exit(0);
}
for(i=0;i<layernumber;i++)
/* for(j=0;j<midianumber;j++)*/
for(k=0;k<ndnumber;k++)
fscanf(pf,"%lf",&inputMode[i][k]);
/////////////////////////////////////////////////////
//用于測試是否能夠正確讀入輸入模式
char character[26];
for(i=0;i<26;i++)
character[i]=(65+i);
ofstream mode ("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\向量模式.txt",ios::out) ;
for(i=0;i<layernumber;i++)
{
mode<<character[i]<<"\n"<<endl;
/* for(j=0;j<midianumber;j++)*/
for(k=0;k<ndnumber;k++)
mode<<inputMode[i][k]<<" ";
mode<<"\n";
mode<<"\n\n\n";
}
////////////////////////////////////////////////////
//權(quán)值初始化,采用隨機賦值的方法
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
for(k=0;k<InputLayerNum;k++)
weight[i*OutputLayerColumn+j][k]=(double(rand()%101))/100.0;
/////////////////////////////////////////////////////
//用于測試是否能夠正確初始化權(quán)值
ofstream quan ("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\初始的權(quán)值.txt",ios::out) ;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
{
quan<<"\n\n\n"<<"Node["<<i+1<<"]["<<j+1<<"]"<<"\n";
for(k=0;k<InputLayerNum;k++)
{
if(k%5==0)
quan<<"\n";
quan<<setprecision(6)<<setiosflags(ios::fixed)<<weight[i*OutputLayerColumn+j][k]<<" ";
}
quan<<"\n\n\n";
}
////////////////////////////////////////////////////
}
/*****************************************************************/
//對網(wǎng)絡進行訓練
/*****************************************************************/
void CSomnet::Treatnet()
{
int l,i,j,k;
int iteration_numbers[26];
int total_num=0;
char character[26];
for(l=0;l<layernumber;l++)
{
iteration_numbers[l]=0;
character[l]=(65+l); //65的ASC碼為'a',character數(shù)組中儲存了26個英文字母
}
initialize();
for(l=0;l<layernumber;l++) //26是指26個英文字母
{
winner_row=OutputLayerRow/2;
winner_column=OutputLayerColumn/2;
while(current_iteration_num<total_iteration_Num)
{
init_distance();
eula_distance(l);
weight_change(l);
if(min_distance<=error_limit)
break;
++current_iteration_num;
paraChange();
};
iteration_numbers[l]=current_iteration_num+1;
neighbor_width=OutputLayerColumn;
neighbor_height=OutputLayerRow;
study_efficiency=efficiency;
current_iteration_num=0;
min_distance=1000.0;
}
for(l=0;l<layernumber;l++)
total_num+=iteration_numbers[l];
ofstream iteration_num("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\迭代次數(shù).txt",ios::out);
for(l=0;l<layernumber;l++)
{
iteration_num<<character[l]<<" 迭代"<<iteration_numbers[l]<<"次!\n"<<endl;
if(l==(layernumber-1))
iteration_num<<"整個訓練過程共迭代"<<total_num<<"次!\n"<<endl;
}
ofstream all_weight("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\訓練后所有權(quán)值.txt",ios::out ) ;
ofstream winner_weight("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\訓練后勝出權(quán)值.txt",ios::out) ;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
{
printf("\n\n\n");
all_weight<<"\n\n\n"<<"Node["<<i+1<<"]["<<j+1<<"]"<<"\n";
for(k=0;k<InputLayerNum;k++)
{
if(k%5==0)
{
printf("\n");
all_weight<<"\n";
}
/* /////////////////////////////////////////////////
if(weight[i*OutputLayerColumn+j][k]>0.9999999)
weight[i*OutputLayerColumn+j][k]=1.0;
if(weight[i*OutputLayerColumn+j][k]<0.0000001)
weight[i*OutputLayerColumn+j][k]=0.0;
*/ /////////////////////////////////////////////////
printf("%f ",weight[i*OutputLayerColumn+j][k]);
all_weight<<setprecision(8)<<setiosflags(ios::fixed)<<weight[i*OutputLayerColumn+j][k]<<" ";
}
}
ofstream winner_node("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\獲勝節(jié)點.txt",ios::out) ;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
{
if(flag[i][j]==100)
{
printf("\n\n\n");
winner_weight<<"\n\n\n"<<"Node["<<i+1<<"]["<<j+1<<"]"<<"\n";
for(k=0;k<InputLayerNum;k++)
{
if(k%5==0)
{
printf("\n");
winner_weight<<"\n";
}
/* /////////////////////////////////////////////////
if(weight[i*OutputLayerColumn+j][k]>0.9999999)
weight[i*OutputLayerColumn+j][k]=1.0;
if(weight[i*OutputLayerColumn+j][k]<0.0000001)
weight[i*OutputLayerColumn+j][k]=0.0;
*/ /////////////////////////////////////////////////
printf("%f ",weight[i*OutputLayerColumn+j][k]);
winner_weight<<setprecision(8)<<setiosflags(ios::fixed)<<weight[i*OutputLayerColumn+j][k]<<" ";
}
winner_node<<"Node["<<i+1<<"]["<<j+1<<"]"<<endl;
}
}
printf("\n");
}
/************************************************************/
//利用標準數(shù)據(jù)測試訓練后的網(wǎng)絡
/************************************************************/
void CSomnet::test_netWork_1()
{
int i,j,l,n;
ofstream test1("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\標準測試.txt",ios::out);
char character[26];
for(i=0;i<26;i++)
character[i]=(65+i);
for(l=0;l<layernumber;l++)
{
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
distance[i*OutputLayerColumn+j]=0.0;
min_distance=1000;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
{
/* for(m=0;m<midianumber;m++)*/
for(n=0;n<ndnumber;n++)
distance[i*OutputLayerColumn+j]+=(long double)
pow(((long double)inputMode[l][n]-(long double)weight[i*OutputLayerColumn+j][n]),2);
if(distance[i*OutputLayerColumn+j]<min_distance)
{
min_distance=distance[i*OutputLayerColumn+j];
temp_row=i;
temp_column=j;
}
}
test1<<character[l]<<"'s winner is Node["<<temp_row+1<<"]["<<temp_column+1<<"]"<<endl<<endl;
}
}
/************************************************************/
//利用非標準數(shù)據(jù)測試訓練后的網(wǎng)絡
/************************************************************/
void CSomnet::test_netWork_2()
{
int i,j,k,l,n;
ofstream test2("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\非標準測試.txt",ios::out);
char character[26];
if(inputMode!=NULL)
delete []inputMode;
inputMode=new double*[layernumber];
for(i=0;i<layernumber;i++)
*(inputMode+i)=new double[ndnumber];
if(inputMode!=NULL)
{
FILE *pf1=fopen("相關(guān)數(shù)據(jù)\\輸入數(shù)據(jù)\\非標準數(shù)據(jù)測試.txt","a+");
if(pf1==NULL)
{
cout<<"Can not open input file!\n";
exit(0);
}
for(i=0;i<layernumber;i++)
/* for(j=0;j<midianumber;j++)*/
for(k=0;k<ndnumber;k++)
fscanf(pf1,"%lf",&inputMode[i][k]);
for(i=0;i<26;i++)
character[i]=(65+i);
// char character[26];
// for(i=0;i<26;i++)
// character[i]=(65+i);
ofstream mode1 ("相關(guān)數(shù)據(jù)\\輸出數(shù)據(jù)\\非標準向量模式.txt",ios::out) ;
for(i=0;i<layernumber;i++)
{
mode1<<character[i]<<"\n"<<endl;
/* for(j=0;j<midianumber;j++)*/
for(k=0;k<ndnumber;k++)
mode1<<inputMode[i][k]<<" ";
mode1<<"\n";
mode1<<"\n\n\n";
}
for(l=0;l<layernumber;l++)
{
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
distance[i*OutputLayerColumn+j]=0.0;
min_distance=1000;
for(i=0;i<OutputLayerRow;i++)
for(j=0;j<OutputLayerColumn;j++)
{
/* for(m=0;m<midianumber;m++)*/
for(n=0;n<ndnumber;n++)
distance[i*OutputLayerColumn+j]+=(long double)
pow(((long double)inputMode[l][n]-(long double)weight[i*OutputLayerColumn+j][n]),2);
if(distance[i*OutputLayerColumn+j]<min_distance)
{
min_distance=distance[i*OutputLayerColumn+j];
temp_row=i;
temp_column=j;
}
}
test2<<character[l]<<"'s winner is Node["<<temp_row+1<<"]["<<temp_column+1<<"]"<<endl<<endl;
}
}
}
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -