[go: up one dir, main page]

Next Article in Journal
Mitigation of Suppressive Interference in AMPC SAR Based on Digital Beamforming
Previous Article in Journal
Interrupted-Sampling Repeater Jamming Countermeasure Based on Intrapulse Frequency–Coded Joint Frequency Modulation Slope Agile Waveform
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigating Intra-Pulse Doppler Frequency Coupled in the Radar Echo Signal of a Plasma Sheath-Enveloped Target

1
Key Laboratory of Equipment Efficiency in Extreme Environment, Ministry of Education, Xidian University, Xi’an 710071, China
2
School of Aerospace Science and Technology, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(15), 2811; https://doi.org/10.3390/rs16152811
Submission received: 22 May 2024 / Revised: 10 July 2024 / Accepted: 29 July 2024 / Published: 31 July 2024
(This article belongs to the Section AI Remote Sensing)
Graphical abstract
">
Figure 1
<p>(<b>a</b>) PSh-enveloped target; (<b>b</b>) Spatial distribution characteristics of corresponding parameters of the plasma sheath.</p> ">
Figure 2
<p>Diagram of the transmitted signal.</p> ">
Figure 3
<p>Diagram of the echo signal.</p> ">
Figure 4
<p>Flow chart of intra-pulse Doppler frequency extraction.</p> ">
Figure 5
<p>Schematic diagram demonstrating computation results of the XOY slice of the flow field.</p> ">
Figure 6
<p>Extraction process of I-D frequency.</p> ">
Figure 7
<p>(<b>a</b>) Time-domain waveform of the transmitted signal; Echo signal at the ranges of (<b>b</b>) 80 km, (<b>c</b>) 100 km, and (<b>d</b>) 120 km.</p> ">
Figure 8
<p>(<b>a</b>) Echo signal; (<b>b</b>) Frequency-domain echo signal; (<b>c</b>) Result of the FRFT of the echo signal; (<b>d</b>) Wigner–Ville transform result of the echo signal.</p> ">
Figure 9
<p>(<b>a</b>) Time-frequency-domain of the echo after dechirp; (<b>b</b>) Frequency-domain of the echo after dechirp.</p> ">
Figure 10
<p>Electron density distribution in simulated flow field at different flight altitudes: (<b>a</b>) H = 30 km; (<b>b</b>) H = 50 km; (<b>c</b>) H = 60 km; (<b>d</b>) H = 70 km.</p> ">
Figure 11
<p>Velocity distribution of plasma sheath at different flight altitudes: (<b>a</b>) H = 30 km; (<b>b</b>) H = 50 km; (<b>c</b>) H = 60 km; (<b>d</b>) H = 70 km.</p> ">
Figure 11 Cont.
<p>Velocity distribution of plasma sheath at different flight altitudes: (<b>a</b>) H = 30 km; (<b>b</b>) H = 50 km; (<b>c</b>) H = 60 km; (<b>d</b>) H = 70 km.</p> ">
Figure 12
<p>Reflection coefficient of plasma sheath to electromagnetic wave at an altitude of 30 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Figure 13
<p>The Doppler spectrum of the echo at an altitude of 30 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Figure 14
<p>The Doppler spectrum of the echo at an altitude of 50 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Figure 14 Cont.
<p>The Doppler spectrum of the echo at an altitude of 50 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Figure 15
<p>The Doppler spectrum of the echo at an altitude of 60 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Figure 16
<p>The Doppler spectrum of the echo at an altitude of 70 km: (<b>a</b>) <span class="html-italic">f</span><sub>c</sub> = 3.3 GHz; (<b>b</b>) <span class="html-italic">f</span><sub>c</sub> = 5.8 GHz; (<b>c</b>) <span class="html-italic">f</span><sub>c</sub> = 9.5 GHz; (<b>d</b>) <span class="html-italic">f</span><sub>c</sub> = 18 GHz.</p> ">
Versions Notes

Abstract

:
In detecting hypersonic vehicles, the radar echo signal is coupled with an intra-pulse Doppler frequency (I-D frequency) component caused by relative motion of a plasma sheath (PSh) and the vehicle, which can induce the phenomenon of a ghost target in a one-dimensional range profile. In order to investigate the I-D frequency generated by the relative motion of a PSh, this study transforms a linear frequency modulated (LFM) signal into a single carrier frequency signal based on echo signal equivalent time delay-dechirp processing and realizes high resolution and fast extraction of the I-D frequency coupled with the frequency-domain echo signal. Furthermore, by relying on the computation of the surface flow field of the RAMC-II Blunt Cone Reentry Vehicle, the coupled I-D frequency in the radar echo signal of a PSh-enveloped target under circumstances of typical altitudes and carrier frequencies is extracted and further investigated, revealing the variation law of I-D frequency. The key findings of this study provide a novel approach for suppressing anomalies in radar detection of PSh-enveloped targets as well as effective detecting and as robust target tracking.

Graphical Abstract">

Graphical Abstract

1. Introduction

During high-velocity motion of near-space hypersonic vehicles, the high-velocity incoming flow causes energy conversion, during which enormous kinetic energy is transformed into internal energy that causes high temperature. According to the “real effect of high-temperature gas”, the vehicle surface is enveloped in a PSh [1,2], which changes the target electromagnetic (EM) characteristics of hypersonic vehicles. When the EM wave propagates in the PSh, charged particles and various ablation products will cause certain complicated EM effects on EM waves, resulting in a poor EM environment of the target and severe distortion of the EM wave [3,4].
The existing research on PSh-enveloped targets primarily attached great importance to the impact of a exerted on EM wave transmission characteristics, analyzing amplitude attenuation and phase distortion of EM waves induced by the PSh [5,6,7,8,9,10,11]. Researchers have investigated plasma ground generation devices to analyze the impact mechanism of plasma on EM waves based on ground tests. As various hypersonic vehicles emerge, robust radar detection on PSh-enveloped targets has gradually been carried out [2,12,13,14,15]. Ding et al. constructed a time-domain radar echo model of a PSh-enveloped target and jointly analyzed the characteristics of a radar echo signal in multiple domains, revealing the cause of ghost targets [16]. Later, based on parallel compensation processing of I-D frequency, they improved the gain of the energy ratio of real target to ghost target and effectively suppressed the ghost target [17].
Owing to the fact that a linear frequency modulated (LFM) signal belongs to the scope of a large time-bandwidth product signal, the effective extraction of I-D frequency remains under the influence of signal bandwidth. The relative motion of a PSh with fluid characteristics can cause coupling amongst different I-D frequency components in radar echo signals that are relatively close. The existing I-D frequency extraction methods for a LFM echo signal usually adopt time-frequency analytical approaches and use time-frequency processing to extract the frequency components of a coupled radar echo signal [18,19,20,21]. Nevertheless, the existing time-frequency approaches have obvious drawbacks. When the Wigner–Ville distribution (WVD) transform is used to extract I-D frequency, it is susceptible to a multi-component signal, and the cross term is formed in the time-frequency domain, causing errors in I-D frequency extraction results. To solve the above problem, some scholars proposed pseudo WVD transform and smooth WVD transform. Despite the fact that the aforementioned approaches can eliminate the influence of the cross term to some degree, the frequency resolution is significantly reduced; for this reason, the relatively close I-D frequency cannot be effectively extracted. Although the short-time Fourier transform (SFT) is capable of frequency analysis under variable time scales, the frequency resolution is very limited. The fractional Fourier transform (FT) has low frequency resolution and can extract a close I-D frequency. However, it requires large amount of computation while consuming more computing resources, which hinders it from being applied in real-time extraction of I-D frequency. During the process of suppressing a ghost target of an echo signal, the abovementioned deficiencies not only cause the degeneration of the extraction accuracy of the I-D frequency but also affect the suppression effect. Additionally, there are few studies on the coupled I-D frequency component in the echo of a PSh-enveloped target, and the influence law of flight altitude and carrier frequency on I-D frequency remains to be obtained.
In this study, by relying on equivalent echo time delay estimation, we transformed the problem of frequency extraction into one of estimating time-domain signal parameters. Based on equivalent time delay-deviation processing, the influence of bandwidth is eliminated, and the I-D frequency high-resolution extraction of the echo signal can be effectively conducted. Furthermore, by adopting our proposed approach, the variation law regarding I-D frequency is revealed according to the computation with respect to the flow field. The key contributions of this study are summarized as follows.
Firstly, by relying on the technology of estimating pulse parameters, the time delay of the radar echo signal is equivalent. By performing dechirp-Fourier transform, we realized high-resolution extraction of the coupled I-D frequency component in a radar echo signal.
Secondly, according to the computation results of the surface flow field of the RAMC-II Blunt Cone Reentry Vehicle, the spatial distribution characteristics of the PSh velocity field along the region from the stagnation point to the tail end are analyzed. Our proposed method is capable of extracting the I-D frequency under typical altitude and carrier frequency, thereby revealing the variation law of the coupled I-D frequency in the echo signal of a PSh-enveloped target (different flight altitudes and carrier frequencies). The key findings of this study provide a novel approach for suppressing anomalies in radar detection of PSh-enveloped targets, in effective target detection as well as in robust target tracking.
The remainder of this paper is organized as follows. Section 2 outlines the extraction process of the I-D frequency coupled in the echo signal of the PSh-enveloped target. Section 3 verifies the feasibility of our proposal via simulation, revealing the variation law of I-D frequency under circumstances of typical flight altitudes and carrier frequencies. Section 4 summarizes the simulation and theoretical research of this study.

2. I-D Frequency Extraction of PSh-Enveloped Target Echo Signal

2.1. Radar Echo Model of the Target Enveloped in a PSh

When the velocity of a hypersonic vehicle exceeds Mach 10, its surface will be enveloped in a PSh, which is demonstrated as follows.
Assume the radar signal as a far-field plane wave when the PSh-enveloped target locates at far field. The echo signals are generated at different positions of the target reference points (the regions of the stagnation point, the middle part, and the tail end). Figure 1b shows that L reference points are selected on the vehicle surface. Assume the electron density of the i-th layer at the l-th reference point is Ne(l,i). Due to the spatial distribution of the characteristic parameters involving the electron density Ne(l,i), the collision frequency Ve(l,i), and the velocity V(l,i) of the PSh, the plasma’s multi-dimensional modulation effect exerted on the echo signals differs at each position [16]. Therefore, the characteristic frequency fp(l,i) of this layer can be described by
f p ( l , i ) = 1 2 π N e ( l , i ) · e 2 m e · ε 0 ,
where e is the unit charge, me is the electron mass, and ε0 is the vacuum permittivity. Correspondingly, the dielectric constant ε(l,i) can be expressed as
ε ( l , i ) = ( 1 2 π f p 2 ( l , i ) ( 2 π f c ) 2 + v e 2 ( l , i ) j 4 π 2 f p 2 ( l , i ) v e ( l , i ) / ( 2 π f c ) ( 2 π f c ) 2 + v e 2 ( l , i ) ) · ε 0
where fc denotes the carrier frequency of EM wave, and j denotes the unit vector of the imaginary part.
By adopting the above Equation (2) and μ0 (vacuum permeability), the propagation constant k(l,i) and the intrinsic wave impedance Z(l,i) of the i-th layer at the l-th position of plasma can be obtained, for which the specific expression can be described by
k ( l , i ) = 2 π f c μ 0 ε ( l , i )
Z ( l , i ) = μ 0 / ε ( l , i ) .
By comparing the carrier frequency fc of the radar signal with the characteristic frequency fp(l,i) at each position, the reflection coefficient at different positions and the coupled velocity components in the radar echo signal can be determined according to the above Equations (3) and (4).
When both the characteristic frequency fp(l,i) of the plasma sheath at the l-th reference point of the i-th layer and the carrier frequency fc of the radar echo signal altogether satisfy the following three conditions
f c > f p ( l , i 1 ) f c f p ( l , i ) f c < f p ( l , i + 1 ) ,
the incident depth of the EM wave at the l-th reference position of the PSh is i, the reflection coefficient of the i-th layer is taken as the reflection coefficient at the reference point position, and the flow-field velocity of the i-th layer is the velocity of the reflected wave coupling at the reference point.
The reflection coefficient R(l) and the coupled velocity v(l) of the echo signal wave at the l-th reference position can be expressed as
R ( l ) = A ( l , i ) + B ( l , i ) / Z ( l , i + 1 ) Z 0 C ( l , i ) + D ( l , i ) / Z ( l , i + 1 ) A ( l , i ) + B ( l , i ) / Z ( l , i + 1 ) + Z 0 C ( l , i ) + D ( l , i ) / Z ( l , i + 1 )
v ( l ) = V ( l , i ) ,
where Z0 denotes the wave impedance in vacuum; Z(l, i + 1) denotes the wave impedance at the l-th reference point of the (i + 1)-th layer of the PSh; and A(l), B(l), C(l), and D(l) denote the transmission parameters at the l-th reference point, which are expressed as
A ( l ) B ( l ) C ( l ) D ( l ) = A 1 ( l ) B 1 ( l ) C 1 ( l ) D 1 ( l ) A i 1 ( l ) B i 1 ( l ) C i 1 ( l ) D i 1 ( l ) A i ( l ) B i ( l ) C i ( l ) D i ( l ) .
In the above Equation (8), the symbol “[.]” denotes the transmission matrix of the i-th layer, which can be rewritten as
A i ( l ) B i ( l ) C i ( l ) D i ( l ) = cosh ( j k ( l , i ) d i ) Z ( l , i ) · sinh ( j k ( l , i ) d i ) sinh ( j k ( l , i ) d i ) / Z ( l , i ) cosh ( j k ( l , i ) d i ) ,
where di denotes the thickness of the i-th layer of the PSh.
When both the characteristic frequency fp(l,i) of the PSh at the l-th reference point of the i-th layer and the carrier frequency fc of the radar signal altogether satisfy the following condition
f c > { f p ( l , i ) i = 1 , 2 I } ,
it means that the EM wave is capable of penetrating into the PSh at this position and will be reflected on the target surface. The above Equations (6) and (7) can be simplified into
R ( l ) = B ( l , i ) Z 0 D ( l , i ) B ( l , i ) + Z 0 D ( l , i )
v ( l ) = V 0 ,
where V0 denotes the actual moving velocity of the target (i.e., the radial velocity).
The LFM pulse signal model transmitted by the radar can be described by [22,23,24]
S ( t ) = A r e c t ( t T p ) exp ( j 2 π ( f c t + u 2 t 2 ) ) ,
where “rect(.)” denotes the gate function, the symbol “A” denotes the radar signal amplitude, u denotes the modulation frequency, and Tp denotes the pulse width.
According to Table 1 presented in the following Section 3.1, the reflection coefficient and the I-D frequency at different points of the PSh-enveloped target are computed, and the radar echo model in time domain is established. It is worth noting that since this paper mainly discusses narrowband radar, and the pulse width is not very large, we adopt the go-stop model when building the radar echo model, which means that the time delay is approximately a fixed value in the pulse time. Therefore, the radar echo model in the time domain can be described as
R p l a ( t ) = l = 1 L R ( l ) · S ( f c + f d ( l ) , t τ 0 ) .
In the above Equation (14), τ0 denotes the time delay, and the term fd(l) = 2v(l) × fc/c denotes the I-D frequency corresponding to the coupling velocity of the radar echo signal at the position of the l-th reference point (c denotes the speed of light). The above Equation (14) can be rewritten as [19]
R p l a ( t ) = l = 1 L R ( l ) · r e c t ( t τ 0 T p ) exp ( j π u ( t τ 0 ) 2 ) · ϕ 0 ( l ) ,
where φ0(l) represents the phase term that includes the I-D frequency, which can be expressed as
φ 0 ( l ) = exp ( j 2 π ( f c + f d ( l ) ) τ 0 ) · exp ( j 2 π f d ( l ) t ) .
It can be deduced from the above Equations (15) to (16) that the relative motion of the PSh causes the echo signals to couple with different I-D frequency components. The I-D frequency component causes abnormal processing of the matched filter, thus causing a “ghost target” phenomenon on the one-dimensional range profile.

2.2. Equivalent Time Delay Derivation

When the PSh-enveloped target is detected by radar, the target-to-radar range causes a time delay, which further leads to the envelope offset of the radar signal on the time axis.
The transmitted radar signal is demonstrated as follows.
When the transmitted signal is reflected on the target surface, its echo will be coupled with the time delay that is caused by the target-to-radar range, thereby inducing the envelope offset phenomenon. The schematic diagram demonstrating the echo signal is presented in the following figure.
It can be observed that the range causes an envelope offset of the echo signal. In Figure 2, the start time of the transmitting signal is 0, and the end time is the pulse width Tp of the signal. In Figure 3, the envelope of the echo is offset to the right side along the time axis, and its start time correlates with the delay of the echo signal coupling. Combining Equations (13) and (14), we can further conclude that the envelope offset of the echo is τ0, while its time delay τ0 can be equivalent by extracting its start time.
The time delay τ0 corresponding to the echo’s start time can be calculated by
τ 0 = T O A = 2 × R 0 c ,
where R0 denotes the target-to-radar range.

2.3. Estimating Pulse Parameter of the Echo

It has been discussed in Section 2.2 that the time delay can be equivalently obtained by estimating the echo signal’s start time. In this subsection, we estimate the echo signal’s start/end time according to the difference of the boxes (DOB) filter.
Let the length of the DOB filter h(n) be N (the order is 2N + 1), so that h(n) can be described by
h ( n ) = 1 n [ N , 1 ] 0 n = 0 1 n [ 1 , N ] .
Based on Equation (18), convolution processing is conducted with respect to the echo signal to obtain the output signal y(n), which can be calculated by
y ( n ) = R p l a ( t ) * h ( n ) = x R p l a ( t ) × h ( n + 1 x ) .
By substituting Equation (15) into Equation (19), the calculation process can be subdivided into the following seven situations:
  • ① When 1 n N , Equation (19) can be rewritten as
    y n = p = 1 n R p l a p ;
  • ② When n = N + 1 , Equation (19) can be rewritten as
    y n = p = 2 N + 1 R p l a p ;
  • ③ When N + 2 n 2 N + 1 , Equation (19) can be rewritten as
    y n = p = n N + 1 n R p l a p q = 1 n N 1 R p l a q ;
  • ④ When 2 N + 2 n M 1 , Equation (19) can be converted into
    y n = p = n N + 1 n R p l a p q = n 2 N n N 1 R p l a q ;
  • ⑤ When M n M + N 1 , Equation (19) can be rewritten as
    y n = p = n N + 1 M R p l a p q = n 2 N n N 1 R p l a q ;
  • ⑥ When n = M + N , Equation (19) can be rewritten as
    y n = p = M N M 1 R p l a p ;
  • ⑦ When M + N + 1 n M + 2 N , Equation (19) can be rewritten as
    y n = p = n 2 N M R p l a p .
Combining Equations (20)–(26), the recursive expression of y(n) can be described by
y n + 1 = y n + R p l a n ,   1 n N 1 y n R p l a n N + R p l a n N + 1 + R p l a n + 1 ,   N + 1 n 2 N y n + R p l a n 2 N + 1 R p l a n N R p l a n N + 1 + R p l a n + 1 ,   2 N + 1 n M 2 y n + R p l a n 2 N R p l a n N R p l a n N + 1 ,   M 1 n M + N 2 y n + R p l a n 2 N ,   M + N 1 n M + 2 N 1 ,
where y(n) is updated by sliding recursion, forming a maximum peak and a minimum peak at the echo’s start time and end time, respectively.
By extracting the peak position, the start time TOA and the end time TOE of the echo signal can be estimated by
T O A = n max N F s T O E = n min N F s ,
where nmax denotes the coordinates of the maximum value, and nmin denotes the coordinates of the minimum value.

2.4. Extracting I-D Frequency

By adopting the echo signal’s prior parameters (modulation frequency μ, pulse width Tp, and Sampling frequency Fs), the dechirp signal D(t) can be constructed, which is expressed as
D ( t ) = exp ( j π μ ( 0 : 1 / F s : T p ) 2 ) .
According to the estimation result of TOA and TOE described in Equation (28), when substituting the dechirp Equation (29) to the echo signal’s start time, we obtain
D Z ( T O A : T O E ) = D ( t ) .
When the echo signal is processed according to Equation (30), then we obtain
R D ( t ) = R p l a ( t ) × D z ( t ) .
The above Equation (31) can further be rewritten as
R D ( t ) = l = 1 L R e f ( n ) · r e c t ( t τ 0 T p ) · φ 0 ( l ) .
When we substitute Equation (16) into the above Equation (32), we obtain
R D ( t ) = l = 1 L R ( l ) · r e c t ( t τ 0 T p ) exp ( j 2 π ( f c + f d ( l ) ) τ 0 ) · exp ( j 2 π f d ( l ) t ) .
From Equation (33), we observe that the echo signal after being dechirp-processed does not contain the time quadratic term, and the coefficient with the time term is the I-D frequency.
When we perform FT with respect to Equation (33), we obtain
R D ( f ) = t = 1 T [ R D ( t ) ] exp ( j 2 π f ( t 1 ) T ) ,
where T denotes the total sampling point numbers of the echo signal. By substituting Equation (33) into Equation (34), the echo signal represented in the frequency domain after being dechirp-processed can be rewritten as
R D ( f ) = l = 1 L R ( l ) s a ( f f d ( l ) ) exp ( j 2 π f τ 0 ) exp ( j 2 π f d ( l ) τ 0 ) exp ( j 2 π f c τ 0 ) ,
where the sa(.) function produces a frequency-domain peak, and different I-D frequencies generate different frequency-domain peaks. The peak coordinates corresponding to fd(l) can be obtained through the peak search, and the frequency corresponding to the peak coordinates can be solved to realize the extraction of I-D frequencies.
For the abovementioned Equations (17)–(35), the following Figure 4 demonstrates the flow chart of the intra-pulse Doppler frequency extraction method proposed in this study.

2.5. Analyzing Variation Law of I-D Frequency

According to the aforementioned computation results, the I-D frequency of the echo signal of a PSh-enveloped target under circumstances of typical altitudes and carrier frequencies is extracted, for which the variation law regarding the I-D frequency is analyzed.
Assume the attack angle of the vehicle in the flow-field calculation is 0 and the PSh on the vehicle surface is symmetric regarding planes passing through the X-axis. The slice processing of the XOY plane is conducted according to the computation results of the of plasma flow field, for which the schematic diagram is presented in Figure 5.
It can be observed that by extracting the data at the position of the reference points, the numerical computation results (electron density, collision frequency, and velocity component) of the characteristic parameters in different areas of the PSh can be obtained, for which the specific process is demonstrated in Figure 6.
The I-D frequency in the echo signal can be extracted via conducting the specific steps presented in the above Figure 6. By using the obtained computation results of the flow field, we investigated the coupled I-D frequency in the echo signal characterized by different flight altitudes and carrier frequencies and revealed the variation law of the I-D frequency.

3. Simulation Analysis of I-D Frequency Extraction

3.1. Pulse Parameter Extraction of Radar Echo Signal

Based on the plasma flow-field data of the blunt cone aircraft at the height of 40 km, the maximum reflection intensity and coupled Doppler frequency of the radar echo are calculated at each reference point after it is perpendicular to the incident plasma flow field on the surface of the aircraft. The plasma flow-field parameters and radar signal parameters are shown in Table 1. By referring to the aforementioned Equations (13) and (15), the transmitted signal and the time-domain radar echo signal are simulated, and a convolution operation is conducted with respect to the echo signal based on the DOB filter to extract the echo’s start time. Considering the target-to-radar range, we also analyzed the relationship between the start time and the time delay.
As shown in Figure 7, the red dotted line indicates results from the echo signal output through the DOB filter. It can be observed that the three target-to-radar ranges differ, and the offset caused by the time-domain echo signal differs accordingly. The transmitted signal possesses 500 μs pulse width, 0 μs start time, and 500 μs end time, respectively.
When the target-to-radar range is 80 km, the echo’s time delay is computed to be 533.3 μs using Equation (17). The DOB filter processes the echo to form two peaks in the time domain. The start time of the echo measured by the method proposed in this paper is 533.3 μs, which is consistent with the theoretical value. The echo’s end time is 1033.4 μs, and its pulse width is estimated to be 20.01 μs throughout the start and end times, in which the error with the actual pulse width is 0.02%. When the target-to-radar range is 100 km, the time delay generated by the target is 666.7 μs. The output maximum peak position based on the DOB filter is 666.7 μs, which is consistent with the time delay caused by the target, in which the error between the echo’s pulse width and the actual pulse width is 0.02%. When the target-to-radar range is 120 km, the target time delay is 800 μs. The DOB output result of the echo is consistent with the actual delay of the target, and the error between the estimated pulse width of the echo and the actual pulse width is 0.02%.
It can be acknowledged that the start time of the echo can be extracted based on DOB filter processing, and the echo’s time delay can be obtained.

3.2. I-D Frequency Extraction of Echo Signals

The existence of the difference in spatial distribution of the PSh suggests that the echo and the multi-I-D frequency components will be coupled. In this section, based on the simulated radar echo signal in Section 3.1, we analyze the inadequate conventional time-frequency processing method (WVD, fractional Fourier transform (FRFT)) to extract the I-D frequency. According to Equations (14), (29), and (35), the simulated radar echo signal is processed, and the frequency extraction method was verified. In this section, the target distance of the simulated echo is 100 km, and other simulation parameters are consistent with those in Section 3.1.
It can be observed from Figure 8a that since the I-D frequencies of the echo (the target echo signal at the same range) differ, the time-domain echo signal is distorted, causing a “periodic dispersion” phenomenon. Additionally, we observe from Figure 8b that a significant modulation phenomenon occurs in the frequency-domain echo signal, and a serious “oscillation” phenomenon appears in the bandwidth.
We observe from Figure 8c that the energy of the FRFT transformation results distributes on a straight line. The FRFT transform of the echo is limited by the frequency resolution. When the I-D frequencies are closely coupled in the echo, energy aliasing will occur, under which circumstance the peak position will expand, preventing the I-D frequency from being extracted effectively.
Furthermore, Figure 8d demonstrates that some slant lines occur in a two-dimensional time-frequency plane of Wigner–Ville transform results of the echo signals, which have the same slope and intercepts (the same modulation frequency and different carrier frequencies). Because of the existence of multiple I-D frequency components in the echo, the time-frequency curve contains many cross terms, causing errors in the extraction results of the I-D frequency.
Figure 9a suggests that when the echo is dechirp-processed, its time-frequency curve is a straight line perpendicular to the frequency axis. According to Equation (32), the dechirp processing can effectively eliminate the modulation frequency u of the echo, for which reason the signal only contains the carrier frequency information. Figure 9b demonstrates that the frequency-domain characteristics of the dechirp-processed echo can be obtained by Fourier transform. Since the modulation frequency of the echo is eliminated, no influence will be exerted on the bandwidth in the frequency domain (the echo signal is transformed from a LFM signal to a single carrier frequency signal). There are several significant peaks in the frequency spectrum, and the peak positions correspond to the signal carrier frequency.
Table 2 shows that under the circumstance of the same parameters, our proposed method is capable of extracting the minimum frequency interval, which is only 0.0657 MHz. Since the I-D frequency is usually low and the intra-pulse frequency components are relatively close to each other, our proposed method can effectively extract the I-D frequency of the coupled echo signal. Additionally, the method exhibits very low computational complexity, and its operation time is much less than the other two methods under the same circumstances. The critical SNR, for which this method can effectively extract the minimum frequency interval, is −9 dB, whereas the critical SNRs of the other two methods are −12 dB and −10 dB, respectively. The critical SNR of our proposed method in this study is a bit higher than the other two methods. One reasonable explanation is that our proposed method needs to estimate the pulse parameters. A SNR that is too low reduces the estimation accuracy of the pulse parameters and causes similar I-D frequencies to be easily aliased due to noise.

3.3. Velocity Distribution Characteristics of the PSh

Based on the simulated flow-field data of the conical target at different heights, this subsection analyzes the distribution characteristics of the velocity of the plasma flow field under circumstances of typical aircraft flight velocity (V = 24 Mach, 1 Mach = 340 m/s) and different flight altitudes (H = 30 km~70 km). The electron density distribution of the simulated flow field of the cone aircraft at different altitudes is shown in Figure 10.
Figure 11a–d respectively represent the velocity distribution characteristics of the plasma flow field at four altitudes. The red dashed line represents the vehicle surface. The plasma is produced within the stagnation point area of the vehicle. Influenced by the action of high-velocity incoming flow, the plasma moves along the stagnation point area to the tail-end area and forms a PSh enveloping the vehicle. As a result, the spatial distribution of the plasma flow velocity varies significantly within different reference areas.
At 30 km altitude, the innermost flow-field velocity of the PSh is similar to the actual velocity of the vehicle. Vertically outward along the vehicle surface, the flow-field velocity declines mildly. The velocity of the outer flow field decreases faster than that of the inner flow field, and the velocity of the outer flow field is far lower than the vehicle velocity. Within the area of the stagnation point, the inner flow-field velocity has an approximately uniform distribution, but the outer flow-field velocity sharply decreases. From the area of the stagnation point to that of the tail end, the velocity distribution interval of the plasma flow field first increases and then decreases, while the velocity descending gradient decreases accordingly.
As shown in Figure 10, with the increase in flight altitude, the electron density of the flow field decreases, but the thickness of the plasma sheath also increases accordingly. At the position close to the vehicle surface, the range of the flow field whose velocity is close to the target velocity gradually increases. The descending gradient of the flow-field velocity within the stagnation point area is approximately constant, the descending gradient of the flow-field velocity in the middle area and the tail-end area decreases.

3.4. Variation Law of the Coupled I-D Frequency in an Echo Signal under Typical Parameters

According to the aforementioned computation results, this subsection analyzes the coupled I-D frequency in the echo signal of the PSh-enveloped target and reveals the variation law of the I-D frequency through different carrier frequencies (fc = 3.3 GHz, 5.8 GHz, 9.5 GHz, and 18 GHz) and different flight altitudes (H = 30 km, 50 km, 60 km, and 70 km).
At the flight altitude of 30 km:
Figure 12. Reflection coefficient of plasma sheath to electromagnetic wave at an altitude of 30 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Figure 12. Reflection coefficient of plasma sheath to electromagnetic wave at an altitude of 30 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Remotesensing 16 02811 g012
Figure 13. The Doppler spectrum of the echo at an altitude of 30 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Figure 13. The Doppler spectrum of the echo at an altitude of 30 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Remotesensing 16 02811 g013
It can be observed that the relative motion of the PSh causes the echo signal to couple with different I-D frequencies. As can be seen from Figure 12, at an altitude of 30 km, with the increase in carrier frequency, the strongest reflection positions of electromagnetic waves in the plasma sheath are all in the outer layer of the PSh. At this time, the plasma sheath has a strong shielding effect on electromagnetic waves. As can be seen from Figure 11a, when the flight altitude of the aircraft is 30 km, the velocity outside the PSh is more than 10 Mach lower than the target flight velocity. Therefore, the Doppler frequency coupled with the outside of the PSh is much lower than the Doppler frequency corresponding to the target velocity. Figure 13a suggests that when the carrier frequency is 3.3 GHz, the maximum I-D frequency coupled in the echo signal is 0.1244 MHz, under which circumstance the corresponding velocity is 16.6 Mach. Because of the extremely high electron density at this flight altitude (the peak characteristic frequency of the plasma flow field within the tail-end area is 23 GHz, which exceeds the radar signal’s carrier frequency), the wave incident depth remains approximately unchanged as the carrier frequency (5.8 GHz, 9.5 GHz, and 18 GHz) increases, under which circumstance the velocity corresponding to the maximum I-D frequency coupled in the echo signal approximates 17 Mach.
At the flight altitude of 50 km:
Figure 14. The Doppler spectrum of the echo at an altitude of 50 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Figure 14. The Doppler spectrum of the echo at an altitude of 50 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Remotesensing 16 02811 g014aRemotesensing 16 02811 g014b
At 50 km altitude, the electron density decreases, the PSh shielding effect imposed on the wave degrades, and the incident depth increases. As shown in Figure 14, when the carrier frequency is 3.3 GHz, the maximum I-D frequency coupled in the echo signal is 0.1595 MHz, for which the corresponding velocity is 21.3 Mach, and the maximum peak energy is 91 dB. When the carrier frequency increases to 5.8 GHz and 9.5 GHz, the maximum I-D frequency coupled in the echo signal increases. This is because the incident depth increases, while the velocity corresponding to the maximum I-D frequency is 22.5 Mach. As the incident depth continues to increase, the secondary energy absorption effect of the PSh exerted on the echo signal intensifies, and the maximum peak energy is reduced by more than 5 dB. As the carrier frequency increases further to 18 GHz, the incident depth continues to increase, and the maximum I-D frequency coupled in the echo signal is 0.9316 MHz, to which the corresponding velocity is 23 Mach. The secondary energy absorption effect reduces the maximum peak energy to 75.6 dB.
At the flight altitude of 60 km:
Figure 15. The Doppler spectrum of the echo at an altitude of 60 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Figure 15. The Doppler spectrum of the echo at an altitude of 60 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Remotesensing 16 02811 g015
When the altitude is 60 km, the electron density of the PSh is further reduced, under which circumstance the PSh shielding effect exerted on the wave is alleviated.
It can be observed from Figure 15a,b that the carrier frequencies are 3.3 GHz and 5.8 GHz, respectively, for which the velocity corresponding to the maximum I-D frequency coupled in the echo signal approximates 21.7 Mach. Furthermore, the secondary energy absorption effect of the PSh imposed on the echo signal lowers the maximum peak energy from 91.5 dB to 82.8 dB. As the carrier frequency increases to 9.5 GHz, the incident depth increases, penetrating through the PSh in some areas (the wave is reflected on the vehicle surface), and the I-D frequency component coupled in the echo signal increases from two to three. An increase in the number of Doppler peaks means an increase in the number of major Doppler components of the coupling, which is due to an increase in strong reflection locations. From Figure 15c, we observe that the maximum I-D frequency coupled in the echo signal is 0.518 MHz, which approximates the Doppler frequency corresponding to the target velocity, and the maximum peak energy is 79.5 dB. As the carrier frequency increases to 18 GHz, the EM wave component penetrating the PSh gradually increases, under which circumstance the energy of the I-D frequency corresponding to the target velocity reaches its maximum value of 86.25 dB.
At a flight altitude of 70 km:
Figure 16. The Doppler spectrum of the echo at an altitude of 70 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Figure 16. The Doppler spectrum of the echo at an altitude of 70 km: (a) fc = 3.3 GHz; (b) fc = 5.8 GHz; (c) fc = 9.5 GHz; (d) fc = 18 GHz.
Remotesensing 16 02811 g016
At the altitude of 70 km, the electron density is far less than that at 30 km, 50 km, and 60 km. The PSh in the vehicle’s middle and tail-end areas exert no shielding effect on the wave, for which reason the wave can penetrate the PSh and then be reflected on the vehicle surface. In the above scenario, the carrier frequency is 3.3 GHz, and the maximum I-D frequency component coupled in the echo signal is 0.1808 MHz, approximating the Doppler frequency corresponding to the target velocity. Under the circumstance of 3.3 GHz carrier frequency, despite the fact that the EM wave penetrates through the PSh, the penetration area is limited, so the secondary absorption effect of the flow field causes the minimum peak energy of the maximum I-D frequency. As the carrier frequency increases to 5.8 GHz and further to 9.5 GHz, the I-D frequency energy corresponding to the target velocity increases from 85 dB to 90 dB. As the carrier frequency increases to 18 GHz, the incident depth increases within the areas of the stagnation point and the near-stagnation point, under which circumstance the echo signal couples more frequency components. At 18 GHz carrier frequency, the wave intensity of the echo signal of the target surface increases, and the I-D frequency energy corresponding to the target velocity in the echo signal increases to 92 dB, which is much higher than other I-D frequency energies.
Table 3. The velocity corresponding to the maximum Doppler frequency.
Table 3. The velocity corresponding to the maximum Doppler frequency.
Maximum Velocityfc = 3.3 GHzfc = 5.8 GHzfc = 9.5 GHzfc = 18 GHz
H = 30 km16.63 Ma16.77 Ma16.68 Ma16.77 Ma
H = 50 km22.32 Ma22.8 Ma22.71 Ma22.83 Ma
H = 60 km21.72 Ma21.53 Ma24.06 Ma23.99 Ma
H = 70 km24.17 Ma23.85 Ma23.92 Ma23.99 Ma
Combining Figure 13, Figure 14, Figure 15 and Figure 16, the simulation results are listed in the Table 3. It can be ascertained that at relatively low flight altitude, due to the extremely strong EM shielding effect of the PSh, the EM wave is reflected outside the PSh, and the I-D frequency coupled in the echo signal is far less than the Doppler frequency corresponding to the target velocity. As the carrier frequency increases, the corresponding velocity of the I-D frequency coupled in the echo signal remains approximately unchanged. As the flight altitude continues to increase, the electron density decreases, the EM shielding effect of the PSh degrades, the incident depth increases, and the I-D frequency component coupled in the echo signal increases. It must be noted that the secondary energy absorption effect of the PSh exerted on the EM wave results in the overall downward trend of the echo signal energy. When the flight altitude increases to 60 km and 70 km, the shielding effect of the PSh sharply degrades. When the carrier frequency remains at no less than 5.8 GHz, the EM wave is capable of penetrating through the PSh and will be reflected on the target surface, and then the echo signal couples with the I-D frequency corresponding to the target velocity. As the carrier frequency increases, the energy absorption effect of the PSh exerted on the reflection wave of the target surface degrades, whereas the energy of the I-D frequency peak corresponding to the target velocity increases slowly.

4. Conclusions

In this study, we proposed an I-D frequency analysis method for the radar echo of a PSh-enveloped target, revealing the variation law of I-D frequency according to the aforementioned numerical computation results cited in the references.
Specifically, we adopted a difference filter to conduct parameter estimation with respect to the echo signal, which is then dealt with based on equivalent time delay theory and dechirp-Fourier transform, thereby realizing high-resolution extraction of I-D frequency. Furthermore, we formulated the spatial distribution characteristics of the velocity field of the PSh along the area of the stagnation point to that of the tail end. By performing our proposed method in this study, we extracted the I-D frequency coupled in the echo signal under circumstances of typical altitudes (H = 30 km, 50 km, 60 km, and 70 km) and typical carrier frequencies (fc = 3.3 GHz, 5.8 GHz, 9.5 GHz, and 18 GHz). We also analyzed the influence exerted by different altitudes and carrier frequencies on the peak position, peak energy, and peak number of I-D frequency components, through which the variation law of I-D frequency was investigated.
Despite the abovementioned research findings, certain positive factors or potential shortcomings must be pointed out concerning related parameter settings, experimental setup, etc., which may facilitate or hinder the practical application of extracting I-D frequency components. Our future research interests will be focusing on suppressing abnormal phenomena in radar detection of PSh-enveloped targets, robust target detection methods, and reliable target tracking, laying a theoretical foundation for improving the generality and robustness of radar detection approaches.

Author Contributions

Conceptualization, B.B. and B.P.; methodology, B.B.; software, B.P.; formal analysis, B.B., B.P., K.Z., Y.Y., X.L. and Y.L.; investigation, B.B., X.L. and Y.L.; writing—original draft preparation, B.P.; writing—review and editing, B.B.; supervision, B.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported partially by the National Natural Science Foundation of China under Grant Nos. 62171349, 92371205, and 62001388 and partially by the Natural Science Basic Research Plan in Shaanxi Province of China under Grant No. 2022KJXX-37.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors would like to thank the editors and the anonymous reviewers for their careful reading and helpful remarks.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hartunian, R.A.; Stewart, G.E.; Fergason, S.D.; Curtiss, T.J.; Seibold, R.W. Causes and Mitigation of Radio Frequency (RF) Blackout during Reentry of Reusable Launch Vehicles; Contractor Rep. ATR-2007(5309)-1; Aerospace Corporation: El Segundo, CA, USA, 2007. [Google Scholar]
  2. Dix, D.M. Typical Values of Plasma Parameters Around a Conical Re-Entry Vehicle; Scientific rept. AD295429; Aerospace Corporation: El Segundo, CA, USA, 1962. [Google Scholar]
  3. Bai, B.; Li, X.; Liu, Y.; Xu, J.; Shi, L.; Xie, K. Effects of hypersonic Plasma Sheath on the Polarization Properties of Obliquely Incident EM Waves. IEEE Trans. Plasma Sci. 2014, 42, 3365–3372. [Google Scholar] [CrossRef]
  4. Sha, Y.X.; Zhang, H.L.; Guo, X.Y.; Xia, M.Y. Analyses of Electromagnetic Properties of a Hypersonic Object With Plasma Sheath. IEEE Trans. Antennas Propag. 2019, 67, 2470–2481. [Google Scholar] [CrossRef]
  5. Rybak, J.P.; Churchill, R.J. Progress in Reentry Communications. IEEE Trans. Aerosp. Electron. Syst. 1971, AES-7, 879–894. [Google Scholar] [CrossRef]
  6. Shi, L.; Liu, Y.; Fang, S.; Li, X.; Yao, B.; Zhao, L.; Yang, M. Adaptive multistate Markov channel modeling method for reentry dynamic plasma sheaths. IEEE Trans. Plasma Sci. 2016, 44, 1083–1093. [Google Scholar] [CrossRef]
  7. Kojima, T.; Higashi, T. Reflection and transmission of electromagnetic waves obliquely incident upon a moving compressible plasma slab. IEEE Trans. Antennas Propag. 1972, 20, 398–400. [Google Scholar] [CrossRef]
  8. Chen, X.Y.; Li, K.X.; Liu, Y.Y.; Zhou, Y.G.; Li, X.P.; Liu, Y.M. Study of the influence of time-varying plasma sheath on radar echo signal. IEEE Trans. Plasma Sci. 2017, 45, 3166–3176. [Google Scholar] [CrossRef]
  9. Laroussi, M.; Roth, J.R. Numerical Calculation of the Reflection, Absorption, and Transmission of Microwaves by a Non-Uniform Plasma Slab. IEEE Trans. Plasma Sci. 1993, 21, 366–372. [Google Scholar] [CrossRef]
  10. Kalluri, D.K.; Goteti, V.R. WKB Solution for Wave Propagation in a Time-Varying Magneto plasma Medium: Longitudinal Propagation. IEEE Trans. Plasma Sci. 1993, 21, 70–76. [Google Scholar] [CrossRef]
  11. Lee, J.H.; Kalluri, D.K. Three-Dimensional FDTD Simulation of Electromagnetic Wave Transformation in a Dynamic Inhomogeneous Magnetized Plasma. IEEE Trans. Antennas Propag. 1999, 47, 1146–1151. [Google Scholar]
  12. Close, S.; Oppenheim, M.; Hunt, S.; Dyrud, L. Scattering characteristics of high-resolution meteor head echoes detected at multiple frequencies. J. Geophys. Res. Space Phys. 2002, 107, SIA-9-1–SIA-9-12. [Google Scholar] [CrossRef]
  13. Luo, C.R.; Ding, C.L.; Duan, L.B. 12th Five-Year Plan Textbooks: Electrodynamics; Publishing House of Electronics Industry: Beijing, China, 2016; pp. 150–176. [Google Scholar]
  14. Cao, C.Q. Basic Series of Modern Physics: Classical Electrodynamics; Science Press: Beijing, China, 2009; pp. 167–180. [Google Scholar]
  15. Kero, J.; Szasz, C.; Wannberg, G. On the meteoric head echo radar cross section angular dependence. Geophys. Res. Lett. 2008, 35, 154–162. [Google Scholar] [CrossRef]
  16. Ding, Y.; Bai, B.; Gao, H.; Niu, G.; Shen, F.; Liu, Y.; Li, X. An Analysis of Radar Detection on a Plasma Sheath Covered Reentry Target. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 4255–4268. [Google Scholar] [CrossRef]
  17. Ding, Y.; Bai, B.; Gao, H.; Liu, Y.; Li, X.; Zhao, M. Method of Detecting a Target Enveloped by a Plasma Sheath Based on Doppler Frequency Compensation. IEEE Trans. Plasma Sci. 2020, 48, 4103–4111. [Google Scholar] [CrossRef]
  18. Mudukutore, A.S.; Chandrasekar, V.; Keeler, R.J. Pulse compression for weather radars. IEEE Trans. Geosci. Remote Sens. 1998, 36, 125–142. [Google Scholar] [CrossRef]
  19. Su, J.; Xing, M.; Wang, G. High-velocity multi-target detection with narrow band radar. IET Radar Sonar Navig. 2010, 4, 595–603. [Google Scholar] [CrossRef]
  20. Huang, P.; Liao, G.; Yang, Z.; Xia, X.G.; Ma, J.T.; Zhang, X. A Fast SAR Imaging Method for Ground Moving Target Using a Second-Order WVD Transform. IEEE Trans. Geosci. Remote Sens. 2016, 54, 1940–1956. [Google Scholar] [CrossRef]
  21. Stankovic, L. On the realization of the polynomial Wigner-Ville distribution for multi-component signals. IEEE Signal Process. Lett. 1998, 5, 157–159. [Google Scholar] [CrossRef]
  22. Zhan, M.; Huang, P.; Liu, X.; Liao, G.; Zhang, Z.; Wang, Z.; Hou, Q. Space Maneuvering Target Integration Detection and Parameter Estimation for a Spaceborne Radar System with Target Doppler Aliasing. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 3579–3594. [Google Scholar] [CrossRef]
  23. Li, X.; Sun, Z.; Yi, W.; Cui, G.; Kong, L. Radar Detection and Parameter Estimation of High-velocity Target Based on MART-LVT. IEEE Sens. J. 2019, 19, 1478–1486. [Google Scholar] [CrossRef]
  24. Chen, J.; He, C.; Liang, X.; Bai, L.; Bai, X.; Zhu, W.; Jin, R. Direction Finding of Linear Frequency Modulation Signal in Time Modulated Array with Pulse Compression. IEEE Trans. Antennas Propag. 2020, 68, 509–520. [Google Scholar] [CrossRef]
Figure 1. (a) PSh-enveloped target; (b) Spatial distribution characteristics of corresponding parameters of the plasma sheath.
Figure 1. (a) PSh-enveloped target; (b) Spatial distribution characteristics of corresponding parameters of the plasma sheath.
Remotesensing 16 02811 g001
Figure 2. Diagram of the transmitted signal.
Figure 2. Diagram of the transmitted signal.
Remotesensing 16 02811 g002
Figure 3. Diagram of the echo signal.
Figure 3. Diagram of the echo signal.
Remotesensing 16 02811 g003
Figure 4. Flow chart of intra-pulse Doppler frequency extraction.
Figure 4. Flow chart of intra-pulse Doppler frequency extraction.
Remotesensing 16 02811 g004
Figure 5. Schematic diagram demonstrating computation results of the XOY slice of the flow field.
Figure 5. Schematic diagram demonstrating computation results of the XOY slice of the flow field.
Remotesensing 16 02811 g005
Figure 6. Extraction process of I-D frequency.
Figure 6. Extraction process of I-D frequency.
Remotesensing 16 02811 g006
Figure 7. (a) Time-domain waveform of the transmitted signal; Echo signal at the ranges of (b) 80 km, (c) 100 km, and (d) 120 km.
Figure 7. (a) Time-domain waveform of the transmitted signal; Echo signal at the ranges of (b) 80 km, (c) 100 km, and (d) 120 km.
Remotesensing 16 02811 g007
Figure 8. (a) Echo signal; (b) Frequency-domain echo signal; (c) Result of the FRFT of the echo signal; (d) Wigner–Ville transform result of the echo signal.
Figure 8. (a) Echo signal; (b) Frequency-domain echo signal; (c) Result of the FRFT of the echo signal; (d) Wigner–Ville transform result of the echo signal.
Remotesensing 16 02811 g008
Figure 9. (a) Time-frequency-domain of the echo after dechirp; (b) Frequency-domain of the echo after dechirp.
Figure 9. (a) Time-frequency-domain of the echo after dechirp; (b) Frequency-domain of the echo after dechirp.
Remotesensing 16 02811 g009
Figure 10. Electron density distribution in simulated flow field at different flight altitudes: (a) H = 30 km; (b) H = 50 km; (c) H = 60 km; (d) H = 70 km.
Figure 10. Electron density distribution in simulated flow field at different flight altitudes: (a) H = 30 km; (b) H = 50 km; (c) H = 60 km; (d) H = 70 km.
Remotesensing 16 02811 g010
Figure 11. Velocity distribution of plasma sheath at different flight altitudes: (a) H = 30 km; (b) H = 50 km; (c) H = 60 km; (d) H = 70 km.
Figure 11. Velocity distribution of plasma sheath at different flight altitudes: (a) H = 30 km; (b) H = 50 km; (c) H = 60 km; (d) H = 70 km.
Remotesensing 16 02811 g011aRemotesensing 16 02811 g011b
Table 1. Simulation parameters of flow field and echo.
Table 1. Simulation parameters of flow field and echo.
Parameter TypeParameter SettingValues
Parameters of target and flow fieldMaximum plasma densityNemax = 2.67 × 1020/m3
HeightH = 40 km
Shape of aircraftBlunt cone
Velocity of vehicleVtar = 4521 m/s
Radar signal parametersCarrier frequencyfc = 3.3 GHz
BandwidthB = 5 MHz
Pulse widthTp = 500 μs
Distance of targetR0 = 80 km/100 km/120 km
Table 2. Performance of the proposed method.
Table 2. Performance of the proposed method.
FRFTWVDProposed Method
Minimum frequency intervalΔf = 0.14 MHzΔf = 0.183 MHzΔf = 0.0657 MHz
Computation timeT = 6.5 sT = 7.8 sT = 1.4 s
Critical SNRSNR = −12 dBSNR = −10 dBSNR = −9 dB
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

Bai, B.; Pu, B.; Zhang, K.; Yang, Y.; Li, X.; Liu, Y. Investigating Intra-Pulse Doppler Frequency Coupled in the Radar Echo Signal of a Plasma Sheath-Enveloped Target. Remote Sens. 2024, 16, 2811. https://doi.org/10.3390/rs16152811

AMA Style

Bai B, Pu B, Zhang K, Yang Y, Li X, Liu Y. Investigating Intra-Pulse Doppler Frequency Coupled in the Radar Echo Signal of a Plasma Sheath-Enveloped Target. Remote Sensing. 2024; 16(15):2811. https://doi.org/10.3390/rs16152811

Chicago/Turabian Style

Bai, Bowen, Bailiang Pu, Ke Zhang, Yilin Yang, Xiaoping Li, and Yanming Liu. 2024. "Investigating Intra-Pulse Doppler Frequency Coupled in the Radar Echo Signal of a Plasma Sheath-Enveloped Target" Remote Sensing 16, no. 15: 2811. https://doi.org/10.3390/rs16152811

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