1. Introduction
Many radar systems use frequency- or phase-modulated signals as their typical waveforms. Each modulated waveform enjoys its own advantages, along with its disadvantages [
1,
2]. A linear frequency modulation (LFM) and bi-phase-coded signal are good examples of classical radar waveforms. On one hand, LFM has a high Doppler tolerance feature and high compression ratio, but it suffers from time-sidelobe effects. On the other hand, bi-phase-coded signals are known for their low sidelobes, but they have a Doppler-intolerant waveform with a limited compression ratio [
3,
4,
5].
The Doppler effect is a critical aspect of radar measurements, as it hinders the effectiveness of radar signals [
6,
7,
8]. Moreover, especially in the case of high-speed targets, it may distort the signal phase and, therefore, no target information would be extracted. This happens for most bi-phase codes, including Barker codes. If we are able to eliminate or even suppress the Doppler effect, we can improve the quality of the range resolution of phase-coded radar signals.
We believe that our paper is the first to discuss both time-domain Doppler estimation and the signal recovery of bi-phase modulation signals in a radar background. In fact, we performed a thorough literature review, carefully examining existing research, from the oldest to the most recent, to ensure that we did not miss any relevant research in this area.
In [
9], the authors proposed a decision-directed Doppler estimation method based on the frequency average for phase-shift keying (PSK) signals for underwater acoustic communications. This estimation is achieved by averaging the frequency estimates from a bank of correlators. A robust estimate can be further computed by increasing the observation intervals at the expense of the computational cost. This method, however, is limited to underwater targets with very slow speeds.
Doppler estimation is achieved in [
10] using the dual alternating-direction method (ADM) with standard statistical tools and a multiple-hypothesis testing procedure using a few pulses and without any frequency or phase modulation. Although the method introduced a good range and Doppler estimation, it gave only approximate values due to cell accuracy. In [
11], the authors attempted to estimate the real-time Doppler effect using up–down and up–mute–down hyperbolic frequency modulation methods. The estimation is based on finding the relationship between the Doppler scaling factor (the rate of change in the delay of a signal in a time-varying channel) and the time delay of the signal at the receiver. Unfortunately, the method was found inadequate for low-speed targets.
The authors of [
12] managed to compensate for the Doppler effect on complementary phase modulation signals. By creating a new complementary code of length
N and multiplying each sample at index
n by the conjugate of the adjacent sample at
in the new code, the method is able to produce the desired compensated code. Although this approach is straightforward for short codes, it becomes increasingly complex and computationally intensive as the code length increases. Additionally, this method was only tested on targets moving at moderate speeds. Furthermore, after compensation, the Doppler information is lost, preventing the determination of the target velocity.
Another method for Doppler estimation was proposed in [
13], and it relied on a combination of hyperbolic frequency modulation and speed spectrum scanning. Although the method reduces interference in speed spectrum peak detection and improves the accuracy of the estimation, the success of this method requires high sampling rates. In addition, the signal suffers from high time sidelobes, which can affect nearby targets.
Yet another Doppler estimation method can be found in [
14], and it was primarily developed for non-terrestrial networks. The method involves using reference signals at different frequency positions in an OFDM carrier to estimate the Doppler shift for high-speed targets. Using this method, a system can calculate the required frequency adjustment value and use it to pre-compensate uplink transmissions and to maintain multi-access orthogonality at the receiver. Unfortunately, the range of the errors in the estimated Doppler shifts is between
kHz and
kHz, which is not acceptable in many radar applications. Moreover, it is not adequate for radar systems that use a single carrier frequency. In [
15], a Doppler cancellation module that integrates data from multiple signals to eliminate Doppler effects from Doppler-intolerant signals, such as phase-modulated signals, is introduced. This allows the use of high-resolution waveforms to detect high-speed targets, even if they are Doppler-intolerant. However, this method does not provide Doppler value estimation and dictates that one of the integrated signals is used for Doppler estimation.
In [
16], the authors presented an enhanced orthogonal matching pursuit (OMP) algorithm to directly estimate the velocity and the acceleration of the target. They rely on pre-constructing a dictionary for sparse motion parameters of the target. This renders the Doppler frequency migration unnecessary. While effective for estimating the velocity and the acceleration, this method is iterative, relies on the greedy nature of OMP, and can be very expensive, especially for large dictionaries of the parameter space.
The authors of [
17] presented another method for recovering radar waveforms and estimating Doppler values for high-speed targets using two types of waveforms: phase-modulated and unmodulated waveforms. This technique is primarily effective for high-speed targets. For low-speed targets, the authors recommend using traditional Doppler estimation methods, thereby incorporating the proposed method as a sub-module within a typical radar system. Doppler values are estimated using the FFT of the unmodulated signal. Consequently, the Doppler ambiguity is constrained by the repetition time and the resolution of FFT.
It is worth highlighting that, to estimate the Doppler frequency, most of the methods mentioned above rely on DFT, which makes them unable to resolve the exact Doppler frequency due to DFT’s inherent Doppler ambiguity. Moreover, many of the available techniques for estimating Doppler frequency are not suitable for bi-phase signals.
In this paper, we present an efficient and inexpensive algorithm that accurately estimates the Doppler value by relying on just one pulse. The algorithm enhances the robustness of its estimate by using an ensemble averaging of Doppler estimators based on the next few successive pulses. To increase the accuracy even further, we rely on an iterative technique to progressively remove any Doppler residue. Since the estimation is performed in the time domain, the proposed algorithm is immune to all Doppler ambiguities. Thereby, the Doppler ambiguity that occurs in the conventional Doppler calculations does not affect our algorithm. Furthermore, the proposed algorithm segregates the Doppler effect from the signal, removes the Doppler effect from any bi-phased waveform using a sequence of subsystems, and enables the return of the waveform to its well-known shape at low signal-to-noise-ratios (SNRs) to produce a reliable range resolution. By doing so, we decrease the time and hardware complexities from which the conventional radar systems suffer, as they rely on expensive MTD processors with a large coherent pulse interval (CPI). In addition, we circumvent the main bottleneck problem in using bi-phase codes for moving target detection, especially for high-speed targets.
Figure 1 shows a block diagram of our proposed modules in a conventional radar system.
As the block diagram depicts, the proposed module can be viewed as an add-on component that sits between the ADC and the rest of the conventional radar receiver modules. Instead of feeding the received signal directly to the matched filter (MF), we first feed it to a sequence of three subsystems (Doppler estimator, phase corrector, and waveform recovery) that form our proposed add-on module. The Doppler estimator applies a nonlinear transformation to the received signal to remove the bi-phase waveform and calculate an estimate for the Doppler value. The phase corrector uses this estimate to reduce the Doppler effect from the received signal and then check if any remainder of the Doppler effect needs to be suppressed. If that is the case, the adjusted signal from the phase corrector is re-fed to the Doppler estimator for further processing. This treatment is repeated until an accurate estimate of the Doppler is obtained. The waveform recovery subsystem uses the boosted estimates of the Doppler value along with the received signal to detect the target speed and recover the original waveform with a minimal Doppler effect. The recovered waveform is then fed to the conventional range radar detector.
It is worth mentioning that the proposed module does not affect the transmitter and is powerful enough to enable the designer to optionally disregard MTD, which requires many pulses to obtain an accurate Doppler estimate without resolving the Doppler ambiguity.
This paper is structured as follows: in
Section 2, the algorithm is described analytically, providing an overview of its key components. The simulation results of the proposed algorithm are presented in
Section 3, highlighting its performance in a simulated environment.
Section 4 investigates the performance of the system and analyzes its effectiveness in practical scenarios. The validation of the theoretical model and the implementation of the algorithm using software-defined radio (SDR) are discussed in
Section 5. Furthermore,
Section 6 identifies the limitations of the current algorithm and offers potential solutions to overcome them. Finally, in
Section 7, the paper concludes by summarizing the key findings and discussing future work that can be pursued in the radar field.
2. Analytical Description of the Proposed Algorithm
In this section, the theory of the proposed algorithm is introduced. A block diagram of the proposed module is shown in
Figure 2 and is used to illustrate the main stages of our algorithm.
First, we generate a complex baseband bi-phase signal and denote it by
. The bi-phase signal is then sent through the traditional modules of a radar transmitter. At the receiver side and after demodulation, the return signal has the following form:
where
is the target delay,
is the target Doppler effect,
is the noise, and
is the delayed bi-phase waveform. As a pre-processing step, the proposed module uses median filtering, a nonlinear digital filtering technique, to reduce noise and improve the performance of the subsequent processing. This filter is used as it effectively preserves edges during the denoising process [
18,
19]. With an appropriate window size, the effect of this filter on the signal and noise can be assumed to have the following form:
where
and
are real numbers with
. The filtered signal is then subjected to a nonlinear transformation
, which enables us to eliminate the effect of the bi-phase and isolate the Doppler effect in the time domain. Although nonlinear transformations were also used in [
20], it was primarily to compensate for very low-frequency offsets. In our case, we apply it for the blind estimation of larger Doppler shifts, starting with a coarse estimate, which is then refined in subsequent iterations. Given that we use a bi-phase waveform, we choose
to be an even power function, in particular, a quadratic one, which we apply in all subsequent iterations of the algorithm.
To illustrate the effect of the transformation, we choose
to be a quadratic function and apply it to
, as shown below.
where
. From this equation, we can see that the nonlinear transformation
serves two objectives: (1) it cancels the bi-phase waveform from the signal, as the term
represents the bi-phase waveform, which oscillates between
and
depending on its phase. When we apply the nonlinear transformation, the term
always results in 1, cancelling the bi-phase signal, and (2) it turns the multiplicative nature of the Doppler effect into an additive form. In other words, at the receiver, the received signal consists of the delayed signal multiplied by the Doppler effect, in addition to the noise. After transformation, the signal can be represented by three terms, as Equation (
3) shows. The first term corresponds to the scaled Doppler signal, the second term represents the pure scaled noise effect, and the last term accounts for the scaled signal multiplied by both noise and the Doppler effect. Given that
, after filtering and transformation, the effect of the last two terms is reduced by
and can be neglected. This shows the importance of introducing a filter before applying the nonlinear transformation. By combining the last two terms and carrying out a few manipulations, the effect of the nonlinear transformation can also be written as
Therefore, the signal is effectively converted from a Doppler multiplicative form into a pure Doppler additive form. This simplifies our task of extracting Doppler information from the return signal. Although the second term,
, has been reduced at the output of the nonlinear denoising filter, we still apply a smoothing filter (to be precise, a moving average filter) before the Doppler estimator to mitigate the effect of the nonlinear transformation process on the undesirable term,
, while reducing the number of noisy, close-by peaks in the first term.
After the smooth filtering, the transformed signal can be considered to have the following form:
where
and
are
with
. The signal is now primed for Doppler estimation and bi-phase signal recovery, as we discuss next.
2.1. Doppler Estimation
The estimation of the Doppler value starts with the local minmax block (see
Figure 2). The local minmax algorithm detects the locations of the minima and maxima and assigns them to the arrays
and
, respectively. This is expressed as follows:
Let us denote by
and
the sizes of
and
, respectively. The quantities
and
derived from the minmax algorithm reflect the density of local minima and maxima in the signal. As the target velocity increases, the frequency modulation induced by the Doppler effect results in a higher density of waveform changes. Consequently, the algorithm detects more minima and maxima within a given time interval, leading to an increase in the size of the
and
arrays. Three cases present themselves.
Case (1):
and
. This case can be divided into four sub-cases:
and
,
and
,
and
, and
and
. For the first three sub-cases, the Doppler frequency is calculated as follows:
where
n is the power exponent of the nonlinear transformation, and
is the estimated Doppler effect from the first incoming pulse. If the indices are used, the expression becomes
where
is the sampling frequency, and
and
are the sampling indices of the elements in the arrays
and
, respectively.
An interesting sub-case is when both
and
are strictly greater than 1. In this case, we can improve the accuracy of our estimated value based on the ensemble of the estimates as follows:
where
.
Case (2):
and
, or
and
. Since we only have a single minmax location, we apply the nonlinear transformation again to increase the number of minmax locations and then apply the above case using Equation (
9). In this case, the estimate becomes
Case (3): and . This indicates either that the original signal has a very low Doppler value and the algorithm should directly feed it to the MF or that the Doppler residue is so small that we can terminate the iterations and proceed to signal recovery before the range calculation. In this case, the Doppler estimate is set to 0.
All of the above cases rely on a single pulse. To increase the accuracy of the estimate, an ensemble average can be taken over a few successive pulses as follows:
where
N is the number of pulses used, and
is the Doppler estimate for the
ith pulse according to the formulas in the cases above.
2.2. Doppler Checker, Phase Corrector, and Bi-Phase Signal Recovery
In each iteration, the estimated Doppler value,
, is passed to the Doppler checker to decide whether to continue the iterations or proceed to range and velocity calculations, as shown in
Figure 2. If
, the Doppler checker passes
to the accumulator, which aggregates the Doppler estimates from the previous iterations.
Since the Doppler estimate can exceed the actual value, the aggregation is not always obvious. In fact, the Doppler checker has to be novel enough to compute the proper sign of before passing the signed to the accumulator.
The accumulated Doppler value
is delivered to the conjugate of the Euler form and Hadamard product blocks to remove the estimated Doppler effect from the signal. Mathematically, the conjugate of the Euler form block computes
, while the Hadamard product, an element-wise multiplication between two signals, outputs the following:
where
is the phase effect of the target Doppler delay. Note that
is a constant, as the Doppler effect and range of targets are unchanged from scan to scan. Therefore, this phase does not affect the local minmax detection.
If
is 0, the Doppler checker knows that the current Doppler residue is small enough to terminate the iterations and deliver the signal for range and velocity calculations, as shown in
Figure 2 with yellow arrows. The ultimate effect of each iteration is, therefore, to reduce the overall Doppler effect
by an amount equal to the accumulated Doppler values
.
Interestingly enough, the expression of
in Equation (
12) is analogous to that of
in Equation (
2) but with a reduced Doppler effect. Therefore, we can recursively repeat the same processing that we used for
on
until we obtain a Doppler estimate value
.
3. Simulation Results of the Proposed Algorithm
In our simulations, we use the parameters shown in
Table 1. The inner workings of the different modules depicted in
Figure 1 are detailed in the flowchart presented in
Figure 3. On the receiver side, the incoming signal is down-converted by the RF front-end and then digitized by the ADC. A median filter is then applied to the digitized signal.
Figure 4a depicts a single-pulse Barker-13 with a coherent pulse interval (CPI) at an SNR of
dB. After applying the median filter, we obtain the signal shown in
Figure 4b. The figure illustrates the effectiveness of the median filter and its ability to significantly decrease the noise while preserving the edges, including the sharp ones at the points where the signal undergoes a phase flip of
. Performing this step is critical for both achieving accurate estimates of the Doppler frequency and minimizing the number of required iterations.
In the first iteration, a nonlinear transformation is applied to the filtered waveform. In our algorithm, we chose a quadratic function to remove bi-phase modulation from the signal and to increase the Doppler cycles. The output from this stage is shown in
Figure 5. It is worth pointing out that although the median filter was effective in reducing noise and preserving edges, it was not able to maintain sharp transitions in the bi-phase waveform. This was due to the nature of the median filter, which approximated sharp transitions differently. This, in turn, caused distortions in the Euler form at the output of the nonlinear transformation. The effect of these distortions was minor, as the Doppler estimator relied on ensemble averaging in most cases.
A moving average filter was employed before the local minmax detector to increase its robustness by smoothing the noisy, close-by peaks, as shown in
Figure 6. While it may appear that the MA filter does not make a substantial difference, it is important to note that before applying the moving average in
Figure 5, the signal exhibits noticeable oscillations. If these are mistakenly considered as peaks by the proposed system, it would significantly affect the Doppler estimation. Therefore, the function of the moving average is to smooth these oscillations and ensure a more accurate Doppler estimation.
The local minmax detector is analogous to the traditional peak-finding algorithm but is also able to detect valleys. The indices from the local minmax detector are then used for Doppler estimation and phase correction, as discussed in
Section 2.1 and
Section 2.2.
Figure 7 shows the signal
obtained from the Hadamard block in the first iteration. By using Equation (
12),
Figure 7 illustrates how the signal undergoes modulation by both fast and slow signals, which is, of course, caused by the Doppler effect. According to the second term in Equation (
12), the noise is modulated by a fast signal due to the estimated Doppler value. This value is very large, as our algorithm was able to accurately approximate the actual Doppler value just from the first iteration. The enlarged portion of the image over a
ms time interval shows the modulated noise. The median filter used in the proposed algorithm can effectively remove this kind of noise, as illustrated by the red curve in the figure.
The figure also shows a second slow modulation associated with the first term in Equation (
12). This slow modulation was due to the Doppler residue after removing the estimated Doppler effect. After filtration, the resulting signal still retains the residual Doppler value, as shown in the same figure using a red color. To show the effect of the second iteration, this signal is also plotted in
Figure 8a.
In the second iteration,
becomes the signal that is shown in
Figure 8b. The new accumulated value of
is
KHz, which is very close to the actual Doppler value of 15 KHz. In the third iteration, our local minmax detector did not find any local minima or maxima; therefore, it terminated the iterations and fed the estimated Doppler value to the block performing the waveform recovery and computing the target velocity.
The novelty of the proposed algorithm is that it converges in just a few iterations. In our case, it only took two iterations to achieve an accurate estimate of the Doppler value.
Figure 9, the output of a Matlab workspace, concludes the above description.
The MF output from the recovered waveform, compared to that of the conventional bi-phased pulse compression (PC) radar system, is shown in
Figure 10. The output demonstrates the superiority of our algorithm in estimating the exact Doppler value and perfectly recovering the signal compared to the conventional system.
We evaluated the robustness of our new module by testing it at very low Doppler values, such as 150 Hz, and comparing it to the conventional PC radar system. The Matlab workspace results, shown in
Figure 11, indicate that our module estimated the Doppler value to be 154 Hz, with a difference of just 4 Hz from the actual Doppler frequency.
Although the conventional system was able to detect the target with a constant false alarm rate (CFAR), the signal was degraded both in the mainlobe and the near sidelobes, as shown in
Figure 12b, compared to the MF output of our module depicted in
Figure 12a.
4. Performance Evaluation
To assess the performance of our proposed Doppler estimation and signal recovery system, two types of evaluations are discussed, namely, Doppler and detection evaluations. The first evaluation measures the SNR against different Doppler values, and the accuracy of the module is then calculated. The second evaluation measures the system SNR against the probability of detection. These evaluations offer a valuable understanding of how effective the proposed module is in practical and real-world applications.
Prior to evaluating our proposed algorithm, it is important to explain the reasoning behind categorizing Doppler values into low, medium, and high intervals, as well as the criteria used to select these intervals. Note that all calculations are performed in the S-band.
4.1. Selection Criteria for Doppler Value Intervals
The algorithm that we developed is capable of handling a wide range of Doppler values, from hundreds of Hz to several thousand Hz. We categorized these values based on the target speed as follows:
Low-Doppler targets: This category includes various types of drones, such as consumer drones, racing drones, and some types of helicopters. Consumer drones typically have speeds ranging from 14 to 20 m/s, with advanced models reaching up to 30 m/s. Racing drones, designed for speed and agility, can exceed speeds of 45 m/s. Some types of helicopters, which are slower and used for specific tasks, also fall into this category. Consequently, the Doppler values for this category range from 100 Hz to 2000 Hz [
21,
22].
Medium-Doppler targets: This category covers fixed-wing aircraft and commercial airplanes. Their average speeds range from 100 m/s to 350 m/s, resulting in Doppler values between 2000 Hz and 7000 Hz [
22,
23].
High-Doppler targets: This category is dedicated to military aircraft, which have speeds ranging from 500 m/s to 750 m/s. The corresponding Doppler values for this category range from 10,000 Hz to 15,000 Hz [
24].
To cover the full Doppler interval from 0 to 15 KHz, we choose the following Doppler ranges:
Low-Doppler Range: 100–2000 Hz;
Medium-Doppler Range: 2000–10,000 Hz;
High-Doppler Range: 10,000–15,000 Hz.
4.2. Selected Threshold for Doppler Values
The Doppler value threshold should be selected such that the radar receiver can extract the target information without the need for our proposed algorithm. Therefore, if the proposed module is activated and the residue is less than the specific threshold, we can still obtain target information.
Based on the parameters listed in
Table 1, the pulse width of the phase-coded waveform equals
ms. Hence, the Doppler completes one cycle during the pulse duration at 770 Hz. To avoid signal distortion, the Doppler value should be less than half of this value. To ensure good quality of the signal in the presence of the Doppler effect at the receiver end, the threshold is set to one-third of this value, approximately 256 Hz. Therefore, the Doppler value threshold is chosen to be 250 Hz.
4.3. Accuracy Evaluation
A cumulative distribution function (CDF) plot was selected to gain a deeper insight into the distribution of estimation errors across various Doppler values and SNR levels.
Figure 13 shows the CDF plot of these estimation errors. The CDF was calculated based on randomly generated Doppler values from the low (100 Hz–2 kHz), medium (2–10 kHz), and high (10–15 kHz) ranges, with SNRs spanning from −12 dB to 0 dB in 1 dB increments, and the calculations were repeated 1000 times each for a total of 39,000 trials.
A steeper slope in the CDF curve indicates more occurrences of smaller errors, particularly for values less than 100 Hz, with a 99.5% probability. This reflects the algorithm’s ability to achieve precise Doppler estimation. Additionally, the plot reveals instances of larger errors occurring less frequently, with values of up to 1000 Hz having a 0.5% chance of occurring. This, however, can be addressed via an ensemble average over a few extra pulses.
To further gauge the accuracy of the proposed module, we rely on box plots. These plots allow us to examine the central tendency and spread of Doppler estimation errors at low, medium, and high Doppler values with a fixed SNR of −3 dB. In each case, the Doppler frequency to estimate is randomly selected and then tested 100 times.
Figure 14 shows the results obtained. Remarkably, our algorithm consistently exhibited excellent performance across all three categories, with the absolute median estimation error not exceeding 30 Hz. It is important to note that the proposed algorithm has the ability to detect both positive and negative differences in all cases. The algorithm is designed to consistently check for these differences to prevent subtracting a value from a negative result, which could otherwise lead to an increase in the output error.
Table 2 provides a quantitative representation of the overall accuracy achieved by incorporating the proposed module into the radar system. The accuracy assessment was conducted based on the predetermined threshold value of 250 Hz (see
Section 4.2), representing the maximum allowable frequency deviation to ensure the extraction of the main signal features without distortion. Additionally, the module was tested across a range of SNR values between −12 dB to 0 dB to provide a more comprehensive evaluation of the system’s performance. The obtained results demonstrated notable accuracy levels, with a 100% accuracy rate being achieved at low Doppler values, 98.23% at medium Doppler values, and 96.32% at high Doppler values. These findings indicate the efficacy of the proposed module in maintaining high accuracy levels across a wide range of Doppler scenarios and SNR values, thereby validating its suitability for the intended application.
4.4. Probability of Detection Evaluation
The performance of the radar system with and without the proposed module is illustrated by the probability of detection (
) vs. SNR curves, which are depicted in
Figure 15. To maintain a probability of false alarms (
) of
, a typical value in radar systems, we used a CFAR processor. For each SNR value, we calculate the
value 500 times and plot the average. Doppler values from the low, medium, and high ranges were chosen to show the robustness of the proposed module. In fact, the module exhibited similarly excellent responses for all of these Doppler values. For an SNR greater than −12 dB, the three curves converge and practically achieve
detection. As the SNR decreases, the three curves diverge slightly, with medium and high Doppler values being better than the lower ones.
Comparing the detection curves of the proposed algorithm with the reference curve without the Doppler effect (indicated by the rhombus marker) reveals a degradation of 3 dB when compared to the medium-Doppler-value curve (indicated by the circle marker) and 5 dB when compared to the low-Doppler-value curve (indicated by the cross marker). It is important to note that this degradation is considered acceptable when compared to the curve for which the Doppler effect (indicated by the pentagon marker) was applied but without using our module; the curve shows no detection at any given SNR values.
5. Proof of Concept
Traditional radar systems typically rely on dedicated hardware components for signal processing, waveform generation, and data analysis. However, by leveraging software-based processing in radar systems, such as in software-defined radar (SDR), multiple radar functions can be consolidated onto a single hardware platform, reducing equipment and maintenance costs [
25,
26,
27,
28]. This offers increased flexibility and adaptability, allowing for the rapid development and deployment of new radar functionalities. SDR platforms can be easily reprogrammed to support different radar modes, waveforms, and signal processing algorithms, enabling radar systems to be highly configurable [
29,
30,
31].
In this section, we present a proof of concept and a prototype for our proposed module in the radar system using SDR. The system uses two SDRs, one for transmitting the radar signal and holding the selected Doppler value and the other for receiving the transmitted signal. We also incorporated a synchronizer in our system by adding a marker signal, such as LFM or hyperbolic frequency modulation (HFM), to accurately estimate the beginning of the bi-phase signal at each CPI. We chose the LFM signal as our synchronizer because of its Doppler-tolerant nature. It is worth noting that we did not extract any information from the LFM signal; it was only used as a synchronization marker.
We provide the setup parameters of our experiment in
Table 1, with the following adjustments based on the specifications of the SDRs used:
We evaluated the performance of our module at the three different Doppler values of 8, 10, and 15 KHz, and we obtained the results shown in
Figure 16. Note that we selected these specific values to test the system under challenging conditions within the medium- and high-Doppler-frequency ranges. The results indicated that our module accurately estimated the Doppler values with minor differences from the actual values. More specifically, for the Doppler value of 8 KHz, the estimated value was 7.987 KHz, with a difference of only 12 Hz; for 10 KHz, the estimated value was 10.031 KHz, with a residue of 31 Hz; for 15 KHz, the estimate was 15.02 KHz, with a difference of 20 Hz. These results align with the simulation results recorded in
Figure 13. Additionally, our module successfully recovered the well-known shape of the signal with all Doppler values, as observed in the MF output shown in
Figure 16b,d,f. This result reflects the robustness and effectiveness of our proposed module in real-world applications, making it a valuable addition to any radar system.
7. Conclusions
A novel and low-computational-complexity algorithm for Doppler estimation and signal recovery in bi-phase radar systems has been proposed. The algorithm is designed to accurately estimate the Doppler frequency and remove the Doppler effect from the received signal, allowing for precise range and velocity calculations.
The algorithm consists of several sub-modules, including a median filter, nonlinear transformation, local minmax detector, Doppler calculation module, Doppler checker, phase corrector, and bi-phase signal recovery. These sub-modules work together to iteratively estimate and reduce the Doppler effect from the received signal until the Doppler residue is zero or sufficiently small.
The excellent performance of the proposed system in estimating the Doppler effect and recovering the signal at different SNRs and Doppler values has been demonstrated. In addition, the implementation using SDR technology showed great potential.
Investigating the possibility of dynamically adapting the radar parameters, expanding the algorithm to handle scenarios with multiple targets simultaneously, and testing the algorithm under different background effects would be valuable extensions that will be the focus of our work in the near future.