BuildMatrices: Building k-box model matrices

In donaldcummins/EBM: Estimate Parameters of Stochastic Energy Balance Models

Building k-box model matrices

Description

BuildMatrices takes physical parameters of the k-box model differential equations and constructs continuous- and discrete-time model matrices.

Usage

BuildMatrices(gamma, C, kappa, epsilon, sigma_eta, sigma_xi)

Arguments

gamma	stochastic forcing correlation parameter.
C	vector of box heat capacities.
kappa	vector of heat transfer coefficients.
epsilon	deep ocean heat uptake efficacy factor.
sigma_eta	stochastic forcing standard deviation parameter.
sigma_xi	standard deviation of stochastic temperature disturbances.

Value

BuildMatrices returns a list containing k-box model matrices: A, B, Q, Ad, Bd, Qd, Gamma0, Cd.

See Also

SimStepData, FitKalman.

Examples

# set physical parameters
parameters <- list(
  gamma = 2.2,
  C = c(7.0, 80.0),
  kappa = c(1.2, 0.75),
  epsilon = 1.2,
  sigma_eta = 0.54,
  sigma_xi = 0.72
)

# build matrices
matrices <- with(parameters, {
  BuildMatrices(gamma, C, kappa, epsilon, sigma_eta, sigma_xi)
})

# print output
print(matrices)

donaldcummins/EBM documentation built on Oct. 15, 2024, 6:17 a.m.