Simulation and Feature Extraction of the Dynamic Electromagnetic Scattering of a Hypersonic Vehicle Covered with Plasma Sheath

Based on the detection platform of space-borne radar (SBR), the dynamic electromagnetic (EM) scattering characteristics of a hypersonic vehicle covered with plasma sheath are calculated and analyzed by the physical optics (PO) method. Using the angle-time conversion method, the dynamic radar cross section (DRCS) of the vehicle in horizontal flying path and different flying states is solved in the range of L-band to X-band. The relevant effective statistical features are extracted from the time-scale matrix of DRCS based on continuous wavelet transform (CWT). The numerical results show that: At different flying altitudes, the backscattering radar cross section (BRCS) of the hypersonic vehicle covered with plasma sheath has irregular fluctuation and large range reduction in the head region of the vehicle. Generally, the higher the incident EM wave frequency is, the more detailed the fluctuation characteristics of the DRCS can be reflected. Within the same radar illuminate angle range, the time consumed by the vehicle DRCS is inversely proportional to the flying height and flying speed and is directly proportional to the detection altitude of the SBR. Besides, the different kinds of statistical eigenvalues extracted from DRCS will play an important role in hypersonic targets classification and recognition.


Introduction
Detection and recognition of hypersonic vehicle by radar remote sensing involves the study of dynamic EM scattering of metal targets combined with the plasma flow field [1,2]; especially without considering the influence of near-space clutter and noise, how to obtain reliable DRCS data of hypersonic target becomes more important. We can get different DRCS data according to different flying states of the target spacecraft and establish a complete classification and recognition database. However, when the aircraft is flying at ultra-high speed in near space, the gas molecules around the vehicle will be excited to ionize, and the outer surface of the vehicle will be covered with a thick layer of plasma, which is called "plasma sheath" in the field of radar remote sensing [3,4]. The plasma sheath absorbs and attenuates the detection radar wave emitted by the SBR, and even makes the EM wave deviate from the original propagation track, resulting in irregular and disordered changes of the scattered signals received by the radar, and thus the main information of the target is lost, which brings great challenges to radar remote sensing, tracking and identifying hypersonic targets [5][6][7][8].
In addition, the detection frequency range of the active space observation radar is mainly concentrated in the L-X band, and the main detection method is to emit multi-frequency continuous plane EM waves. Most EM scattering technology research does not consider the use of pulse compression [9]. Besides, compared with small targets such as hypersonic vehicle, most existing space observation radars are low-resolution radars. Therefore, it is of great practical significance to explore and study the dynamic EM scattering and feature extraction of hypersonic vehicles in the low-resolution radar system [10,11].
At present, the main research targets in hypersonic recognition field can be divided into RAM-C and HTV-2 vehicles, both of which are electrically large targets. NASA has carried out experimental research on RAM-C blunt cone reentry vehicle and published a large number of reliable data with extremely high value [12,13], but the cost is very expensive. Nowadays, with the development of computer technology and some reliable EM numerical methods, the EM simulation technology has gradually come into people's vision. In more cases, the real experiments on site are gradually replaced by computer simulation technology. In recent years, many research reports have studied the propagation and scattering of EM waves at different frequency bands in the plasma sheath by using the finite-difference time-domain (FDTD) method [14][15][16][17], the discontinuous Galerkin time-domain method (DGTD) [18], the method of moments (MOM) based on the volume-surface integral equations (VSIE) [19,20], the physical optics (PO) method [21][22][23] based on the inhomogeneity zonal medium model (IZMM) [24] and the Monte Carlo method [25]. On the other hand, some research work is devoted to the measurement and calculation of the DRCS of the target model and the detection and recognition of simple unmanned aerial vehicles (UAVs), fighters and satellites with different physical shape characteristics [26][27][28][29][30]. However, to the authors' knowledge, there are few reports on the dynamic scattering characteristics of hypersonic vehicles flying in near space based on the spaceborne detection platform, and the interaction between EM wave and plasma sheath is not involved.
In general, an important part of hypersonic targets recognition for SBR is to extract the effective statistical features of the DRCS; there are mainly two methods: (1) extract mathematical features directly from time domain, such as extracting cyclic distribution features of DRCS time series or focusing on the establishment and solution of probability density function (PDF) and cumulative distribution function (CDF) [31,32]; and (2) extract features from transform domain (Fourier transform, wavelet transform, Merlin transform, etc.) [33,34]. Since wavelet transform does not introduce the "cross term" in the time-frequency plane as with the Wigner-Ville transform, and due to the introduction of the concept of multi-scale, the wavelet transform has different resolution at different positions in the time-scale (frequency) plane, which makes the wavelet transform have good local properties in both time domain and frequency domain [35].
In this paper, the parallel PO method is used to establish a more reliable theoretical model to simulate and analyze the dynamic EM scattering characteristics of hypersonic vehicles. To highlight the difference of hypersonic targets in different flying states and reduce the feature dimension, CWT is first needed to extract the features from the DRCS time series of the vehicle, and then the effective statistical features are extracted. The focus of the paper is to provide a database for accurate comparison and reference for SBR remote sensing, detection and recognition of hypersonic targets. The main contributions are as follows: 1.
Based on the flow field data of a hypersonic vehicle covered with plasma sheath, the EM scattering characteristics of target spacecraft at different flying altitudes are calculated and analyzed. On this basis, the horizontal flying path is established; the angle-time conversion method is used to convert the BRCS of the target spacecraft into the DRCS; and the relevant influence of target spacecraft at different flying altitudes, flying speeds and the detected altitude of the space-borne radar on the DRCS is studied. Finally, a hypersonic target tracking and identification database based on space-borne radar in near space is established. 2.
Using the continuous wavelet transform (CWT), the features of DRCS of hypersonic vehicle with different flying states are extracted in the time-scale matrix, and the 10 effective and stable statistical features extracted are taken as the identification features. This conclusion lays a data foundation for space-borne radar to recognize hypersonic targets.
The paper is arranged as follows. In Section 2, the PO method based on the multilayer lossy medium model is briefly introduced, and the angle-time conversion method is used to convert the Remote Sens. 2020, 12, 2740 3 of 20 BRCS of hypersonic vehicles in different flying states into the DRCS. Section 3 calculates and analyzes the BRCS of the target spacecraft at different flying altitudes. On this basis, the influence of the flying altitude, flying speed and space-borne radar detection altitude on DRCS is discussed, and the DRCS database for hypersonic vehicles is established. In Section 4, wavelet transform is carried out for the DRCS observation sequence of the target spacecraft, and effective statistical eigenvalues are extracted from the time-scale matrix. Section 5 summarizes the paper and draws some useful conclusions.

Building the Hypersonic Target Model
In this paper, the hypersonic target is modeled by referring to NASA RAM C-II vehicle with rotationally symmetric structure [13]. Figure 1 shows the hypersonic vehicle physical model. The radius of the ball head is R = 0.1524 m, the total length of the body is L = 0.1524 + 1.295 = 1.4474 m and the half cone angle is α = 9 • . The incident wave electric field is polarized along the X-axis (namely, VV polarization: electric field direction parallel to the XOZ plane (page plane) and perpendicular to the incident direction; all the calculation results of EM scattering in the following sections adopt this polarization mode). The plane EM wave incident along the direction of the head of the blunt cone is the same as the positive direction of the Z-axis. In this paper, based on the electrically large scale model of hypersonic vehicle and considering the main frequency range of spaceborne radar detection (L-band to X-band), the PO method based on the inhomogeneous zonal medium model (IZMM) is used to analyze the dynamic EM scattering of hypersonic vehicle [24]. In particular, the plasma fluid field data of the hypersonic vehicle is taken from the hypersonic vehicle plasma sheath and EM characteristics data manual [36].

Physical Optics (PO) Method Based on the Multilayer Lossy Medium Model (MLMM)
In this study, the BRCS of a hypersonic vehicle is calculated by using the PO method based on the Stratton-Chu scattering field integral equation. Here, the dielectric coefficient of the plasma sheath can be expressed by the Drude model [37] and we solve the collision frequency based on the piecewise calculation method proposed in [38]. According to Stratton-Chu equation and the geometric structure diagram of the hypersonic vehicle in Figure 1, after far-field approximation processing in the case Remote Sens. 2020, 12, 2740 4 of 20 of an incident plane EM wave, the electric and magnetic scattering fields of external point P can be calculated by Equations (1) and (2) [21]: where S l represents the bright area of the scatterer; E T and H T are the total electric field and the total magnetic field on the boundary, respectively; k and η are the propagation constant in free space and the intrinsic impedance, respectively;R is the unit vector of the scattered wave; r is the unit vector of the position of the scatterer unit; r is the position vector of point P (as shown in Figure 1); and r = → r and r = → r .n is the normal unit vector of the surface element at r . In addition, θ i is the incident angle of the plane EM wave,ê ⊥ vector is the polarization direction of the incident and reflected electric fields perpendicular to the incident plane and E ⊥ and E // are the field components of incident electric field inê ⊥ andê i // directions, respectively. R ⊥ and R // are the reflection coefficients of the dielectric surface at vertical and parallel polarization, respectively. Here, the Stepping method is adopted to solve it [22]. In addition, the incident EM wave is usually approximated by the plane EM wave: According to Equations (1)- (6), under far-field scattering, the computational equation of the BRCS of a hypersonic vehicle is as follows:

Angle-Time Conversion Method Based on Horizontal Flying Path
Dynamic RCS (DRCS) can show not only the scattering ability of the target spacecraft to the EM wave radiated by radar, but also the flying altitude, flying speed, radar line-of-sight angle and other feature information of the target spacecraft. In this paper, the hypersonic vehicle is set as the maneuver mode of side station level flying (i.e., the lift of the vehicle is equal to gravity and the thrust is equal to the drag); in the process of side station level flying, the azimuth angle of radar line of sight to the target remains unchanged (default the included angle is 0 • in calculation), the roll angle is 0 • and only the pitch angle θ is changed. This section describes the use of a quasi-static method to obtain the dynamic RCS sequence of the target spacecraft, namely the hypersonic vehicle RAM C-II as the research object.
In addition, the target spacecraft is researched along a horizontal air route path when flying at different altitudes and speeds. The dynamic RCS characteristics along with the change of time and frequency response characteristics of radar detection and identification of the hypersonic vehicle are the basis of the simulation. Here, the angle-time conversion method is used to solve the dynamic RCS sequence of the vehicle. The flying time of the hypersonic target varies unevenly with the unit angle during the distance from the space-borne radar. As shown in Figure 2, H = 200, 000 m is the height from the ground of the space-borne radar mounted on the non-synchronous satellite,H 1 and H 2 are the different flying heights of the hypersonic vehicle and H 0 is the shortest vertical distance between the vehicle and the space-borne radar. that is, when the angle between the incident direction of the space-borne radar EM wave and the axis of the target spacecraft is 2.5 • , the initial time at this time is recorded as 0 s. Here, the first half of the flying path of the vehicle is taken as the calculation interval, the initial value is θ = 89.5 • , the change in unit angle is ∆θ =1 • and the range of change times of the angle in the first half is 0 ≤ n ≤ 87. Then, the interval in which the angle changes in the first half flying of the hypersonic vehicle is: Further, if v is the average speed of the vehicle on the horizontal flying path, then Remote Sens. 2020, 12, 2740 6 of 20 For the first half of the flying path, the cumulative time of the vehicle from infinite distance to each fixed angle is: T n max +1 = n max n=0 t 88−n (0 ≤ n max ≤ 87) (13) Similarly, in the second half of the flying path, the time taken for the vehicle to fly to each fixed angle is symmetrical to that of the first half. Then, the vehicle flies from the vertical shortest distance point (θ = 90 • ) of the space-borne radar to the infinite distance. The accumulated time at each fixed angle is: where m is the number of times the angle changes in the second half, and the accumulated time in the first and second half is the total time taken by the hypersonic vehicle to follow the horizontal flying path at a fixed altitude. Therefore, the accumulated time taken by the vehicle to reach each radar incident angle on the horizontal flying path can be calculated. In addition, when the hypersonic vehicle flies to a fixed position along the flying path, in other words, when it reaches a certain angle θ between the flying direction and the incident direction of space-borne radar incident, it takes different times for an vehicle at different altitudes to reach the same position at the same speed v. As shown in Figure 2, when the flying altitude of the vehicle is H 2 < H 1 , in both cases, the flying direction of the vehicle and the detection direction of space-borne radar are at the same included angle θ, the horizontal distance x 1 > x 1 , namely, t 1 > t 1 . In other words, when the hypersonic vehicle is flying at different altitudes, the accumulated flying time is different when the target spacecraft and the radar irradiation direction form an included angle. The lower is the flying altitude, the longer is the required flying time.
To sum up, this paper takes hypersonic vehicle as the research object and simulates the dynamic EM scattering characteristics of its flying on a set path. The specific steps are as follows: 1.
Obtain the plasma flow field data of hypersonic vehicles in different flying states and establish a BRCS EM scattering research model for hypersonic vehicles.

2.
The physical optics (PO) method based on IZMM is used to simulate the EM scattering characteristics of hypersonic vehicles at different radar line-of-sight angles.

3.
According to the performance parameters of the vehicle and the basic theory of aerodynamics, the flying path of the hypersonic target is set, and the angle-time conversion method is used to convert the BRCS of the vehicle into the DRCS.

Influence of Radar Observation Angle on EM Scattering Properties of Hypersonic Vehicle
This section calculates the BRCS of hypersonic vehicle in the range of 0.5 • to 179.5 • at different flying altitudes. Physical modeling and EM simulation calculations were conducted on a PC with 8 GB RAM and eight 2.6 GHz Intel i7-6700 HQ CPUs. The number of surface subsections of the outer surface of the target aircraft varied according to the flying states of the vehicle covered with the plasma sheath. Table 1 lists the number of triangular patches that subdivide the geometric model (Taking 1/8 of 12 GHz EM wavelength as the dividing criterion) and the time required to calculate the BRCS of the hypersonic targets at 20 Ma and different flying altitudes. Figure 3 shows the electron density distribution of hypersonic vehicle at different flying altitudes [36]. It can be found that, with the increase of flying altitude of the hypersonic target, the plasma sheath thickness around the vehicle becomes thinner gradually, and the electron density begins to decline. This is mainly because, as the air at high altitude is relatively thin, the air ionization is less, and the electron density of plasma coated around the vehicle surface is low. When the flying altitude drops from 70 to 30 km, the low altitude atmosphere is relatively dense, the air on the vehicle surface has large ionization and the electron density of plasma increases gradually.  Figure 3 shows the electron density distribution of hypersonic vehicle at different flying altitudes [36]. It can be found that, with the increase of flying altitude of the hypersonic target, the plasma sheath thickness around the vehicle becomes thinner gradually, and the electron density begins to decline. This is mainly because, as the air at high altitude is relatively thin, the air ionization is less, and the electron density of plasma coated around the vehicle surface is low. When the flying altitude drops from 70 to 30 km, the low altitude atmosphere is relatively dense, the air on the vehicle surface has large ionization and the electron density of plasma increases gradually.  Compared with the curve fluctuation of pure metal blunt cone BRCS, the existence of plasma sheath makes the BRCS of the hypersonic target have a large irregular fluctuation and attenuation in most of the incident angle directions, especially in the illuminate angle range of the vehicle head area. This is due to the high electron density of plasma in the head area of the vehicle; however, the sheath is very thin (as the plasma flow field distribution in Figure 3) and the plasma will produce large collision and absorption of the incident EM wave. In addition, as the flying altitude of the hypersonic target continues to decrease, the BRCS gradually decreases at all angles, especially at the incident angle of 81°. The peak value of the BRCS at the flying altitude of the vehicle at 70 km from the ground is larger than at the other four flying altitudes. Similarly, at an altitude of 30 km, the BRCS of the  Figure 4 shows the BRCS comparison of hypersonic vehicles covered with plasma sheath at different flying altitudes. The EM wave of the detection radar is incident from the head of the blunt cone with an initial angle of θ = 0 • . When the incident angle is θ = 81 • and θ = 261 • , the projection area of the blunt cone spacecraft along the direction of radar irradiation is the largest, and the BRCS reaches a maximum peak point at this point. Similarly, BRCS reaches a maximum value at the tail of the blunt cone. Compared with the curve fluctuation of pure metal blunt cone BRCS, the existence of plasma sheath makes the BRCS of the hypersonic target have a large irregular fluctuation and attenuation in most of the incident angle directions, especially in the illuminate angle range of the vehicle head area. This is due to the high electron density of plasma in the head area of the vehicle; however, the sheath is very thin (as the plasma flow field distribution in Figure 3) and the plasma will produce large collision and absorption of the incident EM wave. In addition, as the flying altitude of the hypersonic target continues to decrease, the BRCS gradually decreases at all angles, especially at the incident angle of 81 • . The peak value of the BRCS at the flying altitude of the vehicle at 70 km from the ground is larger than at the other four flying altitudes. Similarly, at an altitude of 30 km, the BRCS of the hypersonic vehicle covered with a plasma sheath is reduced in almost all incident directions (see the red line in Figure 4). hypersonic vehicle covered with a plasma sheath is reduced in almost all incident directions (see the red line in Figure 4).

Dynamic EM Scattering Analysis of Hypersonic Vehicle
In this section, the space-borne radar is located on a non-synchronous satellite at a height of 200-300 km, and the EM wave bands irradiated by the detection radar are the L-X bands. We studied the effects of incident EM wave frequency, space-borne radar detection altitude, flying altitude and flying speed of the target spacecraft on the dynamic RCS time series of a hypersonic vehicle covered with a plasma sheath when the range of the angle between the target spacecraft and the radar irradiation direction is  oo 2.5 θ 177.5 . Our objective was to establish a relatively complete dynamic RCS target recognition database for hypersonic vehicle of low-resolution radar. The simulation results are shown in Figure 5a-d.

Dynamic EM Scattering Analysis of Hypersonic Vehicle
In this section, the space-borne radar is located on a non-synchronous satellite at a height of 200-300 km, and the EM wave bands irradiated by the detection radar are the L-X bands. We studied the effects of incident EM wave frequency, space-borne radar detection altitude, flying altitude and flying speed of the target spacecraft on the dynamic RCS time series of a hypersonic vehicle covered with a plasma sheath when the range of the angle between the target spacecraft and the radar irradiation direction is 2.5 • ≤ θ ≤ 177.5 • . Our objective was to establish a relatively complete dynamic RCS target recognition database for hypersonic vehicle of low-resolution radar. The simulation results are shown in Figure 5a-d. hypersonic vehicle covered with a plasma sheath is reduced in almost all incident directions (see the red line in Figure 4).

Dynamic EM Scattering Analysis of Hypersonic Vehicle
In this section, the space-borne radar is located on a non-synchronous satellite at a height of 200-300 km, and the EM wave bands irradiated by the detection radar are the L-X bands. We studied the effects of incident EM wave frequency, space-borne radar detection altitude, flying altitude and flying speed of the target spacecraft on the dynamic RCS time series of a hypersonic vehicle covered with a plasma sheath when the range of the angle between the target spacecraft and the radar irradiation direction is  oo 2.5 θ 177.5 . Our objective was to establish a relatively complete dynamic RCS target recognition database for hypersonic vehicle of low-resolution radar. The simulation results are shown in Figure 5a-d.  The technical approach we follow is to convert the backscattering RCS of the target spacecraft with the change in radar incident angle into the dynamic RCS time series as a function of time. The dynamic RCS curve in Figure 5a clearly shows that, with the passage of time, for a target spacecraft flying at a constant speed along the horizontal path, when the distance between the space-borne radar and the target spacecraft changes from far to near and then to far, the fluctuations of the dynamic RCS curve change from gentle to violent and then to gentle. In addition, the distribution of scattering data points of the hypersonic target is from sparse to dense and then to sparse (as shown in Figure 5a, where the curve changes in the three dotted boxes).
When the hypersonic vehicle is located at a great distance from the space-borne radar, its dynamic RCS curve fluctuates gently, and the scattered data points on the curve are distributed sparsely. However, as the target spacecraft approaches the detection radar, until the vehicle flies directly below the space-borne radar (the line between the two centers is perpendicular to the horizontal flying path, i.e., o θ = 90 ), the scattering points on the dynamic RCS curve are the densest with time, and the curve oscillations and jitters are the most intense. This is because the target vehicle is traveling at a constant speed, and, in the process of reducing the distance from the space-borne radar from far to near, the flying time varies unevenly with the unit angle. In other words, as the distance between the vehicle and the detection radar gradually shortens, the flying time of the vehicle gradually decreases when the target spacecraft is scanned by the detection radar at a unit angle. Similarly, while the hypersonic vehicle is flying away from the satellite-borne radar at constant speed, when the detection radar illuminates the target spacecraft and sweeps through the unit angle, the flying time of the vehicle gradually increases from short to long, and the fluctuation of the dynamic RCS curve gradually becomes gentle.
The numerical distribution and fluctuation change of the dynamic RCS curve are related to the flying altitude and flying speed. According to the flying path and the relative position of the space-borne radar and hypersonic vehicle, it can be concluded that, as shown in Figure 5a, the higher is the flying altitude The technical approach we follow is to convert the backscattering RCS of the target spacecraft with the change in radar incident angle into the dynamic RCS time series as a function of time. The dynamic RCS curve in Figure 5a clearly shows that, with the passage of time, for a target spacecraft flying at a constant speed along the horizontal path, when the distance between the space-borne radar and the target spacecraft changes from far to near and then to far, the fluctuations of the dynamic RCS curve change from gentle to violent and then to gentle. In addition, the distribution of scattering data points of the hypersonic target is from sparse to dense and then to sparse (as shown in Figure 5a, where the curve changes in the three dotted boxes).
When the hypersonic vehicle is located at a great distance from the space-borne radar, its dynamic RCS curve fluctuates gently, and the scattered data points on the curve are distributed sparsely. However, as the target spacecraft approaches the detection radar, until the vehicle flies directly below the space-borne radar (the line between the two centers is perpendicular to the horizontal flying path, i.e., θ = 90 • ), the scattering points on the dynamic RCS curve are the densest with time, and the curve oscillations and jitters are the most intense. This is because the target vehicle is traveling at a constant speed, and, in the process of reducing the distance from the space-borne radar from far to near, the flying time varies unevenly with the unit angle. In other words, as the distance between the vehicle and the detection radar gradually shortens, the flying time of the vehicle gradually decreases when the target spacecraft is scanned by the detection radar at a unit angle. Similarly, while the hypersonic vehicle is flying away from the satellite-borne radar at constant speed, when the detection radar illuminates the target spacecraft and sweeps through the unit angle, the flying time of the vehicle gradually increases from short to long, and the fluctuation of the dynamic RCS curve gradually becomes gentle. The numerical distribution and fluctuation change of the dynamic RCS curve are related to the flying altitude and flying speed. According to the flying path and the relative position of the space-borne radar and hypersonic vehicle, it can be concluded that, as shown in Figure 5a, the higher is the flying altitude (H x = 30 km, 40 km, 50 km, 60 km and 70 km), that is, the closer is the hypersonic vehicle to the space-borne radar, the shorter is the time required for the radar to sweep through the fuselage unit angle at the same flying Mach number (v = 20 Mach). In other words, the faster does the dynamic RCS fluctuation of the target spacecraft change in the same time range, the shorter is the time for the curve to change from sparse to dense and then to sparse.
Similarly, as shown in Figure 5b, at the same flying altitude (H x = 60 km), the faster is the hypersonic vehicle flying (v = 15 Mach, 20 Mach and 25 Mach), the shorter is the time required for the radar to sweep through the fuselage unit angle, that is, the target dynamic RCS curve can oscillate to the highest peak point in a short time. For example, the target spacecraft flew at an altitude of 60 km and a speed of 25 Mach (as shown in Figure 5b, the pink curve) when the dynamic RCS curve of the target spacecraft reached the densest center point. That is, when the connection between the vehicle and the space-borne radar center is perpendicular to the horizontal flying path, the dynamic RCS curve of the blunt cone spacecraft with flying altitude of 60 km and flying speed of 15 Mach (red curve in Figure 5b) is still in the gentle and sparse area.
Next, we studied the change in the dynamic RCS of the target spacecraft when the space-borne radar flies at different altitudes (H = 200 km, 250 km and 300 km). In short, the variation rule of dynamic RCS is similar to changes in the flying altitude of the target spacecraft, which is to change the distance between the target spacecraft and the detection radar. The further is the distance between the target spacecraft and the space-borne radar, the longer it takes for the radar to sweep through the fuselage unit angle at the same flying speed (v = 20 Mach) and the slower is the dynamic RCS fluctuation of the target spacecraft in the same time range, as shown by the blue curve of the target dynamic RCS in Figure 5c.
Then, we studied the extent to which the characteristics of the law of the dynamic RCS sequence of the target spacecraft change with the incident EM wave frequency when the frequency range of the incident EM wave is in the L-X bands and the flying state is 60 km and 20 Mach, as shown in Figure 5d. An increase in the frequency of incident EM waves, which would be more obvious, highlights details of the features of the physical structure of the target spacecraft, compared with the target at the low-frequency 1.25 GHz dynamic RCS sequence. When the incident EM wave frequency increases to 8.25 GHz, the dynamic RCS time series curve of the target fluctuates and oscillates more dramatically, and the peak point of the curve becomes sharper when θ = 81 • . In addition, EM waves at higher frequencies can partially or completely penetrate the plasma sheath covering the surface of the target spacecraft, and the energy attenuation is greater because EM waves have a longer propagation path inside the plasma. However, when the hypersonic vehicle covered with a plasma sheath is illuminated by the EM wave at a lower frequency, the EM wave cannot penetrate the plasma sheath, and most of the EM waves are reflected from the surface of the plasma sheath, resulting in the high value of the dynamic RCS curve of the target spacecraft in a large time range. Therefore, effectively reducing the frequency of EM waves emitted by space-borne radar would be more helpful for detecting and tracking a hypersonic vehicle covered with a plasma sheath. Finally, we calculated the frequency of incident EM wave as f = 8.25 GHz; the incident electric field and scattering electric field are in the same polarizations mode (VV/HH), keep the same flying altitude (60 km) and at the different flying speeds. The DRCS time series of hypersonic vehicle are shown in Figure 6a,b (where the detection altitude of space-borne radar remains unchanged, H = 200 km). Figure 6a,b shows that the polarization mode of radar antenna transmitting and receiving EM wave does have an impact on the DRCS time series of the hypersonic targets, which is mainly reflected in the slight difference of DRCS curve waveform trend and oscillation mode. Although these changes are irregular and fuzzy, they do have certain impact on hypersonic targets recognition. In addition, for the two polarization modes (VV/HH), the peak values of DRCS time series at θ = 81 • and θ = 180

Mathematical Sensitivity Function of DRCS
Due to the existence of plasma sheath and the maneuverability of hypersonic targets in space and time, the DRCS of target spacecrafts in different flying states tend to be quite different; thus, the change of DRCS is sensitive, which makes it possible for us to detect and identify hypersonic vehicles effectively. Based on the mathematical characteristics and distribution form of hypersonic targets DRCS, we use probability density function (PDF) and cumulative distribution function (CDF) to solve the probability distribution form of hypersonic targets DRCS coefficient under different flying altitudes, the same flying speed (20 Ma)

Mathematical Sensitivity Function of DRCS
Due to the existence of plasma sheath and the maneuverability of hypersonic targets in space and time, the DRCS of target spacecrafts in different flying states tend to be quite different; thus, the change of DRCS is sensitive, which makes it possible for us to detect and identify hypersonic vehicles effectively. Based on the mathematical characteristics and distribution form of hypersonic targets DRCS, we use probability density function (PDF) and cumulative distribution function (CDF) to solve the probability distribution form of hypersonic targets DRCS coefficient under different flying altitudes, the same flying speed (20 Ma) (Figure 7a,c). In general, in the same detection period of space-borne radar (SBR), when the flying speed remains constant (20 Ma), the higher is the altitude of hypersonic vehicle, the more obvious is the broadening of PDF curve of DRCS, and the peak height of probability density is lower. This is because the higher is the flying altitude, the shorter is the relative distance between the target spacecraft and the SBR. In the same radar detection time range, the illuminate angle range of the radar scanning the hypersonic target is larger, and the distribution range of DRCS coefficient in the range [−60 dBsm, 30 dBsm] becomes wider; accordingly, the peak value of DRCS coefficient distribution probability is lower. Similarly, at the same flying altitude (60 km), the larger is the flying Mach number, the greater is the broadening degree of PDF curve and the lower is the peak probability density, and the principle are the same.
In addition, different flying states of hypersonic targets affect seriously the slope of the CDF curve and the range of DRCS coefficient covered by CDF in the range of [−60 dBsm, 30 dBsm] (Figure 7b,d). Maintaining a certain flying speed (20 Ma), the higher is the flying altitude, the smaller is the slope of CDF curve and the more gentle is its distribution. Similarly, at the same flying altitude, the larger is the flying Mach number, the wider is the DRCS coefficient range covered by the CDF curve, and the trend of curve change will be more gentle.   (Figure 7a,c). In general, in the same detection period of space-borne radar (SBR), when the flying speed remains constant (20 Ma), the higher is the altitude of hypersonic vehicle, the more obvious is the broadening of PDF curve of DRCS, and the peak height of probability density is lower. This is because the higher is the flying altitude, the shorter is the relative distance between the target spacecraft and the SBR. In the same radar detection time range, the illuminate angle range of the radar scanning the hypersonic target is larger, and the distribution range of DRCS coefficient in the range  [ 60 dBsm, 30 dBsm] becomes wider; accordingly, the peak value of DRCS coefficient distribution probability is lower. Similarly, at the same flying altitude (60 km), the larger is the flying Mach number, the greater is the broadening degree of PDF curve and the lower is the peak probability density, and the principle are the same. In addition, different flying states of hypersonic targets affect seriously the slope of the CDF curve and the range of DRCS coefficient covered by CDF in the range of  [ 60 dBsm, 30 dBsm]

Continuous Wavelet Transform and Time-Scale Plane
In this subsection, the dynamic RCS time series of the target spacecraft is extracted by using continuous wavelet transform (CWT) of the joint time-scale plane description signal, where the scale corresponds to the frequency. Based on the different frequency information, the processing of the dynamic RCS observation sequence of low-resolution space-borne radar can clearly show the details of the RCS observation sequence in the time domain and frequency domain, which would be more conducive to extracting the features of the hypersonic vehicle for different flying states.
Scale factor a and time factor b are discretized to realize CWT. The observation sequence of dynamic RCS of space-borne radar is set as DRCS(t), t = 1, 2, . . . , T b , where T b is the total length of dynamic RCS observation data. Then, the wavelet transform of DRCS(t) is: where b = 1, 2, . . . , T b and a = 1, 2, . . . , S a . To enable the final results to be compared for different flying states of the target spacecraft, during the flying process, the range of angle variation with which the space-borne radar illuminates the target spacecraft was always 2.5 • ≤ θ ≤ 177.5 • ; thus, the total time T b in DRCS(t) of the target spacecraft for different flying states was different. In the simulations, the time interval was T s = 1 s, the maximum value of scale factor was S a = 600 and the Symlet(symN) wavelet (compactly supported wavelets with least asymmetry and highest number of vanishing moments for a given support width) was used as the mother wavelet, where N was calculated as 4. As shown in Figure 8a-d, when the flying altitude of spaceborne radar is H = 200 km, for four different flying states and incident EM wave frequencies, the time-scale diagram of dynamic RCS time series of the hypersonic vehicle is obtained after continuous wavelet transform.
To enable the final results to be compared for different flying states of the target spacecraft, during the flying process, the range of angle variation with which the space-borne radar illuminates the target spacecraft was always  oo 2.5 θ 177.5 ; thus, the total time b T in DRCS(t) of the target spacecraft for different flying states was different. In the simulations, the time interval was s T = 1 s , the maximum value of scale factor was a S = 600 and the Symlet(symN) wavelet (compactly supported wavelets with least asymmetry and highest number of vanishing moments for a given support width) was used as the mother wavelet, where N was calculated as 4. As shown in Figure 8a- Figure 8a shows the time-scale diagram of the target spacecraft at 60 km and 20 Mach. The abscissa is the time factor, indicating the total flying time of the target spacecraft, and the ordinate is the scale factor, representing the scale range of the wavelet transform. Here, the scale factor corresponds to the frequency, i.e., the smaller is the scale factor, the greater is the frequency, and vice versa. In other words, when the scale factor a is smaller, more wavelet functions are needed to fit the dynamic RCS time series in this time period, which can more accurately reflect the fluctuation details of the dynamic RCS sequence of the target spacecraft over time. The color intensity in Figure 8a represents the size of the absolute value of the wavelet coefficient, i.e., the brighter is the color, the greater is the value of the wavelet coefficient. For the dynamic RCS time series of the hypersonic vehicle, when the curve fluctuates violently and jitters and sudden changes occur, the inner product of the dynamic RCS function and the sym4 wavelet function shows that the result value is the largest; that is, the value of the wavelet coefficient is the largest. In other words, the higher is the similarity between the target dynamic RCS time series and the sym4 wavelet of the same length in  Here, the scale factor corresponds to the frequency, i.e., the smaller is the scale factor, the greater is the frequency, and vice versa. In other words, when the scale factor a is smaller, more wavelet functions are needed to fit the dynamic RCS time series in this time period, which can more accurately reflect the fluctuation details of the dynamic RCS sequence of the target spacecraft over time. The color intensity in Figure 8a represents the size of the absolute value of the wavelet coefficient, i.e., the brighter is the color, the greater is the value of the wavelet coefficient. For the dynamic RCS time series of the hypersonic vehicle, when the curve fluctuates violently and jitters and sudden changes occur, the inner product of the dynamic RCS function and the sym4 wavelet function shows that the result value is the largest; that is, the value of the wavelet coefficient is the largest. In other words, the higher is the similarity between the target dynamic RCS time series and the sym4 wavelet of the same length in this period, the brighter is the representative part in the time-scale diagram.
Because the fluctuation change recorded by the wavelet function is 0 by default at the initial position, the initial value of the target dynamic RCS curve is less than 0. For example, in Figure 8b, the value of the dynamic RCS curve of the target spacecraft at the initial position is −32 dBm; hence, the wavelet transform determines by default that the initial position of the dynamic RCS sequence of the target suddenly changes. In other words, it is shown in the time-scale diagram that relatively bright longitudinal stripes will appear in the initial time period, as shown in the area enclosed by a red dashed line in Figure 8b.
In addition, when the target spacecraft is flying in the area covered by the space-borne radar during the period of 2000-3000 s, the target dynamic RCS curve fluctuates violently, and the corresponding time-scale diagram shows bright longitudinal stripes, as shown in the area enclosed by the yellow dashed line in Figure 8b. During the period of 3000-5000 s, the dynamic RCS curve fluctuates gently. At the small scale, the wavelet coefficient is almost 0, and the time-scale image is darker. Compared with the larger scale range, the dynamic RCS curve at this time exhibits weak fluctuation, the wavelet coefficient value is small and the time-scale image has low brightness, as shown in the area enclosed by the green dashed line in Figure 8b. In Figure 8c,d, the dynamic RCS curves of the target spacecraft under these two calculation conditions fluctuate wildly. Therefore, the time-scale image after wavelet transformation has relatively dense longitudinal bright stripes, which can clearly distinguish the characteristics of the target dynamic RCS curve for different flying states.

Feature Extraction of Dynamic RCS on Time-Scale Plane
Here, we set A = W DRCS (a, b) . Because the volume of data of Matrix A after wavelet transformation is very large with extensive information redundancy, it is obviously unrealistic to take it directly as a feature, which means we need to extract features with fewer stable dimensions from Matrix A. Because the dynamic RCS observation sequence of the target has random stability, it can be considered to extract mathematical statistical features from Matrix A. After many simulation experiments and using a method based on [39,40], we extracted 10 valid statistical features. They are the maximum singular value feature, mean value feature, maximum value feature, variance feature, scale barycenter feature, central moment feature and invariant moment feature, respectively. The specific calculation formulas appear in Appendix A. We can combine these 10 effective statistical features in a template for target recognition, which is expected to lay a solid foundation for the recognition of hypersonic vehicle in the future. Table 2 lists the 10 relatively stable statistical features extracted from the time-scale plane for different flying states and calculation conditions.
In addition, we examined whether the influence of the plasma sheath on the eigenvalues of the dynamic RCS sequence of a hypersonic vehicle could be more convincingly explained. Suppose that along the horizontal flying path of the target spacecraft, a pure metal blunt cone spacecraft model is placed at each unit angle incident position of space-borne radar. The radar illumination conditions are H = 200 km and f = 3.25 GHz. When the target spacecraft is at 60 km and 20 Ma, the change law of time as a function of angle is taken as the calculation condition to solve the dynamic RCS of the metallic blunt cone spacecraft. The ten effective eigenvalues mentioned above were extracted on the scale-time plane and compared with those extracted by the hypersonic vehicle covered with a plasma sheath at a flying altitude of 60 km and a speed of 20 Ma under the same radar incident conditions, as listed in Table 3.
The eigenvalue data shown in Tables 2 and 3 indicate that, for features f 1 , f 3 and f 4 , which differ from the other features, the eigenvalue of pure metal targets is larger than that of the other targets covered with a plasma sheath, which plays an important role in the classification and recognition of target spacecraft. At the same time, it also shows that the contributions of various original eigenvalues extracted from the dynamic RCS samples of the target spacecraft for classification are different, and the degree of importance of different eigenvalues to samples also varies. In addition, for hypersonic vehicles, the existence of a plasma sheath can greatly change the size of the ten effective statistical eigenvalues, by differentiating and changing them depending on the radar illumination conditions and flying states of the target spacecraft. These feature sets are expected to play an important role in the classification and recognition of hypersonic vehicles in near space.

Conclusions
In this paper, the parallel PO method based on the IZMM is used to simulate and analyze the dynamic EM scattering characteristics of the hypersonic vehicle covered with plasma sheath under the condition of controllable computing time. By using the angle-time conversion method in the side station level flying maneuver mode, we analyze the DRCS of a hypersonic target for different flying states and under different calculation conditions and establish a relatively complete simulation database. Then, the CWT is used to extract ten valid statistical characteristic values on the time-scale plane, which lays a data foundation for the target recognition of a hypersonic vehicle. The results of our calculations enable us to draw the following conclusions: 1.
According to the EM scattering characteristics of hypersonic vehicle at different flying altitudes, we can find that the plasma sheath can greatly change the fluctuation trend of the BRCS of the hypersonic target, making its oscillation law irregular and difficult to predict. In the large illuminate angle range, the BRCS has a reduction effect, which is due to the collision absorption of the incident EM wave in the plasma environment with high electron density. In addition, when the hypersonic target is flying at different altitudes, the EM parameters of the plasma sheath coated on the surface of the hypersonic vehicle are also different, resulting in the BRCS of the vehicle extremely sensitive to the change of flying state, which is a great challenge for hypersonic target identification.

2.
Based on the flying path along the horizontal trajectory, we calculated the DRCS of the hypersonic vehicle for different flying states, and the influence of flying altitude and velocity on the recognition of hypersonic target using DRCS is quantitatively analyzed by using mathematical sensitivity functions PDF and CDF. The results show that the frequency of the incident EM wave, detection altitude and polarizations of SBR, the flying altitude of the target spacecrafts and the flying speed all affect the DRCS time series of a hypersonic vehicle covered with a plasma sheath, but it is possible to detect and identify hypersonic vehicles. In addition, CWT is used to extract the characteristics of the DRCS sequence of the vehicle, and a relatively complete DRCS database for hypersonic vehicle recognition is established.

3.
Finally, based on the urgent needs and existing bottlenecks of hypersonic targets recognition capability, the research results of this paper explore and reveal the dynamic EM scattering characteristics of plasma sheath flow field of hypersonic vehicle. By showing the mathematical and transformation characteristics of DRCS time series, the EM scattering characteristics of hypersonic vehicle covered with plasma sheath are effectively connected with the field of target recognition, and, as a key link, the technical methods are effectively applied to the actual research. The research work published in this paper provides theoretical analysis basis and sample identification database for hypersonic vehicle detection technology based on satellite platform and provides reliable data support and sample reference for target detection, imaging and recognition in hypersonic field in the future.