Compute Discrete Prolate Spheroidal (Slepian) Sequences for use as time-based smoother. This approach uses the tridiagonal method and exploits symmetry. Note the odd order tapers are normalized so that the slope at the centre is positive, in accordance with Slepian (1978) and Thomson (1982). This differs from Percival and Walden (1993). This code follows Chapter 8.3 of Percival and Walden (1993) using LAPACK function calls, Anderson (1999).

1 | ```
.dpss(n,k,nw)
``` |

`n` |
A positive integer, typically the length of the time series-regression data set. |

`k` |
A positive integer, the number of basis vectors for the smoother, often 2*nw - 2. |

`nw` |
A positive double-precision number, the time-bandwidth parameter. The frequency domain analogue of the maximum period of interest. |

`v` |
A 'n' by 'k' matrix of basis vectors of class 'dpss'. Each column is the appropriate Slepian sequence of order 'k-1'. |

Anderson, E. (1999).
*LAPACK Users' guide (Vol. 9).*
SIAM.

Percival, D.B. and Walden, A.T. (1993)
*Spectral analysis for physical applications.*
Cambridge University Press.

Slepian, D. (1978)
Prolate spheroidal wave functions, Fourier analysis, and
uncertainty. V–The discrete case. *Bell System Technical Journal*
Volume **57**, pp. 1371–1430

Thomson, D.J (1982)
Spectrum estimation and harmonic analysis. *Proceedings of the IEEE*
Volume **70**, number 9, pp. 1055–1096.

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.