MaxForce.m

Go to the documentation of this file.

namespace models{
namespace motorunit{
namespace experiments{


/* (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.
 */

function  MaxForce() {

file = fullfile(fileparts(mfilename(" fullpath "))," maxforces ");

/*  % Dynamic IC: yes, SinglePeak: no, OutScaling: no
 * m = models.motorunit.Shorten(true, false, false);
 * 
 * m.T = 2000;
 * m.dt = .1;
 * m.EnableTrajectoryCaching = false;
 * 
 * % s = m.Sampler;
 * % s.Samples = 1000;
 * % s.Seed = 1;
 * % m.off1_createParamSamples;
 * % p = m.Data.ParamSamples;
 * % % Sort samples by fibre type
 * % [~, idx] = sort(p(1,:));
 * % p = p(:,idx);
 * % n = size(p,2);
 * 
 * % Directly set parameter set
 * n = 100;
 * p = [linspace(0,1,n); 10*ones(1,n)];
 * 
 * maxvals = zeros(1,n);
 * maxpos = maxvals;
 * pi = ProcessIndicator('Gathering max forces',n);
 * parfor k=1:n
 *     [t, y] = m.simulate(p(:,k),1);%#ok
 *     [maxvals(k), pos] = max(y(2,:));
 *     maxpos(k) = pos;
 *     if pos == m.T+1
 *         fprintf('Max at final time for param %d\n',k);
 *     end
 *     pi.step;%#ok
 * end
 * pi.stop;
 * save(file,'m','p','maxvals','maxidx','n'); */

load(file);
pm = PlotManager;
pm.LeaveOpen= true;

/*  Load mean current limiting polynomial */
s = load(models.motoneuron.Model.FILE_UPPERLIMITPOLY);
eff_mean_current = polyval(s.upperlimit_poly,p(1,:));

h = pm.nextPlot(" all "," fibre type "," values ");
plot(h,p(1,:),maxvals," r ",p(1,:),eff_mean_current," b ");

/*  Polynomial fit */
hold(h," on ");
coeff = polyfit(p(1,:),maxvals,4);
ax = 0:.01:1;
ay = polyval(coeff,ax);
plot(h,ax,ay," r-- ");
legend(" Max activation "," Effective mean current "," Approx max activation ");
pm.done;
/*  save(file,'coeff','-APPEND'); */
file = fullfile(fileparts(mfilename(" fullpath "))," maxforces_poly ");
save(file," coeff ");
}
};
};
};