/*神经元模型*/
#include < math.h>
#include < stdio.h>
#include < time.h>
#include < stdlib.h>
#include < conio.h>
#define AND1 .571388
#define AND2 .560999
#define OR1 1.149045
#define OR2 1.131896
#define SITA 1
#define NIU .1
#define OMIGA 50
float w[2];
float terg[4][3];
main() {
int m_result = 1, s;
float w1[2];
while ( m_result != 5 ) {
m_result = menu ( m_result );
date();
switch ( m_result ) {
case ( 1 ) : /* 学习*/
case ( 3 ) :
randomize(); /*重新学习*/
w[0] = random ( 80 ) +10;
w[0] = w[0]/100.0;
w[1] = random ( 80 ) +10;
w[1] = w[1]/100.0;
study();
break;
case ( 2 ) :
practic ( w ); /*实践新模型*/
break;
case ( 4 ) :
printf ( " what do you want to see, OR or AND?(1/2) " ); if ( s==1 ) {w1[0] = OR1; w1[1] = OR2; } else if ( s==2 ) {w1[0] = AND1; w1[1] = AND2; } else break;
practic ( w1 );
break; /*演示已有模型*/
default: {}
} /*End Switch*/
} /*End While */
clrscr();
}
menu ( int r ) {
int r1;
clrscr();
printf ( " \t\t1\tStudy a Model\n" ); /* 学习 */ printf ( " \t\t2\tPratice the model\n" ); /*重新学习 */ printf ( " \t\t3\tDefail the MOdel\n" ); /*实践新模型 */ printf ( " \t\t4\tPrepared Model\n" ); /*演示已有模型*/ printf ( " \t\t5\tEnd\n" ); /*结束 */ printf ( " \n\n\t\tWhat do you want to choice?(%d)" , r
); r1 = getche();
if ( r1!=13 ) r = r1-48;
return ( r );
}
date() { /*显示模型*/
int i, j;
clrscr();
for ( i = 0; i < 4; i++ )
for ( j = 0; j< 3; j++ ) {
gotoxy ( 20+j*4, i+5 );
printf ( "%.0f", terg
[i
][j
] ); }
printf ( "\n\n\tThe quarities:%f, %f \n", w
[0], w
[1] ); }
study() { /*学习子函数*/
float x, y, width;
int wrong = 1, i, total = 0;
input(); /*andinput(); */ /*orinput(); */
while ( wrong ) {
wrong = 0;
for ( i = 0; i< 4; i++ ) {
x = OMIGA* ( w[0]*terg[i][0]+w[1]*terg[i][1]-SITA );
y
= 1.0/ ( 1.0+exp ( -x
) ); width = terg[i][2]-y; /*误差 = "教师信号-Y" */
if ( fabs ( width
) >=.01 ) { wrong = 1;
total++;
w[0] = w[0]* ( 1+width*NIU ); /*改变权值0*/
w[1] = w[1]* ( 1+width*NIU ); /*改变权值1*/
if ( total> =200 ) /*防止程序进入死循环*/
{error(); return 1; }
}
}
}/*若i的四次循环均满足" fabs(width)> =.01" 则wrong=0, 结束循环*/
}
practic ( quarity ) /*实践新模型或演示已有模型, 由quarity[]值决定*/
float quarity[2];
{
int k, y;
float x, width, x1, x2;
for ( k = 0; k < 4; k++ ) {
gotoxy ( 0, 4+k*3 );
printf ( "\nPlease input :x1, x2 " ); scanf ( "%f, %f", &x1
, &x2
); x = OMIGA* ( quarity[0]*x1+quarity[1]*x2-SITA );
y
= 1.0/ ( 1.0+exp ( -x
) ) +.01; /*精度为0.01*/ printf ( "\t\t\ty = "%d
" , y); getch();
}
}
input() /*当选择"学习"功能时输入事件*/
{int i;
printf("Please input the 4 groups data:\n\n");
for(i = 0; i< 4; i++){
printf("\tThe %d group x1, x2, teacher's point: ", i+1);
scanf("%f, %f, %f", &terg[i][0], &terg[i][1], &terg[i][2]);
}
clrscr();
}
andinput()
{terg[0][0] = 1; terg[0][1] = 1; terg[0][2] = 1;
terg[1][0] = 1;
terg[1][1] = 0;
terg[1][2] = 0;
terg[2][0] = 0;
terg[2][1] = 1;
terg[2][2] = 0;
terg[3][0] = 0;
terg[3][1] = 0;
terg[3][2] = 0;
}
orinput()
{terg[0][0] = 1;
terg[0][1] = 1;
terg[0][2] = 1;
terg[1][0] = 1;
terg[1][1] = 0;
terg[1][2] = 1;
terg[2][0] = 0;
terg[2][1] = 1;
terg[2][2] = 1;
terg[3][0] = 0;
terg[3][1] = 0;
terg[3][2] = 0;
}
error() /*初始权值选择不当或模型无法实现(如异或模型)时退出*/
{clrscr();
printf("\n\n\n\tBecause the firse QUARITIES's choise is unavailble, \n");
printf("\tIt can't get the right result, forveve!");
}