ICIAMExperiment.m

Go to the documentation of this file.

namespace models{
namespace iciam2011{


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

        dim;
    public: /* ( setObservable ) */


        ICIAMExperiment(dims) {
            this.registerProps(" dim ");

            this.dim= dims;

            this.Sampler= [];/* sampling.GridSampler; */


            /*  This class implements a fake Approx subclass to allow access
             * to the this.Ma property for the error estimator. */
            this.Approx= [];

            /* % System settings */
            this.System= models.iciam2011.ICIAMSystem(this);

            this.T= 20;
            this.dt= 0.05;

            /* s = solvers.MLWrapper(@ode45); */
            s = solvers.ExplEuler;
            s.MaxStep= [];
            /* s = solvers.Heun; */
            this.ODESolver= s;

            this.System.addParam(" mu1 ", .5, " Range ", [0 1], " Desired ", 4);
            this.System.addParam(" mu2 ", 1);
            this.Sampler= sampling.GridSampler;

            e1 = [ones(dims,1) zeros(dims,1)];
            e2 = [zeros(dims,1) ones(dims,1)];

            x0 = dscomponents.AffineInitialValue;
            x0.addMatrix(" mu(2) ", ones(dims,1));
            /* x0.addMatrix('sum(mu)', e1); */
            this.System.x0= x0;

            fexp = this.System.f.Expansion;
            fexp.Ma= repmat(-exp(-fexp.Centers.xi(1,:)/15),dims,1);

            /* this.System.Inputs{1} = @(t)[.04*sin(t/3); exp(-abs(10-t/2))]; */
            this.System.Inputs[1] = @(t)[.04*sin(t/3);.5*exp(-(12-t).^2)];
            /*              m.System.Inputs{2} = @(t)exp(-abs(10-t/2)); % Gutes Beispiel!
             *             m.System.Inputs{3} = @(t).5*exp(-(12-t).^2); % Sehr interessante dynamik
             *             m.System.Inputs{4} = @(t)(t>10)*.1; % okay, aber nicht so spannend */
            B = dscomponents.AffLinInputConv;
            B.addMatrix(" mu(1) ", e1);
            B.addMatrix(" 1-mu(1) ", e2);
            /*              B.addMatrix('-mu(2)', e2); */
            this.System.B= B;

            this.TrainingInputs= 1;

            this.System.C= dscomponents.LinearOutputConv(ones(1,dims)/dims);

            V = ones(dims,1)/sqrt(dims);
            s = spacereduction.ManualReduction(V,V);
            /*              s = spacereduction.PODReducer;
             *             s.Value = 3;
             *             s.Mode = 'abs';
             *             s.UseSVDS = dims > 10000;
             *             this.SpaceReducer = s; */

            s = spacereduction.RotationDecorator(s);
            s.Dims= 100;
            s.Degree= 0.05;
            this.SpaceReducer= s;

            e = error.IterationCompLemmaEstimator;
            e.UseTimeDiscreteC= true;
            this.ErrorEstimator= e;
        }

    
    public: /* ( Static ) */ /* ( setObservable ) */


        static function [d , r , m ] = CreatePlots(dim) {
            
            m = ICIAMExperiment(dim);
            s = m.SpaceReducer;

            d = EstimatorAnalyzer;
            d.EstimatorIterations= [1 2 5];
            d.EstimatorVersions= [1 1 0 0 1 0 0 1 1];
            d.SingleFigures= true;

            d.setModel(m);
            r = d.ReducedModel;

            for mu1 = [0 1]
                d.start([mu1; 1], 1);

                /*  1 = full, 7 = LSLE TD */
                md = d.ModelData(end);
                absmax = 30*abs(md.ErrT(7)-md.ErrT(1));
                relmax = 30*abs(md.RelErrT(7)-md.RelErrT(1));

                f = d.Figures[1];
                a = gca(f);
                axis(a,[0 m.T md.MinErr*.9 min(1e4,absmax)]);
                set(legend(a)," Location "," NorthEast ");
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_errors ", mu1, s.Degree)," fig ");
                title(a,);
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_errors ", mu1, s.Degree));

                f = d.Figures[2]; a = gca(f);
                axis(a,[0 m.T md.MinRelErr*.9 min(1,relmax)]);
                set(legend(a)," Location "," NorthEast ");
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_relerr ", mu1, s.Degree)," fig ");
                title(a,);
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_relerr ", mu1, s.Degree));

                f = d.Figures[3]; a = gca(f);
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_ctimes ", mu1, s.Degree)," fig ");
                title(a,);
                Utils.saveFigure(f,sprintf(" ICIAM_mu1_%d_deg%f_ctimes ", mu1, s.Degree));

            end

            PlotParamSweep(r,[1; 1],1,1,-.1:.05:1.1);
            Utils.saveFigure(gcf," mu1_sweep "," png ");

            d.createStatsTables;
        }



};
}
}