KerMor
0.9
Model order reduction for nonlinear dynamical systems and nonlinear approximation
|
Shorten: Model for a muscle motor unit composed of motoneuron and a sarcomere. More...
Shorten: Model for a muscle motor unit composed of motoneuron and a sarcomere.
The global time unit for this model is milliseconds [ms]. This model has been copied from the MuscleFibreModel in order to obtain a current snapshot without changing the multi-sarcomere adoptions done for the whole fibre model.
This model is intended for research regarding the kernel-approximation of the motor unit w.r.t. the activation and fibre type relation to the force development.
This class is part of the framework
Public Member Functions | |
Shorten (dynamic_ic,logical singlepeakmode, outputscaling) | |
Creates a new motor unit model. More... | |
function pm = | plotMotoSacroLinkFactorCurve () |
function | plotOutputForceScaling (colvec< double > x) |
function pm = | plotState (double t,colvec< double > x, pm) |
function pm = | plot (double t,matrix< double > y, pm) |
Public Member Functions inherited from models.BaseFullModel | |
BaseFullModel (char name) | |
Creates a new instance of a full model. More... | |
function | delete () |
function | off1_createParamSamples () |
Offline phase 1: Sample generation. More... | |
function | off2_genTrainingData () |
Offline phase 2: Snapshot generation for subspace computation. More... | |
function | off3_computeReducedSpace () |
Offline phase 3: Generate state space reduction. More... | |
function | off4_genApproximationTrainData () |
Generates the training data \(x_i\) for the \(\hat{f}\)-approximation and precomputes \(f(x_i)\) values. More... | |
function | off5_computeApproximation () |
Offline phase 5: Core function approximation. More... | |
function | off6_prepareErrorEstimator () |
Prepares offline data for a possible error estimator. More... | |
function | offlineGenerations () |
Performs all large offline computations for model reduction. More... | |
function [
models.ReducedModel reduced , double time ] = | buildReducedModel (varargin) |
Builds a reduced model from a full model. More... | |
function [
t , x , time , cache ] = | computeTrajectory (colvec< double > mu, inputidx) |
Computes a solution/trajectory for the given mu and inputidx in the SCALED state space. More... | |
function [
double t ,
matrix< double > y , colvec < double > mu , in , ct ] = | getSimCacheTrajectory (nr) |
function | plotTrajectoryNr (nr) |
function [
V , W ] = | assembleProjectionMatrices (target_dim) |
function file = | save (char directory,char filename) |
Saves this instance inside the data.ModelData.DataDirectory folder using the model's SaveTag if set or "model.mat". More... | |
Public Member Functions inherited from models.BaseModel | |
BaseModel () | |
function | delete () |
function | initDefaultParameter () |
Reads the default values of the System's ModelParam list and initializes the BaseModel.DefaultMu with it. More... | |
function [
rowvec< double > t , matrix < double > y , double sec , x ] = | simulate (colvec< double > mu,integer inputidx) |
Simulates the system and produces the system's output. More... | |
function [ handle f , handle ax ] = | plot (rowvec< double > t,matrix< double > y, varargin) |
Plots the results of the simulation. Override in subclasses for a more specific plot if desired. More... | |
function [ handle f , handle ax ] = | plotState (rowvec t,matrix< double > x, varargin) |
Plots the results of the simulation. Override in subclasses for a more specific plot if desired. More... | |
function [
f , ax ] = | plotSingle (double t,colvec< double > y, varargin) |
Plots a single solution. Override in subclasses for specific plot behaviour. More... | |
function | plotInputs (pm) |
function [
double t ,
colvec< double > x , double ctime ] = | solveStatic (colvec< double > mu,integer inputidx) |
Solves the linear system \(A(t,\mu)*x + f(t,\mu) + B(t,\mu)*u(t) = 0\). Spatial dependence of f is neglected in the current implementation (since this would require solving a nonlinear system). More... | |
function [
rowvec< double > t , matrix < double > x , double ctime ] = | computeTrajectory (colvec< double > mu,integer inputidx) |
Computes a solution/trajectory for the given mu and inputidx in the SCALED state space. More... | |
function matrix < double > mu = | getRandomParam (integer num,integer seed) |
Gets a random parameter sample from the system's parameter domain P. More... | |
Public Member Functions inherited from KerMorObject | |
KerMorObject () | |
Constructs a new KerMor object. More... | |
function | display () |
disp(object2str(this)); More... | |
function bool = | eq (B) |
Checks equality of two KerMor objects. More... | |
function bool = | ne (B) |
Checks if two KerMorObjects are different. More... | |
function cn = | getClassName () |
Returns the simple class name of this object without packages. More... | |
Public Member Functions inherited from DPCMObject | |
DPCMObject () | |
Creates a new DPCM object. More... | |
DPCMObject () | |
Static Public Member Functions | |
static function res = | test_Shorten () |
Static Public Member Functions inherited from models.BaseFullModel | |
static function | test_BaseModels () |
static function res = | test_BareModel () |
static function | test_LinearModel () |
static function | test_LinearModelNoProj () |
static function | test_LinearModelNoApprox () |
static function | test_LinearModelParams () |
static function | test_LinearModelInputs () |
static function | test_LinearModelParamsInput () |
static function | test_NonlinearModel () |
static function | test_NonlinearModelParams () |
static function | test_NonlinearModelInputs () |
static function | test_NonlinearModelParamsInput () |
static function | test_TimeDependentOutput () |
static function res = | test_FileSystemTidyness () |
Plain small model. More... | |
Public Attributes | |
UseNoise | |
logical | SinglePeakMode = false |
Flag that determines if this system is to be run so that only ever one peak/signal will be issued. More... | |
logical | DynamicInitialConditions = true |
Set this flag to true (in constructor) if you want parameter-dependent initial conditions that have been obtained using long-time simulations of the model for different fibre types and then using the end state as "stable" initial condition. More... | |
logical | SingleTwitchOutputForceScaling = true |
Flag to determine if the output force (e.g. calcium concentration) should be re-scaled so that the forces using a single twitch are the same maximal forces (=1) for each fibre type. More... | |
Public Attributes inherited from models.BaseFullModel | |
data.ModelData | Data |
The full model's data container. Defaults to an empty container. More... | |
export.JKerMorExport | JKerMorExport |
Export instance for possible export of this model to JKerMor. More... | |
sampling.BaseSampler | Sampler |
The sampling strategy the Model uses. More... | |
spacereduction.BaseSpaceReducer | SpaceReducer |
The reduction algorithm for subspaces. More... | |
approx.BaseApprox | Approx = "[]" |
The approximation method for the CoreFunction. More... | |
function_handle | preApproximationTrainingCallback |
Advanced property. Must be a function handle taking the current model instance. More... | |
function_handle | postApproximationTrainingCallback |
Advanced property. Must be a function handle taking the current model instance. More... | |
logical | EnableTrajectoryCaching = false |
Flag that enables caching of computed trajectories in a simulation cache stored in KerMor's TempDirectory folder. More... | |
error.BaseEstimator | ErrorEstimator = "[]" |
The associated error estimator for this model. More... | |
logical | ComputeTrajectoryFxiData = false |
Flag that determines whether fxi data should be computed along with the trajectories or not. More... | |
logical | ProjectApproxTrainingData = false |
Flag to determine if the data collected by the approx.Approx.TrainDataSelector should be projected into the previously computed projection space V (see offline phase 3). If so, the \(f(x_i)\) values are also re-computed for the updated \(x_i\) samples. More... | |
logical | AutoSave = false |
Flag to enable automatic saving of the model after each individual offline phase step and at other locations prone to data loss due to lengthy computations. More... | |
char | SaveTag = "" |
A custom tag that can be used as a prefix to files for corresponding model identification. More... | |
integer | TrainingInputs |
The indices of inputs to use for training data generation. Uses the DefaultInput if not set (and DefaultInput is not empty) More... | |
rowvec< integer > | TrainingParams |
The indices of the model parameters to use for training data generation. More... | |
integer | TrainingInputCount |
Gets the number of inputs used for training. More... | |
rowvec< double > | OfflinePhaseTimes = "[]" |
The computation times for all phases of the last offline generation. More... | |
TrajectoriesCompleted = "[]" | |
Contains a flag for each input/parameter combination which has been successfully computed completely during the offline phase. More... | |
Public Attributes inherited from models.BaseModel | |
models.BaseFirstOrderSystem | System = "[]" |
The actual dynamical system used in the model. More... | |
char | Name = "" |
The name of the Model. More... | |
matrix< double > | G = 1 |
The custom scalar product matrix \(\vG\). More... | |
double | RealTimePlottingMinPause = .1 |
Minimum pause between successive steps when RealTimePlotting is enabled. More... | |
logical | isStatic = false |
Determines if the model is time dependent or static. More... | |
integer | DefaultInput = "[]" |
The default input to use if none is given. More... | |
Times | |
Evaluation points \(\{0=t_0,\ldots,t_n=T\}\) of the model. More... | |
rowvec< double > | scaledTimes |
The time steps Times in scaled time units \(\tilde{t_i} = \frac{t_i}{\tau}\). More... | |
double | Tscaled |
The scaled end time \(\tilde{T} = \frac{T}{\tau}\). More... | |
double | tau |
Time scaling factor \(\tau\). More... | |
double | T |
The final timestep \(T\) up to which to simulate. More... | |
double | dt |
The desired time-stepsize \(\Delta t\) for simulations. More... | |
solvers.BaseSolver | ODESolver |
The solver to use for the ODE. Must be an instance of any solvers.BaseSolver subclass. More... | |
logical | RealTimePlotting |
Determines if the simulation should plot intermediate steps during computation. More... | |
colvec< double > | DefaultMu = "[]" |
The default parameter value if none is given. More... | |
dtscaled = .1 | |
The scaled timestep \(\tilde{\Delta t} = \frac{\Delta t}{\tau}\). More... | |
char | gitRefOnSave = "" |
Contains the GIT revision of this model when it was last saved. More... | |
Public Attributes inherited from DPCMObject | |
WorkspaceVariableName = "" | |
The workspace variable name of this class. Optional. More... | |
ID = "[]" | |
An ID that allows to uniquely identify this DPCMObject (at least within the current MatLab session/context). More... | |
PropertiesChanged = "[]" | |
The Dictionary containing all the property settings as key/value pairs. More... | |
Public Attributes inherited from handle | |
addlistener | |
Creates a listener for the specified event and assigns a callback function to execute when the event occurs. More... | |
notify | |
Broadcast a notice that a specific event is occurring on a specified handle object or array of handle objects. More... | |
delete | |
Handle object destructor method that is called when the object's lifecycle ends. More... | |
disp | |
Handle object disp method which is called by the display method. See the MATLAB disp function. More... | |
display | |
Handle object display method called when MATLAB software interprets an expression returning a handle object that is not terminated by a semicolon. See the MATLAB display function. More... | |
findobj | |
Finds objects matching the specified conditions from the input array of handle objects. More... | |
findprop | |
Returns a meta.property objects associated with the specified property name. More... | |
fields | |
Returns a cell array of string containing the names of public properties. More... | |
fieldnames | |
Returns a cell array of string containing the names of public properties. See the MATLAB fieldnames function. More... | |
isvalid | |
Returns a logical array in which elements are true if the corresponding elements in the input array are valid handles. This method is Sealed so you cannot override it in a handle subclass. More... | |
eq | |
Relational functions example. See details for more information. More... | |
transpose | |
Transposes the elements of the handle object array. More... | |
permute | |
Rearranges the dimensions of the handle object array. See the MATLAB permute function. More... | |
reshape | |
hanges the dimensions of the handle object array to the specified dimensions. See the MATLAB reshape function. More... | |
sort | |
ort the handle objects in any array in ascending or descending order. More... | |
Public Attributes inherited from IParallelizable | |
ComputeParallel = false | |
Flag whether the code should be run in parallel or not. More... | |
Static Protected Member Functions | |
static function this = | loadobj () |
Static Protected Member Functions inherited from models.BaseFullModel | |
static function this = | loadobj (this, sobj) |
Static Protected Member Functions inherited from models.BaseModel | |
static function this = | loadobj (this, s) |
Static Protected Member Functions inherited from DPCMObject | |
static function obj = | loadobj (obj, from) |
Re-register any registered change listeners! More... | |
static function obj = | loadobj (obj, from) |
Additional Inherited Members | |
Protected Member Functions inherited from models.BaseModel | |
function this = | saveobj () |
Store the current GIT branch in the object. More... | |
function | plotstep (src, ed) |
Callback for the ODE solvers StepPerformed event that enables during-simulation-plotting. More... | |
Protected Member Functions inherited from KerMorObject | |
function | checkType (obj, type) |
Object typechecker. More... | |
Protected Member Functions inherited from DPCMObject | |
function | registerProps (varargin) |
Call this method at any class that defines DPCM observed properties. More... | |
function | registerProps (varargin) |
models.motorunit.Shorten.Shorten | ( | dynamic_ic, | |
logical | singlepeakmode, | ||
outputscaling | |||
) |
Creates a new motor unit model.
singlepeakmode | A flag that determines if this model should only ever produce one force peak Default: false |
Definition at line 109 of file Shorten.m.
References models.BaseFullModel.Data, models.BaseModel.DefaultInput, models.BaseModel.DefaultMu, models.BaseModel.dt, DynamicInitialConditions, models.BaseFullModel.EnableTrajectoryCaching, data.ModelData.ModelData(), models.BaseModel.Name, models.BaseModel.ODESolver, models.BaseFullModel.Sampler, models.BaseFullModel.SaveTag, SinglePeakMode, SingleTwitchOutputForceScaling, models.BaseModel.System, models.BaseModel.T, models.BaseFullModel.TrainingInputs, and data.ModelData.useFileTrajectoryData().
|
staticprotected |
Definition at line 214 of file Shorten.m.
References models.BaseFirstOrderSystem.mu, and models.BaseModel.System.
Referenced by plotMotoSacroLinkFactorCurve(), plotOutputForceScaling(), and plotState().
function pm = models.motorunit.Shorten.plotMotoSacroLinkFactorCurve | ( | ) |
Definition at line 171 of file Shorten.m.
References models.BaseFirstOrderSystem.f, PlotManager.LeaveOpen, plot(), and models.BaseModel.System.
Definition at line 185 of file Shorten.m.
References plot(), polyval(), and models.BaseModel.System.
Definition at line 196 of file Shorten.m.
References models.BaseFirstOrderSystem.mu, plot(), and models.BaseModel.System.
|
static |
models.motorunit.Shorten.DynamicInitialConditions = true |
Set this flag to true (in constructor) if you want parameter-dependent initial conditions that have been obtained using long-time simulations of the model for different fibre types and then using the end state as "stable" initial condition.
The actual values have been computed using a polynomial fit of degree seven on each dimension for the parameter range [0,1]
Default: true
SetAccess = Private, GetAccess = Public
Definition at line 68 of file Shorten.m.
Referenced by Shorten().
models.motorunit.Shorten.SinglePeakMode = false |
Flag that determines if this system is to be run so that only ever one peak/signal will be issued.
Default: false
SetAccess = Private, GetAccess = Public
Definition at line 54 of file Shorten.m.
Referenced by Shorten().
models.motorunit.Shorten.SingleTwitchOutputForceScaling = true |
Flag to determine if the output force (e.g. calcium concentration) should be re-scaled so that the forces using a single twitch are the same maximal forces (=1) for each fibre type.
Default: true
SetAccess = Private, GetAccess = Public
Definition at line 89 of file Shorten.m.
Referenced by Shorten().