MEX file that gets the derivative of a fuzzy model with respect to its consequents.
MATLAB help:
% FUZDERCONSEQ Gets the derivative of a fuzzy model respect its consequents.
%
% 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_dconseq = fuzderconseq(Fuzzy_Model, Point)
%
% dh_dconseq = fuzderconseq(Fuzzy_Model, Point, 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.
%
% out -> Output for select the rule for the parameter.
%
% rule -> Rule for the output 'out' for the parameter.
%
% parameter -> Parameter of the consequent (0 -> affine consequent)
%
% See also activation, antec2mat, aproxjac, aproxlinear, conseq2mat, fis2txt,
% fuz2mat, fuzcomb, fuzderantec, fuzderparam, fuzeval, fuzjac, fuzlinear,
% fuzprint, mat2antec, mat2conseq, mat2fuz, fuzsubsystem, txt2fis
Source code:
void mexFunction (int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
size_t output, rule, parameter;
int error=0;
double *Point, *dModel_dconseq;
if(nrhs!=2 && nrhs!=5)
if (nlhs>1)
for (int i=0;i<nrhs;i++)
{
if (mxIsEmpty(prhs[i]) || mxIsNaN(*mxGetPr(prhs[i])))
}
if (mxGetN(prhs[1])!=1)
if (n != mxGetM(prhs[1]))
Point = mxGetPr(prhs[1]);
if (nrhs==5)
{
if ( (mxGetN(prhs[2])!=1) || (mxGetM(prhs[2])!=1))
if ((mxGetN(prhs[3])!=1) || (mxGetM(prhs[3])!=1))
if ( (mxGetN(prhs[4])!=1) || (mxGetM(prhs[4])!=1))
output = ((size_t)*mxGetPr(prhs[2]))-1;
if (output >= m)
rule = ((size_t)*mxGetPr(prhs[3]))-1;
if (rule >= S.
rules(output))
parameter = ((size_t)*mxGetPr(prhs[4]));
plhs[0] = mxCreateDoubleScalar(1);
if (!plhs[0])
dModel_dconseq = mxGetPr(plhs[0]);
(*dModel_dconseq) =
derconseq(S, Point, output, rule, parameter, error);
if (error)
{
mxDestroyArray(plhs[0]);
}
}
else
{
TNT::Array2D<double> dS_dconseq =
derconseq(S, Point);
if (dS_dconseq.dim1()==0 || dS_dconseq.dim2()==0 )
{
mxDestroyArray(plhs[0]);
}
plhs[0] = mxCreateDoubleMatrix(dS_dconseq.dim1(), dS_dconseq.dim2(),mxREAL);
if (!plhs[0])
dModel_dconseq = mxGetPr(plhs[0]);
for (size_t j=0;j<dS_dconseq.dim2();j++)
{
for (size_t i=0;i<dS_dconseq.dim1();i++)
{
*dModel_dconseq = dS_dconseq[i][j];
dModel_dconseq++;
}
}
}
}