Documentation

ESPRIT DOA

ESPRIT direction of arrival (DOA) estimator for ULA

Library

Direction of Arrival (DOA)

phaseddoalib

Description

TheESPRIT DOAblock estimates the direction of arrival of a specified number of narrowband signals incident on a uniform linear array using the ESPRIT algorithm.

Parameters

Signal Propagation speed (m/s)

Specify the propagation speed of the signal, in meters per second, as a positive scalar. You can use the functionphysconstto specify the speed of light.

Operating frequency (Hz)

Specify the operating frequency of the system, in hertz, as a positive scalar.

Number of signals

Specify the number of signals as a positive integer scalar.

Spatial smoothing

Specify the amount of averaging,L, used by spatial smoothing to estimate the covariance matrix as a nonnegative integer. Each increase in smoothing handles one extra coherent source, but reduces the effective number of elements by one. The maximum value of this parameter isN – 2, whereNis the number of sensors.

Type of least squares method

Specify the least squares method used for ESPRIT as one ofTLSorLSwhereTLSrefers to total least squares andLSrefers to least squares.

Forward-backward averaging

Select this parameter to use forward-backward averaging to estimate the covariance matrix for sensor arrays with a conjugate symmetric array manifold.

Row weighting factor

Specify the row weighting factor for signal subspace eigenvectors as a positive integer scalar. This parameter controls the weights applied to the selection matrices. In most cases higher value are better. However, the value can never be greater than(N-1)/2whereNis the number of elements of the array.

Simulate using

Block simulation method, specified asInterpreted ExecutionorCode Generation。如果你想让你的块使用MATLAB®interpreter, chooseInterpreted Execution。If you want your block to run as compiled code, chooseCode Generation。Compiled code requires time to compile but usually runs faster.

Interpreted execution is useful when you are developing and tuning a model. The block runs the underlying System object™ in MATLAB. You can change and execute your model quickly. When you are satisfied with your results, you can then run the block usingCode Generation。长在我模拟跑得比他们快nterpreted execution. You can run repeated executions without recompiling. However, if you change any block parameters, then the block automatically recompiles before execution.

When setting this parameter, you must take into account the overall model simulation mode. The table shows how theSimulate usingparameter interacts with the overall simulation mode.

When the Simulink®model is inAcceleratormode, the block mode specified usingSimulate usingoverrides the simulation mode.

Acceleration Modes

Block Simulation Simulation Behavior
Normal Accelerator Rapid Accelerator
Interpreted Execution The block executes using the MATLAB interpreter. The block executes using the MATLAB interpreter. Creates a standalone executable from the model.
Code Generation The block is compiled. All blocks in the model are compiled.

For more information, seeChoosing a Simulation Mode(Simulink).

Array Parameters

Specify sensor array as

Specify a ULA sensor array directly or by using a MATLAB expression.

Types

Array (no subarrays)
MATLAB expression

Number of elements

Specifies the number of elements in the array as an integer.

Element spacing

Specify the spacing, in meters, between two adjacent elements in the array.

Array axis

This parameter appears when theGeometryparameter is set toULAor when the block only supports a ULA array geometry. Specify the array axis asx,y, orz。All ULA array elements are uniformly spaced along this axis in the local array coordinate system.

Taper

Tapers, also known aselement weights, are applied to sensor elements in the array. Tapers are used to modify both the amplitude and phase of the transmitted or received data.

Specify element tapering as a complex-valued scalar or a complex-valued 1-by-Nrow vector. In this vector,Nrepresents the number of elements in the array. IfTaperis a scalar, the same weight is applied to each element. IfTaperis a vector, a weight from the vector is applied to the corresponding sensor element. A weight must be applied to each element in the sensor array.

Expression

A valid MATLAB expression containing a constructor for a uniform linear array, for example,phased.ULA

Sensor Array Tab: Element Parameters

Element type

Specify antenna or microphone type as

  • Isotropic Antenna

  • Cosine Antenna

  • Custom Antenna

  • Omni Microphone

  • Custom Microphone

Exponent of cosine pattern

This parameter appears when you setElement typetoCosine Antenna

Specify the exponent of the cosine pattern as a scalar or a 1-by-2 vector. You must specify all values as non-negative real numbers. When you setExponent of cosine patternto a scalar, both the azimuth direction cosine pattern and the elevation direction cosine pattern are raised to the specified value. When you setExponent of cosine patternto a 1-by-2 vector, the first element is the exponent for the azimuth direction cosine pattern and the second element is the exponent for the elevation direction cosine pattern.

Operating frequency range (Hz)

This parameter appears whenElement typeis set toIsotropic Antenna,Cosine Antenna, orOmni Microphone

Specify the operating frequency range, in hertz, of the antenna element as a 1-by-2 row vector in the form[LowerBound,UpperBound]。The antenna element has no response outside the specified frequency range.

Operating frequency vector (Hz)

This parameter appears whenElement typeis set toCustom AntennaorCustom Microphone

Specify the frequencies, in Hz, at which to set the antenna and microphone frequency responses as a 1-by-Lrow vector of increasing values. UseFrequency responsesto set the frequency responses. The antenna or microphone element has no response outside the frequency range specified by the minimum and maximum elements ofOperating frequency vector (Hz)

Frequency responses (dB)

This parameter appears whenElement typeis set toCustom AntennaorCustom Microphone

Specify this parameter as the frequency response of an antenna or microphone, in decibels, for the frequencies defined byOperating frequency vector (Hz)。SpecifyFrequency responses (dB)as a 1-by-Lvector matching the dimensions of the vector specified inOperating frequency vector (Hz)

Azimuth angles (deg)

This parameter appears whenElement typeis set toCustom Antenna

Specify the azimuth angles at which to calculate the antenna radiation pattern as a 1-by-Prow vector.Pmust be greater than 2. Angle units are in degrees. Azimuth angles must lie between –180° and 180° and be in strictly increasing order.

Elevation angles (deg)

This parameter appears when theElement typeis set toCustom Antenna

Specify the elevation angles at which to compute the radiation pattern as a 1-by-Qvector.Qmust be greater than 2. Angle units are in degrees. Elevation angles must lie between –90° and 90° and be in strictly increasing order.

Radiation pattern (dB)

This parameter appears when theElement typeis set toCustom Antenna

The magnitude in db of the combined polarized antenna radiation pattern specified as aQ-by-Pmatrix or aQ-by-P-by-Larray. The value ofQmust match the value ofQspecified byElevation angles (deg)。The value ofPmust match the value ofPspecified byAzimuth angles (deg_。The value ofLmust match the value ofLspecified byOperating frequency vector (Hz)

Polar pattern frequencies (Hz)

This parameter appears when theElement typeis set toCustom Microphone

Specify the measuring frequencies of the polar patterns as a 1-by-Mvector. The measuring frequencies lie within the frequency range specified byOperating frequency vector (Hz)。Frequency units are in Hz.

Polar pattern angles (deg)

This parameter appears whenElement typeis set toCustom Microphone

Specify the measuring angles of the polar patterns, as a 1-by-Nvector. The angles are measured from the central pickup axis of the microphone, and must be between –180° and 180°, inclusive.

Polar pattern (dB)

This parameter appears whenElement typeis set toCustom Microphone

Specify the magnitude of the microphone element polar pattern as anM-by-Nmatrix.Mis the number of measuring frequencies specified inPolar pattern frequencies (Hz)Nis the number of measuring angles specified inPolar pattern angles (deg)。Each row of the matrix represents the magnitude of the polar pattern measured at the corresponding frequency specified inPolar pattern frequencies (Hz)and all angles specified inPolar pattern angles (deg)。Assume that the pattern is measured in the azimuth plane. In the azimuth plane, the elevation angle is 0° and the central pickup axis is 0° degrees azimuth and 0° degrees elevation. Assume that the polar pattern is symmetric around the central axis. You can construct the microphone’s response pattern in 3-D space from the polar pattern.

Baffle the back of the element

This check box appears only when theElement typeparameter is set toIsotropic AntennaorOmni Microphone

Select this check box to baffle the back of the antenna element. In this case, the antenna responses to all azimuth angles beyond ±90° frombroadsideare set to zero. Define the broadside direction as 0° azimuth angle and 0° elevation angle.

Ports

Note

The block input and output ports correspond to the input and output parameters described in thestepmethod of the underlying System object. See link at the bottom of this page.

Port Description Supported Data Types
In

Input signals.

The size of the first dimension of the input matrix can vary to simulate a changing signal length. A size change can occur, for example, in the case of a pulse waveform with variable pulse repetition frequency.

Double-precision floating point
Ang

Estimated broadside DOA angles.

Double-precision floating point

Introduced in R2014b