RandomModelEstimatorAnalyzer.m

Go to the documentation of this file.

namespace demos{


/* (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 RandomModelEstimatorAnalyzer
  :public EstimatorAnalyzer {
    public:

        Dims = 500;
        NumCenters = 10;
        PositiveExpansion = false;
        UniformExpansion = false;
    public:


        RandomModelEstimatorAnalyzer(dims) {

            /* % Model settings */
            fm = models.BaseFullModel;
            fm.Name= " Estimator Demo Model ";

            fm.T= 1;
            fm.dt= .025;

            fm.Approx= [];
            fm.Sampler= [];

            /* this.ODESolver = solvers.MLWrapper(@ode45); */
            fm.ODESolver= solvers.ExplEuler(fm.dt);
            /* fm.ODESolver = solvers.Heun(fm.dt); */

            /* % Core function */
            cf = dscomponents.ParamTimeKernelCoreFun;
            kexp = cf.Expansion;
            kexp.TimeKernel= kernels.NoKernel;
            kexp.ParamKernel= kernels.NoKernel;
            kexp.Centers.ti= [];
            kexp.Centers.mui= [];

            /* % System settings */
            sys = models.BaseFirstOrderSystem(fm);
            sys.MaxTimestep= fm.dt;
            sys.f= cf;

            fm.System= sys;
            this.Model= fm;

            if nargin == 1
                this.Dims= dims;
            else
                this.newCoeffs;
            end
        }
        function  setup() {
            k = kernels.GaussKernel(15);
            k.G= 1;
            this.Model.System.f.Expansion.Kernel= k;
            x0 = rand(this.Dims,1);
            if this.PositiveExpansion
                base = linspace(0, 40, this.NumCenters);
                this.Model.System.x0= dscomponents.ConstInitialValue(x0);
            else
                base = linspace(-20, 20, this.NumCenters);
                this.Model.System.x0= dscomponents.ConstInitialValue(x0-.5);
            end
            this.Model.System.f.Expansion.Centers.xi= repmat(base,this.Dims,1);

            if  this.UniformExpansion
                V = ones(this.Dims,1)*sqrt(1/this.Dims);
                s = spacereduction.ManualReduction(V,V);
            else
                s = spacereduction.PODReducer;
                s.Mode= " abs ";
                s.Value= 1;
                s.UseSVDS= false;
            end
            this.Model.SpaceReducer= s;

            /* % Generation */
            this.setModel(this.Model);
        }


        function  newCoeffs() {
            offset = .5;
            if this.PositiveExpansion
                offset = 0;
            end
            /*  Create coefficients */
            if this.UniformExpansion
                ai = (rand(1,this.NumCenters)-offset);
                this.Model.System.f.Expansion.Ma= repmat(ai,this.Dims,1);
            else
                this.Model.System.f.Expansion.Ma= (rand(this.Dims,this.NumCenters)-offset);
            end
            this.setup;
        }
#if 0 //mtoc++: 'set.Dims'
function  Dims(value) {
            this.Dims= value;
            this.newCoeffs;
            /* this.setup; %#ok */
        }

#endif


        
#if 0 //mtoc++: 'set.NumCenters'
function  NumCenters(value) {
            this.NumCenters= value;
            this.newCoeffs;
            /* this.setup; %#ok */
        }

#endif


        
#if 0 //mtoc++: 'set.PositiveExpansion'
function  PositiveExpansion(value) {
            this.PositiveExpansion= value;
            /* this.setup; %#ok */
            this.newCoeffs;
        }

#endif


        
#if 0 //mtoc++: 'set.UniformExpansion'
function  UniformExpansion(value) {
            this.UniformExpansion= value;
            this.newCoeffs;
            /* this.setup; %#ok */
        }

#endif


};
}