# Application of Compressive Sensing in Solving Monostatic Scattering Problems

Authors: Shui-Rong Chai; Li-Xin Guo; Long Li
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

In this paper, a new CS-FMM method that conjugates compressive sensing (CS) with the fast multipole method (FMM) is proposed and validated to efficiently solve monostatic scattering from an arbitrary conducting target. The far zone scattered fields are viewed as the signal of interest. CS is introduced to reduce the number of computations. A new set of incident sources has been generated according to CS. By solving the matrix equations under the new set of incident sources and calculating the related far zone scattered fields, the measurements of the aforementioned signal can be derived. Then, the CS inversion is employed to reconstruct the desired monostatic far zone scattered fields by finding the smallest possibleℓ1 norm solution. Monostatic radar cross section (RCS) from several conducting targets is studied by CS-FMM and by the traditional FMM. And the results are compared with each other to illustrate the accuracy and efficiency of the proposed method.


---

## Body

## 1. Introduction

Electromagnetic scattering from an arbitrary target has attracted a lot of attention over the past few decades for a large number of applications, such as target recognition and tracking, stealth target design, and radar surveillance. There exists many numerical methods that are proposed to analyze the electromagnetic scattering problems. For example, the method and moment (MoM) [1, 2], the shooting and bouncing ray (SBR) [3], and generalized method of moments [4, 5]. However, if monostatic scattering problem, which has wide incident angles, is interested, the aforementioned numerical methods need to obey the so-called Shannon-Nyquist sampling theorem to get a reasonable result. That is to say, the number of times to calculate the induced currents and the scattered fields is equal to the number of incident angles. This is time-consuming and typically leads to a huge amount of computational cost. Therefore, this is a very intuitive thinking that, in order to improve the efficiency of the aforementioned methods in solving the monostatic scattering problems while maintaining sufficient accuracy, one needs to reduce the sampling rate of the incident angles.There exists a new sampling paradigm named compressive sensing (CS), which goes against the Shannon-Nyquist sampling theorem. The CS promises to measure the informative part of the signal directly at a rate significantly below the Nyquist rate and accurately reconstruct the signal from its lower dimensional measurements. Obviously, CS meets the requirement of reducing the number of calculations in monostatic scattering analysis. Therefore, several studies have been published in the literature to exploit the ideas from CS in solving the wide-angle electromagnetic scattering problems.Carin et al. [6] are the first researchers that applied CS in solving wide incident angle scattering problems. He and his team of researchers at Duke University developed an in situ CS to reduce the number of computation required for target characterization by introducing the heterogeneous medium. D’Ambrosio [7] develops the ideas of Carin and integrates CS with the commercial software FEKO. Accordingly, D’Ambrosio’s method has both the benefits of CS and the deeply optimization benefits of FEKO. Chen et al. [8] and Chai and Guo [9, 10] attempt to combine the CS with the MoM for the purpose of efficiently analyzing two and three dimensional electromagnetic monostatic scattering problems. In [11], the triple hybrid approach that integrates the CS, the MoM, and the adaptive cross approximation (ACA) is proposed to further improve the efficiency of MoM in solving monostatic scattering problems. In [12] the hybrid method that conjugates CS with the multilevel fast multipole algorithm (MLFMA) is proposed, where the CS is utilized to reduce the number of incident angles. Wang et al. [13] propose a hybrid method that combines the CS with the characteristic basis function method (CBFM) for fast monostatic scattering analysis. The Bayesian compressive sensing together with the MoM is also utilized to solve monostatic scattering problems in [14].In this paper, the hybrid method that conjugates CS with the FMM is proposed to efficiently analyze monostatic scattering problems. Different from the aforementioned methods, the new CS-FMM method proposed in this paper has its own advantages. First, the induced currents (or the Fourier transforms of the induced currents as in [6, 7]) on one patch over all incident angles are viewed as the proceeded signal, while the signal represents the monostatic scattered fields in the proposed CS-FMM. Second, there exists only one signal in our method, but the number of signals is n in the methods mentioned above, with n denotes the number of unknowns. Third, the new incident source is a linear combination of all plane waves in all of the aforementioned methods, while the new incident source in the newly proposed CS-FMM is a plane wave that derived by randomly extracting from the original incident waves. Fourth, the Gaussian random matrix is chosen to be the measurement matrix in all of the methods mentioned above except [12]. And the measurement matrix is the uniform sampling matrix in CS-FMM, which has a more concise form than that of the Gaussian random matrix. Fifth, the time expenditure on reconstructing the signal of CS-FMM is much smaller than that of the previous methods.


## 2. Theoretical Model

### 2.1. Traditional Fast Multipole Method

The results and discussion may be presented separately, or in one combined section, and may optionally be divided into headed subsections.The purpose of solving a monostatic scattering problem is to evaluate the far zone monostatic scattered fields by(1)Escaθ=−ikηexpikρ04πρ0k^sca×k^sca×∫sexp−ikρ′⋅k^scaJρ′,θds′,where ρ0 is the distance between the observation point and the original point, ρ′ represents the position vector of the source point, k^sca is the unit vector of the scattered wave, θ denotes the incident and the scattering angle in the monostatic scattering analyzing, and J is the induced currents on the surface of the target, which can be described as
(2)Jρ,θ=∑i=1nαiθfiρ.To computeJ in (1) with traditional FMM, the following matrix equation can be derived:
(3)ATDIθ=Vθ,where A, T, and D denote the aggregation, transition, and disaggregation matrix, respectively. This three matrix can be constituted regardless of the choice of the incident angles, V is the known incident field vector, and I represents the unknown vector that formed by grouping the coefficients αiθ.By solving (3) under a particular incident angle θi, one can get the related unknown coefficients vector Iθi, and from Iθi the induced current Jρ,θi. Finally, the desired far zone scattered fields Escaθi can be computed by substituting the computed induced current Jρ,θi to (1).It should be pointed out that the discussion in the previous paragraph is the calculating process of computing the scattered fields under one incident angle. Obviously one needs to repeat the previous processN time, if the monostatic scattering at hand has N incident angles. This is time-consuming and leads to a huge amount of computational cost. Therefore, CS is introduced to improve the efficiency of the FMM in solving monostatic scattering problems by reducing the number of times required to solve (3), (2), and (1).


### 2.2. CS-FMM

To introduce CS into the monostatic scattering problems, the far zone scattered fields are viewed as the signal of interest, which can be described as a complex column vector of lengthN in ℂN. It is clear that there exists just one such signal.To facilitate the analysis, we rewrite (1) into the following matrix form:
(4)EN×1=diagGN×nJn×N,where E and J denote the vector (or matrix) form of the far zone scattered fields and the induced currents, respectively. G is the matrix that generated by gathering the rest part on the right of (1). Definitely, E is the interested signal in CS-FMM.CS attempts to integrate sensing and compression and therefore measure the informative part of the signal directly. Mathematically, the measurement process can be described by(5)sm×1=Φm×NEN×1=Φm×NdiagGN×nJn×N,where Φ is an m×N matrix called the measurement matrix, m denotes the number of measurements, and s is the vector that represents the samples of the signal.The Gaussian matrix, which is a fully populated matrix, is arranged as the measurement matrixΦ in all of the previous works except [12], where the measurement matrix has a more concise expression. However, both the Gaussian matrix and the measurement matrix in [12] make the elements of s a linear combination of the elements of the signal E. This makes it impossible to convert the measurement of the far zone scattered fields to the samples of the induced currents. Therefore, we introduce another type of measurement matrix named uniform sampling matrix, which is constructed as follows:
(6)Φm×N=01⋯00000⋯100⋮⋮⋱⋮⋮⋮10⋯000,where each row has only one non-zero element that equals to one, and the position of the non-zero element is drawn from a discrete uniform distribution over 1 to N. Each row in the uniform sampling matrix can be interpreted as the selection of one incident angle, and there are totally m incident angles that are selected.It can be found that, if the uniform sampling matrix is chosen as the measurement matrix,s can be generated by randomly constructing m elements of E. And (5) can be rewritten as
(7)sm×1=E~m×1=diagG~m×nJ~n×m,where G~ denotes the matrix that derived by randomly constituting m rows of G. J~ represents the induced currents on the surface of the target under the selected m incident angles, which can be derived by solving the associated (3) m times.Another point that should be emphasized is that the CS exploits the fact that the signal of interest is compressible itself or it is compressible in some sparse transform basis. In [15] we have examined that the electric field is compressible in FFT basis. Therefore, in the examples presented below, we choose fast Fourier transform (FFT) as the sparse transform basis.Finally, the CS inversion is employed to recover the original far zone scattered fields by solving the followingℓ1 norm optimization problem
(8)E=Ψy^;y^=argminy′1s.t.ΦΨy^=s,where Ψ is the sparse transform matrix, and y represents the corresponding coefficients vector. The orthogonal matching pursuit (OMP) [16] is applied to solve the ℓ1 norm optimization problem described by (8) in this paper.


## 2.1. Traditional Fast Multipole Method

The results and discussion may be presented separately, or in one combined section, and may optionally be divided into headed subsections.The purpose of solving a monostatic scattering problem is to evaluate the far zone monostatic scattered fields by(1)Escaθ=−ikηexpikρ04πρ0k^sca×k^sca×∫sexp−ikρ′⋅k^scaJρ′,θds′,where ρ0 is the distance between the observation point and the original point, ρ′ represents the position vector of the source point, k^sca is the unit vector of the scattered wave, θ denotes the incident and the scattering angle in the monostatic scattering analyzing, and J is the induced currents on the surface of the target, which can be described as
(2)Jρ,θ=∑i=1nαiθfiρ.To computeJ in (1) with traditional FMM, the following matrix equation can be derived:
(3)ATDIθ=Vθ,where A, T, and D denote the aggregation, transition, and disaggregation matrix, respectively. This three matrix can be constituted regardless of the choice of the incident angles, V is the known incident field vector, and I represents the unknown vector that formed by grouping the coefficients αiθ.By solving (3) under a particular incident angle θi, one can get the related unknown coefficients vector Iθi, and from Iθi the induced current Jρ,θi. Finally, the desired far zone scattered fields Escaθi can be computed by substituting the computed induced current Jρ,θi to (1).It should be pointed out that the discussion in the previous paragraph is the calculating process of computing the scattered fields under one incident angle. Obviously one needs to repeat the previous processN time, if the monostatic scattering at hand has N incident angles. This is time-consuming and leads to a huge amount of computational cost. Therefore, CS is introduced to improve the efficiency of the FMM in solving monostatic scattering problems by reducing the number of times required to solve (3), (2), and (1).


## 2.2. CS-FMM

To introduce CS into the monostatic scattering problems, the far zone scattered fields are viewed as the signal of interest, which can be described as a complex column vector of lengthN in ℂN. It is clear that there exists just one such signal.To facilitate the analysis, we rewrite (1) into the following matrix form:
(4)EN×1=diagGN×nJn×N,where E and J denote the vector (or matrix) form of the far zone scattered fields and the induced currents, respectively. G is the matrix that generated by gathering the rest part on the right of (1). Definitely, E is the interested signal in CS-FMM.CS attempts to integrate sensing and compression and therefore measure the informative part of the signal directly. Mathematically, the measurement process can be described by(5)sm×1=Φm×NEN×1=Φm×NdiagGN×nJn×N,where Φ is an m×N matrix called the measurement matrix, m denotes the number of measurements, and s is the vector that represents the samples of the signal.The Gaussian matrix, which is a fully populated matrix, is arranged as the measurement matrixΦ in all of the previous works except [12], where the measurement matrix has a more concise expression. However, both the Gaussian matrix and the measurement matrix in [12] make the elements of s a linear combination of the elements of the signal E. This makes it impossible to convert the measurement of the far zone scattered fields to the samples of the induced currents. Therefore, we introduce another type of measurement matrix named uniform sampling matrix, which is constructed as follows:
(6)Φm×N=01⋯00000⋯100⋮⋮⋱⋮⋮⋮10⋯000,where each row has only one non-zero element that equals to one, and the position of the non-zero element is drawn from a discrete uniform distribution over 1 to N. Each row in the uniform sampling matrix can be interpreted as the selection of one incident angle, and there are totally m incident angles that are selected.It can be found that, if the uniform sampling matrix is chosen as the measurement matrix,s can be generated by randomly constructing m elements of E. And (5) can be rewritten as
(7)sm×1=E~m×1=diagG~m×nJ~n×m,where G~ denotes the matrix that derived by randomly constituting m rows of G. J~ represents the induced currents on the surface of the target under the selected m incident angles, which can be derived by solving the associated (3) m times.Another point that should be emphasized is that the CS exploits the fact that the signal of interest is compressible itself or it is compressible in some sparse transform basis. In [15] we have examined that the electric field is compressible in FFT basis. Therefore, in the examples presented below, we choose fast Fourier transform (FFT) as the sparse transform basis.Finally, the CS inversion is employed to recover the original far zone scattered fields by solving the followingℓ1 norm optimization problem
(8)E=Ψy^;y^=argminy′1s.t.ΦΨy^=s,where Ψ is the sparse transform matrix, and y represents the corresponding coefficients vector. The orthogonal matching pursuit (OMP) [16] is applied to solve the ℓ1 norm optimization problem described by (8) in this paper.


## 3. Numerical Results and Discussion

To demonstrate the accuracy and efficiency of the proposed CS-FMM, the radar cross section (RCS) from several conducting objects is evaluated. And the results are compared with those derived from traditional FMM. In performing the calculations some parameters are given as follows: the operating frequency isf=300MHz, the incident waves are assumed to be within the xoz plane, and the elevation angle of the incident wave changes from 1∘ to 360∘, with an interval of 1∘. Accordingly, the traditional FMM has to solve (3), (2), and (1) 360 times to derive the wanted far zone scattered fields. The uniform sampling matrix and the FFT are chosen as the measurement matrix and sparse basis functions, respectively. All the simulations were performed on a computer with a 2.67 GHz processor (Intel Core i5 CPU) and 8.0 GB memory.


### 3.1. Cylinder

At first, the CS-FMM is applied to evaluate the monostatic RCS of a cylinder, which has a radius of 0.5 m and a length of 1.0 m. The cylinder axis is assumed to coincide with thex axis. During the simulation, the cylinder is dispersed into 1866 small triangles, which generate 2799 total unknowns, and the number of measurements is set as m=50. In other words, the number of times required to solve (3), (2), and (1) is reduced to 50 by the introduction of CS theory. RCS distribution derived by CS-FMM and traditional FMM is depicted in Figure 1. It can be found that the two methods obtain almost the same RCS results.Figure 1
Monostatic RCS of the cylinder alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


### 3.2. Square Plate

To further verify the accuracy and efficiency of the proposed CS-FMM, the RCS of a square plate with an edge length of 4 m is computed by CS-FMM and by FMM. The tackled plate is located on thexoy plane and is partitioned into 2918 small triangles. The number of unknowns and measurements is 4305 and 95, respectively. RCS results of the FMM and CS-FMM along xoz plane are demonstrated in Figure 2. As shown in Figure 2, the CS-FMM yields approximately accurate results.Figure 2
Monostatic RCS of the square plate alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


### 3.3. Missile

In this case, the proposed CS-FMM and the FMM are utilized to simulate the RCS distribution of a missile, which has a length of 6.7 m, width 0.52 m, and wingspan 3.18 m (as shown in Figure3). Comparisons of RCS distribution derived from FMM and CS-FMM are shown in Figure 4. During the simulation, the missile is dispersed into 2748 small triangles, which generate 4122 total unknowns. In performing the simulation with CS-FMM, the number of measurements is chosen as m=180. As can be observed in Figure 4, there is no visible difference between the two methods.Figure 3
Missile model.Figure 4
Monostatic RCS of the missile alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


### 3.4. Ship

Finally, the angular distribution of the monostatic RCS from a ship (as shown in Figure5) is simulated. The ship has a length of 9.8 m, width 1.2 m, and height 2.65 m. The number of unknowns in this case is 21129. And the number of measurements is set as m=180 in CS-FMM. Monostatic RCS results are presented in Figure 6. It is easy to be observed that the two methods agree very well with each other.Figure 5
Ship model.Figure 6
Monostatic RCS of the ship alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)The simulating time of the FMM and CS-FMM is summarized in Table1. It is observed that the CS-FMM can save significantly over traditionally FMM method.Table 1
Summary of simulating time of four models.
Model
Polarization
FMM (s)
CS-FMM (s)
Cylinder
HH
263
84
VV
210
77
Square plate
HH
420
232
VV
596
221
Missile
HH
747
437
VV
606
393
Ship
HH
3087
1642
VV
3728
2313To provide a quantitative comparison of the number of measurements on the accuracy of the proposed method, the relative root mean square error (R-RMSE) of the far zone scattered fields is defined as follows:(9)R‐RMSE=ECS‐FMM−EFMM2EFMM2,where ECS‐FMM and EFMM denote the far zone scattered fields derived by the proposed CS-FMM and the traditional FMM, respectively.In the following, several functions (FFT, discrete cosine transform (DCT), Legendre polynomials, and the first order Chebyshev polynomials) are utilized as the sparse basis functions in the proposed CS-FMM. The detailed descriptions of the Legendre polynomials and the first order Chebyshev polynomials can be found in [17].The justification for choosing FFT as the sparse basis functions in CS-FMM are presented in Figure7. In Figure 7, the R-RMSE as a function of number of measurements (changed from 10 to 100) for different types of sparse basis functions is depicted. The aforementioned cylinder is chosen as the tackled target. Analogously, the uniform sampling matrix is set as the measurement matrix. Considering the random quality of the measurement matrix, for a given number of measurements, R-RMSE results for 1000 experiments are calculated and averaged. From Figure 7, one can find that the average R-RMSE obtained by FFT basis function is smaller than the other three types of basis functions for both HH and VV polarization. This means that the FFT provides a more accurate result than the other three types of basis functions for a given number of measurements. In another word, the FFT basis function can provide a sparsest representation of the signal. This is the reason why FFT is chosen as the basis function in the cylinder case simulation. Meanwhile, m=50 is chosen in the cylinder case to balance simulating time and accuracy. Analogously, similar results can be derived for plate, missile, and ship cases.Figure 7
Average R-RMSE as a function of number of measurements for different types of sparse basis functions. (a) HH polarization; (b) VV polarization.
(a)
(b)


## 3.1. Cylinder

At first, the CS-FMM is applied to evaluate the monostatic RCS of a cylinder, which has a radius of 0.5 m and a length of 1.0 m. The cylinder axis is assumed to coincide with thex axis. During the simulation, the cylinder is dispersed into 1866 small triangles, which generate 2799 total unknowns, and the number of measurements is set as m=50. In other words, the number of times required to solve (3), (2), and (1) is reduced to 50 by the introduction of CS theory. RCS distribution derived by CS-FMM and traditional FMM is depicted in Figure 1. It can be found that the two methods obtain almost the same RCS results.Figure 1
Monostatic RCS of the cylinder alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


## 3.2. Square Plate

To further verify the accuracy and efficiency of the proposed CS-FMM, the RCS of a square plate with an edge length of 4 m is computed by CS-FMM and by FMM. The tackled plate is located on thexoy plane and is partitioned into 2918 small triangles. The number of unknowns and measurements is 4305 and 95, respectively. RCS results of the FMM and CS-FMM along xoz plane are demonstrated in Figure 2. As shown in Figure 2, the CS-FMM yields approximately accurate results.Figure 2
Monostatic RCS of the square plate alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


## 3.3. Missile

In this case, the proposed CS-FMM and the FMM are utilized to simulate the RCS distribution of a missile, which has a length of 6.7 m, width 0.52 m, and wingspan 3.18 m (as shown in Figure3). Comparisons of RCS distribution derived from FMM and CS-FMM are shown in Figure 4. During the simulation, the missile is dispersed into 2748 small triangles, which generate 4122 total unknowns. In performing the simulation with CS-FMM, the number of measurements is chosen as m=180. As can be observed in Figure 4, there is no visible difference between the two methods.Figure 3
Missile model.Figure 4
Monostatic RCS of the missile alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)


## 3.4. Ship

Finally, the angular distribution of the monostatic RCS from a ship (as shown in Figure5) is simulated. The ship has a length of 9.8 m, width 1.2 m, and height 2.65 m. The number of unknowns in this case is 21129. And the number of measurements is set as m=180 in CS-FMM. Monostatic RCS results are presented in Figure 6. It is easy to be observed that the two methods agree very well with each other.Figure 5
Ship model.Figure 6
Monostatic RCS of the ship alongxoz plane as functions of incident angles. (a) HH polarization; (b) VV polarization.
(a)
(b)The simulating time of the FMM and CS-FMM is summarized in Table1. It is observed that the CS-FMM can save significantly over traditionally FMM method.Table 1
Summary of simulating time of four models.
Model
Polarization
FMM (s)
CS-FMM (s)
Cylinder
HH
263
84
VV
210
77
Square plate
HH
420
232
VV
596
221
Missile
HH
747
437
VV
606
393
Ship
HH
3087
1642
VV
3728
2313To provide a quantitative comparison of the number of measurements on the accuracy of the proposed method, the relative root mean square error (R-RMSE) of the far zone scattered fields is defined as follows:(9)R‐RMSE=ECS‐FMM−EFMM2EFMM2,where ECS‐FMM and EFMM denote the far zone scattered fields derived by the proposed CS-FMM and the traditional FMM, respectively.In the following, several functions (FFT, discrete cosine transform (DCT), Legendre polynomials, and the first order Chebyshev polynomials) are utilized as the sparse basis functions in the proposed CS-FMM. The detailed descriptions of the Legendre polynomials and the first order Chebyshev polynomials can be found in [17].The justification for choosing FFT as the sparse basis functions in CS-FMM are presented in Figure7. In Figure 7, the R-RMSE as a function of number of measurements (changed from 10 to 100) for different types of sparse basis functions is depicted. The aforementioned cylinder is chosen as the tackled target. Analogously, the uniform sampling matrix is set as the measurement matrix. Considering the random quality of the measurement matrix, for a given number of measurements, R-RMSE results for 1000 experiments are calculated and averaged. From Figure 7, one can find that the average R-RMSE obtained by FFT basis function is smaller than the other three types of basis functions for both HH and VV polarization. This means that the FFT provides a more accurate result than the other three types of basis functions for a given number of measurements. In another word, the FFT basis function can provide a sparsest representation of the signal. This is the reason why FFT is chosen as the basis function in the cylinder case simulation. Meanwhile, m=50 is chosen in the cylinder case to balance simulating time and accuracy. Analogously, similar results can be derived for plate, missile, and ship cases.Figure 7
Average R-RMSE as a function of number of measurements for different types of sparse basis functions. (a) HH polarization; (b) VV polarization.
(a)
(b)


## 4. Conclusions

CS together with the FMM is introduced into monostatic scattering problems, and a new method named CS-FMM is proposed. The basic idea of CS-FMM is the sparse property of the far zone scattered fields. The CS is applied to exploit this sparsity and to reduce the number of computations. The FMM is worked as an electromagnetic solver. Numerical simulations show that the CS-FMM can improve the efficiency of traditional FMM dramatically in solving monostatic scattering problems.



---

*Source: 8595104-2019-01-17.xml*