| RMkolmogorov {RandomFields} | R Documentation |
RMkolmogorov corresponds to a vector-valued random field with
covariance function
γ_{ij}(h) = |h|^{2/3}(4/3 * δ_{ij} - 1/3 * h_i h_j / \|h\|^2)
RMkolmogorov(var, scale, Aniso, proj)
var,scale,Aniso,proj |
optional arguments; same meaning for any
|
RMkolmogorov returns an object of class RMmodel.
Martin Schlather, schlather@math.uni-mannheim.de, http://ms.math.uni-mannheim.de
The above formula is from eq. (6.32) of section 6.2 in
Pope, S.B. (2000) Turbulent Flows. Cambridge: Cambridge University Pess.
RMmodel, RMcurlfree,
RMdivfree,
RMvector.
RFoptions(seed=0) ## *ANY* simulation will have the random seed 0; set
## RFoptions(seed=NA) to make them all random again
x <- y <- seq(-2, 2, len=20)
model <- RMkolmogorov()
plot(model, dim=3, MARGIN=1:2, fixed.MARGIN=1)
simu <- RFsimulate(model, x, y, z=0)
plot(simu, select.variables=list(c(1,2)), col=c("red"))