Passive Radar Source Localization Based on Cuckoo Search NVTDOA

In order to meet the requirement of passive radar source localization in electronic warfare, the concept of the virtual time diﬀerences of arrival (VTDOA) is proposed by taking advantage of the characteristics of the same UAV in diﬀerent positions at diﬀerent times and the periodic rotation of radar pulse signal. The VTDOAs are the passive localization information deﬁned as the time diﬀerences of the radar pulse transmission from the radar position to diﬀerent virtual receivers. Firstly, a nonlinear VTDOA (NVTDOA) localization model is constructed. Moreover, suﬃcient conditions for accurately calculating the periodic integers in the model are analyzed, and the observability conditions of the localization model determined are deduced. Secondly, the convergence solution of the NVTDOA localization equation is obtained by Cuckoo search algorithm; thus, passive radar source localization is realized. Finally, the performance of the proposed method is veriﬁed by comparing with the existing methods.


Introduction
Passive localization is a fundamental problem that usually accomplishes localization task using the radiation signals of the target without radiating electromagnetic wave signals outwards, which has been widely used in many fields such as radar [1], sonar [2], sensor network [3], and wireless communication [4]. e passive localization method belongs to a two-step localization system. In the first step, some signal parameters in the target source transmission are extracted, such as received signal strength [5], time difference of arrival (TDOA) [6,7], and angle of arrival (AOA) [8,9]. In the second step, the target position is estimated by solving a set of nonlinear localization equations.
For a long time, the passive localization method based on TDOA has been an important research topic [10]. For noncooperative radiation sources, it is particularly important to study TDOA localization technology for the single station because the time of signal emission is unknown. ere are three main calculation methods of TDOA for noncooperative radiation sources. e first is to take the rising edge of the first pulse as the arrival time [11]. e pulse amplitude jump points from scratch are chosen as the reference point of the method, so it is easier and more accurate to determine the reference point.
However, the signal amplitude of the first pulse is small, and the estimation error becomes larger when the pulse loss is more serious or the amplitude is smaller. e second is to take the average arrival time of all received pulses as the arrival time of radar signals [12]. e method has high stability but high computational complexity. If the received pulses overlap, the estimated performance will deteriorate. e third is to select the arrival time of the pulse with the highest amplitude in the main lobe pulse as the receiving time of radar signal [13]. When the main lobe beam is narrow, the arrival time will have higher stability. But if the main lobe beam is wider, there will be more maximum amplitude values in the pulse, which will degrade the estimation performance. In order to overcome the shortcomings of the above TDOA estimation, a norm approximation method was proposed by Wan [14]. e approximate mathematical expression of the main lobe beam pattern of the radar antenna is obtained at the receiving end, and the time point corresponding to the maximum value is used as the arrival time of radar signal to complete the estimation of VTDOA for radar signal.
With the more and more extensive application of the unmanned aerial vehicle (UAV) in battlefield, there are some related research studies on the algorithm of using a UAV receiver as moving single station to locate a radar emitter. A passive location method of pulse information sources in the battlefield by using a single UAV was proposed by Nickens and Darko [15]. e method takes the rising time of pulse signals received by the UAV as the signal arrival time and builds localization equations based on UAV's movement trajectory. However, when the pulse amplitude of the source is small or far away, the judgment error of the receiver on the rising moment of the pulse will be larger, and the localization accuracy will be worse. Kim et al. [16] proposed a way for a single drone to use a known database to provide rough location services for targets in the battlefield when the satellite positioning system is not working properly. To improve the performance of the location system in scenarios with a small number of receivers or when there are technical restrictions to deploy them in LoS condition, in reference [17][18][19][20], a localization algorithm based on time difference of arrival (TDOA) virtual sensors is proposed, where the reflection points are regarded as TDOA virtual sensors. And the algorithm is applied to localization of UAV in urban scenarios, ground vehicles, emitters, and RF emitters. Based on this, Hao et al. [21], Wan [14], and Zhu et al. [22] put forward a radar radiation source localization method based on passive synthetic aperture antenna array. e method utilizes the characteristics of a single-station small-load UAV in different positions at different times and the periodic rotation characteristics of radar pulse signals and uses VTDOA to complete the localization task. However, the observability of the localization model is not discussed, and the applicable conditions of the model are not analyzed in [14,21,22]. Moreover, it is assumed that the two radar signal waves received by the UAV are the same plane waves to obtain the closed solution of the positioning equation, and the radar rotation speed at the first received signal is used to calculate the radar antenna rotation time difference. erefore, if the UAV runs at a high speed or the localization time difference is long, the localization accuracy of the target radar will become very poor.
For the deficiency of the above algorithm, based on the motion characteristics of UAV and the periodic rotation characteristics of radar pulse signal, a new method of passive localization for the radar emitter based on Cuckoo search NVTDOA is proposed. In conclusion, the Cuckoo search NVTDOA method has made contributions in two aspects compared with the literature. First, the NVTDOA model is not sensitive to the speed of UAV or the localization time difference, so that it has more applicability. However, the approximate calculation of the rotation time difference in the VTDOA model makes its accuracy more sensitive to these, so VTDOA is only applicable to the positioning of UAV at low speed and short stroke. Second, the Cuckoo search algorithm is used instead of the closed-form solution algorithm to solve the localization equation, which finally improves the localization accuracy of the algorithm. In Section 2, the traditional VTDOA localization model is introduced. In Section 3, first, the shortcomings of the VTDOA localization model are analyzed in the first paragraph, so an improved NVTDOA localization model is proposed considering the integral characteristics of radar antenna rotation time difference. Second, the NVTDOA localization model is described according to the traditional two-step localization method. e first step, the derivation of the NVTDOA localization model, is given in Section 3.1, and the second step, the estimation method of virtual time difference in the NVTDOA localization model, is introduced in Section 3.2. erefore, it provides a localization model for the construction of the following localization equation. In Section 4, the observability of the NVTDOA model is introduced to verify the rationality of the model, and the theoretical basis is provided for the approximate calculation of the localization equationin the followi. Based on this, the construction and solution of the location equation are presented in Section 5. To realize the location of a radar emitter, the location models of two flight segments are established in Section 5.1, and the two models are combined to generate the location equations. According to the nonlinear characteristics of the equations, in Section 5.2, the Cuckoo search algorithm is proposed to solve the equations and finally realize the location of the radar emitter. In Section 6, VTDOA and NVTDOA are used to analyze the location of three different passive radars to verify the high accuracy and antinoise of the NVTDOA algorithm. In Section 7, the theoretical derivation and simulation are summarized and analyzed.

VTDOA Localization Model
e angular velocity of the target radar antenna rotating uniformly in counterclockwise direction is Ω, and the periodic rotation angle is Π. Taking the t i moment of the UAV as the origin, a simple motion model of a UAV and the target radar is shown in Figure 1. e angle between the main lobe pulse direction of the far-field radiator and the UAV is φ i . If the UAV moves in a straight line, it can be regarded as a constant.
Assuming that the radar rotation period T � Π/Ω, which can be measured, the UAV receives the radar main lobe pulse signal at position S i at time t i , and after a period of time, the UAV moves to position S i+1 at time t i+1 and receives the radar main lobe pulse signal. Let P be the shortcut point of radiation source S 0 from the UAV track. e length from S 0 to P is D, and the length from S i to P is L. e distance from the initial point S i of the UAV to S 0 is d i and that from S i+1 to S 0 is d i+1 . Similarly, the distance from S i+2 to S 0 is d i+2 . e distance from location S i+1 to S i is l i , that from S i+2 to S i+1 is l i+1, and so on. Of these quantities, l i and l i+1 are known quantities. It can be seen from Figure 1 that, for the fixed radiation source S 0 , it can be characterized by X � D, L { } T . en, the measurement equation of the variables to be solved can be written as follows: where c � 3.0 × 10 8 m/s, which is the propagation speed of radar radio waves. Δt i+1,i is the difference of TOA measurement noise between two locations of UAV, which is modeled as Gaussian white noise with zero mean value and standard deviation of σ t . K i+1,i represents the integer of radar rotation period, and it means that the time interval between two UAV locations is K i+1,i radar period. Δψ is the rotation time difference of the radar antenna from S i to S i+1 .

International Journal of Antennas and Propagation
In the VTDOA localization model, it is assumed that the linear velocity of the radar antenna during the rotation from S i to S i+1 is the same as v ri � Ω · d i . ( en, the VTDOA localization model is that

NVTDOA Localization Model Based on Rotation Time
Difference of Radar Antenna. In [14,17,18], in order to find the closed solution of the localization equation, it is assumed that the two radar signal waves received by UAV are the same plane waves. e rotation time difference of the radar antenna is calculated by using the radar rotation speed at the first received signal. erefore, if the UAV runs at a high speed or the localization time difference is long, the localization accuracy of the target radar will become very poor. So, in this paper, an algorithm is presented to obtain the time difference of the radar antenna rotation based on definite integration of the rotation angle. As shown in Figure 2, Δψ can be obtained by integrating the rotation angle θ of the radar antenna. e rotation time difference of the radar from S θ to any small angle Δθ is Δψ θ ; then So Δψ can be expressed as According to the sine theorem, Δx can be obtained by where Substituting equation (8) into equation (5) yields where Δθ is an infinitesimal amount; then, sin Δθ ≈ Δθ, sin(Δθ + φ i + θ) ≈ sin(φ i + θ), and equation (9) can be reduced to Radar antenna rotates counterclockwise with an angular velocity of Ω R a d a r p u l s e w a v e  International Journal of Antennas and Propagation Substituting equation (10) into equation (6), we have Substituting equation (12) into equation (1) yields

Estimation
Method of NVTDOA. When the radar system conducts periodic scanning, assuming that the UAV starts to collect the data at time t 0 , the UAV will not always be within the range of radar antenna fan beam irradiation. At the ith rotation period of radar, with the antenna rotating, the UAV will enter the range of the antenna main lobe fan beam at time t i and leave the main lobe irradiation area after receiving N P radar main lobe pulses continuously. And the center time of the kth main lobe pulse received in the ith rotation period is t i,k . Under the influence of beam gain directionality G(t), N P received pulse signals have different amplitudes. So the baseband signal waveform of the radar radiation source pulse signal received by UAV is shown in Figure 1, and PRI � (1/PRF) which is the pulse time interval, where PRF is the radar pulse repetition frequency. So the baseband signal of the linear frequency modulated (LFM) radar pulse received by the UAV can be expressed as follows [22]: where ξ(t) is additive random noise, P T is the transmitting power of the radar antenna, G R is the transceiver gains, and φ c is the initial phase of the signal. As shown in Figure 3, the aircraft will receive N P main lobe pulses in each rotation cycle of the radar. It is necessary to mine the hidden useful information from the received signal for radar location. e hidden information used in the paper is defined as NVTDOA. NVTDOA is defined as the difference of transmission time required for the radar signal to reach the same receiver at two different positions. In this paper, NVTDOA is represented by Ψ j,i , and t i and t j represent the initial time for the same UAV to enter the ith and jth scanning periods of the main lobe of the radar antenna, respectively. So In practical application, it is difficult to accurately estimate t i and t j by x(t) of equation (14). As shown in Figure 3 where t i, max and t j, max are the peak times of G(t) of the ith and jth rotation period, respectively. In order to improve the estimation accuracy of NVTDOA, the estimation method of t i, max and t j, max in this paper is referred to [21], which will not be described here.

Observability Analysis of Localization Model
e localization model is obtained in the previous section; the next problem is in which case the equation is solvable, namely, the observability problem of the model [23,24]. In this paper, the observability conditions of the localization model are obtained by identifying the parameters of two systems.

Calculation Conditions for Integers of Radar
Rotation Period

Calculation Conditions with the Radar Rotation Period
Known. When the radar rotation period is known, the period integer can be calculated by the following equation, taking S i and S i+1 as examples: where round(·) denotes the rounding integer operation. Substituting equation (13) into equation (16) yields Sufficient conditions for accurate estimation which can be obtained by equation (17) are as follows: Condition (1a) has two meanings. e first is that d i+1 − d i ≪ cT. e radar rotation period is generally in the order of seconds, so condition (1a) is evidently valid. e second is that Δt i+1,i ≪ T, which is usually satisfied in the radar because the TOA measurement noise is typically on the order of ns and the radar radiation source period is typically on the order of s. e meaning of condition (2a) is that the angle θ i+1,i between two measurements of the UAV and target is less than 0.5Π, so that the integer operation can be avoided to estimate |K i+1,i | � ± i, where i is an integer. e rotation angle period of most target radars is 2π; that is, Π � 2π. When a UAV does not move radially relative to the target, we will have that θ i+1,i < π � 0.5Π, which satisfies the conditions.
It can be seen from conditions (1a) and (2a) that when the radar rotation period is known, K can be accurately estimated.

Calculation Conditions with the Radar Rotation Period
Estimated.
e radar antenna rotation period is unknown but can be estimated by the received radiation source pulse. It can be recorded as T � T + ΔT, and ΔT is the estimation error. We have e integral fraction expansion is as follows: According to equation (19), the sufficient conditions for accurately estimating K i+1,i are as follows: e meaning of conditions (1b) and (2b) is consistent with conditions (1a) and (2a); the meaning of condition (3b) is that the estimation deviation of the rotation period of the radar antenna is negligible compared to the period itself, which is also the basic condition for the accuracy of the period estimation.

Observability Conditions of the Localization Model.
e NVTDOA equation of the target radar position is nonlinear. In most papers, the observability theorem of Lee and Dum [25] is used to solve the observability condition of the equation. When the theorem is applied to discrete nonlinear systems, the equivalent conclusion is to judge the rank of the Gramer matrix. For a nonlinear system as shown in equation (20), let the system set be S and the Gramer matrix of the state variables x * i0 in the set is as in equation (21): where H i � (zy(n)/zx(n)), which is the Jacobian matrix, and Φ is a state transition matrix. If there is an integer N such that rank Γ(i, i + N − 1) � n (state dimension), then the system is completely observable for S [26].

Observability Analysis with the Radar Rotation Period
Known. When the radar period is known, the state transition matrix Φ i,i+1 � I 2×2 and the measurement equation of the variable X � D, L { } T can be simplified as e Jacobi matrix of NVTDOA j,i at X is From Figure 1, it can be seen that

International Journal of Antennas and Propagation
Under far-field radiation source settings, the radar source is a two-dimensional target. Let the position of the UAV receiving the radiation source pulse be S i , S i+1 , and S i+2 , and Gramer matrix based on NVTDOA by three measurements is where e determinant of equation (25) is that where Discussions about M i are as follows: (1) If l i � 0, then M 1 ≠ 0 when l i+1 ≠ 0 and d i ≠ d i+1 . But l i � 0 and d i ≠ d i+1 are contradictory, so this condition is not valid (2) If l i+1 � 0, then M 1 ≠ 0 when l i ≠ 0 and d i+1 ≠ d i+2 . But l i+1 � 0 and d i+1 ≠ d i+2 are contradictory, so this condition is not valid erefore, from (1) and (2), we can know M 1 ≠ 0 when l i ≠ 0, l i+1 ≠ 0, and (5) and (27), it can be seen that when D ≠ 0, It shows that D and L can be estimated by three NVTDOA measurements when UAV must move but not radially. So the target radar position can be determined.

Observability Analysis with the Radar Rotation Period
Estimated. When the period T is unknown, the period estimate X * � D, L, T { } T and the state transition matrix Φ i,i+1 � I 3×3 need to be added to analyze the observability. Without affecting observability, the period integer of two adjacent measurements of UAV is assumed to be 1:

International Journal of Antennas and Propagation
Similarly, the Gramer matrix based on NVTDOA by four measurements is From the calculation in the previous section, we can see the following: In summary, the observability conditions of equation (30) can be obtained; that is, at is to say, D, L, and Tcan be estimated by four NVTDOA measurements when UAV must move but not radially. So the target radar position can be determined.

Construction of NVTDOA Localization Equation.
When the UAV moves, we have e value of Ψ j,i can be estimated by Section 3.2. Referring to equation (13), it can be concluded that where K j,i � round(Ψ j,i /T). e distance from S i to S 0 is d i at t i time, the distance from S j to S 0 is d j , and the distance from S j to S i is l i . It should be noted that from the conclusion of Section 4.2, at low speed the observability condition of the model is related to the change of the position of the UAV and has nothing to do with its speed. If the position of the UAV is recorded during the movement, the error will be caused because the initial time t i is not well determined. In order to locate S i accurately, the UAV hovers approximately before t i to obtain the exact position of the main lobe of the receiving radar antenna during the ith scanning period. After receiving the main lobe pulse in this period, the UAV flies to position S j until receiving the main lobe pulse in the jth scanning period of the radar antenna.
As shown in Figure 1, the geometric relationship diagram of periodic emitter location based on NVTDOA is obtained. For convenience of representation, let S i+1 � S j and S i+2 � S l . In the Cartesian rectangular two-dimensional coordinate system, the radar emitter is located at It is known from the Beidou localization system that its error can be ignored. Based on the accurate estimation of Ψ j,i , equation (33) can be reduced to and Ω � (Π/T). e radar system rotation period angle is Π � 360 ∘ . e approximation method used in the equation is the same as condition (1a) in Section 4. Substituting d i , d j , and l i into θ j,i , we have Similarly, when the UAV flies from S i � [x j , y j ] T to S l � [x l , y l ] T , there will be

International Journal of Antennas and Propagation
By solving the nonlinear equations consisting of equations (36) and (37), the localization results can be obtained that S 0 � [x 0 , y 0 ] T .

Solution of Nonlinear Localization Equations Based on
Cuckoo Search Algorithm 5.2.1. Cuckoo Search Algorithm. e cuckoo search (CS) algorithm is generated by the simulation study of cuckoo's nesting behavior by Yang Xin-she and DEB Suash. e individuals in birds flocking select nests by Levy flight and abandon by discovered probability P a , so the birds population is updated. e paths and locations of cuckoo nests are updated as follows: best is the optimal solution of the tth generation, Levy(λ) ∼ ϕ × u/|v| 1/λ is the dot product, n is the number of nests, and Levy(λ) ∼ ϕ × u/|v| 1/λ is the Levy random search path whose moving step size obeys levy probability distribution as where u and v follow normal distribution and ϕ � ((Γ(1 + λ)sin(λπ/2))/(Γ(1 + λ/2) × 2 (1+λ/2) λ)), where Γ is the standard Gamma function. Combining equations (38) and (39), in CS algorithm, the new solution x t+1 i is generated by After discarding a part of solutions according to the discovery probability, a new solution is generated by where r is the random number of uniform distribution in (0, 1) interval and x t j and x t e are random solutions of the tth generation.

Solution of Nonlinear Localization Equations Based on Cuckoo Search Algorithm.
e localization equations consisting of equations (36) and (37) are recorded as When equation (42) is transformed into the unconstrained optimization problem, the objective function of the cuckoo search algorithm is where w � (x, y) T . So solving equation (42) is to find the minimum solution of equation (43). To sum up, the solution process of CS is as follows: (1) Initialization. e population number of the cuckoo algorithm is m, the maximum number of iterations is K, and the discovered probability is P a . Moreover, m nest positions are randomly generated and the optimal nest position of the current generation,w 0,best (best ∈ 1, 2, . . . , m { }), and the optimal fitness value, F min , are found. In the simulation experiment of the paper, m � 20, K � 2000, and P a � 0.25.
(2) Update of Levy Flight Path (Loop Body). e optimal value of the previous generation, w k− 1,best , is reserved. K represents the current evolutionary generation. And a new nest location is generated by equation (40). en we compare the optimal fitness value of the new nest location with that of the previous nest location and keep the solution with the smaller optimal fitness value which is expressed as I k � [w k,1 , w k,2 , . . . , w k,m ] to the next nest location.

(3) Update of the Discovered Probability (Loop Body).
For each individual in q k � [w k,1 , w k,2 , · · · , w k,m ], r is randomly generated as the probability that the bird's egg is found by the nest owner, that is, the discovered probability. If r < P a , the nest value is unchanged; otherwise, the nest is updated according to equation (41). e new nest is compared with q k � [w k,1 , w k,2 , · · · , w k,m ], and then the better solution which is expressed as q k � [w k,1 , w k,2 , . . . , w k,m ] is kept to the next generation. (4) Judgment Body. e optimal nest position w k,best in q k and the fitness value F min of the solution are calculated. At the same time, it is determined whether the termination condition of the algorithm (the generation is equal to K) is satisfied. If it is satisfied, the global  International Journal of Antennas and Propagation optimal solution is obtained. Otherwise, step 2 is returned, and the iterative operation continues until the termination condition is satisfied.

Simulation Parameter Setting.
e American under AN/SPS49 air-to-air search radar is taken as the simulation object in this paper. e main parameters of the radar are shown in Table 1.
e horizontal pattern of the antenna main lobe beam is of Gaussian type, and the classical mathematical model as equation (44) is used in the simulation: Without loss of generality, the main performance parameters of UAV and receiving loads are set in Table 2.
In this section, the localization accuracy for the radar antenna by the above algorithm is verified. e localization accuracy is calculated by MSE(S 0 ) � L l�1 ‖S l 0 − S 0 ‖ 2 /L, where S l 0 is the estimated value of the lth Monte Carlo simulation of the radar antenna localization, L is the number of simulation, and 1000 is selected here.

Sensitivity of S 0 Estimation to Pulse Repetition Frequency and Received SNR.
e radar pulse energy received by UAV is inversely proportional to the quadratic of d i . According to the radar average power, there will also be a higher SNR. erefore, the radar radiation source is simulated at different positions, and the three target positions are selected as   [20,40] dB. When the radar pulse repetition frequency is PRF � 1000 Hz,PRF � 2000 Hz, the pulse signal time width is τ � 125 μs, the UAV speed is v � 50 m/s, UAV flight time for the single localization ΔT j,i � ΔT j,l � 60 s, and the localization of three radars are S 0,1 � [60, 80] T km, S 0,2 � [90, 90] T km, and S 0,3 � [160, 145] T km. ree radars are positioned separately. at is, the UAV receives only one radar signal at a time, and there is no signal aliasing. In order to verify the accuracy and effectiveness of the proposed algorithm, Hao's method in [21] is compared because Hao's method and the method in this paper both give the solution of the VTDOA localization algorithm. MSE (S 0 ) with different SNRs and PRFs is shown in Figure 4.
As shown in Figure 4, with the increase of SNR, the localization accuracy of Hao and the proposed algorithm also become higher. When the noise power becomes smaller, the influence of Δt i+1,i becomes smaller, which is more in line with the neglect condition of Δt i+1,i . So the localization accuracy is higher. When the PRF is increased, the UAV can obtain more useful information in the same time. So the localization accuracy is higher. As shown in Figure 4, with the increase of PRF, MSE (S 0 ) is reduced and the localization accuracy is higher, which is consistent with the theoretical analysis. By comparing and analyzing the simulation curve of the proposed algorithm and Hao's algorithm, it can be seen that the localization accuracy is improved under different PRFs and SNRs with the proposed algorithm in this paper. However, for the same target radar, PRF � 1000 Hz and PRF � 2000 Hz, the proposed algorithm has little improvement compared with Hao's algorithm. is shows that with the increase of PRF, the improvement of accuracy with the proposed algorithm is not obvious compared with that of Hao's algorithm. is is the result of the estimation of Ψ j,i in this algorithm is the same as that in Hao's. However, with the difference of d i , the improvement of Hao's accuracy will be different which will be analyzed and explained below.
When PRF � 2000 Hz, the curves of localization accuracy of three radars versus SNR are shown in Figure 5. As shown in Figure 5, when PRF and SNR are the same, the MSE (S 0 ) becomes lower as d i increases. Because as d i increases, SNB accepted by UAV increases. But the localization accuracy of this algorithm is higher than that of Hao's algorithm. With different d i the improvement of localization accuracy by the proposed algorithm compared with Hao's algorithm is shown in Figure 6. It can be seen that the improvement becomes better as d i decreases. So the proposed algorithm works better in the relatively close localization area, which is also the requirement of the small aircraft to detect the target at close range.
6.3. Sensitivity of S 0 Estimation to l i . By changing ΔT j,i or v of UAV, the distance l i between two UAVs in a single measurement is changed, where l i � ΔT j,i v. In the simulation, v � 50 m/s , which is kept unchanged, l i is changed by changing ΔT j,i . e accuracy varies with ΔT j,i at PRF � 1000 Hz, SNR � 30 dB, and τ � 125 μs as shown in Figure 7. When ΔT j,i increases from 10 s to 150 s, the localization accuracy of the three radars have the same trend. In Hao's algorithm, when ΔT j,i ∈ [50 s, 100 s], the best localization performance is obtained. When ΔT j,i < 50 s, with the increase of ΔT j,i , the accuracy is greatly improved. Because the aperture of the VTDOA localization system increases in a reasonable range, it will help to improve accuracy. But when ΔT j,i > 100 s, the accuracy becomes worse.
is is due to that the wavefront of radar signal received by UAV twice will not be the same plane wave because of the excessive l i . erefore, the calculation method of rotation time difference of the radar antenna in Hao's  International Journal of Antennas and Propagation algorithm brings more errors. In this paper, the NVTDOA localization equations are proposed by using the improved rotation time difference of the radar antenna calculation method. e Newton iteration is used to obtain the convergence solution, which avoids the approximate calculation in Hao's method. And ΔMSE(S 0 ) is not less than 100.
By the proposed algorithm, with the increase of ΔT j,i , the accuracy increases, and there is no case that the accuracy deteriorates when ΔT j,i is too large. In order to further illustrate the advantages of the proposed algorithm, the improvement of localization accuracy by the algorithm compared with Hao's algorithm is shown in Figure 8. With ΔT j,i increasing, the accuracy increases, which is consistent with the above analysis. e accuracy with smaller D i is better, so the accuracy of close-range target detection by UAV is better improved in this paper.

Conclusion
When the distance between localization points of UAV is large, the localization error of the VTDOA algorithm increases. First, based on the shortcomings of VTDOA, the NVTDOA localization model of periodic radiation source based on single-station small-load UAV is established considering the integral characteristics of radar antenna rotation time difference. Second, the NVTDOA localization equation is constructed and the convergence solution of the NVTDOA localization equation is obtained by the Cuckoo search algorithm to realize the localization of the radar radiation source. Last, the simulation and comparative analysis of the influence of positioning distance, pulse repetition frequency, SNR, and UAV linear distance on location precision are carried out. e simulation results show that Cuckoo search NVTDOA can improve the positioning accuracy and some conclusions are as follows: (1) Comparing the localization accuracy of three different position radars, we can see that, in the relatively close localization area, the improvement effect is better. So the NVTDOA algorithm has higher practical value in the close-range target detection by UAV.
(2) With different SNRs, the localization accuracy of NVTDOA is better than that of VTDOA. erefore, NVTDOA has stronger noise robustness. (3) e simulation experiments of different motion time of the UAV are carried out. e results show that the smaller the linear distance between two positioning of UAV is, the larger the localization error of the VTDOA algorithm is. But localization accuracy of NVTDOA has little dependence on UAV localization distance, which shows that NVTDOA is more suitable for localization of UAV with high speed and long motion time.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.