On the numeric stability of the SFA implementation sfa-tk

12/06/2009
by   Wolfgang Konen, et al.
0

Slow feature analysis (SFA) is a method for extracting slowly varying features from a quickly varying multidimensional signal. An open source Matlab-implementation sfa-tk makes SFA easily useable. We show here that under certain circumstances, namely when the covariance matrix of the nonlinearly expanded data does not have full rank, this implementation runs into numerical instabilities. We propse a modified algorithm based on singular value decomposition (SVD) which is free of those instabilities even in the case where the rank of the matrix is only less than 10 that an alternative way of handling the numerical problems is to inject a small amount of noise into the multidimensional input signal which can restore a rank-deficient covariance matrix to full rank, however at the price of modifying the original data and the need for noise parameter tuning.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
10/07/2013

Singular Value Decomposition of Images from Scanned Photographic Plates

We want to approximate the mxn image A from scanned astronomical photogr...
research
06/07/2017

Fast Eigen Decomposition for Low-Rank Matrix Approximation

In this paper we present an efficient algorithm to compute the eigen dec...
research
02/22/2017

On the Power of Truncated SVD for General High-rank Matrix Estimation Problems

We show that given an estimate A that is close to a general high-rank po...
research
05/07/2015

Filter characteristics in image decomposition with singular spectrum analysis

Singular spectrum analysis is developed as a nonparametric spectral deco...
research
10/01/2017

Efficient Preconditioning for Noisy Separable NMFs by Successive Projection Based Low-Rank Approximations

The successive projection algorithm (SPA) can quickly solve a nonnegativ...
research
02/20/2020

Efficiently updating a covariance matrix and its LDL decomposition

Equations are presented which efficiently update or downdate the covaria...
research
06/02/2018

Fast Rigid 3D Registration Solution: A Simple Method Free of SVD and Eigen-Decomposition

A novel solution is obtained to solve the rigid 3D registration problem,...

Please sign up or login with your details

Forgot password? Click here to reset