1. Introduction
In the last two decades, modern wireless communication system technologies have undergone rapid growth in response to the corresponding demand [
1]. The multipath effect, a prevalent phenomenon in wireless communication, occurs when radio signals arrive at the receiving antenna via multiple paths. This effect can arise from factors such as atmospheric ducting, ionospheric reflection, and reflections from surfaces like water bodies, mountains, and buildings. Interference and phase shifts can ensue when the same signal follows different paths to the receiver. Destructive interference resulting from this can cause signal fading, which may render a signal too weak to be reliably received in certain areas. This phenomenon is one of the most important causes of system fading. On the other hand, if multipath effects are utilized effectively, communication efficiency can also be improved. For example, using drones as relay stations can create intentional multipath effects, enabling transmission across obstacles in complex environments [
2,
3] or ensuring reliable long-distance communication [
4,
5].
Consequently, a comprehensive understanding of multipath fading channels is crucial for the efficient design of wireless communication systems [
6]. In statistical channel models, distribution parameters are typically derived from extensive field measurement data [
7], which are used to characterize typical scenarios. Several classical models have been established, such as the 3GPP/ITU 2D/3D channel models and the 3D millimeter-wave channel model [
8,
9]. However, statistical models frequently fail to precisely capture channel performance at specific locations, leading to the widespread adoption of the ray-tracing method. As shown in
Figure 1, the ray-tracing channel model simplifies the wavefront into multiple small ray tubes that propagate along their respective incident directions [
10]. It is a deterministic model based on high-frequency electromagnetic simulation technology and uses Maxwell’s equations and appropriate boundary conditions to develop a deterministic model that can ascertain the details of multipath propagation and derive the impulse response of multipath fading for specific environmental conditions [
11,
12,
13].
Despite its advantages in providing accurate and detailed modeling of multipath propagation, ray tracing has several inherent drawbacks. One major limitation is its high computational complexity. Additionally, ray-tracing models are highly dependent on detailed and accurate environmental information. Constructing a precise geometric model of the propagation environment, including the shapes, positions, and material properties of all relevant obstacles, can be labor-intensive and time-consuming. Furthermore, while ray tracing can provide detailed insights into a signal’s propagation paths, it does not inherently account for all wave phenomena, such as wave polarization and material anisotropy, which might be relevant in certain scenarios.
Given these challenges, improving the accuracy of multipath channel modeling remains a significant concern in the wireless communication and radar community. Some researchers are attempting to improve the ray-tracing channel model. In [
14], Yi Chen et al. explored combining ray-tracing techniques with terahertz communication to enhance multipath channel modeling. Their hybrid approach outperformed conventional models by better capturing the temporal–spatial characteristics of multipath channels and the role of ray tracing in accurate multipath propagation. However, it still faced challenges related to the high computational demands of ray tracing. In [
15], Yuanzhe Wang et al. introduced Massive MIMO Communications, emphasizing the need to account for reflection, scattering, and diffraction to obtain more precise results. But they were unable to overcome the challenge of modeling multipath channels in complex environments.
Some other researchers hope to improve the accuracy of multipath channel models by using more novel techniques. Anh Hong Nguyen et al. developed a higher-precision model for specular multipath components (SMCs) using gaussian processes, thereby improving the accuracy of multipath channel modeling for ultra-wideband (UWB) and mm-wave signals [
16]. However, in some complex scenarios, Non-Line-of-Sight (NLoS) conditions can severely distort UWB signals, thereby affecting modeling accuracy [
17]. To address this issue, Qiu Wang et al. introduced a 1D-ConvLSTM-Attention network (1D-CLANet) enhanced by Squeeze-and-Excitation (SE). This SE-enhanced 1D-CLANet can identify NLoS conditions during modeling. They will use this 1D-CLANet in future research to reduce the modeling errors caused by NLoS conditions [
18].
Among the work of researchers, Abdelfattah Fawky et al. creatively demonstrated the importance of incorporating radar cross sections (RCSs) into channel modeling, creating a high-precision 3D deterministic model that provided more accurate predictions of signal behavior in complex environments [
19]. This suggests that by accounting for the scattering processes occurring along the signal path in more detail, including the precise RCS values of various scatterers, the overall channel model can be made significantly more accurate. In our paper, we aim to obtain higher-precision RCS values by introducing bistatic testing, and we hope that this high-precision RCS can contribute to developing multipath channel modeling methods that are accurate, have low computational complexity, and are less affected by complex environments.
The RCS of a target represents the effective area perceived by the radar. The formula for defining the RCS is shown in Equation (1) [
20]. It is the ratio of the power reflected by a target towards the emission source within a unit solid angle to the plane wave power density incident on the target from a given direction.
In this formula,
denotes the RCS, typically measured in
or
;
represent the distances between the target and the radar’s transmitting and receiving antennas, respectively;
is the electric field strength of the scattered wave at the receiving antenna; and
is the strength of the incident wave at the target. The RCS is the hypothetical area that intercepts a specific amount of power, which, when scattered uniformly in all directions, produces an echo at the radar station equivalent to that of the actual target. The RCS is the hypothetical area that intercepts a quantity of power such that, when scattered uniformly in all directions, it generates an echo at the radar station equivalent to that of the target. In wireless communication, the radar’s transmitting antenna can be equivalent to the radio wave’s transmitting station, and the radar receiving antenna can be equivalent to the radio wave’s receiving station. Thus, the RCS can characterize the scattering of radio waves by obstacles along the propagation path. The RCS values of some common objects are shown in
Table 1 below.
In practical engineering, the RCS definition formula is generally not used. Instead, the radar equation, as shown in Equation (2), is used for RCS measurement [
22]. In this formula,
represent the radar’s transmitted power and received power, respectively, while
correspond to the gains of the transmitting and receiving antennas, respectively.
denotes the radar’s operating wavelength.
Traditional RCS testing systems typically use a monostatic radar configuration, where the transmitting and receiving antennas are located at the same radar station. However, with continuous technological advancements, bistatic radar RCS measurements are attracting attention due to their higher accuracy and efficiency [
23]. The RCS testing system must ensure coherence between the transmitted and received signals. In monostatic radar systems, where the transmitter and receiver are co-located at the same base station, time and frequency synchronization is readily achievable through the use of a shared clock source, ensuring signal coherence. In bistatic radar systems, where the transmitter and receiver are spatially separated, each must use its own clock. The differences between these clocks make it difficult to coordinate transmission and reception. There are two common methods: One consists of equipping both bases with high-precision atomic clocks. The coherent population trapping (CPT) atomic clock, widely used in synchronization communications, has a frequency stability of approximately
, while the highest-precision rubidium atomic clocks can achieve a frequency stability of
[
24]. However, the production and maintenance costs of high-precision atomic clocks are extremely high. Another method consists of extracting a synchronization signal from one base and transmitting it to another base, thereby achieving synchronization between the two bases. The key to this method is ensuring that the synchronization signal maintains phase stability after long-distance transmission between the two bases.
This article presents a signal transmission system designed for a bistatic RCS testing system to address the issue of time synchronization. As the testing system is an electromagnetic measurement system and the distance between the two radar stations is over 300 m, using wireless electrical signal transmission could interfere with the scattered waves, affecting the test results. On the other hand, using wired electrical signal transmission over such distances would lead to significant attenuation. Therefore, our design adopts an optoelectronic conversion approach. Before transmission, the reference signal is converted into an optical signal, transmitted via optical fiber to the other base station, and then demodulated back into an electrical signal. The optoelectronic conversion system is linear and does not disrupt the coherence of the entire testing system. The linearity of this system is a key indicator in determining whether RCS measurements can be performed. In the following
Section 3, the phase stability of the system is specifically tested, and in
Section 4, the negative impacts that system nonlinearity may have on the overall RCS measurements are analyzed.
2. System Principle and Simulation
The entire testing system is an inter-pulse frequency-hopping radar, and its operating principle is illustrated in
Figure 2. A local oscillator signal provides the time base for two agile frequency signal sources, ensuring time synchronization between them. Agile frequency signal source 1 generates a signal
with a frequency
, which is split into two paths. One path is used as the transmission signal, which, after modulation, amplitude control, and power amplification, is transmitted towards the target by the antenna. The other path serves as the intermediate frequency (IF) for down-conversion to obtain the signal
. The signal
generated by agile frequency signal source 2, with a frequency
, is also split into two paths. One path is used as the reference local oscillator and mixed with the processed echo signal
to obtain signal
. The other path serves as the RF input for down-conversion to obtain the signal
. Thus, the frequency of signal
should be
. Finally, the signal
this time acts as the IF, and the other signal,
, undergoes a series of processes, including sampling, digital phase locking, quadrature digital down-conversion, and baseband signal post-processing. This process allows one to obtain various measurements, such as the target’s RCS, one-dimensional image, and two-dimensional image.
Additionally, the signal transmission system can theoretically be connected at any point in the entire testing system to increase the extent of signal transmission. A schematic diagram of this system is shown in
Figure 3.
The modulated signal, , is input into an electro-optical modulator, where it modulates the optical carrier signal, , generated within the modulator, forming the modulated optical signal, . After long-distance optical fiber transmission, undergoes a phase shift and attenuation, resulting in the optical signal . then enters a photodetector and is converted back into the electrical signal, .
2.1. Selection of an Electro-Optic Modulator
Optical modulators encode an electrical signal onto an optical wave, known as the optical carrier. This encoding can manipulate the amplitude, phase, frequency, polarization, or any combination of these properties of the optical carrier [
25]. To meet the demand for long-distance and high-capacity service in modern optical transmission networks, high speed, compactness, and low power have become the mainstream development themes for future electro-optical modulators [
26,
27]. Among the various electro-optical modulators, the Mach–Zehnder modulator has become one of the most widely used types in recent years due to its advantages: high speed, low cost, and low complexity [
28].
The Mach–Zehnder modulator (MZM) is an interferometric device constructed from a material exhibiting a strong electro-optic effect (such as
). Lithium niobate (
)-based MZMs are commonly used for modulation in communication systems with a capacity of up to
. However, their application is constrained by high drive voltages, a large physical footprint, and a lack of integrability [
29]. In contrast, indium phosphide (
)-based MZMs offer higher 3 dB bandwidths and require lower drive voltages [
30]. But they are also more expensive than the former. The structure of the MZM is illustrated in
Figure 4. The optical input
is divided between the upper and lower modulator arms, where phase modulation occurs via two phase shifters driven by the electrical signals
. The modulated signals are then recombined to produce the optical output
.
The MZM has three operating modes: Maximum Transmission Point (MATP), Quadrature Point (QTP), and Minimum Transmission Point (MITP). By setting an appropriate DC bias voltage, the operating mode of the MZM can be selected. When the bias voltage is set to 0, the MZM is in the MATP mode, and the output optical intensity is at its maximum. In this mode, since the bias voltage is 0, the refractive indices of the upper and lower arms do not change, and the optical signals in both arms are not modulated. When the bias voltage is set to half of the modulator’s half-wave voltage, the MZM is in the QTP mode, and the output optical intensity is half of the maximum optical intensity. The externally applied voltage alters the refractive indices of the upper and lower arms, causing the optical signals in each arm to undergo varying levels of phase modulation. This process induces a phase shift between the optical signals in the two arms, resulting in a coherent superposition within the optical coupler. When the bias voltage is set to the modulator’s half-wave voltage, the MZM is in the MITP mode, and the output optical intensity is 0. In this mode, when the DC bias voltage is equal to the MZM’s half-wave voltage, the optical amplitudes in the upper and lower arms cancel each other out in the optical coupler, resulting in an output optical field of 0.
As the MZM is one of the more mature technologies among the various optoelectronic modulation methods, with relatively simple principles and moderate hardware costs, we employed an MZM operating in single-drive push–pull mode, which modulates the signal in terms of amplitude. In the push–pull-mode MZM, there is only one radio frequency (RF) electrode and one direct current (DC) electrode, as depicted in
Figure 5. Two RF signals of equal amplitude and opposite phases are applied to the upper and lower arms of the MZM. The MZM operates in the QTP state to avoid the need for carrier and phase recovery [
31]. In this state, only the intensity modulation of the input optical signal occurs. The output of the modulator in this operating mode can be represented as follows:
in which
represents the optical carrier signal input to the MZM,
denotes the insertion loss of the MZM,
is the half-wave voltage of the modulator, and
represent the amplitude, angular frequency, and initial phase of the modulating signal
, respectively.
represents the DC bias voltage applied to the upper and lower arms.
2.2. Derivation of the System Formulas
Let us assume the modulating signal
in
Figure 2 is a standard cosine signal:
is the amplitude of the modulating signal.
and
represent the angular frequency (the following will be referred to as “frequency”) and initial phase of the modulating signal, respectively. Since the MZM operates in QTP mode,
. Therefore, the modulated optical signal
obtained after the modulating signal
passes through the electro-optic converter can be represented as follows:
In this formula,
represents the optical carrier signal:
represents the amplitude of the optical carrier signal, and represent the frequency and phase of the optical carrier signal. denotes the insertion loss of the modulator, and is approximately 6 V in this design. is referred to as the modulation index.
According to the Jacobi–Anger expansion,
By substituting Equation (7) into Equation (5), we obtain the following:
Here,
represents the first-kind
n-order Bessel function. In the case of small signal modulation, only the DC component and the first-order component are considered. Therefore, the above equation can be written as follows:
The signal undergoes transmission through a single-mode fiber (SMF) with a length of
, wherein the transfer function of the SMF is given by the following:
represent the length, propagation constant, and amplitude attenuation coefficient of the SMF, respectively. By performing a Taylor expansion on
and using
to denote the propagation constant at point
and the values of the first and second derivatives, we can obtain the following:
From Equations (9)–(11), we can obtain the expression for the modulated optical signal obtained after it propagates through the SMF with a length of
as follows:
Finally, the photocurrent detected by the photodetector is given by the following:
Here,
is the responsivity of the photodetector, and
represents the power attenuation coefficient of the optical wave during transmission in the fiber. Considering the Taylor series expansion of
, the photocurrent can be expressed as follows:
Here, is referred to as the group delay of the carrier, , where , and represent the fiber’s dispersion coefficient, optical wavelength, and the speed of light in a vacuum, respectively. By comparing the input and output of the photoelectric conversion system, it can be observed that, after considering the negative sign in the phase of the first harmonic term of , it can be simplified to . When comparing the phase of the modulating signal , the difference between the two is . The parameter represents the group delay of the entire system, which is solely determined by the system itself and the length of the optical fiber, making a fixed constant.
Moreover, as shown by Equation (12), it is evident that the primary factor determining the intensity of the n-th harmonic is the n-th order Bessel function value of the modulation index
. As shown in
Figure 6, in the case of small-signal modulation, where m is relatively small,
for
becomes significantly smaller than
and
; thus, it can be considered negligible.
2.3. System Simulation
In the previous section, we concluded through theoretical derivation that the system fundamentally meets the phase stability requirements. In this stage, we used Matlab 2020a to conduct simulations to gain an intuitive understanding of the waveform and phase variations at each node of the system as well as observe the impact of various component parameters on system gain and waveform.
Firstly, we set the simulation parameters for the system. Since the carrier used in the photoelectric conversion is typically a laser with a wavelength of 1550 nm, which corresponds to a frequency of approximately
, and the frequency range of the modulated electrical signal is
, it follows that in practical applications,
. Thus, for the modulating signal, we set the parameters as follows: amplitude
, frequency
, and initial phase
. For the optical carrier signal, we set the amplitude
, carrier signal frequency
to satisfy the condition
for amplitude modulation, and carrier signal initial phase
. The half-wave voltage of a
MZM is generally between 3 and 10 V. A lower half-wave voltage can enhance the performance of the MZM by increasing the modulation speed, reducing the heat generation, and lowering the drive voltage requirements, thereby improving overall efficiency. However, MZMs with lower half-wave voltages require more complex fabrication processes and incur higher production costs. After comprehensive consideration, the half-wave voltage of the MZM was set to
, which is also the voltage used for the actual product. Thus, the modulation index
. The insertion loss is
. A waveform comparison between the modulated signal
and the modulating signal
is shown in
Figure 7.
In the figure above, the blue and green curves represent the complete modulated signal and the signal with only the DC and first-order components retained, respectively. The red curve represents the modulating signal . It is evident from the figure that the envelopes of the blue and green curves almost entirely match the red curve. Unlike traditional amplitude modulation, the envelope of the modulated signal is no longer perfectly aligned with the phase of the modulated signal but instead undergoes a phase inversion. This is due to the inherent principles of the MZM, as shown in Equation (9). Comparing the blue and green curves shows that they overlap almost completely, with only slight differences around the peaks, confirming the approximation mentioned in Equation (9) that “In the case of small signal modulation, only the DC component and the first-order component are considered”.
Regarding the case where the modulation index is increased to
, a comparison between the true values and the approximate values is shown in
Figure 8. At this point, there is a significant difference between the true and approximate values, and neither envelope matches the modulated signal anymore, indicating that the intensity modulation of a single push–pull MZM is only applicable when the modulating signal amplitude is less than the half-wave voltage of the modulator. Therefore, in practical applications, it is essential to control the amplitude of the modulating signal.
For the optical fiber in the system, the simulation parameters were set as follows: the attenuation coefficient of the fiber
, the refractive index
, the operating wavelength
, and the dispersion coefficient
. The length was set to
, and a rough calculation revealed that the group delay is
. For the photodetector, the responsivity is a positive number that is always less than 1. For convenience in calculations, it was set to
. A comparison between the system’s output signal
and input signal
is shown in
Figure 9:
In the figure, it can be observed that the demodulated signal has an added DC component and a higher amplitude than the modulating signal, while the frequency remains almost identical to the input signal. Due to the group delay
being only 2420 ns, the phase difference
in
caused by the group delay is negligible, resulting in a phase difference of
between the two signals. Performing a discrete Fourier transform with 18,850 points on the demodulated signal yields the results shown in
Figure 10.
Upon comparing the DFT results of the demodulated and modulating signals, it becomes evident that the proportion of the second harmonic component in the demodulated signal is much smaller than that of the first harmonic and the DC component. This also confirms that neglecting the higher-order components of the demodulated signal, as mentioned in Equation (14), is a reasonable approximation.
The simulation results corroborate the theoretical derivations, with both confirming that the system is fundamentally feasible. These findings provided a solid foundation for the subsequent practical measurements. Moving forward, real-world testing will further validate the system’s feasibility and help identify any potential issues that may not have been captured in the theoretical analysis.
4. System Application
The actual structure of the bistatic testing system after integrating the signal transmission system into the entire testing setup is shown in
Figure 18. Except for the transmitting antenna and the photodetector in the signal transmission system located at the Transmitting Base, all the other components of the testing system are situated at the Receiving Base. The Core contains the local oscillator clock source and the two agile frequency signal sources. It also performs mixing to obtain the intermediate frequency reference signal, as well as the final sampling, quadrature demodulation, and other processing tasks.
Via analysis in conjunction with
Figure 2, it can be gleaned that for a series of n frequency incrementation pulses with a carrier frequency of
, the transmission of a specific pulse can be represented as follows:
where
is the base frequency, which can be set according to the testing requirements, with the test system covering a frequency range of 1–40 GHz;
is the inter-pulse frequency increment, set to 1 MHz;
is the pulse width, determined by the maximum size of the target; and
is the pulse repetition period, determined by the distance to the target. The target echo signal can be represented as follows:
where
,
is the target distance, and
is the speed of light. The signal generated by agile frequency signal source 2 is as follows:
where one path, after mixing with signal
, yields the reference intermediate frequency (IF) signal
, while the other path, after mixing with the echo signal
, produces signal
:
After undergoing a series of processes including sampling, digital phase locking, quadrature digital down-conversion, and baseband signal post-processing, signals
and
yield the processed result
:
Each pulse corresponds to a result, , and the system emits 30 pulses per measurement. As a result, a series of 30 echo information samples was obtained. Performing an IDFT on the echo data yielded the target’s distance and position information, i.e., the one-dimensional image. Next, the region of interest was centered around the target position, with zero-padding applied before performing a DFT. The frequency domain data were then used to compute the power spectrum of the echo, which allowed for the calculation of the target’s RCS value.
According to the system specifications, the transmission power of the system needs to reach at least
to meet the testing requirements. Before introducing the long-distance signal transmission system, the system’s achieved transmission power was
. Therefore, as long as the loss introduced by the signal transmission system is within
, the requirements will be met. As indicated by the previous test results, the system not only did not undergo attenuation but also provided a gain of at least
. Thus, its integration would not affect the overall transmission power of the system. Moreover, the results shown in
Figure 19 indicate that the positional information in the one-dimensional image is only obscured when the introduced random phase error reaches around
. In the test results, the phase error of the signal transmission system is less than 3 degrees (
), which meets the required specifications.
Finally, after integrating the optoelectronic conversion system into the entire bistatic RCS testing system, we completed the measurement of the target’s RCS. Due to the lack of a physical vehicle or model, we selected a standard metal sphere as the test target. The metal sphere had a radius of approximately
meters, and the RCS value for a sphere of this size is around
. The testing environment is shown in
Figure 20.
The frequency chosen for the test was 10 GHz, with elevation angles of 0°, 45°, and 90° and horizontal angles ranging from −45° to 45°. The final measurement results are shown in
Figure 21. From the figure, it can be gleaned that the measurement results differ from the theoretical value of
by no more than 0.1 dB, indicating that the introduction of the optoelectronic conversion system did not affect the accuracy of the RCS measurement.