ParamSweep.m

Go to the documentation of this file.

/* (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 [pm , Y , E ] = ParamSweep(rmodel,mu,inputidx,param,paramvals,pm) {


/*  Validity checks  */
if ischar(param)
    pidx = rmodel.getParamIndexFromName(param);
    pname = param;
elseif isposintscalar(param)
    pidx = param;
    pname = rmodel.System.Params(pidx).Name;
else
    error(" param must be either a char array (parameter name) or positive parameter index ");
end
npar = length(paramvals);
if npar == 0
    error(" Parameter values to sweep must be given. ");
end
if nargin < 6
    pm = PlotManager;
end

haveerrest = ~isempty(rmodel.ErrorEstimator) && rmodel.ErrorEstimator.Enabled;

/*  Iterate over parameter values */
Y = zeros(npar,length(rmodel.Times));
E = [];
if haveerrest
    E = Y;
end
pi = ProcessIndicator(" Starting parameter sweep with %d parameter sets for %s... ",npar,false,npar,pname);
for idx = 1:npar
    mu(pidx) = paramvals(idx);
    [t, y] = rmodel.simulate(mu, inputidx);
    if size(y,1) > 1
        y = Norm.L2(y);
    end
    Y(idx,:) = y;

    if haveerrest
        E(idx,:) = rmodel.ErrorEstimator.OutputError;/* #ok */

    end
    pi.step;
end
pi.stop;

/* % Prepare plot */

[T,MU] = meshgrid(t,paramvals);
ustr = ;
if ~isempty(inputidx)
    ustr = sprintf(" , u_%d ",inputidx);
end
tit = sprintf(" Outputs for parameter sweep of '%s' (idx:%d) from %1.3f to %1.3f, base \\mu = [%s]%s ",...
    pname,param,paramvals(1),paramvals(end),num2str(mu^t),ustr);
h = pm.nextPlot(" param_sweep ",tit," t ",sprintf(" Parameter %d: %s ",param,pname));
/* % Upper bound */
if haveerrest
    obj = surf(h,T,MU,Y+E," EdgeColor "," none "," FaceColor "," red ");
    alpha(obj,.3);
/*      mesh(T,MU,Y+E,'EdgeColor','none','FaceColor','red'); */
    hold on;
end

/* % y plot */
surf(h,T,MU,Y," EdgeColor "," none ");
colormap jet;

/* % Lower bound */
if haveerrest
    obj = surf(h,T,MU,Y-E," EdgeColor "," none "," FaceColor "," red ");
    alpha(obj,.3);
/*      mesh(T,MU,Y-E,'EdgeColor','none','FaceColor','red'); */
    hold off;
end
if nargout < 1
    pm.done;
end
}