-----------------------------------------------------------------------------
-- |
-- Module      :  DSP.Estimation.Spectral.KayData
-- Copyright   :  (c) Matthew Donadio 2002
-- License     :  GPL
--
-- Maintainer  :  m.p.donadio@ieee.org
-- Stability   :  experimental
-- Portability :  portable
--
-- Test vectors from Kay, /Modern Spectral Estimation/
--
-----------------------------------------------------------------------------

module DSP.Estimation.Spectral.KayData (xc,xr) where

import Data.Array (Array, array)
import Data.Complex (Complex((:+)))

-- | Complex test data

xc :: Array Int (Complex Double)
xc = array (0,31) [
   (0,  (( 6.3307)    :+ (-0.174915))),
   (1,  ((-1.33539)   :+ (-0.03044))),
   (2,  (( 3.61896)   :+ (-0.260459))),
   (3,  (( 1.87513)   :+ (-0.323974))),
   (4,  ((-1.08561)   :+ (-0.136055))),
   (5,  (( 3.99114)   :+ (-0.101864))),
   (6,  ((-4.10184)   :+ ( 0.130571))),
   (7,  (( 1.55399)   :+ ( 0.0977916))),
   (8,  ((-2.1258)    :+ (-0.306485))),
   (9,  ((-3.27873)   :+ (-0.0544436))),
   (10, (( 0.241218)  :+ ( 0.0962379))),
   (11, ((-5.74708 )  :+ ( 0.0186908))),
   (12, ((-0.0165977) :+ ( 0.237493))),
   (13, ((-3.28921)   :+ (-0.188478))),
   (14, ((-1.31227)   :+ (-0.120636))),
   (15, (( 0.745251)  :+ (-0.0679575))),
   (16, ((-1.77199)   :+ (-0.416229))),
   (17, (( 2.56419)   :+ (-0.270373))),
   (18, (( 0.21325)   :+ (-0.232544))),
   (19, (( 2.23409)   :+ ( 0.236383))),
   (20, (( 2.2949)    :+ ( 0.173061))),
   (21, (( 1.09186)   :+ ( 0.140938))),
   (22, (( 2.29353)   :+ ( 0.442044))),
   (23, (( 0.695823)  :+ ( 0.509325))),
   (24, (( 0.759858)  :+ ( 0.417967))),
   (25, ((-0.354267)  :+ ( 0.506891))),
   (26, ((-0.594517)  :+ ( 0.39708))),
   (27, ((-1.88618)   :+ ( 0.649179))),
   (28, ((-1.39041)   :+ ( 0.867086))),
   (29, ((-3.06381)   :+ ( 0.422965))),
   (30, ((-2.0433)    :+ ( 0.0825514))),
   (31, ((-2.1628)    :+ (-0.0933218))) ]
-- | Real test data

xr :: Array Int Double
xr = array (0,31) [
   (0,   6.46768),
   (1,  -1.28024),
   (2,   3.74788),
   (3,   1.96092),
   (4,  -0.768349),
   (5,   4.14569),
   (6,  -4.05277),
   (7,   1.65836),
   (8,  -2.06405),
   (9,  -3.33397),
   (10,  0.085145),
   (11, -6.06562),
   (12, -0.411658),
   (13, -3.61831),
   (14, -1.53352),
   (15,  0.481522),
   (16, -1.93653),
   (17,  2.35532),
   (18,  0.145624),
   (19,  2.21991),
   (20,  2.25884),
   (21,  1.07373),
   (22,  2.26531),
   (23,  0.685007),
   (24,  0.762859),
   (25, -0.501008),
   (26, -0.640518),
   (27, -1.99263),
   (28, -1.60416),
   (29, -3.22751),
   (30, -2.21946),
   (31, -2.42246) ]