{"id":1818,"date":"2014-10-24T22:39:00","date_gmt":"2014-10-25T03:39:00","guid":{"rendered":"http:\/\/linuxitos.com\/main\/sistema-de-control-takagi-sugeno-en-c\/"},"modified":"2014-10-24T22:39:00","modified_gmt":"2014-10-25T03:39:00","slug":"sistema-de-control-takagi-sugeno-en-c","status":"publish","type":"post","link":"https:\/\/linuxitos.com\/main\/sistema-de-control-takagi-sugeno-en-c\/","title":{"rendered":"Sistema de control Takagi-Sugeno en C"},"content":{"rendered":"
\n
\/*
 * Autor: A. L. I. Fernando Merino
 * Nombre del archivo: sistemaTakgi.c
 * Descripci\u00f3n: Simula un sistema difuso, usando el m\u00e9todo de Takgi-Sugeno
 * genera un archivo con los  puntos en pares (x,y) para graficar en gnuplot,
 * si desean graficar adicionalmente los puntos que est\u00e9n utilizando generen
 * un archivo extra con los puntos y grafiquenlos, para comparar resultados
 * 24-octubre-2014*\/

#include <stdlib.h>
#include <stdio.h>
#include <unistd.h>
#include <math.h>

\/*define el total de comandos que se enviar\u00e1 a gnuplot*\/
#define NUM_COMANDOS 5

void seleccionGaussiana();
double gaussiana(double c, double gama, double x);
double gaussDerecha(double c, double gama, double x);
double gaussIzquierda(double c, double gama, double x);
void graficarPuntosDesdeArchivo(char *cad, char *archivo, char *archivo1, char *archivo2);

int main(){
 printf(\"ntFunciones de Takgin-Sugenon\");
 seleccionGaussiana();
    return 0;
}

void seleccionGaussiana(){
 FILE * puntGI = fopen(\"gaussIzquierda.txt\", \"w\");
 FILE * puntGD = fopen(\"takagiSugeno.txt\", \"w\");
 FILE * archivoPuntos = fopen(\"gaussiana.txt\", \"w\");
 double x,cenGaus,anchoGaus,inicio, fin;
 double cenIzq,anchoIzq,cenDer,anchoDer;
 double ug=0,ugi=0,ugd=0,y1=0,y2=0,y3=0,y;
 \/*datos que se pueden configurar:
  * muestreo [0,20]
  * las variables a y b se pueden leer desde teclado, a tiene que ser menor a b
  * opcionalmente tambi\u00e9n se puede especificar el perido de muestreo*\/
 printf(\"Ingrese los siguientes datos:n\");
 printf(\"Inicio de la muestra: \");
 scanf(\"%lf\",&inicio);
 printf(\"Fin de la muestra: \");
 scanf(\"%lf\",&fin);
 printf(\"nConfigure la Gaussiana Izquierda:\");
 printf(\"nIngresa el centro: \");
 scanf(\"%lf\",&cenIzq);
 printf(\"Ingresa ancho a los lados del centro: \");
 scanf(\"%lf\",&anchoIzq);
 
 printf(\"nConfigure la Gaussiana central:\");
 printf(\"nIngresa el centro: \");
 scanf(\"%lf\",&cenGaus);
 printf(\"Ingresa ancho a los lados del centro: \");
 scanf(\"%lf\",&anchoGaus);
 
 printf(\"nConfigure la Gaussiana derecha:\");
 printf(\"nIngresa el centro: \");
 scanf(\"%lf\",&cenDer);
 printf(\"Ingresa ancho a los lados del centro: \");
 scanf(\"%lf\",&anchoDer);
 
 for(x=inicio;x<=fin;x+=0.01){
  \/*obtiene el valor de las membres\u00edas de cada valor de x*\/
  ugi = gaussIzquierda(cenIzq,anchoIzq,x);
  ug = gaussiana(cenGaus,anchoGaus,x);
  ugd = gaussDerecha(cenDer,anchoDer,x);
  \/*obtiene las y*\/
  y1=0.214285714*x+2.642857143;
  y2=-0.416666667*x+0.75;
  y3=0.342857143*x-1.528571429;
  \/*se susitituye en la f\u00f3rmula final*\/
  y=((ugi*y1)+(ug*y2)+(ugd*y3))\/(ugi+ug+ugd);
  \/*se guarda en un arhivo de texto con los puntos (x,y)*\/
  fprintf(puntGD, \"%lf %lf n\", x, y);
 }
 fclose(archivoPuntos);
 fclose(puntGD);
 fclose(puntGI);
 graficarPuntosDesdeArchivo(\"set title \"Graficando funciones\"\",
        \"plot \"takagiSugeno.txt\" using 1:2 with lines\",
        \"replot \"puntosAnalizar.txt\" using 1:2 with lines\",
        \"replot \"gaussIzquierda.txt\" using 1:2 with lines\");
}

\/*funciones finales de membresia*\/
double gaussiana(double c, double gama, double x){
 double ux = 0, expo=0;
 expo = pow((x-c)\/gama,2);
 ux = pow(2.7182818284,-0.5*expo);
 return ux;
}

double gaussIzquierda(double c, double gama, double x){
 double ux = 0, expo=0;
 expo = pow((x-c)\/gama,2);
 ux = pow(2.7182818284,-0.5*expo);
 if(x<=c){
  ux=1;
 }
 return ux;
}

double gaussDerecha(double c, double gama, double x){
 double ux = 0, expo=0;
 expo = pow((x-c)\/gama,2);
 ux = pow(2.7182818284,-0.5*expo);
 if(x&gtc){
  ux=1;
 }
 return ux;
}

void graficarPuntosDesdeArchivo(char *cad, char *archivo, char *archivo1, char *archivo2){
 int i;
 \/*lista de comandos para ejecutar y configurar la visualizaci\u00f3n que tendr\u00e1n
     * los puntos en la gr\u00e1fica con gnuplot*\/
    char * configGnuplot[] = {cad, 
                                \"set ylabel \"----grado de membres\u00eda---&gt\"\",
                                \"set xlabel \"----Universo X---&gt\"\",
                                archivo,
                                archivo1,
                                archivo2
                               };
    \/*Se crea una archivo de tipo popen, es una teber\u00eda IPC que se usa, para
     * ejecutar gnuplot y enviarle el archivo a graficar*\/
    \/\/FILE * ventanaGnuplot = popen(\"gnuplot -persist\", \"w\");
    FILE * ventanaGnuplot = popen(\"gnuplot -persist\", \"w\");
    \/\/ Executing gnuplot commands one by one
    for (i=0;i<NUM_COMANDOS;i+ +){
  fprintf(ventanaGnuplot, \"%s n\", configGnuplot[i]);
 }
}
<\/pre>\n<\/div>\n","protected":false},"excerpt":{"rendered":"
\/* * Autor: A. L. I. Fernando Merino * Nombre del archivo: sistemaTakgi.c * Descripci\u00f3n: Simula un sistema difuso, usando el m\u00e9todo de Takgi-Sugeno * genera un archivo con los puntos en pares (x,y) para graficar en gnuplot, * si desean graficar adicionalmente los puntos que est\u00e9n utilizando generen * un archivo extra con los […]<\/p>\n","protected":false},"author":2,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[76,1436,1437,1435,4,1438,1439],"tags":[],"class_list":["post-1818","post","type-post","status-publish","format-standard","hentry","category-c","category-code-c-takagi-sugeno","category-control","category-difuso","category-linux","category-sugeno","category-takagi"],"_links":{"self":[{"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/posts\/1818","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/comments?post=1818"}],"version-history":[{"count":0,"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/posts\/1818\/revisions"}],"wp:attachment":[{"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/media?parent=1818"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/categories?post=1818"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/linuxitos.com\/main\/wp-json\/wp\/v2\/tags?post=1818"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}