|
- UID
- 852722
|
/*神经元模型*/
#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) " );
scanf("%d" , & s);
if(s==1){w1[0]= OR1; w1[1] = OR2; }
else if(s==2){w1[0] = AND1; w1[1] = AND2; }
elsebreak;
practic(w1);
break; /*演示已有模型*/
default:{}
} /*End Switch*/
} /*End While */
clrscr();
printf(" Good Bye!" );
}
menu(int r)
{
int r1;
clrscr();
printf(" \n\n\n" );
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 tochoice?(%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[j]);
}
printf("\n\n\tThe quarities:%f, %f \n", w[0], w[1]);
}
study() /*学习子函数*/
{
float x, y,width;
intwrong = 1, i, total = 0;
input(); /*andinput(); */ /*orinput(); */
while(wrong)
{
wrong= 0;
for(i= 0; i< 4; i++)
{
x= OMIGA*(w[0]*terg[0]+w[1]*terg[1]-SITA);
y= 1.0/(1.0+exp(-x));
width= terg[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("Pleaseinput 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[0], &terg[1], &terg[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\tBecausethe firse QUARITIES's choise is unavailble, \n");
printf("\tItcan't get the right result, forveve!");
} |
|
