sofmsimu.c

nnToolKit 神经网络工具包是基于 MATLAB 神经网络工具箱自行开发的一组神经网络算法函数库

/*
 * MATLAB Compiler: 3.0
 * Date: Sun May 13 16:47:40 2007
 * Arguments: "-B" "macro_default" "-O" "all" "-O" "fold_scalar_mxarrays:on"
 * "-O" "fold_non_scalar_mxarrays:on" "-O" "optimize_integer_for_loops:on" "-O"
 * "array_indexing:on" "-O" "optimize_conditionals:on" "-M" "-silentsetup" "-d"
 * "d:/MATLAB6p5/work/nnToolKit/src" "-B" "csglcom:nnToolKit,nnToolKit,2.0"
 * "-B" "sgl" "-m" "-W" "main" "-L" "C" "-t" "-T" "link:exe" "-h"
 * "libmmfile.mlib" "-W" "mainhg" "libmwsglm.mlib" "-t" "-W"
 * "comhg:nnToolKit,nnToolKit,2.0" "-T" "link:lib" "-h" "libmmfile.mlib" "-i"
 * "-i" "D:/MATLAB6p5/work/nnToolKit/lmnet/LmSimu.m"
 * "D:/MATLAB6p5/work/nnToolKit/lmnet/LmTrain.m"
 * "D:/MATLAB6p5/work/nnToolKit/sofm/SofmSimu.m"
 * "D:/MATLAB6p5/work/nnToolKit/sofm/SofmTrain.m" 
 */
#include "sofmsimu.h"
#include "libmatlbm.h"
#include "libmmfile.h"
#include "libmwsglm.h"
#include "nbgrid.h"
#include "nntwarn.h"
#include "simusm.h"

static mxChar _array1_[3] = { 'O', 'F', 'F' };
static mxArray * _mxarray0_;
static mxArray * _mxarray2_;

static mxChar _array4_[5] = { 'w', '%', 's', '%', 's' };
static mxArray * _mxarray3_;

static mxChar _array6_[4] = { '.', 'd', 'a', 't' };
static mxArray * _mxarray5_;

static mxChar _array8_[1] = { 'r' };
static mxArray * _mxarray7_;

static mxChar _array10_[2] = { '%', 'f' };
static mxArray * _mxarray9_;
static mxArray * _mxarray11_;
static mxArray * _mxarray12_;

static mxChar _array14_[3] = { 'a', 'l', 'l' };
static mxArray * _mxarray13_;

void InitializeModule_sofmsimu(void) {
    _mxarray0_ = mclInitializeString(3, _array1_);
    _mxarray2_ = mclInitializeDouble(-1.0);
    _mxarray3_ = mclInitializeString(5, _array4_);
    _mxarray5_ = mclInitializeString(4, _array6_);
    _mxarray7_ = mclInitializeString(1, _array8_);
    _mxarray9_ = mclInitializeString(2, _array10_);
    _mxarray11_ = mclInitializeDouble(1.0);
    _mxarray12_ = mclInitializeDouble(0.0);
    _mxarray13_ = mclInitializeString(3, _array14_);
}

void TerminateModule_sofmsimu(void) {
    mxDestroyArray(_mxarray13_);
    mxDestroyArray(_mxarray12_);
    mxDestroyArray(_mxarray11_);
    mxDestroyArray(_mxarray9_);
    mxDestroyArray(_mxarray7_);
    mxDestroyArray(_mxarray5_);
    mxDestroyArray(_mxarray3_);
    mxDestroyArray(_mxarray2_);
    mxDestroyArray(_mxarray0_);
}

static mxArray * Msofmsimu(int nargout_,
                           mxArray * ModelNo,
                           mxArray * NetPara,
                           mxArray * SimuData,
                           mxArray * DataDir);

_mexLocalFunctionTable _local_function_table_sofmsimu
  = { 0, (mexFunctionTableEntry *)NULL };

/*
 * The function "mlfSofmsimu" contains the normal interface for the "sofmsimu"
 * M-function from file "d:\matlab6p5\work\nntoolkit\sofm\sofmsimu.m" (lines
 * 1-33). This function processes any input arguments and passes them to the
 * implementation version of the function, appearing above.
 */
mxArray * mlfSofmsimu(mxArray * ModelNo,
                      mxArray * NetPara,
                      mxArray * SimuData,
                      mxArray * DataDir) {
    int nargout = 1;
    mxArray * retstr = NULL;
    mlfEnterNewContext(0, 4, ModelNo, NetPara, SimuData, DataDir);
    retstr = Msofmsimu(nargout, ModelNo, NetPara, SimuData, DataDir);
    mlfRestorePreviousContext(0, 4, ModelNo, NetPara, SimuData, DataDir);
    return mlfReturnValue(retstr);
}

/*
 * The function "mlxSofmsimu" contains the feval interface for the "sofmsimu"
 * M-function from file "d:\matlab6p5\work\nntoolkit\sofm\sofmsimu.m" (lines
 * 1-33). The feval function calls the implementation version of sofmsimu
 * through this function. This function processes any input arguments and
 * passes them to the implementation version of the function, appearing above.
 */
void mlxSofmsimu(int nlhs, mxArray * plhs[], int nrhs, mxArray * prhs[]) {
    mxArray * mprhs[4];
    mxArray * mplhs[1];
    int i;
    if (nlhs > 1) {
        mlfError(
          mxCreateString(
            "Run-time Error: File: sofmsimu Line: 4 Column:"
            " 1 The function \"sofmsimu\" was called with m"
            "ore than the declared number of outputs (1)."),
          NULL);
    }
    if (nrhs > 4) {
        mlfError(
          mxCreateString(
            "Run-time Error: File: sofmsimu Line: 4 Column:"
            " 1 The function \"sofmsimu\" was called with m"
            "ore than the declared number of inputs (4)."),
          NULL);
    }
    for (i = 0; i < 1; ++i) {
        mplhs[i] = NULL;
    }
    for (i = 0; i < 4 && i < nrhs; ++i) {
        mprhs[i] = prhs[i];
    }
    for (; i < 4; ++i) {
        mprhs[i] = NULL;
    }
    mlfEnterNewContext(0, 4, mprhs[0], mprhs[1], mprhs[2], mprhs[3]);
    mplhs[0] = Msofmsimu(nlhs, mprhs[0], mprhs[1], mprhs[2], mprhs[3]);
    mlfRestorePreviousContext(0, 4, mprhs[0], mprhs[1], mprhs[2], mprhs[3]);
    plhs[0] = mplhs[0];
}

/*
 * The function "Msofmsimu" is the implementation version of the "sofmsimu"
 * M-function from file "d:\matlab6p5\work\nntoolkit\sofm\sofmsimu.m" (lines
 * 1-33). It contains the actual compiled code for that M-function. It is a
 * static function and must only be called from one of the interface functions,
 * appearing below.
 */
/*
 * %此为Sofm网络处理程序
 * %自组织特征映射模型(Self-Organizing feature Map)，认为一个神经网络接受外界输入模式时，将会分为不同的区域，各区域对输入模式具有不同的响应特征，同时这一过程是自动完成的。各神经元的连接权值具有一定的分布。最邻近的神经元互相刺激，而较远的神经元则相互抑制，更远一些的则具有较弱的刺激作用。自组织特征映射法是一种无教师的聚类方法。
 * %完成分类训练,并保存权值和分类后各类别下的像素矩阵
 * function retstr = SofmSimu(ModelNo,NetPara,SimuData,DataDir)
 */
static mxArray * Msofmsimu(int nargout_,
                           mxArray * ModelNo,
                           mxArray * NetPara,
                           mxArray * SimuData,
                           mxArray * DataDir) {
    mexLocalFunctionTable save_local_function_table_
      = mclSetCurrentLocalFunctionTable(&_local_function_table_sofmsimu);
    mxArray * retstr = NULL;
    mxArray * c = NULL;
    mxArray * a = NULL;
    mxArray * m = NULL;
    mxArray * count = NULL;
    mxArray * w = NULL;
    mxArray * frw = NULL;
    mxArray * ClassifyNum = NULL;
    mxArray * olddir = NULL;
    mxArray * ans = NULL;
    mclCopyArray(&ModelNo);
    mclCopyArray(&NetPara);
    mclCopyArray(&SimuData);
    mclCopyArray(&DataDir);
    /*
     * NNTWARN OFF
     */
    mlfNntwarn(_mxarray0_);
    /*
     * retstr=-1;
     */
    mlfAssign(&retstr, _mxarray2_);
    /*
     * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     * %ModelNo='1';
     * %DataDir='.';
     * % 网络参数
     * %NetPara(1)=1;      %输入层节点数
     * %NetPara(2)=5;      %分类数
     * %SimuData=0.5882;
     * %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     * %保留原目录
     * olddir=pwd;
     */
    mlfAssign(&olddir, mlfPwd());
    /*
     * %进入数据所在目录
     * cd(DataDir);
     */
    mclAssignAns(&ans, mlfNCd(0, mclVa(DataDir, "DataDir")));
    /*
     * 
     * ClassifyNum=NetPara(2);    %分类数
     */
    mlfAssign(&ClassifyNum, mclIntArrayRef1(mclVa(NetPara, "NetPara"), 2));
    /*
     * % 读入权值文件
     * frw=fopen(sprintf('w%s%s',ModelNo,'.dat'),'r');
     */
    mlfAssign(
      &frw,
      mlfFopen(
        NULL,
        NULL,
        mlfSprintf(
          NULL, _mxarray3_, mclVa(ModelNo, "ModelNo"), _mxarray5_, NULL),
        _mxarray7_,
        NULL));
    /*
     * [w,count]=fscanf(frw,'%f',[ClassifyNum,1]);
     */
    mlfAssign(
      &w,
      mlfFscanf(
        &count,
        mclVv(frw, "frw"),
        _mxarray9_,
        mlfHorzcat(mclVv(ClassifyNum, "ClassifyNum"), _mxarray11_, NULL)));
    /*
     * fclose(frw);
     */
    mclAssignAns(&ans, mlfFclose(mclVv(frw, "frw")));
    /*
     * 
     * m=nbgrid(ClassifyNum);
     */
    mlfAssign(
      &m, mlfNbgrid(mclVv(ClassifyNum, "ClassifyNum"), NULL, NULL, NULL, NULL));
    /*
     * a =simusm(SimuData,w,m,0);
     */
    mlfAssign(
      &a,
      mlfSimusm(
        mclVa(SimuData, "SimuData"),
        mclVv(w, "w"),
        mclVv(m, "m"),
        _mxarray12_));
    /*
     * c=find(a);
     */
    mlfAssign(&c, mlfFind(NULL, NULL, mclVv(a, "a")));
    /*
     * 
     * cd(olddir);
     */
    mclAssignAns(&ans, mlfNCd(0, mclVv(olddir, "olddir")));
    /*
     * 
     * retstr=c;
     */
    mlfAssign(&retstr, mclVv(c, "c"));
    /*
     * close all;
     */
    mclAssignAns(&ans, mlfNClose(0, _mxarray13_, NULL));
    mclValidateOutput(retstr, 1, nargout_, "retstr", "sofmsimu");
    mxDestroyArray(ans);
    mxDestroyArray(olddir);
    mxDestroyArray(ClassifyNum);
    mxDestroyArray(frw);
    mxDestroyArray(w);
    mxDestroyArray(count);
    mxDestroyArray(m);
    mxDestroyArray(a);
    mxDestroyArray(c);
    mxDestroyArray(DataDir);
    mxDestroyArray(SimuData);
    mxDestroyArray(NetPara);
    mxDestroyArray(ModelNo);
    mclSetCurrentLocalFunctionTable(save_local_function_table_);
    return retstr;
}