# Incoherent Integration Detection Method of Airborne Phased Array Radar in a Multipath Environment

Authors: Yili Hu; Yongbo Zhao; Sheng Chen
Journal: International Journal of Antennas and Propagation
 (2021)

---

## Abstract

Airborne phased array radar (PAR) suffers from multipath problems when flying over a calm sea surface. The existence of a multipath phenomenon will cause the electromagnetic echo of the same target to be reflected back to the airborne PAR from two paths, namely, direct path (DP) and multipath. Compared with the ground-based radar, the target echo received by airborne PAR in the multipath environment has two important characteristics: one is that the DP signal and the multipath signal exist in different range bins, and the other is that the radar cross section (RCS) in the DP direction may be smaller than that in the multipath direction. Considering these two characteristics, this paper first proposes a target pairing algorithm for matching the DP range and multipath range of the same target in signal detection, and then, combined with the cell-averaging constant false alarm rate (CA-CFAR) detection model, an incoherent integration detection method for airborne PAR in the multipath environment is proposed. In the target pairing process, the geometric structure relationship of the airborne PAR model can be fully utilized. After a successful target pairing process, the energy of the multipath signal will be incoherently accumulated into the corresponding DP range bin, so as to improve the probability of DP range bin data passing the detection threshold. In essence, the proposed method makes full use of multipath energy to improve the detection capability of airborne PAR in the multipath environment. Finally, the detection probability of the proposed method is given, and the detection performance is analyzed.


---

## Body

## 1. Introduction

Airborne phased array radar (PAR) is a special radar system that places the antenna elements on the airborne platform. It is tactically important because of its wide field of vision, good manoeuvrability, and simultaneous detection capability of ground and air targets. Compared with a ground-based radar, the airborne PAR suffers from more complex clutter and interference, but when combined with the system characteristics of the airborne PAR, many solutions have been proposed to deal with these problems [1–6].Multipath is caused by the high reflectivity of the sea surface at microwave frequencies. Multipath causes multiple propagation paths between the radar and the target, which usually implies that the target is illuminated from two directions and that the radar receives the return energy from two directions [7–9]. Airborne PAR is prone to multipath phenomena in special reflection scenarios such as calm sea surfaces [10–12]. When airborne PAR suffers from multipath problems, most research studies are focused on multipath suppression of clutter and interference [13–16] and parameter estimation [17–20]. In the airborne PAR, the direct path (DP) signal and the multipath signal are not in the same range bin. Therefore, the data in the DP range bin can be directly used for signal detection [21–23]. Obviously, these methods result in a waste of multipath signal energy and cannot fully utilize the targets’ multipath energy for signal detection. In [22], the radar detection problem is devised as a constrained generalized likelihood ratio test in a multipath diffuse reflection environment, and the proposed method still has the property of the constant false alarm rate (CFAR). In [23], an adaptive constrained generalized likelihood ratio test method is proposed to find range diversity targets under diffuse reflection, and the main process is that in each range bin of the original data, the echo data is transformed into the covariance of the secondary data sample, but the detection performance of the method in heterogeneous echo data needs to be improved. In general, the methods in [22, 23] do not make full use of the energy of multipath, and the radar cross section (RCS) characteristics of different path echoes of airborne PAR in a multipath environment are not analyzed.In this paper, the synthetic multipath reflection coefficient (SMRC) of airborne PAR is analyzed. Based on the analysis of SMRC, this paper systematically summarizes the power relationship between the DP signal and the multipath signal of airborne PAR and clearly emphasizes a special multipath phenomenon; that is, the multipath power of airborne PAR in the multipath environment may be greater than the DP power. Furthermore, aiming at the multipath model of airborne PAR, this paper proposes a new signal detection process to improve the detection probability by incoherent integration of multipath energy. Specifically, a target pairing algorithm for matching the DP range and multipath range of the same target is proposed first, and then the corresponding incoherent integration detection method is carried out after a successful target pairing process. Finally, the detection probability expression of the proposed incoherent integration detection method is given.The remainder of this paper is organized as follows: the analysis of the multipath model and the SMRC of airborne PAR in a multipath environment are given in Section2; the proposed incoherent integration detection method is described in Section 3; and Section 4 analyzes the performance of the proposed method. Then, summarizing remarks are given in Section 5.


## 2. Signal Model Analysis of Airborne Phased Array Radar in a Multipath Environment

### 2.1. Analysis of Multipath Model

The echo model of the airborne PAR in a multipath environment is shown in Figure1.Figure 1
The signal model of airborne PAR in a multipath environment.In Figure1, Re is the radius of the earth, H is the platform height of airborne PAR, Ht is the height of the target, θd is the angle between the horizontal and the received DP, θi is the angle between the horizontal and the received multipath, Rd is the DP range, Ri=R1+R2 is the multipath range, and hb is the height of the reflection area. Note that the multipath phenomenon of airborne PAR is mainly caused by strong specular reflection when flying over a calm sea surface, so only the multipath of the elevation direction needs to be considered [24].According to the echo characteristics of the airborne radar, the range difference between the DP signal and the multipath signal is very large, causing the DP signal and multipath signal to not be within the same range bin. Thus, from Figure1, the angle θB between the DP signal and the multipath signal is large, so that the RCS of the DP signal is possible to be smaller than that of the multipath signal, resulting in the echo power of the multipath signal being greater than that of the DP signal. If the power of the multipath signal is not used in signal detection, the multipath return energy will be wasted.Note that in the ground-based radar, the angleθB is very small (usually less than 0.5°), so it is often considered that the RCS values of the DP signal and the multipath signal in the ground-based radar model are approximately equal. The relationship between the RCS values of the DP signal and the multipath signal is an important difference of echo energy between airborne PAR and ground-based radar in the multipath environment.The simplified signal model of Figure1 is shown in Figure 2.Figure 2
The simplified signal model in a multipath environment.From Figure2, the multipath range Ri is greater than the DP range Rd. The range difference between Ri and Rd is a function of the radar’s waveform parameters and the geometry model of the reflection surface. Generally, as long as the platform height H of the airborne PAR model is not very low (i.e., H≥1km), Rd and Ri are in different range bins. Therefore, this paper defaults that the DP range and multipath range are in different range bins.


### 2.2. The Modeling and Analysis of the Synthetic Multipath Reflection Coefficient

Due to the large angleθB of the airborne PAR multipath model shown in Figure 2, the RCS values of the DP signal and the multipath signal are no longer approximately equal. Therefore, it is important to study the echo power of the DP signal and the multipath signal.To analyze the power of signals from different paths of the same target in a multipath environment, the concept of the synthetic multipath reflection coefficient (SMRC) will be introduced. The SMRC reflects the amplitude ratio of the multipath signal to the DP signal.The echo powerPo_d of the DP signal and the echo power Po_i of the multipath signal can be expressed as follows:(1)Po_d=Pt,dGt,dGr,dλ2σd4π3Rd4,(2)Po_i=Pt,iGt,iGr,iλ2σi4π3Ri4,where λ is wavelength, σd and σi are the RCS of the DP signal and the multipath signal, respectively, transmitting power Pt,d=Pt,i, transmitting gain Gt,d=Gt,i, and the reflection attenuation coefficient ε≤1, that is, Gr,i=εGr,d. The SMRC ρ can be represented as follows:(3)ρ2=Po_iPo_d=εσiRd4σdRi4.According to equation (3), the SMRC ρ is related to the RCS, reflection attenuation coefficient ε, and ranges of the DP signal and the multipath signal.In a ground-based radar, the multipath phenomenon is easy to occur when the ground-based radar measures the height of a low-elevation target [25, 26], where σi/σd and Rd4/Ri4 are approximately equal to 1. Therefore, according to equation (3), the SMRC of the ground-based radar is mainly determined by the reflection attenuation coefficient ε. However, in the airborne PAR multipath model, σi/σd and Rd4/Ri4 are no longer approximately equal to 1, and their specific ratio values should be considered in combination with the actual radar system parameters. Note that according to equation (3), the SMRC amplitude ρ may be greater than 1, which means that compared with the ground-based radar, airborne PAR has a special multipath phenomenon; that is, multipath power may be greater than the DP power.In the process of signal detection, the DP rangeRd and the multipath range Ri can be obtained by pulse compression (PC) operation. Through the range difference ΔR=Ri−Rd between Ri and Rd, the relative time delay difference Δτ can be expressed as follows:(4)Δτ=τii−τdd=Ri−Rdc,where τdd=2Rd/c, τii=Rd+Ri/c, and c=3×108m/s. The large range difference ΔR between the DP signal and the multipath signal results in large relative time delay Δτ, that is, Δτ>1/2fs, where fs is the sampling frequency. Therefore, the complex envelopes of the DP signal and the multipath signal will be misalignment, leading to the failure of the joint steering vector class method [27] for signal detection.


## 2.1. Analysis of Multipath Model

The echo model of the airborne PAR in a multipath environment is shown in Figure1.Figure 1
The signal model of airborne PAR in a multipath environment.In Figure1, Re is the radius of the earth, H is the platform height of airborne PAR, Ht is the height of the target, θd is the angle between the horizontal and the received DP, θi is the angle between the horizontal and the received multipath, Rd is the DP range, Ri=R1+R2 is the multipath range, and hb is the height of the reflection area. Note that the multipath phenomenon of airborne PAR is mainly caused by strong specular reflection when flying over a calm sea surface, so only the multipath of the elevation direction needs to be considered [24].According to the echo characteristics of the airborne radar, the range difference between the DP signal and the multipath signal is very large, causing the DP signal and multipath signal to not be within the same range bin. Thus, from Figure1, the angle θB between the DP signal and the multipath signal is large, so that the RCS of the DP signal is possible to be smaller than that of the multipath signal, resulting in the echo power of the multipath signal being greater than that of the DP signal. If the power of the multipath signal is not used in signal detection, the multipath return energy will be wasted.Note that in the ground-based radar, the angleθB is very small (usually less than 0.5°), so it is often considered that the RCS values of the DP signal and the multipath signal in the ground-based radar model are approximately equal. The relationship between the RCS values of the DP signal and the multipath signal is an important difference of echo energy between airborne PAR and ground-based radar in the multipath environment.The simplified signal model of Figure1 is shown in Figure 2.Figure 2
The simplified signal model in a multipath environment.From Figure2, the multipath range Ri is greater than the DP range Rd. The range difference between Ri and Rd is a function of the radar’s waveform parameters and the geometry model of the reflection surface. Generally, as long as the platform height H of the airborne PAR model is not very low (i.e., H≥1km), Rd and Ri are in different range bins. Therefore, this paper defaults that the DP range and multipath range are in different range bins.


## 2.2. The Modeling and Analysis of the Synthetic Multipath Reflection Coefficient

Due to the large angleθB of the airborne PAR multipath model shown in Figure 2, the RCS values of the DP signal and the multipath signal are no longer approximately equal. Therefore, it is important to study the echo power of the DP signal and the multipath signal.To analyze the power of signals from different paths of the same target in a multipath environment, the concept of the synthetic multipath reflection coefficient (SMRC) will be introduced. The SMRC reflects the amplitude ratio of the multipath signal to the DP signal.The echo powerPo_d of the DP signal and the echo power Po_i of the multipath signal can be expressed as follows:(1)Po_d=Pt,dGt,dGr,dλ2σd4π3Rd4,(2)Po_i=Pt,iGt,iGr,iλ2σi4π3Ri4,where λ is wavelength, σd and σi are the RCS of the DP signal and the multipath signal, respectively, transmitting power Pt,d=Pt,i, transmitting gain Gt,d=Gt,i, and the reflection attenuation coefficient ε≤1, that is, Gr,i=εGr,d. The SMRC ρ can be represented as follows:(3)ρ2=Po_iPo_d=εσiRd4σdRi4.According to equation (3), the SMRC ρ is related to the RCS, reflection attenuation coefficient ε, and ranges of the DP signal and the multipath signal.In a ground-based radar, the multipath phenomenon is easy to occur when the ground-based radar measures the height of a low-elevation target [25, 26], where σi/σd and Rd4/Ri4 are approximately equal to 1. Therefore, according to equation (3), the SMRC of the ground-based radar is mainly determined by the reflection attenuation coefficient ε. However, in the airborne PAR multipath model, σi/σd and Rd4/Ri4 are no longer approximately equal to 1, and their specific ratio values should be considered in combination with the actual radar system parameters. Note that according to equation (3), the SMRC amplitude ρ may be greater than 1, which means that compared with the ground-based radar, airborne PAR has a special multipath phenomenon; that is, multipath power may be greater than the DP power.In the process of signal detection, the DP rangeRd and the multipath range Ri can be obtained by pulse compression (PC) operation. Through the range difference ΔR=Ri−Rd between Ri and Rd, the relative time delay difference Δτ can be expressed as follows:(4)Δτ=τii−τdd=Ri−Rdc,where τdd=2Rd/c, τii=Rd+Ri/c, and c=3×108m/s. The large range difference ΔR between the DP signal and the multipath signal results in large relative time delay Δτ, that is, Δτ>1/2fs, where fs is the sampling frequency. Therefore, the complex envelopes of the DP signal and the multipath signal will be misalignment, leading to the failure of the joint steering vector class method [27] for signal detection.


## 3. The Proposed Incoherent Integration Detection Method

### 3.1. The Proposed Method

The CFAR is a common signal detection method. In this section, the cell-average CFAR (CA-CFAR) [28] detection model will be used. Figure 3 shows a flowchart of the proposed incoherent integration detection method for airborne PAR in a multipath environment.Figure 3
The proposed incoherent integration detection model.In Figure3, system parameters mainly refer to the platform height H and the elevation wave position in the fast time data process when the radar transmits a burst of pulses. Sampling echo data refers to the discretized of echo data sampled by the receiving antennas after passing through A/D components. The MTD refers to the moving target detection, and it is also a classic Doppler filtering operation [29], which can effectively suppress interference and sea clutter. The PC refers to the pulse compression which can obtain the target ranges. The process of target pairing is to ensure that the ranges of the DP signal and the multipath signal from the same target can be paired. Note that, in the process of data from PC to target pairing, the operation of plot centroid [24] is required to ensure that the targets are roughly found. However, when the noise amplitude is much larger than the signal amplitude, the rough found target cannot be obtained by the plot centroid operation. In order to avoid the noise becoming the false target and affecting the success probability of signal detection, the proposed method in Figure 3 is not suitable for application at a very low signal-to-noise ratio (SNR).The processes of target pairing and incoherent integration are the most important differences between the proposed method in Figure3 and existing methods [13–16, 21–23]. Target pairing is an important process in Figure 3. The DP range Rd and the multipath range Ri can be obtained after PC, but when we are detecting targets, we cannot directly confirm that Rd and Ri correspond to the two different path ranges of the same target. Therefore, after PC, we need to perform the target pairing process to match the DP range and the multipath range of the same target.In the airborne PAR system, some priori information can be used in the target pairing process, especially the platform heightH and the elevation wave position θ˜d.The specific target pairing process is as follows: firstly, the system parametersH and θ˜d are obtained, and then, the range set of the rough found targets R=R1,…,Rk,…,RK can be obtained after PC, where k=1,2,…,K, K is the number of rough found targets. Finally, according to the system parameters and one of the range value Rk∈R, the multipath range Rk2 can be calculated by using the geometric structure relationship in Figure 1.Furthermore, the specific calculation process of the multipath rangeRk2 can be as follows: the revised platform height is H¯=H−hb, the revised radius of the earth is R¯e=Re+hb, and then, assuming that the DP range is Rk in the range set R for target pairing, the target height Ht can be represented as follows:(5)Ht=Rk2+R¯e+H¯2−2RkR¯e+H¯cosθ˜d+π2−R¯e.The angleϕ in Figure 1 can be obtained as follows:(6)ϕ=cos−1R¯e+H¯2+R¯e+Ht2−Rk22R¯e+H¯R¯e+Ht.The lengthG of the curve in Figure 1 can be represented as follows:(7)G=R¯eϕ.And then,ϕ1 and ϕ2 can be represented as follows:(8)ϕ1=G1R¯e,(9)ϕ2=G−G1R¯e,where(10)G1=G2−ηsinξ3,(11)η=23R¯eH¯+Ht+G24,(12)ξ=sin−12R¯eGHt−H¯η3.Finally, the multipath rangeRk2 can be represented as follows:(13)Rk2=Rk_1+Rk_2,Rk_1=R¯e2+R¯e+H¯2−2R¯eR¯e+H¯cosϕ1,Rk_2=R¯e2+R¯e+Ht2−2R¯eR¯e+Htcosϕ2.IfRk2 is also in the target range set R, it can be considered that the DP range Rk and the multipath range Rk2 are successfully paired. If Rk2 is not in the target range set R, it means that the target at the range Rk has no multipath echo or the multipath echo energy is too weak to be detected.Once the DP range of the target is successfully paired with its multipath range, we can perform the incoherent integration process. If the snapshot data of the DP signal and the multipath signal areXd and Xi, respectively, then the snapshot data Xdi at the DP range bin after incoherent integration is as follows:(14)Xdi2=Xd2+Xi2.Thus, the power of the DP range bin after incoherent integration isPo_di=1+ρ2Po_d, where Po_d is the power of the DP signal before incoherent integration.Note that if the platform height is very low and the airborne PAR system is oversampling, the DP range and multipath range may also be in the same range bin. But the proposed incoherent integration method can still use equations (5)–(13) to calculate the multipath range in the target pairing process. If the calculated range difference between the multipath range and the DP range is within a range bin, the incoherent integration of equation (14) is still valid, and the proposed method can still ensure that the multipath echo energy can be fully utilized for signal detection.


### 3.2. Analysis of Detection Probability of the Proposed Method

The decision criteria of CA-CFAR is as follows:(15)D>TZ,H1,D<TZ,H0,where H1 indicates that the target is detected, H0 means that no target is detected, D stands for detection statistics value, T is the normalization factor, and Z represents the estimated power of interference and noise. It is assumed that only Gaussian white noise exists in the echo data of airborne PAR except DP signals and multipath signals, and interference and sea clutter are suppressed during the MTD operation in Figure 1.In the process of CA-CFAR [28], the relationship between false alarm probability Pfa and normalization factor T is as follows:(16)Pfa=∫0∞e−s/μfSsds=∫0∞e−s/μμ−NsN−1e−s/μds=1+TN−N.Therefore, when thePfa is fixed, the normalization factor T can be represented as follows:(17)T=NPfa−1/N−1,where N=2n is the number of reference cell.The electromagnetic echo signal can be regarded as a sinusoidal sequence [24]. The detection probability expression Pd of the proposed incoherent integration detection method can be represented as follows:(18)Pd=∫S′∞12μzλ′n′−2/4e−z+λ′/2μIn′/2−1zλ′μdz=QM2λ′,2S′+e−S′+λ′S′λ′1/2I12λ′S′,where S′=T/N∑i=1NXi≈Tμ is the detection threshold, μ is the background noise power, n′=2N, λ′=Po_di/μ=1+ρ2Po_d/μ2, QM⋅ is the Q function of Marcum [28], and I1⋅ is the modified Bessel function of first order. The derivation of equation (18) is detailed in the Appendix.


## 3.1. The Proposed Method

The CFAR is a common signal detection method. In this section, the cell-average CFAR (CA-CFAR) [28] detection model will be used. Figure 3 shows a flowchart of the proposed incoherent integration detection method for airborne PAR in a multipath environment.Figure 3
The proposed incoherent integration detection model.In Figure3, system parameters mainly refer to the platform height H and the elevation wave position in the fast time data process when the radar transmits a burst of pulses. Sampling echo data refers to the discretized of echo data sampled by the receiving antennas after passing through A/D components. The MTD refers to the moving target detection, and it is also a classic Doppler filtering operation [29], which can effectively suppress interference and sea clutter. The PC refers to the pulse compression which can obtain the target ranges. The process of target pairing is to ensure that the ranges of the DP signal and the multipath signal from the same target can be paired. Note that, in the process of data from PC to target pairing, the operation of plot centroid [24] is required to ensure that the targets are roughly found. However, when the noise amplitude is much larger than the signal amplitude, the rough found target cannot be obtained by the plot centroid operation. In order to avoid the noise becoming the false target and affecting the success probability of signal detection, the proposed method in Figure 3 is not suitable for application at a very low signal-to-noise ratio (SNR).The processes of target pairing and incoherent integration are the most important differences between the proposed method in Figure3 and existing methods [13–16, 21–23]. Target pairing is an important process in Figure 3. The DP range Rd and the multipath range Ri can be obtained after PC, but when we are detecting targets, we cannot directly confirm that Rd and Ri correspond to the two different path ranges of the same target. Therefore, after PC, we need to perform the target pairing process to match the DP range and the multipath range of the same target.In the airborne PAR system, some priori information can be used in the target pairing process, especially the platform heightH and the elevation wave position θ˜d.The specific target pairing process is as follows: firstly, the system parametersH and θ˜d are obtained, and then, the range set of the rough found targets R=R1,…,Rk,…,RK can be obtained after PC, where k=1,2,…,K, K is the number of rough found targets. Finally, according to the system parameters and one of the range value Rk∈R, the multipath range Rk2 can be calculated by using the geometric structure relationship in Figure 1.Furthermore, the specific calculation process of the multipath rangeRk2 can be as follows: the revised platform height is H¯=H−hb, the revised radius of the earth is R¯e=Re+hb, and then, assuming that the DP range is Rk in the range set R for target pairing, the target height Ht can be represented as follows:(5)Ht=Rk2+R¯e+H¯2−2RkR¯e+H¯cosθ˜d+π2−R¯e.The angleϕ in Figure 1 can be obtained as follows:(6)ϕ=cos−1R¯e+H¯2+R¯e+Ht2−Rk22R¯e+H¯R¯e+Ht.The lengthG of the curve in Figure 1 can be represented as follows:(7)G=R¯eϕ.And then,ϕ1 and ϕ2 can be represented as follows:(8)ϕ1=G1R¯e,(9)ϕ2=G−G1R¯e,where(10)G1=G2−ηsinξ3,(11)η=23R¯eH¯+Ht+G24,(12)ξ=sin−12R¯eGHt−H¯η3.Finally, the multipath rangeRk2 can be represented as follows:(13)Rk2=Rk_1+Rk_2,Rk_1=R¯e2+R¯e+H¯2−2R¯eR¯e+H¯cosϕ1,Rk_2=R¯e2+R¯e+Ht2−2R¯eR¯e+Htcosϕ2.IfRk2 is also in the target range set R, it can be considered that the DP range Rk and the multipath range Rk2 are successfully paired. If Rk2 is not in the target range set R, it means that the target at the range Rk has no multipath echo or the multipath echo energy is too weak to be detected.Once the DP range of the target is successfully paired with its multipath range, we can perform the incoherent integration process. If the snapshot data of the DP signal and the multipath signal areXd and Xi, respectively, then the snapshot data Xdi at the DP range bin after incoherent integration is as follows:(14)Xdi2=Xd2+Xi2.Thus, the power of the DP range bin after incoherent integration isPo_di=1+ρ2Po_d, where Po_d is the power of the DP signal before incoherent integration.Note that if the platform height is very low and the airborne PAR system is oversampling, the DP range and multipath range may also be in the same range bin. But the proposed incoherent integration method can still use equations (5)–(13) to calculate the multipath range in the target pairing process. If the calculated range difference between the multipath range and the DP range is within a range bin, the incoherent integration of equation (14) is still valid, and the proposed method can still ensure that the multipath echo energy can be fully utilized for signal detection.


## 3.2. Analysis of Detection Probability of the Proposed Method

The decision criteria of CA-CFAR is as follows:(15)D>TZ,H1,D<TZ,H0,where H1 indicates that the target is detected, H0 means that no target is detected, D stands for detection statistics value, T is the normalization factor, and Z represents the estimated power of interference and noise. It is assumed that only Gaussian white noise exists in the echo data of airborne PAR except DP signals and multipath signals, and interference and sea clutter are suppressed during the MTD operation in Figure 1.In the process of CA-CFAR [28], the relationship between false alarm probability Pfa and normalization factor T is as follows:(16)Pfa=∫0∞e−s/μfSsds=∫0∞e−s/μμ−NsN−1e−s/μds=1+TN−N.Therefore, when thePfa is fixed, the normalization factor T can be represented as follows:(17)T=NPfa−1/N−1,where N=2n is the number of reference cell.The electromagnetic echo signal can be regarded as a sinusoidal sequence [24]. The detection probability expression Pd of the proposed incoherent integration detection method can be represented as follows:(18)Pd=∫S′∞12μzλ′n′−2/4e−z+λ′/2μIn′/2−1zλ′μdz=QM2λ′,2S′+e−S′+λ′S′λ′1/2I12λ′S′,where S′=T/N∑i=1NXi≈Tμ is the detection threshold, μ is the background noise power, n′=2N, λ′=Po_di/μ=1+ρ2Po_d/μ2, QM⋅ is the Q function of Marcum [28], and I1⋅ is the modified Bessel function of first order. The derivation of equation (18) is detailed in the Appendix.


## 4. Simulation Results

The SMRC characteristics of the DP signal and the multipath signal are analyzed, and the effectiveness of the proposed incoherent integration detection method is verified. Finally, the influence of the SMRCρ and false alarm probability Pfa on detection probability is discussed, respectively.The transmitting and receiving of airborne PAR systems are in phased array mode. For convenience, the amplitude-phase error and mutual coupling effect of the airborne PAR system are not considered. All the results in this section are simulated by MATLAB software.Table1 shows the main system parameters of airborne PAR.Table 1
System parameters.
DescriptionSymbolValueWavelengthλ0.03 mNumber of receiving antennasM15Platform heightH1.5 kmReflection attenuation coefficientε1Earth radiusRe6375 kmDirect path directionθd50°Multipath directionθi−54.11°Sampling frequencyfs2 MHzIn Table1, the receiving antenna M refers to the number of antennas in the elevation direction. Since the reflection area is on the calm sea surface, only the strong specular multipath phenomenon needs to be considered, and the height of reflection area is hb=20m. According to Table 1, a range bin represents 75 m.Simulation 1: the DP rangeRd belongs to 10,40km, and Figure 4 shows the relationship between the angle θB and the DP range Rd.Figure 4
The relationship between the angleθB and the DP range Rd.In Figure4, when DP range belongs to 10,40km and H is 1 km, 1.5 km, 2 km, and 2.5 km, respectively, the minimum range bin differences between the multipath signal and the DP signal are 16, 24, 32, and 40, respectively.Taking the ground-based radar as a comparative experiment, it is assumed that the platform height of the ground-based radar is only 50 m, and other parameters of the ground-based radar are the same as those in Table1. When the DP range belongs to 10,40km, angle θB of the ground-based radar does not exceed 0.22 degree and the maximum range bin difference between the multipath signal and the DP signal is less than 1. Obviously, the angle θB of the airborne PAR is much larger than that of the ground-based radar, and the multipath signal and the DP signal of the airborne PAR are in different range bins.According to the measured RCS data which has been provided in [24], it can be seen that the SMRC amplitude ρ is possible to be more than 1 at such a large angle θB of airborne PAR, which means that the SMRC characteristics introduced previously are in line with the actual situation. Therefore, it is important to utilize the energy of multipath signals for signal detection, which is often ignored in current detection methods [21–23].Simulation 2: since the reflection surface heighthb is an estimated value based on the topography of the reflection surface, the estimated reflection surface height will inevitably lead to a deviation of multipath range estimation in the target pairing process. Assume that the actual height and the calculated height of the reflection surface are h¯b=20m and h˜b, respectively. ΔR=R˜MP−R¯MP, Δh=h˜b−h¯b, where R¯MP is the corresponding multipath range in the condition of h¯b=20m, and R˜MP is the corresponding multipath range when the height of the reflection surface is h˜b. Figure 5 simulates the multipath range deviation when the height deviation between the calculated reflection surface and the actual reflection surface changes.Figure 5
The relationship between multipath range deviation and reflection surface deviation.In Figure5, as Δh increases, ΔR will also increase. But in the case of the same Δh, ΔR will gently decrease with the increase of Rd. It can be seen from Figure 5 that ΔR is always less than one range bin when Δh≤40m. For the multipath model with the reflection area on the calm sea surface, Δh is difficult to exceed 40m, indicating that the reflection surface height deviation has a limited influence on the estimation of multipath range bin.Simulation 3: the airborne PAR transmits LFM signal, the transmitting time width is40us, the bandwidth is 1.5MHz, and the pulse repetition frequency is fr=2000Hz. Assuming that there are three targets with DP ranges of 24.25 km, 25 km, and 25.9 km, respectively, we set the SNR after PC to SNR=12dB, the number of guard cells taken is 1, and the reference cell as follows: N=16, Pfa=10−6, and ρ=−1.2. Figure 6 compares the output gain of conventional CA-CFAR with that of the proposed incoherent integration detection method.Figure 6
Comparison of output gain amplitude. (a) The gain of CA-CFAR. (b) The gain of the proposed incoherent integration detection method.
(a)(b)In Figure6, the DP range bins of three targets are at the 323th, 333th, and 345th sampling data, and the method in Figure 6(b) is the proposed method. The target output gains of CA-CFAR are 32.84 dB, 31.67 dB, and 30 dB, respectively, and the target output gains of the proposed method are 36.19 dB, 36.46 dB, and 35.57 dB, respectively. From Figure 6, the target output gains of the proposed incoherent integration method are larger than the CA-CFAR, and the proposed method can successfully pass the detection threshold at the DP range bins, but the CA-CFAR method cannot. In fact, according to the previous analysis, the proposed method makes use of the return energy of multipath data, so the proposed method is more likely to pass the detection threshold.Simulation 4: the CA-CFAR and GLRT [28, 30] methods are used as comparison results. Supposing Rd=25km and σi/σd∈0.8,4, we can calculate the SMRC belongs to ρ∈−1.67,−0.75. Figure 7(a) shows the influence of the SMRC ρ=−1.67,−0.75 when Pfa=10−6, and Figure 7(b) shows the influence of false alarm probability Pfa=10−6,10−7,10−8 when ρ=−1.2.Figure 7
Comparison of detection probability. (a) Detection probability of differentρ, when Pfa=10−6. (b) Detection probability of different Pfa, when ρ=−1.2.
(a)(b)In Figure7, the CA-CFAR and GLRT methods directly detect the signal of the DP range bin. From Figures 7(a) and 7(b), it is obvious that the proposed incoherent integration detection method makes full use of the energy of the multipath signal, so it has good detection probability performance.In summary, Figure4 demonstrates the special multipath phenomenon of airborne PAR, that is, multipath power of airborne PAR may be greater than the DP power, Figure 5 analyzes the influence of the reflection surface height deviation on the estimation multipath range deviation, Figure 6 demonstrates that the proposed method can work with multiple targets, and Figure 7 shows the detection probability performance of the proposed method.


## 5. Conclusions

This paper analyzes the SMRC characteristics of airborne PAR in a multipath environment. Based on these characteristics, an incoherent integration detection method is proposed, which makes full use of the energy of multipath signals. At the same time, the detection probability of incoherent integration in the multipath environment is given. The proposed detection process and simulation conditions are in line with the practical application. Finally, simulation results show that the proposed method can improve the detection performance.However, accurate prior information is required in the target pairing process. When the prior information is inaccurate, the success probability of the target pairing process will decrease. Therefore, solving the dependence of the target pairing process on prior information is the focus of future research.



---

*Source: 8699512-2021-12-15.xml*