# High Precise Scattering Centers Models for Cone-Shaped Targets Based on Induced Currents

Authors: Qi-Feng Li; Kun-Yi Guo; Xin-Qing Sheng; Tian-Shu Liu
Journal: International Journal of Antennas and Propagation
 (2017)

---

## Abstract

Based on induced currents of cone-shaped targets, the high precise scattering center model is derived in this paper. The distribution characteristics of the induced currents and their relationships with the attributes of scattering centers are investigated in detail; the high precise scattering center model is obtained. In order to validate the scattering center model, numerical simulations of two cone-shaped targets are presented, and the accuracies of the models are validated through comparing the range profiles simulated by the models with those by a full-wave numerical method. The validation results demonstrate that this model is superior to the existing model in precisely characterizing the scattering centers induced by creeping waves.


---

## Body

## 1. Introduction

Scattered fields of radar targets at sufficiently high frequencies can be approximated as the sum of scattered fields from individual scattering sources, and these scattering sources are generally called scattering centers [1]. The term of scattering center has been expanded for now. Scattering phenomena under cases of lower frequency, creeping waves, for example, can also be described in terms of scattering centers and these scattering centers are also significant scattering components in some conditions [2, 3]. The characteristics of radar returns from a target are primarily determined by the attributes of its scattering centers, such as locations and scattering amplitudes, and their dependencies on aspect angle, frequency, and polarization [4]. Therefore, parametric scattering center models that concisely and precisely characterize radar returns have gained considerable attention in many radar applications, such as automatic target recognition (ATR) [5, 6] and radar image interpretation [7, 8].Parametric scattering center model is the combination of a series of analytical functions that describe the changes of location and scattering amplitude of scattering centers with frequency and aspect angle of line of sight (LOS) of radar. The existing scattering center models for monostatic radar configuration include the damped exponential model [9, 10], the geometrical theory of diffraction (GTD) based model [11, 12], the attributed scattering center model [13–15], and the sliding scattering center model [16]. These scattering center models are all developed for scattering problems of high frequency, such as specular reflection by large-scale smooth surface or diffraction by large-scale straight edge. The precision of these models is greatly challenged for describing the scattering problem of lower frequency, such as diffraction by small-scale edges and rescattered waves of creeping waves and travelling waves on smooth surface. And even worse, there are no generalized analytical solutions for these kinds of scattering problems.To tackle this problem, the models of the induced currents on the target are built at first; then the scattering center model is derived by the electric field integral expression based on the current model. The modelling of the induced currents on the target is easier than the scattering center modelling, for the precise numerical results of currents can be acquired by full-wave numerical computation. The well validated full-wave numerical method, namely, the hybrid finite element-boundary integral-multilevel fast multipole algorithm (FE-BI-MLFMA) [17–20] is used in this paper. It is found that for cone-shaped target, its scattering characteristics are mainly determined by one-dimensional (1D) currents on the lines in the plane constructed by the LOS and the rotational axis of the cone. Therefore, the modelling of the induced currents on the target is greatly simplified.Based on the built model of the 1D currents, the relationship between the characteristics of currents and the attributes of scattering centers is investigated; the formation processes of scattering centers by its contributing currents are analyzed. It is found that one parameter of current models is sometimes related to several scattering centers, and vice versa. In other words, the coupling of multiple scattering sources has been included in the current model. Moreover, the scattering phenomena of creeping waves are also included in the current model, and the corresponding scattering center can be easily acquired based on the current model. In order to validate the presented method, numerical simulations of two cone-shaped targets are investigated in detail, and the accuracies of the built models are validated through comparing the range profiles (RPs) simulated by the models with the real RPs of the two targets computed by FE-BI-MLFMA, as well as the RPs simulated by the attributed scattering center model obtained by traditional method.The remainder of the paper is organized as follows: Section2 presents the modelling of the induced currents on a blunt-nosed cone target. Section 3 presents the derived scattering center model based on the current model. Section 4 presents the validation of the built scattering center models. Finally, Section 5 provides the conclusions of this paper.


## 2. Parametric Expressions of Scattering Centers

The induced currents on a cone-shaped target are investigated in this section. The cone-shaped target is a perfect electric conductor (PEC). The geometry of the cone and the relative position of the target and radar are shown in Figure1. The body-fixed Cartesian coordinate system (x,y,z) is used to describe the geometry, its origin is at the bottom center of the target, and the rotationally symmetric axis of the cone is located along the z-axis, as shown in Figure 1. The geometrical parameters of the target are set as R=0.732m, r=0.061m, h=3.721m, and α=10.4°. The radar parameters are set as follows: the frequency range isf= 1~2 GHz, with sampling interval of Δf=12.5MHz.Figure 1
The geometry of the cone-shaped target and the radar.Suppose that the distance between the radar and the target satisfies the far-field condition (R0>2L2/λ, where R0 is the distance between the radar and the target and L is the size of the target). In far-field condition, the incident wave can be taken as plane wave; then the induced currents on the target illuminated by the plane waves from the monostatic radar are computed by FE-BI-MLFMA. The characteristics of currents on the flank and the bottom are investigated separately in the following.


### 2.1. Characteristics of Induced Currents on the Flank

As the cone is a rotationally symmetrical object, the change of the LOS in azimuth angle (denoted byϕ) does not affect the characteristics of induced currents. Therefore, without loss of generality, we set the change of LOS to be in the xoz plane, and the elevation angle of the LOS is denoted by θ. Then, the wave vector of incident wave can be expressed as k→i=2π/λ(-cos⁡θz^-sin⁡θx^).When the incident electric field is under vertical polarization, the magnetic field can be expressed asH→i=-Hmϕ^=-Hmy^. The induced currents on the cone are presented in Figure 2, where the used parameters of the incident waves are f=1GHz, θ=8°. Figure 2(a) shows the amplitude of the induced currents and Figure 2(b)L1 is the intersection line of the flank of the cone with the plane constructed by the axis of cone and the direction of the incident electric field, and L2 is the intersection line of the flank with the plane constructed by the axis of cone and the direction of the incident magnetic field. In our geometry setting, the generatrix in xoz plane of lit region is L1 and the generatrix in yoz is L2.Figure 2
Induced currents on the flank of the cone. (a) Amplitude; (b) phase.
(a)
(b)As shown in Figure2(a), the maximal values of currents are on L1 and the minimal values of currents are on L2. The conclusion agrees with the analytical results based on physical optics method (PO) [21–24]. The normals of the flank on L1 and L2 are n^L1=cos⁡αx^+sin⁡αz^ and n^L2=cos⁡αy^+sin⁡αz^. Based on PO, induced current can be approximately expressed as J=2n^×Hi, so the currents on L1 and L2 can be derived as JL1=2Hm(-cos⁡αz^+sin⁡αx^) and JL2=2Hm(sin⁡αx^). The incident electric field is generally set as 1 V/m in computation; then the incident magnetic field is Hm=1/Z, where Z is the wave impedance of free space.Optical currents onL1 computed by PO method are compared with the values computed by FE-BI-MLFMA (called full-wave currents for short), as presented in Figure 3 (θ=8°). Figure 3(a) shows the amplitude (A/m2) of the currents, and Figure 3(b) shows the unwrapped phase (rad) of the currents. The range of horizontal axis indicates the distance between a point on L1 and the top of the coneFigure 3
Amplitude and phase of the currents onL1. (a) Amplitude; (b) unwrapped phase.
(a)
(b)Figure3(a) shows that there are significant differences between the optical currents and the full-wave currents. The full-wave currents along L1 are more complex than that of the optical currents. Near the two ends of L1, the amplitude of the full-wave currents is obviously larger than that of the optical currents. Within the internal region of L1, the amplitude fluctuates as a combination form of periodicity and linearity. Figure 3(b) shows that the two results agree well with each other, and the phase of induced current linearly changes along L1.The generatrix inxoz plane in shadow region is denoted by L1′. The induced currents on L1′ are shown in Figure 4. The fluctuation feature and the values of the induced currents on L1′ are similar to the case of L1, which means the scattering components of creeping waves are the denominated scattering contributions when the elevation angle of the LOS is small (θ=8°<α). The range of horizontal axis indicates the distance between a point on L1′ and the top of the cone.Figure 4
Amplitude and phase of the induced currents onL1′ in shadow region. (a) Amplitude; (b) unwrapped phase.
(a)
(b)Based on the characteristics of the induced currents onL1 and L1′ under different LOSs, the full-wave currents are described by the following expressions:(1)Jρ,f,θ=Amρ,f,θexp⁡jϕmρ,f,θ(2)Amρ,f,θ=am1sin⁡am2θfρ+am3θρ+am4θf+am5θ+am6ρ+am7(3)ϕmρ,f,θ=am8θfρ+am9θρ,where (1) describes the complex value of induced current; its amplitude and phase are indicated by Am and ϕm. Equation (2) is proposed according to the amplitude characteristics of the induced currents on L1 and L1′ under different LOSs. Equation (3) is proposed according to the phase characteristics of the induced currents on L1 and L1′. ρ indicates the position of the current element on L1m=1 or L1′m=2. Am(·) and ϕm(·) are the amplitude and the phase of the currents, respectively. amn (n=1,2,…,9) indicates the parameters to be estimated. amn is irrelevant to frequency, but am2, am3, am4, am5, am8, and am9 are linear-circular correlated with θ.To show characteristics of the currents on the whole conical surface, we slice the cone along lineL1 and unfold the flank to a fan surface, as shown in Figure 5. The currents on the fan and the bottom are shown in Figure 6. The half-angle of the fan is γ=πR/l, where l is the outer radius of the fan.Figure 5
Unfolded fan.Figure 6
Currents on the fan and the bottom.Based on PO method, the two-dimensional currents can be expressed asJ(ρ,β)=2n^(β)×Hi=2Hmsin⁡αx^-cos⁡αcos⁡β/γπz^, where (ρ,β) is the coordinates of the individual current on the fan. The component of J(ρ,β) along vertical direction (v^=-cos⁡αz^+sin⁡αx^) can be expressed as (4). Therefore, the change of currents with β can be described by the currents on L1 and the dependant function Fβ.(4)Jvρ,β=Jvρ,0Fβ(5)Fβ=cos2⁡αcos⁡βγπ+sin2⁡α.Equations (4) and (5) describe the fact that the change of induced currents with β can be expressed by the currents on L1 and the dependent function F(·).


### 2.2. Characteristics of Induced Currents on the Bottom

Similarly, the induced currents on the bottom are investigated here. The amplitudes and phases of induced currents on the bottom are presented in Figures7(a) and 7(b), respectively. The parameters of the incident waves are the same as those of Figure 2. The line in xoz plane along +x axis is denoted by I1, the line in xoz plane along -x axis is denoted by I1′, and the line in yoz plane along +y axis is denoted by I2. The maximal values of currents are on I1 and the minimal values of currents are on I2. The scattering field from the bottom is mainly contributed by the currents on I1 and I1′.Figure 7
Currents on the bottom. (a) Amplitude; (b) phase.
(a)
(b)The currents onI1 are presented in Figure 8. The range of horizontal axis indicates the distance between a point on I1 and the point on the curved edge. From Figure 8(a), it can be seen that the amplitude of full-wave currents fluctuates as a combination form of periodicity and linearity. From Figure 8(b), it can be seen that the phase of induced current changes linearly along I1. The characteristics of the currents on I1′ are similar to I1.Figure 8
Currents onI1. (a) Amplitude; (b) unwrapped phase.
(a)
(b)These characteristics of the currents on the bottom are similar to the case of currents on the flank. Therefore, the currents on the bottom are also described by the expression of (1), only with different parameters. The currents on I1 and I1′ are indicated by the case of m=3. Similarly, the two-dimensional currents on the bottom are also expressed by (4).


### 2.3. Estimation of the Current Model Parameters

The coefficients in the presented formula,am1~am9, are estimated through matching the currents computed by the model with the currents computed by FE-BI-MLFMA. The used method for optimal estimation is genetic algorithm (GA) [25–28]. The currents on the flank computed by the model with estimated parameters (the simulated currents for short) are presented in Figures 9(a) and 10(a). For comparison, the currents computed by FE-BI-MLFMA (the real currents for short) are presented in Figures 9(b) and 10(b). The used parameters in computation aref = 1~2 GH and ∆f=0.125GHz, θ=8°. The range of horizontal axis indicates the distance between a point on L1 and the top of the cone.Figure 9
Amplitudes of the currents. (a) Simulated currents; (b) real currents.
(a)
(b)Figure 10
Unwrapped phases of the currents. (a) Simulated currents; (b) real currents.
(a)
(b)The errors between the simulated currents and the real currents are examined: the correlation coefficient of the amplitudes is92.81% and the root mean square error (RMSE) is 6.01×10-5V/m2; the correlation coefficient of the unwrapped phases is 98.72% and the RMSE is 2.99 rad.Other cases of incident angle (θ = 1°~10°) are also investigated. The parameters of currents are estimated with the same method. The errors between the simulated currents and the real currents are the correlation coefficient of the amplitudes is larger than 90% and the RMSE is less than 10-4V/m2; the correlation coefficient of the unwrapped phases is larger than 95% and the RMSE is less than 4 rad.


## 2.1. Characteristics of Induced Currents on the Flank

As the cone is a rotationally symmetrical object, the change of the LOS in azimuth angle (denoted byϕ) does not affect the characteristics of induced currents. Therefore, without loss of generality, we set the change of LOS to be in the xoz plane, and the elevation angle of the LOS is denoted by θ. Then, the wave vector of incident wave can be expressed as k→i=2π/λ(-cos⁡θz^-sin⁡θx^).When the incident electric field is under vertical polarization, the magnetic field can be expressed asH→i=-Hmϕ^=-Hmy^. The induced currents on the cone are presented in Figure 2, where the used parameters of the incident waves are f=1GHz, θ=8°. Figure 2(a) shows the amplitude of the induced currents and Figure 2(b)L1 is the intersection line of the flank of the cone with the plane constructed by the axis of cone and the direction of the incident electric field, and L2 is the intersection line of the flank with the plane constructed by the axis of cone and the direction of the incident magnetic field. In our geometry setting, the generatrix in xoz plane of lit region is L1 and the generatrix in yoz is L2.Figure 2
Induced currents on the flank of the cone. (a) Amplitude; (b) phase.
(a)
(b)As shown in Figure2(a), the maximal values of currents are on L1 and the minimal values of currents are on L2. The conclusion agrees with the analytical results based on physical optics method (PO) [21–24]. The normals of the flank on L1 and L2 are n^L1=cos⁡αx^+sin⁡αz^ and n^L2=cos⁡αy^+sin⁡αz^. Based on PO, induced current can be approximately expressed as J=2n^×Hi, so the currents on L1 and L2 can be derived as JL1=2Hm(-cos⁡αz^+sin⁡αx^) and JL2=2Hm(sin⁡αx^). The incident electric field is generally set as 1 V/m in computation; then the incident magnetic field is Hm=1/Z, where Z is the wave impedance of free space.Optical currents onL1 computed by PO method are compared with the values computed by FE-BI-MLFMA (called full-wave currents for short), as presented in Figure 3 (θ=8°). Figure 3(a) shows the amplitude (A/m2) of the currents, and Figure 3(b) shows the unwrapped phase (rad) of the currents. The range of horizontal axis indicates the distance between a point on L1 and the top of the coneFigure 3
Amplitude and phase of the currents onL1. (a) Amplitude; (b) unwrapped phase.
(a)
(b)Figure3(a) shows that there are significant differences between the optical currents and the full-wave currents. The full-wave currents along L1 are more complex than that of the optical currents. Near the two ends of L1, the amplitude of the full-wave currents is obviously larger than that of the optical currents. Within the internal region of L1, the amplitude fluctuates as a combination form of periodicity and linearity. Figure 3(b) shows that the two results agree well with each other, and the phase of induced current linearly changes along L1.The generatrix inxoz plane in shadow region is denoted by L1′. The induced currents on L1′ are shown in Figure 4. The fluctuation feature and the values of the induced currents on L1′ are similar to the case of L1, which means the scattering components of creeping waves are the denominated scattering contributions when the elevation angle of the LOS is small (θ=8°<α). The range of horizontal axis indicates the distance between a point on L1′ and the top of the cone.Figure 4
Amplitude and phase of the induced currents onL1′ in shadow region. (a) Amplitude; (b) unwrapped phase.
(a)
(b)Based on the characteristics of the induced currents onL1 and L1′ under different LOSs, the full-wave currents are described by the following expressions:(1)Jρ,f,θ=Amρ,f,θexp⁡jϕmρ,f,θ(2)Amρ,f,θ=am1sin⁡am2θfρ+am3θρ+am4θf+am5θ+am6ρ+am7(3)ϕmρ,f,θ=am8θfρ+am9θρ,where (1) describes the complex value of induced current; its amplitude and phase are indicated by Am and ϕm. Equation (2) is proposed according to the amplitude characteristics of the induced currents on L1 and L1′ under different LOSs. Equation (3) is proposed according to the phase characteristics of the induced currents on L1 and L1′. ρ indicates the position of the current element on L1m=1 or L1′m=2. Am(·) and ϕm(·) are the amplitude and the phase of the currents, respectively. amn (n=1,2,…,9) indicates the parameters to be estimated. amn is irrelevant to frequency, but am2, am3, am4, am5, am8, and am9 are linear-circular correlated with θ.To show characteristics of the currents on the whole conical surface, we slice the cone along lineL1 and unfold the flank to a fan surface, as shown in Figure 5. The currents on the fan and the bottom are shown in Figure 6. The half-angle of the fan is γ=πR/l, where l is the outer radius of the fan.Figure 5
Unfolded fan.Figure 6
Currents on the fan and the bottom.Based on PO method, the two-dimensional currents can be expressed asJ(ρ,β)=2n^(β)×Hi=2Hmsin⁡αx^-cos⁡αcos⁡β/γπz^, where (ρ,β) is the coordinates of the individual current on the fan. The component of J(ρ,β) along vertical direction (v^=-cos⁡αz^+sin⁡αx^) can be expressed as (4). Therefore, the change of currents with β can be described by the currents on L1 and the dependant function Fβ.(4)Jvρ,β=Jvρ,0Fβ(5)Fβ=cos2⁡αcos⁡βγπ+sin2⁡α.Equations (4) and (5) describe the fact that the change of induced currents with β can be expressed by the currents on L1 and the dependent function F(·).


## 2.2. Characteristics of Induced Currents on the Bottom

Similarly, the induced currents on the bottom are investigated here. The amplitudes and phases of induced currents on the bottom are presented in Figures7(a) and 7(b), respectively. The parameters of the incident waves are the same as those of Figure 2. The line in xoz plane along +x axis is denoted by I1, the line in xoz plane along -x axis is denoted by I1′, and the line in yoz plane along +y axis is denoted by I2. The maximal values of currents are on I1 and the minimal values of currents are on I2. The scattering field from the bottom is mainly contributed by the currents on I1 and I1′.Figure 7
Currents on the bottom. (a) Amplitude; (b) phase.
(a)
(b)The currents onI1 are presented in Figure 8. The range of horizontal axis indicates the distance between a point on I1 and the point on the curved edge. From Figure 8(a), it can be seen that the amplitude of full-wave currents fluctuates as a combination form of periodicity and linearity. From Figure 8(b), it can be seen that the phase of induced current changes linearly along I1. The characteristics of the currents on I1′ are similar to I1.Figure 8
Currents onI1. (a) Amplitude; (b) unwrapped phase.
(a)
(b)These characteristics of the currents on the bottom are similar to the case of currents on the flank. Therefore, the currents on the bottom are also described by the expression of (1), only with different parameters. The currents on I1 and I1′ are indicated by the case of m=3. Similarly, the two-dimensional currents on the bottom are also expressed by (4).


## 2.3. Estimation of the Current Model Parameters

The coefficients in the presented formula,am1~am9, are estimated through matching the currents computed by the model with the currents computed by FE-BI-MLFMA. The used method for optimal estimation is genetic algorithm (GA) [25–28]. The currents on the flank computed by the model with estimated parameters (the simulated currents for short) are presented in Figures 9(a) and 10(a). For comparison, the currents computed by FE-BI-MLFMA (the real currents for short) are presented in Figures 9(b) and 10(b). The used parameters in computation aref = 1~2 GH and ∆f=0.125GHz, θ=8°. The range of horizontal axis indicates the distance between a point on L1 and the top of the cone.Figure 9
Amplitudes of the currents. (a) Simulated currents; (b) real currents.
(a)
(b)Figure 10
Unwrapped phases of the currents. (a) Simulated currents; (b) real currents.
(a)
(b)The errors between the simulated currents and the real currents are examined: the correlation coefficient of the amplitudes is92.81% and the root mean square error (RMSE) is 6.01×10-5V/m2; the correlation coefficient of the unwrapped phases is 98.72% and the RMSE is 2.99 rad.Other cases of incident angle (θ = 1°~10°) are also investigated. The parameters of currents are estimated with the same method. The errors between the simulated currents and the real currents are the correlation coefficient of the amplitudes is larger than 90% and the RMSE is less than 10-4V/m2; the correlation coefficient of the unwrapped phases is larger than 95% and the RMSE is less than 4 rad.


## 3. Approach of Scattering Center Modelling

There are two classical methods for scattering center modelling. Scattering centers can be extracted directly from radar images by peak or pattern detection, or they can be described by parametric models and the unknown parameters of the model then be acquired through optimal estimation. The second method is not limited by radar resolutions and could provide more precise and elaborate descriptions of scattering centers, yet it is under the strict precondition of the accuracy and application scope of parametric model. The commonly used parametric models are proposed according to the solutions of several simple targets based on the physical optics method (PO) or geometrical theory of diffraction (GTD). Due to limitations of these asymptotic methods, the existing parametric models fail in some cases, for example, creep waves. In this paper, the parametric model is derived from the induced currents which are computed by a full-wave numerical method. All the scattering contributes are included in the derived model.Based on the electric field integral expression under the far-field condition in the free space [29], the scattered field can be expressed as(6)Evs=-jωμNv=-jkZNv,where Z is wave impedance of free space, Nv is the component of the vector N along v^ direction.(7)N=e-jkr4πr∫Jr′ejkr′·r^dS′,where r^=sin⁡θcos⁡ϕx^+sin⁡θsin⁡ϕy^+cos⁡θz^ and r′=ρsin⁡αcos⁡β/γπx^+ρsin⁡αsin⁡β/γπy^+h-ρcos⁡αz^.Substituting (4) into (7), we can derive the results of the currents in lit region (β = 0~γ).(8)Nv=e-jkr4πr∫l0∫0γJvρ,0Fβexp⁡jkcos⁡θh-ρcos⁡α+ρsin⁡θsin⁡αcos⁡βγπρdρdβ.The integral part related toβ can be derived based on the principle of stationary phase [30], as follows:(9)Gρ=∫0γFβexp⁡jkρsin⁡θsin⁡αcos⁡βγπdβ≈2πjkρsin⁡θsin⁡αexp⁡jkρsin⁡θsin⁡α.Then, the scattered field can be expressed as(10)Evs=-jkZe-jkr4πr∫l0Jvρ,0exp⁡jkcos⁡θh-ρcos⁡αGρρdρ=-jkZe-jkr4πr∫l0Jvρ,0exp⁡-jkr′·r^Vρ1/2dρ,where V=2π/jksin⁡θsin⁡α and r′ is the location vectors of currents on L1.From (10) we can see that the characteristics of the scattering fields from the flank in lit region are mainly determined by the currents on L1. Similarly, the characteristics of the scattering fields from the flank in shadow region are mainly determined by the currents on L1′, and the characteristics of the scattering fields from the bottom are mainly determined by the currents on I1 and I1′. The total fields are the supposition of the three components, as given below.(11)Evs=∑m=13Ema1m,…,a9m,where E1 indicates the electric field induced by the currents in lit region; E2 indicates the electric field induced by the currents in shadow region; E3 indicates the electric field induced by the currents on the bottom.Based on (1) and (10), we can derive the analytic expression of the scattering center model at VV polarization. VV means the incident wave and reflected wave are all in vertical polarization. The analysis of the scattering centers corresponding to the currents on L1 is provided. The analysis processes are similar for the cases of L1′, I1, and I1′, so the repeated contents are ignored here.Substituting the expressions of currents onL1, L1′, I1, and I1′ into (10), we can derive the total scattering fields, as given by (A.1) in the Appendix. Then through the regular computation of the integral along ρ in (A.1), we can get 5 exponential terms as shown in (12); the five functions are the five corresponding parametric models of scattering centers contributed by the currents on L1.(12)pSC1=eja13+a19l+a15a11l2a12+a18f-2kr′·r^+a112ja12+a18-kr′·r^2f2eja12+a18l-kr′·r^+a14fpSC2=e-ja13-a19l+a15a11l2a12-a18f+2kr′·r^-a112ja12-a18+kr′·r^2f2e-ja12-a18l+kr′·r^+a14fpSC3=eja15a112ja12+a18-kr′·r^2f2eja14fpSC4=-e-ja15a112ja12-a18+kr′·r^2f2e-ja14fpSC5=eja19lja18f-kr′·r^-la17-l2a16+a17+2la16ja18f-kr′·r^-2a16ja18-kr′·r^f2eja18lf-kr′·r^.According to the numerical results of the parameters, the components,pSC1, pSC3, and pSC4, whose amplitudes are quite small compared with the others, can be ignored in the following analysis. The components, pSC2 and pSC5, include the contributions of four scattering centers. The four dominated scattering centers, denoted by SC1, SC2, SC3, and SC4, are illustrated in Figure 11. The relationship of the parameters of the model with the amplitudes and locations of scattering centers is listed in Table 1.Table 1
Associations between the current model parameters and scattering centers parameters.
Scattering center
Position
Amplitude
S
C
1
a
12, a14, a18, a22, a24, a28
a
11, a12, a13, a15, a18, a19, a21, a22, a23, a25, a28, a29
S
C
2
a
18
a
16, a17, a18, a19
S
C
3
a
28
a
26, a27, a28, a29
S
C
4
a
32, a34, a38
a
31, a32, a33, a35, a38, a39Figure 11
Positions of four dominated scattering centers on the target.S
C
1 locates at the top of the cone, which is mainly contributed by the components pSC2 of L1 and L1′; SC2 locates at edges of the bottom in lit region, which is mainly contributed by the components pSC5 of L1 and I1; SC3 locates at edges of the bottom in shadow region, which is mainly contributed by the components pSC5 of L1′ and I1′; SC4 locates outside the geometry of the cone, which is mainly contributed by the components pSC2 of I1 and I1′, that is, the creeping waves in the shadow region.


## 4. Validation of the Model

According to the principle of RP, the intensity and position of the peak reflect to a certain extent the scattering amplitude and downrange location of the corresponding scattering center. The accuracy of the derived scattering center model is validated through comparing the RP simulated by the model (the simulated RP for short) and the data simulated by the scattered fields computed by FE-BI-MLFMA (the real RP for short). Two cone-shaped targets are investigated here. The geometry of target I is the same as the case of Figure1. The radiuses of bottom disks and spherical tops of target II are the same as target I. The height of target II is h=5.4m. Two cases with different parameters setting are investigated.Case 1 (θ=0, ϕ=0, f=1～2 GH, △f=0.125GHz).
The simulated RP and the real RP are presented in Figure12(a). For demonstration of the superiority of this model, the RP obtained by the attributed scattering center model through traditional method are also compared with the simulated RP, as presented in Figure 12(b). The errors of amplitude and locations of the scattering centers are listed in Table 2.Table 2
Errors of amplitude and locations of the scattering centers.
Case1
Position errors
Amplitude errors
Parameters
S
C
1
S
C
2
S
C
3
S
C
1
S
C
2
S
C
3
Current model
0.19
λ
0.13
λ
0.03
λ
0.14 dB
0.02 dB
0.04 dB
Attributed scattering center model
0.19
λ
0.13
λ
Fail
0.15 dB
0.04 dB
FailFigure 12
Comparison with the RP. (a) Comparison of the simulated RP with the real RP; (b) comparison of the RP obtained by the attributed scattering center model with the simulated RP by this model.
(a)
(b)The comparison results demonstrate that the simulated RPs by this model agree well with the real RP. Superior to the exiting model, the scattering centerSC4 induced by creeping waves is precisely characterized by this model.Case 2 (θ=0～10∘, Δθ=0.1°, ϕ=0, f=1～2 GH, Δf=0.125 GHz).
The two-dimensional image of the RPs under the continuously changedθ (the RP history for short) is investigated. The simulated RP history and the real RP history of target I are presented in Figures 13(a) and 13(b). The correlation coefficient of the two images is 86.3%. The RMSE is 2.56 dB.Figure 13
RP history of target I. (a) Simulated RP history of target I; (b) real RP history of target I.
(a)
(b)The simulated RP history and the real RP history of target II are presented in Figures14(a) and 14(b). The correlation coefficient of the two images is 88.7%. The RMSE is 2.21 dB.Figure 14
RP history of target II. (a) Simulated RP history of target II; (b) real RP history of target II.
(a)
(b)From the images of the RP history of two targets, it can be seen that the simulated results by this model agree well with the results obtained by the full-wave method. The four dominated scattering centers under different aspect angle of LOS are all precisely characterized. Whenθ is within half conical angle, the SC4 induced by creeping waves in shadow region is a dominated scattering center; its intensity is even larger than the SC2 and SC3 on the edges. Therefore, in these cases, SC4 should be considered in scattering center modelling.


## 5. Conclusions

The induced currents on cone-shaped targets are investigated in detail in this paper. It is found that scattering characteristics of the cone-shaped targets are mainly determined by one-dimensional currents on the lines in the plane constructed by the LOS and the rotational axis of the cone. The parametric model of the currents on the lines is built based on their distribution characteristics. Based on the induced currents computed by FE-BI-MLFMA, the parameters of the model are estimated using GA. Then, based on the current model, the scattering center model is derived according to electric field integral expression under far-field condition. This model is suitable for cone-shaped targets with different cone angle and height and can precisely describe scattered waves under arbitrary incident directions.The precision of this model is validated by the RPs under different aspect angles of LOS. The RPs simulated by this model agree well with the RP simulated by the scattered fields computed by FE-BI-MLFMA. The location errors of scattering centers are less than half wavelength, the amplitude errors are less than 1dB. The comparison results between this model and the attributed scattering center model show that this model is of higher accuracy, especially in describing the scattering centers contributed by creeping waves.



---

*Source: 7482895-2017-03-07.xml*