# A Fast Algorithm for Electromagnetic Scattering from One-Dimensional Rough Surface

Authors: Jing Jing Wang; An Qi Wang; Zhi Xiang Huang; Tie Zhen Jiang
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

In this paper, the Adaptive Modified Characteristic Basis Function Method (AMCBFM) is proposed for quickly simulating electromagnetic scattering from a one-dimensional perfectly electric conductor (PEC) rough surface. Similar to the traditional characteristic basis function method (CBFM), Foldy-Lax multiple scattering equations are applied in order to construct the characteristic basis functions (CBFs). However, the CBFs of the AMCBFM are different from those of the CBFM. In the AMCBFM, the coefficients of the CBFs are first defined. Then, the coefficients and the CBFs are used to structure the total current, which is used to represent the induced current along the rough surface. Moreover, a current criterion is defined to adaptively halt the order of the CBFs. The validity and efficiency of the AMCBFM are assessed by comparing the numerical results of the AMCBFM with the method of moments (MoM). The AMCBFM can effectively reduce the size of the matrix, and it costs less than half the CPU time used by the MoM. Moreover, by comparing it with the traditional CBFM, the AMCBFM can guarantee the accuracy, reduce the number of iterations, and achieve better convergence performance than the CBFM does. The second order of the CBFs is set in the CBFM. Additionally, the first order of the CBFs of the AMCBFM alone is sufficient for this result.


---

## Body

## 1. Introduction

Electromagnetic (EM) scattering from rough surface has been widely studied and applied in the research areas of marine communication, target detection, and stealth technology [1–5]. The method of moments (MoM) [6–8] has high accuracy and is widely used to simulate the EM scattering from rough surface. However, the dimensions of the solution matrix of the MoM are N×N, where N represents the number of unknowns. When solving this problem using the MoM, one has to simulate a sufficiently long rough surface, which results in many unknowns N and time-consumption. Several numerical methods have been developed, such as the fast multipole method (FMM) [9] and the characteristic basis function method (CBFM) [10–12]. In the FMM, an additional algorithm is used to deal with the kernel function. In addition, the iterative method is applied to solve the matrix equation. However, the FMM is an iterative method, which is limited by the convergence of the matrix equation. In the CBFM, the Foldy-Lax multipath scattering equations [13] are applied in order to structure the characteristic basis functions (CBFs). The self-interaction is first considered in order to generate the primary characteristic basis functions (PCBFs). Moreover, the mutual coupling effect of the subcells is considered in order to build the secondary characteristic basis functions (SCBFs). In previous work of our teams [14], the CBFM is used to simulate the EM scattering from a dielectric rough surface. In this work, the orders of the SCBFs were set by us and were determined through several tries. In the previous work of others [15], the following criterion was used to determine the order of the SCBFs:(1)esp=V-Z·Jt2V2×100%where Z is the impedance matrix, V is the known vector that is excited by the incident wave, Jt represents the induced current that is obtained by the CBFM, and ·2 is the L2 norm. From (1), it is readily seen that Z must be refilled and the matrix vector product is calculated in order to determine the order of the SCBFs. Therefore, determining the order of the SCBFs using (1) will be time-consuming and require more memory.In this paper, the Adaptive Modified Characteristic Basis Function Method (AMCBFM) [16] with a new current criterion is proposed. Similar to the traditional CBFM, the Foldy-Lax multipath scattering equations are also used to structure the CBFs. The PCBFs, such as the CBFM, are first obtained. However, the coefficients of the PCBFs are defined. Then, the primary total current of each cell is obtained by combining the coefficients and the PCBFs. Then, the coefficients are taken with the PCBFs in order to get the first order of the SCBFs. The coefficients of the first order of the SCBFs and the first total current of each cell are obtained. Like the primary total current, the first total current is used to get the second order of the SCBFs. Higher orders of the SCBFs are similarly obtained. Moreover, the new current criterion that is defined in this paper is used to adaptively halt the order of the SCBFs.The remainder of this paper is organized as follows. In Section2, the theoretical formulations of the AMCBFM are exhibited in detail. In Section 3, the validity and efficiency of the AMCBFM are assessed and compared with the MoM and the CBFM. Section 4 presents a summary of this paper and further research plans on this topic.


## 2. Theoretical Formulations

### 2.1. EM Scattering from a Rough Surface Using the MoM

Assume that there is a tapered plane waveψinc incident upon a PEC rough surface, as shown in Figure 1, where z=f(x) is a Gaussian rough surface with a simulated length L. θi and θs are the incident angle and the scattered angle, respectively. Space Ω0 with the relative permittivity ε0 and the relative permeability μ0 is the incident space. Here, space Ω0 is assumed to be the free space; i.e., ε0=1, and μ0=1. ψsca represents the scattered wave in space Ω0. The position vector is r=xx^+zz^.Figure 1
Geometric model of the EM scattering from a PEC rough surface.For a Horizontal (H) polarized incident wave, the surface integral equation for this scattering problem is [17](2)ψincr=∫sds′Gr,r′n^′·∇ψtr′where S denotes the interface of the PEC rough surface, n^′ is the unit normal vector of the rough surface, Gr,r′=i/4H01k0r-r′ is the two-dimensional Green’s function in the space Ω0, and ψtr′=ψincr′+ψscar′ represents the total wave including the incident wave and the scattered wave in space Ω0.For the Vertical (V) case, the surface integral equation is as follows [17].(3)ψincr=12ψtr-∫Sds′ψtr′n^′·∇Gr,r′By using the MoM with the pulse basis functions and the point matching technique, the integral equations (2) and (3) can be discretized into the following matrix equations.(4)A-·U=b(5)B-·I=VThe elements of the impedance matrix A- and B- in detail are(6)Amn=Δx14H01k0r-r′,m≠niΔx41+i2πln⁡γk04eΔlm,m=n(7)Bmn=-Δx14H01k0r-r′,m≠n12-Δx4πf′′xm1+f′xm2,m=nwhere k0 is the wavenumbers in space Ω0 and Δlm=Δx1+(f′(xm))2. The subscripts m and n are the field point and source point along the rough surface, respectively. The elements of the right vector b and V are ψincrm, and the left vectors U and I are the unknown current vectors.


### 2.2. Theoretical Formulations of AMCBFM

Analyzing (4) and (5), the two equations can be rewritten as(8)Z-·J=Vwhere Z- is the impedance matrix with the dimensions N×N, the dimensions of the right vector V are N×1, the dimensions of the unknown current vector J are N×1, and N is the number of the discrete segments along the rough surface.In the AMCBFM, theN segments of the rough surface are first divided into M cells; every cell contains Ni(i=1,2,…,M) segments, and N1+N2+…+NM=N. Equation (8) can be rewritten as [18](9)Z-11Z-12⋯Z-1MZ-21Z-22⋯Z-2M⋮⋮⋯⋮Z-M1Z-M2⋯Z-MMJ1TotJ2Tot⋮JMTot=V1V2⋮VMwhere Z¯iji=1,2,…,M;j=1,2,…,M is the mutual impedance matrix of the ith cell, and the jth cell with the dimensions Ni×Nj. JiTot represents the total current of the ith cell, which represents the induced current along the ith cell. Vi indicates the incident wave upon the ith cell.(  1) The Primary Total Current J i T o t ( p ). For the ith(i=1,2,……,M) cell, only the self-interaction is considered to generate the PCBFs JiP such as the CBFM(10)Z-ii·Jip=Vii=1,2,…,Mwhere Z-ii is the self-interaction impedance matrix of the ith cell, and Jip and Vi are the PCBFs and the incident wave vector of the ith cell, respectively. Using the lower upper (LU) decomposition technique, one can obtain the PCBFs Jip. The lower and the upper triangular matrices of Z-ii are obtained and saved.Assuming the total currentJiTot is only from the PCBFs, it is called the primary total current JiTot(p), and it is calculated as(11)JiTotp=apipJipwhere apip is the coefficient of the PCBFs.Taking (11) into (9), one can get the following.(12)Z-11Z-12⋯Z-1MZ-21Z-22⋯Z-2M⋮⋮⋯⋮Z-M1Z-M2⋯Z-MMJ1TotpJ2Totp⋮JMTotp=V1V2⋮VMThen, by multiplying (12) by JipH, where JipH is the conjugate transpose of Jip, one can obtain the following.(13)J1pHZ-11J1pJ1pHZ-12J2p⋯J1pHZ-1MJMpJ2pHZ-21J1pJ2pHZ-22J2p⋯J2pHZ-2MJMp⋮⋮⋯⋮JMpHZ-M1J1pJMpHZ-M2J2p⋯JMpHZ-MMJMpap1pap2p⋮apMp=J1pHV1J2pHV2⋮JMpHVMUsing the LU decomposition technique, the PCBFs coefficients ap1(p)ap2(p)⋯apM(p)T and the primary total current JiTot(p) are acquired. It is obviously found that the mutual coupling effect of the subcells is neglected, and thus JiTot(p) is not accurate. However, the coefficients ap1(p)ap2(p)⋯apM(p)T can be used to get the SCBFs.(  2) The First Total Current J i T o t ( 1 ). Considering the mutual coupling effect of the subcells when building the SCBFs, for the ith(i=1,2,……,M) cell, the first order of the SCBFs JiS1 is defined as (14)Z-iiJiS1=-∑i=1,j=1,j≠iMapjpZ-ijJjpwhere JiS1is the first order SCBFs. In addition, the coefficients ap1(p)ap2(p)⋯apM(p)T are obtained using (13).The first total currentJiTot(1) is defined as(15)JiTot1=api1Jip+a1i1JiS1where api1 and a1i1 are the coefficients of the first order of SCBFs.Taking (15) into (9), one can get the following.(16)Z-11J1pZ-12J2p⋯Z-1MJMpZ-11J1S1Z-12J2S1⋯Z-1MJMS1Z-21J1pZ-22J2p⋯Z-2MJMpZ-21J1S1Z-22J2S1⋯Z-2MJMS1⋮⋮⋯⋮⋮⋮⋯⋮Z-M1J1pZ-M2J2p⋯Z-MMJMpZ-M2J1S1Z-M2J2S1⋯Z-MMJMS1ap11ap21⋮apM1a111a121⋮a1M1=V1V2⋮VMIt is readily seen that (16) has 2M unknowns and M equations. In order to solve this problem, (16) is multiplied by JipH and JiS1H, which results in the following.(17)J1pHZ-11J1p⋯J1pHZ-1MJMpJ1pHZ-11J1S1⋯J1pHZ-1MJMS1⋮⋯⋮⋮⋯⋮JMpHZ-M1J1p⋯JMpHZ-MMJMpJMpHZ-M1J1S1⋯JMpHZ-MMJMS1J1S1HZ-11J1p⋯J1S1HZ-1MJMpJ1S1HZ-11J1S1⋯J1S1HZ-1MJMS1⋮⋯⋮⋮⋯⋮JMS1HZ-M1J1p⋯JMS1HZ-MMJMpJMS1HZ-M1J1S1⋯JMS1HZ-MMJMS1ap11⋮apM1a111⋮a1M1=J1pHV⋯JMS1HVMJ1S1HV⋯JMS1HVMTThe CBFs coefficients ap1(1)⋯apM(1)a11(1)⋯a1M(1)T and the secondary total current JiTot(1) are obtained.(3) The n t h Total Current J i T o t ( n ) . Similar to the first order of SCBFs, the nth order of the SCBFs is defined as(18)Z-iiJiSn=-∑i=1,j=1,i≠jMan-1jn-1Z-ijJjSn-1where JjSn-1 is the (n-1)th order of the SCBFs, and a(n-1)j(n-1) is the coefficient of JjSn-1.Thenth total current JiTot(n) is defined as follows.(19)JiTotn=apinJip+a1inJiS1+⋯+aninJiSn(  4) The Current Criterion.The current criterion that determines whether the higher-order SCBFs are needed is defined as [16](20)ε=JTotn+12-JTotn2JTotn+12×100%where ·2 is the L2 norm. It is readily seen that the order of the SCBFs can be adaptively halted.(  5) The Percentage Error.The percentage error of the AMCBFM or the CBFM is defined as(21)PE=JAMCBFM/CBFM-JMoM2JMoM2×100%where JAMCBFM/CBFM represents the induced current that is obtained using the AMCBFM or the CBFM, JMoM denotes the induced current that is calculated using the MoM, and ·2 is the L2 norm.


## 2.1. EM Scattering from a Rough Surface Using the MoM

Assume that there is a tapered plane waveψinc incident upon a PEC rough surface, as shown in Figure 1, where z=f(x) is a Gaussian rough surface with a simulated length L. θi and θs are the incident angle and the scattered angle, respectively. Space Ω0 with the relative permittivity ε0 and the relative permeability μ0 is the incident space. Here, space Ω0 is assumed to be the free space; i.e., ε0=1, and μ0=1. ψsca represents the scattered wave in space Ω0. The position vector is r=xx^+zz^.Figure 1
Geometric model of the EM scattering from a PEC rough surface.For a Horizontal (H) polarized incident wave, the surface integral equation for this scattering problem is [17](2)ψincr=∫sds′Gr,r′n^′·∇ψtr′where S denotes the interface of the PEC rough surface, n^′ is the unit normal vector of the rough surface, Gr,r′=i/4H01k0r-r′ is the two-dimensional Green’s function in the space Ω0, and ψtr′=ψincr′+ψscar′ represents the total wave including the incident wave and the scattered wave in space Ω0.For the Vertical (V) case, the surface integral equation is as follows [17].(3)ψincr=12ψtr-∫Sds′ψtr′n^′·∇Gr,r′By using the MoM with the pulse basis functions and the point matching technique, the integral equations (2) and (3) can be discretized into the following matrix equations.(4)A-·U=b(5)B-·I=VThe elements of the impedance matrix A- and B- in detail are(6)Amn=Δx14H01k0r-r′,m≠niΔx41+i2πln⁡γk04eΔlm,m=n(7)Bmn=-Δx14H01k0r-r′,m≠n12-Δx4πf′′xm1+f′xm2,m=nwhere k0 is the wavenumbers in space Ω0 and Δlm=Δx1+(f′(xm))2. The subscripts m and n are the field point and source point along the rough surface, respectively. The elements of the right vector b and V are ψincrm, and the left vectors U and I are the unknown current vectors.


## 2.2. Theoretical Formulations of AMCBFM

Analyzing (4) and (5), the two equations can be rewritten as(8)Z-·J=Vwhere Z- is the impedance matrix with the dimensions N×N, the dimensions of the right vector V are N×1, the dimensions of the unknown current vector J are N×1, and N is the number of the discrete segments along the rough surface.In the AMCBFM, theN segments of the rough surface are first divided into M cells; every cell contains Ni(i=1,2,…,M) segments, and N1+N2+…+NM=N. Equation (8) can be rewritten as [18](9)Z-11Z-12⋯Z-1MZ-21Z-22⋯Z-2M⋮⋮⋯⋮Z-M1Z-M2⋯Z-MMJ1TotJ2Tot⋮JMTot=V1V2⋮VMwhere Z¯iji=1,2,…,M;j=1,2,…,M is the mutual impedance matrix of the ith cell, and the jth cell with the dimensions Ni×Nj. JiTot represents the total current of the ith cell, which represents the induced current along the ith cell. Vi indicates the incident wave upon the ith cell.(  1) The Primary Total Current J i T o t ( p ). For the ith(i=1,2,……,M) cell, only the self-interaction is considered to generate the PCBFs JiP such as the CBFM(10)Z-ii·Jip=Vii=1,2,…,Mwhere Z-ii is the self-interaction impedance matrix of the ith cell, and Jip and Vi are the PCBFs and the incident wave vector of the ith cell, respectively. Using the lower upper (LU) decomposition technique, one can obtain the PCBFs Jip. The lower and the upper triangular matrices of Z-ii are obtained and saved.Assuming the total currentJiTot is only from the PCBFs, it is called the primary total current JiTot(p), and it is calculated as(11)JiTotp=apipJipwhere apip is the coefficient of the PCBFs.Taking (11) into (9), one can get the following.(12)Z-11Z-12⋯Z-1MZ-21Z-22⋯Z-2M⋮⋮⋯⋮Z-M1Z-M2⋯Z-MMJ1TotpJ2Totp⋮JMTotp=V1V2⋮VMThen, by multiplying (12) by JipH, where JipH is the conjugate transpose of Jip, one can obtain the following.(13)J1pHZ-11J1pJ1pHZ-12J2p⋯J1pHZ-1MJMpJ2pHZ-21J1pJ2pHZ-22J2p⋯J2pHZ-2MJMp⋮⋮⋯⋮JMpHZ-M1J1pJMpHZ-M2J2p⋯JMpHZ-MMJMpap1pap2p⋮apMp=J1pHV1J2pHV2⋮JMpHVMUsing the LU decomposition technique, the PCBFs coefficients ap1(p)ap2(p)⋯apM(p)T and the primary total current JiTot(p) are acquired. It is obviously found that the mutual coupling effect of the subcells is neglected, and thus JiTot(p) is not accurate. However, the coefficients ap1(p)ap2(p)⋯apM(p)T can be used to get the SCBFs.(  2) The First Total Current J i T o t ( 1 ). Considering the mutual coupling effect of the subcells when building the SCBFs, for the ith(i=1,2,……,M) cell, the first order of the SCBFs JiS1 is defined as (14)Z-iiJiS1=-∑i=1,j=1,j≠iMapjpZ-ijJjpwhere JiS1is the first order SCBFs. In addition, the coefficients ap1(p)ap2(p)⋯apM(p)T are obtained using (13).The first total currentJiTot(1) is defined as(15)JiTot1=api1Jip+a1i1JiS1where api1 and a1i1 are the coefficients of the first order of SCBFs.Taking (15) into (9), one can get the following.(16)Z-11J1pZ-12J2p⋯Z-1MJMpZ-11J1S1Z-12J2S1⋯Z-1MJMS1Z-21J1pZ-22J2p⋯Z-2MJMpZ-21J1S1Z-22J2S1⋯Z-2MJMS1⋮⋮⋯⋮⋮⋮⋯⋮Z-M1J1pZ-M2J2p⋯Z-MMJMpZ-M2J1S1Z-M2J2S1⋯Z-MMJMS1ap11ap21⋮apM1a111a121⋮a1M1=V1V2⋮VMIt is readily seen that (16) has 2M unknowns and M equations. In order to solve this problem, (16) is multiplied by JipH and JiS1H, which results in the following.(17)J1pHZ-11J1p⋯J1pHZ-1MJMpJ1pHZ-11J1S1⋯J1pHZ-1MJMS1⋮⋯⋮⋮⋯⋮JMpHZ-M1J1p⋯JMpHZ-MMJMpJMpHZ-M1J1S1⋯JMpHZ-MMJMS1J1S1HZ-11J1p⋯J1S1HZ-1MJMpJ1S1HZ-11J1S1⋯J1S1HZ-1MJMS1⋮⋯⋮⋮⋯⋮JMS1HZ-M1J1p⋯JMS1HZ-MMJMpJMS1HZ-M1J1S1⋯JMS1HZ-MMJMS1ap11⋮apM1a111⋮a1M1=J1pHV⋯JMS1HVMJ1S1HV⋯JMS1HVMTThe CBFs coefficients ap1(1)⋯apM(1)a11(1)⋯a1M(1)T and the secondary total current JiTot(1) are obtained.(3) The n t h Total Current J i T o t ( n ) . Similar to the first order of SCBFs, the nth order of the SCBFs is defined as(18)Z-iiJiSn=-∑i=1,j=1,i≠jMan-1jn-1Z-ijJjSn-1where JjSn-1 is the (n-1)th order of the SCBFs, and a(n-1)j(n-1) is the coefficient of JjSn-1.Thenth total current JiTot(n) is defined as follows.(19)JiTotn=apinJip+a1inJiS1+⋯+aninJiSn(  4) The Current Criterion.The current criterion that determines whether the higher-order SCBFs are needed is defined as [16](20)ε=JTotn+12-JTotn2JTotn+12×100%where ·2 is the L2 norm. It is readily seen that the order of the SCBFs can be adaptively halted.(  5) The Percentage Error.The percentage error of the AMCBFM or the CBFM is defined as(21)PE=JAMCBFM/CBFM-JMoM2JMoM2×100%where JAMCBFM/CBFM represents the induced current that is obtained using the AMCBFM or the CBFM, JMoM denotes the induced current that is calculated using the MoM, and ·2 is the L2 norm.


## 3. Numerical Simulations and Discussions

The validity and efficiency of the AMCBFM are first assessed by comparing its results with the results that were obtained from the CBFM and the MoM. Both the H and V polarized incident waves are considered. In Figure2, the bistatic scattering coefficients that are obtained using these three methods are shown. The induced currents along the rough surface that are obtained by the AMCBFM and the MoM are given in Figure 3. The frequency of the incident wave is f=1.5GHz, and the incident angle is θi=30∘. The rms height and the correlation length of the simulated rough surface are δ=0.1λ and l=1.0λ, respectively. The length of the simulated Gaussian rough surface is L=409.6λ with a discrete interval of Δx=0.1λ. In other words, there are N=4096 segments along the rough surface. Thus, the dimension of the MoM matrix is 4096×4096. The rough surface is equally divided into 4 cells. Thus, the dimensions of the CBFM matrix and the AMCBFM matrix are only 1024×1024. It is clearly determined that the dimensions of the dealing matrix have been significantly reduced in the AMCBFM and the CBFM. The orders of the SCBFs of the CBFM are set as n=2. The current criterion of the AMCBFM is ε=0.1 and the halted order of the SCBFs is n=1. From Figures 2 and 3, it is readily seen that the AMCBFM and the CBFM agree well with the MoM in both the H and V polarization. However, the halted order of SCBFs of AMCBFM is n=1 through the criterion ε=0.1. The orders of the SCBFs of the CBFM are set as n=2.Figure 2
The bistatic scattering coefficients obtained by the different methods. (a) H polarization and (b) V polarization.
(a) (b)Figure 3
The currents obtained by the AMCBFM and the CBFM. (a) H polarization and (b) V polarization.
(a) (b)Moreover, in Table1, the CPU times of the AMCBFM, the CBFM, and the MoM are shown. It is obviously found that the AMCBFM costs less than half the CPU time of the MoM. In addition, Table 1 compares the percentage errors of the AMCBFM and the CBFM. It is readily seen that the AMCBFM with the order n=1 is more efficient than the CBFM with the order n=2. The results are implemented on a computer with an Intel Core i7-6700 processor, 4 GB of memory, and the Microsoft Windows 10 operating system.Table 1
The CPU time and the percentage error of the AMCBFM and the CBFM.
polarization method CPU Time (seconds) Percentage error H AMCBFM 6.132 5.198 × 1 0 - 3 CBFM 6.689 8.611 × 1 0 - 3 MoM 14.016 ——  V AMCBFM 6.243 1.987 × 1 0 - 3 CBFM 6.647 1.305 × 1 0 - 3 MoM 14.138 ——One may note that only one sample of the rough surface is used in Figures2 and 3. In order to further test the efficiency of the proposed algorithm, the results obtained with δ=0.5λ and l=1.5λ are shown in Table 2. These results are based on the rough surface with a simulated length of L=409.6λ, a discrete interval of Δx=0.1λ, and an incident angle of θi=30∘. The CPU time of the MoM is 12.995 seconds for H polarization and 13.01 seconds in the V case. In the AMCBFM, the rough surface is equally divided into 4 cells. Thus, the dimensions of the AMCBFM matrix are only 1024×1024. The current criteria of AMCBFM are set as ε=1×10-1, ε=1×10-2, ε=1×10-3, ε=1×10-4, and ε=1×10-5. From Table 2, it is readily seen that the order of the SCBFs can be adaptively halted in the AMCBFM. Furthermore, the current criterion ε=1×10-1 is sufficient for retaining the accuracy of the induced current. By comparing the CPU times between the AMCBFM and the MoM, it is easily determined that the AMCBFM with ε=1×10-1 costs less than half the CPU time of the MoM.Table 2
The results of the AMCBFM with different error thresholds.
polarization Error threshold Orders of SCBFs Percentage error CPU Time (seconds) H 1 × 1 0 - 1 1st 1.206 × 1 0 - 2 6.022 1 × 1 0 - 2 1st 1.206 × 1 0 - 2 6.068 1 × 1 0 - 3 1st 1.206 × 1 0 - 2 6.131 1 × 1 0 - 4 2nd 3.168 × 1 0 - 3 9.625 1 × 1 0 - 5 6th 6.008 × 1 0 - 4 43.368  V 1 × 1 0 - 1 1st 1.145 × 1 0 - 2 6.164 1 × 1 0 - 2 1st 1.145 × 1 0 - 2 6.287 1 × 1 0 - 3 1st 1.145 × 1 0 - 2 6.302 1 × 1 0 - 4 2nd 1.386 × 1 0 - 3 9.188 1 × 1 0 - 5 2nd 1.386 × 1 0 - 3 9.172


## 4. Conclusions

In this paper, the AMCBFM was used to simulate the EM scattering from a Gaussian rough surface. Good agreement between the AMCBFM and the MoM was found. In addition, the AMCBFM costs less than half the CPU time of the MoM. Furthermore, compared with the traditional CBFM, the lower order of the CBFs of the AMCBFM is sufficient for guaranteeing the accuracy. Note that the scattering model is only a 2D model. In the future, this method will be expanded to solve the more realistic 3D scattering problem. In addition, this method can also be applied to calculate the dielectric rough surface with a cylinder located above it. Finally, the AMCBFM can be combined with the compressing sensing technology in order to solve the EM scattering over a frequency band.



---

*Source: 7985421-2019-04-11.xml*