matlab_utilities/fuzderantec.cpp

MEX file that gets the derivative of a fuzzy model with respect to its antecedents.

MATLAB help:

% FUZDERANTEC Gets the derivatives of a fuzzy model respect its antecedents.
%
%   It is assumed that the output of the system is the same as the parameter,
%   otherwise it should be considered derivative equal to 0.
%
%   dh_dantec = fuzderconseq(Fuzzy_Model, Point)
%
%   dh_dantec = fuzderconseq(Fuzzy_Model, Point, input, out, rule, parameter)
%
% Arguments:
%
%   Fuzzy_Model -> This fuzzy model could be a '.txt' file, a '.fis' file,
%                  or a 'FIS' variable from MATLAB Workspace.
%
%   Point -> Point where the derivative will be calculate.
%
%   input -> Input for the parameter.
%
%   out -> Output for select the rule for the parameter.
%
%   rule -> Rule for the output 'out' for the parameter.
%
%   parameter -> Parameter of the antecedent.
%
% See also activation, antec2mat, aproxjac, aproxlinear, conseq2mat, fis2txt,
%          fuz2mat, fuzcomb, fuzderconseq, fuzderparam, fuzeval, fuzjac,
%          fuzlinear, fuzprint, mat2antec, mat2conseq, mat2fuz, subsystem,
%          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"
#include <flt/derivatives.hpp>

using namespace FLT;
using namespace TNT;

void mexFunction (int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
  size_t input, output, rule, param, m;
  int error=0;
  double *Point, *dModel_dantec;
  
  if(nrhs!=2 && nrhs!=6)
    ERRORMSG(E_NumberArgIn)
  if (nlhs>1)
    ERRORMSG(E_NumberArgOut)
  for (int i=0;i<nrhs;i++)
  {
    if (mxIsEmpty(prhs[i]) || mxIsNaN(*mxGetPr(prhs[i])))
      ERRORMSG(E_ArgNoValid)
  }
  System S;
  if (readModel(prhs[0],S))
    ERRORMSG(E_Model)

  if (1 != mxGetN(prhs[1]))
    ERRORMSG(E_Column)
  if (S.inputs() != mxGetM(prhs[1]))
    ERRORMSG(E_PointCoherent)
  Point = mxGetPr(prhs[1]);

  if (nrhs==6)
  {
    if ((mxGetN(prhs[2])!=1) || (mxGetM(prhs[2])!=1))
      ERRORMSG(E_In_No_Scalar)
    if ((mxGetN(prhs[3])!=1) || (mxGetM(prhs[3])!=1))
      ERRORMSG(E_Out_No_Scalar)
    if ((mxGetN(prhs[4])!=1) || (mxGetM(prhs[4])!=1))
      ERRORMSG(E_R_No_Scalar)
    if ((mxGetN(prhs[5])!=1) || (mxGetM(prhs[4])!=1))
      ERRORMSG(E_Param_No_Scalar)
      
    input = ((size_t)*mxGetPr(prhs[2]))-1;
    if (input >= S.inputs())
      ERRORMSG(E_InputExceded)

    output = ((size_t)*mxGetPr(prhs[3]))-1;
    if (output >= S.outputs())
      ERRORMSG(E_OutputExceded)

    rule = ((size_t)*mxGetPr(prhs[4]))-1;
    if (rule >= S.rules(output))
      ERRORMSG(E_RulesExceded)
      
    param = ((size_t)*mxGetPr(prhs[5]))-1;
    Membership *M = S.readRule(output, rule)->readFunction(input);
    if (M->type() == ANYMF) // ANYMF doesn't have parameters
    {
      plhs[0]=mxCreateDoubleMatrix(0,0,mxREAL);
      if (!plhs[0])
        ERRORMSG(E_NumberArgOut)
      return;
    }
    if (param >= M->num_params())
      ERRORMSG(E_ParamExceded)

    plhs[0] = mxCreateDoubleScalar(1);
    if (!plhs[0])
      ERRORMSG(E_NumberArgOut)
    dModel_dantec = mxGetPr(plhs[0]);

    (*dModel_dantec) = derantec(S, Point, input, output, rule, param, error);
    if (error)
    {
      mxDestroyArray(plhs[0]);
      ERRORMSG(O_GeneralError)
    }
  }
  else // nrhs=2
  {
    m = S.outputs();
    plhs[0] = mxCreateDoubleMatrix(m,S.NumberOfAntecedents(),mxREAL);
    if (!plhs[0])
      ERRORMSG(E_NumberArgOut)
    dModel_dantec=mxGetPr(plhs[0]);
    TNT::Array2D<double> dS_dantec = derantec(S, Point);
    if (dS_dantec.dim1()==0 || dS_dantec.dim2()==0 )
    {
      mxDestroyArray(plhs[0]);
      ERRORMSG(O_GeneralError)
    }
    for (size_t j=0;j<dS_dantec.dim2();j++)
    {
      for (size_t i=0;i<dS_dantec.dim1();i++)
      {
        *dModel_dantec = dS_dantec[i][j];
        dModel_dantec++;
      }
    }
  }
}