singermeas

Measurement function for Singer acceleration motion model

Description

example

measurements = singermeas(states) returns the measurements in rectangular coordinates for the Singer motion model based on the current states.

measurements = singermeas(states,frame) specifies the measurement output coordinate system, frame.

measurements = singermeas(states,frame,sensorpos,sensorvel) also specifies the sensor position, sensorpos, and the sensor velocity, sensorvel.

measurements = singermeas(states,frame,sensorpos,sensorvel,laxes) specifies the local sensor axes orientation, laxes.

measurements = singermeas(states,measurementParameters) specifies the measurement parameters, measurementParameters.

Examples

collapse all

Define a state for a 2-D Singer acceleration motion.

state = [1;10;3;2;20;5];

Obtain the measurement in a rectangular frame.

measurement = singermeas(state)
measurement = 3×1

1
2
0

Obtain the measurement in a spherical frame.

measurement = singermeas(state, 'spherical')
measurement = 4×1

63.4349
0
2.2361
22.3607

Obtain the measurement in a spherical frame relative to a stationary sensor located at [1;-2;0].

measurement = singermeas(state, 'spherical', [1;-2;0], [0;0;0])
measurement = 4×1

90
0
4
20

Obtain the measurement in a spherical frame relative to a stationary sensor located at [1;-2;0] that is rotated by 90 degrees around the z axis relative to the global frame.

laxes = [0 -1 0; 1 0 0; 0 0 1];
measurement = singermeas(state, 'spherical', [1;-2;0], [0;0;0], laxes)
measurement = 4×1

0
0
4
20

Obtain measurements from multiple 2D states in a rectangular frame.

states = [1 2 3; 10 20 30; 2 4 5; 20 30 40; 5 6 11; 1 3 1.5];
measurements = singermeas(states)
measurements = 3×3

1     2     3
20    30    40
0     0     0

Input Arguments

collapse all

Current states, specified as a real-valued 3N-by-1 vector or a real-valued 3N-by-M matrix. N is the spatial degree of the state, and M is the number of states.

The state vector in each column takes different forms based on its spatial dimensions.

Spatial DegreesState Vector Structure
1-D[x;vx;ax]
2-D[x;vx;ax;y;vy;ay]
3-D[x;vx;ax;y;vy;ay;z;vz;az]

For example, x represents the x-coordinate, vx represents the velocity in the x-direction, and ax represents the acceleration in the x-direction. If the motion model is in one-dimensional space, the y- and z-axes are assumed to be zero. If the motion model is in two-dimensional space, values along the z-axis are assumed to be zero. Position coordinates are in meters. Velocity coordinates are in meters/second. Acceleration coordinates are in m/s2.

Example: [5;0.1;0.01;0;-0.2;-0.01;-3;0.05;0]

Measurement output frame, specified as 'rectangular' or 'spherical'. When the frame is 'rectangular', a measurement consists of x, y, and z Cartesian coordinates. When specified as 'spherical', a measurement consists of azimuth, elevation, range, and range rate.

Data Types: char

Sensor position with respect to the navigation frame, specified as a real-valued 3-by-1 column vector. Units are in meters.

Data Types: double

Sensor velocity with respect to the navigation frame, specified as a real-valued 3-by-1 column vector. Units are in m/s.

Data Types: double

Local sensor coordinate axes, specified as a 3-by-3 orthogonal matrix. Each column specifies the direction of the local x-, y-, and z-axes, respectively, with respect to the navigation frame. That is, the matrix is the rotation matrix from the global frame to the sensor frame.

Data Types: double

Measurement parameters, specified as a structure or an array of structures. For more details, see Measurement Parameters.

Data Types: struct

Output Arguments

collapse all

Measurement vector, returned as an N-by-1 column vector of scalars or an N-by-M matrix of scalars. The form of the measurement depends upon which syntax you use.

• When the syntax does not use the measurementParameters argument, the measurement vector is [x,y,z] when the frame input argument is set to 'rectangular' and [az;el;r;rr] when the frame is set to 'spherical'.

• When the syntax uses the measurementParameters argument, the size of the measurement vector depends on the values of the Frame, HasVelocity, and HasElevation fields in the measurementParameters structure.

FrameMeasurement
'spherical'

Specifies the azimuth angle, az, elevation angle, el, range, r, and range rate, rr of the measurements.

Spherical Measurements

HasElevation
falsetrue
HasVelocityfalse[az;r][az;el;r]
true[az;r;rr][az;el;r;rr]

Angle units are in degrees, range units are in meters, and range rate units are in m/s.

'rectangular'

Specifies the Cartesian position and velocity coordinates of the measurements.

Rectangular Measurements

 HasVelocity false [x;y;z] true [x;y;z;vx;vy;vz]

Position units are in meters and velocity units are in m/s.

Data Types: double

collapse all

Azimuth and Elevation Angle Definitions

Define the azimuth and elevation angles used in the toolbox.

The azimuth angle of a vector is the angle between the x-axis and its orthogonal projection onto the xy plane. The angle is positive in going from the x axis toward the y axis. Azimuth angles lie between –180 and 180 degrees. The elevation angle is the angle between the vector and its orthogonal projection onto the xy-plane. The angle is positive when going toward the positive z-axis from the xy plane. Measurement Parameters

The MeasurementParameters property consists of an array of structures that describe a sequence of coordinate transformations from a child frame to a parent frame or the inverse transformations (see Frame Rotation). If MeasurementParameters only contains one structure, then it represents the rotation from one frame to the other. If MeasurementParameters contains an array of structures, then it represents rotations between multiple frames.

The fields of MeasurementParameters are shown here. Not all fields have to be present in the structure.

 Field Description Frame Enumerated type indicating the frame used to report measurements. When detections are reported using a rectangular coordinate system, Frame is set to 'rectangular'. When detections are reported in spherical coordinates, Frame is set to 'spherical' for the first structure. OriginPosition Position offset of the origin of the child frame relative to the parent frame, represented as a 3-by-1 vector. OriginVelocity Velocity offset of the origin of the child frame relative to the parent frame, represented as a 3-by-1 vector. Orientation 3-by-3 real-valued orthonormal frame rotation matrix. The direction of the rotation depends on the IsParentTochild field. IsParentToChild A logical scalar indicating whether Orientation performs a frame rotation from the parent coordinate frame to the child coordinate frame. If false, Orientation performs a frame rotation from the child coordinate frame to the parent coordinate frame instead. HasElevation A logical scalar indicating if the measurement includes elevation. For measurements reported in a rectangular frame, if HasElevation is false, the measurements are reported assuming 0 degrees of elevation. HasAzimuth A logical scalar indicating if the measurement includes azimuth. HasRange A logical scalar indicating if the measurement includes range. HasVelocity A logical scalar indicating if the reported detections include velocity measurements. For measurements reported in the rectangular frame, if HasVelocity is false, the measurements are reported as [x y z]. If HasVelocity is true, measurements are reported as [x y z vx vy vz].

 Singer, Robert A. "Estimating optimal tracking filter performance for manned maneuvering targets." IEEE Transactions on Aerospace and Electronic Systems 4 (1970): 473-483.

 Blackman, Samuel S., and Robert Popoli. "Design and analysis of modern tracking systems." (1999).

 Li, X. Rong, and Vesselin P. Jilkov. "Survey of maneuvering target tracking: dynamic models." Signal and Data Processing of Small Targets 2000, vol. 4048, pp. 212-235. International Society for Optics and Photonics, 2000.