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

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