/*
 Fuzzy string matching (1996)
by Cristobal Brahm lbrahm@cec.uchile.cl
(theoretical method given by Maurice.Clerc@WriteMe.com)
*/

#include <conio.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <alloc.h>
#define ALPHA		1
#define BETA		-1.2
#define GAMMA		2
#define LARGO_MAX	5
#define LARGO_ALF	4

float confusion[LARGO_ALF][LARGO_ALF]={
{1,0.2,0,0},
{0.2,1,0,0},
{0,0,1,0},
{0,0,0,1}
};

char alfabeto[LARGO_ALF]={'A','B','C','D'};

int letra_numero(char letra)
{
int a;
for (a=0;a<LARGO_ALF;a++)
{
	if (alfabeto[a]==letra)
		return(a);
}
return(-1);
}


float fuzzy_law(int d)
{
return (ALPHA*exp(BETA*pow(d,GAMMA)));
}

float similarity(char a, char b)
{
if (confusion[letra_numero(a)][letra_numero(b)]>=0)
	return(confusion[letra_numero(a)][letra_numero(b)]);
else
	return(0);
}

char string1[LARGO_MAX]="ABADC";
char string2[LARGO_MAX]="BADCA";
char string3[LARGO_MAX]="BADC";
char string4[LARGO_MAX]="ABAC";
char string5[LARGO_MAX]="AADC";
char string6[LARGO_MAX]="ABAD";


calcula_fuzzy(char string[LARGO_MAX], float fuzzy_vector[LARGO_ALF*LARGO_MAX])
{
int x,y,z,d,pos;
float fl,fuzzy,s;
for (x=0;x<LARGO_ALF;x++)
{
	for (y=0;y<LARGO_MAX;y++)
	{
		fuzzy=0;
		for (z=0;z<LARGO_MAX;z++)
		{
			d=abs(y-z);
			fl=fuzzy_law(d);
			s=similarity(alfabeto[x],string[z]);
			fuzzy=fuzzy+fl*s;
		}
		fuzzy_vector[LARGO_MAX*x+y]=min(fuzzy,1);
	}
}
}

float fuzzy_cos(float f1[LARGO_ALF*LARGO_MAX], float f2[LARGO_ALF*LARGO_MAX])
{
int x;
float punto,mod1,mod2,coseno;
punto=0;
mod1=0;
mod2=0;
for (x=0;x<LARGO_ALF*LARGO_MAX;x++)
{
	punto=punto+f1[x]*f2[x];
	mod1=mod1+f1[x]*f1[x];
	mod2=mod2+f2[x]*f2[x];
}
mod1=sqrt(mod1);
mod2=sqrt(mod2);
if (mod1*mod2!=0)
	coseno=punto/(mod1*mod2);
else
	coseno=0;
return(coseno);
}

main()
{
int x;
float fv1[LARGO_ALF*LARGO_MAX];
float fv2[LARGO_ALF*LARGO_MAX];
float fv3[LARGO_ALF*LARGO_MAX];
float fv4[LARGO_ALF*LARGO_MAX];
float fv5[LARGO_ALF*LARGO_MAX];
float fv6[LARGO_ALF*LARGO_MAX];
float coseno;
clrscr();
calcula_fuzzy(string1,fv1);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv1[x]);
printf("\n\n");
getch();
calcula_fuzzy(string2,fv2);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv2[x]);
printf("\n\n");
getch();
calcula_fuzzy(string3,fv3);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv3[x]);
printf("\n\n");
getch();
calcula_fuzzy(string4,fv4);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv4[x]);
printf("\n\n");
getch();
calcula_fuzzy(string5,fv5);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv5[x]);
printf("\n\n");
getch();
calcula_fuzzy(string6,fv6);
for(x=0;x<LARGO_ALF*LARGO_MAX;x++)
	printf(" %f ",fv6[x]);
printf("\n\n");
getch();
coseno=fuzzy_cos(fv1,fv1);
printf("Coseno(1,1) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv1,fv2);
printf("Coseno(1,2) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv1,fv3);
printf("Coseno(1,3) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv1,fv4);
printf("Coseno(1,4) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv1,fv5);
printf("Coseno(1,5) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv1,fv6);
printf("Coseno(1,6) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv2,fv2);
printf("Coseno(2,2) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv2,fv3);
printf("Coseno(2,3) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv2,fv4);
printf("Coseno(2,4) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv2,fv5);
printf("Coseno(2,5) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv2,fv6);
printf("Coseno(2,6) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv3,fv3);
printf("Coseno(3,3) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv3,fv4);
printf("Coseno(3,4) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv3,fv5);
printf("Coseno(3,5) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv3,fv6);
printf("Coseno(3,6) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv4,fv4);
printf("Coseno(4,4) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv4,fv5);
printf("Coseno(4,5) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv4,fv6);
printf("Coseno(4,6) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv5,fv5);
printf("Coseno(5,5) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv5,fv6);
printf("Coseno(5,6) : %f\n",coseno);
getch();
coseno=fuzzy_cos(fv6,fv6);
printf("Coseno(6,6) : %f\n",coseno);
getch();
}