matlab_utilities/mat2fuz.cpp

MEX file that changes all data of a fuzzy model.

MATLAB help:

% MAT2FUZ Changes all data of a fuzzy model.
%
%   Final_Model = mat2fuz(Initial_Model,Vector)
%
%   Final_Model = mat2fuz(Initial_Model,Vector,Output,Rule)
%
%   [Final_Model,length] = mat2fuz(Initial_Model,Vector)
%
%   [Final_Model,length] = mat2fuz(Initial_Model,Vector,Output,Rule)
%
% Arguments:
%
%   Final_Model -> Changed fuzzy model.
%
%   Output -> Use to select the output.
%
%   Rule -> Use to select the rule for given output.
%
%   length -> The length of Vector, 0 if error.
%
%   Initial_Model -> Fuzzy model to change. This model could be a '.txt' file,
%                    a '.fis' file, or a 'FIS' variable from MATLAB Workspace.
%
%   Vector -> Vector with data for fuzzy model.
%
% See also activation, antec2mat, aproxjac, aproxlinear, conseq2mat, fis2txt,
%          fuz2mat, fuzcomb, fuzeval, fuzjac, fuzlinear, fuzprint, mat2antec,
%          mat2conseq, fuzsubsystem, txt2fis

Source code:

/*  Copyright (C) 2004-2011
  ANTONIO JAVIER BARRAGAN, antonio.barragan@diesia.uhu.es
  http://uhu.es/antonio.barragan

  Collaborators:
  JOSE MANUEL ANDUJAR, andujar@diesia.uhu.es
  MARIANO J. AZNAR, marianojose.aznar@alu.uhu.es

  DPTO. DE ING. ELECTRONICA, DE SISTEMAS INFORMATICOS Y AUTOMATICA
  ETSI, UNIVERSITY OF HUELVA (SPAIN)

  For more information, please contact with authors.
  
  This software is free software: you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.

  This software is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/

#include "aux_matlab.hpp"

using namespace FLT;

void mexFunction (int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
  size_t i, length, out, rule;
  double *data;
  
  if(nrhs<2)
    ERRORMSG(E_NumberArgIn)
  if (nrhs>2)
  {
    if (nrhs!=4)
      ERRORMSG(E_NumberArgIn)
    if (mxGetM(prhs[2])>1 || mxGetM(prhs[3])>1)
      ERRORMSG(E_RuleOutputFileVector)
    if (mxGetN(prhs[2])!=mxGetN(prhs[3]))
      ERRORMSG(E_NumberArg)
  }
    
  if (nlhs>2)
    ERRORMSG(E_NumberArgOut)
  if (mxIsEmpty(prhs[0]) || mxIsNaN(*mxGetPr(prhs[0])))
    ERRORMSG(E_ArgNoValid)
    
  System S;
  if (readModel(prhs[0],S))
    ERRORMSG(E_Model)
  
  size_t NoC = S.inputs() + 1;

  // Save data in S
  data = mxGetPr(prhs[1]);
  int error = 0;
  if (nrhs==2)
  {
    error = S.setAntecedents(data);
    length = S.NumberOfAntecedents();
    S.setConsequents(data + length);
    length += S.NumberOfConsequents();
  }
  else
  {
    length = 0;
    for (i=0;i<mxGetN(prhs[2]);i++)
    {     
      out = ((size_t)*(mxGetPr(prhs[2])+i))-1; // MATLAB is index-1 but C/C++ is index-0
      rule = ((size_t)*(mxGetPr(prhs[3])+i))-1;
      
      Rule *R = S.readRule(out,rule);
      
      size_t numAntec = R->NumberOfAntecedents();
      error += R->setAntecedents(data);
      data += numAntec;
      
      R->setConsequents(data);
      data += NoC;

      length += numAntec + NoC;
    }

  }
  
  if (error>0)
  {
    ERRORMSG(E_Store)
    plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL); // Empty matrix
    return;
  }
  
  if (nlhs==2)
    plhs[1] = mxCreateDoubleScalar(length);
    

   // Create the FIS variable
  const char* names[] = {"name","type","andMethod","orMethod","defuzzMethod","impMethod","aggMethod","input","output","rule"};
  mwSize dim[] = {1,1};
  plhs[0] = mxCreateStructArray(2,dim,10,names);
  
  if(System2FIS(S, plhs[0]))
  {
    mxDestroyArray(plhs[0]);
    ERRORMSG(E_CreateFIS)
  }
}