get_l_rough {PEIP}

One-D Roughening

Description

Returns a 1D differentiating matrix operating on a series with n points.

Usage

get_l_rough(n, deg)

Arguments

n

integer, number of data points

deg

order of the derivative to approximate

Details

Used to get first and 2nd order roughening matrices for 1-D problems

Value

Matrix:discrete differentiation matrix

Author(s)

Jonathan M. Lees<[e-mail]>

References

Aster, R.C., C.H. Thurber, and B. Borchers, Parameter Estimation and Inverse Problems, Elsevier Academic Press, Amsterdam, 2005.

Examples

### first order roughening matrix for a 10 by 10 model: a sparse matrix
N = 10
L1 = get_l_rough(10,1);

### second order roughening matrix for a 10 by 10 model
N = 10
L2 = get_l_rough(10,2);