A Frequency Selection Method Based on the Pole Characteristics

Due to the heavy jamming band of high frequency, frequency selecting strategies are serious issues for the system designed to achieve its best performance. Pole is independent of the direction and polarization of the incident wave, but the residue corresponding to the pole is related to the direction and polarization of the incident wave. And the value of residue is proportional to the value of the pole. This paper chooses the frequency which can maximize the residue in the high-frequency band as the optimal frequency for accurate extraction.The simulation result of a large number of ship targets shows remarkable rise in average recognition rate by using this method, compared with the average recognition rate of randomly selected frequency.


Introduction
With the extensive application of high-frequency radar, radar frequency selection has become an important research topic [1][2][3][4][5].The propagation characteristics of the high-frequency electromagnetic waves are not only a function of the distance but also a function of frequency [4].In the case of the same distance, the higher the frequency, the faster the attenuation of electromagnetic waves.The frequency selection will seriously affect the performance of the high-frequency surface wave radar on target detection.It is necessary to optimize the selection of operating frequency of the high-frequency radar in order to get excellent detection performance of high-frequency radar.This paper uses pole characteristic to select frequency, explore criteria for frequency selection optimization, and get the selected frequency optimization method.Finally, the result of simulation for a variety of ship target detection verified the correctness and validity of the method.

Frequency Optimization Criterion
To get an optimal operating frequency, the first thing to consider is the frequency optimization criterion.
A high-frequency radar operating frequency can be obtained through scientific quantitative calculation in accordance with this criterion so that the high-frequency radar can achieve optimal detection performance.Radar is able to find the goal, that is to say, the average correct recognition rate to meet certain requirements.Paper [6] puts forward that radar working frequency optimization selection is investigated for an optimization selection method based on the minimum classification error rule in order to improve airplane target identification performance and that MER can change to the principle of average correct recognition rate maximum.Optimal operating frequency is the one which makes the average correct identification rate maximum.This also can be called as the principle of average correct recognition rate maximum.The average correct recognition rate can be got from dividing the number of correctly identified test by the total number of randomized trials.
Pole is independent of the direction and polarization of the incident wave, but the residue corresponding to the pole is related to the direction of polarization of the incident wave.And the value of residue is proportional to the value of the pole.Considering the above two factors, select the frequency which corresponds to the maximum value of residue as the optimal frequency.

Calculation of Residue
Matrix Pencil Method [7] is proposed on the basis of Pencil of Function [8] and has a higher computational efficiency.It only needs two Hankel matrices constructed by time-domain data sequence of the target.Then seeking the generalized eigenvalue of these two Hankel matrix we can calculate the poles of target.Because Matrix Pencil Method needs not to calculate the root of the characteristic equation as Prony method, Matrix Pencil Method has high antinoise performance.To improve the anti-noise performance of Matrix Pencil Method, Sarkar first applied Matrix Pencil Method in matrix singular value decomposition [9] and gave a method to decide the number of poles, through the matrix low rank approximation which could well inhibit the effects of noise, and the efficiency of the algorithm is improved greatly.But Matrix Pencil Method depends on the selection of relevant parameters; Sarkar studies the beam parameters selection [10] and the method for determining the number of target poles [9].Matrix Pencil Method has high computational efficiency in certain SNR and good operation steadiness, which is very suitable for the extraction of the target pole.
From previous research we can see that the properties of Matrix Pencil Method are superior than others among a large number of pole extraction methods [6].But there are some disadvantages, such as the unknown number of poles that will to a large extent affect the matrix method of pole extraction results.We use the improved the overall minimum square Matrix Pencil Method to extract the pole of complex targets [5,9].

Principles of Matrix Beam Method.
In 1971, after a lot of electromagnetic pulses tests, Baum found that, when objects are excited by electromagnetic pulse, their late response is a series of superposition of attenuation sinusoidal oscillation, and that attenuation factor and oscillation frequency composed the complex natural resonant frequency, which is the pole of target.Baum puts forward the theory of singularity expansion method on the basis of the above phenomena [11,12].
are eigenvectors corresponding to eigenvalue   .Equation (10) can also be written as where superscript " †" represents the generalized inverse of the matrix.It can be seen from ( 11) that it is possible to obtain the required target pole when eigenvalue of Y † 1 Y 2 is a known value.
And then, target complex residue can be obtained through the least squares method: where

The Choice of Pole Number.
Advanced transient response of fully conducting targets can expand into an infinite number of decaying sinusoidal vibrations and forms.And formula (1) constrained the complex attenuation frequency on the .
It is One of the important questions about extracting pole of complex targets is how to select the appropriate  [13].If the  value is selected too small, it will miss the target's true pole, and the pole will have a large deviation; if the  value is selected too large, it will not only produce false pole but also will lead to true pole position moved.
Pole number of simple targets within a certain frequency band can be easy to get; but for the aircraft and ship target, the pole number  contains in the signal is unknown.For this issue, this paper uses the selection method based on timedomain signal reconstruction error minimized to decide the number of poles [5,14].
In order to get better results of pole extraction, we can apply the singular value decomposition technique to pole extraction problem with matrix beam methods.Let us suppose formula (1) is accurate; if we ignore deviation and noise which come from modeling and calculations, matrix Y will become a full rank matrix.
By singular value of formula (3), we can get this: 1 ,  2 , . . .,   is the biggest  singular values of Y.If there is no noise,  +1 ,  +2 , . . .will be zero.Take low rank of Y to approximation of matrix Ŷ: And we can remove the last column of Ŷ, then we will get matrix Ŷ1 , and if we remove the first column, we will get matrix Ŷ2 like this Ŷ1 , Ŷ2 singular values can be decomposed as removes the last column of matrix V H will get V H 1 , and if remove the first column, we'll get V H 2 .Theorem [15] If matrix Y's rank greater than or equal , then there exists the only matrix of the same dimensions Ŷ, and the rank is equal to , which makes the norm of ‖Y − Ŷ‖  Frobenius minimum.The extent of best approximation is described by ‖Y − Ŷ‖  = ∑ − =+1   , where ‖ ⋅ ‖  indicates Frobenius norm.
Theorem explains that the use of Ŷ is the optimal approximation to Y under F-norm in theory.We will get Ŷ1 and Ŷ2 through Ŷ; then the eigenvalue of Ŷ † 1 Ŷ2 will be got.This method can effectively inhibit the effects of noise on pole extraction.This singular value decomposition technique applied to a method of matrix beam method is called total least square matrix beam method [5].
Theorem in the assumption of formula ( 1) is accurate and on the condition of  is known.In the actual pole extraction,  is unknown, and it cannot be determined only through the distribution of singular values  1 ,  2 , . . .,   ,  +1 ,  +2 , . ... Because of the distribution of singular value, neither  is rendered when you select a value, then the singular value is far smaller than the preceding singular value, nor  appears when you select a value, followed by singular value of size are essentially the same.Therefore you need to focus on pole number  of selection method.
It can judged by pole extraction results to decide the parameters selection is good or not.We will use the deviations between theoretical and extraction as criteria to constrain parameters selection and to get the best parameters .But the most target theoretical pole is unknown, so it needs to give a reasonable and actionable criterion.
If the target number of pole  is known, pole   can be achieved by matrix beam method, and the complex residue of target   can be got by arithmetic expression (4).It also can reconstruct the time domain signal  rec () of target, then we will use the error between the reconstructed domain signal  rec () and domain signal  cal () As a criterion min Choosing the appropriate parameter  opt .
Cause  can meet And usually target scattering data length  is not too long; the value of  is generally around /3.So the value of  is an integer from 1 to /3.We can directly choose all parameters  within the interval range to calculate () select minimal parameter () which is corresponding to  as the most appropriate parameters  opt .Figure 1 shows that the ship target error () and pole number determine the results of process variation of () with parameter .It can be seen that as  is increasing, errors it first shows shock reducing, then minimum is reached, and at last it will be larger again.
We can see from Figure 1 that in this paper, the best number of pole is  = 46.This approach can reduce the impact because models are inaccurate, and there are deviations in the data and the noise, extracting the precise targets pole.

The Rationality of the Residue Selected Frequency.
To prove the rationality of the residue selected frequency, we need to verify that the residue is a main contribution to RCS.Because the residue is obtained from the late response, we need to compare the contribution of the early and late response to the RCS for the first step.
where  ear (),  lat () represent the early and late responses.And  lat () is the () in formula () = ∑  =1       ; we can rewrite formula () = ∑  =1       as In theory, no matter what form of target transient response  all () is in formula  all () =  ear () +  lat (), it can be expanded into an infinite sum number of complex exponential where R is the coefficient of expansion corresponding to exponential term  s  .Basis function set {s 1 , s2 , . ..} definitely contains the target pole set { 1 ,  2 , . . .,   }.Although using only the pole set { 1 ,  2 , . . .,   } to expand the transient response  all () may exist certain error, but that error is admissible.
where R can be called the new residue corresponding   to different from   .

The Comparison Simulation.
Here we get the original RCS by using software FEMAP to model the ship target and FEKO to calculate the RCS.The target we choose 49.05 meters long, 12.40 meters high, and 56.39 meters wingspan airplane, the model of airplane, is shown in Figure 2. The comparison result is shown in Figures 3 and 4, red line represents recovered RCS only by late time response and blue line shows the original RCS, which proves that the late response is a very important contribution to RCS.The above conclusion implies that the residue calculated from late response is also very important to RCS.Next, we will test and verify the above conclusion by estimating the contribution of residue to RCS.
To make the simulation contain common complex goal, this time I use software FEMAP and FEKO to model one ship and calculate its RCS, the model is shown in Figure 5. Then calculating the residue according to the method I described in the paper, incident angle ranges from 0 ∘ to 180 ∘ at a 10 ∘ interval.The frequency selection is from 0.2 MHz to 20 MHz at a 0.2 MHz interval.The number of the pole is 46.The result is shown in Figure 6, We can see that the residue changes with frequency and angle.It shows peak appeared when angle is 90 ∘ and frequency is 9.6 MHz.
We can set that maximum residue as zero and then recover the RCS.Comparing the variation of RCS before and after.Figure 7 shows that when maximum residue was replaced by 0, RCS declined to 8 dB-10 dB; the above simulation proves the influence of residue on the value of RCS.The simulation results provide a theoretical basis for using residue to select frequency.high-frequency band into  frequency points with a fixed interval, obtain the RCS by FEKO calculating in all angle, we can get a target database in high-frequency band:

Compare the Average
Extracted  frequency points from  frequency points, are equivalent to extracted  column vectors from  column vectors in above matrix, to constitute a submatrix: D  is a radar target feature vector; name D as radar target feature space composite of  types target feature vector.
Obtain the RCS on  frequency points by actual measurement, write it as name it radar target measurements eigenvectors.
According to the existing research Ksienski and Lin show that the nearest neighbor classification method for the class of complex radar target aircraft has better recognition results [11].The Euclidean distance between measuring characteristic feature and the  types target in target feature space is According to nearest neighbor classification, if target is  type.

Simulation.
In order to study the effect of the maximum number residue frequency-selective method of ship target also verifies that the maximum residue frequency selection method can effectively improve the recognition effect.I do a large number of random tests to verify the correctness of the theory and the effectiveness of the method.
The simulation selected 6 typical classes of ship target as a target to be classified from the library of high-band ship target characteristics; the type of ship targets and basic dimensions are shown in   Assume that select the  target, add Gaussian white noise with zero mean; mean square deviation is , and signal-tonoise ratio is calculated using Here does not consider the situation that the target is not included in the target class to be classified.1000 randomized tests are conducted to study the changes of the correct identification rate in different signal-to-noise ratios with the different frequency selection method.Figure 8 shows average correct recognition rate comparison between maximum residue frequency selection method and random frequency selection method of 6 goals.Apparently, the average correct recognition rate of maximum residue selection frequency method is better than random frequency selection method.The former method is about six percent higher than the latter.The average correct recognition rate of maximum residue selection frequency method can be as high as 0.8.
To make the simulation much more closer to real conditions, we need to consider that there may be deviations between the measurement RCS and calculated RCS from simulation.That deviations may come from the bug of software or the noise that affects the measurement of RCS.In this simulation, a larger amount of uniformly distributed random value is added to calculated RCS as that deviation.Taking into account that the RCS deviation may occur in the reality as  dB, to  target's  frequency, we add uniformly distributed random value in the range of [−Δ  , Δ  ], satisfying  = 10 ⋅ log (Δ  /  ).At the same time, Gaussian white noise we discuss in Figure 6 is still added (calculated RCS with noise).Figure 9 shows the recognition results of 6 ship targets at two operating frequencies 13.4 MHz and 15 MHz when taking consideration of the deviations between the measure RCS and calculated RCS.
From Figure 8, Compare the average correct recognition rate of the above method and the random frequency selection method when angle is 90 ∘ ; average correct identification rate of my method is 82.76% when SNR is 13 dB, while the average correct identification rate of the random frequency selection mode is 61.33% when SNR is 13 dB, which verifies the effectiveness of the method.
In Figure 9, we can see that the average correct recognition rate declined quickly when measuring deviation becomes larger.When measuring deviation is ±1 dB, the average correct recognition rate can reach 0.87.When measuring deviation is ±2 dB, the average correct recognition rate can reach 0.68.Comparing Figure 9 with Figure 8 we can that the measuring deviation does have an influence on the average correct recognition rate.But it can be overcome by adding more frequency.The simulation of 4 frequencies 13.4 MHz, 15 MHz, 16.6 MHz, and 20 MHz is shown in Figure 10, from which we can see that the average correct recognition rate can reach 0.85.
From above comparison we can get the conclusion that the measuring deviation has a great influence on recognition International Journal of Antennas and Propagation  rate.But this problem can be overcome by adding more frequency even when the measuring deviation exists, and the average recognition rate of the method adopted in this paper is significantly larger than random selected frequency method, which indicated that the robustness of the method is absolutely better than random selected frequency method.And the method is feasible.

Conclusion
This paper studies the frequency selection method in highfrequency band radar target recognition.Based on the average correct recognition rate maximization principle, considering the residues corresponding to pole is related to the incident wave direction and polarization; the value of residue is proportional to the value of the pole, a maximum residue frequency selection method was proposed, and satisfactory recognition result is obtained with the simulation experiment on multiclass ship targets.Considering the difference between real measured RCS, and software calculated RCS, we add measurement deviation in software calculated RCS, and the ship target recognition effect under the conditions of different measurement deviations is analyzed.This verifies the correctness and validity of the frequency optimization method.

Figure 1 :
Figure 1: The error is changing with the value of  in different.

4. 1 . 1 .
How to Recover RCS by Using Late Response Only.Target time-domain transient response and frequency response are the Fourier transform for each other.If we can recover the time-domain transient response, the frequency domain response of the target can be obtained naturally.Target timedomain transient response is composed of early and late responses, which can be written as  all () =  ear () +  lat () ,

Figure 3 :
Figure 3: The comparison between original RCS and the RCS recovered only by late time response in frequency domain.(a) The angle between the incident and the bow is 0 ∘ .(b) The angle between the incident and the bow is 180 ∘ .

Figure 4 :
Figure 4: The comparison between original RCS and the RCS recovered only by late time response in time domain.(a) The angle between the incident and the bow is 0 ∘ .(b) The angle between the incident and the bow is 180 ∘ .

Figure 7 :
Figure 7: Comparison of RCS between original and residue = 0.

Figure 8 :
Figure 8: Average correct recognition rate changes with the signalto-noise ratio by different frequency selection methods for 6 classes of target.

Figure 9 :
Figure 9: The average correct identification rate changes trend of 6 targets with different SNR when deviations exist.(a) The measuring deviation is ±1 dB.(b) The measuring deviation is ±2 dB.

Table 1 ,
Including aircraft carriers, cruisers, destroyers, frigates, hunting submarines, missile boats, and different types of ships.Here added noise to calculated target RCS to simulate the measured data, cause many of them are secret we cannot get the measured data.Here we can name this RCS calculated RCS with noise to different from below.