Numeric.EMD

Empirical Mode Decomposition in pure Haskell.

Main interface is emd, with defaultEO. A tracing version that outputs a log to stdout is also available, as emdTrace. This can be used to help track down a specific IMF that might be taking more time than desired.

This package uses "sized vectors" as its main interface, to ensure:

1. The resulting EMD contains IMFs that are all the same length as the input vector
2. We provide a vector of size of at least one.

There are many functions to convert unsized vectors to sized vectors in Data.Vector.Sized and associated modules, including toSized (for when you know the size at compile-time) and withSized (for when you don't).

Empirical Mode Decomposition

emd :: (Vector v a, KnownNat n, Fractional a, Ord a) => EMDOpts a -> Vector v (n + 1) a -> EMD v (n + 1) a

EMD decomposition of a given time series with a given sifting stop condition.

Takes a sized vector to ensure that:

1. The resulting EMD contains IMFs that are all the same length as the input vector
2. We provide a vector of size of at least one.

emdTrace :: (Vector v a, KnownNat n, Fractional a, Ord a, MonadIO m) => EMDOpts a -> Vector v (n + 1) a -> m (EMD v (n + 1) a)

emd, but tracing results to stdout as IMFs are found. Useful for debugging to see how long each step is taking.

emd' :: (Vector v a, KnownNat n, Fractional a, Ord a, Applicative m) => (SiftResult v (n + 1) a -> m r) -> EMDOpts a -> Vector v (n + 1) a -> m (EMD v (n + 1) a)

emd with a callback for each found IMF.

data EMD v n a

An EMD v n a is an Empirical Mode Decomposition of a time series with n items of type a stored in a vector v.

Constructors

EMD

Fields

• emdIMFs :: ![Vector v n a]
• emdResidual :: !(Vector v n a)

Instances

Show (v a) => Show (EMD v n a)

Defined in Numeric.EMD

Methods

showsPrec :: Int -> EMD v n a -> ShowS

show :: EMD v n a -> String

showList :: [EMD v n a] -> ShowS

data EMDOpts a

Options for EMD composition.

Constructors

EO

Fields

• eoSiftCondition :: SiftCondition a

  stop condition for sifting

• eoSplineEnd :: SplineEnd a

  end conditions for envelope splines

• eoBoundaryHandler :: Maybe BoundaryHandler

  process for handling boundary

Instances

Eq a => Eq (EMDOpts a)

Defined in Numeric.EMD

Methods

(==) :: EMDOpts a -> EMDOpts a -> Bool

(/=) :: EMDOpts a -> EMDOpts a -> Bool

Ord a => Ord (EMDOpts a)

Defined in Numeric.EMD

Methods

compare :: EMDOpts a -> EMDOpts a -> Ordering

(<) :: EMDOpts a -> EMDOpts a -> Bool

(<=) :: EMDOpts a -> EMDOpts a -> Bool

(>) :: EMDOpts a -> EMDOpts a -> Bool

(>=) :: EMDOpts a -> EMDOpts a -> Bool

max :: EMDOpts a -> EMDOpts a -> EMDOpts a

min :: EMDOpts a -> EMDOpts a -> EMDOpts a

Show a => Show (EMDOpts a)

Defined in Numeric.EMD

Methods

showsPrec :: Int -> EMDOpts a -> ShowS

show :: EMDOpts a -> String

showList :: [EMDOpts a] -> ShowS

defaultEO :: Fractional a => EMDOpts a

Default EMDOpts

data BoundaryHandler

Constructors

BHClamp

Clamp envelope at end points (Matlab implementation)

BHSymmetric

Extend boundaries symmetrically

Instances

Eq BoundaryHandler

Defined in Numeric.EMD

Methods

(==) :: BoundaryHandler -> BoundaryHandler -> Bool

(/=) :: BoundaryHandler -> BoundaryHandler -> Bool

Ord BoundaryHandler

Defined in Numeric.EMD

Methods

compare :: BoundaryHandler -> BoundaryHandler -> Ordering

(<) :: BoundaryHandler -> BoundaryHandler -> Bool

(<=) :: BoundaryHandler -> BoundaryHandler -> Bool

(>) :: BoundaryHandler -> BoundaryHandler -> Bool

(>=) :: BoundaryHandler -> BoundaryHandler -> Bool

max :: BoundaryHandler -> BoundaryHandler -> BoundaryHandler

min :: BoundaryHandler -> BoundaryHandler -> BoundaryHandler

Show BoundaryHandler

Defined in Numeric.EMD

Methods

showsPrec :: Int -> BoundaryHandler -> ShowS

show :: BoundaryHandler -> String

showList :: [BoundaryHandler] -> ShowS

data SiftCondition a

Stop conditions for sifting process

Constructors

SCStdDev !a

Stop using standard SD method

SCTimes !Int

Stop after a fixed number of sifting iterations

SCOr (SiftCondition a) (SiftCondition a)

One or the other

SCAnd (SiftCondition a) (SiftCondition a)

Stop when both conditions are met

Instances

Eq a => Eq (SiftCondition a)

Defined in Numeric.EMD

Methods

(==) :: SiftCondition a -> SiftCondition a -> Bool

(/=) :: SiftCondition a -> SiftCondition a -> Bool

Ord a => Ord (SiftCondition a)

Defined in Numeric.EMD

Methods

compare :: SiftCondition a -> SiftCondition a -> Ordering

(<) :: SiftCondition a -> SiftCondition a -> Bool

(<=) :: SiftCondition a -> SiftCondition a -> Bool

(>) :: SiftCondition a -> SiftCondition a -> Bool

(>=) :: SiftCondition a -> SiftCondition a -> Bool

max :: SiftCondition a -> SiftCondition a -> SiftCondition a

min :: SiftCondition a -> SiftCondition a -> SiftCondition a

Show a => Show (SiftCondition a)

Defined in Numeric.EMD

Methods

showsPrec :: Int -> SiftCondition a -> ShowS

show :: SiftCondition a -> String

showList :: [SiftCondition a] -> ShowS

defaultSC :: Fractional a => SiftCondition a

Default SiftCondition

data SplineEnd a

End condition for spline

Constructors

SENotAKnot

SENatural

SEClamped a a

Instances

Eq a => Eq (SplineEnd a)

Defined in Numeric.EMD.Internal.Spline

Methods

(==) :: SplineEnd a -> SplineEnd a -> Bool

(/=) :: SplineEnd a -> SplineEnd a -> Bool

Ord a => Ord (SplineEnd a)

Defined in Numeric.EMD.Internal.Spline

Methods

compare :: SplineEnd a -> SplineEnd a -> Ordering

(<) :: SplineEnd a -> SplineEnd a -> Bool

(<=) :: SplineEnd a -> SplineEnd a -> Bool

(>) :: SplineEnd a -> SplineEnd a -> Bool

(>=) :: SplineEnd a -> SplineEnd a -> Bool

max :: SplineEnd a -> SplineEnd a -> SplineEnd a

min :: SplineEnd a -> SplineEnd a -> SplineEnd a

Show a => Show (SplineEnd a)

Defined in Numeric.EMD.Internal.Spline

Methods

showsPrec :: Int -> SplineEnd a -> ShowS

show :: SplineEnd a -> String

showList :: [SplineEnd a] -> ShowS

Internal

sift :: (Vector v a, KnownNat n, Fractional a, Ord a) => EMDOpts a -> Vector v (n + 1) a -> SiftResult v (n + 1) a

Iterated sifting process, used to produce either an IMF or a residual.

data SiftResult v n a

The result of a sifting operation. Each sift either yields a residual, or a new IMF.

Constructors

SRResidual !(Vector v n a)

SRIMF !(Vector v n a) !Int

number of sifting iterations

envelopes :: (Vector v a, KnownNat n, Fractional a, Ord a) => SplineEnd a -> Maybe BoundaryHandler -> Vector v (n + 1) a -> Maybe (Vector v (n + 1) a, Vector v (n + 1) a)

Returns cubic splines of local minimums and maximums. Returns Nothing if there are not enough local minimum or maximums to create the splines.