Function File: irsa_act.m

[rxp,ryp,err] = irsa_act (xp, yp, fmax, [Nout], [itmax],[thresh])

Reconstructs a band-limited signal from an irregular sampled timeseries using Adaptive weights, Conjugate gradients and Toeplitz matrices.

Input:

xp : Vector – sampling points

yp : Vector – corresponding sampling values

fmax : Scalar – maximum (band-limiting) frequency

Nout : Vector – Number of sampling points for the representation of the reconstruced signal. Default is the shortest possible representation.

itmax : Scalar – maximum number of reconstruction iterations. Default is twice the number of the frequencies.

thresh: Scalar – threshold for conjugate gradients (stop criterion). Default is machine precision.

Output:

rfyp : Columnvector – reconstructed spectrum

err : Columnvector – error for each iteration step