KernelTestFun.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 KernelTestFun
  :public dscomponents.ParamTimeKernelCoreFun {
    public: /* ( setObservable ) */

        svNum;
    public: /* ( setObservable ) */

        KernelTestFun(sys,pos_flag) {
            
            this = this@dscomponents.ParamTimeKernelCoreFun(sys);
            this.registerProps(" svNum ");

            if nargin == 1
                pos_flag = false;
            end

            /*  Sample bases */
            this.svNum= 20;
            kexp = kernels.ParamTimeKernelExpansion;
            kexp.Centers.xi= repmat(linspace(-24,24,this.svNum),sys.Model.dim,1);
            kexp.Centers.ti= .5*(1:this.svNum);
            kexp.Centers.mui= rand(2,this.svNum);

            /*  Function coefficients */
            offset = .5;
            /* offset = 0; */
            if pos_flag
                offset = 0;
            end
            ai = (rand(1,this.svNum)-offset);

            kexp.Ma= repmat(ai,sys.Model.dim,1);

            /* % KernelApprox settings
             *this.Kernel = kernels.GaussKernel(4*sqrt(dims)); */
            kexp.Kernel= kernels.GaussKernel(sqrt(sys.Model.dim));
            kexp.Kernel.G= 1;
            kexp.TimeKernel= kernels.GaussKernel(.4);
            kexp.TimeKernel.G= 1;
            kexp.ParamKernel= kernels.GaussKernel(.4);
            kexp.ParamKernel.G= 1;
            this.Expansion= kexp;
        }


        function  showBaseFun() {
            figure;
            x = repmat(linspace(-50,50,max(this.svNum*5,100)),this.System.Model.dim,1);
            fx = this.evaluate(x,0,[]);
            plot(x(1,:),fx(1,:)," r ");
            xlabel(" x "); ylabel(" f(x) ");
            title(" KernelTest base function ");
        }
};
}
}