Abstract
A broadband radar cross section (RCS) calculation approach is proposed based on the characteristic basis function method (CBFM). In the proposed approach, the desired arbitrary frequency band is adaptively divided into multiple subband in consideration of the characteristic basis functions (CBFs) number, which can reduce the universal characteristic basis functions (UCBFs) numbers after singular value decomposition (SVD) procedure at lower subfrequency band. Then, the desired RCS data can be obtained by splicing the RCS data in each subfrequency band. Numerical results demonstrate that the proposed method achieve a high accuracy and efficiency over a wide frequency range.
1. Introduction
With the development of detection technology, the topic of efficient and accurate solution of the electromagnetic scattering problems over a broad frequency band has drawn increasing attention. The broadband electromagnetic scattering is very important in many fields such as modern radar target recognizing, microwave imaging, and microwave remote sensing. We can obtain more comprehensive information from broadband RCS data, such as volume information, contour information, details information, and so on, which has important significance in stealth and antistealth fields.
The traditional method of moments (MoM) can be applied for accurate solution of broadband electromagnetic problems. To obtain the RCS data over a wide band range using MoM, one has to repeat the calculation processes at each frequency point over the desired band, and this is inefficient when the object size as well as the number of frequency sampling points is large. Recently, several techniques have been proposed to circumvent this problem. In [1], the model-based parameter estimation based on rational function approximation is used to reduce the number of frequency points at which solutions or samples are required in broadband RCS calculation. In the papers [2, 3], asymptotic waveform estimation (AWE) technique has been successfully applied to obtain the broadband RCS data. But AWE is accurate only around the frequency of expansion; its accuracy deteriorates beyond a certain bandwidth. In [4], the impedance matrix is computed at relatively large frequency intervals and then interpolated to approximate its values. However, all the above techniques need resort to iterative methods and thus the target’s monostatic RCS cannot be calculated very efficiently. In the papers [5–7], an algorithm that utilizes the characteristic basis function method (CBFM) [8, 9] has been presented for computing the broadband RCS data. This approach takes the characteristic basis functions (CBFs) at the highest frequency after singular value decomposition (SVD) procedure as universal characteristic basis functions (UCBFs). However, it should be noted that the number of UCBFs is higher than necessary and the computational complexity will be increased when applying the UCBFs to the lower frequency points. In many practical application of engineering, one often need to analyze the electromagnetic scattering characteristics over arbitrary frequency band. In this paper, we introduce a new adaptive construction method that can calculate the broadband RCS over any given frequency band. This method leads to relatively small number of UCBFs at lower subfrequency band, so the fill process of the impedance matrix has been sped up and the memory requirements have been reduced significantly.
The rest of the paper is organized as follows: the detail implementation process of the CBFM is presented in Section 2. Section 3 derives the adaptive construction method of the UCBFs. In Section 4, numerical results for two test targets are presented and compared with the UCBF method over the bandwidth. Conclusions are drawn in Section 5.
2. Universal Characteristic Basis Functions Method
The algorithm begins by modeling the target at the highest frequency point of desired frequency band, using conventional MoM formulation, the electric field integral equation (EFIE) can be transformed into a dense, complex linear system of the form (1) where represents the generalized impedance matrix of dimension , and are vectors, and is the number of unknowns.
In this paper, we apply the CBFM based on SVD to generate the characteristic basis functions (CBFs) by a series of illuminating wave. Let us suppose that we calculate the CBFs at the highest frequency point , the CBFM-SVD begins with dividing the geometry of the object into blocks, characterizing these parts alone by using CBFs, which represent not only the self-interactions but also account for the mutual coupling effects between the self-blocks and the remaining blocks. We consider totally plane waves, for examples , here, the number of plane wave angles is deliberately overestimated. The CBFs on the subblock are constructed by the equation where, represents a matrix containing the CBFs stored column-wise; represents a matrix containing the excitation vector used to illuminate the subblock ; and is a impedance matrix of the subblock , is the number of unknowns in the subblock . Typically, the number of plane waves we have used to generate the CBFs would exceed the number of degrees of freedom (DoFs) associated with the block, to eliminate the redundant information in caused by overestimation, an orthogonalization procedure based on SVD method is used to reduce the final number of CBFs, that also serves to improve the condition number of the reduced matrix.
We arranged the CBFs obtained from multiple incident conditions in a rectangular matrix and retained only those whose singular values are above a certain threshold (typically chosen to be 1.0e-3, associated with the first and the most influential singular values). For the sake of simplicity, we assume that all of the blocks contain the same number of CBFs after SVD for each of the blocks. The solution to the entire problem is then expressed as a linear combination of the CBFs as follows: In (3), represents the th CBFs of the block after SVD. represents the unknown weight coefficients and the solution is obtained at frequency . Substitute formula (3) into MoM equation (1), we can generate the reduced matrix where, represents the reduced impedance matrix of dimension , each element of which can be expressed as the elements of can be expressed as Finally, the reduced system (4) is solved for the unknown weight coefficients and the CBFs are obtained at frequency . If these CBFs can be used in the whole frequency band, then we called them the UCBFs. The UCBFs capture the electromagnetic behaviour of the lower frequencies as well, enable one to solve the scattering for any frequency sample in the band without going through the time-consuming process to generate the CBFs anew.
3. Adaptive Generation of UCBFs
In the proposed approach, the any given frequency band is adaptively divided into multiple subbands in consideration of the UCBFs number, then we can calculate the RCS utilizing the UCBFs generated at the highest frequency point in each subband, and the whole RCS data can be obtained by splicing the RCS data of each frequency subband. For a given arbitrary frequency band, the adaptive method named after adaptive UCBF method (Adaptive UCBFM) is described as follows:(1)setting the initial frequency point and ;(2)generating the UCBFs in frequency point and , the total UCBFs number of all subblocks is In (7), represent the UCBFs number of subblock at frequency point , and represent the ones at frequency point . Given a threshold level , typically chosen to be 0.5), if , the RCS over the whole frequency band can be calculated using the UCBFs generating at repeatedly, else take , repeat step 2) in and until ;(3)getting the wide RCS data over a given arbitrary frequency band by splicing the RCS data in each frequency subband.
4. Numerical Results
In order to validate the accuracy and efficiency of the proposed approach; we present the results for the problem of scattering by a PEC plate and a NASA almond. A conventional triangular patch segmentation and RWG basis functions [10] is employed. The results have been compared with those derived by using UCBF and MoM (FEKO). All the results are computed on an Intel(R) Core (TM)2 Duo PC with 3.40 GHz processor and 4 GB RAM, only one core is used. The percentage error of Adaptive UCBFM is defined as
As a validation example, let us start by considering the problem of scattering by a PEC plate over a frequency range of 0.3 GHz to 3 GHz, which is 25 cm × 25 cm in size located on the plane. The geometry is divided into 1176 triangular patches with an average length of at 3 GHz, and this gives rise to 1764 unknowns. The geometry is then subdivided into 4 blocks, each block has been extended by in all directions and analyzed for a spectrum of plane waves incident from and , with , . The UCBFs after SVD retained at each subband were shown in Table 1.
The object is then excited by a normally incident (°, °, polarization angle = 90°) plane wave. Then the monostatic RCS results (21 frequency sampling points) using the method proposed in this paper over the bandwidth are shown in Figure 1. The bistatic RCS calculated for HH polarization at 1.5 GHz is shown in Figure 2. The results obtained by using the commercial MoM (FEKO) and the UCBF approaches are also included in Figure 2 for the sake of comparison.
From Figures 1 and 2, we can see that the Adaptive UCBFM has a better match with classic MoM (FEKO) than UCBF without a compromise of the accuracy.
To further illustrate the accuracy of the method, we also present the results for the problem of scattering by a PEC almond shown in Figure 3, which is described in [11].
The total length of almond is 25 cm along axis . The discretisation in triangular patches is carried out at 3.1 GHz with a mean edge length of , this leads to a number of unknowns of 4110. The structure is divided into 4 blocks, as shown in Figure 3. Each block has been extended by in all directions. The number of UCBFs after SVD at each frequency subband was shown in Table 2.
The object is excited by a normally incident (°, °, polarization angle = 90°) plane wave. In Figure 4, the monostatic RCS values (19 frequency sampling points) have been obtained using the Adaptive UCBFM over a frequency range of 0.3 GHz to 3.1 GHz.
The 2.0 GHz bistatic RCS results are shown in Figure 5, the solution using FEKO and UCBF has been presented also for validation. The results of NASA PEC almond thus demonstrate that the proposed approach is able to accurately compute the RCS of target with sharp corners structure.
It should be noted that the number of UCBFs at lower frequency band derived by using the Adaptive UCBFM is evidently smaller than that obtained using the UCBF method, and thus serves to accelerate the fill procedure of reduced matrix at lower frequencies. Table 3 lists the entire CPU time (which includes the time to fill the reduced impedance matrix at each frequency point, calculates the UCBFs at the highest frequency point, and calculates the current coefficient on the surface of the target) and RAM requirements for broadband RCS calculation process. It can be seen that the Adaptive UCBFM has compressed much more CPU time and random memory demand than conventional MoM and UCBF.
5. Conclusion
This paper has presented a method for analyzing broadband electromagnetic behaviors of metallic objects. The proposed technique is based on the adaptive construction of the UCBFs over a wide frequency band. This method leads to relatively small number of UCBFs at lower frequencies, so the reduced matrix size can be decreased at lower frequency band, which leads to a significant reduction of solver time with respect to conventional CBFM procedure. The results thus demonstrate that the proposed approach is able to accurately compute the RCS in arbitrary frequency band, resulting in a substantial time-saving, without a compromise of the accuracy.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China under Grant no. 61172020 and the Specialized Research Fund for the Doctoral Program of Higher Education of China under Grant no. 20123401110006.