MathMODSystem.m

Go to the documentation of this file.

namespace models{
namespace mathmod2012{


/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

class MathMODSystem
  :public models.BaseFirstOrderSystem {
    public: /* ( setObservable ) */

        svNum;
    public: /* ( setObservable ) */


        MathMODSystem(models.BaseFullModel model) {
            
            this = this@models.BaseFirstOrderSystem(model);

            this.registerProps(" svNum ");

            /* % System settings */
            dims = model.dim;
            this.NumStateDofs= dims;

            this.MaxTimestep= [];

            this.addParam(" input shift ",1);
            this.addParam(" expansion param ",5," Range ",[0 10]);
            this.addParam(" initial value ",.5," Range ",[0 1]);
            this.DependentParamIndices= 2;

            /*  Sample bases */
            this.svNum= 20;

            f = dscomponents.ParamTimeKernelCoreFun(this);
            fe = f.Expansion;
            this.f= f;

            /*  Choose the kernel with such that each kernel vanishes (<
             * kerneleps) after kernelsupport support vectors.  */
            kernelsupport = 2;
            kerneleps = 0.00001;

            /*  Space kernel */
            space = linspace(0,50,this.svNum);
            fe.Centers.xi= repmat(space,dims,1);
            d = sqrt(dims)*space(2);
            fe.Kernel= kernels.GaussKernel;
            fe.Kernel.setGammaForDistance(kernelsupport*d,kerneleps)
            fe.Kernel.G= 1;

            /*  Time kernel */
            fe.TimeKernel= kernels.NoKernel;
            fe.Centers.ti= [];

            /*  Param kernel */
            pk = kernels.GaussKernel;
            pspace = linspace(0,10,this.svNum);
            /*  Only the 2nd entry is used, so fill rest with zeros */
            pk.P= 2;
            pk.G= 1;
            pk.setGammaForDistance(sqrt(this.ParamCount)*pspace(2)*20,kerneleps);
            fe.ParamKernel= pk;
            fe.Centers.mui= [zeros(1,this.svNum); pspace; zeros(1,this.svNum)];

            this.updateDimensions;
            this.updateSparsityPattern;
        }


};
}
}