Rotation.m

Go to the documentation of this file.

namespace models{
namespace synth{


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

        Rotation() {
            this.T= 10;
            this.dt= .025;

            /* s = sampling.RandomSampler;
             *s.Samples = 10;
             *this.Sampler = s; */
            this.Sampler= sampling.GridSampler;

            this.System= models.synth.RotationDynSys(this);
            this.Name= " Synthetic rotation model ";

            this.SpaceReducer= [];
/*              
             * Core Approximation
 *             a = approx.algorithms.Componentwise;
 *             a.CoeffComp = general.regression.KernelLS;
 *             a.TimeKernel = kernels.GaussKernel;
 *             a.Kernel = kernels.GaussKernel(2);
 *             a.ParamKernel = kernels.LinearKernel;
 *             a.lambda = 2; */

/*              a = approx.algorithms.Componentwise;
 *             a.CoeffComp = general.regression.ScalarEpsSVR;
 *             a.TimeKernel = kernels.GaussKernel;
 *             a.Kernel = kernels.GaussKernel(2);
 *             a.ParamKernel = kernels.LinearKernel;
 *             a.eps = .05;
 *             a.C = 100;
 *             this.Approx = a; */

            a = approx.KernelApprox(this.System);
            a.Expansion.TimeKernel= kernels.GaussKernel;
            a.Expansion.Kernel= kernels.GaussKernel(2);
            /* a.ParamKernel = kernels.LinearKernel; */
            a.Expansion.ParamKernel= kernels.GaussKernel(2);
            this.Approx= a;

            this.ODESolver= solvers.MLWrapper(@ode45);
        }


        function  plot(double t,matrix<double> y) {
            pm = PlotManager(false,1,2);
            ax = pm.nextPlot(" 3d "," 3D plot over time (z-axis) "," x_1 "," x_2 ");
            plot3(ax,y(2,:),y(1,:),t," r ");
            ax = pm.nextPlot(" 1d "," Plot over time "," t "," x_1,x_2 ");
            plot(ax,t,y(2,:)," r ",t,y(1,:)," b ");
            pm.LeaveOpen= true;
            pm.done;
        }


    
    public: /* ( Static ) */

        static function  test_Rotation_Simulation() {
            m = models.synth.Rotation;
            for k = 1:4
                mu = m.getRandomParam;
                [t,y] = m.simulate(mu);
                m.plot(t,y);
            end
        }

    

};
}
}