SarcoScaling.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  SarcoScaling() {

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

/*  m = models.motorunit.Shorten(true, true);
 * 
 * m.T = 80;
 * 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);
 * 
 * % Directly set parameter set
 * p = [linspace(0,1,60); linspace(2,3,60)];
 * 
 * n = size(p,2);
 * maxvals = zeros(1,n);
 * maxidx = [];
 * pi = ProcessIndicator('Gathering max forces',n);
 * parfor k=1:n
 *     [t, y, ct, x] = m.simulate(p(:,k),1);%#ok
 *     [maxvals(k), pos] = max(x(59,:));
 *     if pos == m.T+1
 *         fprintf('Max at final time for param %d\n',k);
 *         maxidx = [maxidx k];
 *     end
 *     pi.step;%#ok
 * end
 * pi.stop;
 * save(file,'m','p','maxvals','maxidx'); */

load(file);
pm = PlotManager(false,2,2);
tri = delaunay(p(1,:),p(2,:));
h = pm.nextPlot(" all ");
trisurf(tri, p(1,:),p(2,:),maxvals," Parent ",h, " FaceColor "," interp "," EdgeColor "," interp ");

/*  m.Sampler.Domain = MotoneuronParamDomain; */
[~, idx] = m.Sampler.Domain.filter(p);
/*  tri = delaunay(p(1,idx),p(2,idx));
 * h = pm.nextPlot('valid');
 * trisurf(tri, p(1,idx),p(2,idx), maxvals(idx), 'Parent',h, 'FaceColor','interp','EdgeColor','interp'); */

x = p(1,idx);
y = 1./maxvals(idx);
h = pm.nextPlot(" all ");
plot(h, x, y);
axis(h," tight ");

/*  Polynomial fit */
h = pm.nextPlot(" all ");
plot(h,x,y);
hold(h," on ");
c = polyfit(x,y,12);
ax = 0:.01:1;
ay = polyval(c,ax);
plot(h,ax,ay," r ");
axis(h," tight ");
save(file," c "," -APPEND ");
}
};
};
};