# A Truncated Singular Value Decomposition Enhanced Nested Complex Source Beam Method

Authors: Lianning Song; Zaiping Nie
Journal: International Journal of Antennas and Propagation
 (2017)

---

## Abstract

This work presents a novel matrix compression algorithm to improve the computational efficiency of the nested complex source beam (NCSB) method. The algorithm is based on the application of the truncated singular value decomposition (TSVD) to the multilevel aggregation, translation, and disaggregation operations in NCSB. In our implementation, the aggregation/disaggregation matrices are solved by the truncated far-field matching, which is based on the directional far-field radiation property of the complex source beams (CSBs). Furthermore, the translation matrices are obtained according to the beam width of CSBs. Due to the high directivity of the radiation patterns of CSBs, all the far-field related interaction matrices are low-ranked. Therefore, TSVD can be employed and a new set of equivalent sources can be constructed by a linear combination of the original CSBs. It is proved that the radiation power of the new sources is proportional to the square of the corresponding singular values. This provides a theoretical guideline to drop the insignificant singular vectors in the calculation. In doing so, the efficiency of the original NCSB method can be much improved while a reasonably good accuracy is maintained. Several numerical tests are conducted to validate the proposed method.


---

## Body

## 1. Introduction

Integral equation (IE) method has been intensively studied in the analysis of electromagnetic (EM) radiation and scattering in recent years. For a numerical solution, the IE is usually discretized into a matrix equation by the method of moments (MoM) [1]. In order to conduct large-scale simulations, a variety of numerical techniques have been developed based on MoM. These techniques are carried out efficiently either by reduction of the far-field interaction exploiting the physical or mathematical properties of the IE or by using well-designed basis/test functions to reduce the total number of unknowns.In the first category, the multilevel fast multipole algorithm (MLFMA) [2] is widely used, where the far-field interactions can be accelerated by aggregation, translation, and disaggregation operations. Besides it, the low rank property of the MoM far-field submatrices leads to a series of matrix decomposition methods, such as the multilevel matrix decomposition algorithm (MLMDA) [3] and the adaptive cross approximation (ACA) [4, 5]. For the second category, the higher order basis functions (HOBF) [6], phase extracted basis functions (PE) [7, 8], characteristic basis function (CBF) [9, 10], or body of revolution (BoR) [11] MoM are proposed to reduce the number of unknowns in a given problem.Following the first category, a complex source beam-method of moments (CSB-MoM) is recently proposed to accelerate the far-field interactions of MoM [12]. In this method, the object is divided into groups and complex source beams (CSBs) [13] are used to expand the fields of the basis functions residing in each group [14, 15]. Hence, the far-field interactions of these basis/test functions can be accounted for by their equivalent CSB expansions. Due to the directional nature of CSBs [16], the interactions usually involve only a small portion of the total CSBs. A multilevel version of this method is developed in [17] to further improve the computational efficiency. The branch cut issue of CSBs, which might degrade the accuracy in the evaluation of group interactions, can be avoided by a proper choice of the CSB parameters.Unfortunately, in this method, the CSB expansions at any level are calculated directly from the basis functions contained in the finest level. This operation is very computationally expensive. Hence, the application of this method to electrically large problems is prohibited. To overcome this difficulty, a nested complex source beam (NCSB) method is proposed by utilizing an equivalent relationship between adjacent levels [18]. This relationship is built by treating CSBs in the child group as new sources and applying the far-field matching to get their CSB expansions in the parent group. In doing this, the computational complexity of NCSB can be reduced to O(Nlog⁡N), where N is the number of unknowns. However, it should be noted that the CSBs involved in the translation process have to be determined empirically. Moreover, the aggregation and disaggregation matrices are still in a dense format.To fully exploit the directional property of CSBs, a truncated singular value decomposition (TSVD) method is applied to compress the aggregation, translation, and disaggregation matrices of NCSB in this paper. After SVD, a set of equivalent sources can be obtained by a linear combination of the original CSBs. A theoretical proof reveals that the radiation power of the new sources is proportional to the square of the corresponding singular values. This provides a guideline for the truncation in the calculation. Thereby, the proposed method not only leads to a significant improvement of the computational efficiency but also provides a flexible compromise between accuracy and computational cost.


## 2. Formulations

Given a 3D perfectly electrical conducting (PEC) body defined by its surface, a MoM matrix equation can be obtained as(1)ZI=ZNI+ZFI=V,where I is the unknown vector containing the expansion coefficients of the current, V is the excitation vector, and Z is the MoM impedance matrix. After grouping, this dense matrix can be separated into the near-field interaction part ZN and the far-field interaction counterpart ZF.In CSB-MoM, the far-field interactions between different groups are carried out by a series of CSBs launched on a complex equivalence surface enclosing each group [14]. The far-field part of the matrix-vector product (MVP) in CSB-MoM can be represented as(2)ZFI=Zmm′Im′=∑p,p′=θ,ϕYmpTL,mm′pp′Wm′p′Im′,where m and m′ denote the observation and source groups and L indicates the finest level (single level in CSB-MoM). Wm′p′ and Ymp are, respectively, the expansion matrix and local expansion matrix for both θ and ϕ components. The far-field matching technique can be used for constructing Wm′p′ and Ymp in electric field integral equation (EFIE), magnetic field integral equation (MFIE) [19], and PMCHWT integral equation [20]. TL,mm′pp′ is the translation matrix, of which the elements are expressed as(3)TL,mm′pp′qq′=p^q·G̿r~L,mq∣r~′L,m′q′·p^q′′,where p^q and p^q′′ denote the unit vectors of qth and q′th CSB, respectively. The complex position vectors r~L,mq and r~′L,m′q′ are the launch points of CSBs. G̿r~mq∣r~′m′q′ is the dyadic Green function with complex arguments.By using this representation, the CSB expansion coefficients for the source groupm′ can be expanded from the surface currents as(4)SL,m′=SL,m′θSL,m′ϕ=Wm′θWm′ϕIm′.On the other hand, the CSB expansion coefficients RL,m′ for the receiving group m are translated from the transmitting group m′:(5)RL,m′=TL,mm′pp′SL,m′.As has been mentioned previously, the directional property of CSBs can be used to reduce the computational cost in the translation procedure. The translation window can be set based on numerical experiments for different group sizes and CSB parameters, as is done in [18].


### 2.1. NCSB Formulations

Motivated by the idea of MLFMA [2], CSB-MoM can be naturally extended to a multilevel version if a similar octree data structure is adopted. However, the multilevel CSB or MLCSB proposed in [17] calculates the CSB expansion coefficients for each level directly from the basis functions in the finest level. The lack of a proper mechanism of interlevel aggregation and disaggregation prohibits this method from analyzing large-unknown problems. To end this, the NCSB is proposed in [18] to introduce the aggregation and disaggregation operations between every two adjacent levels through a proper far-field matching.In the following, we obtain the aggregation matrix in a slightly different way by further utilizing the directional property of the radiating far field of CSBs. Moreover, the symmetry is fully exploited so that only one aggregation matrix is needed for all the eight child groups of one parent group. First, a linear system is set up to build the equivalent relationship of CSBs between two adjacent levels:(6)ZlAlc=Fl+1,cc=1,2,…,8,where Zl is the matching matrix that connects the equivalent CSB sources with the far fields Fl+1,c. The expression of the matching matrix is the same as the matrix used for CSB expansion with electric current type equivalent sources in CSB-MoM [12, 14]. Different from CSB-MoM of single level, the equivalent CSB sources in Zl are launched in level l, and the far fields here are radiated by CSBs in child level l+1 with both θ and ϕ components:(7)Fl+1,c=El+1,cθθEl+1,cθϕEl+1,cϕθEl+1,cϕϕ,where the elements in each block matrix are(8)El+1,cp′pqlql+1=p^ql′G̿rql∣r~ql+1,c′p^ql+1.It is noted that, for a parent group in the octree, there are eight different child groups at most. Hence, eight corresponding aggregation matrices are required in the calculation. However, by fully exploiting the symmetry property, we discover that only one aggregation matrix is needed and others can be easily deduced from it. In the following, the subscriptc will be omitted for simplicity.In our method, the radiating far fields on the right-hand-sides of (6) can be fast calculated by truncating CSB directional fields within the paraxial regions. Once Zl and Fl+1,c are assembled, the aggregation matrix for level l(9)Al=A~lθθA~lθϕA~lϕθA~lϕϕ2Ql×2Ql+1can be numerically solved. Here Ql and Ql+1 are the numbers of CSBs in levels l and l+1. In this work, we apply the least square method to improve the stability of the solution. In (9), the element A~lp′p′′qlql+1 in the block aggregation matrix is the mapping coefficient between the ql+1th CSB in the child group and the qlth CSB in its parent group. By using this relationship, the disaggregation matrix can be easily obtained from the transpose of the aggregation matrix. The CSB expansion coefficients of the parent group in level l can be obtained efficiently from its child groups in level l+1 with the aggregation matrix:(10)Sl,m′=Sl,m′θSl,m′ϕ=∑n′∈childm′AlSl+1,n′.Similar to MLFMA, the CSB expansion coefficients of a receiving group in levell+1 is obtained from both of the translation in the same level and the disaggregation from its parent level l. Therefore, we have(11)Rl+1,m=Rl+1,mθRl+1,mϕ=AlTRl,nθRl,nϕ=AlTRl,n.As has been elaborated in the above, the NCSB algorithm can be constructed in much the same way as MLFMA. However, all the operation matrices involved in NCSB are dense matrices, in contrast with the diagonal matrices in MLFMA. Improvement can be made by utilizing the unique property of CSBs. Specifically, all the CSBs are directional; that is to say, a CSB can only directionally interact with another CSB in the far-field region. Hence, the dense aggregation/disaggregation and translation matrices are usually low-ranked. To take advantage of the directional property of CSBs, TSVD will be applied to compress these matrices in the following sections.


### 2.2. TSVD in Aggregation/Disaggregation Process

For the aggregation matrix in (9), the SVD factorization can be applied as(12)Al=UlAΣlAVlAH,where UlA2Ql×2Ql and VlA2Ql+1×2Ql+1 are unitary matrices and ΣlA=diag⁡{σ1,σ2,…,σ2Ql+1}2Ql×2Ql+1 is a diagonal matrix whose elements are the nonnegative real singular values listed in a descending order. The columns of VlA represent a new set of orthonormal child CSB modes in level l+1, which are obtained as linear combinations of the original CSBs. Matrix UlA defines a complete set of orthonormal parent CSB modes associated with the child modes. Diagonal matrix ΣlA maps the child modes to the corresponding parent modes with weighting of the singular values.The physical meaning of the above SVD can be interpreted as follows. One set of the child CSBs is coupled to another set of parent CSBs, with the coupling strength governed by the corresponding singular value. Hence, after SVD, the aggregation procedure is converted from beam-to-beam coupling into mode-to-mode coupling. Furthermore, we can prove that the radiated power by each child CSB mode is proportional to the square of the singular value (see Appendix). This power-related mode provides a helpful guideline for extracting the most significant part in the aggregation matrix.Ifrl is the effective rank of Al for a prescribed threshold ϵ, the SVD decomposition in (12) can be truncated and approximated as(13)Al≃U~lAΣ~lAV~lAH,where Σ~lA is formed by the first (largest) rl singular values in ΣlA satisfying σrl+1<ϵσ1. U~lA and V~lA denote the corresponding submatrices consisting of the first rl columns of UlA and VlA, respectively. Hence, V~lA contains all the child CSB modes having radiated power higher than ϵ2σ12. According to the prescribed accuracy, this reduced set of CSBs in child group suffices to represent the far field of the parent group. Since disaggregation is implemented as the transpose of aggregation, the TSVD procedure of disaggregation matrix can be implemented by fully utilizing U~lA and V~lA in (13).


### 2.3. TSVD in Translation Process

As the CSB decays rapidly in the directions orthogonal to the beam propagation axis, in the translation process, the interactions between far-field groups are carried out by CSBs whose beam axis are alined with the line connecting the group centers. These CSBs can be easily selected by setting a conical truncation window. To avoid amount of experiments needed to determine the minimal conical angle, in this paper, we first set the truncation angle when the field of a CSB decays to-20 dB of the aforementioned prescribed accuracy times the maximum value. According to this criterion, the CSB selection matrices Ll,mm′′ and Ll,mm′ are introduced to pick up the most significant CSBs in the transmitting and receiving groups; namely,(14)Tl,mm′pp′Ql×Ql=Ll,mm′Ql×KlT~l,mm′pp′Kl×KlLl,mm′′Kl×Ql,where only one element in each row (column) of Ll,mm′′ (Ll,mm′) is 1, if the corresponding CSB is within the conical truncation window; otherwise they are set to 0.After this truncation based on the beam width of CSBs, the dimension of the translation matrixT~l,mm′pp′ is reduced from Ql×Ql to Kl×Kl. However, it is still rank-deficient. The TSVD procedure is then applied to this reduced translation matrix to gain further compression:(15)T~l,mm′pp′≃U~l,mm′T,pp′Σ~l,mm′T,pp′V~l,mm′T,pp′H,where U~l,mm′T,pp′, Σ~l,mm′T,pp′, and V~l,mm′T,pp′ are defined similarly as in (13), except that the threshold is set as ϵ2 instead of ϵ to account for the effects from both of the transmitting and receiving groups.The proposed TSVD procedure is effective by providing new set of CSB modes based on the singular vectors. The significance in the total radiated power is governed by the descending singular values. The SVD process fully exploits the spatial window property of CSBs, which is independent of the geometry, current distribution, and types of IE. By defining an error threshold, redundant modes which contribute negligibly to the radiated power are excluded. Therefore, the aggregation, translation, and disaggregation matrices can be represented in a very compact form. Furthermore, the balance between accuracy and computational efficiency can be easily controlled by adjusting the truncation threshold.


## 2.1. NCSB Formulations

Motivated by the idea of MLFMA [2], CSB-MoM can be naturally extended to a multilevel version if a similar octree data structure is adopted. However, the multilevel CSB or MLCSB proposed in [17] calculates the CSB expansion coefficients for each level directly from the basis functions in the finest level. The lack of a proper mechanism of interlevel aggregation and disaggregation prohibits this method from analyzing large-unknown problems. To end this, the NCSB is proposed in [18] to introduce the aggregation and disaggregation operations between every two adjacent levels through a proper far-field matching.In the following, we obtain the aggregation matrix in a slightly different way by further utilizing the directional property of the radiating far field of CSBs. Moreover, the symmetry is fully exploited so that only one aggregation matrix is needed for all the eight child groups of one parent group. First, a linear system is set up to build the equivalent relationship of CSBs between two adjacent levels:(6)ZlAlc=Fl+1,cc=1,2,…,8,where Zl is the matching matrix that connects the equivalent CSB sources with the far fields Fl+1,c. The expression of the matching matrix is the same as the matrix used for CSB expansion with electric current type equivalent sources in CSB-MoM [12, 14]. Different from CSB-MoM of single level, the equivalent CSB sources in Zl are launched in level l, and the far fields here are radiated by CSBs in child level l+1 with both θ and ϕ components:(7)Fl+1,c=El+1,cθθEl+1,cθϕEl+1,cϕθEl+1,cϕϕ,where the elements in each block matrix are(8)El+1,cp′pqlql+1=p^ql′G̿rql∣r~ql+1,c′p^ql+1.It is noted that, for a parent group in the octree, there are eight different child groups at most. Hence, eight corresponding aggregation matrices are required in the calculation. However, by fully exploiting the symmetry property, we discover that only one aggregation matrix is needed and others can be easily deduced from it. In the following, the subscriptc will be omitted for simplicity.In our method, the radiating far fields on the right-hand-sides of (6) can be fast calculated by truncating CSB directional fields within the paraxial regions. Once Zl and Fl+1,c are assembled, the aggregation matrix for level l(9)Al=A~lθθA~lθϕA~lϕθA~lϕϕ2Ql×2Ql+1can be numerically solved. Here Ql and Ql+1 are the numbers of CSBs in levels l and l+1. In this work, we apply the least square method to improve the stability of the solution. In (9), the element A~lp′p′′qlql+1 in the block aggregation matrix is the mapping coefficient between the ql+1th CSB in the child group and the qlth CSB in its parent group. By using this relationship, the disaggregation matrix can be easily obtained from the transpose of the aggregation matrix. The CSB expansion coefficients of the parent group in level l can be obtained efficiently from its child groups in level l+1 with the aggregation matrix:(10)Sl,m′=Sl,m′θSl,m′ϕ=∑n′∈childm′AlSl+1,n′.Similar to MLFMA, the CSB expansion coefficients of a receiving group in levell+1 is obtained from both of the translation in the same level and the disaggregation from its parent level l. Therefore, we have(11)Rl+1,m=Rl+1,mθRl+1,mϕ=AlTRl,nθRl,nϕ=AlTRl,n.As has been elaborated in the above, the NCSB algorithm can be constructed in much the same way as MLFMA. However, all the operation matrices involved in NCSB are dense matrices, in contrast with the diagonal matrices in MLFMA. Improvement can be made by utilizing the unique property of CSBs. Specifically, all the CSBs are directional; that is to say, a CSB can only directionally interact with another CSB in the far-field region. Hence, the dense aggregation/disaggregation and translation matrices are usually low-ranked. To take advantage of the directional property of CSBs, TSVD will be applied to compress these matrices in the following sections.


## 2.2. TSVD in Aggregation/Disaggregation Process

For the aggregation matrix in (9), the SVD factorization can be applied as(12)Al=UlAΣlAVlAH,where UlA2Ql×2Ql and VlA2Ql+1×2Ql+1 are unitary matrices and ΣlA=diag⁡{σ1,σ2,…,σ2Ql+1}2Ql×2Ql+1 is a diagonal matrix whose elements are the nonnegative real singular values listed in a descending order. The columns of VlA represent a new set of orthonormal child CSB modes in level l+1, which are obtained as linear combinations of the original CSBs. Matrix UlA defines a complete set of orthonormal parent CSB modes associated with the child modes. Diagonal matrix ΣlA maps the child modes to the corresponding parent modes with weighting of the singular values.The physical meaning of the above SVD can be interpreted as follows. One set of the child CSBs is coupled to another set of parent CSBs, with the coupling strength governed by the corresponding singular value. Hence, after SVD, the aggregation procedure is converted from beam-to-beam coupling into mode-to-mode coupling. Furthermore, we can prove that the radiated power by each child CSB mode is proportional to the square of the singular value (see Appendix). This power-related mode provides a helpful guideline for extracting the most significant part in the aggregation matrix.Ifrl is the effective rank of Al for a prescribed threshold ϵ, the SVD decomposition in (12) can be truncated and approximated as(13)Al≃U~lAΣ~lAV~lAH,where Σ~lA is formed by the first (largest) rl singular values in ΣlA satisfying σrl+1<ϵσ1. U~lA and V~lA denote the corresponding submatrices consisting of the first rl columns of UlA and VlA, respectively. Hence, V~lA contains all the child CSB modes having radiated power higher than ϵ2σ12. According to the prescribed accuracy, this reduced set of CSBs in child group suffices to represent the far field of the parent group. Since disaggregation is implemented as the transpose of aggregation, the TSVD procedure of disaggregation matrix can be implemented by fully utilizing U~lA and V~lA in (13).


## 2.3. TSVD in Translation Process

As the CSB decays rapidly in the directions orthogonal to the beam propagation axis, in the translation process, the interactions between far-field groups are carried out by CSBs whose beam axis are alined with the line connecting the group centers. These CSBs can be easily selected by setting a conical truncation window. To avoid amount of experiments needed to determine the minimal conical angle, in this paper, we first set the truncation angle when the field of a CSB decays to-20 dB of the aforementioned prescribed accuracy times the maximum value. According to this criterion, the CSB selection matrices Ll,mm′′ and Ll,mm′ are introduced to pick up the most significant CSBs in the transmitting and receiving groups; namely,(14)Tl,mm′pp′Ql×Ql=Ll,mm′Ql×KlT~l,mm′pp′Kl×KlLl,mm′′Kl×Ql,where only one element in each row (column) of Ll,mm′′ (Ll,mm′) is 1, if the corresponding CSB is within the conical truncation window; otherwise they are set to 0.After this truncation based on the beam width of CSBs, the dimension of the translation matrixT~l,mm′pp′ is reduced from Ql×Ql to Kl×Kl. However, it is still rank-deficient. The TSVD procedure is then applied to this reduced translation matrix to gain further compression:(15)T~l,mm′pp′≃U~l,mm′T,pp′Σ~l,mm′T,pp′V~l,mm′T,pp′H,where U~l,mm′T,pp′, Σ~l,mm′T,pp′, and V~l,mm′T,pp′ are defined similarly as in (13), except that the threshold is set as ϵ2 instead of ϵ to account for the effects from both of the transmitting and receiving groups.The proposed TSVD procedure is effective by providing new set of CSB modes based on the singular vectors. The significance in the total radiated power is governed by the descending singular values. The SVD process fully exploits the spatial window property of CSBs, which is independent of the geometry, current distribution, and types of IE. By defining an error threshold, redundant modes which contribute negligibly to the radiated power are excluded. Therefore, the aggregation, translation, and disaggregation matrices can be represented in a very compact form. Furthermore, the balance between accuracy and computational efficiency can be easily controlled by adjusting the truncation threshold.


## 3. Numerical Results

In this section, several numerical tests are conducted to demonstrate the efficiency and validity of the proposed method.Firstly, the low rank property of the aggregation matrices of NCSB method is studied. Figure1 shows the normalized singular values of the aggregation matrices for level 2 (with electrical size of 8λ and aggregating from level 3) and level 3 (with electrical size of 4λ and aggregating from level 4). The numbers of CSBs for levels 2, 3, and 4 are 5738, 1758, and 682, as listed in Table 1. Considering both θ and ϕ components of CSBs, the dimensions of the aggregation matrices for level 2 and level 3 are 11476×3516 and 3516×1364, respectively. As shown in Figure 1, the singular values for both matrices decrease rapidly, indicating a quick decay of the radiated power from individual child CBS mode. Due to this fact, the aggregation matrices are truncated by thresholds of ϵ=10-2, 10-3, and 10-4 for comparison. Figure 2 depicts the relative errors of the radiated far field (an effective measure of the aggregation error) with respect to the observation angles. In particular, this relative error is defined as(16)eϕ=Eϕ-Eaggϕ2Eϕ2and here E(ϕ) is the reference field, Eagg(ϕ) is the field obtained via aggregations by TSVD, and ϕ is the observation angle which is set on a azimuth circle. In Figure 2, the reference field is radiated directly by the CSBs of a group with the size of 2λ in the fourth level. Then two aggregation steps with TSVD are performed to this group until the second level. Finally, the field Eagg(ϕ) is calculated by the CSBs in the second level. Figure 2 shows that the relative aggregation error can be effectively controlled by the TSVD threshold ϵ.Table 1
Computational statistics.
Method
Q
l (level 2–4)
Aggregation TSVDrank (level 2-3)
K
l (level 2–4)
Translation TSVDrank (level 2–4)
Aggregation &Translation Mem.
Setup time
Solving time
Total time
NCSB
5738/1758/682
—
1213/681/289
—
19230.5 MB
511 s
3864 s
4375 s
ϵ
=
1
0
-
2
635/509
59/48/66
1429.9 MB
936 s
1178 s
2114 s
ϵ
=
1
0
-
3
1154/730
119/131/108
3119.3 MB
958 s
1345 s
2303 s
ϵ
=
1
0
-
4
1600/919
206/199/142
5548.7 MB
992 s
1488 s
2480 sFigure 1
Normalized singular values of the aggregation matrix.Figure 2
Relative aggregation errors measured by the radiated far field.Next, we investigate the TSVD based translation process. As shown in Table1, after the aforementioned window selection, the number of CSBs for the reduced translation (Kl) drops to 1213, 681, and 289 for levels 2, 3, and 4, respectively. However, the information of the reduced translation matrices is still redundant, and the effective rank can be further compressed. Figure 3 shows the normalized singular values of the reduced translation matrices, where a rapid decrease can again be observed for individual CSB mode. As has been explained, the truncation thresholds for translation matrices should be set as ϵ2, namely, 10-4, 10-6, and 10-8, in this example. The relative translation error(17)eq=Rqtsvd-RqRqq=1,2,…,Klis depicted in Figure 4, where Rq and Rqtsvd are the CSB expansion coefficients in level 2, translated by the original translation and the TSVD counterpart, respectively. Even though the error is relatively higher for small indices of CSBs, the corresponding normalized expansion coefficients Rqnor=Rq/maxjRj are also relatively small. Hence, these CSBs will contribute to the total field insignificantly. In other words, TSVD maintains reasonably good translation accuracy for both the strongly coupled CSBs (with larger indices) and the weakly coupled CSBs (with smaller indices).Figure 3
Normalized singular values of the reduced translation matrix.Figure 4
Relative translation errors.Finally, to demonstrate the overall performance of the proposed method, the scattering problem of a diameter 32λ sphere is presented. It is illuminated by an incident plane wave (θi=0∘, φi=0∘) at 300MHz and solved by combined field integral equation (CFIE). Figure 5 shows the bistatic RCS results compared with the analytical solution of Mie series from 120∘ to 180∘. We can see a very good agreement is obtained for ϵ=10-3 and ϵ=10-4. Furthermore, a good agreement is still obtained for ϵ=10-2. To quantify the accuracy, the RMS value of the error E=δNCSB-TSVD-δMie/δMiemax has been calculated, δ being the RCS (m2). RMS errors of 0.18%, 0.0806%, and 0.0768% are found for the thresholds 10-2, 10-3, and 10-4. At the same time, the computational efficiency can be significantly improved by the proposed method. The statistics are summarized in Table 1. It can be observed that a great computational improvement is obtained in terms of CPU time and memory consumption. It is also shown that a good compromise between accuracy and the computational efficiency can be easily achieved, by a proper choice of the truncation parameter.Figure 5
Bistatic RCS of a PEC sphere with a diameter32λ, by NCSB-TSVD and Mie-series solution.


## 4. Conclusion

In this paper, we have proposed an efficient implementation of NCSB. Different from previous work, the aggregation process is first constructed by the truncated far-field matching. The dimension of the translation matrix is then reduced based on the beam width of the CSBs. Finally, TSVD is used to compress the aggregation, translation, and disaggregation matrices by fully exploiting the directional property of CSBs. It is shown that the radiated power of the new sources from SVD is proportional to the square of the corresponding singular values. This power-dependent mode provides a theoretical guideline for extracting the most significant CSB contributions in the calculation. Therefore, the desired balance between accuracy and computational efficiency can be easily controlled by adjusting the truncation threshold.



---

*Source: 9072421-2017-08-01.xml*