dpsssave - Discrete prolate spheroidal or Slepian sequence database

Syntax

Description

dpsssave(time_halfbandwith,dps_seq,lambda) creates a database of discrete prolate spheroidal (DPSS) or Slepian sequences and saves the results in dpss.mat. The time half bandwidth producttime_halfbandwith is a real-valued scalar determining the frequency concentration of the Slepian sequences in dps_seq. dps_seq is a NxK matrix of Slepian sequences where N is the length of the sequences. lambda is a 1xK vector containing the frequency concentration ratios of the Slepian sequences in dps_seq.

If the database dpss.mat exists, subsequent calls to dpsssave append the Slepian sequences to the existing file.

status = dpsssave(time_halfbandwith,dps_seq,lambda) returns a 1 if the database operation was successful or a 0 if unsuccessful.

Definitions

Discrete Prolate Spheroidal Sequences

The discrete prolate spheroidal or Slepian sequences derive from the following time-frequency concentration problem. For all finite-energy sequences index limited to some set , which sequence maximizes the following ratio:

where Fs is the sampling frequency . In other words, which index-limited sequence has the largest proportion of its energy in the band [-W,W]. For index-limited sequences, the ratio must satisfy the inequality . The sequence maximizing the ratio is the first discrete prolate spheroidal or Slepian sequence. The second Slepian sequence maximizes the ratio and is orthogonal to the first Slepian sequence. The third Slepian sequence maximizes the ratio of integrals and is orthogonal to both the first and second Slepian sequences. Continuing in this way, the Slepian sequences form an orthogonal set of band limited sequences.

Time Half Bandwidth Product

The time half bandwidth product is NW where N is the length of the sequence and [-W,W] is the effective bandwidth of the sequence. In constructing Slepian sequences, you choose the desired sequence length and bandwidth 2W. Both the sequence length and bandwidth affect how many Slepian sequences have concentration ratios near one. As a rule, there are 2NW-1 Slepian sequences with energy concentration ratios approximately equal to one. Beyond 2NW-1 Slepian sequences, the concentration ratios begin to approach zero. Common choices for the time half bandwidth product are: 2.5, 3, 3.5, and 4.

You can specify the bandwidth of the Slepian sequences in Hz by defining the time half bandwidth product as NW/Fs where Fs is the sampling frequency.

Examples

Create Slepian sequence database in current directory:

seq_length=512;
time_halfbandwidth=2.5;
num_seq=4;
[dps_seq,lambda]=dpss(seq_length,time_halfbandwidth);
% Create databased dpss.mat in current working directory
status=dpsssave(time_halfbandwidth,dps_seq,lambda);
% status should equal 1

References

Percival, D.B., and A.T. Walden. Spectral Analysis for Physical Applications. Cambridge: Cambridge University Press, 1993.

See Also

dpss | dpssclear | dpssdir | dpssload

© 1984-2012- The MathWorks, Inc.    -   Site Help   -   Patents   -   Trademarks   -   Privacy Policy   -   Preventing Piracy   -   RSS