3 % This is the
default implementation
for the calculation of the gamma value
5 % Andreas Schmidt, 2016
7 function data = gen_detailed_data(rmodel, model_data, detailed_data)
8 % Store some additional data
for the online calculation
9 % Actually, we need nothing!
13 function fn = gen_reduced_data(rmodel, model_data, ...
14 detailed_data, reduced_data)
15 fn = @(rmodel, reduced_data, rbsim) -1;
18 function gamma = gamma(this, model, model_data, dsim)
19 % CALCULATE Calculate the value of gamma through a high
20 % dimensional solution of the
Lyapunov equation
22 [E,A,B,~,~,R] = model.assemble(model_data);
23 Y = A - B*inv(R)*B'*dsim.Z*dsim.Z'*E;
25 gamma = normest(lyap( Y', -eye(model.n), [], E));
LYAPUNOV This is the default implementation for the calculation of the gamma value.
BASE Basis class for the preparation and calculation of gamma Any subclass should also overwrite the...
Implementation of the parametric algebraic Riccati equation.