Force.m

Go to the documentation of this file.

namespace models{
namespace golf{


/* (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 Force
  :public dscomponents.ACoreFun {
    public:


        
    public:

        Force(sys) {
            this = this@dscomponents.ACoreFun(sys);
            this.xDim= 4;
            this.fDim= 2;
            this.JSparsityPattern= sparse(true(2,4));
        }


        function dv = evaluate(x_v,double t) {
            m = this.System.Model;
            mu = this.mu;

            gr = m.gradgreen(x_v(1:2));
            /*  Normale berechnen */
            normale = [gr; -1]/norm([gr; -1]);
            /*  Normalenkraft */
            Fn = normale * -mu(7)*normale(3);
            /*  Hangabtriebskraft */
            Fb = [0;0;-mu(7)] - Fn;
            /*  Dreht den Fb-Vektror in die x,y-Richtung (kleiner Optimierungsgedanke)
             *Fb = -gr/norm(gr)*norm(Fb); */

            /* % MIT Reibung
             * Reibungskonstante * |F_n| * normierter Geschwindigkeitsvektor * zeiteinheit
             *Fr = mu*norm([gr;1]) * (v / norm(v)); %ALT */
            Fr = mu(8)*norm(Fn) * (x_v(3:4) / norm(x_v(3:4)));
            /*  Gesamtkraft; wir lassen die z-Komponente der Hangabtriebskraft
             * ausser acht, da diese auf einem Putting-Green erwartungsgemäß
             * klein ist. */
            dv = Fb([1 2])-Fr;

            /* % OHNE Reibung
             *dv = Fb; */
        }


        function  evaluateCoreFun() {
            error(" evaluate is overridden directly. ");
        }

    

};
}
}