Next Article in Journal
An Improved Artificial Potential Field Method for Ship Path Planning Based on Artificial Potential Field—Mined Customary Navigation Routes
Next Article in Special Issue
Response of Shallow-Water Temperature and Significant Wave Height to Sequential Tropical Cyclones in the Northeast Beibu Gulf
Previous Article in Journal
Hydrodynamic Interactions and Enhanced Energy Harnessing amongst Many WEC Units in Large-Size Wave Parks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Doppler Characteristics of Maritime Airborne Corner Reflector

College of Weaponry Engineering, Naval University of Engineering, Wuhan 430033, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2024, 12(5), 727; https://doi.org/10.3390/jmse12050727
Submission received: 24 March 2024 / Revised: 21 April 2024 / Accepted: 26 April 2024 / Published: 27 April 2024
(This article belongs to the Special Issue Ocean Observations)

Abstract

:
The maritime airborne corner reflector (ACR) is a radar reflector that can measure wind speed in an unknown sea area in real time over a long distance. To improve our understanding of how the ACR works, we investigated the Doppler characteristics of the ACR for the first time from a dynamic perspective. First, we constructed a radar echo signal model of the ACR. Then, we obtained the dynamic Doppler characteristics through pulse Doppler processing and discussed the special phenomenon of Doppler broadening. Finally, we proposed a rectangular window decomposition method to analyze the inner principle of the Doppler broadening phenomenon in more detail. In conclusion, this study provides valuable insights into the Doppler characterization of an ACR from a dynamic viewpoint, which contributes to enriching the basic theory of this equipment.

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
s t ( t ) = rect t T p A t exp j 2 π f c t + μ 2 t 2
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:
rect t T p = 1 t / T p 1 / 2 0 t / T p > 1 / 2
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
s r ( t ) = rect t t d T p A r exp j 2 π f c ( t t d ) + μ 2 ( t t d ) 2
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:
t d = 2 d c
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]
P r = P t G 2 λ 2 σ 4 π 3 d 4 L s L α d
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:
A r = A t 2 G 2 λ 2 σ 4 π 3 d 4 L s L α d   = A t G λ 8 d 2 σ π 3 L s L α d
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:
x = d 0 2 z 0 2 cos ϕ 0 + v x t y = d 0 2 z 0 2 sin ϕ 0 + v y t z = z 0 + v z t d = x 2 + y 2 + z 2 ϑ = π 2 + arcsin z d ϕ = ϕ 0 + w t
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]:
σ = 4 π A e 2 / λ 2
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°,
A e = 2 exp ϑ 2 / ξ 2 · L 2
(2)
When ϑ approaches 90°, ϕ approaches 0°, or ϕ approaches 90°,
A e = 2 exp ϑ 90 2 / ξ 2 · L 2 · sin ϑ 0 < ϕ 45 ° cos ϑ 45 ° < ϕ < 90 °   if   ϑ   c l o s e   t o   90 ° A e = 2 exp ϕ 2 / ξ 2 · L 2 · sin ϕ 0 < ϑ 45 ° cos ϕ 45 ° < ϑ < 90 °         if   ϕ   c l o s e   t o   0 °   A e = 2 exp ϕ 90 2 / ξ 2 · L 2 · sin ϕ 0 < ϑ 45 ° cos ϕ 45 ° < ϑ < 90 °   if   ϕ   c l o s e   t o   90 °
(3)
When 2° < ϑ < 88° and 2° < ϕ < 88°,
A e = L 2 · 2 sin 2 ϑ · sin 2 ϕ + 90 ° / f ϑ , ϕ ϑ < ϑ 1 f ϑ , ϕ 2 / f ϑ , ϕ     ϑ 1 < ϑ < ϑ 2 2 sin ϕ + 45 ° · sin 2 ϑ / f ϑ , ϕ ϑ > ϑ 2 ϑ 1 = arccot 2 sin ϕ + 45 ° ϑ 2 = arccot 2 cos ϕ + 45 ° f ϑ , ϕ = sin ϑ · cos ϕ + sin ϕ + cos ϑ
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]:
s r ( t , t m ) = rect t t m T p A r exp j 2 π f c ( t t m ) + μ 2 ( t t m ) 2
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
t m = 2 d v d m T r c 1 m M
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:
χ ( t , m ) = A r T p sin c B t 2 d v d m T r c exp j 2 π f c 2 d v d m T r c
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:
χ m = A r T p exp j 2 π f c 2 d c sin c B t 2 d c + m v d T r Δ d exp j 2 π f c 2 v d c m T r
where Δd = c/2B 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
A r T p exp j 2 π f c 2 d c sin c B t 2 d c + m v d T r Δ d X ω
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:
χ m Y ω = X ω 2 π f c 2 v d c T r
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
ω 0 = 2 π f c 2 v d c T r = 2 π f d f r
The Doppler frequency shift equation can be derived from the following equation:
f d = f c 2 v d c
To further analyze the spectrum of the slow-time series χ[m], Equation (15) is simplified by considering the general case Δd << vdTr:
χ m = A r T p exp j 2 π f c 2 d c sin c B t d Δ d exp j 2 π f d m T r           = A exp j 2 π f d m T r           m = 0 , , M 1
The spectrum of χ[m] is obtained by applying a DTFT to the above equation:
Y f = m = + χ m exp j 2 π f m T r         = A sin π f f d M T r sin π f f d T r exp j π M 1 f f d T r         P R F 2 f < P R F 2
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:
K [ t ] = 1 M m = 1 M Y [ n d , m ] Y a v e 4 ρ 4
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
Q [ t ] = 0 w h e n   max Y [ n d , m ] > T h r 1 w h e n   max Y [ n d , m ] < T h r
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:
A m = A 1 m + A 2 m + A 3 m + A 4 m m = 0 , , M 1 A 1 m = a 1 m = 0 , , M 1 1 0   e l s e A 2 m = a 2 m = M 1 , , M 2 1 0 e l s e A 3 m = a 3 m = M 2 , , M 3 1 0 e l s e A 4 m = a 4 m = M 3 , , M 1 0 e l s e
By substituting the above equations into Equation (20), the slow-time series is
χ m = χ 1 m + χ 2 m + χ 3 m + χ 4 m   = A 1 m exp j 2 π f d m T r + A 2 m exp j 2 π f d m T r     + A 3 m exp j 2 π f d m T r + A 4 m exp j 2 π f d m T r     m = 0 , , M 1
The following can be obtained according to the linear nature of the Fourier transform:
Y f = DTFT χ 1 m + χ 2 m + χ 3 m + χ 4 m = Y 1 f + Y 2 f + Y 3 f + Y 4 f = a 1 sin π f f d M 1 T r sin π f f d T r exp j π M 1 1 f f d T r + a 2 sin π f f d M 2 M 1 T r sin π f f d T r exp j π M 1 + M 2 1 f f d T r + a 3 sin π f f d M 3 M 2 T r sin π f f d T r exp j π M 2 + M 3 1 f f d T r + a 4 sin π f f d M M 3 T r sin π f f d T r exp j π M + M 3 1 f f d T r P R F 2 f < P R F 2
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/2q rectangular window functions with the following expressions:
A m = i = 0 M 2 q 1 A i m m = 0 , , M 1 A i m = χ 2 q i m = 2 q i , , 2 q i + 2 q 0 e l s e
Then, the slow-time series can be re-expressed as
χ m = i = 0 M 2 q 1 A i m exp j 2 π f d m T r m = 0 , , M 1
DTFT on the new slow-time series χ′[m] yields
Y f = DTFT i = 0 M 2 q 1 A i m exp j 2 π f d m T r = i = 0 M 2 q 1 χ 2 q i sin π f f d 2 q T r sin π f f d T r exp j π 4 q i + 2 q 1 f f d T r = sin π f f d 2 q T r sin π f f d T r i = 0 M 2 q 1 χ 2 q i exp j π 4 q i + 2 q 1 f f d T r P R F 2 f < P R F 2
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 ffd, 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.

5. Conclusions

Based on the constructed ACR radar echo signal model, the dynamic Doppler characteristics are investigated through pulse Doppler processing. The effects of various ACR condition parameters on the severity and frequency of Doppler broadening are analyzed. Finally, a simplified rectangular window method is proposed for in-depth analysis, which helps to identify the intrinsic causes of the Doppler broadening phenomenon. The detailed analysis of the dynamic Doppler characteristics of ACRs enhances the understanding of their working mechanisms, laying a theoretical foundation for future in-depth research on the application, detection, and identification of this type of equipment.
This research was conducted during the steady dropping phase of an ACR. Analyzing the Doppler characteristics of ACRs in the initial stage after launch poses greater challenges, which will be the focus of our future endeavors.

Author Contributions

Conceptualization, L.W.; methodology, L.W. and S.H.; software, L.W. and S.H.; validation, Z.L.; formal analysis, L.W. and Y.L.; investigation, L.L. and Y.L.; resources, L.W. and C.F.; data curation, L.W.; writing—original draft preparation, L.W.; writing—review and editing, L.W.; visualization, L.W.; supervision, C.F.; project administration, S.H.; funding acquisition, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cardillo, E.; Ferro, L. Multi-Frequency Analysis of Microwave and Millimeter-Wave Radars for Ship Collision Avoidance. In Proceedings of the Microwave Mediterranean Symposium (MMS), Pizzo Calabro, Italy, 9–13 May 2022. [Google Scholar] [CrossRef]
  2. Hu, H.; Zhang, L.; Zhang, X. A study on equipment developments and operational using of ship born gas-filled anti-missile multi-corner reflector. Def. Technol. Rev. 2018, 39, 74–77. [Google Scholar]
  3. Gan, L.; Sun, G.; Feng, D.; Li, J. Characteristics of an Eight-Quadrant Corner Reflector Involving a Reconfigurable Active Metasurface. Sensors 2022, 22, 4715. [Google Scholar] [CrossRef]
  4. Kubicke, G.; Bourlier, C.; Saillard, J. Comparison of bistatic signatures of octahedral and icosahedral reflectors in high-frequency domain. In Proceedings of the European Microwave Conference, Munich, Germany, 9–12 October 2007. [Google Scholar]
  5. Kubicke, G.; Bourlier, C.; Saillard, J. Polarimetric bistatic signature of a faceted octahedron in high-frequency domain. Prog. Electromagn. Res. 2007, 71, 173–209. [Google Scholar] [CrossRef]
  6. Kubicke, G.; Bourlier, C.; Saillard, J. High-Frequency Bistatic Scattering by Depolarizing, Nearly Omnidirectional Reflectors: Higher Order Polyhedral Reflectors. IEEE Trans. Antennas Propag. 2008, 56, 3029–3035. [Google Scholar] [CrossRef]
  7. Zhang, J.; Hu, S.; Wang, P.; Xueman, F. RCS model construction and analysis for corner reflector based on PO/AP algorithm. Syst. Eng. Electron. 2018, 40, 1478–1485. [Google Scholar]
  8. Cong, Z.; He, Z.; Taiping, S.; Ding, D.; Chen, R. A local current based iterative physical optics method for computing the RCS of target above sea surface. Eng. Anal. Bound. Elem. 2023, 146, 318–326. [Google Scholar] [CrossRef]
  9. Luo, Y.; Guo, L.; Zuo, Y. Investigations of Effects of Geometric Characteristics on RCS for Corner Reflectors. In Proceedings of the International Conference on Computational Electromagnetics, Shanghai, China, 20–22 March 2019. [Google Scholar]
  10. Wu, L.; Hu, S.; Xu, J.; Liu, Z. Dynamic sequential radar cross section properties of airborne corner reflector in array. IET Radar Sonar Navig. 2023, 17, 1405–1419. [Google Scholar] [CrossRef]
  11. Wu, L.; Hu, S.; Zhang, J.; Weng, Y.; Li, S.; Yang, J.; Yi, H. Fast RCS Estimation Method for Bipyramid Corner Reflector. Tactical Missile Technol. 2021, 5, 29–35. [Google Scholar]
  12. Bai, J.; Qin, P.; Li, H.; Zhang, R.; Zhang, X. An Improved Adaptive Pulse Compression Algorithm Based on Linear Frequency Modulation Signal. In Proceedings of the IEEE International Conference on Signal, Information and Data Processing (ICSIDP), Chongqing, China, 11–13 December 2019. [Google Scholar] [CrossRef]
  13. Wu, L.; Xu, J.; Hu, S.; Liu, Z. High-frequency backscattering properties of quasi-omnidirectional corner reflector: The great-icosahedral-like reflector. AIP Adv. 2022, 12, 105225. [Google Scholar] [CrossRef]
  14. Ding, L.; Geng, F.; Chen, J. Radar Principles, 6th ed.; Publishing House of Electronics Industry: Beijing, China, 2020; pp. 47–65. [Google Scholar]
  15. Richards, M.A. Fundamentals of Radar Signal Processing, 2nd ed.; McGraw-Hill Education: New York, NY, USA, 2005; pp. 34–45. [Google Scholar]
  16. Fan, X.; Hu, S.; He, J. High-frequency method for the evaluation of the radar cross section of corner reflectors. Chin. J. Radio Sci. 2016, 31, 331–335+362. [Google Scholar]
  17. Bjorklund, S.; Nelander, A.; Pettersson, M.I. Fast-Time and Slow-Time Space-Time Adaptive Processing for bistatic radar interference suppression. In Proceedings of the IEEE National Radar Conference, Arlington, VA, USA, 10–15 May 2015. [Google Scholar]
  18. Mahafza, B.R. Radar Systems Analysis and Design Using MATLAB, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2016; pp. 316–330. [Google Scholar]
  19. Dwivedi, Y.; Rao, S.S. A test for second order stationarity of a time series based on the Discrete Fourier Transform (Technical Report). J. Time 2015, 32, 68–91. [Google Scholar] [CrossRef]
  20. Alabaster, C. Pulse Doppler Processing; IET Digital Library: London, UK, 2012; pp. 35–61. [Google Scholar]
  21. Zhang, J.; Cheng, G.; Tang, J.; Wu, H.; Tian, Z. A Subaperture Motion Compensation Algorithm for Wide-Beam, Multiple-Receiver SAS Systems. J. Mar. Sci. Eng. 2023, 11, 1627. [Google Scholar] [CrossRef]
  22. Ermoshkin, A.V.; Kosteev, D.A.; Ponomarenko, A.A.; Razumov, D.D.; Salin, M.B. Surface Waves Prediction Based on Long-Range Acoustic Backscattering in a Mid-Frequency Range. J. Mar. Sci. Eng. 2022, 10, 722. [Google Scholar] [CrossRef]
  23. Yang, Y.; Fang, S. Improved Velocity Estimation Method for Doppler Sonar Based on Accuracy Evaluation and Selection. J. Mar. Sci. Eng. 2021, 9, 576. [Google Scholar] [CrossRef]
  24. Skrai, K.; Petrovi, J.; Pale, P. Classification of Low- and High-Entropy File Fragments Using Randomness Measures and Discrete Fourier Transform Coefficients. Vietnam J. Comput. Sci. 2023, 10, 433–462. [Google Scholar] [CrossRef]
  25. Li, P.; Wu, Y.; Guo, W.; Cao, C.; Ma, Y.; Li, L.; Leng, H.; Zhou, A.; Song, J. Striation-Based Beamforming with Two-Dimensional Filtering for Suppressing Tonal Interference. J. Mar. Sci. Eng. 2023, 11, 2117. [Google Scholar] [CrossRef]
  26. Varbanets, R.; Fomin, O.; Píštěk, V.; Klymenko, V.; Minchev, D.; Khrulev, A.; Zalozh, V.; Kučera, P. Acoustic Method for Estimation of Marine Low-Speed Engine Turbocharger Parameters. J. Mar. Sci. Eng. 2021, 9, 321. [Google Scholar] [CrossRef]
  27. Liu, S.; Guo, H.; Qian, Z.; Li, J.; Wang, X.; Sun, W.; Zhang, L.; Zhang, A. A Study on a Novel Inverted Ultra-Short Baseline Positioning System and Phase Difference Estimation. J. Mar. Sci. Eng. 2023, 11, 952. [Google Scholar] [CrossRef]
  28. Ye, R.; Xing, H.; Zhou, X. Sea-Surface Small Target Detection Based on Improved Markov Transition Fields. J. Mar. Sci. Eng. 2024, 12, 582. [Google Scholar] [CrossRef]
  29. Zhang, F.; Zhang, W.; Cheng, C.; Hou, X.; Cao, C. Detection of Small Objects in Side-Scan Sonar Images Using an Enhanced YOLOv7-Based Approach. J. Mar. Sci. Eng. 2023, 11, 2155. [Google Scholar] [CrossRef]
  30. Zhao, D.; Xing, H.; Wang, H.; Zhang, H.; Liang, X.; Li, H. Sea-Surface Small Target Detection Based on Four Features Extracted by FAST Algorithm. J. Mar. Sci. Eng. 2023, 11, 339. [Google Scholar] [CrossRef]
  31. Foss, T.; Joreskog, K.G.; Olsson, U.H. Testing structural equation models: The effect of kurtosis. Comput. Stat. Data Anal. 2011, 55, 2263–2275. [Google Scholar] [CrossRef]
  32. Zhang, X.; Yang, P. Back Projection Algorithm for Multi-Receiver Synthetic Aperture Sonar Based on Two Interpolators. J. Mar. Sci. Eng. 2022, 10, 718. [Google Scholar] [CrossRef]
  33. Flamarion, M.V.; Pelinovsky, E. Evolution and Statistical Analysis of Internal Random Wave Fields within the Benjamin–Ono Equation. J. Mar. Sci. Eng. 2023, 11, 1853. [Google Scholar] [CrossRef]
  34. Dai, Z.; Li, H.; Wang, C.; He, Y. Backscattering Statistics of Labeled Sentinel-1 Wave Mode Imagettes for Ten Geophysical Phenomena. J. Mar. Sci. Eng. 2022, 10, 1594. [Google Scholar] [CrossRef]
Figure 1. Typical modality of ACR. (a) Structure of Octahedral Reflector. (b) Parachute stagnation assist.
Figure 1. Typical modality of ACR. (a) Structure of Octahedral Reflector. (b) Parachute stagnation assist.
Jmse 12 00727 g001
Figure 2. Pulsed Doppler radar workflow diagram.
Figure 2. Pulsed Doppler radar workflow diagram.
Jmse 12 00727 g002
Figure 3. Processing of pulsed Doppler data.
Figure 3. Processing of pulsed Doppler data.
Jmse 12 00727 g003
Figure 4. Fast-time–slow-time data graph.
Figure 4. Fast-time–slow-time data graph.
Jmse 12 00727 g004
Figure 5. RD data graph.
Figure 5. RD data graph.
Jmse 12 00727 g005
Figure 6. Sequential Doppler spectrogram.
Figure 6. Sequential Doppler spectrogram.
Jmse 12 00727 g006
Figure 7. Target false dismissal caused by Doppler spectrum broadening.
Figure 7. Target false dismissal caused by Doppler spectrum broadening.
Jmse 12 00727 g007
Figure 8. Time-sequential kurtosis diagram.
Figure 8. Time-sequential kurtosis diagram.
Jmse 12 00727 g008
Figure 9. Time-sequential false dismissal judgement diagram.
Figure 9. Time-sequential false dismissal judgement diagram.
Jmse 12 00727 g009
Figure 10. Doppler broadening phenomenon under different parameter schemes. (a) Doppler broadening severity for different parameter values. (b) Doppler broadening frequency for different parameter values.
Figure 10. Doppler broadening phenomenon under different parameter schemes. (a) Doppler broadening severity for different parameter values. (b) Doppler broadening frequency for different parameter values.
Jmse 12 00727 g010
Figure 11. Broadening spectrum and its original slow-time series. (a) Amplitude fluctuation in the original slow-time series. (b) Doppler spectrum with broadening phenomena.
Figure 11. Broadening spectrum and its original slow-time series. (a) Amplitude fluctuation in the original slow-time series. (b) Doppler spectrum with broadening phenomena.
Jmse 12 00727 g011
Figure 12. Simplified rectangular window for slow-time series amplitude.
Figure 12. Simplified rectangular window for slow-time series amplitude.
Jmse 12 00727 g012
Figure 13. Doppler spectrum of slow-time series with rectangular window.
Figure 13. Doppler spectrum of slow-time series with rectangular window.
Jmse 12 00727 g013
Figure 14. Omnidirectional RCS of ACR.
Figure 14. Omnidirectional RCS of ACR.
Jmse 12 00727 g014
Table 1. Selected values of ACR condition parameters.
Table 1. Selected values of ACR condition parameters.
No.Size (m)Rotation Rate (rad/s)Dropping Velocity (m/s)Height (m)Range (m)
10.20.10.220200
20.40.20.440400
30.60.30.660600
40.80.40.880800
51.00.51.01001000
61.20.61.21201200
71.40.71.41401400
81.60.81.61601600
91.80.91.81801800
102.01.02.02002000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wu, L.; Hu, S.; Feng, C.; Luo, Y.; Liu, Z.; Lin, L. Dynamic Doppler Characteristics of Maritime Airborne Corner Reflector. J. Mar. Sci. Eng. 2024, 12, 727. https://doi.org/10.3390/jmse12050727

AMA Style

Wu L, Hu S, Feng C, Luo Y, Liu Z, Lin L. Dynamic Doppler Characteristics of Maritime Airborne Corner Reflector. Journal of Marine Science and Engineering. 2024; 12(5):727. https://doi.org/10.3390/jmse12050727

Chicago/Turabian Style

Wu, Lingang, Shengliang Hu, Chengxu Feng, Yasong Luo, Zhong Liu, and Li Lin. 2024. "Dynamic Doppler Characteristics of Maritime Airborne Corner Reflector" Journal of Marine Science and Engineering 12, no. 5: 727. https://doi.org/10.3390/jmse12050727

APA Style

Wu, L., Hu, S., Feng, C., Luo, Y., Liu, Z., & Lin, L. (2024). Dynamic Doppler Characteristics of Maritime Airborne Corner Reflector. Journal of Marine Science and Engineering, 12(5), 727. https://doi.org/10.3390/jmse12050727

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop