55 this.
Name=
" MathMOD 2012 Experiments ";
65 s = solvers.ExplEuler;
71 this.
System= models.mathmod2012.MathMODSystem(
this);
77 this.
Sampler= sampling.GridSampler;
79 e1 = [ones(dims,1) zeros(dims,1)];
80 e2 = [zeros(dims,1) ones(dims,1)];
82 x0 = dscomponents.AffineInitialValue;
83 x0.addMatrix(
" mu(3) ", ones(dims,1));
88 f.Ma= repmat(-exp(-f.Centers.xi(1,:)/15),dims,1);
91 this.
System.
Inputs[2] = @(
t)[.5*sin(
t/2); 4*exp(-7*(12-
t).^2)-.5*exp(-(5-
t).^2)];
92 B = dscomponents.AffLinInputConv;
93 B.addMatrix(
" mu(1) ", e1);
94 B.addMatrix(
" 1-mu(1) ", e2);
102 V = ones(dims,1)/sqrt(dims);
103 s = spacereduction.ManualReduction(V,V);
110 s = spacereduction.RotationDecorator(s);
122 m = models.mathmod2012.MathMODExperiment(dim);
127 d.EstimatorVersions= [1 1 0 0 1 0 0 1 1];
128 d.SingleFigures=
true;
134 d.start([mu1; 1], 1);
137 md = d.ModelData(end);
138 absmax = 30*abs(md.ErrT(7)-md.ErrT(1));
139 relmax = 30*abs(md.RelErrT(7)-md.RelErrT(1));
143 axis(a,[0 m.T md.MinErr*.9 min(1e4,absmax)]);
144 set(legend(a),
" Location ",
" NorthEast ");
145 fi = fullfile(
KerMor.
App.DataDirectory,
" MathMOD ",sprintf(
" MathMOD_mu1_%d_deg%f_ ", mu1, s.Degree));
150 f = d.Figures[2]; a = gca(f);
151 axis(a,[0 m.T md.MinRelErr*.9 min(1,relmax)]);
152 set(legend(a),
" Location ",
" NorthEast ");
157 f = d.Figures[3]; a = gca(f);
164 PlotParamSweep(r,[1; 1],1,1,-.1:.05:1.1);
char Name
The name of the Model.
Collection of generally useful functions.
The base class for any KerMor detailed model.
double dt
The desired time-stepsize for simulations.
integer TrainingInputs
The indices of inputs to use for training data generation. Uses the DefaultInput if not set (and Defa...
sampling.BaseSampler Sampler
The sampling strategy the Model uses.
models.BaseFirstOrderSystem System
The actual dynamical system used in the model.
dscomponents.AInputConv B
The input conversion.
function registerProps(varargin)
Call this method at any class that defines DPCM observed properties.
rowvec EstimatorIterations
How many estimator iterations should be performed?
dscomponents.AInitialValue x0
Function handle to initial state evaluation.
solvers.BaseSolver ODESolver
The solver to use for the ODE. Must be an instance of any solvers.BaseSolver subclass.
static function saveFigure(fig, filename, ext)
Opens a matlab save dialog and saves the given figure to the file selected.
Inputs
The system's possible input functions. A cell array of function handles, each taking a time argument ...
approx.BaseApprox Approx
The approximation method for the CoreFunction.
double T
The final timestep up to which to simulate.
Analysis class for the error estimators.
spacereduction.BaseSpaceReducer SpaceReducer
The reduction algorithm for subspaces.
dscomponents.LinearOutputConv C
The output conversion Defaults to an LinearOutputConv instance using a 1-matrix, which just forwards ...
dscomponents.ACoreFun f
The core f function from the dynamical system.
Numerical experiments class for Paper ICIAM.
Global configuration class for all KerMor run-time settings.
static function KerMor theinstance = App()
The singleton KerMor instance.
double MaxTimestep
The maximum timestep allowed for any ODE solvers.
dim
The system's dimension.
static function [ d , r , m ] = CreatePlots(dim)