Pool.m

Go to the documentation of this file.

namespace models{
namespace motorunit{


/* (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 Pool
  :public handle {
    public:

        FibreTypes;

        Version;

        ShortenInput = @("t")3;

        
    public:

        fugle;

        shorten;

        alpha;

        
    public:


        Pool(version) {
            if nargin < 1
                version = 1;
            end
            this.Version= version;
            this.fugle= models.motorunit.Fuglevand;
            this.shorten= models.motorunit.Shorten;
        }


        function  prepareSimulation(T,meancurrent) {
            ft = this.FibreTypes;
            nf = length(ft);
            a = zeros(nf,T+300);
            this.fugle.Seed= 1;
            if this.Version == 1
                for fidx = 1:nf
                    a(fidx,:) = this.fugle.getForces([ft(fidx); meancurrent],T+300);
                end
                startpos = find(sum(a,1),1)-10;
                if startpos > 0
                    a(:,1:startpos) = [];
                    a(:,end-300+startpos+1:end) = [];
                end
            elseif this.Version == 2
                m = this.shorten;
                m.T= T;
                m.System.Inputs[1] = this.ShortenInput;
                for fidx = 1:nf
                    a(fidx,:) = this.shorten.simulate([ft(fidx); meancurrent],1);
                end
            end
            this.alpha= a;
            mval = max(a(:));
            if mval > 0
                this.alpha= a/mval;
            end
        }


        function alpha = getActivation(double t) {
            a = this.alpha;
            alpha = a(:,min(size(a,2),round(t+1)));
        }

    

};
}
}