Rotating Target Detection Using Commercial 5G Signal

: Passive radar detection emerges as a pivotal method for environmental perception and target detection within radar applications. Through leveraging its advantages, including minimal electromagnetic pollution and e ﬃ cient spectrum utilization, passive radar methodologies have garnered increasing interest. In recent years, there has been an increasing selection of passive radar signal sources, and the emerging 5G has the characteristics of a high-frequency band, high band-width, and a large number of base stations, which give it signi ﬁ cant advantages for use in passive radar. Therefore, in this paper, we introduce a passive radar target detection method based on 5G signals and design a rotating target speed measurement experiment. In the experiment, this paper validated the method of detecting rotating targets using 5G signals and evaluated the measurement accuracy, providing a research foundation for passive radar target detection using 5G signals and detecting rotating targets such as drone rotors.


Introduction
With the continuous development of communication technology and the integration of sensing and communication systems, the types of space radiation signal sources are increasing, improving passive radar detection technology.Radio signal candidates range from analog signals, such as television and frequency modulation (FM) signals, to digital communication signals, such as the global system for mobile (GSM) communications and wireless fidelity (WiFi) signals.As the signal bandwidth and carrier frequency gradually increase, detection performance also continues to improve.Moreover, with the development and application of the integration of fifth-generation mobile communication technology (5G), the use of 5G signals as passive radar signal sources for target detection also has research significance.
In recent years, the emergence of signal sources such as GSM, WiFi, and long-term evolution (LTE) signals has led to rapid development in the short-range application of passive radar.Although they have low power, they are widely present in the environment and have a large signal bandwidth, making them suitable for short-range detection.Therefore, these signals are applied to passive radar systems for target detection in different scenarios [1][2][3][4][5].Additionally, the bandwidth and carrier frequency of communication signals are gradually increasing, which will improve the detection performance of passive radar, indicating the direction for its development.The rapidly developing 5G is also being applied to passive radar [6][7][8].
At present, 5G is in a stage of rapid development, fostering widespread accessibility to its corresponding resources.These resources can be used for passive radar systems, and many passive radar researchers often focus on larger targets such as vehicles and people.When the target is small, targets such as drones or drone rotors, whether passive radars based on 5G can correctly detect the target is a worthwhile research question.Therefore, based on existing research, this paper proposes a method for detecting small targets over short distances using the Channel State Information-Reference Signal (CSI-RS) in 5G.Additionally, a rotating target experimental model employing a stepper motor is constructed to accurately simulate target movement scenarios.To conduct experimental measurements, the CSI-RS in the 5G signal is applied in different experimental scenarios.Then, by comparing the measurement results with the set parameters of the experimental system, the correctness of the method of measuring rotating targets using 5G signals is verified.In a word, this paper provides an approach for small target detection based on 5G passive radar.

Related Work
This chapter provides a detailed introduction to the related work mentioned in the first chapter.This work is classified into the following three categories based on signals.
GSM and LTE: Samczynski et al. proposed the concept of GSM passive radar for road vehicle monitoring and verified its feasibility, providing an example for other GSM-based applications [1].As for LTE, the feasibility of using LTE signals for target detection has been explored [4].
WiFi: Colone et al. designed a passive radar data processing method based on WiFi signals and verified the feasibility of the system through experiments [2].In 2017, Colone et al. applied the WiFi passive radar system to the detection and monitoring of targets such as aircraft and vehicles at small airports, proving the potential of the system in such scenarios [3].In terms of human perception, Xie et al. proposed a respiratory rate monitoring system based on WiFi signals, verifying the ability of passive radar based on WiFi [5].These studies reveal the potential of passive radar for short-range applications and also provide different signal-processing methods.For example, Fabiola Colone et al. obtained target information by setting signal monitoring channels and reference channels, while Xie et al. obtained WiFi channel state information ratio(CSIR) for target detection.
5G: Kanhere conducted a simulation analysis on the method of using 5G signals for target localization and evaluated the performance of applying 5G in radar systems [6].Księżyk et al. provided a detailed introduction to the basic principles of 5G sensing and conducted concept validation, explaining the feasibility and limitations of using 5G signals for passive radar and providing guidance for further research on passive radar based on 5G [7].Abratkiewicz et al. also studied the use of Synchronization Signaling Block (SSB) signals for target detection in addition to classical passive radar signal processing methods, and validated this method with cars as experimental targets, providing a new approach for target detection [8].

5G Signal Overview
The 5G system is the evolution of the 4G system and meets the needs for higher data transmission rates, greater data traffic, low latency, and large-scale connections.Additionally, 5G new radio (NR) is a key technology in 5G, which is a global 5G standard based on OFDM's new air interface design.Table 1 presents a comparison of the main parameters of 4G and 5G NR, including waveform, bandwidth, and subcarrier spacing.  1 that both the uplink and downlink signals of 5G signals can use cyclic prefix orthogonal frequency division multiplexing (CP-OFDM) waveforms and the physical layer uplink can also use discrete Fourier transform-spread orthogonal frequency division multiplexing (DFT-S-OFDM) waveforms.OFDM technology converts high-speed information streams into lower-rate parallel data streams through serial-toparallel conversion and then modulates the parallel data with multiple orthogonal subcarriers.The 5G signal used for detection in this paper is the CP-OFDM waveform.When the cyclic prefix is not considered, the OFDM-modulated signal can be expressed as follows: ( ) where OFDM T is the symbol period, N is the number of subcarriers in the OFDM modulation process, c f is the carrier frequency, and i d is the amplitude of the i-th modulation symbols.When the carrier frequency is not considered, the baseband signal waveform ( ) x t can be obtained from (1), which can be represented as follows: ( ) Additionally, the sampled signal can be expressed as follows: It can be seen that ; thus, OFDM modulation can be achieved using Inverse Discrete Fourier Transform.

Features of 5G Signal Physical Layer
The method used in this paper is to obtain channel response through CSI-RS signals in 5G for object detection.Therefore, it is necessary to understand the main parameter settings of 5G and CSI-RS signals.This chapter introduces the main parameter configuration of 5G and CSI-RS signals.

Cyclic prefix (CP)
CP-OFDM waveforms can be used for uplink and downlink transmission of 5G signals.After obtaining the OFDM modulation waveform, a segment of data at the end of each symbol is copied to the front end to obtain the CP-OFDM waveform, which is the cyclic prefix.The cyclic prefix prevents the previous symbol from falling into the range of this symbol due to multipath effects and does not destroy the orthogonality between subcarriers.When inter-symbol interference occurs, the presence of the cyclic prefix causes a cyclic shift in the signal, resulting in a phase shift in the demodulated symbols.The 5G system has flexible subcarrier settings and corresponding cyclic prefix types, as listed in Table 2 [9], where μ is the subcarrier parameter.
When the subcarrier spacing is 30 kHz, each subframe contains two slots, that is, each slot is 0.5 ms [11].Additionally, when the type of cyclic prefix is a normal cyclic prefix, a slot contains 14 OFDM symbols; when it is an extended cyclic prefix, a slot contains 12 OFDM symbols [11].In the 5G protocol, many time-frequency resource units are defined, such as the resource element (RE), resource block (RB), physical resource block (PRB), and time-frequency resource grid.The RE is the smallest resource unit in 5G, occupying one OFDM symbol in the time domain and representing one subcarrier in the frequency domain.An RB contains 12 subcarriers in the frequency domain and is not defined in the time domain.PRB is a physical RB within the bandwidth part (BWP).BWP is a partial continuous bandwidth configured for UE in the transmission channel.Within each BWP, PRB is independently numbered, and, in order to manage different BWPs, a common resource block (CRB) is defined in 5G.The schematic diagram of BWP is shown in Figure 1 [10], where

CSI-RS Signal
This paper uses the CSI-RS in 5G to detect rotating targets.CSI-RS [11][12][13][14] is a reference signal used in 5G signals to obtain channel state information.It is similar to a pilot signal.In 4G, the cell-specific reference signal (CRS) is used to obtain the downlink channel status between the base station and the user.And CRS can be configured with 1, 2, and 4 antenna ports.Compared to 4G, the functions of CSI-RS in 5G have been expanded, including time-frequency tracking, channel state information obtaining, and more flexible configuration of CSI-RS parameters such as the number of antenna ports, signal period, and multiplexing method.Additionally, the parameter configuration of CSI-RS is very important for using it for target detection.
The basic parameters of CSI-RS signals mainly include the number of antenna ports, period, and density.The number of antenna ports supported by a CSI-RS signal includes 1, 2, 4, 8, 12, 16, 24, and 32.The period of CSI-RS can be 4, 5, 8, 10, 16, 20, 32, 40, 64, 80, 160, 320, and 640 time slots, and the length of the time slot is related to the length of the CP.CSI-RS has three densities: 0.5, 1, and 3.This paper obtains channel response through the CSI-RS signal, so it is necessary to know the position of the CSI-RS signal in the resource grid.The position of the CSI-RS signal in this experiment is mainly determined by period, density, symbol location, and subcarrier location.The period determines the frequency of CSI-RS occurrence in the time domain.Density determines how many subcarriers CSI-RS occupies within a PRB, symbol location determines the time starting position of CSI-RS in a slot, and subcarrier location determines the subcarrier starting position of CSI-RS within a PRB. Figure 2 shows a CSI-RS RE mapping, which includes a slot in the time domain and a PRB in the frequency domain.When the subcarrier spacing is 30 kHz, a slot contains 14 OFDM symbols.The number of antenna ports is one and the density is three.Three CSI-RS symbols occupy three REs in a PRB.In Figure 2, the symbol location is five, starting from the fifth OFDM symbol, and the subcarrier location is one, starting from the first subcarrier.

The Bistatic Radar Model
Utilizing 5G signals for rotating target detection essentially employs bistatic radar principles for target identification.The transmitting antenna and receiving antenna of the bistatic radar are separated from each other, and their structural model is shown in Figure 3 [15].Where v is the velocity of the target, L is the distance between the receiving and transmitting antennas, T L is the distance between the target and the transmitting an- tenna, R L is the distance between the target and the receiving antenna, β is the double base angle, and δ is the angle of the target at the velocity direction relative to the bisector of the double base angle.Assuming the wavelength of the signal is λ , the Doppler frequency shift generated by the moving target can be expressed as [16] ( )  The signal received by the receiving antenna can be represented as follows:

Unilateral rotation target
The frequency offset curve of the target in different rotation directions is shown in Figure 5b when the speed of the unilateral rotating target is 2 rps.Except for two special cases where the target rotation center is on the centerline of the receiving and transmitting antenna line and the double base angle is zero, the Doppler frequency deviation curve generated by a uniformly rotating target resembles a sine curve, and its maximum and minimum values are not equal in size.Their size relationship is related to the target area and rotation direction.If the target area is known, the change in the target rotation direction can be determined based on it.

Bilateral rotating target
The bilateral rotating target model can be established as shown in Figure 6.Set the position relationship as shown in Figure 3 with a speed of 2 rps.When the radius of both sides of the target is 0.3 m, the target frequency offset curve is as shown in Figure 7a.When the bilateral target radii are 0.3 m and 0.15 m, the target frequency offset curve is as shown in Figure 7b.

Target Detection Processing Method
For passive radar signal processing, the time-frequency cross-correlation method is often used to detect targets [17].Generally, two signal-receiving channels are set up, namely, a reference signal channel for receiving direct waves and a monitoring channel for receiving target echoes.Then, adaptive filtering or extended cancellation algorithms are used to filter out direct waves and multipath clutter [18][19][20].This paper mainly studies methods of target detection using channel response.The CSI-RS in the 5G signal can be used to estimate the channel between the transmitting and receiving antennas and obtain the Doppler frequency shift of the rotating target based on the channel response.When the 5G signal configuration is known, the corresponding reference signal can be generated locally to obtain the transmission signal.Rotating the target causes the channel between the transmitting and receiving antennas to change.The channel response is obtained by transmitting and receiving signals, and the information about the rotating target is obtained from the channel response.

Channel Estimation
The frequency domain least squares (LSs) method is a commonly used channel estimation method in OFDM systems.For OFDM systems that rely on a pilot symbol, the transmitted pilot symbol vector is P X and the received pilot symbol vector is ( ) To obtain the minimum value of ( ) , set its partial derivative to 0, i.e., ( ) , and the estimated channel response value is [21].
The flow chart of target frequency offset estimation based on 5G signals is shown in Figure 8 [22].( ) ( ) ( ) where T is period which includes the basic symbol length and the cyclic prefix length, ( ) ( ) ( ) According to the LS method, the obtained channel response can be expressed as follows: In this paper, CSI-RS is applied to obtain channel responses, and the Equation ( 11) can be represented as follows: where CSI T is the period of the CSI-R, and is the interval between adjacent subcarriers of the CSI-RS.The dimension of the channel response matrix is 1 1

N M
× , where 1 N is the number of subcarriers occupied by the CSI-RS in the frequency domain, and 1 M is the number of OFDM symbols occupied in the time domain.

Doppler Frequency Offset Estimation
According to Equation ( 12), there will be a phase shift caused by a Doppler frequency shift between different symbols on the same subcarrier.Therefore, the Doppler frequency shift produced by the rotating target can be calculated through the phase shift between symbols.Then, the rotation speed of the rotating target can be estimated based on the obtained Doppler frequency deviation versus the time curve.
In this experiment, the carrier frequency is much larger than the signal bandwidth, so the target Doppler frequency shifts obtained by different subcarriers can be approximately equal.To obtain the target Doppler frequency offset, the channel response obtained at different single subcarrier positions can be used for processing, that is, a timefrequency analysis of ( ) . However, due to frequency selectivity, the channel response obtained on a single subcarrier may have lower energy and greater interference.It is also possible to stack the channel responses at multiple subcarrier positions in the frequency domain before conducting a time-frequency analysis, i.e., conducting a time-frequency analysis on ( ) . The time-frequency analysis serves to perform a short-time Fourier transform on the obtained channel response to obtain frequency domain information and time domain information at the same time.The time-frequency analysis schematic is shown in Figure 9.In data processing, it is necessary to choose an appropriate window length to achieve better results.The channel response extracted using the n-th subcarrier can be represented as n H , and the result of the time-frequency analysis on n H can be represented as follows: The size of the matrix ( ) × , where the column vector represents the result of the FFT transformation at the corresponding time j .
The global maximum value index is obtained by detecting the amplitude of the FFT result corresponding to a certain time j , and the frequency corresponding to the index is the target Doppler frequency at that time, which can be expressed as follows: arg max abs , where ( ) ⋅ represents the index corresponding to the maximum of the k-th along the vector index i .When the number of targets is known, the local maximum values can be directly sorted to obtain the frequency offset of all targets.Inversely, if the number of targets is unknown, the algorithm of constant false alarm rate [23,24] detection can be applied to generate an adaptive threshold; thus, the local maximum value exceeding the threshold will be determined as a target.

Interference Suppression 1. Direct wave and multipath clutter suppression
The signals received by the receiving antenna include target echoes, direct waves from the transmitting antenna, and multipath signals reflected by other obstacles in the environment.After the signal is synchronized, the delay and frequency offset of the direct wave are both 0. According to Equation ( 12), the direct wave will add a constant to the corresponding symbol and subcarrier position in the matrix

( )
, H n m .The delay of the multipath signal relative to the direct wave is τ , and the frequency offset is 0, which is expressed as a constant complex number in ( ) , H n m .So, a strong peak will appear in the spectrum at zero frequency.Therefore, it is necessary to suppress the zero-frequency component of the signal to obtain the Doppler frequency offset of the target, which can be achieved by removing the average value of the signal or using a filter.
Using the rotation model shown in Figure 5, the target echo intensity is set to 0 dB, the direct wave intensity is set to 20 dB, and the signal-to-noise ratio is 15 dB.In the simulation, the period of the CSI-RS signal is set to five time slots, that is, 50 ms.Figure 10a,b shows the time-frequency analysis results and extracted frequency offset curve after removing the signal mean value.Figure 10c shows the error between the actual measurement value and the theoretical value.It can be seen that the error is small and random, which has a small impact on the measurement of the frequency offset variation period.For the filter, the equal ripple design method is used to design a bandpass filter with an order of 226.After filtering, delay compensation is performed on the signal.The timefrequency analysis results and the extracted frequency offset curve are shown in Figure 11a-c and show the error between the actual measurement value and the theoretical value.The reason why the starting position in Figure 11b is 0 for tens of milliseconds is that the filter initially acts on the signal, resulting in poor filtering performance.In data processing, it is possible to discard data with poor filtering performance at the forefront of the wave.

Random phase noise suppression
In the experiment, the signal-receiving and transmitting antennas are separated, which would lead to the problem of clock desynchronization.Clock desynchronization may result in random phase shifts that need to be processed to obtain more accurate Doppler shifts.Multiple antennas can be used to receive signals in the experiment.If two receiving antennas share a crystal oscillator, the clock offset of the two antennas at the same time will be the same.Taking advantage of this characteristic, the cross-antenna crosscorrelation (CACC) method or the cross-antenna signal ratio (CASR) method can be used to eliminate random phase shifts [25,26].In this paper, we use the CASR method to eliminate random phase shifts.
Taking the channel response extracted from a subcarrier as an example, the channel response can be divided into static path components and dynamic path components.For the i-th antenna, the obtained channel response can be written as follows [25]: where ( ) i s h t is the static path component in the i-th antenna channel response, ( ) ( ) h t is a constant in the dynamic path component, and ( ) t θ is a random phase offset.For antennas i and 1 i + using the same oscillator, the signal ratio can be expressed as follows [26]: For a rotating target, it is approximately fixed in a short period, during which channel states π have the same periodicity, so the CASR method can remove phase noises while utilizing the periodic changes in phase to obtain Doppler frequency shift information.

Summary of Detection Methods
In summary, the overall process of using 5G signals for rotating target detection can be obtained as shown in Figure 12.When a 5G base station is used for measurement, it is necessary to obtain the configuration parameters of the CSI-RS signal.Firstly, demodulate the 5G signal obtained by the receiving antenna to obtain the resource grid, and then extract the CSI-RS signal based on the CSI-RS configuration parameter to obtain the channel response.The next step is to suppress random phase noise and zero frequency components in the channel response.The signal after interference suppression is subjected to time-frequency analysis, and the target frequency offset is extracted from the results of the time-frequency analysis.After smoothing the frequency offset curve, the rotating target speed is obtained through FFT.The key to the method proposed in this paper is to understand the parameters of CSI-RS signals and obtain the parameters of 5G base stations.If the configuration parameters of 5G base stations cannot be directly queried, other methods need to be found to obtain them.
To make this method more widely used, it is necessary to develop applications that can obtain configuration information of any 5G base station and perform signal demodulation in order to further integrate it with the signal processing method in this paper.

Experiments and Data Analysis
In experiments, the CSI-RS signal period is 20 ms, and the maximum unambiguous Doppler frequency is 50 Hz.The measurable Doppler frequency shift range is [−25 Hz, 25 Hz].For a unilateral rotation target, when the signal carrier frequency is 3.55 GHz, the rotation radius is 0.3 m, and the double base angle is 0, the maximum measurable speed is 0.56 rps.The presence of the double base angle makes the maximum measurable speed slightly greater than 0.56 rps.The experiment is conducted in two measurement scenarios.

5G Experimental Base Station Testing
In the first measurement scenario, an experiment is conducted using a 5G base station in a 5G laboratory.The CSI-RS signal parameters of 5G base stations in the laboratory are shown in Table 4.The measurement scenario is shown in Figure 13a, and the model is shown in Figure 13b.

Unilateral Target
The radius of the unilateral rotation target is 0.3 m, and the data with a speed of 0.5 rps and clockwise rotation are analyzed as an example.To obtain better display results, the initial 1 s of data are discarded.Figure 14a shows the time-frequency analysis results of the channel response without interference suppression, Figure 14b shows the time-frequency analysis after random phase noise suppression, and Figure 14c shows the timefrequency analysis after suppressing the zero frequency component.Figure 14d shows the smoothed frequency offset curve extracted by peak detection, in which the blue solid line represents the actual measurement results, and the red dashed line represents the theoretical values derived from the rotating target model.Figure 15a shows the error between the measured Doppler frequency and the theoretical value.The initial large error is due to the poor filtering effect of the filter when it first acts on the data.It can be seen that the overall error is small, and the experimental results are consistent with the theoretical model.The errors are mainly caused by model measurement errors and the FFT spectral resolution.From the spectrum graph Figure 15b, it can be observed that the measured speed is 0.500488 rps, with an error of 0.1%.For the frequency offset curves corresponding to different rotation directions at the same position, as shown in Figure 15c, it can be seen that the relationship between the maximum and minimum values of the frequency offset curves corresponding to different rotation directions is different.

Bilateral Rotating Target
In the bilateral rotating target experiment, data measurements are conducted on two types of bilateral targets as shown in Figure 16a,b.Taking the counterclockwise speed of 0.375 rps as an example for analysis, the time-frequency analysis results are shown in Figure 17a when the model has the same radius.The radii of the two targets are both 0.3 m.The result of frequency offset extraction is shown in Figure 17c and the spectrum is shown in Figure 17e.The obtained speed is 0.757 rps, which is twice the actual speed.The reason for this is that, when the radii are the same, the frequency offset curves of the two targets are the same, with only half a cycle of time delay.When the radii of the two targets are 0.3 m and 0.15 m, respectively, the time-frequency analysis results are shown in Figure 17b,d, and the corresponding spectrum is shown in Figure 17f.The obtained speed is 0.378 rps, with an error of 0.8%.

5G Commercial Base Station Testing
Based on the 5G experimental base station experiment, a measurement experiment using a 5G commercial base station is also designed.In this measurement scenario, the accuracy of target speed estimation is verified using a commercial 5G base station.Unilateral rotating target models with radii of 0.3 m and 0.15 m are tested.The measurement scenario is shown in Figure 18a, and the obtained CSI-RS signal parameters of the 5G base station are shown in Table 5.The measurement results of the unilateral target rotating model with a radius of 0.3 m are shown in the red line in Figure 18b and Table 6.The measurement results of the unilateral target rotating model with a radius of 0.15 m are shown in the blue line in Figure 18b and Table 7.When the rotation radius of the unilateral target is 0.3 m and the theoretical speed is 0.75 rps, the reason for the error of 100% is that the Doppler frequency exceeds the measurement range, resulting in Doppler blur.

Analysis of Computational Load and Processing Time
This paper uses Matlab for signal processing, and the computer CPU used is an Intel(R) Core(TM) i5-9300H CPU.According to the processing flowchart shown in Figure 12, Table 8 shows the computational complexity of each step.

Number of Floating Point Multiplication Description
Random phase noise suppression The computational load and processing time are mainly determined by the amount of data used in data processing.The data processing in this paper generally takes several seconds.Taking the data in Section 5.1.1 as an example, 1 400 N = , , the data processing takes 5.891 s.

Conclusions
WiFi, 4G LTE, and other signals are widely used in passive radar.The rise of 5G has led scholars to study passive radar based on 5G.In addition to the classic processing methods of passive radar, some scholars have studied signal processing methods based on certain signals within 5G, such as signal processing methods based on SSB signals.This paper focuses on whether other signals in 5G can be used for target detection and proposes a method of using CSI-RS signals of 5G for channel estimation and extracting Doppler frequency offset from channel responses.Meanwhile, compared to other studies using passive radars based on 5G for vehicle detection or personnel localization, the method in this paper focuses more on detecting weak and small targets at short distances.
In order to verify the correctness of the method proposed in this paper, the rotating target was chosen as the detection object, as the actual motion state of the rotating target, such as velocity, can be more easily obtained.In this paper, a rotating target echo model is established, including bilateral rotating targets and unilateral rotating targets.Additionally, the frequency offset curve characteristics of different types of rotating targets are analyzed.The impact of interference signals and receiver phase noise in actual measurements is also analyzed.At the same time, a rotating target model for actual measurement is established.This paper conducted measurement experiments in different environments and parameter configurations, and compared the measurement results with theoretical values, providing data for research on passive radar based on 5G.According to the experiment results, the speed estimation error does not exceed 5% within the measurable range.The results have verified the feasibility and correctness of the detection method for rotating target detection and velocity estimation using CSI-RS signals.This paper can provide a foundation and direction for further research and application of passive radar based on 5G.The limitations of this method and its application in other scenarios can be further studied in the future.
of the cell containing this BWP relative to CRB0, of BWP relative to CRB0.The resource grid contains two dimensions, the time domain and frequency domain, and it is a method of signal representation.Each grid represents a subcarrier in the frequency domain and an OFDM symbol in the time domain.

Figure 1 .
Figure 1.The schematic diagram of BWP.

Figure 3 .
Figure 3. Bistatic radar model.TX and RX are the transmitting and receiving antennas, and the speed of the target is v .β is the double base angle, and δ is the angle between the direction of the target velocity and the bisector of the double base angle.
) 4.1.2.The Model of the Received Signal In the actual environment, the receiving antenna will receive the direct wave signal from the base station and the multipath interference caused by the reflection of other obstacles in the environment while receiving the target echo.The transmission diagram of the 5G signal from the transmitting base station to the receiving antenna is shown in Figure 4.

Figure 4 .
Figure 4. Signal propagation diagram.The 5G base station is the signal transmission location, and the receiving antenna is the signal arrival location.The direct wave is the signal that directly reaches the receiving antenna from the 5G base station.Multipath signals are signals reflected by static obstacles.The target echo is reflected by the target.

.
) where ( ) s t is the transmitted signal, 0 t is the time delay of the signal propagating in a straight line between the receiving antenna and the base station, 0 m t τ + is the time delay generated by the transmitted signal reflecting from the m-th stationary obstacle to the receiving antenna, 0 n t τ + is the time delay generated by the transmitted signal reflecting from the n-th moving target to the receiving antenna, n d f is the Doppler frequency shift caused by the radial velocity generated by the motion of the n-th moving target, and ( ) n t is noises in the environment.4.1.3.Rotating Target Measurement Scene The rotating target measurement scene established based on the bistatic radar model is shown in Figure 5a.TX is the 5G signal-transmitting antenna, and RX is the signal-receiving antenna.0 T L is the projected distance from the target rotation center to the trans- mitting antenna on the antenna connection line, and 0 R L is the projected distance from the target rotation center to the receiving antenna on the antenna connection line.r L is the distance from the rotation center to the antenna connection, and r is the target rotation radius.In the simulation, the carrier frequency is 3.45 GHz, The frequency offset curves of different types of rotating targets are discussed below.

Figure 5 .
Figure 5. (a) Rotating target measurement scene.TX is the signal-transmitting antenna, and RX is the signal-receiving antenna.The target rotates at a constant speed of radius r under the drive of the pole.T L , R L , L , r L , 0 T L , and 0 R L are distances between different points; (b) unilateral ro- tating target Doppler frequency offset curve.

Figure 6 . 2 rFigure 7 .
Figure 6.Bilateral rotating target model.Target 1 and Target 2 are two targets with rotation radii 1 r N , where p N is the noise vector.The objective function of the method can be expressed as follows:

Figure 8 .
Figure 8. Flow chart of target frequency offset estimation based on 5G signal.The transmission resource grid of 5G signals can be represented as TX G , where TX G is a matrix of size N M × , and M is the number of OFDM modulation symbols contained in the signal.The baseband modulation signal can be represented as follows: is the frequency of the n-th subcarrier in the baseband signal, and ( ) rect ⋅ is the rectangular window function.Assuming the time delay of the target echo signal is τ and the Doppler frequency shift is d f , the received signal can be represented as follows: is the transmission resource grid of the signal,
f j is the target Doppler frequency offset index at the corresponding time j , ( ) abs ⋅ represents the modulus value, and ( ) arg max ⋅ represents the index corresponding to the global maximum.For multiple targets, the global maximum cannot detect all targets.Instead, local maximum detection can be used to extract the frequency offset of multiple targets by sorting the local maximum or setting a threshold, which can be represented as follows: the Doppler frequency offset index of the k-th target at time j , and

Figure 10 .
Figure 10.(a) Time-frequency analysis results after removing the signal mean value; (b) frequency offset curve after removing the signal mean value; (c) the error between the actual measurement value and the theoretical value.

Figure 11 .
Figure 11.(a) Time-frequency analysis results after filter processing; (b) frequency offset curve after filter processing; (c) the error between the actual measurement value and the theoretical value.

Figure 14 .Figure 15 .
Figure 14.(a)Time-frequency analysis results of channel response without interference suppression; (b) time-frequency analysis results after suppressing random phase noise; (c) time-frequency analysis results after suppressing zero frequency component; (d) smoothed frequency offset curve of unilateral target.

Figure 16 .Figure 17 .
Figure 16.(a) Bilateral rotating target model with the same radius; (b) bilateral rotating target models with a different radius.

Figure 18 .
Figure 18.(a) 5G commercial base station testing measurement scenario; (b) results of unilateral target with a radius of 0.3 m and unilateral target with a radius of 0.15 m.

1 4M
Add up the channel responses obtained from 1 N subcarriers; the channel re- Signal has J real sampling points; the order of the LPF filter is 31

Table 4 .
CSI-RS parameter table for 5G base station in 5G laboratory.

Table 5 .
CSI-RS signal parameters of the 5G commercial base station.

Table 6 .
Results of unilateral targets with a radius of 0.3 m.

Table 7 .
Results of unilateral targets with a radius of 0.15 m.

Table 8 .
Analysis of computational complexity for each processing step.