




Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Main points of this file are Create, Cumulative Sum, Differences of Order, Running Median of Length, The T4253H Smoothing Function, Fast Fourier Transform, Inverse Fast Fourier Transform of Two Series
Typology: Study notes
1 / 8
This page cannot be seen from the preview
Don't miss anything!





1
CREATE produces new series as a function of existing series.
The following notation is used throughout this chapter unless otherwise stated: Existing Series X 1 , K,X (^) n New Series Y 1 , K,Yn
Y (^) j X (^) i j n i
j = = =
โ 1
Define
Z (^) j 1 6k = Z (^) j 1 k โ 16 โ Z (^) jโ 11 k โ 16 k = 1 , K, m j = k + 1 , K,n
with
Z (^) j 1 6 0 = X (^) j j = 1 , K,n
then
Y (^) j = %&Z^ jm^ j^ =^ m^ + n '
1 6 1,^ K, SYSMIS otherwise
X j m n j (^) j m = j^ m =^ + =
% & '
Y X^ j^ n^ m j (^) j n m n = j^ m =^ โ = โ +
% & '
If m is odd, define
q = mโ^1 2
then
Y (^) j X^ j^ k m^ j^ q^ n^ q k q
q = =^ +^ โ
% &
KK
'
KK
=โ
โ 1,^ K, SYSMIS otherwise
If m is even, define q = m 2 and
Z (^) j X (^) j k m j q n q k q
q = (^) + = โ =โ +
โ 1
then
Y (^) j = Z^ j^ +^ Z^ j j^ =^ q^ +^ n^ โq
% &
K 'K^
(^3) โ 1 8 2 1 ,^ K, SYSMIS otherwise
where
Z (^) j 1 6 0 = X (^) j j = 1 , K,n
then
Y (^) j = Z (^) j1 6m j = mp +1, K,n
The original series is smoothed by a compound data smoother based on Velleman (1980). The smoother starts with:
Z (^) j + 1 2 = median 3 X (^) j โ 1 , X (^) j , X (^) j + 1 , X (^) j+ 2 8 j = 2 , K,nโ 2
and
Z (^) n X (^) n X (^) n X (^) n X (^) n Z (^) n Xn
0 5
1 1 1 5^1 1 2 12 1
(^1) 1 2 1 12 1 1 1 2
.
( ) ( ).^ median^ ,
median ,
1 6 1 6
0 5 1 6 1 6 0 5
Z 1 1 6^1 = Z (^) 0 5. Z (^) n1 6^1 = Zn+1 2
and
Z 1 6j 1 = 21 4 Z (^) j โ1 2 + Z (^) j+1 2 9 j = 2 , K, nโ 1.
Let Z 1 6^2 be the resulting series, then
n n
n n n n
1
2 1
(^1 2 )
2
2 1
1 2
1 3
1
1
2 2
1 1
1 1
0 5 0 5 0 5 0 5
0 5 0 5 0 5 0 5
0 5 0 5 0 5 0 5
4 9
4 9
median , ,
median , ,
and
Z 1 6j 2 = median Z 1 6j 1 โ 2 , Z 1 6j 1 โ^1 , Z 1 6j 1 , Z 1 6j 1 + 1 , Z 1 6j^1 +^2 j = 3 , K,nโ 2
Let Z 1 6^3 be the resulting series, then
Z Z Z Z j n
j j j j
n n n n n
3 1
2 2 1
2
13 23 33 12 23
(^3313 )
0 5 0 5 0 5 0 5
0 5 0 5 0 5 0 5 0 5
0 5 0 5 0 5 0 5 0 5
4 9
4 9
4 9
โ +
โ โ โ
median , , , , ,
median , ,
median , ,
Z Z Z Z j n
j j j j
n n
(^4 ) 4 1
(^3 ) 2
(^3 ) 4 1
3
14 13 4 3
0 5 0 5 0 5 0 5 2 1
0 5 0 5 0 5 0 5
Thus a, b are two sequences generated by FFT and they are called real and imaginary, respectively.
a n
X f t k r
b n
X f t k r
k t k t
n
k t k t
n
=
=
1
1
cos , ,
sin , ,
ฯ
ฯ
1 0 56
1 0 56
where
r n^ n n n
= %& โ '
1 6 if^ is odd if is even
and
a X b n
X (^) t t t
n 0 0 1
=
โ cos^2 ฯ^1
The inverse Fourier Transform of two series {a, b} is defined as
X (^) t a b t a (^) k f (^) k t b (^) k f (^) kt k
q
k
q = โ โ + โ โ โ
!
"
$
(^0 0) โ= 1 โ= 1 # cos 2 ฯ 1 16 7 2 cos 2 2 ฯ 1 16 7 sin 22 ฯ 1 167
Brigham, E. O. 1974. The fast fourier transform. Englewood Cliffs, N.J.: Prentice- Hall.
Monro, D. M. 1975. Algorithm AS 83: Complex discrete fast Fourier transform. Applied Statistics, 24: 153โ160.
Monro, D. M., and Branch, J. L. 1977. Algorithm AS 117: The chirp discrete Fourier transform of general length. Applied Statistics, 26: 351โ361.
Velleman, P. F. 1980 Definition and comparison of robust nonlinear data smoothing algorithms. Journal of the American Statistical Associations, 75: 609โ615.
Velleman, P. F., and Hoaglin, D. C. 1981. Applications, basics, and computing of exploratory data analysis. Boston: Duxbury Press.