
DPSSX.OMS
Script File:
# Description:
# Computes and plots the first three discrete prolate sequences.
#
# The sequences plotted in this example appear in Fig. 2 of the article
# "Multitaper spectral analysis of high  frequency seismograms"
# by J. Park (JGR vol. 92 no. B12 p. 12675  12684
#
clear
# number of points in the sequences
N = 128
# number of sequences
m = 3
# bandwidth of the sequences
w = 2e0 / N
# compute the sequences
S = dpss(N, m, w)
# two line title for plot
title = { ...
"Prolate Spheroidal Windows", ...
" N = 128, W = 4. / 128" ...
}
# place the title at the top of the plot
gaddtext(title, [.5, 1.], [.5, 1.])
# plot sequences with x axis between 0 and 1
gplot(seq(N) / real(N), S)
Output:

