Using radar targets of opportunity to build a monopulse calibration table

Abstract

In a radar system, a monopulse calibration table is constructed from live targets of opportunity. A center of gravity or weighted average of normalized signals ΔV received at SUM and DIFF channels from a live target are used to determine the target's actual azimuth. Off bore sight angles (OBA) of the target are then determined from the target's actual azimuth. Normalized received signal values of ΔV are converted to nearest-valued integers. The OBA s that correspond to each integer-valued normalized received signal are averaged and can then be plotted as a function of normalized received signal value ΔV. Different tables or plots can be constructed for elevation angles. An equation of a best-fit line the matches or at least closely approximates the plotted data is determined to smooth the actual data.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram of a radar system and a target of opportunity;

FIG. 2 is a chart showing the entire process of generating the calibration table which includes determining a target's off-bore sight location from the strength of signals received at SUM and DIFF antennas of the antenna, the calculation of the table, storing the tables in the arrays, calculating the best fit line and building the final calibration table;

FIG. 3 is a graph of data obtained from an actual target having strong SUM and DIFF channel signals;

FIG. 4 is a graph of data obtained from an actual target having weak SUM and DIFF channel signals;

FIG. 5 is a graph of data obtained from an actual target having noisy SUM and DIFF channel signals that may result from interference;

FIG. 6 is block diagram showing steps of another embodiment of the method disclosed herein; and

FIG. 7 is a graph of an actual calibration table measured data and best fit line superimposed on it calibration table (one out of 5 used) where the calibration angle is plotted as a function of the monopulse discriminant.

FIG. 8 shows a plot of the data of a monopulse calibration table obtained from live data, and a best-fit line approximation of the live data, derived from the live data.

DETAILED DESCRIPTION

FIG. 1 is a depiction of a radar system 10 comprised of a rotating antenna 12 having two input channels or input patterns known in the art as the SUM 13A channel and a DIFFERENCE or DIFF channel 13B. The antenna 12 has a beam width 14, the center line or bore sight of which is identified by reference numeral 16. As the antenna 12 rotates about a vertical axis, the bore sight 16 and a reference direction such as true north form an angle 18 referred to herein as azimuth. In one embodiment using SSR, the target of opportunity 11 has a transponder that emits a signal when the transponder detects a signal from a radar transmitter. When PSR is used, the radar signal returned by the target surface is used for the monopulse process. For each signal or pulse that the transmitter 20 emits at a given azimuth angle θ, the transponder in the target 11 emits a response signal or the target's surface returns the signal. Thus, there are SUM and DIFF signals produced at each of several different azimuth angles θ.

FIG. 2 shows a chart that depicts the steps of the general method 200 of building a monopulse calibration table using a radar target of opportunity 11 with a radar system or radar receiver, such as the system and receiver shown in FIG. 1. FIG. 2 also shows the signals used by the method.

In a first step that is identified in FIG. 2 by reference numeral 202, monopulse data is collected from a target and stored in a computer memory device, such as random access memory (RAM) or magnetic or optical disc. A computer memory storage device is identified in FIG. 1 by reference numeral 24. The collected data for each of several azimuth angles θ provides a monopulse estimate, which is depicted in FIG. 2 by reference numeral 206. The monopulse estimate is defined as:

ΔV=Log(SUM)−Log(DIFF)  (1)

where SUM and DIFF are the amplitudes of the signals received at the sum and difference channels of an antenna respectively.

In addition to storing a monopulse estimate, a beam side indication is also stored. The beam side indication and its storage is identified in FIG. 2 by reference numeral 208.

The antenna's azimuth angle θ is stored, the value and storage of which is identified in FIG. 2 by reference numeral 210. The antenna's azimuth angle θ is the angle of the antenna where the SUM and DIFF signals used in the monopulse estimate were obtained.

If it is available, the target's height is stored. In FIG. 2, the target's height and its storage is identified by reference numeral 212.

The target's range is stored, which is identified by reference numeral 214 and the radar's interrogation mode is also stored, which is identified by reference numeral 216.

In a second step of the method of building a monopulse calibration table from a live target, the target's actual location azimuth is determined and the azimuth readings obtained from a target are normalized. The azimuth angle readings from a target are normalized by subtracting the actual azimuth angle from the azimuth readings where the sum and diff signals were obtained. In FIG. 2, the process of determining off bore sight angles of a live target is identified by reference numeral 220. The sub-process of determining a live target's actual azimuth or angle center is identified by reference numeral 222 and the process of normalizing other off bore sight azimuth angles θ is identified by reference numeral 224.

In building a monopulse calibration table from a live target, the target in one example is swept several times. After the last sweep, the data that was collected in the step identified by reference numeral 202 is first used to calculate the target's azimuth center, which is labeled “angle center” in FIG. 2. The calculation of a target's actual azimuth angle θ_boresight or “angle center” in FIG. 2 can be determined two ways.

In a simple case where the SUM and DIFF signals from a target are strong and the beam side indicator (embodied as a DIFF signal phase change) is unambiguous, a target's actual azimuth angle θ_boresight (identified as the angle center in FIG. 2) can be determined by first computing ΔV for each of several different azimuth angles θ using equation (1) above. As a second step of the process 222 of determining a target's actual azimuth or angle center, the values of ΔV for each azimuth angle are compared to each other. The azimuth angle θ where ΔV is maximum is considered to be the azimuth angle θ_boresight where the target is actually located.

FIG. 3 is a graph of ΔV signals at various azimuth angles that were obtained from an actual live target having strong SUM and DIFF channel signals. As can be seen in FIG. 3, there is a distinct maximum ΔV at an azimuth angle of approximately 165.5 degrees (for the specific target represented in FIG. 3), which is also where the DIFF signal phase angle or beam side indicator (identified by reference numeral 208 in FIG. 2) changes from one to zero, indicating a 180 degree phase change. Once the target's actual azimuth angle θ_boresight is determined, the other values of ΔV for the target can be correlated to the azimuth angles where every other ΔV was calculated. The other values of ΔV are correlated to corresponding off bore sight angles by calculating the arithmetic difference between θ_boresight and the azimuth angle θ where every other ΔV was determined. By plotting the value of ΔV obtained at θ_boresight and plotting the values of ΔV obtained at several other OBA as a function of OBA, a simple monopulse calibration table of ΔV values vs. off bore sight angles (OBAs) can be constructed from signals received from (or returned by) an actual target.

In a radar system where the SUM and DIFF channels are strong, the above-identified simple process of determining an azimuth angle where ΔV is maximum can work effectively as will the above-identified process of determining normalized azimuth values. FIG. 4, however, is a graph of ΔV signals at various azimuth angles that were obtained from an actual target having weak SUM and DIFF channel signals. As shown in FIG. 4, when signals from a target are weak, the graph of ΔV as a function of off bore sight azimuth angles θ can have many ΔV that can each appear to be maxima. Moreover, there may be more than one azimuth angle where the DIFF signal phase change is not distinct, as shown by the multiple beam side indicator transitions.

FIG. 5 shows a graph of actual ΔV signals at various azimuth angles calculated from a garbled or noisy target. As shown in FIG. 5, there does not appear to be any maxima of ΔV at any particular azimuth angle θ.

Together, FIGS. 3, 4 and 5 show that a single set of SUM and DIFF signals might not show an azimuth angle of where a target is actually located. Thus, in one embodiment, finding the target boresight using an actual target of opportunity uses weighted averages of ΔV throughout an entire beam width scan, using several ΔV signal levels at several azimuth angles. In such an embodiment, the process of determining the actual azimuth angle θ_boresight of an actual target (identified by reference numeral 222 in FIG. 2) is determined using a weighted average of all i azimuth angles and corresponding ΔV values. The weighted average value of θ_boresight in this embodiment is expressed as:

$\begin{array}{cc}{\theta }_{\mathrm{boresight}}=\frac{\sum _i{\theta }_i\Delta V_i}{\sum _i\Delta V_i}& \left(2\right)\end{array}$

where θ_i is an azimuth angle where a corresponding ΔVi is determined using equation (1) above.

θ_boresight is more reliably determined using a weighted average of several azimuth angles θ and the corresponding ΔV. After θ_boresight is determined, OBA azimuth readings can be determined for the other azimuth angles where ΔV was determined by subtracting θ_boresight from the azimuth angle where each ΔV was determined from equation (1) above.

FIG. 6 (FIGS. 6A and 6B) shows steps of a method for constructing a monopulse calibration table from ΔV signals received from several targets of opportunity at a single elevation angle, which can of course be repeated for targets at different elevation angles.

At step 602, a normalized received signal ΔVi is calculated for signals that are received at the antenna's SUM and DIFF channels from a target at each azimuth angle θ_i of several different i azimuth angles. At step 604, the target's actual azimuth angle θ_boresight is calculated as the sum of the product of all θ_i·ΔV_i, where θ_i is the azimuth angle where a particular ΔV_i is determined, divided by the sum of all ΔV_i.

In step 606, for each of the i-number of azimuth angles θ where a ΔV_i is determined, the off bore sight angle θ_{off bore sight} where each ΔV_i was determined is calculated as the arithmetic difference between each azimuth angle θ_i, of the first step and the azimuth angle θ_boresight determined in the second step to obtain i, off-bore sight angles, θ_{off bore sight i}, whereat each of the normalized signals ΔV_i of step a) was determined. The normalized received signals ΔV_i obtained in the first step are correlated to corresponding off-bore sight angles (OBA) expressed as θ_{off bore sight}. In step 608, the values of ΔV_i determined at each off bore sight angle are stored and in step 610, a running average of each the off bore sight angles for each ΔV_i is calculated.

As can be seen in from Eq. (2) above, when determining θ_boresight using a weighted average, many azimuth angles are included in the determination of θ_boresight during which, some readings may cross North or be detected as having crossed North because of the antenna beam width. Therefore, when determining θ_boresight using a weighted average, a test is first made to determine if a reading crossed North.

The determination of whether a reading or readings crossed North during the collection of azimuth angle measurements is made by determining whether the arithmetic difference between the last azimuth angle collected (expressed as θ_end) and the first azimuth angle collected (expressed as θ_start) is less than zero. This test can be expressed as:

IF(θend − θstart) < 0 then
target_North_crossing == true;

If the readings crossed North during the collection of azimuth angles, 360 degrees is added to any azimuth angle value measured to be less than 180 degrees, when determining the weighted average of θ_boresight. This test of azimuth angles and the adjustment to measured azimuth angles θ less than 180 degrees when the readings crossed North can be stated as:

IF(target_North_crossing == true && θi < 180)
then θi = θi + 360;

If the azimuth center θ_boresight as determined by using equation (2) above is calculated to be greater than 360 degrees, then 360 degrees needs to be subtracted from the calculated θ_boresight. This step can be stated as:

IF(θboresight > 360)
then θboresight = θboresight − 360;

After the azimuth angles and normalized received signals are accumulated, and after θ_boresight is determined, a table or graph of the normalized received signal values and an average of corresponding off bore sight azimuth angles can be created for each normalized received signal ΔV. Each OBA azimuth angle θ_i to be entered into the table is determined from the expression:

θ_{off boresight}=θ_boresight−θ_i  (3)

with values on the right side of the beam being assigned a positive value and values on the left side of the beam being assigned a negative value. For a reading crossing North, a correction is applied, the value of which is determined by whether θ_i is greater than 180 degrees or less than −180 degrees. The correction is expressed as:

If (θoff boresight > 180); θoff boresight = θoff boresight −360; and
If (θoff boresight < −180); θoff boresight = θoff boresight + 360;

In the course of building or populating a table using data from multiple targets and multiple sweeps, one target might provide say, 20 ΔV values obtained at 20 different azimuth angles. A second target might provide a different set or sets of ΔV values at its corresponding azimuth angles. Third, fourth, fifth, etc. targets might each provide fewer or greater numbers of ΔV values at again different azimuth angles. It is also likely that there will be multiple azimuth angles for a particular value of ΔV. Put another way, for a table with 512 values of ΔV, any particular value of ΔV might be obtained from different OBA. In one embodiment, each OBA entry of a monopulse calibration table is determined as a running average of azimuth angles where each value of ΔV was obtained. The running average θ(n)_count of an OBA for a particular value of ΔV can be calculated in one more than one way but in the above-identified embodiment it is calculated using running average where the new average, is calculated from the last average and the new azimuth input OBA angle θ_i:

$\begin{array}{cc}{\stackrel{\_}{\theta }\left(n\right)}_{\mathrm{Count}+1}=\frac{\mathrm{Count}*{\stackrel{\_}{\theta }\left(n\right)}_{\mathrm{count}}+{\theta }_i}{\mathrm{Count}+1}& \left(7\right)\end{array}$

where θ(n)_Count is a previously-determined arithmetic average of a Count-number of azimuth angles θ, and θ_i is a new or additional azimuth angle to be added into the average θ(n)_count+1 value of a Count+1-number of azimuth angles. Thus each of the off bore sight angle averages for each ΔV can be plotted on a graph or left in a table format.

Moreover, the method of building a monopulse calibration table from live targets in one example includes building multiple tables, with each table being for ΔV values and azimuth angles obtained a target at a particular elevation. In one embodiment, five separate monopulse calibration tables were created. Table 0 was for readings obtained from a target or targets at elevation angles from 0 to 10 degrees; Table 1 was for readings obtained at elevation angles of 10 to 20 degrees; Table 2 was for readings obtained at elevation angles from 20 to 30 degrees; Table 3: 30 to 40 degrees; and Table 5: 40 to 50 degrees elevation.

In SSR radar, elevation angle of a target is provided by the Mode C response, which is identified in FIG. 2 by reference numeral 216. When 3-D radar is used, elevation angle can be extracted from the radar height data (identified in FIG. 2 by reference numeral 212) and the range data (identified in FIG. 2 by reference numeral 214).

Once a reported target reaches its end, the strong or high quality signals are evaluated. In one embodiment, if at least three height replies were within 100 feet of each other, the height is determined to be valid with the first height reply being used to determine the target's elevation. If the height data is not valid or non-existent, the target height is set to zero and all of the signals are processed into Table 0. Once target height is known, the target elevation angle relative to the horizon can be calculated as:

The relative elevation to the antenna's boresight elevation θ_ANT—ELEV is determined as:

θ_Elev=θ_Elev−θ_ANT—ELEV  (4)

The elevation angle θ_Elev of a target is used to determine which of the elevation-angle monopulse calibration tables to update using the processes described above. After all of the data from a particular target is collected, a particular table can be constructed and updated according to the elevation angle where the data is obtained. A target's azimuth center, θ_boresight, is calculated for each table using the equations above. The off boresight angle for each ΔV is also calculated for each table.

Referring again to FIG. 2, at step 228, the normalized received signal values ΔV_i for each table is the index to the memory, with their corresponding running average OBA to form monopulse calibration tables for different elevations.

Referring now to FIG. 8 there is shown a monopulse calibration table, which is constructed by plotting values for ΔV along the horizontal or x-axis with the corresponding off bore sight angle (OBA) running averages plotted on the vertical or y-axis. In constructing the table shown in FIG. 8, before any ΔV values are stored or plotted, the ΔV values are first converted to integers. The integer values of ΔV then become pointers or indices into the calibration tables for the different elevation angles.

The monopulse calibration table thus has a monopulse table index or MonoTableIndex, which is expressed from the monopulse discriminant:

MonoDiscriminant=Log₂SUM−Log₂DIFF+(float)0.5;  (5)

where the term (float)0.5 is a value required to be added or subtracted to convert the Log₂SUM−Log₂DIFF value, to the nearest-value integer. The monopulse calibration table index can then be expressed as:

MonoTableIndex=(int)MonoDiscriminant.  (6)

Put another way, the monopulse calibration table index is equal to ΔV, converted to the nearest-valued integer. For each monopulse calibration table index, there will exist a corresponding off boresight angle θ.

As can be seen in FIG. 8, the horizontal axes of the monopulse calibration tables becomes series of integers starting from either zero or 1, up to a user-specified number of ΔV values which the table will have, when enough data to fill the table has been captured.

A calibration table with a relatively large number of ΔV values from a large number of azimuth angles will provide greater resolution or accuracy in the determination of a target's off bore sight location than will a calibration table of a relatively small number of ΔV values obtained from a relatively small number of azimuth angles. In one embodiment, the monopulse calibration table had 512 values for ΔV, corresponding to 512 the SUM and DIFF signals obtained at different azimuth angles where each ΔV was determined. Half of the ΔV values, i.e., 256, were made at azimuth angles on each side of the antenna's bore sight. Of the 256 values of ΔV on each side of the antenna's bore sight there are corresponding average OBA, which in one example are determined using a running average of the OBA.

The size of a monopulse calibration table is a design choice but a large table of say, 512 entries will provide a monopulse calibration table with greater resolution than will a short table of say 20 or 30 entries. The number of entries in the table should be evenly divided on both sides of an antenna's boresight.

In order to construct a monopulse calibration table of a large number of entries, (e.g., 256 entries on each side of the antenna's boresight) ΔV measurements are needed for each table entry. By way of example, for a table having 512 entries, 512 separate monopulse table indices, with each index corresponding to a ΔV value, are needed.

FIG. 8 also shows a best-fit line approximation of the live data, derived from the live data using the process described above. The horizontal or x-axis of the graph in FIG. 8 has 512 separate data values numbered from 0 to 511, each of which corresponds to a ΔV obtained at a corresponding azimuth angle, θ. Two hundred fifty six (256) values are plotted on each side of the antenna's bore sight azimuth angle θ boresight. The OBA angles where each ΔV was obtained, are plotted along the y-axis.

Those of ordinary skill in the art will recognize that it might not be possible to obtain 512 valid ΔV values from 512 different azimuth angles θ in a single sweep of a single live target. Thus, completion of a relatively long monopulse calibration table can require the accumulation of live data from multiple targets, over multiple scans at perhaps multiple elevation angles (in embodiments where separate tables are prepared for separate elevation angle sectors).

In order to improve accuracy and to fill tables having entries that might not be populated by data obtained during actual measurements, a “best-fit” line can be determined using a second order polynomial. As used herein, the best-fit line is considered to be a line or curve that closely approximates, but which in one example is identical to the shape of the curve plotted through the actual data values of off bore sight angles as a function of ΔV values.

FIG. 8 shows a best-fit line identified by reference numeral 802, which was derived from the running average of OBA and normalized ΔV. As a first step in calculating a best-fit line, an approximation is performed for the data on each side of the θ_boresight. In an embodiment where 512 data points are collected, the approximation is made separately for the 256 data points on each side θ_boresight. As shown by reference numeral 230 in FIG. 2, a “best-fit” line in one example is calculated for each separate elevation table.

In order to determine a best fit line, a polynomial approximation table Az(i) is constructed. For an embodiment of 512 data points, Az(i) for i=0 to i<256 is calculated first as:

$\begin{array}{cc}\mathrm{Az}\left(i\right)=\frac{1}{a_0+a_1i+a_2i^2+a_3i^3}& \left(8\right)\end{array}$

Higher indices, e.g., 256 to 511, can use the same process if the function to be approximated is represented in a reverse order. To approximate this function, a direct polynomial combination needs to be established:

y _i(i)=a₀+a₁i+a₂i²+a₃i³ where  (9)

$\begin{array}{cc}{y}_i\left(i\right)=\frac{1}{{\mathrm{Az}}_i\left(i\right)}& \left(10\right)\end{array}$

The best-fit line coefficients a_i can be obtained by solving the matrix equation that returns a vector A that contains the best fin line coefficients:

[M]{right arrow over (A)}={right arrow over (Y)}  (11)

where

$\begin{array}{cc}\left[M\right]=\left[\begin{array}{cccc}m& \sum _mi& \sum _mi^2& \sum _mi^3\\ \sum _mi& \sum _mi^2& \sum _mi^3& \sum _mi^4\\ \sum _mi^2& \sum _mi^3& \sum _mi^4& \sum _mi^5\\ \sum _mi^3& \sum _mi^4& \sum _mi^5& \sum _mi^6\end{array}\right]\mathrm{and}\overrightarrow{A}=\left[\begin{array}{c}{a}_0\\ a_1\\ a_2\\ a_3\end{array}\right]\mathrm{and}\overrightarrow{Y}=\left[\begin{array}{c}\sum _m\frac{1}{{\mathrm{Az}}_i\left(i\right)}\\ \sum _m\frac{i}{{\mathrm{Az}}_i\left(i\right)}\\ \sum _m\frac{i^2}{{\mathrm{Az}}_i\left(i\right)}\\ \sum _m\frac{i^3}{{\mathrm{Az}}_i\left(i\right)}\end{array}\right]& \left(12\right)\end{array}$

The solution for the vector coefficient matrix {right arrow over (A)} is given as

{right arrow over (A)}=[M]⁻¹{right arrow over (Y)} where the one-over matrix expression

$M^{-1}=\frac{1}{M}\mathrm{adjA}$

and where |M| is the determinant of the square matrix M and adjA is the matrix formed by replacing each element with its cofactor and transposing the result. Once the 256 points best-fit table is calculated by the function, it is loaded back to a 512 array in the proper order.

A best fit line in one example is tested to be valid or invalid. In one embodiment, a best-fit line calculation is considered to be valid for any of the five monopulse tables if a table had at least fifty cells and each cell has at least ten entries. Those of ordinary skill will recognize however that other best-fit line validity criteria could be used in alternate embodiments.

In the above-identified embodiment, if a best-fit line is not valid using the above-identified criteria then the above-identified default table is used for all five elevation tables. If the standard deviation for Table 0 is calculated to be more than a default value of 5, the default table was used for Table 0. If the standard deviation for Table 0 is calculated to be less than or equal to 5, the best-fit line is constructed using a combination of measured data and best-fit line data described above. For tables 1 to 4, if the best-fit line is invalid, or if the best fit line standard deviation is more than 1, the best-fit line for Table 0 is used.

The standard deviation of the best fit line is calculated from the equation:

$\begin{array}{cc}\mathrm{Std}=\sqrt{\frac{\sum _N{\left[\mathrm{Fit}\left(i\right)-\mathrm{Az}\left(i\right)\right]}^2}{N}}& \left(13\right)\end{array}$

Where

Fit(i)=the fit line points as extracted from the best fit table, described below; and

Az(i)=the measured data points, i.e., the data points derived from actual signals.

A comparison of the calculated best-fit line to actual data shows that the best-fit line provides a good approximation of monopulse calibration data where live target data is noisy. When the data quality is high, however, and when noise is minimal, a small bias error in the best-fit line can be observed at the edges of the actual data. Therefore, a best-fit table should be a linear combination of the best-fit line as calculated above and measured data, as can be seen in FIG. 8. In one embodiment, the best fit line replaces all points where the count of data points is less than a default number of 30 points and the measured data is used elsewhere. The default number of 30 points is a design choice and other values could be used in alternate embodiments. Finally, it should be noted that in embodiments using multiple elevation tables, a best-fit line is calculated as a composite of all the tables, as shown in FIG. 2.

As for an apparatus to perform the method set forth herein, FIG. 1 shows a depiction of a radar system 10. The system 10 is comprised of a rotating antenna 12 that receives signals at SUM 13A and difference 13B antenna channels within an angular-shaped beam 14, the center of which is the antenna's bore-sight 16. The angle formed between a reference direction, such as true north and the bore-sight 16 is the azimuth angle θ, which is also identified by reference numeral 18.

A transmitter/receiver 20 is coupled to the antenna 12 and receives signals from the SUM 13A DIFF 13B channels of the antenna 12, which is known to those of ordinary skill in monopulse radar.

A computer 22 is coupled to the transmitter/receiver 20 and receives data from the receiver, which represents the signals that are received at the SUM 13A and DIFF 13B channels of the antenna 12. Computer program instructions are stored in a memory device 24 coupled to the computer 20 though an address and control bus 25. The memory device 24 can be implemented using any appropriate device, examples of which include static or dynamic RAM, EEPROM, magnetic or optical disk storage devices. A display device 26, such as a CRT or LCD, displays locations of the target 11 for an operator.

When the instructions in the memory device 24 are read by the computer 22 and executed, they cause the computer to perform the steps depicted in the FIGS. 2 and 6 and as described above, thereby providing an apparatus by which a monopulse calibration table can be constructed using live, radar targets of opportunity. More particularly, the program instructions stored in memory cause the computer to read SUM and DIFF signal values from the radar receiver, calculate normalized received values thereof and to identify ΔV maxima, perform center of gravity averaging, building a monopulse calibration table, by calculating the running averages of azimuth angles for monopulse calibration table indices and calculating the best-fit line for each of the tables.

The steps or operations and the apparatus described herein are examples of a method and apparatus by which a monopulse calibration table can be constructed using a live radar target of opportunity. There may be many variations to these steps or operations and variations to the apparatus without departing from the spirit of the invention without departing from the spirit of the invention and these are therefore considered to be within the scope of the invention as defined in the following claims.

Claims

1. A method of building a monopulse calibration table from a live target of opportunity using signals received from the target at sum and difference antenna channels of a rotating radar antenna that receives signals within an angular-shaped beam, the center of which is the antenna's bore-sight, the angle formed between a reference direction and the bore-sight being an azimuth angle, the method comprising the steps of: a) at each of a plurality of N different azimuth angles, calculating a normalized received signal from a target of opportunity using a signal received from the target at a sum channel and using a signal received from the target at a difference channel;b) identifying the azimuth angle of where the target of opportunity is actually located as the azimuth angle whereat the magnitude of the normalized received signal obtained in step a) is greatest;c) for each of the other N azimuth angles of step a), calculating the difference between each of the N azimuth angles of step a) and the azimuth angle determined in step b) to obtain N, off-bore sight angles where each of the normalized received signals of step a) was determined;d) storing each of the normalized received signal obtained at step a) with a corresponding off bore sight azimuth angles determined in step c).

2. The method of step 1 wherein the step of calculating a normalized received signal is comprised of the steps of: dividing the magnitude of the signal from the sum channel by the magnitude of the signal from the difference channel to obtain a value ΔV.

3. The method of claim 2 wherein the step of identifying the azimuth angle whereat the magnitude of the normalized received signal obtained in step a) is greatest is further comprised of identifying the azimuth angle where the target center exists by calculating a weighted average of all azimuth angles for all corresponding ΔV values.

4. The method of step 1 wherein the step of calculating a normalized received signal is comprised of the steps of: determining a logarithmic value of the magnitude of the signal from the difference channel;determining a logarithmic value of the magnitude of the signal from the sum channel;subtracting the logarithmic value of the magnitude of the signal from the difference channel from the logarithmic value of the signal from the sum channel.

5. In a radar system, a method of determining off-bore sight angles of a target using a plurality of signals received from the target at sum and difference antenna channels of a rotating radar antenna at each azimuth angle of a plurality of azimuth angles, the method comprising the steps of: a) at each azimuth angle θi of a plurality of i azimuth angles, calculating a normalized received signal ΔVi for signals that are received at the antenna from a target from a Vsum signal and a Vdiff signal;b) calculating the target's actual azimuth angle θboresight as the sum of the product of all θi·ΔVi, divided by the sum of all ΔVi;c) for each of the i azimuth angles θ, of step a), calculating an off bore sight angle θoff bore sight i for each corresponding ΔVi by calculating the difference between each azimuth angle θi, of step a) and the azimuth angle θboresight determined in step b) to obtain i, off-bore sight angles, θoff bore sight i, whereat a corresponding ΔVi of step a) was determined;d) correlating each of the normalized received signals ΔVi obtained at step a) to a corresponding off-bore sight angle θoff bore sight angle determined in step c).

6. The method of claim 5 further including the steps of: converting each value of ΔVi to a nearest-valued integer, int(ΔVi).

7. The method of claim 6 further including the steps of: storing each int(ΔVi) with a corresponding θoff bore sight angle to determine a target's off bore sight angle.

8. The method of claim 5 wherein the step of calculating a normalized received signal is comprised of the steps of: dividing the magnitude of the signal from the sum channel Vsum by the magnitude of the signal Vdiff from the difference channel.

9. The method of claim 5 wherein each normalized received signal is determined by the steps of: determining a logarithmic value of the magnitude of the signal Vdiff from the difference channel;determining a logarithmic value of the magnitude of the signal Vsum from the sum channel;determining the difference between the logarithmic value of the magnitude of the signal from the difference channel and the logarithmic value of the signal from the sum channel.

10. The method of claim 5 wherein the step of identifying the azimuth angle whereat the magnitude of the normalized received signal obtained in step a) is greatest is further comprised of the steps of: storing each of the N azimuth angles whereat a sum and difference signal ratio is determined;storing each of the N sum and difference signal ratios;sorting each of the N sum and difference signal ratios to determine the greatest ratio; andidentifying the corresponding azimuth angle where the greatest sum and difference signal ratio occurred.

11. The method of claim 5 further comprising the step of: calculating a running average of the azimuth angles θ for each corresponding ΔVi.

12. The method of claim 11, wherein the running average is calculated by the equation:

13. The method of claim 5 further comprising the step of: assigning a plurality of signals received from the target at sum and difference antenna channels of a rotating radar antenna at a plurality of azimuth angles to a predetermined vertical angle segment.

14. The method of claim 5 further comprising the step of: determining a best-fit line for a graph of all ΔV and all θoff bore sight angles.

15. In a monopulse radar receiver, a method of building a monopulse calibration table from SUM and DIFF signals received from a plurality of live radar targets, the table having N entries, each of the N entries comprising a monopulse discriminant and a comprising a corresponding azimuth angle whereat signals from the live targets are obtained, the method comprising the steps of: determining a first azimuth angle θ that corresponds to a first monopulse discriminant, from a first live target;determining a second azimuth angle θ that corresponds to the first monopulse discriminant, from the first live target;calculating a running average of the first and second azimuth angles θ for the first monopulse discriminant for the first live target.

16. The method of claim 15 further comprising the step of: identifying an azimuth angle of a second live target from the running average of the azimuth angles for the first monopulse discriminant obtained from azimuth angles obtained from the first live target.

17. The method of claim 15 wherein the running average of azimuth angles is determined by the equation:

18. The method of claim 15 further comprising the step of: receiving signals from a live radar target that identifies the target's elevation as being within an elevation range.

19. The method of claim 18 further comprising the steps of: assigning the running average of the first and second azimuth angles θ for the first monopulse discriminant for the first live target to a first range of elevation angles.

20. The method of claim 15 further comprising the step of: calculating a best-fit line for azimuth angles and normalized received signals.

21. In a monopulse radar receiver, a method of building a monopulse calibration table from signals received from a plurality of live radar targets, the table having N entries, each of the N entries comprising a monopulse discriminant and a comprising a corresponding azimuth angle whereat signals from the live targets are obtained, the method comprising the steps of: determining a first azimuth angle θ that corresponds to a first monopulse discriminant, from a first live target;determining a second azimuth angle θ that corresponds to the first monopulse discriminant, from a second live target;calculating a running average of the first and second azimuth angles θ for the first monopulse discriminant obtained from the first and second live targets.

22. The method of claim 21 further comprising the step of: identifying an azimuth angle of a third live target from the running average of the azimuth angles for the first monopulse discriminant obtained from azimuth angles obtained from the first and second live targets.

23. The method of claim 21 wherein the running average of azimuth angles is determined by the equation:

24. The method of claim 21 further comprising the step of: receiving signals from a live radar target that identifies the target's elevation as being within an elevation range.

25. The method of claim 21 further comprising the steps of: assigning the running average of the first and second azimuth angles θ for the first monopulse discriminant for the first live target to a first range of elevation angles.

26. The method of claim 21 further comprising the step of: calculating a best-fit line for azimuth angles and normalized received signals.

27. A radar system comprised of: a rotating antenna that receives signals at sum and difference antenna channels from targets at different elevations, within an angular-shaped beam, the center of which is the antenna's bore-sight, the angle formed between a reference direction and the bore-sight being an azimuth angle;a receiver coupled to the antenna and receiving signals from the sum and difference channels of the antenna in an elevation angle range;a computer coupled to the receiver and receiving data that represents signals received at the sum and difference channels of the antenna; andmemory, coupled to the computer, said memory storing program instructions, which when executed cause the computer to: a) calculate a normalized received signal from a target of opportunity in an elevation angle range at each of a plurality of i different azimuth angles θ, using a signal received from the target at the antenna's sum channel and using a signal received from the target at the antenna's difference channel;b) identify an azimuth angle θboresight of where the target of opportunity is actually located, as the azimuth angle whereat the magnitude of the normalized received signal obtained in step a) is greatest;c) for each of the other i azimuth angles of step a), calculate the difference between each of the i azimuth angles of step a) and the azimuth angle determined in step b) to obtain i, off-bore sight angles where each of the normalized received signals of step a) was determined;d) correlate each of the normalized received signal obtained at step a) to a corresponding off-bore sight angle determined in step c).

28. The radar system of claim 27, wherein said program instructions cause the computer to: e) determine a second target's off-bore sight angle by correlating a normalized received signal obtained from the second target, to a normalized received signal value determined for the target of opportunity and an associated off-bore sight angle determined in step d).

29. The radar system of claim 27, wherein said program instructions cause the computer to: divide the magnitude of the signal from the sum channel by the magnitude of the signal from the difference channel.

30. The radar system of claim 27, wherein said program instructions cause the computer to: divide the magnitude of the signal from the difference channel by the magnitude of the signal from the sum channel.

31. The radar system of claim 27, wherein said program instructions cause the computer to: determine a logarithmic value of the magnitude of the signal from the difference channel;determine a logarithmic value of the magnitude of the signal from the sum channel;subtract the logarithmic value of the magnitude of the signal from the difference channel from the logarithmic value of the signal from the sum channel.

32. A radar system comprised of: a rotating antenna that receives signals at sum and difference antenna channels within an angular-shaped beam, the center of which is the antenna's bore-sight, the angle formed between a reference direction and the bore-sight being an azimuth angle;a receiver coupled to the antenna and receiving signals from the sum and difference channels of the antenna;a computer coupled to the receiver and receiving data that represents signals received at the sum and difference channels of the antenna; andmemory, coupled to the computer, said memory storing program instructions, which when executed cause the computer to:a) determine at each azimuth angle θi of a plurality of i azimuth angles in an elevation angle range, a normalized received signal ΔVi for signals that are received at the antenna from a target in said elevation angle range;b) calculate a target's actual azimuth angle θboresight in the elevation angle range as the sum of the product of all θi·ΔVi, divided by the sum of all ΔVi;c) for each of the i azimuth angles θ in the elevation angle range of step a), calculate the off bore sight angle θoff bore sight of each ΔVi by calculating the angle difference between each azimuth angle θi, of step a) and the azimuth angle θboresight determined in step b) to obtain i, off-bore sight angles, θoff bore sight i, whereat each of the normalized signals ΔVi of step a) was determined;d) correlate each of the normalized received signals ΔVi obtained at step a) to a corresponding off-bore sight angle θoff bore sight i determined in step c).

33. The radar system of claim 32, wherein said program instructions cause the computer to: determine when the phase angle difference between a signal received at the antenna's sum channel and at the difference channel changes by 180 degrees.

34. The radar system of claim 32, wherein said program instructions cause the computer to calculate a normalized received signal ΔVi by: dividing the magnitude of the signal Vsum from the antenna's sum channel by the magnitude of the signal Vdiff from the antenna's difference channel.

35. The radar system of claim 32, wherein said program instructions cause the computer to: divide the magnitude of the signal from the difference channel Vdiff by the magnitude of the signal Vsum from the sum channel.

36. The radar system of claim 32, wherein said program instructions cause the computer to: calculate a logarithmic value of the magnitude of the signal Vdiff from the difference channel;calculate a logarithmic value of the magnitude of the signal Vsum from the sum channel;calculate the difference between the logarithmic value of the magnitude of the signal from the difference channel and the logarithmic value of the signal from the sum channel.

37. The radar system of claim 32, wherein said program instructions cause the computer to identify the azimuth angle whereat the magnitude of the normalized received signal obtained in step a) is greatest by: storing each of the N azimuth angles whereat a sum and difference signal ratio is determined;storing each of the N sum and difference signal ratios;sorting each of the N sum and difference signal ratios to determine the greatest ratio; andidentifying the corresponding azimuth angle where the greatest sum and difference signal ratio occurred.

38. In a monopulse radar receiver, a method of building a plurality of monopulse calibration tables from signals received from a plurality of live radar targets at a plurality of different elevation ranges, each of the plurality of monopulse calibration tables determining the off bore sight angle location of targets within a corresponding elevation range, each of the tables having N entries, each of the N entries comprising a monopulse discriminant and a comprising a corresponding azimuth angle whereat signals from the live targets are obtained, the method comprising the steps of: receiving a first signal that identifies a first live radar target as being within a first elevation range;determining a first azimuth angle θ that corresponds to a first monopulse discriminant, for the first live target;receiving a second signal that identifies a second live radar target as being within said first elevation range;determining a second azimuth angle θ that corresponds to the first monopulse discriminant, for the second live target;calculating a running average of the first and second azimuth angles θ for the first monopulse discriminant obtained from the first and second live targets within said first elevation range;assigning the running average of the first and second azimuth angles θ for the first monopulse discriminant to a first monopulse calibration table for live targets within said first elevation range.

39. The method of claim 38 further comprising the steps of: receiving a third signal that identifies a third live radar target as being within a second elevation range;determining a third azimuth angle θ that corresponds to a first monopulse discriminant, from the third live target;receiving a fourth signal that identifies a fourth live radar target as being within said second elevation range;determining a fourth azimuth angle θ that corresponds to the first monopulse discriminant, from the fourth live target;calculating a running average of the third and fourth azimuth angles θ for the first monopulse discriminant obtained from the third and fourth live targets within said second elevation range;assigning the running average of the third and fourth azimuth angles θ for the first monopulse discriminant to a second monopulse calibration table for live targets within said second elevation range.

40. The method of claim 38 wherein said first signal is an IFF Mode C response.

41. The method of claim 38 wherein said radar is a 3-D radar and said first signal identifies a target's range and horizontal distance and the target's elevation is calculated by the equation:

42. The method of claim 38, wherein each of the calibration tables corresponds to a range of target elevation angles.

43. The method of claim 42 where the number of calibration tables is five.

44. The method of claim 38 further comprising the steps of: determining an equation for a best fit line for each of the calibration tables.

45. The method of claim 44 further comprising the step of building a monopulse calibration table for live targets using the best fit line for the first monopulse calibration table and the monopulse discriminates and azimuth angles obtained from the first and second live targets.

46. The method of claim 45 further comprising the step of building a monopulse calibration table for live targets using both the best fit lines for the first and second monopulse calibration tables and the monopulse discriminates and azimuth angles obtained from the first, second, third and fourth live targets.

47. The method of claim 44 further comprising the steps of determining whether the best fit line is valid or invalid.

48. The method of claim 47 wherein the determination of a best fit line's validity is determined by calculating a standard deviation of fit line points and measured data points.