This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.


System object: phased.BeamscanEstimator
Package: phased

Plot spatial spectrum


h = plotSpectrum(___)


plotSpectrum(H) plots the spatial spectrum resulting from the last call of the step method.

plotSpectrum(H,Name,Value) plots the spatial spectrum with additional options specified by one or more Name,Value pair arguments.

h = plotSpectrum(___) returns the line handle in the figure.

Input Arguments


Spatial spectrum estimator object.

Name-Value Pair Arguments

Specify optional comma-separated pairs of Name,Value arguments. Name is the argument name and Value is the corresponding value. Name must appear inside single quotes (' '). You can specify several name and value pair arguments in any order as Name1,Value1,...,NameN,ValueN.


Set this value to true to plot the normalized spectrum. Setting this value to false plots the spectrum without normalization.

Default: false


Character vector to use as figure title.

Default: ''


Plot units, specified as 'db', 'mag', or 'pow'.

Default: 'db'


expand all

Estimate the DOA's of two signals received by a 10-element ULA with element spacing of 1 meter. The antenna operating frequency is 150 MHz. The actual direction of the first signal is 10° in azimuth and 20° in elevation. The direction of the second signal is 60° in azimuth and -5° in elevation.

Create the signals and array.

fs = 8000;
t = (0:1/fs:1).';
x1 = cos(2*pi*t*300);
x2 = cos(2*pi*t*400);
array = phased.ULA('NumElements',10,'ElementSpacing',1);
array.Element.FrequencyRange = [100e6 300e6];
fc = 150e6;
x = collectPlaneWave(array,[x1 x2],[10 20;60 -5]',fc);
noise = 0.1*(randn(size(x)) + 1i*randn(size(x)));

Solve for the DOAs.

estimator = phased.BeamscanEstimator('SensorArray',array, ...
[y,doas] = estimator(x + noise);
doas = broadside2az(sort(doas),[20 -5])
doas = 

    9.5829   60.3813

Plot the beamscan spectrum


Was this topic helpful?