PiecewiseLinear.m

Go to the documentation of this file.

namespace general{
namespace functions{


/* (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 PiecewiseLinear
  :public general.functions.AFunGen {
    protected:

        pts;

        vals;

        
    public:

        PiecewiseLinear(colvec<double> x,matrix<double> y) {
            if nargin < 2
                if nargin < 1
                    x = 0:2;
                    y = [0 1 0];
                else
                    y = x;
                    x = 1:size(x,2);
                end
            end
            if length(x) ~= length(unique(x))
                error(" Have to provide disjoint x locations ");
            end
            [x, idx] = sort(x," ascend ");
            y = y(idx);
            this.pts= x;
            this.vals= y;
        }


        function [fhandle , dfhandle ] = getFunction() {
            [x,y] = this.transform(this.pts, this.vals);
            xy = [x;y];
            nf = size(xy,2)-1;
            funs = cell(1,nf);
            dfuns = cell(1,nf);
            for k=1:nf
                xy1 = xy(:,k);
                xy2 = xy(:,k+1);
                cond = sprintf(" (t>=%g & t <%g) ",xy1(1),xy2(1));
                a = (xy2(2)-xy1(2))/(xy2(1)-xy1(1));
                b = xy1(2)-a*xy1(1);
                fun = sprintf(" (%g*t+%g) ",a,b);
                funs[k] = [cond " .* " fun];
                dfuns[k] = [cond " .* " sprintf(" %g ",a)];
            end
            funstr = sprintf(" @(t)%s ",Utils.implode(funs,"  +  "));
            dfunstr = sprintf(" @(t)%s ",Utils.implode(dfuns,"  +  "));
            fhandle = eval(funstr);
            dfhandle = eval(dfunstr);
        }


        function str = getConfigStr() {
            str = sprintf(" Points: %d ",size(this.pts,2));
        }


        function pm = plot(range) {
            if nargin < 2
                x = this.transform(this.pts,this.vals);
                range = [min(x) max(x)];
            end
            pm = plot@general.functions.AFunGen(this, " R ", range);
        }


    
    protected:

        function [colvec<double>x , matrix<double>y ] = transform(colvec<double> x,matrix<double> y) {
        }
    public: /* ( Static ) */

        static function  test_PiecewiseLinear() {
            p = general.functions.PiecewiseLinear;
            p.plot;
            p = general.functions.PiecewiseLinear(rand(1,7));
            p.plot;
            p = general.functions.PiecewiseLinear(rand(1,7),rand(1,7));
            p.plot;
        }

    

};
}
}