psalsa.h File Reference

Baseline estimation using Peaked Signal’s Asymmetric Least Squares Algorithm(psalsa). More...

#include <stdexcept>
#include <vector>

Go to the source code of this file.

Functions
const std::vector< double >	sablib::BaselinePsalsa (std::vector< double > &y, const double lambda, const double p, const double k, const unsigned int s=2, const unsigned int loop=10, const double eps=1e-3)
	Performs baseline estimation using Peaked Signal’s Asymmetric Least Squares Algorithm (psalsa).

Detailed Description

Baseline estimation using Peaked Signal’s Asymmetric Least Squares Algorithm(psalsa).

Author

Izadori

References:

• Oller-Moreno, S.; Pardo, A.; Jiménez-Soto, J. M.; Samitier, J.; Marco, S. "Adaptive Asymmetric Least Squares baseline estimation for analytical instruments" 2014 IEEE 11th International Multi-Conference on Systems, Signals & Devices (SSD14), Barcelona, Spain, 2014, pp. 1-5.

Definition in file psalsa.h.

Function Documentation

BaselinePsalsa()

const std::vector< double > sablib::BaselinePsalsa	(	std::vector< double > &	y,
		const double	lambda,
		const double	p,
		const double	k,
		const unsigned int	s = 2,
		const unsigned int	loop = 10,
		const double	eps = 1e-3 )

Performs baseline estimation using Peaked Signal’s Asymmetric Least Squares Algorithm (psalsa).

Parameters

y	The input data for baseline estimation.
lambda	Smoothing parameter.
p	Weight (asymmetry parameter, typically 0.001 to 0.1).
k	Exponential decay of the weights.
s	The order of the difference (usually s = 1, 2, or 3).
loop	Maximum number of iterations.
eps	Convergence threshold.

Returns

The estimated baseline.

Exceptions

std::invalid_argument

One or more parameters are wrong.

Definition at line 16 of file psalsa.cpp.