Classification Method of Uniform Circular Array Radar Ground Clutter Data Based on Chaotic Genetic Algorithm

The classification and recognition of radar clutter is helpful to improve the efficiency of radar signal processing and target detection. In order to realize the effective classification of uniform circular array (UCA) radar clutter data, a classification method of ground clutter data based on the chaotic genetic algorithm is proposed. In this paper, the characteristics of UCA radar ground clutter data are studied, and then the statistical characteristic factors of correlation, non-stationery and range-Doppler maps are extracted, which can be used to classify ground clutter data. Based on the clustering analysis, results of characteristic factors of radar clutter data under different wave-controlled modes in multiple scenarios, we can see: in radar clutter clustering of different scenes, the chaotic genetic algorithm can save 34.61% of clustering time and improve the classification accuracy by 42.82% compared with the standard genetic algorithm. In radar clutter clustering of different wave-controlled modes, the timeliness and accuracy of the chaotic genetic algorithm are improved by 42.69% and 20.79%, respectively, compared to standard genetic algorithm clustering. The clustering experiment results show that the chaotic genetic algorithm can effectively classify UCA radar’s ground clutter data.


Introduction
Uniform circular array (UVA) radar is a uniform circular array radar composed of isotropic antennas, which has the advantages of space omnidirectional scanning, various search and tracking methods, and flexible beam pointing [1][2][3]. It is especially suitable for the monitoring of low and slow targets under complex terrain background. The signal received by UCA ground-based radar includes the electromagnetic scattering signal generated by the surrounding environment besides the target signal. Ground clutter refers to the echo signal formed after the electromagnetic wave signal emitted by radar is reflected by ground object background [4], which is usually affected by radar system parameters, ground object background type, wind speed and other factors. UCA ground-based radar is easily affected by strong scattering clutter signals during search and tracking missions, which is a crucial factor restricting the performance of UCA radar.
The performance of radar target detection depends on recognizing clutter characteristics and the estimation of relevant statistical parameters to a certain extent. Research on radar clutter signal characteristics can be summarized into four aspects: power spectrum characteristics, non-stationarity and nonlinearity, amplitude distribution characteristics, and correlation. The spectral characteristics of clutter are studied based on measured and low local search efficiency. Introducing chaos theory into the standard PSO algorithm can effectively solve the problem of particle swarm easily falling into the local extreme point [36,37]. Therefore, in the field of unsupervised classification, GA clustering can be used to achieve data classification in different fields.
Based on the classification problem of clutter data collected by UCA ground-based radar, a data classification method based on the Chaotic Genetic Algorithm (CGA) is proposed in this paper. According to the ground clutter data collected in five different environments (microwave anechoic chamber, highway, dry grassland, town, and gravel land), three characteristics of the ground clutter data, including correlation, non-stationarity, and clutter spectrum (range-Doppler domain), were extracted, and the clutter data were classified by the analysis results of characteristic information. In addition, the optimization genetic algorithm of the chaos property is adopted to make the individuals move alternately between chaos and stability and the ergodic property of chaos theory is fully utilized to overcome the problem of precocious convergence of the population and improve the global convergence ability of the population.
In general, the main contributions of this paper are listed below: • Firstly, the characteristics of ground clutter data measured in different UCA groundbased radar scenarios are studied, and the correlation, non-stationary, and statistical characteristics of the range-Doppler domain of clutter data are analyzed. • Secondly, a GA clustering method based on chaos theory is proposed to overcome standard GA's defects, such as premature convergence and weak local optimization ability, and complete the data classification and recognition according to the feature factors extracted from the measured clutter data.
The rest of this paper is organized as follows: In Section 2, the basic structure of the UCA radar system, the process of data preprocessing, the usual analysis method of ground clutter data (correlation analysis, non-stationary analysis, and statistical analysis of range-Doppler map), and the basic principle of clustering algorithm adopted in this paper are introduced. In Section 3, the classification results of ground clutter data in different scenarios and two kinds of beam control modes are presented and evaluate the two clustering algorithms' classification performance. Finally, the conclusions and prospects of this study are summarized in Section 4.

Uniform Circular Array Radar and Experiment Sites
A uniform circular array radar is used for data collection. The radar system consists of a receiver, transmitter, transceiver antenna, signal processing unit, and power module. The antenna part comprises eight waveguide slot antennas with the same array element structure and adopts the operation mode of split transceiver. The UCA radar transceiver antenna is designed by eight conformal array elements with an interval of 45 • along the circumferential azimuth direction, and the azimuth 360 • full coverage circumferential scanning is realized by means of phase weighting. The antenna emission pattern is a multi-main-lobe pattern, and eight antennas receive signals independently. The simulation diagram of azimuth emission in two-beam control modes is presented in Figure 1. When the channels are in phase, the synthetic antenna's emission pattern has eight main lobes, pointing to 0 • , 45 • , 90 • , 135 • , 180 • , 225 • , 270 • , and 315 • respectively (Figure 1a). When each channel is combined at a phase interval of 90 • , the main lobe directions of the antenna emission pattern are 22.5 • , 67.5 • , 112.5 • , 157.5 • , 202.5 • , 247.5 • , 292.5 • , and 337.5 • respectively ( Figure 1b). In addition, the pitching 3 dB beam width of a single element is 17.47 • , and the pitching main beam direction of all transmitting and receiving array element is 90 • . The radar system realizes the whole scene's detection through an alternating transformation of the two beam-controlled modes. The beam control code switches to display three states during regular operation ( Figure 1). The experiment includes microwave anechoic chamber experiment and field experiment ( Figure 2). The microwave anechoic chamber is equipped with uniformly distributed absorbing materials and is equipped with a test stand and probe. Through the automatic control of the rotating platform, the radar system's antenna parameters and the effectiveness of sending and receiving data can be easily tested (Figure 2b). Four scenarios with different environmental complexity were selected for multiple measurement experiments, which were located in Wuxi City, Jiangsu Province, Pinggu District, Beijing, Haidian District, Beijing, and Bengbu City, Anhui Province, China. According to topography, environment, and functional uses, it is classified as highway, dry grassland, town, and sandy land. Trees and shrubs are planted on both sides of the road with green protection belts (Figure 2c). There are a large number of weeds and a few low shrubs in the grassland monitoring area (Figure 2d). In the urban environment, a typical parking lot is selected, and there are several vehicles with strong interference parked in the monitoring area ( Figure 2e). A little gravel and low weeds are scattered on the gravel ground ( Figure  2f). Different background components and materials in the above selected experimental test scenes are the signal sources of ground clutter in the radar data under four different environments. The experiment includes microwave anechoic chamber experiment and field experiment ( Figure 2). The microwave anechoic chamber is equipped with uniformly distributed absorbing materials and is equipped with a test stand and probe. Through the automatic control of the rotating platform, the radar system's antenna parameters and the effectiveness of sending and receiving data can be easily tested (Figure 2b). Four scenarios with different environmental complexity were selected for multiple measurement experiments, which were located in Wuxi City, Jiangsu Province, Pinggu District, Beijing, Haidian District, Beijing, and Bengbu City, Anhui Province, China. According to topography, environment, and functional uses, it is classified as highway, dry grassland, town, and sandy land. Trees and shrubs are planted on both sides of the road with green protection belts ( Figure 2c). There are a large number of weeds and a few low shrubs in the grassland monitoring area (Figure 2d). In the urban environment, a typical parking lot is selected, and there are several vehicles with strong interference parked in the monitoring area (Figure 2e). A little gravel and low weeds are scattered on the gravel ground ( Figure 2f). Different background components and materials in the above selected experimental test scenes are the signal sources of ground clutter in the radar data under four different environments.

Data Pre-Processing
Since each antenna of the radar receives signals independently, a pulse contains e echo data channels, and the data structure is shown in Figure 3. In the distance dimens the received signals are discharged in accordance with the sequence of channels and pulse data of each channel contains frame header information, which is used to mark pulse number and the wave-controlled mode (Figure 3a). Radar parameters such as p repetition frequency (PRF), sampling points of range dimension, and coherent proces interval (CPI) are slightly different among the data collected in different experimental narios. To observe the composition and structure of radar data, the frame header in mation is removed in Figure 3a, and multiple radar pulse data of eight receiving chan are spliced into a data cube (Figure 3b). The longitudinal direction represents the d sampling of the slow time dimension. The UCA radar data cube's structure model is basis for analyzing radar data characteristics and understanding signal processing op tions.

Data Pre-Processing
Since each antenna of the radar receives signals independently, a pulse contains eight echo data channels, and the data structure is shown in Figure 3. In the distance dimension, the received signals are discharged in accordance with the sequence of channels and the pulse data of each channel contains frame header information, which is used to mark the pulse number and the wave-controlled mode (Figure 3a). Radar parameters such as pulse repetition frequency (PRF), sampling points of range dimension, and coherent processing interval (CPI) are slightly different among the data collected in different experimental scenarios. To observe the composition and structure of radar data, the frame header information is removed in Figure 3a, and multiple radar pulse data of eight receiving channels are spliced into a data cube (Figure 3b). The longitudinal direction represents the data sampling of the slow time dimension. The UCA radar data cube's structure model is the basis for analyzing radar data characteristics and understanding signal processing operations.  Before analyzing the characteristics of ground clutter, it is necessary to preprocess the original data collected in different scenes, including three parts ( Figure 4). Data processing steps of UCA radar are listed as follows: (1) Channel data separation: By analyzing the frame header information of eight channels' pulse-echo data, the pulse sequence and beam control modes are marked; (2) Data extraction and reorganization determine the switching state of radar beam pointing, and data collected under the same beam control mode of eight channels are grouped into the same group; and (3) Data preprocessing and analysis, including Hilbert transform, range-direction FFT, denoising, normalization processing, and data statistical analysis.
Since the UCA radar's working mode is the direct sampling of time-domain signals, it is necessary to perform a Hilbert transformation on the sampled signals and obtain two signals for coherent processing.

Characteristic Factors
The correlation, non-stationarity, and range-Doppler characteristics are selected as the characteristic factors of ground clutter and the characteristic differences in five different scenarios and two beam control modes are analyzed.

Correlation Analysis of Ground Clutter Data
Pearson correlation coefficient (PCC) is used to measure the correlation between different data types [38], and the beam-related correlation function, azimuth correlation function, and range-related correlation function are obtained. Beam correlation refers to the correlation between the sampled data obtained by two groups of different beam  Before analyzing the characteristics of ground clutter, it is necessary to preprocess the original data collected in different scenes, including three parts ( Figure 4). Data processing steps of UCA radar are listed as follows: (1) Channel data separation: By analyzing the frame header information of eight channels' pulse-echo data, the pulse sequence and beam control modes are marked; (2) Data extraction and reorganization determine the switching state of radar beam pointing, and data collected under the same beam control mode of eight channels are grouped into the same group; and (3) Data preprocessing and analysis, including Hilbert transform, range-direction FFT, denoising, normalization processing, and data statistical analysis. Before analyzing the characteristics of ground clutter, it is necessary to preprocess the original data collected in different scenes, including three parts ( Figure 4). Data processing steps of UCA radar are listed as follows: (1) Channel data separation: By analyzing the frame header information of eight channels' pulse-echo data, the pulse sequence and beam control modes are marked; (2) Data extraction and reorganization determine the switching state of radar beam pointing, and data collected under the same beam control mode of eight channels are grouped into the same group; and (3) Data preprocessing and analysis, including Hilbert transform, range-direction FFT, denoising, normalization processing, and data statistical analysis.
Since the UCA radar's working mode is the direct sampling of time-domain signals, it is necessary to perform a Hilbert transformation on the sampled signals and obtain two signals for coherent processing.

Characteristic Factors
The correlation, non-stationarity, and range-Doppler characteristics are selected as the characteristic factors of ground clutter and the characteristic differences in five different scenarios and two beam control modes are analyzed.

Correlation Analysis of Ground Clutter Data
Pearson correlation coefficient (PCC) is used to measure the correlation between different data types [38], and the beam-related correlation function, azimuth correlation function, and range-related correlation function are obtained. Beam correlation refers to the correlation between the sampled data obtained by two groups of different beam (3) Since the UCA radar's working mode is the direct sampling of time-domain signals, it is necessary to perform a Hilbert transformation on the sampled signals and obtain two signals for coherent processing.

Characteristic Factors
The correlation, non-stationarity, and range-Doppler characteristics are selected as the characteristic factors of ground clutter and the characteristic differences in five different scenarios and two beam control modes are analyzed.

Correlation Analysis of Ground Clutter Data
Pearson correlation coefficient (PCC) is used to measure the correlation between different data types [38], and the beam-related correlation function, azimuth correlation function, and range-related correlation function are obtained. Beam correlation refers to the correlation between the sampled data obtained by two groups of different beam directions in the same channel, corresponding to the data samples represented by blue and brown in Figure 5a, one pulse has eight beam correlation coefficients. Azimuthal correlation refers to the correlation of clutter data of different beam directions in the same beam control mode, which corresponds to the data samples represented by the same color in Figure 5a. The range correlation is calculated from data samples of the same color in the same channel (single pulse sampling data of a beam). The case of multi-channel and multi-pulse data is shown in Figure 5b. Data samples with this structure can calculate short-term correlation (intra-frame data) and long-term correlation (inter-frame data). In addition, radar range-Doppler domain data can be calculated from single-frame data of one channel.
Sensors 2021, 21, x FOR PEER REVIEW 7 of 20 directions in the same channel, corresponding to the data samples represented by blue and brown in Figure 5a, one pulse has eight beam correlation coefficients. Azimuthal correlation refers to the correlation of clutter data of different beam directions in the same beam control mode, which corresponds to the data samples represented by the same color in Figure 5a. The range correlation is calculated from data samples of the same color in the same channel (single pulse sampling data of a beam). The case of multi-channel and multi-pulse data is shown in Figure 5b. Data samples with this structure can calculate short-term correlation (intra-frame data) and long-term correlation (inter-frame data). In addition, radar range-Doppler domain data can be calculated from single-frame data of one channel.

Radar Beam Correlation
Each channel of the UCA radar antenna can transmit or receive two beams indifferent directions, and the beam correlation is calculated by the data under different beam in the same channel. When extracting the characteristic factor of beam correlation, the data of different beams in the same channel in one frame are firstly taken to calculate the correlation (see Equation (1)) and the calculated results of eight channels acquired are averaged to obtain the finally calculated results for beam correlation of ground clutter data in one frame. The same pulse contains eight beam-correlation coefficients, so the beam-correlation function is shown as follow: where , and , are the distance unit sampling data of two beams in the same channel, and are the mean value of two groups of data, respectively. When extracting the beam correlation feature factor, firstly, the data of different beams in the same frame are taken to calculate the correlation, then the eight-channel calculation results are averaged to eliminate the influence of noise data. Finally, the beam correlation calculation results of one frame of ground clutter data can be obtained.

Radar Beam Correlation
Each channel of the UCA radar antenna can transmit or receive two beams indifferent directions, and the beam correlation is calculated by the data under different beam in the same channel. When extracting the characteristic factor of beam correlation, the data of different beams in the same channel in one frame are firstly taken to calculate the correlation (see Equation (1)) and the calculated results of eight channels acquired are averaged to obtain the finally calculated results for beam correlation of ground clutter data in one frame. The same pulse contains eight beam-correlation coefficients, so the beam-correlation function is shown as follow: where x i,l and y i,l are the distance unit sampling data of two beams in the same channel, x l andŷ l are the mean value of two groups of data, respectively. When extracting the beam correlation feature factor, firstly, the data of different beams in the same frame are taken to calculate the correlation, then the eight-channel calculation results are averaged to eliminate the influence of noise data. Finally, the beam correlation calculation results of one frame of ground clutter data can be obtained.

Azimuth Correlation
Azimuth correlation function refers to the range unit correlation of pulse data between different channels in the same beam control mode. The antenna beam direction is covered at equal intervals of 360 • , realizing the whole scene range's detection. Therefore, the azimuth correlation function can reflect the differences in different scenes, and it is defined as follows: where x i,l and x i,l+j are the pulse distance unit data with two beam pointing intervals of j channels in the same beam control mode. The beam position interval of the same beam control mode covers the azimuthal direction of 360 • , so the range unit data of adjacent spaced beams have a relatively strong correlation. Azimuth correlation depends not only on the beam position interval, but also on the wave control mode, and the latter has a more significant influence on the azimuth correlation. To distinguish the azimuthal correlation of two-beam control modes, the characteristic factor is obtained by averaging the correlation coefficients of different interval beam positions.

Range Correlation
The range correlation function of ground clutter is usually determined by the spatial difference of each range unit's scattering distribution in the test scene. The correlation calculation formula between each distance unit is as follows: where x i and x i+j are the radar sampling data in range direction with the interval of j between two range units, andm is the mean value of the sampling data of the range unit. The range-direction autocorrelation function gradually decreases with the increase of spatial distance sampling interval related to the beam control mode. When the sampling interval is 10, the average value of eight channels is calculated as the characteristic factor.

Recursive Graph
A recurrence plot (RP) is used to analyze the non-stationary difference of ground clutter data under different beam control modes in the same scene of the UCA radar. Assuming the time series {x(1), x(2), · · · , x(N)}, a set of delay m-dimensional vectors can be obtained by the phase space reconstruction technique. The m-dimensional vector → X(l) is as follows: where → X(l) is the delay vector, and the delay time is τ. At this time, L = N − (m − 1) × τ is the number of reconstruction vectors, and the embedding dimension is m. Based on the embedding theorem [39], the appropriate phase space reconstruction parameters are selected for nonlinear time series. The Euclidean distance between any two points in L's phase spaces is defined as: It can be seen from Equation (5) that X(i) and X(j) are any two points in the phase space, respectively. According to the Euclidean distance of each point in the phase space, the definition of the elements in the recursive graph can be obtained as follow: where i is the number of rows, and j is the number of columns. The distance threshold is ε and θ(·) is the Heaviside function. When d i,j is less than the threshold ε, the value R i,j is 1.
The recursion rate (RR) of the recursive graph of the UCA radar's ground clutter sampling signal in range direction is used as the characteristic parameter to classify the beam control mode. RR refers to the ratio of the number of recursive points in a recursive graph to the total number of points. The calculation formula is: The number of coordinate axis points in the recursive graph is L and R i,j is the current position's recursive value. The larger RR is, the higher the recursion rate is. The RP describes the non-stationarity of the one-dimensional range-direction data, which is also related to the beam control mode. The average of the RR of the eight channels data is calculated as the characteristic factor.

The Range-Doppler Maps
To analyze the statistical characteristics of the near-ground clutter region, the range-Doppler map is effectively clipped. Feature factors, including histogram mean and variance, are extracted from the clipped range-Doppler image's statistical histogram. A frame of data includes range-Doppler data of eight channels in two-beam control modes. That is, it has a total of 16 range-Doppler images. The final feature factor of a frame of data is formed by taking the average of these images' histogram statistical results. Figure 6 illustrates the process of range-Doppler clipping and statistical analysis by taking the dry grassland scene as an example. By extracting the clipping region's sample data, the statistical histogram is made, and the mean value and variance of the sample data are calculated. It can be obtained as follows: where c mn is the amplitude value at positions m and n in the range-Doppler map, and c represents the mean value of clutter amplitude in the clipped region, where the number of distance units is M, and the number of Doppler units is N. One frame of data includes a total of 16 range-Doppler images of eight channels in two-beam control modes. Based on the mean value of statistical histogram results, the range-Doppler domain feature factor of frame data is formed. Sensors 2021, 21, x FOR PEER REVIEW 10 of 20

Feature Factor Extraction and Analysis
Correlation characteristic map and recursive analysis results of range sampling signal of UCA radar ground clutter can be seen from Figures A1 and A2. Three different characteristic factors of radar beam correlation, range-Doppler domain mean, and variance are selected to cluster the ground clutter data of different scenes. Three characteristic factors of the azimuth correlation function, range correlation coefficient, and recursive rate are selected to cluster the ground clutter data of different beam control modes in the same scene (Table 1, takes the ground clutter data obtained from the Wuxi highway scene as an example to calculate the value range of six characteristic factors). The UCA radar has eight channels and 2048 range sampling units. Reading five frames of data, each frame of data pulse number is 64, so a total of 320 pulse samples. Two different beam control modes are distinguished, and the sample size is 640.

Clustering Algorithms
Aiming at the defects of standard genetic algorithm, an improved genetic algorithm is proposed to realize the clutter data clustering function by introducing chaotic sequences' characteristics.

Standard Genetic Algorithm Clustering
The standard GA algorithm constructs the fitness function by searching the minimum clustering center of the objective function value:

Feature Factor Extraction and Analysis
Correlation characteristic map and recursive analysis results of range sampling signal of UCA radar ground clutter can be seen from Figures A1 and A2. Three different characteristic factors of radar beam correlation, range-Doppler domain mean, and variance are selected to cluster the ground clutter data of different scenes. Three characteristic factors of the azimuth correlation function, range correlation coefficient, and recursive rate are selected to cluster the ground clutter data of different beam control modes in the same scene (Table 1, takes the ground clutter data obtained from the Wuxi highway scene as an example to calculate the value range of six characteristic factors). The UCA radar has eight channels and 2048 range sampling units. Reading five frames of data, each frame of data pulse number is 64, so a total of 320 pulse samples. Two different beam control modes are distinguished, and the sample size is 640.

Clustering Algorithms
Aiming at the defects of standard genetic algorithm, an improved genetic algorithm is proposed to realize the clutter data clustering function by introducing chaotic sequences' characteristics.

Standard Genetic Algorithm Clustering
The standard GA algorithm constructs the fitness function by searching the minimum clustering center of the objective function value: where x k is the sample and the clustering as center is m j . Firstly, a population is initialized randomly, and the clustering center is encoded by floating point number. The coding of the clustering center is optimized through selection, crossover, and mutation operations, then the result of clustering division is judged based on the sum of eigenvectors of samples to be classified and the corresponding Euclidean distance of the clustering center. Floating-point coding is to regard a chromosome as a string composed of K clustering as centers, that is, the K clustering analysis of D-dimensional feature vectors and the chromosome structure of the clustering center is defined as: where the feature dimension of the data sample is d, and k represents the central cluster, that is, each chromosome is a floating-point code string with the length of k × d.

Chaotic Genetic Algorithm Clustering
Based on the standard genetic algorithm, we propose an improved idea, which introduces chaotic disturbance terms at the initial position and the local convergence position respectively, which is called chaotic genetic algorithm in this paper. The ergodicity of chaos theory can overcome local minima and achieve global optimization. Adding a chaotic sequence when the population is initialized can give full play to the global optimization ability of the genetic algorithm, and improve the convergence speed of the population by using the ergodicity of chaotic particles. The introduction of chaotic disturbance term at the local convergence position can break the monopoly position of individual dominant individuals in the population, and prevent the group from losing its competitiveness, which not leading to the stagnation of the population.
At the initial position and local convergence position of the standard GA, the logistics equation is adopted to introduce the interference term to generate chaotic sequences, which are determined by the following equation: where z (0) ∈ (0, 1) and z (0) / ∈ (0.25, 0.5, 0.75), when µ = 4, the generated sequence is completely chaotic. The population size (labeled as N) is set to be within 100 and the maximum number of iterations (labeled as T) to be 1000. The evolutionary process of the population is ended by judging the termination conditions. The hybridization probability P c and mutation probability P m are 0.9 and 0.01, respectively. The GA based on chaos theory is obtained by combining chaos theory with GA, and its process is shown in Figure 7.
In the initial stage, the chaotic sequence corresponding to the population size is generated by a chaotic operator to initialize the population. The individual elimination in the process of GA Solver adopts the method of survival of the fittest, that is, the dominant individuals replace the poorest individuals. Chaos genetic algorithm involves selection, crossover and mutation operations, respectively, using roulette algorithm, single-point crossover algorithm, and basic bit algorithm. There are two termination conditions of the algorithm, one is to achieve the set number of iterations, the other is the clustering criterion function to achieve the set threshold. In the process of population evolution, the quantitative calculation is introduced to judge whether the population is precocious, as shown in the following formula: where λ 2 s is the population fitness variance of the population, representing the aggregation degree of individuals of the population. When λ 2 s < C (C is a constant), the population  In the initial stage, the chaotic sequence corresponding to the population size is generated by a chaotic operator to initialize the population. The individual elimination in the process of GA Solver adopts the method of survival of the fittest, that is, the dominant individuals replace the poorest individuals. Chaos genetic algorithm involves selection, crossover and mutation operations, respectively, using roulette algorithm, single-point crossover algorithm, and basic bit algorithm. There are two termination conditions of the algorithm, one is to achieve the set number of iterations, the other is the clustering criterion function to achieve the set threshold. In the process of population evolution, the quantitative calculation is introduced to judge whether the population is precocious, as shown in the following formula: where is the population fitness variance of the population, representing the aggregation degree of individuals of the population. When ( is a constant), the population falls into precocity.
is the fitness value of the current iteration and is the current average fitness of the population.
is the normalization factor.

Results
This section mainly introduces the clustering results of UCA radar clutter measured data and uses the criterion function based on Euclidean distance to evaluate GA clustering's performance and chaotic GA clustering proposed in this paper. The feature factor is normalized according to the maximum and minimum value of the feature factor. In the scene data clustering analysis, five clustering centers are set according to different scenes of measured data. In the clustering analysis of beam control mode data, two clustering centers are set.

Clustering of Clutter Data in Different Scene
The analysis results come from five experimental scenes, and the population size is

Results
This section mainly introduces the clustering results of UCA radar clutter measured data and uses the criterion function based on Euclidean distance to evaluate GA clustering's performance and chaotic GA clustering proposed in this paper. The feature factor is normalized according to the maximum and minimum value of the feature factor. In the scene data clustering analysis, five clustering centers are set according to different scenes of measured data. In the clustering analysis of beam control mode data, two clustering centers are set.

Clustering of Clutter Data in Different Scene
The analysis results come from five experimental scenes, and the population size is set to 10, 15, 20, 25, and 30, and the maximum number of iterations is set to 1000. The process of population evolution is terminated when the central cluster reaches the specified minimum Euclidean distance (the value is actually 20) or exceeds the maximum number of iterations. The clustering experiment in each case is performed five times, and the clustering result is the average value of five times ( Table 2). Compared with SGA clustering, the chaotic SGA clustering proposed in this paper has obvious advantages, which are mainly reflected in two aspects: 1.
It has a faster convergence speed, which can save 34.60% of the time.

2.
It has a higher classification accuracy, and the average criterion function value is reduced by 42.82%. As the population size increases, the classification accuracy of both SGA and chaotic SGA both show an upward trend (from 20% and 40% to 60% and 100%, respectively), while the criterion function decreases significantly (from 49.52 and 23.63 to 28.59 and 21.03, respectively), and the convergence speed of the population gradually slows down (the average convergence speed of SGA clustering increased from 29.85 to 80.31, while chaotic SGA clustering increased from 23.63 to 43.31). Since the introduction of chaotic interference SGA clustering can shorten the time for the population to jump out of the local optimal position, the convergence speed of the chaotic SGA decreases slowly. Besides, the average evolution times of SGA clusters have reached or approached the set maximum evolution generation of 1000, which means that the population is at risk of falling into a local optimum.
From the cluster analysis of the population size, we can see that the larger the population size and evolutionary generation, the better the clustering effect, but the slower the convergence speed. Based on the two-way trade-off between classification accuracy and convergence speed, the population size is set to 20 in the experimental analysis of this study.
When the population size of SGA clustering and chaotic SGA clustering is set to 20, the clustering results and the evolution process of corresponding populations under five different experimental scenes are shown in Figure 8. Compared with SGA clustering, the clustering center obtained by chaotic SGA clustering is more consistent with the classification and division of actual characteristic factors, and the corresponding criterion function values are 24.27 and 19.63, respectively. At the same time, it can be seen from the iterative process of the two algorithms that the average fitness curve and the optimal fitness curve of the SGA clustering evolution process are in the state of fitting for a long time. From the evolution process graph of chaotic SGA clustering, it can be seen that the average fitness function curve shows an upward trend after chaotic disturbances are introduced at the initial position and the local optimal position, respectively. Besides, compared with SGA clustering (55.51 s), chaotic SGA clustering takes less time to find the clustering center (21.04 s).

Clutter Data Clustering of Two-Beam Control Modes in The Same Scene
Taking the highway scene as an example, the clustering analysis results of ground clutter data with different beam control modes are presented in Table 3. The population size is set as 5, 10, 15, 20, and 25, and the maximum number of iterations is set as 500. When the cluster center reaches the specified minimum Euclidean distance (which is actually 75) or exceeds the maximum number of iterations, the population evolution process is terminated. The experiments are repeated five times for each population size, and the clustering performance is evaluated by taking the mean value of the five results.
As can be seen from Table 3, compared with SGA clustering, the chaotic SGA clustering process has a faster convergence speed and smaller criterion function value, with the convergence speed increased by 42.69% and criterion function decreased by 20.79%, while the number of evolutions is less (decreased by 41.69%). The results of several experiments show that the classification effect is good only when the value of the criterion function is lower than 75, while the SGA clustering algorithm is greater than this value even when the population size is 25. Therefore, the SGA clustering algorithm cannot well realize the classification of the two-beam control modes' data. Besides, with the increase of population size from 10 to 25, the average value of criterion function of chaotic SGA is between 73 and 74, and the classification accuracy does not significantly improve with the

Clutter Data Clustering of Two-Beam Control Modes in The Same Scene
Taking the highway scene as an example, the clustering analysis results of ground clutter data with different beam control modes are presented in Table 3. The population size is set as 5, 10, 15, 20, and 25, and the maximum number of iterations is set as 500. When the cluster center reaches the specified minimum Euclidean distance (which is actually 75) or exceeds the maximum number of iterations, the population evolution process is terminated. The experiments are repeated five times for each population size, and the clustering performance is evaluated by taking the mean value of the five results. As can be seen from Table 3, compared with SGA clustering, the chaotic SGA clustering process has a faster convergence speed and smaller criterion function value, with the convergence speed increased by 42.69% and criterion function decreased by 20.79%, while the number of evolutions is less (decreased by 41.69%). The results of several experiments show that the classification effect is good only when the value of the criterion function is lower than 75, while the SGA clustering algorithm is greater than this value even when the population size is 25. Therefore, the SGA clustering algorithm cannot well realize the classification of the two-beam control modes' data. Besides, with the increase of population size from 10 to 25, the average value of criterion function of chaotic SGA is between 73 and 74, and the classification accuracy does not significantly improve with the increase of the population size. However, the number of iterations decreased significantly from 340 to 238. Compared with SGA clustering, chaotic SGA clustering shows a better superiority in both local and global search capabilities.
When the population size is greater than 10, chaotic SGA clustering can realize the effective classification of data under two beam control modes. Set the population size as 15, and the clustering results of SGA and chaotic SGA are shown in Figure 9. The average fitness curve of chaotic SGA increased steadily, and due to the influence of the retraction of the average fitness curve, the stagnation period of the optimal fitness curve in the evolutionary process diagram is relatively short, and the population jumps out of the optimal local position and has relatively strong local search ability (the right of Figure 9b). Although both SGA and chaotic SGA clustering can obtain classification results that match the actual beam control mode, the criterion function values of the clustering center are obviously different (SGA is 90.02, and chaotic SGA is 72.99). At the same time, compared with SGA clustering (21.5373 s, 500), the convergence speed of the chaotic SGA clustering process is faster (12.71 s), and the number of evolution generation is smaller (223). Therefore, chaotic SGA clustering uses a shorter evolutionary generation than SGA clustering to complete the global optimal position search, which has stronger global searchability.

Conclusions
The ground clutter data collected in five field test environments of UCA radar are analyzed, and the effective classification of data in different scenes and two types of beam control modes is completed. In this paper, statistical characteristics (mean value, variance) of the range-Doppler domain and beam correlation feature factors are extracted as scene classification of clutter data. In addition, the azimuth correlation, range autocorrelation, and recursive plot recursive feature factors are extracted as beam control mode classification of clutter data. Finally, the chaotic genetic algorithm proposed in this paper is used

Conclusions
The ground clutter data collected in five field test environments of UCA radar are analyzed, and the effective classification of data in different scenes and two types of beam control modes is completed. In this paper, statistical characteristics (mean value, variance) of the range-Doppler domain and beam correlation feature factors are extracted as scene classification of clutter data. In addition, the azimuth correlation, range autocorrelation, and recursive plot recursive feature factors are extracted as beam control mode classification of clutter data. Finally, the chaotic genetic algorithm proposed in this paper is used for data clustering analysis under five different scenes and two different beam control modes. The conclusions are as follows: 1.
Compared with SGA clustering, the clustering center obtained by chaotic SGA clustering is more consistent with the classification and division of actual characteristic factors. From the scene data, the criterion function values of SGA and chaotic SGA clustering corresponding to scene classification are 38.03 and 21.74, respectively, and the time consumed is 55.57 and 36.34 s, respectively. From the beam control mode of data classification, the criterion functions are 93.99 and 74.45, respectively, and the convergence speeds are 17.47 and 10.01 s, respectively.

2.
Chaotic SGA clustering has high local search ability and global searchability, realizing the effective classification of data samples.

3.
The effective classification and analysis of ground clutter data can improve UCA radar adaptability to clutter environments to enhance target detection performance.
In the future, more ground clutter data of measurement experiments will be carried out to improve the establishment of a ground clutter database in different scenes. It realizes the subtle perception of ground clutter characteristics of UCA radar and provides more space modeling information of ground clutter for the application of auxiliary adaptive detection technology.