get_l_rough: One-D Roughening

In PEIP: Geophysical Inverse Theory and Optimization

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);

PEIP documentation built on Aug. 21, 2023, 9:10 a.m.