Airborne direction finding method based on Doppler-phase measurement

Abstract

A new direction finding (DF) method, in which the high-accuracy measuring can be realized only with single baseline, is presented used for airborne based on Doppler-phase measurement. The analysis discovers that the integer of wavelength in radial distance can be directly derived compositely, making use of the velocity vector equation and Doppler shift, as well as Doppler changing rate equation. From this, the integer difference of wavelength in path length difference of radial distance between two adjacent antenna elements can be obtained. As soon as the value less than a wavelength in path length difference is determined by phase difference measurement, the direction angle of target can be obtained. As compared with now existing interferometry first determining phase difference, this sort of direction finding method combining Doppler with phase difference first by determining path length difference does not have phase ambiguity nor require restricting base length. By simple mathematical identity transformation, we can prove that the equation derived in this paper is equivalent to an existing one from phase interferometry. The new method presented in this paper will certainly increase new developing force for the research and development of airborne single station direction finding system.

Key words： phase interferometer; direction finding; Doppler frequency, airborne single station location

Introduction

The phase interferometry is a direction finding (DF) method with better measurement accuracy. It is widely used for active and passive detection system [1-3]. However, for single baseline phase interferometry, there is a contradiction between the accuracy of direction finding and maximum unambiguity angle [4,5]. To solve this problem, the existing method is to utilize multibaseline system including the method combining long baselines with short ones [6-8] and algorithm resolving phase ambiguity with multibaseline [9,10].

In actual application, the method combining long baselines with short ones have two limitations [11,12]. In fact, corresponding baselines will also become extremely small since wavelength is very short for high-frequency signal. At this moment, not only must the antenna element be made very small, but also very high demand is put forward for antenna arrangement. It will bring about coupled between antennas and bring down antenna gain. At the same time, higher demand will be required for measurement accuracy of interferometer. For algorithm resolving phase ambiguity with multibaseline, the computing amount is heavy due to demanding multidimensional integer search [13,14].

The study shows that airborne single baseline interferometer will realize high-accuracy direction finding without phase ambiguity after combining with Doppler information. At the same time, the baseline length can be arbitrarily selected only from the standpoint of measurement principle. This paper presents the basic measurement principle and detailed analytic derivation method of Doppler-phase interference measurement.

Single baseline Doppler-phase difference DF principle

Existing DF expression

For contrasting, we first give the existing expression. Provided that incident wave is approximatively plane wave, the phase difference between two adjacent antennas takes the form

(1)\mathrm{2}\mathrm{\pi }{N}_{\mathrm{0}}+\varphi ={\displaystyle \frac{\mathrm{2}\mathrm{\pi }L}{\lambda }}\sin \theta ,

where N_{\mathrm{0}} is an integer, \lambda is the wavelength, L is the distance between two antennas, \varphi is the phase difference, and \theta is the incidence angle of target signal.

Specific value of Doppler change ratio

Figure 1 shows the geometric relationship of single baseline direction-finding antenna array applying for flying vehicle based on Doppler shift-phase difference measurement. Provided that target is static or low speed, when flying vehicle is in uniform motion, the expressions of Doppler change ratio in two antenna elements are

(2)\stackrel{·}{f_{\mathrm{d}\mathrm{1}}}={\displaystyle \frac{v_{\mathrm{t}\mathrm{1}}^{\mathrm{2}}}{\lambda r_{\mathrm{1}}}},

(3)\stackrel{·}{f_{\mathrm{d}\mathrm{2}}}={\displaystyle \frac{v_{\mathrm{t}\mathrm{2}}^{\mathrm{2}}}{\lambda r_{\mathrm{2}}}},

where r_i(i=\mathrm{1},\mathrm{2}) is the radial distance, and v_{\mathrm{t}i} is the tangential velocity.

The specific value of Eqs. (2) and (3) is

(4)q={\displaystyle \frac{\stackrel{·}{f_{\mathrm{d}\mathrm{2}}}}{\stackrel{·}{f_{\mathrm{d}\mathrm{1}}}}}={\displaystyle \frac{r_{\mathrm{1}}}{r_{\mathrm{2}}}}{\displaystyle \frac{v_{\mathrm{t}\mathrm{2}}^{\mathrm{2}}}{v_{\mathrm{t}\mathrm{1}}^{\mathrm{2}}}}.

According to sine theorem, the specific value of radial distance corresponding with two antennas can be also obtained:

(5){\displaystyle \frac{r_{\mathrm{2}}}{r_{\mathrm{1}}}}={\displaystyle \frac{\sin {\beta }_{\mathrm{1}}}{\sin {\beta }_{\mathrm{2}}}}={\displaystyle \frac{v\sin {\beta }_{\mathrm{1}}}{v\sin {\beta }_{\mathrm{2}}}}={\displaystyle \frac{v_{\mathrm{t}\mathrm{1}}}{v_{\mathrm{t}\mathrm{2}}}}.

It shows that the specific value of radial distance equals to the one of tangential velocity corresponding with two antennas when flying vehicle is uniform motion. Substituting Eq. (5) into Eq. (4) gives

(6)q={\displaystyle \frac{v_{\mathrm{t}\mathrm{2}}^{\mathrm{3}}}{v_{\mathrm{t}\mathrm{1}}^{\mathrm{3}}}}.

Fig.1 Geometric schematic for airborne DF with single baseline

Full size|PPT slide

Integer value of radial distance

According to velocity resolution, the identical relation of velocity can be written as

(7)v^{\mathrm{2}}=v_{\mathrm{r}\mathrm{1}}^{\mathrm{2}}+v_{\mathrm{t}\mathrm{1}}^{\mathrm{2}}=v_{\mathrm{r}\mathrm{2}}^{\mathrm{2}}+v_{\mathrm{t}\mathrm{2}}^{\mathrm{2}}.

where v_{\mathrm{r}i} is the radial velocity.

After rearrangement, we have

(8)v_{\mathrm{r}\mathrm{1}}^{\mathrm{2}}-v_{\mathrm{r}\mathrm{2}}^{\mathrm{2}}=v_{\mathrm{t}\mathrm{2}}^{\mathrm{2}}-v_{\mathrm{t}\mathrm{1}}^{\mathrm{2}}.

Respectively substituting Doppler shift equation and Doppler change ratio as well as specific value into Eq. (8) gives

(9)\lambda (f_{\mathrm{d}\mathrm{1}}^{\mathrm{2}}-f_{\mathrm{d}\mathrm{2}}^{\mathrm{2}})=r_{\mathrm{1}}\stackrel{·}{f_{\mathrm{d}\mathrm{1}}}(u-\mathrm{1}),

(10)\lambda (f_{\mathrm{d}\mathrm{1}}^{\mathrm{2}}-f_{\mathrm{d}\mathrm{2}}^{\mathrm{2}})=r_{\mathrm{2}}\stackrel{·}{f_{\mathrm{d}\mathrm{2}}}\left(\mathrm{1}-{\displaystyle \frac{\mathrm{1}}{u}}\right),

where u=\sqrt[\mathrm{3}]{q^{\mathrm{2}}}.

Thus, we can obtain integer values at about two radial distances:

(11)N_{\mathrm{1}}=\mathrm{int}\left[{\displaystyle \frac{r_{\mathrm{1}}}{\lambda }}\right]=\mathrm{int}\left[{\displaystyle \frac{f_{\mathrm{d}\mathrm{1}}^{\mathrm{2}}-f_{\mathrm{d}\mathrm{2}}^{\mathrm{2}}}{\stackrel{·}{f_{\mathrm{d}\mathrm{1}}}(u-\mathrm{1})}}\right],

(12)N_{\mathrm{2}}=\mathrm{int}\left[{\displaystyle \frac{r_{\mathrm{2}}}{\lambda }}\right]=\mathrm{int}\left[{\displaystyle \frac{(f_{\mathrm{d}\mathrm{1}}^{\mathrm{2}}-f_{\mathrm{d}\mathrm{2}}^{\mathrm{2}})u}{\stackrel{·}{f_{\mathrm{d}\mathrm{2}}}(u-\mathrm{1})}}\right].

DF equation based on Doppler shift-phase difference

Provided that the phase measured by descriminator is respectively {\varphi }_{\mathrm{1}} and {\varphi }_{\mathrm{2}} responding with two radial distances, radial distance can be represented by

(13)r_{\mathrm{1}}=\lambda \left(N_{\mathrm{1}}+{\displaystyle \frac{{\varphi }_{\mathrm{1}}}{\mathrm{2}\mathrm{\pi }}}\right),

(14)r_{\mathrm{2}}=\lambda \left(N_{\mathrm{2}}+{\displaystyle \frac{{\varphi }_{\mathrm{2}}}{\mathrm{2}\mathrm{\pi }}}\right).

Corresponding path difference is

(15)\mathrm{\Delta }r=\lambda \left(\mathrm{\Delta }N+{\displaystyle \frac{\mathrm{\Delta }\varphi }{\mathrm{2}\mathrm{\pi }}}\right),

where \mathrm{\Delta }N=N_{\mathrm{1}}-N_{\mathrm{2}} is the integer difference that has been determined by Doppler shift as well as changing rate, and \mathrm{\Delta }\varphi ={\varphi }_{\mathrm{1}}-{\varphi }_{\mathrm{2}} is the phase difference between two antennas that can be determined by phase difference measurement.

Due to sine theorem, target azimuth is

(16)\sin \theta ={\displaystyle \frac{\mathrm{\Delta }r}{L}}={\displaystyle \frac{\lambda }{L}}\left(\mathrm{\Delta }N+{\displaystyle \frac{\mathrm{\Delta }\varphi }{\mathrm{2}\mathrm{\pi }}}\right).

If N_{\mathrm{0}}=\mathrm{\Delta }N and \varphi =\mathrm{\Delta }\varphi, then we can obtain the entire same result with existing phase interferometry Eq. (1). The different is that there is no phase ambiguity in Doppler shift-phase difference equation, wherein the measurement is perfectly independent of baseline length.

Conclusion

From the point of view of theory analysis, derived equation is only suitable to direction finding when flying vehicle is uniform motion because quoted Doppler changing ratio is only in uniform motion. However, this limitation may be solved by improving arithmetic step by step. At the same time, the uniform flight will be normal for practical aerial reconnaissance.

Since the integer value of radial distance can be first solved by use of Doppler information, there is no phase ambiguity in measurement and restriction for select of baseline based on the method presented in this paper. Hence, there is no contradiction for single baseline phase interferometry between the accuracy of direction finding and maximum unambiguity angle. This means that sense finding of high accuracy can be realized only by single baseline. This method has following advantages: 1) mounting demand is low; 2) detecting time is short; and 3) it is suitable for carrying out detection operation in a broadband.

References

1 Zhao G Q. Radar Countermeasures Principle. Xi’an: Xidian University Press, 1999, 54–56 (in Chinese)

2 Messer H, Singal G. On the achievable DF accuracy of two kinds of active interferometers. IEEE Transactions on Aerospace and Electronic Systems, 1996, 32(3): 1158–1164

3 Ji X G, Gao X G. An airborne passive location algorithm-interferometer locating. Fire Contro & Command Control, 2008, 33(11): 158–161 (in Chinese)

4 Li Y, Zhao G W, Li T. Algorithm of solving interferometer phase difference ambiguity by airborne single position. Chinese Journal of Sensors and Actuators, 2006, 19(6): 2600–2602 (in Chinese)

5 Jacobs E, Ralston E W. Ambiguity resolution in interferometry. IEEE Transactions on Aerospace and Electronic Systems, 1981, AES-17(6): 766–780

6 An X J. Research on DF based on improved phase interferometer. Radio Engineering of China, 2009, 39(3): 59–61 (in Chinese)

7 Wei X, Wan J W, Huang F K. Study of passive location system based on multi base-line interferometers. Modern Radar, 2007, 29(5): 22–25, 35 (in Chinese)

8 Wang Z R. On high precision DOA estimation. Radio Engineering of China, 2007, 37(11): 24–25 (in Chinese)

9 Lin Y M, Liu Y, Zhang Y N. Algorithm of direction finding for broadband digital signal. Journal of Nanjing University of Aeronautics & Astronautics, 2005, 37(3): 335–340 (in Chinese)

10 Zhou Y Q, Huang F K. Solving ambiguity problem of digitized multi-baseline interferometer under noisy circumstance. Journal on Communications, 2005, 26(8): 16–21 (in Chinese)

11 Gong X Y, Huang F K, Yuan J Q. A new algorithm for estimation of direction of arrival based on the second-order difference of phase of interferometer array. Acta Electronica Sinica, 2005, 33(3): 444–446 (in Chinese)

12 Sundaram K R, Murthy U M S, Mallik R J. Module conversion method for estimating the direction of arrival. IEEE Transactions on Aerospace and Electronic Systems, 2000, 36(4): 1391–1396

13 Gong X Y, Yuan J Q, Su L H. A multi-pare unwrap ambiguity of interferometer array for estimation of direction of arrival. Journal of Electronics & Information Technology, 2006, 28(1): 55–59 (in Chinese)

14 Zhang G B, Liu Y, Liu Z M. Unwrapping phase ambiguity algorithm based on baseline ratio. Journal of Nanjing University of Aeronautics & Astronautics, 2008, 40(5): 665–669 (in Chinese)