RCLadder.m

Go to the documentation of this file.

namespace models{
namespace circ{


/* (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 RCLadder
  :public models.BaseFullModel {
    public:

        Dims;
    public:


        RCLadder(dims) {
            
            this = this@models.BaseFullModel;

            if nargin == 0
                dims = 30;
            end
            this.Dims= dims;

            this.registerProps(" Dims ");

            this.Name= " RC Ladder circuit model ";

            this.T= 3;
            this.dt= .0025;
            this.tau= 1;
            this.SaveTag= sprintf(" rcladder_d%d_T%g ",dims,this.T);
            this.Data= data.ModelData(this);

            this.System= models.circ.RCLadderSys(this);

            /*  Only train with first input! */
            this.TrainingInputs= [2 3];

            this.Sampler= [];

            app = approx.KernelApprox(this.System);
            a = approx.algorithms.VKOGA;
            a.MaxRelErr= 1e-5;
            a.MaxAbsResidualErr= 1e-3;
            ec = kernels.config.ExpansionConfig;
            ec.StateConfig= kernels.config.GaussConfig(" D ",.3:.1:2);
            a.ExpConfig= ec;
            app.Algorithm= a;

            t = data.selection.TimeSelector;
            t.Size= 12000;
            app.TrainDataSelector= t;
            this.Approx= app;

            s = spacereduction.PODReducer;
            s.Mode= " abs ";
            s.Value= 3;
            s.UseSVDS= false;
            this.SpaceReducer= s;

            s = solvers.SemiImplicitEuler(this);
            s.MaxStep= .05; /*  Stability constraint due to diffusion term */

            this.ODESolver= s;

            this.ErrorEstimator= error.IterationCompLemmaEstimator;
        }


        
#if 0 //mtoc++: 'set.Dims'
function  Dims(value) {
            if ~isposintscalar(value)
                error(" value must be a positive integer scalar ");
            end
            this.Dims= value;
        }

#endif

    
    public: /* ( Static ) */

        static function  test_RCLadderModel() {
            m = models.circ.RCLadder(100);
            m.simulate([], 1);
            [t,y,~,x] = m.simulate([], 2);
            m.plot(t,y);
            m.plotState(t,x);
        }

    
    protected: /* ( Static ) */

        static function this = loadobj() {
            if ~isa(this, " models.circ.RCLadder ")
                s = this;
                this = models.circ.RCLadder(s.Dims);
                /*  field "Dims" is set in constructor */
                this = loadobj@models.BaseFullModel(this, s);
            else
                this = loadobj@models.BaseFullModel(this);
            end
        }

};
}
}