1. Introduction
Corner reflectors, a type of radar-reflecting equipment [
1,
2], redirect incident rays back to their original direction through a three-sided metal corner structure, thereby generating a robust echo signal. Traditional maritime corner reflector devices are typically deployed using carrier fixation and inflatable floating methods. However, in recent years, airborne corner reflectors (ACRs) have emerged as a novel focus in the field of maritime corner reflection. ACRs are launched into the air using conventional carriers for efficient signal reflection. Notable versions of ACRs include [
3] “WIZARD” developed by Rafael Corporation, “OCR” developed by Rheinmetall Corporation, “SEALEM” developed by Lacroix Defense Corporation, and “TORERO” jointly developed by IHI Aerospace Corporation and Chemring Group. These ACRs share similar design characteristics: supplementary parachutes serve as the primary solution for extending their operational duration; Octahedral Reflector (OR) structures constitute the mainstream approach for enhancing reflective capability. This paper primarily focuses on an ACR configuration featuring the typical modality of “parachute + OR”, as illustrated in
Figure 1.
Due to the inherent disparity between ACRs and traditional corner reflectors in terms of modality, scholars have conducted a series of theoretical investigations into their jamming mechanism. Kubicke [
4,
5,
6] proposed a hybrid algorithm based on Physical Optics (PO), Geometric Optics (GO), and the equivalent current method, which was employed to examine the polarization characteristics of OR structures. Zhang [
7] introduced a high-frequency hybrid algorithm grounded in PO and area projection, enabling the establishment of a statistical model of a Radar Cross Section (RCS) [
8] of the OR array. Luo [
9] simulated the RCS of ORs with diverse geometric features using the Shooting and Bouncing Ray (SBR) method, summarizing regulars of its RCS concerning the angles and radius. However, these aforementioned studies primarily concentrate on electromagnetic scattering properties of OR structures under instantaneous static conditions.
One particular application of maritime ACR equipment at sea involves the determination of long-distance sea-surface wind speed. In situations where obtaining accurate ambient wind speed in an unknown sea area is necessary, but telemetry is insufficient to meet accuracy requirements or when it is impractical to directly dispatch a carrier with an anemometer for data collection in close proximity, ACRs offer a unique advantage. By launching an ACR into the desired location at sea, it can steadily move along with the wind [
10]. Then, the velocity of the floating ACR can be measured in real time using radar, based on the Doppler principle. Therefore, studying the Doppler characteristics of ACR targets holds significant practical importance.
The electromagnetic scattering characteristics of airborne motion radar targets, specifically ACRs, have only been analyzed under static conditions, which presents certain limitations and one-sidedness. The existing static target characterization of ACRs does not encompass or support the study of Doppler characteristics. To analyze the dynamic characteristics of ACR targets, a real-time Radar Cross Section (RCS) estimation method for OR structures [
11] and a motion model [
10] of maritime ACRs at sea were proposed in our previous work. These research works laid an important foundation for conducting dynamic Doppler characterization of ACRs. The innovation in this study primarily lies in the following aspects:
Constructing a radar echo signal model of the ACR target based on a pulsed linear frequency modulation (FM) signal;
Analyzing the dynamic Doppler characteristics of ACR targets using the pulse Doppler principle, and finding a special Doppler broadening phenomenon that hinders velocity measurement;
Analyzing the effects of different ACR conditioning parameters on the severity and frequency of Doppler broadening, and finding that the rotational rate of an ACR is the main factor affecting the broadening phenomenon;
A spectrum analysis method based on rectangular window decomposition was proposed to reveal the inherent causes of the Doppler broadening phenomenon and prove that Doppler broadening is indeed closely related to the rotational motion of ACRs.
Section 1 focuses on constructing the radar echo signal model of an ACR.
Section 2 involves obtaining a time series Doppler characteristic image based on pulse Doppler processing and discussing the correlation between multiple condition parameters and Doppler broadening.
Section 3 introduces an analytical method based on rectangular window decomposition and investigates the causes of Doppler broadening.
2. Radar Echo Signal Modeling for ACR
Modern radar technology pursues a greater signal bandwidth. Linear frequency modulation (LFM) signals [
12], commonly used for matched-filter pulse compression, achieve a large time–broadband product by employing nonlinear phase modulation.
The LFM rectangular pulse signal is one of the most commonly used radar-transmitting signals. Its complex expression is
where
Tp is the pulse width,
At is the transmit signal amplitude,
fc is the signal carrier frequency,
μ =
B/
τ is the FM slope,
B is the signal bandwidth, and rect( ) is the rectangular function:
When processing the echo signal from a target, we need to account for various factors like propagation loss, the target reflection coefficient, and the distance traveled. These factors influence both the amplitude and phase of the transmitted signal. The complex expression for the echo signal is
where
Ar is the echo signal amplitude and
td is the echo time delay. The time delay (
td) of the echo in relation to the transmitted signal is directly proportional to the distance (
d) between the radar and the target:
where
c denotes the speed of light.
The center of mass of an ACR target serves as the scattering center [
10]. The distance traveled by all the rays reflected by it equals twice the distance between the radar and the scattering center [
13]. Consequently, the ACR can be regarded as a point-target scatterer. For simple point-target scatterers, the effect of the scattering process on the signal amplitude measured by the radar can be analyzed using the radar range equation.
The radar range equation provides a deterministic model that correlates the received echo power (Pr) to the transmitted signal power (Pt) through various parameters. In the case of high-carrier-frequency radars (X-band), assuming that the adopted LFM bandwidth is much smaller than the signal carrier frequency (B/fc < 1/10), the echo signal can still be considered a narrowband pulse. The received power (Pr) estimated from the range equation can be directly related to the amplitude (Ar) of the received pulse signal.
According to the radar range equation, the total backscattered power obtained by the receiving antenna is [
14,
15]
where
λ =
c/
fc is the signal wavelength,
σ is the RCS value that represents the ability of the target to reflect the radar wave,
G is the antenna power gain,
Ls is the system loss factor, and
Lα(
d) is the atmospheric attenuation factor associated with distance,
d.
Since the radar signal power is proportional to the square of the electric field strength, the amplitude of the target echo signal can be estimated using the following range equation:
Substituting Equations (4) and (6) into Equation (3) yields an echo signal model of an ACR within a pulse repetition interval (PRI).
Therefore, the ACR echo signal is described as a complex model that incorporates multiple parameters. These parameters consist of radar system parameters influenced by factors like the transmit signal amplitude, carrier frequency, bandwidth, and antenna gain. Additionally, there are unknown parameters related to the condition of the ACR target, such as the RCS and distance.
The acquisition of the RCS and distance necessitates additional ACR motion modeling. Reference [
10] provides a comprehensive analysis and introduction of ACR motion characteristics, facilitating the establishment of a simplified ACR motion model. The typical motion process of the ACR consists of two distinct phases: “rapid descent” and “steady descent”. During the rapid descent phase, which lasts only for a few seconds, the ACR descends rapidly with vigorous oscillations and rotations. In this phase, its horizontal speed undergoes changes until it matches that of the wind speed. Subsequently, in the subsequent steady descent phase, the ACR gradually descends at a constant speed while maintaining horizontal rotational motion. At this stage, its horizontal velocity remains consistent with that of the wind speed. Due to the short duration of the rapid descent phase and the complexity of the motion process, a general ACR motion model is established primarily based on the steady descent phase:
where (
x,
y,
z) denotes the real-time position of the ACR; (
vx,
vy) denotes the wind speed in the horizontal direction;
d0 denotes the initial distance between the radar and the ACR;
ϕ0 denotes the initial incidence azimuth;
z0 denotes the initial altitude of the ACR;
vz denotes the vertical descent speed of the ACR, which is a fixed value related to the performance of the parachute and the overall mass of the system;
d denotes the real-time distance between the radar and the ACR;
ϑ denotes the radar incidence pitch;
ϕ denotes the radar incidence azimuth; and
w denotes the angular speed of the ACR’s horizontal rotational motion.
The real-time distance (
d) can be easily obtained from the ACR motion model. After obtaining the longitudinal and transverse angles (
ϑ and
ϕ) of the radar incident wave, the RCS value of the ACR can be further obtained based on the geometrical optics method [
16]:
where
Ae is the equivalent flat plate area, whose formula has different forms at different angles [
11]:
- (1)
When ϑ approaches 0°, or ϑ approaches 90° and ϕ approaches 0° or 90°,
- (2)
When ϑ approaches 90°, ϕ approaches 0°, or ϕ approaches 90°,
- (3)
When 2° < ϑ < 88° and 2° < ϕ < 88°,
L is the length of the inner edge of the ACR, typically representing its size; ξ is the standard deviation of the normal distribution used to model the changing trend of Ae when the aspect angle is approaching the angular split.
3. Analysis of Dynamic Doppler Characteristics
The processing of echo signals by radar typically involves two steps, fast-time processing and slow-time processing [
17], as illustrated in
Figure 2. Fast-time processing, often referred to as pulse compression techniques, focuses on processing a single pulse echo to extract a range of information about the target. However, radars typically emit not just a single pulse but a periodic sequence of pulses. Typically, a coherent processing interval (CPI) of
M pulses is utilized. Within each CPI, slow-time processing can gather velocity information about the target.
The echo signal of multiple pulses can be expressed as follows [
18]:
The above equation contains two independent variables, where
t is the fast-time independent variable and
tm is the slow-time independent variable.
tm denotes the time delay of the
m-th pulse echo, which is expressed as
where
Tr denotes the PRI; the radial velocity (
vd) denotes the projection of the target velocity vector in the radar line-of-sight direction, which is positive when the target is moving in the direction of approaching the radar, and negative when it is moving in the opposite direction.
Pulse compression on each of the
M signals yields
M fast-time samples:
At this point, the slow-time independent variable is transitioned from the time delay (
tm) to
m, denoting the PRI sequence number. The fast-time independent variable (
t) is regarded as a constant during slow-time processing. By further isolating the terms independent of
m to simplify the equation, the slow-time series with respect to
m is obtained:
where Δ
d =
c/2
B denotes the distance resolution.
The spectral form of the slow-time series
χ[
m] is complex but can be analyzed using the Discrete Time Fourier Transform (DTFT) [
19]. Let
be a DTFT pair.
ω denotes the normalized digital domain corner frequency. Despite the complex form of
X(
ω), there is no shifting of the spectrum. The DTFT of the slow-time series
χ[
m] can be obtained from the frequency shift property as follows:
The spectral center of the slow-time series
χ[
m] is shifted to
ω0 = 2
πfc(
2vd/
c)
Tr. In addition, the velocity information of the target is included in the position of the spectral center of the slow-time series. Hence, the extraction of velocity information is accomplished by performing the slow-time Fast Fourier Transform (FFT) to identify the position of the spectral center. The relationship between the digital domain angular frequency (
ω) and the true frequency (
f) is
The Doppler frequency shift equation can be derived from the following equation:
To further analyze the spectrum of the slow-time series
χ[
m], Equation (15) is simplified by considering the general case Δ
d <<
vdTr:
The spectrum of
χ[
m] is obtained by applying a DTFT to the above equation:
Y(f) is an asinc function with a peak at f = fd and an amplitude peak at MA. The Doppler Equation (19) maps the frequency domain to the velocity domain and yields the radial velocity of the target by finding the location of the peak of the spectrum Y(f).
The signal processing process described above is also known as pulsed Doppler processing [
20,
21,
22,
23]. In actual pulsed Doppler processing, the set of fast-time sample series with length
N obtained from each
M PRIs is typically regarded as a two-dimensional fast-time–slow-time data matrix
χ[
n,
m] within a CPI. Here, the dimension (
n) where the range cells are located is termed the fast-time axis, while the dimension
(m) where the number of pulses is located is termed the slow-time axis. Pulse Doppler processing involves extracting the slow-time data series per range cell from the data matrix
χ[
n,
m] and conducting Discrete Fourier Transform (DFT) [
24,
25,
26,
27] on it. Ultimately, the fast-time–slow-time data matrix is transformed into a Range-Doppler (RD) data matrix
Y[
n,
m], as depicted in
Figure 3.
The following parameters for radar operation are configured: carrier frequency:
fc = 10 GHz, wavelength:
λ = 0.03 m, transmit power:
Pt = 35 kW, antenna gain:
G = 60, bandwidth:
B = 60 MHz, pulse width:
Tp = 2 × 10
−4 s, pulse repetition interval:
Tr = 1 × 10
−3 s, coherent processing period:
M = 64, and initial incident azimuthal angle:
ϕ0 = 0°. In the hypothetical case, an ACR target with a size
L = 1 m appears at the moment
t = 0 at a range
d0 = 1000 m from the radar, with an initial height of 100 m in the air, a horizontal rotational rate of 0.5 rad/s, a dropping velocity of 1 m/s, and an ambient wind speed of (5, 0, 0) m/s. At this point, the fast-time–slow-time data matrix is obtained by pulse compression, i.e., the fast-time pulse compression output of all pulse signals within a CPI, as illustrated in
Figure 4.
In the above figure, the horizontal coordinate represents the fast-time (range) axis, while the vertical coordinate represents the slow-time (number of pulses) axis. DFT was separately performed on each row of the slow-time data series, resulting in the RD data plot as illustrated in
Figure 5.
The slow-time series within the range cell where the target is situated is compressed into a spike. Furthermore, the spike’s location is consistently indicated to be at 5 m/s in the velocity dimension, aligning with the actual radial velocity magnitude. It indicates that the pulsed Doppler processing effectively extracted the velocity information of the target.
To analyze the dynamic Doppler characteristics of the ACR target over time, a single CPI was regarded as a frame, wherein each frame consisted of M PRIs. Pulse compression was subsequently performed on the echo signals within each PRI to obtain fast-time–slow-time data matrices. Following this, spectral analysis is performed on the slow-time series of the range cell where the target is located.
With
T0 = 0 s as the initial moment and
Tend = 2 s as the final moment, a total of 31 frames were generated under the set PRI and CPI conditions. The generated dynamic sequential Doppler spectrogram is shown in
Figure 6.
Over time, the peak of the Doppler spectrum consistently close to 5 m/s. This occurs because the flight direction of the ACR aligns with the radar’s line-of-sight direction, and the flight speed of the ACR is very close to the wind speed (5 m/s). An interesting phenomenon arises when the Doppler spectrum of the ACR target exhibits significant broadening at specific moments, such as frames 4, 12, 20, and 28.
The presence of Doppler broadening is detrimental to detecting target velocity. When broadening occurs, the energy that was originally concentrated in the spectrum is dispersed, which perhaps causes false dismissal of an ACR target because the peak value of the spectrum may be lower than the detection threshold [
28,
29,
30], as shown in
Figure 7. Target false dismissal will destabilize target tracking and perplex later data analysis. Therefore, the Doppler broadening phenomenon should be minimized or avoided as far as possible.
In order to quantify the severity of Doppler broadening, the kurtosis [
31,
32,
33,
34] of the spectrum in the
t-th frame is extracted:
where
nd is the number of range cells where the ACR target is located,
ρ is the standard deviation of the spectral amplitude, and
Yave is the average amplitude of the spectrum.
Figure 8 shows time-sequential kurtosis extracted from the time-sequential spectrogram illustrated in
Figure 6.
Kurtosis (
K) is a measure of fat tails for the probability density distribution of the spectrum. A higher
K value, which represents a sharper Doppler spectrum shape, is more conducive to speed detection. In contrast, a lower
K value means that the Doppler broadening is more severe. In
Figure 8, the kurtosis (
K) varies with time-sequential frames. As one can see, the kurtosis is generally high in whole frames, while there is significant attenuation at a few moments, such as frames 4, 12, 20, and 28. The appearances of these kurtosis valleys exactly corresponds to the moments of Doppler broadening in
Figure 6, which proves the rationality of employing kurtosis to measure the severity of Doppler broadening.
Furthermore, in order to quantitatively analyze the frequency of Doppler broadening, we determine whether the false dismissal occurs at frame
t by
where
Q[
t] = 1 represents the occurrence of target false dismissal in the
t-th frame, max( ) is the operation of taking the maximum value from the data, and
Thr is a fixed detection threshold.
Figure 9 shows time-sequential judgements for the occurrence of false dismissal based on the time-sequential spectrogram illustrated in
Figure 6.
It can be seen that there is a false dismissal of the target at frames 4, 12, 20, and 28, respectively, which exactly corresponds to the moments when the Doppler broadening occurs in
Figure 6. By summing
Q[
m], it can be concluded that the total times of false dismissal is four, which is consistent with the occurrence of Doppler broadening illustrated in
Figure 6, proving the rationality of employing the number of false dismissals to measure the frequency of Doppler broadening.
Further analysis was conducted on the effects of different target conditions on the Doppler broadening phenomenon. Ten different values were selected for each of the five types of ACR target condition parameters: size, rotational rate, dropping velocity, height, and range. All values followed a linear increasing trend, as outlined in
Table 1.
Continuing to adopt the simulation parameter settings mentioned in
Figure 4, a therein single parameter would be transformed in turn by one value selected from
Table 1 for each experiment. A total of 50 experiments were conducted to obtain 50 sets of time-sequential Doppler spectra. Then,
K[
m] and
Q[
m] were obtained from the 50 Doppler spectra. Finally, we collected the minimum kurtosis, min(
K[
m]), and the false dismissal frequency, sum(
Q[
m]), representing the severity and frequency of the Doppler broadening phenomenon separately, as shown in
Figure 10.
As one can see in
Figure 10a, the severity of Doppler broadening is mainly related to the rotational rate of ACR. As the rotational rate gradually increases, the minimum kurtosis decreases, indicating that the degree of Doppler broadening becomes more severe. When the other four parameters change, the minimum kurtosis is almost constant, indicating that these parameters have little impact on the severity of the Doppler broadening phenomenon.
In
Figure 10b, it can be seen that the frequency of Doppler broadening is also mainly related to the rotational rate of the ACR. As the rotational rate gradually increases, the number of false dismissals increases, indicating that the appearance of Doppler broadening becomes more frequent. When the other four parameters change, the number of false dismissals remains stable at around four, indicating that these parameters have little impact on the frequency of the Doppler broadening phenomenon.
Based on the above analysis, it can be concluded that the Doppler broadening phenomenon is likely related to the rotational motion of the ACR. The faster the ACR rotates, the more severe and more frequent the Doppler broadening phenomenon appears. Therefore, in order to avoid missed detections and achieve stable tracking of ACR targets, the design of ACRs should try to suppress their rotational motion trend as far as possible.
4. Analysis of The Doppler Broadening Phenomenon Based on Rectangular Window Decomposition
To analyze the cause of the Doppler broadening phenomenon, the slow-time series and its spectrum at a particular moment of Doppler broadening were extracted as shown in
Figure 11.
In
Figure 11a, the horizontal coordinates represent pulse numbers, while the vertical coordinates depict the slow-time series amplitudes. In
Figure 11b, the horizontal coordinates represent velocities, and the vertical coordinates represent spectral amplitudes. According to
Figure 11, the Doppler spectrum at this moment exhibits significant broadening, while the corresponding original slow-time series demonstrates noticeable changes in amplitude oscillation. Conventional targets such as aircraft or ships typically maintain a constant echo amplitude (backscattered energy stabilization) within a CPI due to their stable navigational attitude. Consequently, the target’s velocity information can be extracted by analyzing the phase change information of the slow-time series through DFT. However, the ACR target remains in a spin state for a prolonged period during airborne flight following launch, owing to its rotational rate. This results in a constant change in echo amplitude with attitude angle. Moreover, the slow-time series amplitude fluctuation of the ACR target is highly complex due to the uneven distribution of backscattered energy across the ACR’s angular domain. The amplitude fluctuation in the slow-time series generated by the ACR’s spin can be viewed as applying a specialized window function to the slow-time series. The form of this window function is intricate, and there is no universal expression for arbitrary conditions. Therefore, it necessitates case-by-case analysis.
As depicted in
Figure 12, the slow-time series is simplified to the sum of four rectangular window functions for ease of analysis, with the following expressions:
By substituting the above equations into Equation (20), the slow-time series is
The following can be obtained according to the linear nature of the Fourier transform:
The Doppler spectrum of the simplified rectangular window slow-time series χ’[m] can be obtained from the above equation, as illustrated in
Figure 13.
The variation in the total rectangular window slow-time series spectrum (
Y) in
Figure 13 closely resembles that of the actual spectrum in
Figure 11b. The broadening phenomenon observed in the ACR Doppler spectrum is precisely caused by the variation in the echo amplitude resulting from its spin. Analysis of the spectra of the slow-time series using four sub-rectangular window functions (
Y1,
Y2,
Y3, and
Y4) indicates that the broadening of the spectra is primarily related to the rectangular window A2. The higher amplitude of the rectangular window (A2) and its narrower width accentuate the flat trend of the asinc function in the spectrum.
To enhance generalization, the arbitrary slow-time series amplitude is expressed as the sum of
M/2
q rectangular window functions with the following expressions:
Then, the slow-time series can be re-expressed as
DTFT on the new slow-time series
χ′[
m] yields
The spectrum function has the following characteristics:
When f = fd, the phase of the spectra generated by all rectangular window slow-time series is 0. In this case, the total spectral amplitude equals the sum of all sub-spectral amplitudes. However, when f ≠ fd, the spectra of all rectangular window slow-time series do not interact with each other. Consequently, the amplitude will decay significantly. Therefore, the spectrum always peaks at f = fd regardless of the value of q;
By changing the value of q, different fitting effects can be achieved. A larger q value results in fewer required rectangular windows and a simpler shape of the spectrum. Conversely, a smaller q value results in more required rectangular windows and a more complex shape of the spectrum;
A small q value results in flatter asinc functions. When the amplitude of a specific rectangular window, Ax, is significantly higher than the other rectangular windows, the interference effect of the other rectangular windows on the spectrum of Ax is reduced. This leads to an apparent broadening phenomenon in the total spectrum. Additionally, if there are multiple rectangular windows with prominent amplitudes and small phase differences between their spectra, the total spectrum will also broaden.
Based on the above analysis, the Doppler broadening of ACRs is caused by the unique variation in amplitude in the slow-time series. When the amplitude of a consecutive set of pulses (Mh) in the slow-time series is significantly higher than the others, it leads to Doppler spectral broadening. This broadening phenomenon becomes more noticeable with smaller values of Mh. The variation in amplitude in the slow-time series is determined by the echo amplitude (Ar), which is, in turn, determined by the RCS value of the ACR.
The RCS variation of an ACR with the aspect angle [
11] is depicted in
Figure 14. The rotational motion of the ACR induces continuous changes in the relative incidence angle of the radar beam, consequently leading to fluctuations in its RCS. As shown in
Figure 14, there will be a sharp increase or decrease in ACR backscattered energy near certain angles. The Doppler broadening phenomenon exhibited by the ACR arises precisely due to the fact that the aspect angle aligns with these critical angles within a coherent processing interval (CPI), resulting in a few pulse repetition intervals (PRIs) collecting higher amplitudes and forming distinctive amplitude variations in a slow-time series.
In conclusion, the rotational motion of an ACR leads to a sharp change in its RCS amplitude at certain moments, which is reflected as broadening in the spectrum. This explains why the severity and frequency of Doppler broadening are related to the rotational rate and proves the correctness of the analysis results in
Section 3.