Detecting and localization method of unknown signal using aircraft with ADS-B system

Abstract

A method of detecting an unknown signal and estimating a source location of the unknown signal using aircraft based on an automatic dependent surveillance-broadcast (ADS-B) system is provided. The method includes a first step (S100) for obtaining from a plurality of airborne aircrafts provided with a network system, aircraft signals transmitted to an air traffic control (ATC) and a second step (S200) for detecting, by the ATC, a presence of an unknown signal in the aircraft signals based on one of a time difference of arrival (TDOA) method, a time of arrival (TOA) method, and an angle of arrival (AOA) method. The method further includes a third step (S300) for estimating the source location of the unknown signal and a fourth step (S400) for transmitting unknown signal generation information associated with the unknown signal and the source location of the unknown signal to neighboring aircraft and the ATC.

Claims

1. A method of detecting an unknown signal and estimating a source location of the unknown signal using network connected aircraft, the method comprising: obtaining, in a first step (S100), from a plurality of airborne aircrafts operating in a wireless communication network implementing an automatic dependent surveillance-broadcast (ADS-B) system, aircraft signals transmitted to an air traffic control (ATC), the aircraft signals comprising global navigation satellite system (GNSS) signals transmitted from navigation satellites and an occurrence of disturbance or interference of the GNSS signals, wherein in the ADS-B system, the ATC serves as a wireless access point and the plurality of airborne aircrafts serve as a client; detecting, in a second step (S200), by the ATC, a presence of an unknown signal in the aircraft signals based on the aircraft signals transmitted from the plurality of aircraft in the first step (S100); wherein the airborne aircrafts are located within a physical distance capable of receiving the occurrence of disturbance or interference at a predetermined relevant time period, the disturbance or interference encompassed within a frequency spectrum of the GNSS signals; wherein detecting in the second step (S200) further comprises: dividing the aircraft signals into identical time blocks; converting divided aircraft signals to a plurality of frequency areas; obtaining power statistics from the plurality of frequency areas; applying a spectrum mask of the GNSS to the power statistics; determining a presence of a power of a frequency deviated from a power of a frequency spectrum mask of the GNSS signals; and determining the occurrence of disturbance or interference of the GNSS signals based on the presence of the deviated power of the frequency; estimating, in a third step (S300), the source location of the unknown signal based on one of a TDOA method, a TOA method, and an AOA method and an altitude of the unknown signal when the occurrence of disturbance or interference associated with the unknown signal on the GNSS signals is detected from the aircraft signals in the second step (S200); and transmitting, in a fourth step (S400), a message including a source code indicative of the presence of the occurrence of disturbance or the interference associated with the unknown signal, detected in the second step (S200) and the source location of the unknown signal, estimated in the third step (S300), to one or more neighboring aircrafts through flight information services-broadcast (FIS-B) or a notice to airmen (NOTAM); wherein in the detection step (S200), detection of the presence of the unknown signal is not based on a signal received at a receiver installed at the ATC, thereby enhancing the precision of the detection of the presence of the unknown signal.

2. The method of claim 1, wherein the third step (S300) comprises: estimating the source location of the unknown signal using the TDOA measurement .sub.ri when the presence of the unknown signal is detected based on the aircraft signals in the second step (S200), wherein .sub.ri is a difference of delayed arrival times of two of the aircraft signals.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a diagram illustrating an overview configuration of an automatic dependent surveillance-broadcast (ADS-B) system according to an embodiment of the present invention.

(2) FIGS. 2 and 3 are flowcharts illustrating a method of detecting an unknown signal and detecting a source location of the unknown signal using aircraft based on an ADS-B system according to an embodiment of the present invention.

(3) FIG. 4A is a diagram illustrating an overview of a method of estimating a source location of an unknown signal according to an embodiment of the present invention.

(4) FIG. 4B is a graph illustrating a time difference of arrival (TDOA) measurement when a cross-correlation function is applied according to an embodiment of the present invention.

(5) FIG. 5 is a structure diagram illustrating a method of detecting an unknown signal and detecting a source location of the unknown signal using aircraft based on an ADS-B system according to an embodiment of the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(6) Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the like elements throughout. The embodiments are described below in order to explain the present invention by referring to the figures.

(7) An aspect of the present invention provides a method of detecting an unknown signal and estimating a source location of the unknown signal that rapidly and precisely detects an unknown signal and estimates a source location of the unknown signal using airborne aircraft provided with an automatic dependent surveillance-broadcast (ADS-B) system. The method enables detection of an occurrence of the unknown signal and detection of the source location of the unknown signal from any one of an air traffic control (ATC) and the aircraft. Accordingly, hereinafter, the method will be discussed in two embodiments based on an estimation agent.

Embodiment 1: Estimation Performed by Ground ATC

(8) A method of detecting an unknown signal and estimating a source location of the unknown signal performed by an ATC includes obtaining, in a first step (S100), a signal, detecting, in a second step, (S200) an unknown signal, estimating, in a third step (S300), a source location of the unknown signal, and transmitting, in a fourth step (S400) unknown signal generation information.

(9) The first step (S100), obtaining of the signal includes receiving, by a plurality of airborne aircraft provided with an ADS-B system, aircraft signals comprising global navigation satellite system (GNSS) signals transmitted from navigation satellites and an unknown signal transmitted from a source location of the unknown signal, and transmitting the aircraft signals comprising the GNSS signal and the unknown signal to an ATC through the ADS-B system.

(10) When a jammer to noise (J/N) power ratio measurement is used to detect the unknown signal, the plurality of airborne aircraft transmits the J/N measurements to the ATC on land in lieu of the GNSS signal and the unknown signal through the ADS-B system. Descriptions pertaining to the J/N measurements will be provided later.

(11) In a manner dissimilar to a method of detecting an unknown signal using a signal received from conventional ground receivers through an ATC, the method according to the present embodiment receives, by airborne aircraft, signals transmitted from navigation satellites and an unknown source. Here, the signals are transmitted to the ground ATC through the ADS-B system in the form of aircraft signals. As previously described, since a receiver installed on land such as the ATC may not receive an unknown signal, and receive a signal distorted by refraction or reflection, the method according to the present disclosure is provided to enhance precision of detecting the unknown signal using the signal received from the aircraft rather than the signal received from the ATC.

(12) According to an aspect of the present invention, the unknown signal may be continuously detected by the aircraft provided with the ADS-B system that receives the GNSS signals transmitted from the navigation satellites and the unknown signal transmitted from the unknown source location, and, through aircraft signals, transmits the signals to the ground ATC through the ADS-B system. To this end, according to an aspect of the present invention, a message, for example, a source code, to detect a presence of the unknown signal is included in addition to information on a location, a speed, and an altitude of an aircraft basically collected from the ADS-B system, thus establishing a system for detecting an unknown signal without changing or additionally installing a communication instrument.

(13) The second step (S200), detecting of the unknown signal includes analyzing aircraft signals transmitted from a plurality of aircraft, and detecting a presence of an unknown signal from the aircraft signals.

(14) In the analysis of the signals transmitted from the plurality of aircraft, when the unknown signal is included in the aircraft signals, each of the plurality of aircraft divides the aircraft signals transmitted from the plurality of aircraft into identical time blocks, converts the signals to a plurality of frequency areas, and obtains power statistics from the plurality of frequency areas. When a spectrum mask of a legitimate signal is applied to the power statistics, and a frequency having power significantly deviating from a frequency mask is verified to be present, interference by an unknown signal is determined to be detected, and therefore a presence of the unknown signal is detected.

(15) The third step (S300), estimating of the source location of the unknown signal includes estimating a source location of an unknown signal when the unknown signal is determined to be included in signals transmitted from a plurality of aircraft through the foregoing process.

(16) When the unknown signal is wirelessly propagated from the source location, the unknown signal is emitted to all directions as a spherical wave. In this example, a time delay occurs because the unknown signal, for example, a wireless signal, moves at a predetermined speed, and power decreases due to spherical spreading. In detail, when the unknown signal is emitted from the unknown source, a time delay occurs in proportion to a distance between a transmission point and a reception point, and power is reduced by an amount proportional to the distance between the transmission point and the reception point. As a result, the source location of the unknown signal is estimated based on such characteristics of the time delay and the power reduction.

(17) Accordingly, the method disclosed herein may be classified into a method based on time delay and a method based on power change during propagation of an unknown signal. Examples of the method based on the time delay may include a time of arrival (TOA) method based on an arrival time of a propagated signal, a time difference of arrival (TDOA) method based on a difference of arrival times of propagated signals, and an angle of arrival (AOA) method based on an arrival angle of a signal.

(18) Also, examples of the method based on the power change may include a received signal strength (RSS) method based on strength of a received signal, a power of arrival (POA) method, a power difference of arrival (PDOA) method, and a frequency difference of arrival (FDOA) method based on a frequency change due to Doppler effects caused by a movement of a receiver.

(19) However, the TOA method generally used in a GNSS requires time synchronization between a transmission point in time and a reception point in time. The AOA method, albeit insusceptible to the time synchronization, requires antenna alignment amongst reception sensors and has inferior location estimation precision to the TDOA method. The TDOA method exhibits greater efficiency as a method of estimating a location of an unknown signal than the other methods because the TDOA method may be applied to an instance in which an input signal is unknown.

(20) Hereinafter, the method of estimating the source location of the unknown signal based on the TDOA method will be described.

(21) A TDOA-based location estimation method employs a difference of delayed arrival times of signals includes disposing a reference sensor and a predetermined sensor differing from the reference sensor, and obtaining a location of a signal source using a difference of delayed arrival times of signals received from each of the sensors. In this example, the location of the signal source obtained by the difference of delayed arrival times is represented by a plurality of hyperbolic curves as shown in FIG. 4A. The location of the signal source is estimated by obtaining a point at which the plurality of hyperbolic curves converges.

(22) The difference of delayed arrival times between the two sensors is obtained using a cross-correlation function expressed by Equation 1. When the signals received from the two sensors s.sub.r and s.sub.i, respectively, are cross-correlated, a cross-correlation value R.sub.ri() indicated by a curve of FIG. 4A has a single maximum value. In this example, a delay time until reaching the maximum value is a TDOA measurement, for example, a difference of delayed arrival times denoted by .sub.ri.

(23) $\begin{array}{cc}{R}_{\mathrm{ri}}\left(\right)E\left[s_r\left(t\right)s_i\left(t-\right)\right]=\frac{1}{T}{}_0^T{s}_r\left(t\right)s_i\left(t-\right)d& \left[\mathrm{Equation}1\right]\end{array}$

(24) In Equation 1, R.sub.ri() denotes a cross-correlation value of signals received from the reference sensor s.sup.r and an i-th sensor s.sub.i. s.sub.r(t) and s.sub.i(t) each denotes a signal received from the reference sensor s.sup.r and a signal received from the i-th sensor s.sub.i at a t point in time. T denotes an integral time, and denotes a delay time.

(25) A difference of the delayed arrival times between the reference sensor s.sub.r provided on an aircraft and the i-th sensor s.sub.i is a TDOA measurement .sup.ri, which corresponds to a difference of locations with respect to a speed c of signals. Accordingly, the TDOA measurement .sub.ri is calculated using a function with respect to the source location of the unknown signal and the locations of the two sensors expressed by Equation 4.

(26) $\begin{array}{cc}\begin{array}{c}{}_{\mathrm{ri}}=\left(t_r-t_s\right)-\left(t_i-t_s\right)\\ =\frac{\sqrt{{\left(x-x_r\right)}^2+{\left(y-y_r\right)}^2}-\sqrt{{\left(x-x_r\right)}^2+{\left(y-y_i\right)}^2}}{c}\\ =f_{r,i}\left(x,y,x_r,y_r,x_i,y_i\right)\end{array}& \left[\mathrm{Equation}4\right]\end{array}$

(27) In Equation 4, t.sub.s denotes a point in time at which a signal is transmitted from an unknown source, t.sub.r denotes a point in time at which the reference sensor s.sub.r receives the unknown signal, t.sub.i denotes a point in time at which the i-th sensor receives the unknown signal, (x,y) denotes location coordinates of the unknown signal source, (x.sub.r,y.sub.r) denotes location coordinates of the reference sensor s.sub.r, and c denotes a speed of a signal.

(28) When the TDOA measurement .sub.ri is obtained, the source location (x,y) of the unknown signal is calculated by Equation 4. (x,y) is not calculated directly because Equation 4 is a non-linear equation. Accordingly, Equation 4 is linearized to be represented by Equation 5 by applying Taylor series to Equation 4.

(29) $\begin{array}{cc}\begin{array}{c}{}_{\mathrm{ri}}=f_{r,i}(.Math.)\\ {f_{r,i}(.Math.).Math.}_{\left(x_0,y_0\right)}+\frac{{f_{r,i}(.Math.).Math.}_{\left(x_0,y_0\right)}}{x_0}x+\frac{{f_{r,i}(.Math.).Math.}_{\left(x_0,y_0\right)}}{y_0}y\\ =\frac{1}{c}\sqrt{{\left(x_0-x_r\right)}^2+{\left(y_0-y_r\right)}^2}-\\ \sqrt{{\left(x_0-x_i\right)}^2+{\left(y_0-y_i\right)}^2}+\frac{1}{c}\frac{x_0-x_i}{\sqrt{{\left(x_0-x_i\right)}^2+{\left(y_0-y_i\right)}^2}}-\\ \frac{x_0-x_1}{\sqrt{{\left(x_0-x_1\right)}^2+{\left(y_0-y_1\right)}^2}}+\\ \frac{1}{c}\frac{y_0-y_i}{\sqrt{{\left(x_0-x_i\right)}^2+{\left(y_0-y_i\right)}^2}}-\frac{y_0-y_1}{\sqrt{{\left(x_0-x_1\right)}^2+{\left(y_0-y_1\right)}^2}}\end{array}& \left[\mathrm{Equation}5\right]\end{array}$

(30) In Equation 5, (x.sub.0,y.sub.0) denotes initial location coordinates of the unknown signal source.

(31) When N number of sensors exist, Equation 6 in a matrix type is given by setting a first sensor to be a reference sensor r, and applying each of remaining sensors, for example, a second sensor through an N-th sensor, to an i-th sensor one at a time in a sequence.

(32) $\begin{array}{cc}\frac{1}{c}\left[\begin{array}{cc}{G}_{x21}& G_{y21}\\ G_{x31}& G_{y31}\\ \mathrm{.Math.}& \mathrm{.Math.}\\ G_{\mathrm{xN}1}& G_{\mathrm{yN}1}\end{array}\right]\left[\begin{array}{c}x\\ y\end{array}\right]=\left[\begin{array}{c}{}_{12}-{}_{12}{.Math.}_{\left(x_0,y_0\right)}\\ {{}_{13}-{}_{13}.Math.}_{\left(x_0,y_0\right)}\\ \mathrm{.Math.}\\ {{}_{1N}-{}_{1N}.Math.}_{\left(x_0,y_0\right)}\end{array}\right]G=Z& \left[\mathrm{Equation}6\right]\end{array}$

(33) In Equation 6, G.sub.xi1 is calculated by

(34) $\frac{x_0-\mathrm{xi}}{\sqrt{{\left(x_0-x_i\right)}^2+{\left(y_0-y_i\right)}^2}}-\frac{x_0-x_1}{\sqrt{{\left(x_0-x_1\right)}^2+{\left(y_0-y_1\right)}^2}},$ G.sub.yi1 is calculated by

(35) $\frac{y_0-y_i}{\sqrt{{\left(x_0-x_i\right)}^2+{\left(y_0-y_i\right)}^2}}-\frac{y_0-y_1}{\sqrt{{\left(x_0-x_1\right)}^2+{\left(y_0-y_1\right)}^2}}$ .sub.ri|(x.sub.0,y.sub.0) is calculated by ({square root over ((x.sub.0x.sub.r).sup.2+(y.sub.0y.sub.r).sup.2)}{square root over ((x.sub.0x.sub.i).sup.2+(y.sub.0y.sub.i).sup.2)})/c, denotes a location variation, and denotes a residual, for example, a difference between a measurement and an estimation.

(36) An estimation of a location variation is used because the location variation is difficult to be calculated directly by Equation 6, which is represented by Equation 7 suggested by W. H. Foy in Position-location Solution by Taylor-series Estimation (The Institute of Electrical and Electronics Engineers (IEEE) Transactions on Aerospace and Electronic System, vol. AES-12, no. 2, pp. 187-194, March 1976), which is hereby incorporated by reference herein.
{circumflex over ()}=[G.sup.TQ.sup.1G].sup.1G.sup.TQ.sup.1Z[Equation 7]

(37) In Equation 7, Q denotes a covariance matrix of a measured error.

(38) Accordingly, as represented by Equation 8, estimated location coordinates ({circumflex over (x)},) of the unknown signal source is calculated by a sum of an initial location (x.sub.0,y.sub.0) and the location variation {circumflex over ()} estimated by Equation 7.

(39) $\begin{array}{cc}\left[\begin{array}{c}\hat{x}\\ \hat{y}\end{array}\right]=\left[\begin{array}{c}{x}_0\\ y_0\end{array}\right]=\left[\begin{array}{c}\hat{x}\\ \hat{y}\end{array}\right]& \left[\mathrm{Equation}8\right]\end{array}$

(40) As represented by Equation 8, the estimated location coordinates ({circumflex over (x)},) are largely based on the initial location (x.sub.0,y.sub.0). As a result, a final estimated location coordinates ({circumflex over (x)},) of the unknown signal source are not directly obtained through a single calculation, however, rather obtained by re-applying the estimated location coordinates ({circumflex over (x)},) to the initial location (x.sub.0,y.sub.0) and iterating the foregoing process of Equations 5 through 8. The iterating of the process is performed until a stop condition is satisfied. An example of the stop condition may include an instance in which the location variation is {circumflex over ()} less than a predetermined threshold value (TH).

(41) Thus far, the method of estimating the unknown signal source location based on the TDOA method is provided. However, the estimation is also possible by the other methods, for example, the RSS method or a differential RSS (DRSS) method based on a J/N measurement.

(42) The RSS method using the J/N measurement is a method of estimating a source location of an unknown signal based on a difference of strength of received signals using a J/N measurement. As used herein, the J/N measurement refers to a strength of an unknown signal with respect to noise.

(43) When an unknown signal is included in a signal received by a GNSS receiver, strength of a J/N measurement varies, which enables an easy detection of the unknown signal and estimation of a source location of the unknown signal. The J/N measurement is obtained by a J/N meter provided along with an analog-to-digital (A/D) converter that configures the GNSS receiver.

(44) A plurality of J/N measurements obtained through the foregoing process is transmitted and collected by a neighboring aircraft or a ground ATC. When the collected J/N measurement is greater than a reference value set using a range of J/N measurements in normal circumstances, the unknown signal is determined to be present, and an occurrence of an error in a predetermined receiver is also determined based on the plurality of collected J/N measurements.

(45) Hereinafter, descriptions pertaining to estimation of the source location of the unknown signal based on a variation in the J/N measurement in response to a distance change will be provided.

(46) A difference dJ.sub.ij of strength of unknown signals received by receivers (i, j), respectively, is represented by Equation 9 based on a path loss model.

(47) $\begin{array}{cc}{\mathrm{dJ}}_{\mathrm{ij}}=J_i-J_j=10{\log }_{10}\left(\frac{d_j}{d_i}\right)& \left[\mathrm{Equation}9\right]\end{array}$

(48) In Equation 9, J.sub.i and J.sub.j each denote strengths of unknown signals received by receivers (i,j), denotes a path loss index, d.sub.i and d.sub.j each denote a distance from a source location of an unknown signal to the receivers (i,j).

(49) As represented by Equation 9, the path loss index varies based on a propagation environment, and in an absence of an obstacle between the source location of the unknown signal and the receivers, the path loss index is 2. An approximate value close to 2 is applied to the path loss index according to the method of the present invention because such obstacles are rarely present between a source location of an unknown signal and a GNSS receiver provided in an aircraft.

(50) The difference dJ.sub.ij of strength of unknown signals is calculated by Equation 9 as described above, however, also calculated by Equation 10 using a difference of J/N measurements obtained by the receivers (i,j). In Equation 10, JNR.sub.i and JNR.sub.j each denote J/N measurement in the receivers (i,j). N.sub.i denotes noise of the receiver (i), and since the noise may differ in a plurality of receivers based on design methodologies, a standardization or revision process may be required as necessary.
dJ.sub.ij=JNR.sub.iJNR.sub.j=(J.sub.iN.sub.i)(J.sub.jN.sub.j)J.sub.iJ.sub.j[Equation 10]

(51) When the difference dJ.sub.ij of strength of the unknown signals is calculated by Equation 10, the source location

(52) $X=\left[\begin{array}{c}x\\ y\end{array}\right]$
of the unknown signal is calculated by Equation 11.
AX=B[Equation 11]

(53) In Equation 11,

(54) 0$A=\left[\begin{array}{cc}{x}_{c13}-x_{c12}& y_{c13}-y_{c12}\\ x_{c23}-x_{c12}& y_{c23}-y_{c12}\\ x_{c23}-x_{c13}& y_{c23}-y_{c13}\end{array}\right],B=\frac{1}{2}\left[\begin{array}{c}{R}_{12}^2-R_{13}^2-\left(x_{c12}^2+y_{c12}^2\right)+\left(x_{c13}^2+y_{c13}^2\right)\\ R_{12}^2-R_{23}^2-\left(x_{c12}^2+y_{c12}^2\right)+\left(x_{c23}^2+y_{c23}^2\right)\\ R_{13}^2-R_{23}^2-\left(x_{c13}^2+y_{c13}^2\right)+\left(x_{c23}^2+y_{c23}^2\right)\end{array}\right],x_{\mathrm{cij}}=\left(r_{\mathrm{ij}}^2{x}_j-x_i\right)/\left(r_{\mathrm{ij}}^2-1\right),y_{\mathrm{cij}}=\left(r_{\mathrm{ij}}^2{y}_j-y_i\right)/\left(r_{\mathrm{ij}}^2-1\right),R_{\mathrm{ij}}^2=\left[1/{\left(r_{\mathrm{ij}}^2-1\right)}^2+1/\left(r_{\mathrm{ij}}^2-1\right)\right]\left[{\left(x_i-x_j\right)}^2+{\left(y_i-y_j\right)}^2\right],r_{\mathrm{ij}}=\frac{d_i}{d_j}={10}^{\frac{{\mathrm{dJ}}_{\mathrm{ji}}}{10}},$
and (x.sub.i,y.sub.i) denotes a location of the receiver (i).

(55) Accordingly, estimation of the source location of the unknown signal using the J/N measurements is relatively simple, inexpensive, and convenient when compared to the TDOA method because transmission of massive data is unnecessary.

(56) The fourth step (S400), transmitting S400 of the unknown signal generation information includes, in the case of detecting, by a ground ATC, an occurrence of an unknown signal and thereby estimating a source location of the unknown signal, inputting, by a ground ADS-B system, the information to a computer in a control console of a neighboring aircraft through flight information services-broadcast (FIS-B) or a notice to airmen (NOTAM), informing the aircraft flying in a neighborhood of the occurrence of the unknown signal, and at the same time, converting to a backup system, for example, a conventional radar facility or a conventional navigation system. In this example, flight safety may be significantly increased by a ground system transmitting a location solution to the neighboring aircraft using backup navigation equipment, and continuing the flight in a state in which the unknown signal is removed.

Embodiment 2: Estimation Performed by Aircraft

(57) A method of detecting an unknown signal and estimating a source location of the unknown signal performed by aircraft includes obtaining, in a first step (S100), signals, detecting, in a second step (S200), an unknown signal, estimating, in a third step (S300), a source location of the unknown signal, and transmitting, in a fourth step (S400), unknown signal generation information as shown in FIG. 3. Hereinafter, brief descriptions of each operation will be provided since an overall process of Embodiment 2 is the same as Embodiment 1, barring a difference in that the estimation of Embodiment 2 is performed by the aircraft rather than the ATC.

(58) The first step (S100), obtaining of the signal includes receiving, from aircraft provided with an ADS-B system, an aircraft signal comprising GNSS signals transmitted from navigation satellites and an unknown signal transmitted from a source location of the unknown signal, and transmitting the aircraft signal comprising the GNSS signal and the unknown signal to neighboring aircraft through the ADS-B system.

(59) The second step (S200), detecting of the unknown signal includes analyzing aircraft signals transmitted from neighboring aircraft, and detecting a presence of an unknown signal in the aircraft signals.

(60) The third step (S300), estimating of the source location of the unknown signal includes estimating a source location of an unknown signal using the TDOA measurement .sub.ri obtained in the second step (S200) of the unknown signal. In a manner dissimilar to Embodiment 1, the source location of the unknown signal is estimated by an airborne aircraft, rather than a ground ATC. However, description previously provided in Embodiment 1 may be applied to a remaining process aside from the estimation.

(61) The fourth step (S400), transmitting of the unknown signal generation information includes, in the case of detecting, by an aircraft, an occurrence of an unknown signal and thereby estimating a source location of the unknown signal, inputting, by ADS-B equipment provided in the aircraft, the information to a computer in a control console of aircraft flying in a neighborhood through FIS-B or an NOTAM, informing the neighboring airborne aircraft and an ATC of the occurrence of the unknown signal, and at the same time, converting to a backup system.

(62) Although the forgoing descriptions provide that a source location of an unknown signal is detected using aircraft signals obtained from a plurality of airborne aircraft, according to the method of the present invention, a single or a couple of aircraft flying circuitous routes in a predetermined region may be used to collect J/N measurements using signals transmitted from an unknown signal source andnavigation satellites, and detect the source location of the unknown signal using the collected J/N measurements.

(63) According to embodiments of the present invention, it is possible to rapidly and precisely detect an unknown signal and estimate a source location of the unknown signal by collecting information associated with the unknown signal from aircraft provided with an ADS-B system, detecting the source location of the unknown signal using the information, and thus removing effects of refraction or reflection of a signal.

(64) Although a few embodiments of the present invention have been shown and described, the present invention is not limited to the described embodiments. Instead, it would be appreciated by those skilled in the art that changes may be made to these embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.