spectrum.music

Multiple signal classification spectrum

Syntax

Hs = spectrum.music
Hs = spectrum.music(NSinusoids)
Hs = spectrum.music(NSinusoids,SegmentLength)
Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent)
Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName)
Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName,SubspaceThreshold)
Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName,SubspaceThreshold,InputType)


Description

Note

The use of spectrum.music is not recommended. Use pmusic instead.

Hs = spectrum.music returns a default multiple signal classification (MUSIC) spectrum object, Hs, that defines the parameters for the MUSIC spectral estimation algorithm, which uses Schmidt's eigenspace analysis algorithm. This object uses the following default values.

Default Values

Property Name	Default Value	Description
NSinusoids	2	Number of complex sinusoids
SegmentLength	4	Length of each of the time-based segments into which the input signal is divided.
OverlapPercent	50	Percent overlap between segments
WindowName	'Rectangular'	Window name or 'User Defined' (see window for valid window names). For more information on each window, refer to its reference page). This argument can also be a cell array containing the window name or 'User Defined' and, if used for the particular window, an optional parameter value. The syntax is {wname,wparam}. You can use set to change the value of the additional parameter or to define the MATLAB® expression and parameters for a user-defined window (see spectrum for information on using set).
SubspaceThreshold	0	Threshold is the cutoff for signal and noise separation. The threshold is multiplied by λmin , the smallest estimated eigenvalue of the signal's correlation matrix. Eigenvalues below the threshold (λmin*threshold) are assigned to the noise subspace.
InputType	'Vector'	Type of input that will be used with this spectrum object. Valid values are 'Vector', 'DataMatrix' and 'CorrelationMatrix'.

Hs = spectrum.music(NSinusoids) returns a spectrum object, Hs, with the specified number of sinusoids and default values for all other properties. Refer to the table above for default values.

Hs = spectrum.music(NSinusoids,SegmentLength) returns a spectrum object, Hs, with the specified segment length.

Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent) returns a spectrum object, Hs, with the specified overlap between segments.

Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName) returns a spectrum object, Hs, with the specified window.

Note

Window names must be enclosed in single quotes, such as spectrum.music(3,32,50,'chebyshev') or spectrum.music(3,32,50,{'chebyshev',60})

Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName,SubspaceThreshold) returns a spectrum object, Hs, with the specified subspace threshold.

Hs = spectrum.music(NSinusoids,SegmentLength,... OverlapPercent,WindowName,SubspaceThreshold,InputType) returns a spectrum object, Hs, with the specified input type.

Note

See pmusic for more information on the MUSIC algorithm.

Examples

MUSIC Pseudospectrum of a Sinusoidal Signal

Define a complex signal with three sinusoids, add noise, and estimate its pseudospectrum using the MUSIC algorithm.

n = 0:99;
s = exp(1i*pi/2*n) + 2*exp(1i*pi/4*n) + exp(1i*pi/3*n) + randn(1,100);

Hs = spectrum.music(3,20);

pseudospectrum(Hs,s)

References

[1] Harris, Fredric. J. “On the Use of Windows for Harmonic Analysis with the Discrete Fourier Transform.” Proceedings of the IEEE^®. Vol. 66, January 1978, pp. 51–83.

Version History

Introduced before R2006a