Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

High-resolution imaging of a high-speed target based on a reconfigurable photonic fractional Fourier transformer

Open Access Open Access

Abstract

The previously reported photonics-based radar working with a large bandwidth has the advantages of realizing high-resolution imaging of targets with low velocity. However, the high velocity of a target will introduce Doppler dispersion to the echo signals, which severely deteriorates the imaging resolution. This problem becomes more noticeable as the bandwidth increases. In this paper, we propose a radar receiver based on a reconfigurable photonic fractional Fourier transformer (PFrFTer). The order of the PFrFTer can be reconstructed flexibly by changing the optical transform kernel. When the transform order matches the velocity of the target, the chirp echo signals behave as narrow impulses in the fractional Fourier domain, showing the range information with a high resolution. In the experiment, a PFrFTer is established and applied to process the echo signals with a bandwidth of 12 GHz. A lossless range resolution of 1.4 cm is obtained in range profiles and inverse synthetic aperture radar imaging for high-speed targets. This range resolution is much higher than that in the classical optical de-chirping receiver. These results demonstrate the PFrFTer is immune to the Doppler dispersion effect and is excellent for high-resolution imaging of high-speed target. The introduced technique would be of practical interest in the detection and recognition of targets.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Imaging radar plays an essential role in the detection and recognition of targets. It has many applications such as security check [1], foreign-object debris detection in airport runways [2], satellite monitoring [3] and so on [4]. In order to obtain more accurate identification, the imaging radar needs to have higher resolution capability, which demands the radar to work at a larger bandwidth and a higher frequency. However, due to the electronic bottlenecks, it is difficult for the traditional radar to generate and process such signals [5].

Fortunately, microwave photonics technology offers a promising solution to solve this issue for its unique characters of high frequency, wide bandwidth, immunity to electromagnetic interference, etc [6,7]. Due to these merits, many relevant techniques are proposed and demonstrated to generate, distribute and process radio frequency (RF) signals. As for the generation of RF signals, there are several methods including photonic digital-to-analog converter [8], frequency to time mapping [9], and photonic frequency multiplication [10]. As for dissemination of ultra-stable RF signal over a large area, passive compensation scheme and active compensation scheme have been presented, such as photonic microwave phase conjugation, phase-locked loops [11,12]. In terms of processing RF signals, techniques such as photonic filter [13,14], photonic phase shifter [15], photonic analog-to-digital converter (PADC) [16] and photonic de-chirping [17] have been developed extensively. In recent years, the development of these technologies has led to the emergence of high-performance microwave photonic radar systems, including photonic-based inverse synthetic aperture radar (ISAR) [1721], photonics-based multiple-input and multiple-output radar [2224], and photonics-based distributed radar [25,26]. Most of these radar systems have high–resolution capability. For example, previously, we proposed a W-band photonics-based radar with a bandwidth of 8 GHz [19]. A toy gun with a width of 18 cm is imaged and recognized by this radar. In [20], a small electric fan is imaged with a frame rate of 100 fps by a K-band photonics-based radar. In [21], a small unmanned aerial vehicle with an arm width of 1.9 cm is imaged and the detail structures of legs and arms are distinguished by a Ka-band photonics-based radar. From these research results, it is obvious that microwave photonics technology has taken radar imaging to the next level.

However, until now, the high-resolution imaging results taken by the photonics-based radar with a large bandwidth are obtained from stationary or low-speed targets. For high-speed targets, Doppler dispersion distorts the echo signals severely and deteriorates the imaging resolution, especially when the Doppler dispersion products $2BTv/c$ is larger than 1, where $v$, $B$, $T$ are the velocity, the bandwidth and the pulse width, respectively [27,28]. For example, in [29], when the velocity is 300 m/s, the spectrum of the de-chirped signal is expanded, resulting in a decrease in resolution. Although some compensation algorithms are proposed to solve this problem, they need to be performed in the digital signal processer with complex computation and take some time, which slow down the imaging rate, especially in the scenario where the imaging accumulation time is long and the amount of data is large [30,31]. Thus, the photonics-based radar mentioned before can work with a large bandwidth, but the imaging resolution would deteriorate if the target moves rapidly.

Since Doppler dispersion distorts the echo signals both in the time domain and frequency domain, it is invalid to suppress the Doppler dispersion effect by receiving and processing the echo signals only through methods in one of the domains, such as PADC, photonic de-chirping. Further, using these methods to process the echo signals by waveform design, such as dual-chirp LFM signals, is also invalid because the Doppler dispersion effect will act simultaneously on each chirp LFM signals, instead of cancelling each other out. As a two-dimensional signal processing method in the time-frequency domain, the fractional Fourier transform (FrFT) is an attractive method to suppress the Doppler dispersion effect. FrFT can be interpreted as a rotation in the time-frequency plane. It is a useful tool to analyze signals with the time-dependent spectrum, such as linear frequency modulated wave (LFMW) [32]. For a given LFMW, it behaves as a slant edge in the time-frequency plane. A proper rotation will transform the LFMW to an impulse signal in the fractional Fourier domain (FrFD). The Doppler dispersion effect can only change the chirp rate and center frequency of an LFMW, which still behaves as a slant edge in the time-frequency plane. Therefore, the received LFMW from high-speed targets can still concentrate to narrow pulses in the FrFD when the order of the FrFT matches the velocity, showing the high-resolution range information.

Some methods for realizing photonic fractional Fourier transformer (PFrFTer) have been proposed, and they can be used to process optical signal [33,34] or RF signal [3538]. The method capable of processing RF signal is mainly summarized into two categories. In the first category, the FrFT of a given RF signal is obtained by coherent addition of temporally delayed replicas of the input signal after modulation by a quadratic phase term [35]. Because of the large time delay, the bandwidth of the signal is limited to tens of MHz in the implementation. Peng et al. proposed an approach to extend the bandwidth and verified the effectiveness by numerical simulation, while the architecture is complex and hard to be implemented [36]. In the second category, the FrFT is mainly realized based on the product of the input RF signal and an LFMW [37,38]. The approach can be used to process signals with large bandwidth and the architecture is simple. However, it lacks tunability.

In this paper, we realize high-resolution imaging of high-speed target based on a reconfigurable PFrFTer. A linear chirp light serves as the transform kernel to rotate the time-frequency plane, resulting in the FrFD. The LFMWs received from high-speed targets are projected onto it and behave as narrow impulses when the order of the PFrFTer matches the velocity of targets. The order can be reconstructed by changing the chirp rate of the light, which can be practically achieved by a dispersion module (DM). Besides, the architecture of the PFrFTer is simple and it is able to process radar signals with broad bandwidth up to tens of GHz. Experiments are performed to verify the effectiveness of the PFrFTer. The transmitted signal is an LFMW with a large bandwidth of 12 GHz. The imaging of moving targets in different scenarios is obtained from the PFrFTer receiver and the classical optical de-chirping receiver. Comparing the two sets of results, it can be seen that the imaging obtained from the PFrFTer receiver is clearer than that from the classical optical de-chirping receiver. The imaging resolution from the PFrFTer receiver keeps lossless as 1.4 cm for high-speed targets, which is almost equal to the theoretical value of 1.3 cm. This technique would be of interest for the detection and recognition of targets.

2. Principle

2.1 FrFT of the echo signals from moving targets

The mathematical FrFT expression of a given signal $x(t )$ is written as [32]

$${F_p}\{{x(t )} \}(u )= \int_{ - \infty }^{ + \infty } {x(t ){K_\alpha }({u,t} )} dt\;\;\;\;\alpha \ne m\pi \;({m \in Z} ),$$
where ${K_\alpha }({u,t} )$ is the kernel, which is given by
$${K_\alpha }({u,t} )= \sqrt {1 - j\cot \alpha } \cdot \exp\left[ {j2\pi \left( { - u\csc \alpha t + \frac{{\cot \alpha }}{2}{t^2}\textrm{ + }\frac{{{u^2}\cot \alpha }}{2}} \right)} \right],$$
where u is the fractional frequency, $\alpha$ is the transform angle and $p = 2\alpha /\pi$ is the order. The FrFT can be interpreted as a rotation in the time-frequency plane with the angle $\alpha$, which is decided by the kernel. The square modulus of the FrFT is equal to the projection onto the u axis of the time-frequency distribution of the given signal, seen as Fig. 1(a). This characteristic is particularly valuable in processing the LFMW, whose time-frequency distribution is a slant line. At a proper order, the LFMW will concentrate maximally and behave as an impulse in the FrFD.

 figure: Fig. 1.

Fig. 1. The FrFT in the time-frequency domain. (a) Square modulus of the FrFT of a given signal at a specific order $2\alpha /\pi$. (b) Square modulus of the FrFT of the echo signals from a stationary target and a moving target at proper order $2{\alpha _1}/\pi$ and $2{\alpha _2}/\pi$, respectively.

Download Full Size | PDF

In a radar system, the transmitted signal is an LFMW, the pulse width, the center frequency and the chirp rate of which are T, ${f_1}$ and ${k_1}$, respectively. Assuming the distance and velocity of a target are R and v, the time delay can be expressed as $\tau (t )= \frac{{2({R\textrm{ - }vt} )}}{{c\textrm{ - }v}} \approx \frac{{2R}}{c} - \frac{{2v}}{c}t$, where c is the velocity of the light [39]. Then the received signal of a moving target can be expressed as

$$\begin{aligned} {s_r}(t )&\textrm{ = }rect\left( {\frac{{t - \tau (t )}}{T}} \right)\exp [{j({2\pi {f_1}({t - \tau (t )} )+ \pi {k_1}{{({t - \tau (t )} )}^2}} )} ]\\ &\approx rect\left( {\frac{{t - 2R/c}}{T}} \right)\exp [{j({2\pi {f_2}({t - {\tau_2}} )+ \pi {k_2}{{({t - {\tau_2}} )}^2}} )} ], \end{aligned}$$
where ${f_2} = \left( {1 + \frac{{2v}}{c}} \right){f_1}$, ${k_2} = {\left( {1 + \frac{{2v}}{c}} \right)^2}{k_1}$, ${\tau _2} = \frac{{2R}}{{c + 2v}}$. From Eq. (3), it can be seen that the chirp rate of the echo signal will change with the velocity. The FrFT of the echo signal can be written as
$${F_p}\{{{s_r}(t )} \}(u )\textrm{ = }\int_{ - \infty }^{ + \infty } {{s_r}(t ){K_\alpha }({u,t} )} dt.$$

When the order does not satisfy the formula $p = 2arc\cot \left( { - {{\left( {1 + \frac{{2v}}{c}} \right)}^2}{k_1}} \right)/\pi$, the echo signals will concentrate and behave as broadened impulses in the FrFD according to Eq. (4). Further, when $|{{k_2} + \cot \alpha } |{T^2} \gg 1$, Eq. (4) can be rewritten as

$${F_p}\{{{s_r}(t )} \}(u )\propto rect\left( {\frac{{u\csc \alpha \textrm{ + }{k_2}{\tau_2} - {f_2}}}{{({{k_2} + \cot \alpha } )T}}} \right).$$

The distance is derived as

$$R = ({c + 2v} )({{f_2} - u\csc \alpha } )/({2{k_2}} ).$$

According to Eq. (5) and Eq. (6), the resolution of $u\csc \alpha$ is $|{{k_2} + \cot \alpha } |T$, so the range resolution is approximately calculated as $|{{k_2} + \cot \alpha } |{T^2}c/({2B} )$. The resolution becomes worse as the absolute value of the sum of ${k_2}$ and $\cot \alpha$ is greater. Therefore, in order to get the highest resolution of the imaging, the angle should satisfy $\cot \alpha ={-} {k_2}$, meaning the order should be $p = 2arc\cot \left( { - {{\left( {1 + \frac{{2v}}{c}} \right)}^2}{k_1}} \right)/\pi$. Then the FrFT of the echo signal is calculated as

$${F_p}\{{{s_r}(t )} \}(u )\propto \textrm{sinc} [{T({u\csc \alpha \textrm{ + }{k_2}{\tau_2} - {f_2}} )} ].$$

It can be seen that the echo signals from stationary or moving targets are transformed into sinc-type pulses in the FrFD, as shown in Fig. 1(b). The distance is also derived as Eq. (6). According to Eq. (6) and Eq. (7), the resolution of $u\csc \alpha$ becomes $1/T$, so the range resolution is approximately calculated as , which is determined by the bandwidth of the transmitted signal. It is obvious that the resolution is lossless for targets with any velocity and is equal to the ideal resolution for an LFMW with a bandwidth of $B$[40].

2.2 Model of the proposed radar receiver based on the PFrFTer

Since the target's velocity is varying, the order of the FrFT must be reconstructed to match the velocity for high-resolution imaging. Here, a reconfigurable PFrFTer is proposed to realize it. The schematic diagram of the radar receiver based on the PFrFTer is shown in Fig. 2. The PFrFTer mainly consists of four components: light source, electro-optical conversion, opto-electronic conversion and Fourier transform. The light source includes an optical carrier and a coherent linear chirp light. Its optical field can be expressed as

$$\begin{aligned} {E_1} &\propto {E_0} + {E_0}s(t )\\ &\textrm{ = }\underbrace{{{E_0}}}_{{optical\;\;carrier}}\textrm{ + }\underbrace{{{\beta _1}{E_0}\exp [{j({2\pi {f_s}t + \pi {k_s}{t^2}\textrm{ + }{\varphi_s}} )} ]}}_{{coherent\;linear\;chirp\;light}}, \end{aligned}$$
where ${\beta _1}$ is the intensity ratio of the coherent linear chirp light to the optical carrier, ${f_s}$ is the difference of the center frequency between the optical carrier and the coherent linear chirp light, ${k_s}$ and ${\varphi _s}$ are the chirp rate and the initial phase of the linear chirp light, respectively. This light source can be generated by many methods. Here, it is generated by a laser diode (LD), single-sideband (SSB) modulator and a DM. The SSB modulator is driven by the reference signal, which is a delayed copy of transmitted signal and is written as ${s_{ref}}(t )= exp [{j({2\pi {f_1}t + \pi {k_1}{t^2}} )} ]$. In single lower sideband modulation, the output signal of the DM is the light source as expressed in Eq. (8) and ${f_s}$, ${k_s}$ and ${\varphi _s}$ are calculated as follows [41]
$$\begin{array}{l} {f_s} = c{f_1}/({c\textrm{ + }{k_1}D{\lambda^2}} )\\ {k_s} = c{k_1}/({c\textrm{ + }{k_1}D{\lambda^2}} )\\ {\varphi _s} = \pi f_1^2D{\lambda ^2}/({c\textrm{ + }{k_1}D{\lambda^2}} ), \end{array}$$
where $D ={-} 2\pi c{\beta _2}L/{\lambda ^2}$ (ps/nm) is the total dispersion value of the DM, ${\beta _2}$ is the dispersion coefficient, and $\lambda$ is the wavelength of the optical carrier. The coherent linear chirp light ${E_0}s(t )$ serves as the optical kernel of the PFrFTer. The optical carrier ${E_0}$ provides a way for the received signal ${s_r}(t )$ to turn into the optical domain by electro-optical conversion, resulting in a new component in the light ${E_0}{s_r}(t )$. Through opto-electronic conversion, this component multiplies the kernel, resulting in a mixing signal. After performing Fourier transform (FT) to the mixing signal, the echo signal can be projected to the FrFD, transforming to a narrowband signal, which is
$$i(f )= \int_{ - \infty }^{ + \infty } {{s_r}(t )\exp[{ - j({2\pi {f_s}t + \pi {k_s}{t^2}\textrm{ + }{\varphi_s}} )} ]} \cdot \exp ({ - j2\pi ft} )dt.$$

Making a variable substitution: ${k_s}\textrm{ = } - \cot \alpha ,\;\;f = u\csc \alpha - {f_s}$, Eq. (10) can be written as

$$\begin{aligned} |{i(u )} |&\propto \left|{\int_{ - \infty }^{ + \infty } {{s_r}(t )\exp\left[ {j2\pi \left( { - u\csc \alpha t + \frac{{\cot \alpha }}{2}{t^2}\textrm{ + }\frac{{{u^2}\cot \alpha }}{2}} \right)} \right]} } \right|dt\\ &\propto |{{F^p}[{{s_r}(t )} ](u )} |. \end{aligned}$$

 figure: Fig. 2.

Fig. 2. Schematic diagram of the radar receiver based on the PFrFTer. Ant: antenna; LD: laser diode; SSB: single-sideband; DM: dispersion module; E/O: electro-optical conversion; O/E: opto-electronic conversion. FT: Fourier transform.

Download Full Size | PDF

The order of the PFrFT is derived as

$$\begin{aligned} p &= 2arc\cot ({ - {k_s}} )/\pi \\ &= 2arc\cot \left( { - \frac{c}{{c\textrm{ + }{k_1}D{\lambda^2}}}{k_1}} \right)/\pi . \end{aligned}$$

It can be seen that the order can be reconstructed by changing the chirp rate of the linear chirp light. Specifically, it can be reconstructed by setting the dispersion. According to the discussion in section 2.1, a high-resolution imaging would be obtained when we adjust the order and make it equal to $2arc\cot ({ - {k_2}} )/\pi$. Besides, according to Eq. (6), the distance of the target can be rewritten as

$$R ={-} ({c + 2v} )({f + {f_s} - {f_2}} )/({2{k_2}} ).$$

In order to make the order equal to $2arc\cot ({ - {k_2}} )/\pi$, we only need to set the dispersion to

$$D ={-} \frac{{c({{{({c + 2v} )}^2} - {c^2}} )}}{{{k_1}{\lambda ^2}{{({c + 2v} )}^2}}} \approx{-} \frac{{4v}}{{{k_1}{\lambda ^2}}}.$$

It is worth noting that the previous analysis is also right for the single upper sideband modulation. In this case, the dispersion is needed to be set as

$$D = \frac{{c({{{({c + 2v} )}^2} - {c^2}} )}}{{{k_1}{\lambda ^2}{{({c + 2v} )}^2}}} \approx \frac{{4v}}{{{k_1}{\lambda ^2}}}.$$

Here, we can see that the dispersion with the same value can be used to process the targets with the same absolute velocity in opposite direction by single lower sideband modulation or single upper sideband modulation.

3. Experiment and discussion

To investigate the performance of the PFrFTer in the radar system, an experiment is carried out based on the schematic diagram in Fig. 2. An arbitrary waveform generator (AWG) is performed as the transmitter. Besides, it is also used to generate the simulated echo signals according to Eq. (3). Most studies about the detection of moving targets are all based on such experimental approach [29]. Although distortion will occur in the echo signals due to the electromagnetic scattering characteristic of the targets, transport characteristics of the environment in actual applications, it can be eliminated to a certain extent by some compensation approaches [42]. Besides, it does not affect the effectiveness of the proposed method to suppress the Doppler dispersion effect and realize high-resolution imaging. Therefore, these influences are not considered in the echo signals generated from the AWG. A light centered at 1545 nm with a power of 16 dBm is generated from an LD (keysight 81606A). SSB modulation is implemented by a Mach-Zehnder modulator (MZM) (Fujitsu FTM7938E) and an optical band-pass filter (OBPF) (WaveShaper 1000s). The MZM is biased at the quadrature point, resulting in dual-sideband modulation, and the OBPF is used to filter out one of the sidebands. The DM is realized by the dispersion compensation fiber (DCF). The electro-optical conversion is completed by an MZM (Photline MXAN-LN-40), which is also biased at the quadrature point. The opto-electronic conversion is achieved by a photodiode with a bandwidth of 3 GHz. A digital storage oscilloscope (Keysight UXR0134A) is used to sample the signal generated from the electro-optical conversion and the FT is performed by a digital signal processor.

For better analyzing the effectiveness of the PFrFTer in the radar system, a controlled experiment of classical optical de-chirping processing is also carried out. Most of the high-resolution photonics imaging radars presented before are based on classical optical de-chirping processing technology [1921]. The detail introduction about the optical de-chirping receiver can be found in our previous work [25]. The transmitted signal is an LFMW centered at 14 GHz with 12 GHz bandwidth and 250 µs pulse width, and its spectrum is shown in Fig. 3(a). The OBPF is set to filter out the upper sideband, resulting in single lower sideband modulation. The spectrum of the light is shown in Fig. 3(b).

 figure: Fig. 3.

Fig. 3. (a) The spectrum of the transmitted signal. (b) The spectrum of the optical single-sideband signal.

Download Full Size | PDF

Four targets with different velocities are assumed and the echo signals from them are generated by the AWG. The four targets are all single scattering point targets. In other words, there is only one strong scattering point in each target. Their velocities are 0 m/s, 56.91m/s, 113.45m/s and 151.40 m/s, as shown in Table 1. The optical de-chirping process is performed to the echo signals, and the range profile results of the target are obtained. Figures 4(a)-(d) are the range profiles of the targets 1-4 respectively. As can be seen from the graphs, the −3 dB main-lobe width increases and the peak side-lobe ratio (PSLR) decreases with larger velocity, which means the range resolution deteriorates when the velocity gets larger. These results are caused by the Doppler dispersion effect. For a stationary target, the de-chirped signal is a single frequency signal with a narrow spectrum. However, for a moving target, the de-chirped signal is an LFMW thanks to the Doppler dispersion effect and its spectrum expands with the increase of the velocity, making the deterioration of the range resolution. The PFrFTer receiver is also used to detect targets 1-4. According to Eq. (12), the order of the PFrFTer is reconstructed by setting the dispersion of the DCF to 0 ps/nm, −1986.9 ps/nm, −3960.7 ps/nm and −5285.4 ps/nm, respectively, which is shown in Table 1. The echo signals from the targets are all converted into narrow pulses in the FrFD at the proper order, thus the range profiles are obtained and shown in Fig. 4(e)-(h). From these graphs, it can be seen that the −3 dB main-lobe width keeps 1.1 cm and the PSLR keeps 13.2 dB when the velocity varies. This −3 dB main-lobe width is basically the same as the theoretical value of $B$ is the bandwidth of the transmitted signal). Obviously, the PFrFTer is immune to the Doppler dispersion effect and achieves lossless-resolution range profiles of moving targets.

 figure: Fig. 4.

Fig. 4. The range profiles obtained from the optical de-chirping receiver when the velocities of the target are (a) 0 m/s, (b) 56.91 m/s, (c) 113.45 m/s (d) 151.4 m/s. The range profiles obtained from the PFrFTer receiver when the velocities of the target are (e) 0 m/s, (f) 56.91 m/s, (g) 113.45 m/s (h) 151.4 m/s.

Download Full Size | PDF

Tables Icon

Table 1. The velocities of the targets and the corresponding dispersion values

To further validate the range resolution, it is assumed that there are two targets in space with a velocity of 113.45m /s. The range profiles are obtained from the optical de-chirping receiver and the PFrFTer receiver when the targets are separated by different distances. Figures 5(a) and 5(b) show the range profiles from the optical de-chirping receiver when the targets are separated by 5 cm and 1.4 cm, respectively. It is hard to distinguish the two targets for the optical de-chirping receiver when the distance between the targets is less than 5 cm. Figures 5(c) and 5(d) show the corresponding range profiles from the PFrFTer receiver, in which the dispersion is set to −3960.7 ps/nm. From these graphs, we can see the two targets can be distinguished clearly by the PFrFTer even if they are separated by 1.4 cm. The distances between them are measured to be 4.9 cm and 1.6 cm when they are separated by 5 cm and 1.4 cm, respectively. The measurement error of the distances is mainly caused by the noise of the receiver and the interference between the frequency spectrums of the echo signals from different targets. This experiment shows the resolution of the PFrFTer can be as high as 1.4 cm for high-speed targets, which is basically consistent with the theoretical value of 1.3 cm.

 figure: Fig. 5.

Fig. 5. The range profiles obtained by the de-chirping receiver when the targets are separated by (a) 5 cm and (b) 1.4 cm. The range profiles obtained by the PFrFTer receiver when the targets are separated by (c) 5 cm and (d) 1.4 cm.

Download Full Size | PDF

The two-dimensional imaging capability of PFrFTer is verified by ISAR imaging experiments. The target is a small quadcopter with four scattering centers. These scattering centers are evenly distributed on a circle with a diameter of 12 cm, as shown in Fig. 6(a). First, the distance ${R_0}$ between the target and the radar is 100 m. The target moves at a 45-degree angle to the initial line of the radar's sight with a speed of ${v_1} = 160.44$m/s. The radial velocity relative to the radar is calculated as 113.45 m/s. The accumulation time of ISAR imaging is 0.3 s. According to it, the rotation angle $\delta \theta$ is calculated as 14.25°.

 figure: Fig. 6.

Fig. 6. (a) Schematic diagram of the target shape. (b) The trajectory of the target.

Download Full Size | PDF

Two groups of experiments are implemented based on the optical de-chirping and the PFrFTer. The range-Doppler algorithm is applied to process the data. The range profiles are obtained at each pulse firstly. Figures 7(a) and 7(d) show the range profiles of the first pulse from the optical de-chirping receiver and the PFrFTer receiver, respectively. The latter can clearly distinguish four scattering centers, while the former cannot, which again proves the resolution deterioration of the optical de-chirping. Then range alignment is conducted on the range profiles and the results are shown in Fig. 7(b) and 7(e). Obviously, the range alignment result based on the optical de-chirping is worse than that based on the PFrFTer because of the deteriorated range resolution. Further, the error of the range alignment will affect the azimuth compression. The ISAR imaging based on the optical de-chirping is shown in Fig. 7(c). This image is blurred, and it is difficult to recognize the shape of the target. However, due to the high-resolution range profile based on the PFrFTer, high precision range alignment can be achieved and it further guarantees effective azimuth compression. As a result, clear imaging with high resolution in range and cross-range direction is obtained, as shown in Fig. 7(f). Generally, the image entropy is applied to measure the focus quality of the ISAR image [43]. The smaller the value of the image entropy, the better the focus quality of the ISAR image. Here, the image entropies obtained from the ISAR images showed in Fig. 7(c) and 7(f) are 7.1 and 6.2, respectively. These results show the PFrFTer is a promising technology in the recognition of moving targets.

 figure: Fig. 7.

Fig. 7. (a) The range profile, (b) the result of range alignment and (c) the ISAR imaging from the optical de-chirping receiver. (d) The range profile, (e) the result of range alignment and (f) the ISAR imaging from the PFrFTer receiver.

Download Full Size | PDF

In order to get the imaging of high-speed targets with a lossless and high resolution based on the PFrFTer, the order of the PFrFTer must satisfy the relationship:$p = 2arc\cot \left( { - {{\left( {1 + \frac{{2v}}{c}} \right)}^2}{k_1}} \right)/\pi$. From this property, we can see that the proposed system is more suitable to be applied to the detection and recognition of a large single target or multiple targets with the relatively same speed, such as missiles [3,44]. Besides, we should obtain the velocity of the target firstly before getting the imaging of the target. As for cooperative targets, it is easy to obtain the velocity. As for non-cooperative targets, the velocity measurement can be realized by many simple systems such as lidar [45]. With these systems, the radar system based on the PFrFTer can realize high-resolution imaging of a moving target.

The order of the PFrFTer is reconstructed by changing the chirp rate of the light, which can be practically realized by a DM. From Eq. (14) and Eq. (15), it is known that the required dispersion value is determined by the target velocity, the chirp rate of the transmitted signal and the wavelength of the light. A simulation is conducted to show this relationship. In the simulation, the wavelength of the light is set to a normal value, 1545 nm. Figure 8 shows the required dispersion as a function of the chirp rate of the transmitted signal and the target velocity under single lower sideband modulation. From the graph, we can see that the absolute dispersion value is less than 6×104 ps/nm when the velocity is smaller than Mach 10, and the chirp rate is larger than 0.1 GHz/µs, which are feasible for most scenarios and imaging radar systems. The dispersion can be provided by the tunable dispersion compensation module, DCF or other dispersive devices, the dispersion of which can be as large as −70790 ps/(nm·km) [46].

 figure: Fig. 8.

Fig. 8. Numerically simulated dispersion as a function of the chirp rate of the transmitted signal and the target velocity.

Download Full Size | PDF

In the actual application of the PFrFTer in radar imaging, there are errors between the obtained velocity and the real velocity of the target, as well as between the set dispersion and the calculated ideal dispersion. These errors will deteriorate the resolution. From [27], it is known that the −3 dB main-lobe width of the range profile will hardly expand when the velocity is less than $c/2BT$. It means the Doppler dispersion effect can be ignored when the deviation between the set dispersion and the optimal dispersion is less than $({4/{k_1}{\lambda^2}} )\cdot ({c/2BT} )$ according to Eq. (14) and (15), or in other words, when the set dispersion D and the real velocity v satisfy following equation

$$\left|{D \pm \frac{{4v}}{{{k_1}{\lambda^2}}}} \right|< \left|{\frac{4}{{{k_1}{\lambda^2}}} \cdot \frac{c}{{2BT}}} \right|,$$
where “+” and “−” are established under lower sideband modulation and upper sideband modulation in the SSB modulation, respectively. Because of ${k_1} = B/T$, Eq. (16) can be rewritten as
$$\left|{D \pm \frac{{4v}}{{{k_1}{\lambda^2}}}} \right|< \left|{\frac{{2c}}{{{B^2}{\lambda^2}}}} \right|.$$

In the ISAR experiments, the radial velocity of the target is continuously increasing in the detection processing. During the accumulation time, the velocity and the dispersion satisfy Eq. (17), so the imaging still keeps a high resolution, which can be seen in Fig. 7(d) and 7(f).

4. Conclusions

In conclusion, we have proposed and demonstrated a radar receiver based on a reconfigurable PFrFTer, aiming to process the wideband echo signals, eliminate the Doppler dispersion effect and obtain high-resolution imaging of targets with high speed. A receiver based on the PFrFTer is established and investigated in the experiments. Advantages of the proposed PFrFTer compared with the traditional de-chirping receiver are demonstrated through detection and imaging of different targets with high speed. The proposed technique would be of interest in the detection and recognition of targets.

Funding

National Key Research and Development Program of China (2019YFB2203301); National Natural Science Foundation of China (61690191, 61690192, 61621064).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but maybe obtained from the authors upon reasonable request.

References

1. H. Cheng, S. Li, H. Zheng, H. Jing, and H. Sun, “A W-band auto-focus holographic imaging system for security screening,” IEICE Electron. Express 14(11), 20170347 (2017). [CrossRef]  

2. T. Kawanishi, A. Kanno, N. Yamamoto, N. Yonemoto, N. Shibagaki, and K. Kashima, “Optical fiber network-connected distributed mm-wave radar system,” in 2017 IEEE Photonics Society Summer Topical Meeting Series (IEEE, 2017), pp. 197–198.

3. W. W. Camp, J. T. Mayhan, and R. M. O’Donnell, “Wideband radar for ballistic missile defense and range-Doppler imaging of satellites,” Lincoln Lab. J. 12(2), 267–280 (2000).

4. T. Long, Z. Liang, and Q. Liu, “Advanced technology of high-resolution radar: target detection, tracking, imaging, and recognition,” Sci. China Inf. Sci. 62(4), 40301 (2019). [CrossRef]  

5. J. Vankka and K. A. I. Halonen, Direct digital synthesizers: theory, design and applications (Springer Science & Business Media, 2013).

6. P. Ghelfi, F. Laghezza, F. Scotti, G. Serafino, S. Pinna, D. Onori, E. Lazzeri, and A. Bogoni, “Photonics in radar systems: RF integration for state-of-the-art functionality,” IEEE Microw. 16(8), 74–83 (2015). [CrossRef]  

7. M. Li and N. Zhu, “Recent advances in microwave photonics,” Front. Optoelectron. 9(2), 160–185 (2016). [CrossRef]  

8. J. Liao, H. Wen, X. Zheng, H. Zhang, Y. Guo, and B. Zhou, “Novel 2N bit bipolar photonic digital-to-analog converter based on optical DQPSK modulation coupled with differential detection,” Opt. Lett. 37(9), 1502–1504 (2012). [CrossRef]  

9. C. Wang and J. Yao, “Photonic generation of chirped millimeter-wave pulses based on nonlinear frequency-to-time mapping in a nonlinearly chirped fiber Bragg grating,” IEEE Trans. Microwave Theory Techn. 56(2), 542–553 (2008). [CrossRef]  

10. Y. Tong, D. Han, R. Cheng, Z. Liu, W. Xie, J. Qin, and Y. Dong, “Photonics-based coherent wideband linear frequency modulation pulsed signal generation,” Opt. Lett. 43(5), 1023–1026 (2018). [CrossRef]  

11. X. Wang, Z. Liu, S. Wang, D. Sun, Y. Dong, and W. Hu, “Photonic radio-frequency dissemination via optical fiber with high-phase stability,” Opt. Lett. 40(11), 2618–2621 (2015). [CrossRef]  

12. H. Wang, X. Xue, S. Li, and X. Zheng, “All-optical arbitrary-point stable quadruple frequency dissemination with photonic microwave phase conjugation,” IEEE Photonics J. 10(4), 1–8 (2018). [CrossRef]  

13. D. Marpaung, B. Morrison, M. Pagani, R. Pant, D. Choi, B. Luther-Davies, S. J. Madden, and B. J. Eggleton, “Low-power, chip-based stimulated Brillouin scattering microwave photonic filter with ultrahigh selectivity,” Optica 2(2), 76–83 (2015). [CrossRef]  

14. L. Li, X. Yi, S. Song, S. X. Chew, R. Minasian, and L. Nguyen, “Microwave photonic signal processing and sensing based on optical filtering,” Appl. Sci. 9(1), 163 (2019). [CrossRef]  

15. S. Pan and Y. Zhang, “Tunable and wideband microwave photonic phase shifter based on a single-sideband polarization modulator and a polarizer,” Opt. Lett. 37(21), 4483–4485 (2012). [CrossRef]  

16. Y. Xu, S. Li, X. Xue, X. Xiao, X. Zheng, and B. Zhou, “An interleaved broadband photonic ADC immune to channel mismatches capable for high-speed radar imaging,” IEEE Photonics J. 11(4), 1–9 (2019). [CrossRef]  

17. R. Li, W. Li, M. Ding, Z. Wen, Y. Li, L. Zhou, S. Yu, T. Xing, B. Gao, and Y. Luan, “Demonstration of a microwave photonic synthetic aperture radar based on photonic-assisted signal generation and stretch processing,” Opt. Express 25(13), 14334–14340 (2017). [CrossRef]  

18. P. Ghelfi, F. Laghezza, F. Scotti, G. Serafino, A. Capria, S. Pinna, D. Onori, C. Porzi, M. Scaffardi, and A. Malacarne, “A fully photonics-based coherent radar system,” Nature 507(7492), 341–345 (2014). [CrossRef]  

19. S. Peng, S. Li, X. Xue, X. Xiao, D. Wu, X. Zheng, and B. Zhou, “High-resolution W-band ISAR imaging system utilizing a logic-operation-based photonic digital-to-analog converter,” Opt. Express 26(2), 1978–1987 (2018). [CrossRef]  

20. F. Zhang, Q. Guo, Z. Wang, P. Zhou, G. Zhang, J. Sun, and S. Pan, “Photonics-based broadband radar for high-resolution and real-time inverse synthetic aperture imaging,” Opt. Express 25(14), 16274–16281 (2017). [CrossRef]  

21. A. Wang, J. Wo, X. Luo, Y. Wang, W. Cong, P. Du, J. Zhang, B. Zhao, J. Zhang, and Y. Zhu, “Ka-band microwave photonic ultra-wideband imaging radar for capturing quantitative target information,” Opt. Express 26(16), 20708–20717 (2018). [CrossRef]  

22. F. Zhang, B. Gao, and S. Pan, “Photonics-based MIMO radar with high-resolution and fast detection capability,” Opt. Express 26(13), 17529–17540 (2018). [CrossRef]  

23. S. Maresca, F. Scotti, G. Serafino, L. Lembo, A. Malacarne, F. Falconi, P. Ghelfi, and A. Bogoni, “Coherent MIMO radar network enabled by photonics with unprecedented resolution,” Opt. Lett. 45(14), 3953–3956 (2020). [CrossRef]  

24. F. Berland, T. Fromenteze, D. Boudescoque, P. Di Bin, H. H. Elwan, C. Aupetit-Berthelemot, and C. Decroze, “Microwave Photonic MIMO Radar for Short-Range 3D Imaging,” IEEE Access 8, 107326–107334 (2020). [CrossRef]  

25. X. Xiao, S. Li, X. Xue, L. Xing, S. Peng, X. Zheng, and B. Zhou, “Photonics-assisted broadband distributed coherent aperture radar for high-precision imaging of dim-small targets,” IEEE Photonics J. 11(5), 1–9 (2019). [CrossRef]  

26. J. Dong, F. Zhang, Z. Jiao, Q. Sun, and W. Li, “Microwave photonic radar with a fiber-distributed antenna array for three-dimensional imaging,” Opt. Express 28(13), 19113–19125 (2020). [CrossRef]  

27. W. R. Remley, “Doppler dispersion effects in matched filter detection and resolution,” Proc. IEEE 54(1), 33–39 (1966). [CrossRef]  

28. C. Pang, S. Liu, and Y. Han, “High-speed target detection algorithm based on sparse Fourier transform,” IEEE Access 6, 37828–37836 (2018). [CrossRef]  

29. J. Zhang, W. Jiang, Y. Yu, and X. Zhang, “Photonics-based simultaneous measurement of distance and velocity using multi-band LFM microwave signals with opposite chirps,” Opt. Express 27(20), 27580–27591 (2019). [CrossRef]  

30. C. Zeng, W. Zhu, and X. Jia, “Bistatic ISAR Sparse Imaging Method for High-Speed Moving Target Based on Dechirping Processing,” Int. J. Antenn. Propag. 2019, 1–15 (2019). [CrossRef]  

31. L. Yang, M. Xing, L. Zhang, G. Sun, Y. Gao, Z. Zhang, and Z. Bao, “Integration of rotation estimation and high-order compensation for ultrahigh-resolution microwave photonic ISAR imagery,” IEEE Trans. Geosci. Remote Sensing 59(3), 2095–2115 (2021). [CrossRef]  

32. L. B. Almeida, “The fractional Fourier transform and time-frequency representations,” IEEE Trans. Signal Process. 42(11), 3084–3091 (1994). [CrossRef]  

33. A. W. Lohmann and D. Mendlovic, “Fractional Fourier transform: photonic implementation,” Appl. Opt. 33(32), 7661–7664 (1994). [CrossRef]  

34. C. Cuadrado-Laborde, A. Carrascosa, A. Diez, J. L. Cruz, and M. V. Andres, “Photonic fractional Fourier transformer with a single dispersive device,” Opt. Express 21(7), 8558–8563 (2013). [CrossRef]  

35. C. Schnebelin and H. G. de Chatellus, “Agile photonic fractional Fourier transformation of optical and RF signals,” Optica 4(8), 907–910 (2017). [CrossRef]  

36. D. Peng, H. Jiang, F. Shen, H. Xia, and K. Xie, “Reconfigurable high-resolution photonic fractional Fourier transformer for broadband frequency measurement and large-scale chirp rate characterization of microwave signals,” Opt. Commun. 460, 125157 (2020). [CrossRef]  

37. G. Han, S. Li, X. Xue, and X. Zheng, “Photonic fractional Fourier transformer for chirp radar with ghost target elimination,” Opt. Lett. 45(15), 4228–4231 (2020). [CrossRef]  

38. G. Han, S. Li, X. Xue, and X. Zheng, “Photonic chirp rates estimator for piecewise linear frequency modulated waveforms based on photonic self-fractional Fourier transform,” Opt. Express 28(15), 21783–21791 (2020). [CrossRef]  

39. P. Koivumäki, Triangular and Ramp Waveforms in Target Detection with a Frequency Modulated Continuous Wave Radar (Aalto University, 2017).

40. V. C. Chen and M. Martorella, Inverse Synthetic Aperture Radar Imaging: Principles, Algorithms, and Applications (SciTech Publishing, 2014).

41. Y. Zha, X. Xue, S. Li, X. Zheng, and B. Zhou, “Wideband Radar Signal Distribution With an Idler-Free Photonic Microwave Frequency Shifter,” IEEE Photonics Technol. Lett. 30(22), 1948–1951 (2018). [CrossRef]  

42. P. Hu, P. Hu, S. Xu, and Z. Chen, “Wideband radar system distortion compensation using spherical satellite echo,” in Sixth Asia-Pacific Conference on Antennas and Propagation (IEEE, 2017), pp. 1–3.

43. S. Zhang, Y. Liu, and X. Li, “Fast Entropy Minimization Based Autofocusing Technique for ISAR Imaging,” IEEE Trans. Signal Process. 63(13), 3425–3434 (2015). [CrossRef]  

44. Z. Wang, F. Yan, F. He, and J. Zhu, “Missile target automatic recognition from its decoys based on image time-series,” Pattern Recogn. 43(6), 2157–2164 (2010). [CrossRef]  

45. R. Diaz, S. Chan, and J. Liu, “Lidar detection using a dual-frequency source,” Opt. Lett. 31(24), 3600–3602 (2006). [CrossRef]  

46. B. Wang, W. Zheng, R. Chen, H. Wei, B. Li, and J. Hsu, “Dispersion compensation photonic crystal fibers with ultra-high dispersion coefficient and extremely low loss,” in 2013 6th International Conference on Advanced Infocomm Technology (IEEE, 2013), pp. 97–98.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but maybe obtained from the authors upon reasonable request.

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (8)

Fig. 1.
Fig. 1. The FrFT in the time-frequency domain. (a) Square modulus of the FrFT of a given signal at a specific order $2\alpha /\pi$ . (b) Square modulus of the FrFT of the echo signals from a stationary target and a moving target at proper order $2{\alpha _1}/\pi$ and $2{\alpha _2}/\pi$ , respectively.
Fig. 2.
Fig. 2. Schematic diagram of the radar receiver based on the PFrFTer. Ant: antenna; LD: laser diode; SSB: single-sideband; DM: dispersion module; E/O: electro-optical conversion; O/E: opto-electronic conversion. FT: Fourier transform.
Fig. 3.
Fig. 3. (a) The spectrum of the transmitted signal. (b) The spectrum of the optical single-sideband signal.
Fig. 4.
Fig. 4. The range profiles obtained from the optical de-chirping receiver when the velocities of the target are (a) 0 m/s, (b) 56.91 m/s, (c) 113.45 m/s (d) 151.4 m/s. The range profiles obtained from the PFrFTer receiver when the velocities of the target are (e) 0 m/s, (f) 56.91 m/s, (g) 113.45 m/s (h) 151.4 m/s.
Fig. 5.
Fig. 5. The range profiles obtained by the de-chirping receiver when the targets are separated by (a) 5 cm and (b) 1.4 cm. The range profiles obtained by the PFrFTer receiver when the targets are separated by (c) 5 cm and (d) 1.4 cm.
Fig. 6.
Fig. 6. (a) Schematic diagram of the target shape. (b) The trajectory of the target.
Fig. 7.
Fig. 7. (a) The range profile, (b) the result of range alignment and (c) the ISAR imaging from the optical de-chirping receiver. (d) The range profile, (e) the result of range alignment and (f) the ISAR imaging from the PFrFTer receiver.
Fig. 8.
Fig. 8. Numerically simulated dispersion as a function of the chirp rate of the transmitted signal and the target velocity.

Tables (1)

Tables Icon

Table 1. The velocities of the targets and the corresponding dispersion values

Equations (17)

Equations on this page are rendered with MathJax. Learn more.

F p { x ( t ) } ( u ) = + x ( t ) K α ( u , t ) d t α m π ( m Z ) ,
K α ( u , t ) = 1 j cot α exp [ j 2 π ( u csc α t + cot α 2 t 2  +  u 2 cot α 2 ) ] ,
s r ( t )  =  r e c t ( t τ ( t ) T ) exp [ j ( 2 π f 1 ( t τ ( t ) ) + π k 1 ( t τ ( t ) ) 2 ) ] r e c t ( t 2 R / c T ) exp [ j ( 2 π f 2 ( t τ 2 ) + π k 2 ( t τ 2 ) 2 ) ] ,
F p { s r ( t ) } ( u )  =  + s r ( t ) K α ( u , t ) d t .
F p { s r ( t ) } ( u ) r e c t ( u csc α  +  k 2 τ 2 f 2 ( k 2 + cot α ) T ) .
R = ( c + 2 v ) ( f 2 u csc α ) / ( 2 k 2 ) .
F p { s r ( t ) } ( u ) sinc [ T ( u csc α  +  k 2 τ 2 f 2 ) ] .
E 1 E 0 + E 0 s ( t )  =  E 0 o p t i c a l c a r r i e r  +  β 1 E 0 exp [ j ( 2 π f s t + π k s t 2  +  φ s ) ] c o h e r e n t l i n e a r c h i r p l i g h t ,
f s = c f 1 / ( c  +  k 1 D λ 2 ) k s = c k 1 / ( c  +  k 1 D λ 2 ) φ s = π f 1 2 D λ 2 / ( c  +  k 1 D λ 2 ) ,
i ( f ) = + s r ( t ) exp [ j ( 2 π f s t + π k s t 2  +  φ s ) ] exp ( j 2 π f t ) d t .
| i ( u ) | | + s r ( t ) exp [ j 2 π ( u csc α t + cot α 2 t 2  +  u 2 cot α 2 ) ] | d t | F p [ s r ( t ) ] ( u ) | .
p = 2 a r c cot ( k s ) / π = 2 a r c cot ( c c  +  k 1 D λ 2 k 1 ) / π .
R = ( c + 2 v ) ( f + f s f 2 ) / ( 2 k 2 ) .
D = c ( ( c + 2 v ) 2 c 2 ) k 1 λ 2 ( c + 2 v ) 2 4 v k 1 λ 2 .
D = c ( ( c + 2 v ) 2 c 2 ) k 1 λ 2 ( c + 2 v ) 2 4 v k 1 λ 2 .
| D ± 4 v k 1 λ 2 | < | 4 k 1 λ 2 c 2 B T | ,
| D ± 4 v k 1 λ 2 | < | 2 c B 2 λ 2 | .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.