Next Article in Journal
Proposal and Implementation of a Procedure for Compliance Recognition of Objects with Smart Tactile Sensors
Next Article in Special Issue
Dynamic Gesture Recognition Based on FMCW Millimeter Wave Radar: Review of Methodologies and Results
Previous Article in Journal
New Biosensor for Determination of Neuropilin-1 with Detection by Surface Plasmon Resonance Imaging
Previous Article in Special Issue
A Novel Approach for Simulation of Automotive Radar Sensors Designed for Systematic Support of Vehicle Development
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Double Fourier-Transform Imaging Algorithm for a 24 GHz FMCW Short-Range Radar

Department of Information Engineering, Electronics and Telecommunications, Sapienza University of Rome, 00184 Rome, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(8), 4119; https://doi.org/10.3390/s23084119
Submission received: 10 March 2023 / Revised: 14 April 2023 / Accepted: 16 April 2023 / Published: 19 April 2023
(This article belongs to the Special Issue Recent Advances of FMCW-Based Radar Sensors)

Abstract

:
A frequency-modulated continuous-wave radar for short-range target imaging, assembling a transceiver, a PLL, an SP4T switch, and a serial patch antenna array, was realized. A new algorithm based on a double Fourier transform (2D-FT) was developed and compared with the delay and sum (DAS) and multiple signal classification (MUSIC) algorithms proposed in the literature for target detection. The three reconstruction algorithms were applied to simulated canonical cases evidencing radar resolutions close to the theoretical ones. The proposed 2D-FT algorithm exhibits an angle of view greater than 25° and is five times faster than DAS and 20 times faster than the MUSIC one. The realized radar shows a range resolution of 55 cm and an angular resolution of 14° and is able to correctly identify the positions of single and multiple targets in realistic scenarios, with errors lower than 20 cm.

1. Introduction

In the last years, radar systems have been widely used in defense, automotive, navigation, traffic control, and environmental monitoring sectors, as well as medicine and sports [1]. The main advantages of radars, with respect to ultrasound (sonar) and optical (LIDAR) systems, are their relative insensitivity to weather conditions (rain, snow, dust, etc.), the possibility to accurately identify target movements, and the ability to penetrate non-metallic materials. Various types of radars, such as Doppler [2,3], ultra-wideband [4,5], and frequency-modulated continuous-wave (FMCW) [6,7] radars, were proposed in the literature.
An FMCW radar for the berthing of large ships in inclement weather with poor visibility was proposed in [8]. The radar operates at a central frequency of f0 = 9.25 GHz, with a bandwidth of B = 100 MHz and employs 4 transmitting (TX) and 16 receiving (RX) antennas, synthesizing, in this way, a uniform linear array of 64 virtual antennas. Image reconstruction was achieved by employing the beam-steering technique for different identification distances. A multi-beam and multi-range FMCW radar (f0 = 24 GHz, B = 100 MHz) was presented in [9]. The proposed radar can be used for cruise control and stop-and-go systems. A two-channel TX and four-channel RX structure was designed and connected to series-fed patch antennas. In order to extract angular information with high accuracy, the subspace-based digital beam forming (DBF) with a multiple signal classification (MUSIC) algorithm was used. In [10], a surveillance FMCW radar (f0 = 24 GHz, B = 2 GHz) to monitor car and pedestrian traffic, as well as intruders, was designed. The developed radar system consists of one transceiver with two selectable TX antennas, eight parallel RX antennas, and a back-end module for data logging and control. In the target detection algorithms, the ranges and relative velocities of multiple targets were obtained by means of 2D Range-Doppler-FFT processing, while the angular information was determined by digital beam forming based on the received signals.
Further applications of FMCW radars concerned the monitoring of physical activity of subjects [11,12], the classification of hand movements using convolutional neural networks [13,14], the localization and monitoring of vital signatures of subjects [15], and the recognition of the smoothness of railway tracks [16].
Various types of FMCW transceivers are available on the market that operate in the 24 GHz industrial, scientific, and medical (ISM) band [17,18,19,20,21,22]. Distance2Go is a 24 GHz radar based on the BGT24MTR11 transceiver chipset from Infineon [17]. The radar board, equipped with one TX and one RX antenna, can measure the distance, speed, and direction of movement. These features make the radar suitable for various applications, such as motion detection, presence sensing, proximity sensing, etc. [18,19]. The EV-TINYRAD24G from Analog-Devices is a 24 GHz radar evaluation module that allows for the implementation and testing of radar-sensing applications [20]. This radar has two TX and four RX integrated antenna array in a multiple-input/multiple-output (MIMO) radar scheme. The combined antennas allow the TinyRad module to detect the distance, speed, and angular position of multiple targets simultaneously. The IMST 24 GHz FMCW radar Development Kit DK-sR-1200e is provided with high level human-tracker and range-finder algorithms. The radar has one TX and two RX channels equipped with integrated patch antennas [21] and was used to monitor the evolution of lake ice thickness and snow mass [22].
In the mentioned papers, very complex and integrated radar realizations were presented. Moreover, only one reconstruction algorithm was used, and no comparison between different techniques was reported.
In this work a simple and modular FMCW radar for the control of subjects on the sidewalks of train and subway stations and for monitoring elderly subjects in their home environment was designed and realized. Since, for the abovementioned applications, fast acquisition times are required, a system with one TX and four RX antennas was chosen. A new, fast, and simple reconstruction algorithm based on double-fast Fourier transform (2D-FT) was proposed and compared with already available delay and sum and MUSIC algorithms [23]. To evaluate the overall performance of the designed radar system, both the simulated and experimental results are presented.

2. 2D–FT Reconstruction Algorithm

To generate a two-dimensional image of a scenario by using radar data or numerical simulations, suitable reconstruction algorithms must be used. Five different reconstruction algorithms, used to process data acquired with a frequency step radar (SFR), were presented in [23]. In the radar proposed in [23], the scattering parameters (S21) between M selectable pairs of transmitting (TX) and receiving (RX) antennas (MIMO radar) were measured for a given number (N) of frequencies.
Since FMCW radars use IQ demodulators, the two IF outputs generated by mixing the received and transmitted signals are proportional to the real and the imaginary part of the transmission scattering parameter, which, in turn, are proportional to the target radar cross-section. In conclusion, the reconstruction algorithms presented in [23], with reference to a frequency step radar, can also be applied to the present case concerning an FMCW radar. Two different algorithms, the delay and sum (DAS) and the multiple signal classification (MUSIC) are used in the following. Moreover, since the fast Fourier (FT) algorithm described in [23] does not work properly with few antennas (M = 4 in this case), a new double Fourier transform (2D-FT) intensity function is introduced to overcome this drawback.
The normalized 2D-FT intensity function can be derived by the DAS intensity function given by the following:
I D A S x , y = F D A S x , y F D A S M a x
F D A S x , y = f 0 B 2 f 0 + B 2 L 2 + L 2 Γ ξ , f e + j 4 π f c x ξ 2 + y 2 d ξ d f
where F D A S M a x is the maximum value returned by F D A S x , y in the Region of Interest (ROI), f 0 is the central frequency of the operating band B, and L is the length of the array centered in the origin of the Cartesian reference system (see Figure 1).
A second-order Taylor expansion of the phase term appearing in the integrand function in (2) is introduced to reduce its numerical computational burden. By doing so, as well as by employing the approximation reported in Appendix A, we obtain the following:
F D A S x , y f 0 B 2 f 0 + B 2 L 2 + L 2 Γ ξ , f e + j 4 π f c x 2 + y 2 j 4 π f 0 c x ξ x 2 + y 2 d ξ d f = = f 0 B 2 f 0 + B 2 L 2 + L 2 Γ ξ , f e + j 2 π f τ + j 2 π γ ξ d ξ d f γ = 2 f 0 c x x 2 + y 2 τ = 2 c x 2 + y 2 = F 2 D F T x , y
Equation (3) allows us to define the intensity function, I 2 D F T x , y , based on the 2D inverse Fourier transform through the following analytical link between the variables (γ, τ) and the observation point (x, y):
γ = 2 f 0 c x x 2 + y 2 τ = 2 c x 2 + y 2
thus obtaining
I 2 D F T x , y F 2 D F T x , y F 2 D F T M a x
where F 2 D F T M a x is the maximum magnitude value returned by F 2 D F T x , y in the ROI.
Figure 2 shows the mapping of the ROI in the corresponding region expressed by means of the γ and τ variables, where
γ m i n = 2 f 0 c x m a x x m a x 2 + y m i n 2
γ m a x = 2 f 0 c x m i n x m i n 2 + y m i n 2
τ m i n = 2 y m i n c
τ m a x = 2 d m a x c
and d m a x is the maximum distance of the ROI points from the origin of the Cartesian reference system.
The mapping from the γ and τ domain to the spatial variables x and y is expressed as follows:
x = c τ γ 4 λ 0 y = c τ 2 1 λ 0 γ 2 2
If the following condition is met,
λ 0 γ 2 1
which, according to (4), implies
x m i n y m i n 1 x m a x y m i n 1
then (10) can be adequately approximated as follows:
x = γ λ 0 2 y y = c τ 2
Note that (12) corresponds to a ROI located far from the Radar system array, i.e., a situation of practical interest.
The numerical implementation of the proposed 2D-FT reconstruction algorithm is reported in Appendix B.

3. Radar DESIGN

Figure 3a shows the schematic of the designed modular radar. The EVAL-BGT24MTR11 board by Infineon was used as the transceiver [24]. The BGT24MTR11 chip contains a VCO operating between 24 and 25.5 GHz, power and low-noise amplifiers, and an I-Q mixer (see Figure 3b). The scaled by 16 board output frequency (Q1) is sent to the VCO input of the EVAL ADF4159 board by Analog Devices used as PLL circuit (see Figure 3c) [25]. The charge pump output (CP) of the board is filtered by the loop filter and sent to the VCO control input (FINE/COARSE). By properly setting the registers of the ADF4159, the BGT24MTR11 board generates a chirp signal in the ISM band between 24 and 24.25 GHz, with a fixed time duration. One of the two balanced outputs of the transceiver is sent to the transmitting antenna (TX), while the other (TXX) is terminated on a 50 Ω load. The RF output power of the board is about 11 dBm, with a power consumption of 500 mW in continuous operating mode. The low-noise amplifier noise figure is 12 dB, and the mixer conversion gain is equal to 26 dB.
The signals coming from the 4 receiving antennas are sent to the SP4T Hittite board EVAL HMC1084LC4 [26]. This board can select through an electronic control (control voltage ranging between 0 and −3 V), one of the four outputs that is sent to the BGT24MTR11 input (RFIN). The generated intermediate signals (IFI and IFQ) are acquired through the NI-USB 6361 DAQ card [27] and transferred to a PC to be used as input for the image reconstruction algorithms.
Various design solutions have been proposed in the literature for the radar radiating system, such as horn, Vivaldi, or dielectric antennas [28,29,30,31,32,33]. In this study, a simple radiating system based on a serial patch antenna array [34] was chosen to ensure high gain, light weight, and low cost (see Figure 4).
The realized serial patch array antenna is printed on a Roger RO4350B dielectric substrate (𝜀𝑟 = 3.66) with a thickness of 30 mils (762 μm), and a loss tangent of 0.0037. Since the array currents have to flow in the same direction, the feeding point is located as shown in Figure 4a. In this way, the microstrip line located below the feeding point shifts the current phase by 180° at the antenna central frequency, and therefore, on the two pairs of three patches, the current flows in the same direction. A taper on the patches’ width (along the radiating edge of each patch) was applied to reduce side lobe levels [35,36]. The taper factors are 0.75 and 0.5 for the last two patches located in the opposite directions. To take into account the parasitic effects introduced by the SMA connector-antenna transition, a fine-tuning using the Microwave Studio full-wave software by CST was performed. The optimized dimensions of the antenna microstrip lines are LS = 3.5 mm and WS = 0.5 mm. All the patches have the same length, LP = 3.12 mm, while the widths are WP1 = 4.01 mm, WP2 = 3.0075 mm, and WP3 = 2.005 mm. The antenna radiation pattern, simulated with the full-wave CST software, is reported in Figure 4b and evidences a directivity of 14.2 dBi, with a side-lobe level of −14 dB. This value is 4 dB better than the one achieved simulating a uniform-width patch array and ensuring good rejection of environmental clutter.
The overall antenna system consists of a serial patch array antenna as transmitter (TX in Figure 5) and 4 serial array antennas, one wavelength spaced, as receiver (RX in Figure 5). The antenna was realized starting from the Gerber file generated by the CST software by using a numerically controlled milling table (Quick Circuit 5000).
The scattering parameters of the TX and RX antenna system, simulated by CST and measured using the Agilent PNA E8363C network analyzer, are reported in Figure 6. A good agreement between measurements and numerical results can be observed. In particular, the measurements show an impedance bandwidth (for a −10 dB threshold) of the transmitting antenna of about 600 MHz (between 23.8 GHz and 24.4 GHz) and a TX/RX coupling isolation greater than 35 dB.
The ADSimPLL software [37] was used to design the chirp shape and the loop filter components. The PLL parameters (maximum, minimum, reference and phase detector frequency, VCO tuning law, etc.) were considered to be the ADSimPLL inputs. In particular, since the PLL is controlled by a 16-time-scaled output frequency of the BGT24MTR11 board, the VCO central frequency is 24,125/16 = 1507.8 MHz, and the bandwidth 250/16 = 15.625 MHz. To reach a good compromise between the carrier phase noise and its settling time, a loop bandwidth of 45 kHz with a phase margin of 47° was chosen. To this end, a 3-pole 4th-order loop filter was designed using the ADSimPLL software. The filter schematic is reported in Figure 7a, while its experimental prototype with surface mount devices (SMDs) is shown in Figure 7b. The filter transfer impedance was analyzed by using Microwave Office software by Cadence. Figure 8 shows a comparison between the simulated results and the measurements performed on the filter prototype from which an average error of the order of 0.9 dB, corresponding to 11%, can be extrapolated. An excellent agreement between measured and simulated responses is evidenced.
To start up the radar, the registers of the ADF4159 PLL board are initialized by using the Analog Devices ADF4158-9 PLL software. With this tool, the radar chirp, identified with ADSimPLL software, was synthetized by means of 128 frequency steps, with a deviation of 122 kHz and a time duration of 8.12 μs. In this way, the time chirp duration is 1039 μs, and the total ramp frequency shift is 15.625 MHz.
To verify the actual operation of the PLL circuit used in the radar modular board, the ramp driving the VCO at the output of the loop filter was measured (see Figure 9), showing a chirp repetition frequency of 962 Hz, corresponding to a period of 1039 μs. As can be observed, the agreement between measurement and numerical results achieved from using the ADSimPLL software is good.
The phase noise on the tones synthesized by the PLL was estimated with the ADIsimPLL software. At the scaled output Q1 of the BGT24MTR11 block scheme shown in Figure 3, a Single-Sideband-to-Carrier Ratio (SSCR) of −105 dBc/Hz at 10 kHz offset frequency was obtained.
Finally, the NI-USB 6361 DAQ card was used to select the EVAL HMC1084LC4 switch outputs. A suitable voltage-level adapter is used to adapt the TTL output of the DAQ card to the inputs of the EVAL HMC1084LC4 board. All the operations of the NI-USB 6361 board are controlled by a PC on which a LabVIEW virtual instrument is running. The radar prototype, the used boards, the voltage adapter used to drive the SP4T switch, and the laser distance meter used as reference for the target distance are shown in Figure 10.

4. Numerical and Experimental Results

4.1. Test of the 2D-FT Algorithm on Canonical Cases

A first example of 2D-FT reconstruction test, based on the simulation of a point target located at xs = 0.5 m, ys = 3.0 m and with the ROI located between xmin = −2 m, xmax = 2 m, ymin = 1 m, and ymax = 5 m, employs an array of Na antennas spaced by da = λ0/2 (see Figure 11).
Under the hypothesis of point-like targets located far from the antenna array, the reflection coefficients Γ p q in (A17) take the following form:
Γ p q σ   e x p j 4 π f q r p c
where σ is a parameter linked to the target radar cross-section and to the electrical path of the RF signal, while r p is the target-to-antenna distance (see Figure 1).
To identify how the accuracy of the reconstructed target position depends on the fractional bandwidth and on the radar carrier frequency, three different cases were analyzed. In the first case, a radar operating at the central frequency of f0 = 24.125 GHz with a B = 250 MHz [FBW = 0.01 (1%)] was considered; in the second case, the radar exhibits a B = 2.4 GHz [FBW = 0.1 (10%)] at f0 = 24 GHz; and in the last case, the radar frequency is f0 = 2 GHz and B = 0.2 GHz [FBW = 0.1 (10%)]. Numerical simulations were performed by considering Na = 4 and Na = 16. The achieved error parameters, E Q and E L , defined in Appendix A, are reported in Table 1 and Table 2, respectively.
From the above tables, it appears that a critical condition for the method ( E Q M a x >> 1) is obtained with f0 = 2 GHz, FBW = 0.1, and Na = 16. To establish the suitability of these parameters in identifying the goodness of the method in finding the correct position of the target, the 2D-FT algorithm was applied, and the normalized intensity function versus position was reported in Figure 12a, along with the cross-range direction for the 4 antennas and 16 antenna cases. The figure outlines that, with 16 antennas, the 2D-FT technique is not able to identify the target position. In all other cases, the target is correctly identified. As an example, Figure 12b shows the cross-range cut of the normalized intensity function for the case f0 = 24 GHz, FBW = 0.1. In this case, both with 4 and with 16 antennas, the target position is correctly identified.
The tables also show that the condition based on the maximum average error in the entire ROI is less stringent and better suitable to characterize the performance of the 2D-FT method. Note that the quadratic term in (A2) is responsible for a worsening accuracy when the frequency is lower for the same FBW, while an increase in accuracy is observed when the antenna number decreases.

4.2. Target Identification with Simulated Data

The DAS, 2D-FT, and MUSIC algorithms were compared by considering the simulated scenario reported in Figure 13.
A point target (red dot) is placed 3 m far from the antenna plane. The radar, which is equipped with an array of Na = 4 antennas spaced at da = λ0/2, exhibits a bandwidth of B = 250 MHz at 24.125 GHz.
Figure 14 shows the cuts of the normalized intensity function along the range (a) and cross-range (b) directions performed by using the three considered techniques.
The figures show that all the considered reconstruction techniques are able to correctly localize the target position with errors lower than 1 cm. In particular, the DAS and 2D-FT algorithms present a similar behavior, while the MUSIC, having super resolution properties, better identifies the target position (see Figure 14). For the considered radar, the theoretical resolution in range is δr = c/2B = 60 cm, while the angular resolution is δθ = λ/(2∙da∙(Na − 1)) = 0.33 rad (19°).
The half-power beam width (HPBW) in range and cross-range for DAS and 2D-FT is of 55 cm and 70 cm, respectively, which is very close to the 60 cm and 100 cm of the theoretical resolution (see Figure 14). In contrast with the two considered algorithms, the HPBW in range and cross-range of the MUSIC algorithm is of about 2 cm.
The considered algorithms required a computational time (CT) (MATLAB programs ran on a Pentium i7 processor) equal to the following: 2D-FT, 0.62 s; DAS, 2.94 s; and MUSIC, 14.5 s. Consequently, the 2D-FT algorithm using only two fast Fourier transforms has the best CT performance.
A two-target scenario consisting of two point scatters placed at a distance of 3 m from the antenna plane was considered to better highlight the resolution properties of the considered algorithms. The cross-range distance between the two targets was varied up to reaching, in x = 0, a −3dB value for DAS and 2D-FT algorithms (see Figure 15). This is achieved when the distance between the two targets is 1.2 m. This value is close to the theoretical cross-range resolution of the considered radar system (δcr = 1 m at 3 m distance). Moreover, in this case, the super resolution properties of the MUSIC algorithm are also evidenced.
A further study was performed to evaluate the field of view (FOV) of the three considered techniques. To this end, the point target was moved along a line parallel to the antenna array. This means that different θ angles with respect to the antenna center identify the target (see Figure 13). Figure 16 shows the percentage error on the estimated distance with respect to the real one obtained for the three considered algorithms.
The figure shows that if the considered percentage error is 5%, the FOV is equal to ±16°, ±21°, and ±30° for MUSIC, 2FFT, and DAS algorithms, respectively.
Finally, it is interesting to note that the theoretical unambiguous angular section of the radar for d a = λ 0 2 results in ϑ M A X = ± λ 0 2 d a = ± 28.6 ° , which is very close to the simulated DAS FOV.

4.3. Target Identification with Measured Data

The realized modular radar was used to study a target location in a realistic environment. The considered scenario investigated consists of a 60 × 60 cm2 metallic panel placed at a distance (d) from the antenna plane (see Figure 17).
The signal coming from each antenna, sequentially selected with the help of an electronic switch (see Figure 3), is sent to the receiver, where a DAQ board, operating at a sampling frequency of 200 kHz, takes 208 samples over a chirp period of 1039 μs. Since the sawtooth signal used to generate the chirp introduces a consistent ripple level on the received signal, the first 28 samples are canceled. The same procedure is applied to the signals coming from the other three antennas, whose related data are stored in a 4 × 180 matrix.
Figure 18 and Figure 19 show the range and cross-range cuts of the intensity function, corresponding to targets placed 2.6 m and 5.4 m away from the antenna plane, respectively.
For the panel placed at d = 2.6 m, the achieved HPBW in range for DAS and 2D-FT was 55 cm (see Figure 18a), which is close to the theoretical resolution (60 cm), while for the HPBW in cross-range, values of 60 cm both for 2D-FT and for DAS were achieved. Both these latter values are close to the theoretical resolution of 0.33 × 260 = 85 cm. Finally, the MUSIC results outline a range and cross-range HPBW of about 10 cm.
The same result of 55 cm was obtained for the range HPBW when the metal panel was placed at d = 5.4 m, while values of 130 cm for 2D-FT and for DAS of cross-range HPBW were achieved, at which a theoretical resolution of 0.33 × 540 = 178 cm was related. In the same condition, the MUSIC algorithm gives rise to a range and cross-range HPBW of about 6 cm and 20 cm, respectively.
The achieved results for DAS and 2D-FT confirm that the range HPBW, being related to the signal bandwidth, does not change with the distance, while the cross-range HPBW, depending on the array–target distance, varies for the two considered cases.
To evaluate the effectiveness of the radar spatial discrimination, the experimental scenario consisting of two metal panels, as shown in Figure 20, was considered. Figure 21 shows the map of the intensity function achieved with the three considered reconstruction algorithms: DAS (a), 2D-FT (b), and MUSIC (c).
All the considered algorithms were able to identify the two metal targets with errors less than 20 cm. The hot spot produced by the closest panel in (−0.7 m, 2.8 m) is less intense due to the smaller panel size and lateral displacement, both of which are responsible for a lower back-scattered field toward the radar antennas with respect to the furthest panel. In particular, the radar cross-section of a square metal panel can be evaluated as follows [38]:
R C S = 4 π a 2 λ 2 s i n k a s i n θ k a s i n θ
where “a” is the panel side, k = 2 π λ is the free-space wavenumber, and θ is the angle with respect to the panel normal. With reference to the square panels used in the experimental scenario reported in Figure 20, and taking into account the panel sizes and positions, the RCS seen by the radar is 19 m2 and 950 m2 for Panel 1 and Panel 2, respectively.

5. Conclusions

An FMCW radar, working at 24 GHz for short-range target imaging, realized by assembling the BGT24MTR11, ADF4159, and HMC1084LC4 Eval boards, and using a simple serial patch antenna array, was presented. A new, fast, and simple-to-implement 2D-FT algorithm was developed and tested on canonical scenarios. The proposed algorithm showed that a condition based on the maximum average error over the entire ROI is suitable to characterize its performance. The proposed algorithm was compared with the DAS and MUSIC algorithms. In particular, performances similar to those of the DAS were obtained, while lower resolution but better field of view were observed with respect to the results obtainable using MUSIC. Moreover, the 2D-FT algorithm is about 5 times faster than the DAS and more than 20 times faster than the MUSIC one. Another advantage of the radar proposed is that the number and the technology (printed, dielectric, and metallic) of antennas can be easily modified according to the design requirements. In conclusion the proposed radar system can be suitably employed for short-range real-time applications, such as the control of subjects and objects on the sidewalks of train and subway stations and the monitoring of elderly subjects in their home.

Author Contributions

Conceptualization, R.C., S.P., E.P. and O.T.; Methodology, R.C., S.P. and O.T.; Validation, R.C., S.P., E.P. and O.T.; Writing—original draft, R.C., S.P. and O.T.; Writing—review & editing, E.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Intensity Function Approximations

The second-order Taylor expansion of the phase term in (2) is as follows:
4 π f c x ξ 2 + y 2 4 π f c x 2 + y 2 4 π f c x ξ x 2 + y 2 + 2 π f c y 2 ξ 2 x 2 + y 2 3 / 2
Equation (A1) can be rewritten as follows,
4 π f c x ξ 2 + y 2 4 π f c x 2 + y 2 4 π f 0 c x ξ x 2 + y 2
if it is possible to neglect the quadratic term in ξ, and if the radar system fractional bandwidth is narrow. The validity of these assumptions is described below.
  • Condition for neglecting the quadratic term in (3)
    (a)
    Very stringent condition based on the maximum error over the entire ROI:
    E Q M a x = M a x 4 π f c y 2 ξ 2 2 x 2 + y 2 3 / 2 f 0 B 2 f f 0 + B 2 L 2 ξ L 2 x m i n x x m a x y m i n y y m a x 1
The maximum is reached for the following:
f = f m a x ,   x = 0 ,   y = y m i n ,   ξ = L 2
This then gives the following:
E Q M a x = π f m a x c L 2 2 y m i n = π 2 L 2 λ m i n y m i n 1
(b)
Less stringent condition based on the maximum average error over the entire ROI:
E Q M e a n M a x = M a x 1 L B f 0 B 2 f 0 + B 2 L 2 + L 2 4 π f c y 2 ξ 2 2 x 2 + y 2 3 / 2 d ξ d f x m i n x x m a x y m i n y y m a x 1
The maximum is obtained for the following:
x = 0 ,   y = y m i n
This then gives the following:
E Q M e a n M a x = π f 0 L 2 6   c   y m i n = π 6 L 2 λ 0   y m i n 1
  • Conditions for approximating the linear term in (3)
    (a)
    Very stringent condition based on the maximum error over the entire ROI:
E L M a x = M a x 4 π f f 0 c x ξ x 2 + y 2 f 0 B 2 f f 0 + B 2 L 2 ξ L 2 x m i n x x m a x y m i n y y m a x 1
The maximum is reached for the following:
f f 0 = B 2 ,   y = y m i n ,   ξ = L 2
This then gives the following:
E L M a x = π B L λ 0 f 0 M a x x m i n x m i n 2 + y m i n 2 , x m a x x m a x 2 + y m i n 2 1
Then that is rearranged as the following:
B f 0 λ 0 π L M i n x m i n 2 + y m i n 2 x m i n , x m a x 2 + y m i n 2 x m a x
(b)
Less stringent condition based on the maximum average error over the entire ROI:
E L M e a n M a x = M a x 1 L B f 0 B 2 f 0 + B 2 L 2 + L 2 4 π f f 0 c x ξ x 2 + y 2 d ξ d f x m i n x x m a x y m i n y y m a x 1
The maximum is reached for the following:
y = y m i n
This then gives the following:
E L M e a n M a x = π B L 4 λ 0 f 0 M a x x m i n x m i n 2 + y m i n 2 , x m a x x m a x 2 + y m i n 2 1
Then this is rearranged as the following:
B f 0 4 λ 0 π L M i n x m i n 2 + y m i n 2 x m i n , x m a x 2 + y m i n 2 x m a x
It is interesting to note that these last two conditions express an upper limit to the fractional bandwidth that can be adopted by the radar system used to detect the target.

Appendix B. Double-Fast FT Technique Numerical Implementation

Equation (3) can be implemented efficiently by using the 2D-FFT. Let Γ be the matrix of dimensions N a (even) × N f (odd), whose coefficients are given by
Γ p q = Γ ξ p , f q
where
ξ p = N a 1 2 d a + d a p             p = 0 , , N a 1
f q = f m i n + f q                   q = 0 , , N f 1
and
f = B N f 1
The sampling step of the variables (ξ and f) reported in (A18)–(A19) guarantees the coverage of the domain (γ and τ) shown in Figure 2 if the following conditions are satisfied:
1 2 d a γ m i n 1 2 d a 1 2 N a γ m a x
1 2 f 1 1 N f τ m i n 1 2 f 1 1 N f τ m a x
Equations (A21) and (A22) are derived by using the domain-to-codomain relationship of the 2D-FFT transform.
By taking (∆x and ∆y) the spatial steps employed to sample the ROI with ( N x   a n d   N y ) samples, we obtain the following:
x = x m a x x m i n N x 1
y = y m a x y m i n N y 1
Then, under the hypothesis based on the approximation (13), the extension of zero-padding ( N a , N f ) related to the spatial steps (∆x and ∆y) is obtained as follows:
N a = 2 λ 0 y m a x 4 x d a χ x
N f = 2 c 4 y f χ y + 1
where χ x 1 and χ y 1 are, respectively, the subsampling factors useful to reduce the size of the matrices used in 2D-FFT operations when very small sampling steps are used (note that, for χ x = 1 and χ y = 1 , no subsampling is employed), while · is the floor operator.
Finally, by using (A17)–(A26), the 2D-FFT takes the following form:
Γ r s F F T p = 0 N a 1 q = 0 N f 1 Γ p q e j 2 π s p N a e j 2 π r q N f
where
Γ p q = Γ p q             0 p N a 1             0 q N f 1 0                                                 o t h e r w i s e
(Note that (A27) can be straightforwardly implemented by using the ifft · , N a , N f instruction defined in the MATLAB environment with N a and N f zero-padding.)
Rearranging the coefficients of the matrix, Γ r s F F T , in such a way as to be associated with the sampled transformed variables γ r , τ s results in the following:
γ r = 1 2 d a + r N a d a                       0 r N a 1
τ s = 1 N f 2 1 N f f + s N f f                   0 s N f 1
Then the following holds:
Γ r s F F T = Γ r + N a 2 s + N f + 1 2 F F T           0 r N a 2 1           0 s N f 1 2 1 Γ r + N a 2 s N f 1 2 F F T           0 r N a 2 1           N f 1 2 s N f 1 Γ r N a 2 s + N f + 1 2 F F T           N a 2 r N a 1           0 s N f 1 2 1 Γ r N a 2 s N f 1 2 F F T           N a 2 r N a 1           N f 1 2 s N f 1
(Note that (A31) can be straightforwardly implemented by using the f f t s h i f t · function defined in the MATLAB environment.)
Using (A29)–(A30), the limits of the index r that allow us to cover the interval γ m i n , γ m a x , as well those useful to identify the limits of the index s that allow us to cover the interval τ m i n , τ m a x , are given by the following:
r m i n = γ m i n + 1 2 d a N a d a
r m a x = γ m a x + 1 2 d a N a d a
s m i n = τ m i n 1 N f 2 1 N f f N f f
s m a x = τ m a x 1 N f 2 1 N f f N f f
Finally, using (4), (A29), and (A30), the 2D-FT intensity function at the ROI sampling points x u , y v is
x u = x m i n + x   u             u = 0 , , N x
y v = y m i n + y   v             v = 0 , , N y
and is evaluated by means of the bilinear interpolator procedure reported in (A38).
F 2 F T x u , y v = = N a d a N f f 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p N a d a 2 c x u 2 + y u 2 1 N f 2 1 N f f q N f f Γ p + 1 q + 1 F F T + 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p N a d a 2 c x u 2 + y u 2 1 N f 2 1 N f f q + 1 N f f Γ p + 1 q F F T + 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p + 1 N a d a 2 c x u 2 + y u 2 1 N f 2 1 N f f q N f f Γ p q + 1 F F T + + 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p + 1 N a d a 2 c x u 2 + y u 2 1 N f 2 1 N f f q + 1 N f f Γ p q F F T r m i n p r m a x 1 s m i n q s m a x 1 N a d a 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p N a d a Γ p + 1 q F F T + 2 f 0 c x u x u 2 + y u 2 + 1 2 d a p + 1 N a d a Γ p q F F T r m i n p r m a x 1 q = s m a x N f f 2 c x u 2 + y u 2 1 N f 2 1 N f f q N f f Γ p q + 1 F F T + 2 c x u 2 + y u 2 1 N f 2 1 N f f q + 1 N f f Γ p q F F T p = r m a x s m i n q s m a x 1 Γ p q F F T p = r m a x q = s m a x 0 otherwise
where
p = 2 f 0 c   x u x u 2 + y u 2 + 1 2   d a N a   d a
q = 2 c x u 2 + y u 2 1 N f 2   1 N f   Δ f N f   Δ f
Equation (A38) can be straightforwardly implemented by using the i n t e r p 2 · , l i n e a r function available in the MATLAB environment. Doing so results in the following:
I 2 D F T x u , y v = F 2 D F T x u , y v F 2 D F T M a x
where F 2 D F T M a x is the maximum magnitude value taken by F 2 D F T x u , y v in the ROI.

References

  1. Pisa, S.; Pittella, E.; Piuzzi, E. A survey of radar systems for medical applications. IEEE Aerosp. Electron. Syst. Mag. 2016, 31, 64–81. [Google Scholar] [CrossRef]
  2. Kim, S.; Nguyen, C. On the development of a multifunction millimeter-wave sensor for displacement sensing and low-velocity measurements. IEEE Trans Microw. Theory Tech. 2004, 52, 2503–2512. [Google Scholar] [CrossRef]
  3. Li, C.; Cummings, J.; Lam, J.; Graves, E.; Wu, W. Radar remote monitoring of vital signs. IEEE Microw. Mag. 2009, 49, 47–56. [Google Scholar] [CrossRef]
  4. McEwan, T.E. Ultra Wideband Radar Motion Sensor. U.S. Patent 5,361,070, 1 November 1994. [Google Scholar]
  5. Pisa, S.; Bernardi, P.; Cavagnaro, M.; Pittella, E.; Piuzzi, E. A circuit model of an ultra-wideband impulse radar system for breath-activity monitoring. Int. J. Numer. Model. 2012, 25, 46–63. [Google Scholar] [CrossRef]
  6. Wang, G.; Gu, C.; Inoue, T.; Li, C. A hybrid FMCW-interferometry radar for indoor precise positioning and versatile life activity monitoring. IEEE Trans. Microw. Theory Tech. 2014, 62, 2812–2822. [Google Scholar] [CrossRef]
  7. Sacco, G.; Piuzzi, E.; Pittella, E.; Pisa, S. An FMCW Radar for localization and vital signs measurement for different chest orientations. Sensors 2020, 20, 3489. [Google Scholar] [CrossRef]
  8. Huang, Y.; Brennan, P.V.; Patrick, D.; Weller, I.; Roberts, P.; Hughes, K. FMCW based MIMO imaging radar for maritime navigation. Prog. Electromagn. Res. 2011, 115, 327–342. [Google Scholar] [CrossRef]
  9. Jeong, S.H.; Yu, H.Y.; Lee, J.E.; Oh, J.N.; Lee, K.H. A multi-bean and multi-range Radar with FMCW and digital beam forming for automotive applications. Prog. Electromagn. Res. 2012, 124, 285–299. [Google Scholar] [CrossRef]
  10. Hyun, E.; Jin, Y.S.; Lee, J.H. Design and implementation of 24 GHz multichannel FMCW surveillance Radar with a software-reconfigurable baseband. J. Sens. 2017, 2017, 3148237. [Google Scholar] [CrossRef]
  11. Rizzi Varela, V.G.; Rodrigues, D.V.Q.; Zeng, L.; Li, C. Multitarget Physical Activities Monitoring and Classification Using a V-Band FMCW Radar. IEEE Trans. Instrum. Meas. 2023, 72, 8500910. [Google Scholar] [CrossRef]
  12. Jiang, L.; Wu, M.; Che, L.; Xu, X.; Mu, Y.; Wu, Y. Continuous Human Motion Recognition Based on FMCW Radar and Transformer. J. Sens. 2023, 2023, 2951812. [Google Scholar] [CrossRef]
  13. Smith, J.W.; Thiagarajan, S.; Willis, R.; Makris, Y.; Torlak, M. Improved Static Hand Gesture Classification on Deep Convolutional Neural Networks Using Novel Sterile Training Technique. IEEE Access 2021, 9, 10893–10902. [Google Scholar] [CrossRef]
  14. Choi, J.-W.; Park, C.-W.; Kim, J.-H. FMCW Radar-Based Real-Time Hand Gesture Recognition System Capable of Out-of-Distribution Detection. IEEE Access 2022, 10, 87425–87434. [Google Scholar] [CrossRef]
  15. Mercuri, M.; Sacco, G.; Hornung, R.; Zhang, P.; Visser, H.; Hijdra, M.; Pisa, S.; Liu, Y.-H.; van Liempd, B.; Torfs, T. 2-D localization, angular separation and vital signs monitoring using a SISO FMCW Radar for smart long-term health monitoring environments. IEEE Internet Things J. 2021, 8, 11065–11077. [Google Scholar] [CrossRef]
  16. Li, S.; Ding, J.; Liu, W.; Li, H.; Zhou, F.; Zhu, Z. Railway Track Recognition Based on Radar Cross-Section Statistical Characterization Using mmWave Radar. Remote Sens. 2022, 14, 294. [Google Scholar] [CrossRef]
  17. Infineon Technologies AG. AN543: Distance2Go App; Note BGT24MTR11; Infineon Technologies AG: Neubiberg, Germany, 2019. [Google Scholar]
  18. Cardillo, E.; Li, C.; Caddemi, A. Vital sign detection and Radar self-motion cancellation through clutter identification. IEEE Trans. Microw. Theory Tech. 2021, 69, 1932–1942. [Google Scholar] [CrossRef]
  19. Tavanti, E.; Rizik, A.; Fedeli, A.; Caviglia, D.D.; Randazzo, A. A Short-Range FMCW Radar-Based Approach for Multi-Target Human-Vehicle Detection. IEEE Trans. Geosci. Remote Sens. 2022, 60, 2003816. [Google Scholar] [CrossRef]
  20. Analog Devices UG-1709: Evaluating the TinyRad 24 GHz Demonstration Platform. 2020. Available online: https://www.analog.com/media/en/technical-documentation/user-guides/ev-tinyrad24g-ug-1709.pdf (accessed on 15 April 2023).
  21. IMST GmbH DK-sR-1200e: 24 GHz FMCW-Radar Development Kit. 2017. Available online: https://imst.com/imst/en/support/download/circuits/LR_HZ_DK-sR-1200e_Radar_RZ_17_9_E.pdf (accessed on 15 April 2023).
  22. Pomerleau, P.; Royer, A.; Langlois, A.; Cliche, P.; Courtemanche, B.; Madore, J.B.; Picard, G.; Lefebvre, É. Low cost and compact FMCW 24 GHz Radar applications for snowpack and ice thickness measurements. Sensors 2020, 20, 3909. [Google Scholar] [CrossRef]
  23. Cicchetti, R.; Pisa, S.; Piuzzi, E.; Pittella, E.; D’Atanasio, P.; Testa, O. Numerical and experimental comparison among a new hybrid FT-Music technique and existing algorithms for through-the-wall radar imaging. IEEE Trans. Microw. Theory Tech. 2021, 69, 3372–3387. [Google Scholar] [CrossRef]
  24. Infineon Technologies, Data sheet BGT24MTR11: Silicon Germanium 24 GHz Transceiver MMIC. 2014. Available online: https://www.infineon.com/dgdl/Infineon-BGT24MTR11-DS-v03_01-EN.pdf (accessed on 15 April 2023).
  25. Analog Devices, User Guide, EV-ADF4159EB1Z/EV, UG-383: Frequency Synthesizer for Phase-Locked Loops. 2013. Available online: https://www.analog.com/media/en/technical-documentation/evaluation-documentation/UG-383.pdf (accessed on 15 April 2023).
  26. Analog Devices, Datasheet HMC1084LC4, GaAs MMIC SP4T Reflective Switch 23–30 GHz. 2014. Available online: https://www.analog.com/media/en/technical-documentation/data-sheets/hmc1084.pdf (accessed on 15 April 2023).
  27. National Instrument, Datasheet DAQ NI-USB 6361. 2023. Available online: https://www.ni.com/docs/en-US/bundle/pcie-pxie-usb-6361-specs/page/specs.html (accessed on 15 April 2023).
  28. Caliskan, A.; Yamac, Y.E.; Uysal, H.; Turk, A.S.; Kizilay, A.; Orhan, M. 24 GHz short range radar system measurements for synthetic aperture radar imaging. In Proceedings of the IEEE Microwaves, Radar and Remote Sensing Symposium (MRRS), Kiev, Ukraine, 29–31 August 2017; pp. 285–287. [Google Scholar]
  29. Ali, F.; Urban, A.; Vossiek, M. A Short Range Synthetic Aperture Imaging Radar with Rotating Antenna. Int. J. Electron. Telecommun. 2011, 57, 97–102. [Google Scholar] [CrossRef]
  30. Baldazzi, E.; Al-Rawi, A.; Cicchetti, R.; Smolders, A.B.; Testa, O.; van Coevorden Moreno, C.D.J.; Caratelli, D. A high-gain dielectric resonator antenna with plastic-based conical horn for millimeter-wave applications. IEEE Antennas Wirel. Propag. Lett. 2020, 19, 949–953. [Google Scholar] [CrossRef]
  31. Caratelli, D.; Cicchetti, R.; Cicchetti, V.; Testa, O. A wideband high-gain circularly-polarized dielectric horn antenna equipped with Lamé-axicon stacked-disk lens for remote sensing, air traffic control and satellite communications. IEEE Access 2023, 11, 20912–20922. [Google Scholar] [CrossRef]
  32. Bourqui, J.; Okoniewski, M.; Fear, E.C. Balanced antipodal Vivaldi antenna for breast cancer detection. In Proceedings of the Second European Conference on Antennas and Propagation, Edinburgh, UK, 11–16 November 2007; pp. 1–5. [Google Scholar]
  33. Cicchetti, R.; Cicchetti, V.; Faraone, A.; Foged, L.; Testa, O. A compact high-gain wideband lens Vivaldi antenna for wireless communications and through-the-wall imaging. IEEE Trans. Antennas Propag. 2021, 69, 3177–3192. [Google Scholar] [CrossRef]
  34. Pisa, S.; Pittella, E.; Piuzzi, E. Serial patch array antenna for an FMCW radar housed in a white cane. Int. J. Antennas Propag. 2016, 2016, 9458609. [Google Scholar] [CrossRef]
  35. Yin, J.; Wu, Q.; Yu, C.; Wang, H.; Hong, W. Low-Sidelobe-Level Series-Fed Microstrip Antenna Array of Unequal Interelement Spacing. IEEE Antennas Wirel. Propag. Lett. 2017, 16, 1695–1698. [Google Scholar] [CrossRef]
  36. Yuan, T.; Yuan, N.; Li, L.-W. A Novel Series-Fed Taper Antenna Array Design. IEEE Antennas Wirel. Propag. Lett. 2008, 7, 362–365. [Google Scholar] [CrossRef]
  37. ADIsimPLL, version 5.10; Analog Devices, Inc.: Cambridge, MA, USA, 2019.
  38. Chen, W.-K. The Electrical Engineering Handbook; Elsevier Academic Press: San Diego, CA, USA, 2005. [Google Scholar]
Figure 1. Considered geometry. Antenna array (blue region) and ROI (blue dashed line) are shown in the figure.
Figure 1. Considered geometry. Antenna array (blue region) and ROI (blue dashed line) are shown in the figure.
Sensors 23 04119 g001
Figure 2. Mapping of the ROI in the corresponding region expressed by means of the γ and τ variables.
Figure 2. Mapping of the ROI in the corresponding region expressed by means of the γ and τ variables.
Sensors 23 04119 g002
Figure 3. Block diagram of the realized radar (a). Block diagram of the BGT24MTR11 (b) and the ADF4159 (c).
Figure 3. Block diagram of the realized radar (a). Block diagram of the BGT24MTR11 (b) and the ADF4159 (c).
Sensors 23 04119 g003
Figure 4. Schematic of the realized serial patch array antenna (a) and its radiation pattern (b). Theta = 0° corresponds to the normal at the air–dielectric interface of the array antenna.
Figure 4. Schematic of the realized serial patch array antenna (a) and its radiation pattern (b). Theta = 0° corresponds to the normal at the air–dielectric interface of the array antenna.
Sensors 23 04119 g004
Figure 5. Schematic of the TX/RX antenna array.
Figure 5. Schematic of the TX/RX antenna array.
Sensors 23 04119 g005
Figure 6. Measured (meas) and simulated (sim) scattering parameters of the array antenna system.
Figure 6. Measured (meas) and simulated (sim) scattering parameters of the array antenna system.
Sensors 23 04119 g006
Figure 7. Schematic of the PLL loop (a) and its realization with SMD technology (b).
Figure 7. Schematic of the PLL loop (a) and its realization with SMD technology (b).
Sensors 23 04119 g007
Figure 8. Measured and simulated responses of the loop filter.
Figure 8. Measured and simulated responses of the loop filter.
Sensors 23 04119 g008
Figure 9. Measured ramp signal driving the VCO; and VCO output frequency as a function of time.
Figure 9. Measured ramp signal driving the VCO; and VCO output frequency as a function of time.
Sensors 23 04119 g009
Figure 10. Prototype of the realized modular radar.
Figure 10. Prototype of the realized modular radar.
Sensors 23 04119 g010
Figure 11. Geometry of the scenario: in the x–y plane. Point target and considered ROI (dashed blue line).
Figure 11. Geometry of the scenario: in the x–y plane. Point target and considered ROI (dashed blue line).
Sensors 23 04119 g011
Figure 12. Cross-range cuts of the normalized intensity function for the 4 antennas and 16 antenna cases: (a) f0 = 2 GHz, FBW = 0.1; (b) f0 = 24 GHz, FBW = 0.1.
Figure 12. Cross-range cuts of the normalized intensity function for the 4 antennas and 16 antenna cases: (a) f0 = 2 GHz, FBW = 0.1; (b) f0 = 24 GHz, FBW = 0.1.
Sensors 23 04119 g012
Figure 13. Geometry of the scenario in the x–y plane. Point target 3 m far from the antenna plane and considered ROI (dashed blue line).
Figure 13. Geometry of the scenario in the x–y plane. Point target 3 m far from the antenna plane and considered ROI (dashed blue line).
Sensors 23 04119 g013
Figure 14. Comparison between the three considered techniques in the presence of a single target. The normalized intensity function versus (a) range and (b) cross-range distance is reported.
Figure 14. Comparison between the three considered techniques in the presence of a single target. The normalized intensity function versus (a) range and (b) cross-range distance is reported.
Sensors 23 04119 g014
Figure 15. Comparison between the three considered techniques in the presence of two targets. The normalized intensity function versus cross-range distance is reported.
Figure 15. Comparison between the three considered techniques in the presence of two targets. The normalized intensity function versus cross-range distance is reported.
Sensors 23 04119 g015
Figure 16. Target position error versus identification ϑ angle shown in Figure 13 for the three considered techniques.
Figure 16. Target position error versus identification ϑ angle shown in Figure 13 for the three considered techniques.
Sensors 23 04119 g016
Figure 17. Sketch of the experimental investigated scenario.
Figure 17. Sketch of the experimental investigated scenario.
Sensors 23 04119 g017
Figure 18. The normalized intensity function versus (a) range and (b) cross-range distance is reported. Panel distance is 2.6 m away from the radar antenna.
Figure 18. The normalized intensity function versus (a) range and (b) cross-range distance is reported. Panel distance is 2.6 m away from the radar antenna.
Sensors 23 04119 g018
Figure 19. The normalized intensity function versus (a) range and (b) cross-range distance is reported. Panel distance is 5.4 m away from the radar antenna.
Figure 19. The normalized intensity function versus (a) range and (b) cross-range distance is reported. Panel distance is 5.4 m away from the radar antenna.
Sensors 23 04119 g019
Figure 20. Sketch of the experimental investigated scenario with two square metal panel targets: metal panel 1 (30 × 30 cm) and metal panel 2 (60 × 60 cm).
Figure 20. Sketch of the experimental investigated scenario with two square metal panel targets: metal panel 1 (30 × 30 cm) and metal panel 2 (60 × 60 cm).
Sensors 23 04119 g020
Figure 21. Map of the normalized intensity function in the ROI for the scenario in Figure 20 obtained with DAS (a), 2D-FT (b), and MUSIC (c) algorithms.
Figure 21. Map of the normalized intensity function in the ROI for the scenario in Figure 20 obtained with DAS (a), 2D-FT (b), and MUSIC (c) algorithms.
Sensors 23 04119 g021
Table 1. Phase error produced by the approximation of the square term.
Table 1. Phase error produced by the approximation of the square term.
E Q M a x E Q M e a n M a x
f0 = 24.125 GHz, B = 250 MHz FBW = 0.01, Na = 40.0440.015
f0 = 24 GHz, B = 2.4 GHz FBW = 0.10, Na = 40.0460.015
f0 = 2 GHz, B = 0.2 GHz FBW = 0.10, Na = 40.560.18
f0 = 24.125 GHz, B = 250 MHz FBW = 0.01, Na = 161.110.37
f0 = 24 GHz, B = 2.4 GHz FBW = 0.10, Na = 161.160.37
f0 = 2 GHz, B = 0.2 GHz FBW = 0.10, Na = 1613.94.42
Table 2. Phase error produced by the approximation of the linear term.
Table 2. Phase error produced by the approximation of the linear term.
E L M a x E L M e a n M a x
f0 = 24.125 GHz, B = 250 MHz FBW = 0.01, Na = 40.044 B f 0 = 0.24 M a x 0.011 B f 0 = 0.95 M a x
f0 = 24 GHz, B = 2.4 GHz FBW = 0.10, Na = 40.42 B f 0 = 0.24 M a x 0.10 B f 0 = 0.95 M a x
f0 = 2 GHz, B = 0.2 GHz FBW = 0.10, Na = 40.42 B f 0 = 0.24 M a x 0.10 B f 0 = 0.95 M a x
f0 = 24.125 GHz, B = 250 MHz FBW = 0.01, Na = 160.22 B f 0 = 0.048 M a x 0.055 B f 0 = 0.19 M a x
f0 = 24 GHz, B = 2.4 GHz FBW = 0.10, Na = 162.10 B f 0 = 0.048 M a x 0.53 B f 0 = 0.19 M a x
f0 = 2 GHz, B = 0.2 GHz FBW = 0.10, Na = 162.11 B f 0 = 0.048 M a x 0.53 B f 0 = 0.19 M a x
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cicchetti, R.; Pisa, S.; Piuzzi, E.; Testa, O. A Double Fourier-Transform Imaging Algorithm for a 24 GHz FMCW Short-Range Radar. Sensors 2023, 23, 4119. https://doi.org/10.3390/s23084119

AMA Style

Cicchetti R, Pisa S, Piuzzi E, Testa O. A Double Fourier-Transform Imaging Algorithm for a 24 GHz FMCW Short-Range Radar. Sensors. 2023; 23(8):4119. https://doi.org/10.3390/s23084119

Chicago/Turabian Style

Cicchetti, Renato, Stefano Pisa, Emanuele Piuzzi, and Orlandino Testa. 2023. "A Double Fourier-Transform Imaging Algorithm for a 24 GHz FMCW Short-Range Radar" Sensors 23, no. 8: 4119. https://doi.org/10.3390/s23084119

APA Style

Cicchetti, R., Pisa, S., Piuzzi, E., & Testa, O. (2023). A Double Fourier-Transform Imaging Algorithm for a 24 GHz FMCW Short-Range Radar. Sensors, 23(8), 4119. https://doi.org/10.3390/s23084119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop