Quantcast

Documentation Center

  • Trial Software
  • Product Updates

step

System object: phased.ConstantGammaClutter
Package: phased

Simulate clutter using constant gamma model

Syntax

Y = step(H)
Y = step(H,X)
Y = step(H,STEERANGLE)
Y = step(H,X,STEERANGLE)

Description

Y = step(H) computes the collected clutter return at each sensor. This syntax is available when you set the TransmitSignalInputPort property to false.

Y = step(H,X) specifies the transmit signal in X. Transmit signal refers to the output of the transmitter while it is on during a given pulse. This syntax is available when you set the TransmitSignalInputPort property to true.

Y = step(H,STEERANGLE) uses STEERANGLE as the subarray steering angle. This syntax is available when you configure H so that H.Sensor is an array that contains subarrays and H.Sensor.SubarraySteering is either 'Phase' or 'Time'.

Y = step(H,X,STEERANGLE) combines all input arguments. This syntax is available when you configure H so that H.TransmitSignalInputPort is true, H.Sensor is an array that contains subarrays, and H.Sensor.SubarraySteering is either 'Phase' or 'Time'.

Input Arguments

H

Constant gamma clutter object.

X

Transmit signal, specified as a column vector.

STEERANGLE

Subarray steering angle in degrees. STEERANGLE can be a length-2 column vector or a scalar.

If STEERANGLE is a length-2 vector, it has the form [azimuth; elevation]. The azimuth angle must be between –180 and 180 degrees, and the elevation angle must be between –90 and 90 degrees.

If STEERANGLE is a scalar, it represents the azimuth angle. In this case, the elevation angle is assumed to be 0.

Output Arguments

Y

Collected clutter return at each sensor. Y has dimensions N-by-M matrix. M is the number of subarrays in the radar system if H.Sensor contains subarrays, or the number of sensors, otherwise. When you set the OutputFormat property to 'Samples', N is specified in the NumSamples property. When you set the OutputFormat property to 'Pulses', N is the total number of samples in the next L pulses. In this case, L is specified in the NumPulses property.

Tips

The clutter simulation that ConstantGammaClutter provides is based on these assumptions:

  • The radar system is monostatic.

  • The propagation is in free space.

  • The terrain is homogeneous.

  • The clutter patch is stationary during the coherence time. Coherence time indicates how frequently the software changes the set of random numbers in the clutter simulation.

  • The signal is narrowband. Thus, the spatial response can be approximated by a phase shift. Similarly, the Doppler shift can be approximated by a phase shift.

  • The radar system maintains a constant height during simulation.

  • The radar system maintains a constant speed during simulation.

Examples

expand all

Clutter Simulation of System with Known Power

Simulate the clutter return from terrain with a gamma value of 0 dB. The effective transmitted power of the radar system is 5 kw.

Set up the characteristics of the radar system. This system has a 4-element uniform linear array (ULA). The sample rate is 1 MHz, and the PRF is 10 kHz. The propagation speed is 300,000 km/s, and the operating frequency is 300 MHz. The radar platform is flying 1 km above the ground with a path parallel to the ground along the array axis. The platform speed is 2000 m/s. The mainlobe has a depression angle of 30 degrees.

Nele = 4;
c = 3e8; fc = 3e8; lambda = c/fc;
ha = phased.ULA('NumElements',Nele,'ElementSpacing',lambda/2);

fs = 1e6; prf = 10e3;
height = 1000; direction = [90; 0];
speed = 2000; depang = 30;

Create the clutter simulation object. The configuration assumes the earth is flat. The maximum clutter range of interest is 5 km, and the maximum azimuth coverage is +/– 60 degrees.

Rmax = 5000; Azcov = 120;
tergamma = 0; tpower = 5000;
hclutter = phased.ConstantGammaClutter('Sensor',ha,...
    'PropagationSpeed',c,'OperatingFrequency',fc,'PRF',prf,...
    'SampleRate',fs,'Gamma',tergamma,'EarthModel','Flat',...
    'TransmitERP',tpower,'PlatformHeight',height,...
    'PlatformSpeed',speed,'PlatformDirection',direction,...
    'BroadsideDepressionAngle',depang,'MaximumRange',Rmax,...
    'AzimuthCoverage',Azcov,'SeedSource','Property',...
    'Seed',40547);

Simulate the clutter return for 10 pulses.

Nsamp = fs/prf; Npulse = 10;
csig = zeros(Nsamp,Nele,Npulse);
for m = 1:Npulse
    csig(:,:,m) = step(hclutter);
end

Plot the angle-Doppler response of the clutter at the 20th range bin.

hresp = phased.AngleDopplerResponse('SensorArray',ha,...
    'OperatingFrequency',fc,'PropagationSpeed',c,'PRF',prf);
plotResponse(hresp,shiftdim(csig(20,:,:)),...
    'NormalizeDoppler',true);

Clutter Simulation Using Known Transmit Signal

Simulate the clutter return from terrain with a gamma value of 0 dB. The step syntax includes the transmit signal of the radar system as an input argument. In this case, you do not record the effective transmitted power of the signal in a property.

Set up the characteristics of the radar system. This system has a 4-element uniform linear array (ULA). The sample rate is 1 MHz, and the PRF is 10 kHz. The propagation speed is 300,000 km/s, and the operating frequency is 300 MHz. The radar platform is flying 1 km above the ground with a path parallel to the ground along the array axis. The platform speed is 2000 m/s. The mainlobe has a depression angle of 30 degrees.

Nele = 4;
c = 3e8; fc = 3e8; lambda = c/fc;
ha = phased.ULA('NumElements',Nele,'ElementSpacing',lambda/2);

fs = 1e6; prf = 10e3;
height = 1000; direction = [90; 0];
speed = 2000; depang = 30;

Create the clutter simulation object and configure it to take a transmit signal as an input argument to step. The configuration assumes the earth is flat. The maximum clutter range of interest is 5 km, and the maximum azimuth coverage is +/– 60 degrees.

Rmax = 5000; Azcov = 120;
tergamma = 0;
hclutter = phased.ConstantGammaClutter('Sensor',ha,...
    'PropagationSpeed',c,'OperatingFrequency',fc,'PRF',prf,...
    'SampleRate',fs,'Gamma',tergamma,'EarthModel','Flat',...
    'TransmitSignalInputPort',true,'PlatformHeight',height,...
    'PlatformSpeed',speed,'PlatformDirection',direction,...
    'BroadsideDepressionAngle',depang,'MaximumRange',Rmax,...
    'AzimuthCoverage',Azcov,'SeedSource','Property',...
    'Seed',40547);

Simulate the clutter return for 10 pulses. At each step, pass the transmit signal as an input argument. The software automatically computes the effective transmitted power of the signal. The transmit signal is a rectangular waveform with a pulse width of 2 µs.

tpower = 5000;
pw = 2e-6;
X = tpower*ones(floor(pw*fs),1);
Nsamp = fs/prf; Npulse = 10;
csig = zeros(Nsamp,Nele,Npulse);
for m = 1:Npulse
    csig(:,:,m) = step(hclutter,X);
end

Plot the angle-Doppler response of the clutter at the 20th range bin.

hresp = phased.AngleDopplerResponse('SensorArray',ha,...
    'OperatingFrequency',fc,'PropagationSpeed',c,'PRF',prf);
plotResponse(hresp,shiftdim(csig(20,:,:)),...
    'NormalizeDoppler',true);

More About

Was this topic helpful?