lorentz: Lorentz-function

In SvenjaFischer/Disagg: Disaggregation of daily mean discharges to hourly discharges including the peak

Description

The Lorentz-function given by Q_i(t)=α_{1_i}+2\frac{α_{2_i}}{π}\frac{α_{3_i}}{4(t-α_{4_i})^2+α_{3_i}^2}

with parameters α_{j_i}

Usage

lorentz(t, par = runif(4, 0, 10))

Arguments

t	function argument, here: time step
par	parameters α_{j_i}, with starting values

Value

returns Q_i(t)

Author(s)

Svenja Fischer

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (t, par = runif(4, 0, 10))
{
    par[1] + 2 * par[2]/pi * par[3]/(4 * (t - par[4])^2 + par[3]^2)
  }

SvenjaFischer/Disagg documentation built on May 14, 2019, 10:37 a.m.