US20090091490A1 - Method and system for radar tracking of moving target from moving station - Google Patents

Abstract

A method and system is proposed for use by a moving station (such as a jetfighter) for radar tracking of a moving target (such as an air-to-air missile). The proposed method and system involves the use of a hybrid FSK/LFM (Frequency Shift Keying & Linear Frequency Modulation) scheme for acquiring a collection of raw radar data, a first Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the raw radar data, a trilateration module, and a second Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the trilateration-resulted radar data. These features allow the radar tracking of moving objects to be more fast and accurate.

Description

BACKGROUND OF THE INVENTION
• 1. Field of the Invention
• This invention relates to radar technology, and more particularly, to a method and system for radar tracking of a moving target (such as an air-to-air missile) from a moving station (such as a jetfighter) with improved accuracy and performance.
• 2. Description of Related Art
• A radar system is a radio-based surveilance system that uses electromagnetic waves to identify the range, altitude, direction, or speed of moving or fixed objects such as aircraft, ships, and motor vehicles. In operation, a radar system is capabale of emitting a beam of radio wave toward the area under surveilance, and if a target exists in that area, the emitted radio wave will be reflected back. The reflected radio wave is then received by the radar system and analyzed to determine the position, speed, and acceleration of the target.
• Radar systems are categorized into ground-based type and airborne type. The airborne-type radar systems are installed on airplanes such as jetfighters for detecting air-to-air missiles or nearby airplanes. In this type of application, since the radar system and the target are both moving, the detected radar data represents the relative position, relative speed, and relative acceleration of the target.
• Many research efforts have been conducted on airborne radar systems intended for air-to-air missile detections. For example, the technical paper “WAVEFORM DESIGN PRINCIPLES FOR AUTOMOTIVE RADAR SYSTEMS” authored by H. Rohling and M. M. Meinecke, IEEE Radar, pp. 1-4. October 2001, proposes a radar system that utilizes a hybrid FSK/LFM (Frequency Shift Keying & Linear Frequency Modulation) technique for detection of relative distance and relative speed between the radar system and the target.
• One drawback to the aforementioned radar system, however, is that the hybrid FSK/LFM sensor is only capable of detecting the range and speed of the target, and incapable of detecting the acceleration.
• Moreover, many conventional radar systems utilize trilateration to obtain the range, velocity, and acceleration of the target from a collection of raw radar data. One drawback to the trilateration process, however, is that it would undesirably degrade the S/N (signal-to-nose ratio) of the original radar data.
SUMMARY OF THE INVENTION
• It is therefore an objective of this invention to provide a new radar method and system that utilizes a new and modified FSK/LFM technique capable of detecting the acceleration of the target in addition to the target's range and speed.
• It is another objective of this invention to provide a new radar method and system that can provide enhanced S/N ratio of the detected data for improved accuracy of target tracking.
• Defined as a method, the invention comprises: (M1) emitting a radar signal beam; (M2) sensing an echo of the radar signal beam from a moving target and producing a raw set of radar data about the moving target, including range, velocity, and acceleration; (M3) performing S/N (signal-to-noise ratio) enhancement of each raw set of radar data; (M4) performing trilateration on the S/N-enhanced sets of radar data in reference to a predefined 3-dimenional coordinate system having a first axis, a second axis, and a third axis to obtain three axis-oriented sets of radar data respectively in reference to the first axis, the second axis, and the third axis of the 3-dimenional coordinate system; (M5) performing S/N enhancement on the trilateration-resulted sets of radar data; and (M6) a radar data presentation module, which is capable of presenting the S/N-enhanced sets of radar data produced by the second Gaussian-noise filter array in a human-cognizable data form.
• The moving object radar tracking system according to the invention comprises: (A) a radar signal emitter; (B) a radar signal sensor array including an array of N radar signal sensors that are preferably each implemented with a hybrid FSK/LFM sensors; (C) a first Gaussian-noise filter array including an array of N Gaussian-noise filters that are preferably implemented with one-stage linear Kalman filters; (D) a trilateration module; (E) a second Gaussian-noise filter array, which includes at least 3 Gaussian-noise filters that are preferably implemented with an array of one-stage linear Kalman filters; and (F) a radar data presentation module.
• The method and system for tracking moving target from moving station according to the invention is characterized by the use of a hybrid FSK/LFM (Frequency Shift Keying & Linear Frequency Modulation) scheme for acquiring a collection of raw radar data, a first Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the raw radar data, a trilateration module, and a second Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the trilateration-resulted radar data. These features allow the radar tracking of moving objects to be more fast and accurate.
BRIEF DESCRIPTION OF DRAWINGS
• The invention can be more fully understood by reading the following detailed description of the preferred embodiments, with reference made to the accompanying drawings, wherein:
• FIG. 1 is a schematic diagram showing the application of the moving target radar tracking system according to the invention;
• FIG. 2 is a schematic diagram showing the internal architecture of the moving target radar tracking system according to the invention;
• FIG. 3A is a graph showing the characteristic plot of frequency variation of an LFM signal using a conventional FSK/LFM technique; and
• FIG. 3B is a graph showing the characteristic plot of frequency variation of a pair of LFM signals using a modified FSK/LFM technique in accordance with the invention.
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS
• The moving target radar tracking system of the invention according to the invention is disclosed in full details by way of preferred embodiments in the following with reference to the accompanying drawings.
Application of the Invention
• FIG. 1 is a schematic diagram showing an application example of the moving target radar tracking system according to the invention (which is here encapsulated in a box indicated by the reference numeral100). As shown, the moving target radar tracking system of theinvention100 is designed for installation on a moving station, such as ajetfighter10, for thejetfighter10 to detect the position and motion of a moving target, such as an air-to-air missile (AAM)20 that is in constant changes of position and motion relative to thejetfighter10. During operation, when the AAM20 is in proximity to thejetfighter10, the moving object radar tracking system of theinvention100 will respond to this condition by indicating the position, speed, and acceleration of theAAM20 to the pilot of the moving target radar tracking system of theinvention100.
Architectrue of the Invention
• As shown inFIG. 2, in architecture, the moving target radar tracking system of theinvention100 comprises: (A) aradar signal emitter110; (B) a radarsignal sensor array120, which includes an array of N radar signal sensors that are preferably each implemented with a hybrid FSK/LFM sensors; (C) a first Gaussian-noise filter array130, which includes N Gaussian-noise filters that are preferably each implemented with a one-stage linear Kalman filter; (D) atrilateration module140; (E) a second Gaussian-noise filter array150, which includes at least 3 Gaussian-noise filters151,152,153, that are preferably implemented with an array of one-stage linear Kalman filters; and (F) a radardata presentation module160. The respective functions of the constituent components of the moving target radar tracking system of theinvention100 are described in details in the following.
Radar Signal Emitter110
• Theradar signal emitter110 is used for emitting a beam of radar signal into the space under surveillance by thejetfighter10. If an AAM20 is in the proximity of thejetfighter10, the radar signal beam will hit the AAM20 and reflect as an echoed radar signal back to thejetfighter10.
RadarSignal Sensor Array120
• The radarsignal sensor array120 includes an array of N radar signal sensors that are preferably implemented with a plurality of hybrid FSK/LFM sensors, each of which is based on a combination of the FSK (Frequency Shift Keying) and the LFM (Linear Frequency Modulation) techniques for extracting a set of position-motion data about theAAM20, including the range, radial velocity, and radial acceleration of theAAM20 based on the echoed radar signal. The number N can be 2, 3, or more. The N hybrid FSK/LFM sensors in the radarsignal sensor array120 are respectively used for the sampling of a raw set of position-motion data, represented by {r₁, v₁, a₁}, {r₂, v₂, a₂}, . . . , and {r_N, v_N, a_N}, respectively. Compared to the conventional FSK/LFM technique proposed by H. Rohling and M. M. Meinecke in the technical paper “WAVEFORM DESIGN PRINCIPLES FOR AUTOMOTIVE RADAR SYSTEMS”, it is an important aspect of the invention that the hybrid FSK/LFM sensors in the radarsignal sensor array120 is capable of additionally acquiring the acceleration of theAAM20, rather than just the range and velocity of the AAM20 by this prior art. Each hybrid FSK/LFM sensor in the radarsignal sensor array120 is designed to compute for the range r, radial velocity v, and radial acceleration a of theAAM20 based on a modified FSK/LFM technique as described below.
• FIG. 3A a graph showing the characteristic plot of the frequency variation of an LFM signal using the conventional FSK/LFM technique (in this graph, B_sweeprepresents bandwidth and T_LFMrepresents dwell time); whereasFIG. 3B is a graph showing the characteristic plot of the frequency variation of a pair of LFM signals utilized by the invention. The two LFM signals respectively have a positive slope of S_Aand S_B. If 2N samples are to be acquired during an internal of τ=T_LFM/2N, then the frequency increment over 2τ is ƒ_inc=B_sweep/N. In one typical application of the invention, for example, the parameters B_sweep, ƒ_A(t), ƒ_B(t), T_LFM, and N are respectively B_sweep=150 MHz, ƒ_B(t)−ƒ_A(t)=300 kHz, T_LFM=256 ms, and N=256. The shift in the frequency of the second LFM signal measured at the temporal point t=2n·τ from the frequency of the first LFM signal measured at the temporal point t=(2n−1)·τ is represented by ƒ_shift, and ƒ_shift=ƒ_B(0)−ƒ_A(0)+S·τ, where S=S_Afor the first LFM signal and S=S_Afor the first LFM signal. Further, the difference between the frequency of the echoed radar signal and the frequency of the originally-emitted radar signal beam is Δƒ_A=ƒ_τA(t)−ƒ_A(t) for the first LFM signal, and Δƒ_B=ƒ_τB(t)−ƒ_B(t) for the second LFM signal. If we define S_A=Δƒ_A·T_LFMand S_B=Δƒ_B·T_LFM, then during the time interval 0<t<T_LFM, it can be deduced that:
• $\begin{array}{c}{S}_A=T_{\mathrm{LFM}}f_{\mathrm{rA}}\left(t\right)-f_A\left(t\right)\\ =T_{\mathrm{LFM}}\left\{\left[f_0+s\left(t-\frac{2{\tau }_A}{c}\right)\right]\left(1-\frac{2v_A}{c}\right)-\left(f_0+\mathrm{st}\right)\right\}\\ =T_{\mathrm{LFM}}\left[-\frac{2v_A}{c}f_A\left(t\right)-\frac{2{\mathrm{sr}}_A}{c}\left(1-\frac{2v_A}{c}\right)\right]\\ \simeq T_{\mathrm{LFM}}\left[-\frac{2s}{c}r_A-\frac{2f_A\left(t\right)}{c}v_A\right]\\ S_B=T_{\mathrm{LFM}}\left[f_{\mathrm{rB}}\left(t\right)-f_B\left(t\right)\right]\\ \simeq T_{\mathrm{LFM}}\left[-\frac{2{\mathrm{sr}}_B}{c}-\frac{2v_B}{c}f_B\left(t\right)\right]\\ =T_{\mathrm{LFM}}\left[-\frac{2s}{c}\left(r_A+v_A\tau +a{\tau }^2\text{/}2\right)-\frac{2}{c}\left(v_A+a\tau \right)f_B\left(t\right)\right]\\ =T_{\mathrm{LFM}}\left\{\frac{-2s}{c}r_A-\frac{2}{c}\left[f_B\left(t\right)+s\tau \right]v_A-\frac{2}{c}\left[\tau f_B\left(t\right)+\frac{s{\tau }^2}{2}\right]a\right\}\end{array}$
• where
• τ_A, v_Aare respectively the range and velocity of theAAM20 at t=(2n−1)·τ;
• τ_B, v_Bare respectively the range and velocity of theAAM20 at t=2n·τ.
• Further, the phase difference Δφ between the second LFM signal measured at t=2n·τ and the first LFM signal measured at t=(2n−1)·τ can be obtained from the following equation:
• $\begin{array}{c}\mathrm{\Delta \varphi }={\varphi }_B-{\varphi }_A\\ =\frac{4\pi }{c}\{\left[f_B\left(0\right)+2\mathrm{ns}\tau -\frac{2{\mathrm{sr}}_B}{c}\right]\left(1-\frac{2v_B}{c}\right)r_B-\\ \left[f_A\left(0\right)+\left(2n-1\right)s\tau -\frac{2{\mathrm{sr}}_A}{c}\right]\left(1-\frac{2v_A}{c}\right)r_A\}\\ \simeq \frac{4\pi }{c}\{\left[f_B\left(0\right)+2\mathrm{ns}\tau -\frac{2s{\tau }_B}{c}\right]r_B-\\ \left[f_A\left(0\right)+\left(2n-1\right)s\tau -\frac{2s{\tau }_A}{c}\right]r_A\}\\ \simeq \frac{4\pi }{c}f_{\mathrm{shift}}{\tau }_A+\frac{4\pi }{c}f_B\left(t\right)\left(v_A\tau +\frac{1}{2}a{\tau }^2\right)\end{array}$
• Therefore, it can be deduced that:
• $\left[\begin{array}{c}r\\ v\\ a\end{array}\right]=\left[\begin{array}{c}{r}_A\\ v_A\\ a\end{array}\right]={\stackrel{\_}{J}}^{-1}·\left[\begin{array}{c}-S_AC\text{/}T_{\mathrm{LFM}}\\ -S_BC\text{/}T_{\mathrm{LFM}}\\ \mathrm{\Delta \varphi }C\text{/}\pi \end{array}\right]$
• where
• $\stackrel{\_}{J}=\left[\begin{array}{ccc}2s& 2f_A\left(t\right)& 0\\ 2s& 2\left[f_B\left(t\right)+s·\tau \right]& 2\left[f_B\left(t\right)·\tau +s·{\tau }^2\text{/}2\right]\\ 4f_{\mathrm{shif}}t& 4f_B\left(t\right)·\tau & 2f_B\left(t\right)·{\tau }^2\end{array}\right]$
• Moreover, if we choose ƒ_shift=−ƒ_inc/2, a higher level of accuracy can be achieved for the range and velocity measurement. The output N datasets {r₁, v₁, a₁}, {r₂, v₂, a₂}, . . . , and {r_N, v_N, a_N} from this radarsignal sensor array120 are then transferred respectively to the N Gaussian-noise filters in the first Gaussian-noise filter array130 for S/N enhancement.
First Gaussian-Noise Filter Array130
• The first Gaussian-noise filter array130 includes N Gaussian-noise filters which are preferably implemented with an array of N one-stage linear Kalman filters, each being coupled to one of the N hybrid FSK/LFM sensors in the radarsignal sensor array120.
• These N one-stage linear Kalman filters are capable of reducing the Gaussian-noise in the output datasets {r₁, v₁, a₁}, {r₂, v₂, a₂}, . . . , and {r_N, v_N, a_N} from the N hybrid FSK/LFM sensors in the radarsignal sensor array120. These N one-stage linear Kalman filters are based on a conventional Gaussian-noise filtering technology, so detailed description thereof will not be given in this specification. References about the internal structure and input-output characteristics of the one-stage linear Kalman filtering method can be found, for example, in the book “Adaptive Filter Theory”, fourth edition, authored by S. Haykin and published by Prentice Hall, 2002; and the technical paper entitled “RADAR TRACKING FOR AIR SUREILLANCE IN A STRESSFUL ENVIRONMENT USING A FUZZY-GAIN FILTER” by K. C. C. Chan et al, IEEE Trans. Fuzzy Syst. vol. 5, no. 1, pp. 80-89, June 1997. The None-stage linear Kalman filters in the first Gaussian-noise filter array130 are used in combination to produce a collection of N S/N-enhanced sets of position-motion data about theAAM20, which are represented by {{circumflex over (r)}₁, {circumflex over (v)}₁, â₁}, {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂}, . . . , and {{circumflex over (r)}_N, {circumflex over (v)}_N, â_N}, respectively. These S/N-enhanced datasets {{circumflex over (r)}₁, {circumflex over (v)}₁, â₁}, {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂}, . . . , and {{circumflex over (r)}_N, {circumflex over (v)}_N, â_N} are then transferred to thetrilateration module140 for further processing.
Trilateration Module140
• Thetrilateration module140 is used for processing the N S/N-enhanced datasets {{circumflex over (r)}₁, {circumflex over (v)}₁, â₁}, {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂}, . . . , {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂} from the first Gaussian-noise filter array130 in reference to a predefined 3-dimenional (3-D) coordinate system, preferably a 3-D rectangular coordinate system having an x-axis, a y-axis, and a z-axis, to thereby obtain a collection of 3 axis-oriented sets of position-motion data respectively in reference to the x-axis, the y-axis, and the z-axis of the 3-D rectangular coordinate system. The trilateration process yields a set of position data {{circumflex over (x)}, ŷ, {circumflex over (z)}}, a set of velocity data {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}, and a set of acceleration data {â_x, â_y, â_z} in accordance with the following equations (A1), (A2), and (A3):
• $\begin{array}{cc}\hat{x}=\frac{x_1^2-x_2^2-{\hat{r}}_1^2+{\hat{r}}_2^2}{2·\left(x_1-x_2\right)}\hat{y}=\frac{-p_2\sqrt{p_2^2-4{\mathrm{<figure-callout\; label="p"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">p</figure-callout>}}_1p_3}}{2{\mathrm{<figure-callout\; label="p"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">p</figure-callout>}}_1}\hat{z}=p_4-\frac{y_3}{z_3}·\hat{y}& \mathrm{Equation}\left(A1\right)\\ \left[\begin{array}{c}{\hat{v}}_x\\ {\hat{v}}_y\\ {\hat{v}}_z\end{array}\right]={\stackrel{\_}{L}}^{-1}·\left[\begin{array}{c}{\hat{v}}_1·{\hat{r}}_1\\ {\hat{v}}_2·{\hat{r}}_2\\ {\hat{v}}_3·{\hat{r}}_3\end{array}\right]\mathrm{where}\stackrel{\_}{L}=\left[\begin{array}{ccc}\hat{x}-x_1& \hat{y}& \hat{z}\\ \hat{x}-x_2& \hat{y}& \hat{z}\\ \hat{x}& \hat{y}-y_3& \hat{z}-z_3\end{array}\right]& \mathrm{Equation}\left(A2\right)\\ \left[\begin{array}{c}{\hat{a}}_x\\ {\hat{a}}_y\\ {\hat{a}}_z\end{array}\right]={\stackrel{\_}{L}}^{-1}·\left[\begin{array}{c}{\hat{a}}_1·{\hat{r}}_1\\ {\hat{a}}_2·{\hat{r}}_2\\ {\hat{a}}_3·{\hat{r}}_3\end{array}\right]& \mathrm{Equation}\left(A3\right)\end{array}$
• The above-listed equations (A1), (A2), and (A3) are deduced as follows. Assume that the location of the (i)th hybrid FSK/LFM sensor in the radarsignal sensor array120 is (xi, yi, zi), the originally-emitted radar signal beam hits theAAM20 and reflects back at ({circumflex over (x)}, ŷ, {circumflex over (z)}) (which represented the detected position of the AAM20), and theAAM20 moves at a velocity of {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}). Then, in the case of N=3, the 3 hybrid FSK/LFM sensors in the radarsignal sensor array120 can be used to respectively acquire three samples of range data {{circumflex over (r)}₁, {circumflex over (r)}₂, {circumflex over (r)}₃}, which are mathematically expressed as follows:
• 
  {circumflex over (r)}₁²=({circumflex over (x)}−x₁)²+ŷ²+{circumflex over (z)}²  (1)
• 
  {circumflex over (r)}₂²=({circumflex over (x)}−x₂)²+ŷ²+{circumflex over (z)}²  (2)
• 
  {circumflex over (r)}₃²={circumflex over (x)}²+(ŷ−y₃)²+({circumflex over (z)}−z₃)²  (3)
• Note that y₁=z₁=y₂=z₂=x₃=0. From (1) and (2), it can be deduced that:
• 
  {circumflex over (r)}₁²−({circumflex over (x)}−x₁)²={circumflex over (r)}²−({circumflex over (x)}−x₂)²
• Therefore, it can be deduced that:
• $\hat{x}=\frac{x_1^2-x_2^2-{\hat{r}}_1^2+{\hat{r}}_2^2}{2·\left(x_1-x_2\right)}$$\hat{y}=-\frac{-p_2+\sqrt{p_2^2-4{\mathrm{<figure-callout\; label="p"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">p</figure-callout>}}_1p_3}}{2{\mathrm{<figure-callout\; label="p"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">p</figure-callout>}}_1}$$\mathrm{where}$${\mathrm{<figure-callout\; label="p"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">p</figure-callout>}}_1=1+\frac{y_3^2}{z_3^2}$$p_2=-2p_4·\frac{y_3}{z_3}$$p_3=p_4^2-{\hat{r}}_2^2+{\left(\hat{x}-x_2\right)}^2$$p_4=\frac{r_2^2-{\left(\hat{x}-x_2\right)}^2-{\hat{r}}_3^2+{\hat{x}}^2+y_3^2+z_3^2}{2z_3}$
• Further, it can be deduced that {{circumflex over (v)}₁, {circumflex over (v)}₂, {circumflex over (v)}₃} is mathematically related to {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z} as follows:
• ${\hat{v}}_1=\frac{\hat{x}-x_1}{{\hat{\mathrm{<figure-callout\; label="r\; ^"\; filenames="US20090091490A1-20090409-D00002.png"\; state="\{\{state\}\}">r</figure-callout>}}}_1}·{\hat{v}}_x+\frac{\hat{y}}{{\hat{r}}_1}·{\hat{v}}_y+\frac{\hat{z}}{{\hat{r}}_1}·{\hat{v}}_z$${\hat{v}}_2=\frac{\hat{x}-x_2}{{\hat{r}}_2}·{\hat{v}}_x+\frac{\hat{y}}{{\hat{r}}_2}·{\hat{v}}_y+\frac{\hat{z}}{{\hat{r}}_2}·{\hat{v}}_z$${\hat{v}}_3=\frac{\hat{x}}{{\hat{r}}_3}·{\hat{v}}_x+\frac{\hat{y}-{\mathrm{<figure-callout\; label="y"\; filenames="US20090091490A1-20090409-D00003.png"\; state="\{\{state\}\}">y</figure-callout>}}_3}{{\hat{r}}_2}·{\hat{v}}_y+\frac{\hat{z}-{\mathrm{<figure-callout\; label="z"\; filenames="US20090091490A1-20090409-D00003.png"\; state="\{\{state\}\}">z</figure-callout>}}_3}{{\hat{r}}_2}·{\hat{v}}_z$
• From the above equations, it can be obtained that:
• $\left[\begin{array}{c}{\hat{v}}_x\\ {\hat{v}}_y\\ {\hat{v}}_z\end{array}\right]={\stackrel{\_}{L}}^{-1}·\left[\begin{array}{c}{\hat{v}}_1·{\hat{r}}_1\\ {\hat{v}}_2·{\hat{r}}_2\\ {\hat{v}}_3·{\hat{r}}_3\end{array}\right]\mathrm{where}\stackrel{\_}{L}=\left[\begin{array}{ccc}\hat{x}-x_1& \hat{y}& \hat{z}\\ \hat{x}-x_2& \hat{y}& \hat{z}\\ \hat{x}& \hat{y}-y_3& \hat{z}-z_3\end{array}\right]$
• and similarly, it can be deduced that {â_x, â_y, â_z} is related to {â₁, â₂, â₃} as follows:
• $\left[\begin{array}{c}{\hat{a}}_x\\ {\hat{a}}_y\\ {\hat{a}}_z\end{array}\right]={\stackrel{\_}{L}}^{-1}·\left[\begin{array}{c}{\hat{a}}_1·{\hat{r}}_1\\ {\hat{a}}_2·{\hat{r}}_2\\ {\hat{a}}_3·{\hat{r}}_3\end{array}\right]$
• One problem in the use of this trilateration process, however, is that it would undesirably cause the 3 outputted datasets {{circumflex over (x)}, ŷ, {circumflex over (z)}}, {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}, and {â_x, â_y, â_z} to be degraded in S/N. As a solution to this problem, these 3 datasets {{circumflex over (x)}, ŷ, {circumflex over (z)}}, {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}, {â_x, â_y, â_z} are rearranged into three groups: {{circumflex over (x)}, {circumflex over (v)}_x, â_x}, {ŷ, {circumflex over (v)}_y, â_y}, {{circumflex over (z)}, {circumflex over (v)}_z, â_z}, which are then transferred to the second Gaussian-noise filter array150 for S/N enhancement.
Second Gaussian-Noise Filter Array150
• The second Gaussian-noise filter array150 includes an array of at least 3 Gaussian-noise filters which are preferably implemented with 3 one-stage linear Kalman filters including a first Gaussian-noise filter151, a second Gaussian-noise filter152, and a third Gaussian-noise filter153. The first Gaussian-noise filter151 is used for S/N enhancement of the x-axis oriented dataset {{circumflex over (x)}, {circumflex over (v)}_z, â_x} to thereby obtain a noise-reduced dataset, here represented by {{tilde over (x)}, {tilde over (v)}_x, ã_x}; the second Gaussian-noise filter152 is used for S/N enhancement of the y-axis oriented set of position-motion data {ŷ, {circumflex over (v)}_y, â_y} to thereby obtain a noise-reduced dataset, here represented by {{tilde over (y)}, {tilde over (v)}_y, ã_y}; and the third Gaussian-noise filter153 is used for S/N enhancement of the z-axis oriented set of position-motion data {{circumflex over (z)}, {circumflex over (v)}_z, â_z} to thereby obtain a noise-reduced dataset, here represented by {{tilde over (z)}, {tilde over (v)}_z, ã_z}.
RadarData Presentation Module160
• Theradar presentation module160 is used to present the output S/N-enhanced datasets {{tilde over (x)}, {tilde over (v)}_x, ã_x}, {{tilde over (y)}, {tilde over (v)}_y, ã_y}, {{tilde over (z)}, {tilde over (v)}_z, ã_z} in human-coginzable data format for the purpose of informing the pilot of thejetfighter10 of the position, speed, and acceleration of the threateningAAM 20. First, the output datasets {{tilde over (x)}, {tilde over (v)}_x, ã_x}, {{tilde over (y)}, {tilde over (v)}_y, ã_y}, {{tilde over (z)}, {tilde over (v)}_z, ã_z} from the second Gaussian-noise filter array150 are rearranged into {{tilde over (x)}, {tilde over (y)}, {tilde over (z)}} for position, {{tilde over (v)}_x, {tilde over (v)}_y, {tilde over (v)}_z} for velocity, and {ã_x, ã_y, ã_z} for acceleration. These datasets are then used to drive, for example, a monitor screen (not shown) for visually presenting the position, speed, and acceleration of the threateningAAM 20 to the pilot of thejetfighter10.
Operation of the Invention
• The following is a detailed description of the operation of the moving target radar tracking system of theinvention100. In this application example, it is assumed that the moving target radar tracking system of theinvention100 is installed on ajetfighter10 and, during flight of thejetfighter10, anAAM20 is launched against thejetfighter10.
• During flight of thejetfighter10, the radarsignal sensor array110 is activated to emit a beam of radar signal into the space under surveillance by thejetfighter10. When the emitted radar signal beam hits theAAM20, it will reflect as an echoed radar signal back to thejetfighter10 and which is sampled by the N hybrid FSK/LFM sensors in the radarsignal sensor array120 to produce a collection of N datasets {r₁, v₁, a₁}, {r₂, v₂, a₂}, . . . and {r_N, v_N, a_N}.
• Subsequently, these N datasets {r₁, v₁, a₁}, {r₂, v₂, a₂}, . . . , and {r_N, v_N, a_N} are transferred to the first Gaussian-noise filter array130, which includes an array of N one-stage linear Kalman filters, for S/N enhancement before undergoing trilateration. The S/N-enhanced datasets {{circumflex over (r)}₁, {circumflex over (v)}₁, â₁}, {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂}, . . . , and {{circumflex over (r)}_N, {circumflex over (v)}_N, â_N} are then transferred to thetrilateration module140, which is capable of processing the N S/N-enhanced datasets {{circumflex over (r)}₁, {circumflex over (v)}₁, â₁}, {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂}, . . . , {{circumflex over (r)}₂, {circumflex over (v)}₂, â₂} in reference to a 3-D rectangular coordinate system to thereby obtain a collection of 3 axis-oriented sets of position-motion data respectively in reference to the x-axis, the y- axis, and the z-axis of the 3-D rectangular coordinate system. The trilateration process yields a set of position data {{circumflex over (x)}, ŷ, {circumflex over (z)}}, a set of velocity data {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}, and a set of acceleration data {â_x, â_y, â_z}.
• Since the foregoing trilateration process might increase the level of noise in the resulted data {{circumflex over (x)}, ŷ, {circumflex over (z)}}, {{circumflex over (v)}_x, {circumflex over (v)}_y, {circumflex over (v)}_z}, and {â_x, â_y, â_z}, these data are rearranged into three groups: {{circumflex over (x)}, {circumflex over (v)}_x, â_x}, {ŷ, {circumflex over (v)}_y, â_y}, and {{circumflex over (z)}, {circumflex over (v)}_z, â_z} for further processing by the second Gaussian-noise filter array140 for S/N enhancement.
• In the second Gaussian-noise filter array150, the first Gaussian-noise filter151 is used for S/N enhancement of the x-axis oriented dataset {{circumflex over (x)}, {circumflex over (v)}_x, â_x}; the second Gaussian-noise filter152 is used for S/N enhancement of the y-axis oriented dataset {ŷ, {circumflex over (v)}_y, â_y}; and the third Gaussian-noise filter153 is used for S/N enhancement of the z-axis oriented dataset {{circumflex over (z)}, {circumflex over (v)}_z, â_z}. The resulted S/N-enhanced datasets are represented by {{tilde over (x)}, {tilde over (v)}_x, ã_x}, {{tilde over (y)}, {tilde over (v)}_y, ã_y}, and {{tilde over (z)}, {tilde over (v)}_z, ã_z}, respectively.
• Finally, the S/N-enhanced datasets {{tilde over (x)}, {tilde over (v)}_x, ã_z}, {{tilde over (y)}, {tilde over (v)}_y, ã_y}, {otl z, {tilde over (v)}_z, ã_z} are transferred to the radardata presentation module160, which rearranges these datasets into {{tilde over (x)}, {tilde over (y)}, {tilde over (z)}} for position, {{tilde over (v)}_x, {tilde over (v)}_y, {tilde over (v)}_z} for velocity, and {ã_x, ã_y, ã_z} for acceleration, and then uses these position, velocity, and acceleration datasets {{tilde over (x)}, {tilde over (y)}, {tilde over (z)}}, {{tilde over (v)}_x, {tilde over (v)}_y, {tilde over (v)}_z}, {ã_x, ã_y, ã_z} to drive a radar monitor screen (not shown) for visually presenting the position, speed, and acceleration of the threateningAAM20 to the pilot of thejetfighter10.
• In conclusion, the invention provides a method and system for radar tracking of a moving target (such as an air-to-air missile) from a moving station (such as a jetfighter) with improved accuracy and performance. The proposed method and system is characterized by the use of a hybrid FSK/LFM scheme for acquiring a collection of raw radar data, a first Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the raw radar data, a trilateration module, and a second Gaussian-noise filter array of one-stage linear Kalman filters for S/N-enhancement of the trilateration-resulted radar data. These features allow the radar tracking of moving objects to be more fast and accurate. The invention is therefore more advantageous to use than the prior art.
• The invention has been described using exemplary preferred embodiments. However, it is to be understood that the scope of the invention is not limited to the disclosed embodiments. On the contrary, it is intended to cover various modifications and similar arrangements. The scope of the claims, therefore, should be accorded the broadest interpretation so as to encompass all such modifications and similar arrangements.

Claims (12)

1. A moving target radar tracking method, comprising:

emitting a radar signal beam;

sensing an echo of the radar signal beam from a moving target and producing a raw set of radar data about the moving target, including range, velocity, and acceleration;

performing S/N (signal-to-noise ratio) enhancement of each raw set of radar data;

performing trilateration on the S/N-enhanced sets of radar data in reference to a predefined 3-dimenional coordinate system having a first axis, a second axis, and a third axis to obtain three axis-oriented sets of radar data respectively in reference to the first axis, the second axis, and the third axis of the 3-dimenional coordinate system;

performing S/N enhancement on the trilateration-resulted sets of radar data; and

presenting the S/N-enhanced trilateration-resulted sets of radar data in a human-cognizable data form.

2. The moving target radar tracking method ofclaim 1, wherein the sensing of the echo of radar signal is implemented with a hybrid FSK/LFM scheme.

3. The moving target radar tracking method ofclaim 1, wherein the S/N-enhancement of the raw sets of radar data is implemented with a Kalman filtering process.

4. The moving target radar tracking method ofclaim 3, wherein the Kalman filtering process is based on a one-stage linear scheme.

5. The moving target radar tracking method ofclaim 1, wherein the S/N-enhancement of the trilateration-resulted sets of radar data is implemented with a Kalman filtering process.

6. The moving target radar tracking method ofclaim 5, wherein the Kalman filtering process is based on a one-stage linear scheme.

7. A moving target radar tracking system, comprising:

a radar signal emitter, which is capable of emitting a radar signal beam;

a radar signal sensor array, which includes an array of radar signal sensors, each capable of sensing an echo of the radar signal beam from a moving target and producing a raw set of radar data about the moving target, including range, velocity, and acceleration;

a first Gaussian-noise filter array, which includes an array of Gaussian-noise filters respectively coupled to the radar signal sensors in the radar signal sensor array for S/N (signal-to-noise ratio) enhancement of each raw set of radar data produced by the radar signal sensor array;

a trilateration module for trilaterating the S/N-enhanced sets of radar data outputted from the first Gaussian-noise filter array in reference to a predefined 3-dimenional coordinate system having a first axis, a second axis, and a third axis to obtain three axis-oriented sets of radar data respectively in reference to the first axis, the second axis, and the third axis of the 3-dimenional coordinate system;

a second Gaussian-noise filter array, which includes an array of at least three Gaussian-noise filters respectively used for S/N-enhanced of the set of first-axis oriented radar data, the set of second-axis oriented radar data, and the set of third-axis oriented radar data produced by the trilateration module; and

a radar data presentation module, which is capable of presenting the S/N-enhanced trilateration-resulted sets of radar data produced by the second Gaussian-noise filter array in a human-cognizable data form.

8. The moving target radar tracking system ofclaim 7, wherein the radar signal sensor array is composed of an array of hybrid FSK/LFM sensors.

9. The moving target radar tracking system ofclaim 7, wherein the first Gaussian-noise filter array is composed of an array of Kalman filters.

10. The moving target radar tracking system ofclaim 9, wherein the Kalman filters are one-stage linear type.

11. The moving target radar tracking system ofclaim 7, wherein the second Gaussian-noise filter array is composed of an array of Kalman filters.

12. The moving target radar tracking system ofclaim 7, wherein the Kalman filters are one-stage linear type.

Images