Выделить слова: 

Патент США №

8390508

Автор(ы)

Lee и др.

Дата выдачи

05 марта 2013 г.

Generating radar cross-section signatures



РЕФЕРАТ

In one aspect, a method to generate radar cross section (RCS) signatures, includes determining a spectrum of an object and using the spectrum of the object to generate RCS signatures of a plurality of objects. In another aspect, an apparatus to generate radar cross section (RCS) signatures includes circuitry to determine a spectrum of an object; and use the spectrum of the object to generate RCS signatures of a plurality of objects. In a further aspect, an article includes a machine-readable medium that stores executable instructions to generate radar cross section signatures (RCS). The executable instructions cause a machine to determine a spectrum of an object and use the spectrum of the object to generate RCS signatures of a plurality of objects.

Авторы:

Chul J. Lee (Lexington, MA), Richard A. Gilstrap (Burlington, MA)

Патентообладатель:
ИмяГородШтатСтранаТип

Chul J. Lee
Richard A. Gilstrap
Lexington
Burlington
MA
MA
US
US

Заявитель:

Raytheon Company (Waltham, MA)

ID семейства патентов

47677207

Номер заявки:

12/753,982

Дата регистрации:

05 апреля 2010 г.

Класс патентной классификации США:

342/90; 342/13; 342/175

Класс международной патентной классификации (МПК):

G01S 13/00

Класс совместной патентной классификации:

G01S 7/412 (20130101)

Область поиска:

342/13,175

Использованные источники

[Referenced By]

Патентные документы США

4123165October 1978Brown et al.
4373808February 1983Pell et al.
5086396February 1992Waruszewski, Jr.
5096281March 1992Windebank et al.
H1181May 1993Rihaczek
5227801July 1993Pierce
5305430April 1994Glassner
5317689May 1994Nack et al.
5355442October 1994Paglieroni et al.
5392050February 1995Guerci et al.
5583975December 1996Naka et al.
5588032December 1996Johnson et al.
5594844January 1997Sakai et al.
5616031April 1997Logg
5933146August 1999Wrigley
5953722September 1999Lampert et al.
6005916December 1999Johnson et al.
6031542February 2000Wittig
6212132April 2001Yamane et al.
6750805June 2004Cameron
6750859June 2004Sowizral et al.
6924763August 2005Poullin
6941303September 2005Perrizo
7123548October 2006Uzes
7289118October 2007Schmittler et al.
7348975March 2008Reshetov et al.
7369081May 2008Ganz et al.
7382397June 2008Mottur
7535408May 2009Kuhn
7567205July 2009Lee
7592947September 2009Lee
7602332October 2009Lee
7616151November 2009Lee
7646332January 2010Lee
7652620January 2010Lee
7750842July 2010Lee
2002/0060784May 2002Pack et al.
2002/0075260June 2002Brokenshire et al.
2002/0087858July 2002Oliver et al.
2003/0011519January 2003Breglia et al.
2003/0022395January 2003Olds
2003/0093430May 2003Mottur
2003/0112234June 2003Brown et al.
2005/0138073June 2005Zhou et al.
2006/0059494March 2006Wexler et al.
2006/0160495July 2006Strong
2006/0210169September 2006Colestock et al.
2007/0061619March 2007Nemirovsky et al.
2007/0165042July 2007Yagi
2008/0052334February 2008Yamazoe
2008/0132174June 2008Corral et al.
2009/0040096February 2009Lee et al.
2009/0040098February 2009Lee et al.

Зарубежные патентные документы

WO 2009/021100Feb., 2009WO
WO 2009/038894Mar., 2009WO

Другие источники


US. Appl. No. 12/212,783, filed Sep. 18, 2008, 76 pages. cited by applicant .
Notice of Allowance and Fee(s) Due dated Jul. 14, 2010, U.S. Appl. No. 12/212,783, filed Sep. 18, 2008, 8 pages. cited by applicant .
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, PCT/US2008/072434 dated Jun. 2, 2009, 4 pages. cited by applicant .
The Written Opinion of the International Searching Authority, PCT/US2008/072434 dated Jun. 2, 2009, 9 pages. cited by applicant .
Walker et al.: "Parallel Computation of Time-Domain Integral Equation Analyses of Electromagnetic Scattering and RCS", IEEE Transactions on Antennas and Propagation, IEEE Service Center, Piscataway, NJ, US, vol. 45, No. 4, Apr. 1, 1997, XP011002960, ISSN: 0018-926X. cited by applicant .
Ngoly et al.: "Parallel Post-Processing Techniques for Fast Radar Cross-Section Computation", 2006 12.sup.th Biennial IEEE Conference on Electromagnetic Field Computation, Apr. 30, 2006, Piscataway, NJ, USA, IEEE, pp. 377-377, XP010917223, ISBN: 978-1-4244-0320-2. cited by applicant .
Ling et al.: "Shooting and bouncing rays: calculating the RCS of an arbitrarily shaped cavity", IEEE Transactions on Antennas and Propagation USA, vol. 37, No. 2, Feb. 1989, pp. 194-205, XP002528225, ISSN: 0018-926X. cited by applicant .
Savides et al.: "Radar simulation using the shooting and bouncing ray technique" CCECE 2003, Canadian Conference on Electrical and Computer Engineering, Montreal, Canada, May 4-7, 2003, New York, NY, IEEE, US, vol. 1, May 4, 2003, pp. 307-310, XP010653888, ISBN: 978-0/7803-7781-3, abstract. cited by applicant .
Hu et al.: "RCS computation of jet engine with complex termination based on multiplaten Z-buffer algorithm", Microwave Conference, 1999 Asia Pacific Singapore, Nov. 30-Dec. 3, 1999, Piscataway, NJ, USA, IEEE, US, Nov. 30, 1999, pp. 781-784, XP010374299, ISBN: 978-0-7803-5761-7, abstract. cited by applicant .
Notification Concerning Transmittal of International Preliminary Report on Patentability (Chapter 1 of the Patent Cooperation Treaty), PCT/US2008/072434 dated Feb. 18, 2010, 11 pages. cited by applicant .
Notification of Transmittal of the International Search Report and the Written Opinion of the International Searching Authority, or the Declaration, PCT/US2008/072433 dated Feb. 24, 2009, 5 pages. cited by applicant .
The Written Opinion of the International Searching Authority,PCT/US2008/072433 dated Feb. 24, 2009, 8 pages. cited by applicant .
Lozano et al.: "Improvements in Ray-Tracing Acceleration Techniques to Compute Diffraction Effect and Doubles and Triples Effects in the RCS Prediction of Complex Targets" Antennas and Propagation Society Symposium, 2005. IEEE Washington, DC, Jul. 3-8, 2005, Piscataway, NJ: IEEE, US, vol. 3A, Jul. 3, 2005, pp. 93-96, XP010859931, ISBN: 978-0-7803-8883-3. cited by applicant .
Schmitz et al.:"Zpatch. A high frequency bistatic signature prediction code" Radar conference, 1997, IEEE National Syracuse, NY, USA, May 13-15, 1997, New York, NY, USA, IEEE, US, May 13, 1997, pp. 232-236, XP010224771, ISBN: 978-0-7803-3731-2. cited by applicant .
Yu et al.: "Radar cross section computation and visualization by shooting-and-bouncing ray (SBR) technique" Proceedings of the Antennas and Propagation Society International Symposium (APSIS). Chicago, Jul. 20-24, 1992; [Proceedings of the Antennas and Propagation Society International Symposium (APSIS)], New York, IEEE, US, vol. -, Jul. 18, 1992, pp. 1323-1326, XP010065787, ISBN: 978-0-7803-0730-8. cited by applicant .
Notification Concerning Transmittal of Intrnational Preliminary Report on Patentability (Chapter 1 of the Patent Cooperation Treaty), PCT/US2008/072433 dated Feb. 18, 2010, 11 pages. cited by applicant .
Ozturk: "Implementation of Physical Theory of Diffraction for Radar Cross Section Calculations", The Institute of Engineering and Sciences of Bilkent University, in Partial Fulfillment of the Requirements for the Degree of Master of Science, Jul. 2002, 72 pages. cited by applicant .
Shore et al.: "Application of incremental length diffraction coefficients to calculate the pattern effects of the rim and surface cracks of a reflector antenna", Antennas and Propagation Society International Symposium, 1993, AP-S Digest, pp. 1350-1353. cited by applicant .
Ozdemir et al.: "Fast ASAR Imag Formation Using the Shooting and Bouncing Ray Technique", IEEE Antennas and Propagation Society International Symposium, vol. 4, pp. 2605-2608, Jul. 13-18, 1997. cited by applicant .
Hansen, "Corner Diffraction Coefficients for the Quarter Plate", IEEE transactions on Antennas Propagation, vol. 39, No. 7, pp. 976-984, Jul. 1991. cited by applicant .
Shore et al.: "Application of Incremental Length Diffraction Coefficients to Calculate the Pattern Effects of the Rim and Surface Cracks of a Reflector Antenna", IEEE Antennas and Propagation Society International Symposium, 1993. cited by applicant .
Havran: "Heuristic Ray Shooting Algorithms", PhD thesis, Czech Technical University in Prague, Nov. 2000, 220 pages. cited by applicant .
Antenna Handbook, vol. I Fundamentals and Mathematical Techniques; P.H. Pathak, Ohio State University ElectroScience Laboratory; Edited by Y.T. Lo and S.W. Lee; Chapter 4, pp. 4-1-4-110. cited by applicant .
Ufimtsev, P. Ya., "Diffraction of Plane Electromagnetic Waves by a Thin Cylindrical Conductor", Radio Eng. Electron Phys., vol. 7, pp. 241-249, 1962. cited by applicant .
Pathak et al., "Model, Ray, and Beam Techniques for Analyzing the EM Scattering by Open-Ended Waveguide Cavities", IEEE Transaction on Antennas and Propagation, vol. 37, No. 5, May 1989, pp. 635-647. cited by applicant .
Pathak et al., "An Efficient Hybrid FE-BI-TW-Collective Ray Formulation for Analysis of Large Conformal Arrays", 2008 Union Radio Scientifique Internationale--(URSI) General Assembly, Aug. 7-16, 2008, 4 pages. cited by applicant .
Oguzer et al., "On the Elimination of Infinities in the PO Component of Equivalent Edge Currents", Wave Motion, vol. 18, pp. 1-10, 1993. cited by applicant .
Prabhakar H. Pathak, "High-Frequency Techniques for Antenna Analysis", Proceedings of the IEEE, vol. 80, No. 1, Jan. 1992, pp. 44-65. cited by applicant .
Book, Edited by Andrew S. Glassner, "An Introduction to Ray Tracing", U.S. Edition Published by Academic Press, 1989, 329 pages. cited by applicant .
Steve Kosanovich, "Fundamentals of Xpatch", Training Class Manual, SAIC DEMACO, Apr. 28, 2008. cited by applicant .
Real-time radar cross section of complex targets by physical optics graphical processing; Rius, J.M.; Ferrando, M.; Antennas and Propagation Society International Symposium, 1990, AP-S. Merging Technologies for the 90's'. Digest. May 7-11, 1990 pp. 1280-1283 vol. 3. cited by applicant .
Xpatch 4: the next generation in high frequency electromagnetic modeling and simulation software; Andersh, D.; Moore, J.; Kosanovich, S.; Kapp, D.; Bhalla, R.; Kipp, R.; Courtney, T.: Nolan, A.; German, F.; Cook, J.; Hughes, J.; Radar Conference, 2000. The Record of the IEEE 2000 International. cited by applicant .
3D scattering center extraction from Xpatch, Bhalla, R. Hao Ling, Dept. of Electr. & Comput. Eng., Texas Univ., Austin, TX; Antennas and Propagation Society International Symposium, 1995. AP-S. Digest, Jun. 18-23, 1995, vol. 4, on pp. 1906-1909 vol. 4, Jun. 18-23, 1995. cited by applicant .
Bhalla et al.: "Three-dimensional scattering center extraction using the shooting and bouncing ray technique", Antennas and Propagation, IEEE Transactions, vol. 44, Issue 11, Nov. 1996, pp. 1445-1453. cited by applicant .
Bhalla et al.: "Near-field signature prediction using far-field scattering centers extracted from the shooting and bouncing ray technique"; Antennas and Propagation, IEEE Transactions, vol. 48, Issue 2, Feb. 2000, pp. 337-338. cited by applicant .
U.S. Appl. No. 12/212,787, filed Sep. 18, 2008, 291 pages. cited by applicant .
U.S. Appl. No. 12/138,711, filed Jun. 13, 2008, 313 pages. cited by applicant .
U.S. Appl. No. 12/138,814, filed Jun. 13, 2008, 216 pages. cited by applicant .
U.S. Appl. No. 12/212,779, filed Sep. 18, 2008, 398 pages. cited by applicant .
U.S. Appl. No. 12/212,783, filed Sep. 18, 2008, 339 pages. cited by applicant .
U.S. Appl. No. 12/212,786, filed Sep. 18, 2008, 232 pages. cited by applicant .
U.S. Appl. No. 11/889,197, filed Aug. 9, 2007, 254 pages. cited by applicant .
U.S. Appl. No. 11/889,198, filed Aug. 9, 2007, 219 pages. cited by applicant.

Главный эксперт: Brainard; Timothy A
Уполномоченный, доверенный или фирма: Daly, Crowley, Mofford & Durkee, LLP

Интересы правительства




НАУЧНЫЕ ИССЛЕДОВАНИЯ ФИНАНСИРУЕМЫЕ ПРАВИТЕЛЬСТВОМ



This invention was made with Government support under Contract Number HQ0006-03-C-0047 awarded by the Missile Defense Agency (MDA). The United States Government has certain rights in the invention.

ФОРМУЛА ИЗОБРЕТЕНИЯ



What is claimed is:

1. A method to generate radar cross section (RCS) signatures, comprising: determining a spectrum of a first object using a computer; and using the spectrum of the first object to generate RCS signatures of a plurality of objects, using the computer, by applying the spectrum of the first object to each of the plurality of objects, using the spectrum of the first object comprising: interpolating time; selecting frequencies to match a center frequency and bandwidth of a radar pulse; resampling to frequency sampling of radar to match a range bin size; and shifting and applying trajectory phase and scaling.

2. The method of claim 1 wherein using the spectrum of the first object to generate RCS signatures of a plurality of objects comprises using a lookup table representing the spectrum of the first object.

3. The method of claim 1 wherein determining a spectrum of the first object comprises determining a frequency spectrum of the first object.

4. The method of claim 1 wherein determining a spectrum of the first object comprises determining a spectrum of an object comprising chaff.

5. The method of claim 4 wherein using the spectrum of the first object to generate RCS signatures of a plurality of objects comprises using the spectrum of the first object comprising chaff to generate RCS signatures of a plurality of objects comprising chaff.

6. The method of claim 1, further comprising: determining a spectrum of a first macro-object object using the computer; and using the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects, using the computer, by applying the spectrum of the first macro-object to each of the plurality of macro-objects, wherein using the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects comprises: interpolating time; selecting frequencies to match a center frequency and bandwidth of a radar pulse; resampling to frequency sampling of radar to match a range bin size; and shifting and applying trajectory phase and scaling.

7. The method of claim 6 wherein determining the spectrum of first object comprises determining a spectrum of a chaff puck, and wherein determining the spectrum of a first macro-object comprises determining a spectrum of a chaff dispenser.

8. An apparatus to generate radar cross section (RCS) signatures, comprising: circuitry configured to: determine a spectrum of a first object; and use the spectrum of the first object to generate RCS signatures of a plurality of objects, by applying the spectrum of the first object to each of the plurality of objects, wherein the circuitry configured to use the spectrum of the first object comprises circuitry configured to: interpolate time; select frequencies to match a center frequency and bandwidth of a radar pulse; resample to frequency sampling of radar to match a range bin size; and shift and apply trajectory phase and scaling.

9. The apparatus of claim 8 wherein the circuitry comprises at least one of a processor, a memory, programmable logic and logic gates.

10. The apparatus of claim 8 wherein the circuitry to use the spectrum of the first object to generate RCS signatures of a plurality of objects comprises circuitry to use a lookup table representing the spectrum of the first object.

11. The apparatus of claim 8 wherein the circuitry to determine a spectrum of the first object comprises circuitry to determine a frequency spectrum of the first object.

12. The apparatus of claim 8 wherein the circuitry to determine a spectrum of the first object comprises circuitry to determine a spectrum of an object comprising chaff.

13. The apparatus of claim 12 wherein the circuitry to use the spectrum of the first object to generate RCS signatures of a plurality of objects comprises circuitry to use the spectrum of the first object comprising chaff to generate RCS signatures of a plurality of objects comprising chaff.

14. The apparatus of claim 8 wherein the circuitry is further configured to: determine a spectrum of a first macro-object object using the circuitry; and use the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects, using the circuitry, by applying the spectrum of the first macro-object to each of the plurality of macro-objects, wherein the circuitry configured to use the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects comprises circuitry configured to: interpolate time; select frequencies to match a center frequency and bandwidth of a radar pulse; resample to frequency sampling of radar to match a range bin size; and shift and apply trajectory phase and scaling, wherein the circuitry configured to determine the spectrum of first object comprises circuitry configured to determine a spectrum of a chaff puck, and wherein the circuitry configured to determine the spectrum of a first macro-object comprises circuitry configured to determine a spectrum of a chaff dispenser.

15. An article, comprising: a non-transitory machine-readable medium that stores executable instructions to generate radar cross section signatures (RCS), the executable instructions causing a machine to: determine a spectrum of an object; and use the spectrum of the object to generate RCS signatures of a plurality of objects, by applying the spectrum of the first object to each of the plurality of objects, wherein the instructions causing the machine to use the spectrum of the first object comprises instructions causing the machine to: interpolate time; elect frequencies to match a center frequency and bandwidth of a radar pulse; resample to frequency sampling of radar to match a range bin size; and shift and apply trajectory phase and scaling.

16. The article of claim 15 wherein the instructions causing the machine to use the spectrum of the first object to generate RCS signatures of a plurality of objects comprises instructions causing the machine to use a lookup table representing the spectrum of the first object.

17. The article of claim 15 wherein the instructions causing the machine to determine a spectrum of the first object comprises instructions causing the machine to determine a frequency spectrum of the first object.

18. The article of claim 15 wherein the instructions causing the machine to determine a spectrum of the first object comprises instructions causing the machine to determine a spectrum of an object comprising chaff.

19. The article of claim 18 wherein the instructions causing the machine to use the spectrum of the first object to generate RCS signatures of a plurality of objects comprises instructions causing the machine to use the spectrum of the first object comprising chaff to generate RCS signatures of a plurality of objects comprising chaff.

20. The article of claim 15 further comprising instructions causing the machine to: determine a spectrum of a first macro-object object using the machine; and use the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects, using the machine, by applying the spectrum of the first macro-object to each of the plurality of macro-objects, wherein the instructions causing the machine to use the spectrum of the first macro-object to generate RCS signatures of a plurality of macro-objects comprises instructions causing the machine to: interpolate time; select frequencies to match a center frequency and bandwidth of a radar pulse; resample to frequency sampling of radar to match a range bin size; and shift and apply trajectory phase and scaling, wherein the instructions causing the machine to determine the spectrum of first object comprises instructions causing the machine to determine a spectrum of a chaff puck, and wherein the instructions causing the machine to determine the spectrum of a first macro-object comprises instructions causing the machine to determine a spectrum of a chaff dispenser.


ОПИСАНИЕ




УРОВЕНЬ ТЕХНИКИ



A radar system emits radio waves that are reflected by an object in a form of a reflected signal that is detected by the radar system. Based on the round trip time duration and the angle of the reflected signal, the location of the object may be determined. In training scenarios, instead of using actual objects, it is more practical and cost effective to use simulated radar objects. The simulated radar objects may be generated using radar signature modeling tools that emulate the radar object.

Chaff is used by aircraft to evade hostile forces or by ballistic missiles to mask the objects in the missile complex. For example, the chaff is released by an aircraft and used to create radar objects to draw enemy fire away from the aircraft. Generally, chaff includes metallic structures that are easily detected by radar. To simulate chaff in the past, extensive lookup tables were used to form a single instance of radar data of simulated chaff. These prior art approaches are time consuming, even on modern computers, so that chaff cannot be generated in real-time for radar simulation and Monte Carlo testing.

Referring to FIG. 1A, a missile 12 includes a chaff dispenser 14. The chaff dispenser releases (dispenses) chaff packages 16a, 16b, 16c (called herein "pucks") into the air. Generally, the pucks are dispensed at different times. The pucks 16a-16c disperse a cloud of chaff 18. A radar 20 detects the chaff from the pucks 16a-16c and the missile 12 so that it appears to a user there are more missiles so that more enemy munitions are directed towards the chaff than the missile.


СУЩНОСТЬ



In one aspect, a method to generate radar cross section (RCS) signatures, includes determining a spectrum of an object and using the spectrum of the object to generate RCS signatures of a plurality of objects. In another aspect, an apparatus to generate radar cross section (RCS) signatures includes circuitry to determine a spectrum of an object; and use the spectrum of the object to generate RCS signatures of a plurality of objects. In a further aspect, an article includes a machine-readable medium that stores executable instructions to generate radar cross section signatures (RCS). The executable instructions cause a machine to determine a spectrum of an object and use the spectrum of the object to generate RCS signatures of a plurality of objects.

One or more of the aspects above may include one or more of the following features. Using the spectrum of the object to generate RCS signatures of a plurality of objects may include using a lookup table representing the spectrum of the object. Determining a spectrum of an object may include determining a frequency spectrum of the object. Determining a spectrum of an object may include determining a spectrum of an object that includes chaff. Using the spectrum of the object to generate RCS signatures of a plurality of objects may include using the spectrum of the object that includes chaff to generate RCS signatures of a plurality of objects that includes chaff. Using the spectrum of the object that includes chaff to generate RCS signatures of a plurality of objects that includes chaff may include interpolating time, selecting frequencies to match a center frequency and bandwidth of a radar pulse, resampling to frequency sampling of radar to match a range bin size and shifting and applying trajectory phase and scaling.


ОПИСАНИЕ РИСУНКОВ



FIG. 1A is a view of a missile releasing chaff.

FIG. 1B is a drawing of coordinate systems of importance for this process.

FIG. 1C is a drawing of a radar antenna showing four quadrants and the different line of sight distances from them used in phase comparison monopulse.

FIG. 2 is a flowchart of an example of a process to generate radar cross section (RCS) signatures.

FIG. 3 is a flowchart of an example of a process to generate RCS signatures for multiple pucks of chaff.

FIG. 4 is a diagram depicting range determination.

FIG. 5 is an example of spectrum of one puck of chaff.

FIGS. 6A to 6C are a set of drawings of sequences of radar pulses in the time domain.

FIG. 7 is a block diagram of an example of a computer on which the process of FIG. 4 may be implemented.


ПОДРОБНОЕ ОПИСАНИЕ



Prior art radar cross section (RCS) signature prediction models are used to generate radar chaff signatures for each scenario and stored in lookup tables. A single database for multiple scenarios is not practical to generate all potential scenarios and waveforms for chaff. The alternative approach to a massive database generation for chaff is to predict RCS signature in real-time using a much smaller lookup table. The methods described herein provide a technique to generate chaff RCS signatures, for example, in real-time, while retaining a high-fidelity model.

The Acceleration Method, described herein, is a technique to quickly generate RCS signatures for the multiple puck and dispenser activities. Emulating chaff signatures for the multiple puck/dispenser activities has been achieved by combining corresponding puck signatures, extracted from single puck's signature database, after linear phase shifts according to scenario trajectories and radar beam scheduling.

While the particular examples described herein focus on chaff, the techniques described herein may be applied to other types of objects like reentry vehicles and rocket bodies. The focus is on chaff in this description because of the large numbers of pieces of chaff that can be deployed (e.g., thousands or tens of thousands deployed) whereas for other object types (e.g., reentry vehicles, rocket bodies and so forth) the numbers are typically on the order of ones to tens. To consider implementation for other object types a Final Look-up Table (LUT) of chaff is determined and a chaff cloud trajectory is generated. The Final LUT may be substituted with a target LUT of target scattering data as a function of aspect angle and chaff puck trajectory information substituted with target Center of Gravity (CG) trajectory and rotational orientation information.

Before describing the Final LUT the theory and assumptions are first described.

In signal processing one typically observes a sinusoidal signal, z(t)=Aexp(i2.pi.ft+.PHI.), (i= -1), that depends on time as the independent variable and has a fixed constant frequency, f is used. Fourier transformation provides a frequency spectrum with a peak at f. When working with Linear Frequency Modulation (LFM) waveforms in radar the convention is to use the intermediate radar pulse, prior to Fourier Transformation, as a frequency domain signal. In this convention the radar signal is represented as s(f)=Aexp(i2.pi..tau.f+.PHI.), a sinusoid that depends on frequency as the independent variable and .tau., a fixed constant for a pulse, that represents the target offset from the pulse's time reference point. Typically, the time reference is the beginning or middle of the pulse. A Fourier transform performed on this signal will have a peak at .tau. in the resulting spectrum that represents the target's time offset from the pulse's reference time. The sign convention of the signal varies from radar to radar so that one will see the use of a Fourier transform or an inverse Fourier transform depending on the radar's convention. This final Fourier transform to complete `Pulse Compression` can be a multi-step process with windowing and shifts such that any implementation of this algorithm has to match the specific radar's final Fourier transform process and sign convention.

As used herein, reference is made to `wrapping` or `ambiguous Doppler` after the Fourier transformation. This is a numerical phenomenon that occurs when working with discreet data series and the discreet Fourier transforms. To see the origin of this consider the signal exp(i.omega.n). If .omega.=.DELTA..omega.+2.pi.m for m integer, then the signal looks like exp(i.DELTA..omega.n+i2.pi.mn)=exp(i.DELTA..omega.n)exp(i2.pi.mn)=ex- p(i.DELTA..omega.n)1=exp(i.DELTA..omega.n) because exp(i2.pi.mn)=1 for all m and n.

As described herein a deramped LFM Radar Pulse Model is used. An LFM electromagnetic pulse has a frequency of oscillation that varies linearly with time (f=Bt/T+f.sub.c) so that the signal form is exp(i2.pi.(Bt.sup.2/T+f.sub.ct+.phi..sub.0)). This electromagnetic pulse reflects off the target and returns to the radar. On arrival at the radar the received pulse is mixed with the conjugate of a duplicate of the transmitted pulse that is triggered at the expected time of arrival of the reflected pulse and has used a stable oscillator as carrier frequency reference (f.sub.c). The resulting deramped signal then has the form exp(i2.pi.(B.tau..sub.rofft/T+f.sub.c.tau..sub.rt)) ignoring Doppler effects. In practice the Doppler effect on the frequency slope is pre-compensated on transmit or receive and the residual effects are negligible. The Doppler shift component due to the carrier frequency gives rise to the Range Doppler Coupling (RDC) effect and is treated as a time/range offset in later processing and may be pre-compensated. .tau..sub.rt is the round trip time for the pulse from the radar to the target and back for a coherent radar. .tau..sub.roff is a time offset between the time when the reflected pulse arrived at the radar and the reference pulse was triggered. The time offset, .tau..sub.roff would be zero if the radar could determine exactly where the target was as a function time and choose to trigger the reference pulse appropriately. For fixed sampling rate of the deramped signal t/T=n/(N-1), N=# of samples, and the deramped signal pulse for a scatterer that has not been trajectory phase corrected is: exp(i2.pi.(B.tau..sub.roffn/(N-1)+f.sub.c.tau..sub.rt)), n={-N/2 to N/2}, which is considered to be the LFM frequency domain radar signal and is the final form of the simulated pulse that is added into the radar data.

To enable coherent processing of a sequence of radar pulses it is advantageous to apply a phase correction to each pulse. .tau..sub.rt=2R(t)/c, R(t) is the target center of gravity slant range from the radar as a function of time, c is the speed of light in vacuum. In actual radar data there are tropospheric delays due to the index of refraction in the earths' atmosphere. There are other radar frequency dependent effects also. The actual term is (2R(t)+.delta.R)/c for Phase Comparison Monopulse with .delta.R depending on which quadrant the antenna phase data is processed from and where c is the speed of light. The product f.sub.c.tau..sub.rt=f.sub.c2R(t)/c=2R(t)/.lamda..sub.c typically changes by a large amount from one pulse to another given target range rates and is input to a sinusoidal function. A phase correction is applied for this trajectory motion. Typically, an estimate of target (e.g., chaff puck) center of gravity (CG) location as a function of time pulse to pulse phase variation is reduced to scatterer motion about the CG by this phase correction. Assuming RDC will be applied correctly as an additive effect later and assuming the radar reference pulse was triggered for the target CG (or Chaff puck CG), the phase corrected frequency domain LFM signal to generate the Final LUT is of the form: S.sub.nexp(i2.pi.(B.tau..sub.roffn/(N-1)+f.sub.c.tau..sub.roff)), where the phase term depends on .tau..sub.roff, the same as the frequency varying time offset. This time offset is the time offset of a scatterer/chaff piece from the CG of the target/chaff puck, the round trip time of an EM wave from the scatterer to the CG along the direction defined by the radar line of sight. It has a correspondence to the scatterer range offset from the pulse's reference range of .tau..sub.roff=2.DELTA.R.sub.roff/c. The actual term is .tau..sub.roff=(2.DELTA.R.sub.roff+.delta.R)/c for Phase Comparison Monopulse with .delta.R depending on which quadrant of antenna phase data is processed from.

FIG. 1B shows a radar located on the surface of the earth. There are many variations of coordinate systems one could employ for this problem. Of note are Cartesian coordinate systems centered in the earth, centered at the radar antenna with z axis vertical and centered at the radar antenna center with z axis perpendicular to a plane tangent to the antenna.

The other two coordinate systems of note are common radar coordinate measurement frames. RAE, Range, Azimuth, and Elevation, centered at the radar antenna, Range is the distance from the antenna center to the object, Azimuth is the angle of the object relative to the radar measured clockwise from north and elevation is the angle of the object relative to the radar measured up from the horizontal. RUV, Range, U, V is a phased array radar coordinate system, Range is the distance from the antenna center to the object, U is the direction cosine of the line of sight of the object from the radar with the transverse direction of the antenna. (x-axis of FIG. 1C) V is the direction cosine of the line of sight of the object from the radar with the vertical direction of the antenna. (y-axis of FIG. 1C).

Referring to FIG. 1C, the model described herein uses four phase centers, q.sub.1, q.sub.2, q.sub.3, q.sub.4. These phase centers can be offset from the physical antenna by appropriate relative phasing of transmit and receive processing. In particular, all four quadrants of the radar antenna transmit energy. The resulting signal arrives at the target as if it was transmitted as a point source located at the antenna center, t.sub.c. The radar energy is reflected off the target and returns to the radar antenna. Each of the four quadrants operates as its own receiver channel that is centered in the quadrant. The received signal at a quadrant then has a round trip range dependence that includes the range from the antenna center to the target, R.sub.tc, plus the range from the particular quadrant (quadrant s={1, 2, 3, 4}) to the target, R.sub.tqs=R.sub.tc+.delta.R.sub.tqs(.delta.R.sub.tqs.about..+-.d.sub.qs.d- elta..theta..sub.tq with .delta..theta..sub.tq the target angular offset from the beam center and d.sub.qs the quadrant phase center offset from the antenna center in the horizontal or vertical), so for quadrant s the received signal looks like S.sub.qs=exp(i2.pi.B/c(2.DELTA.R.sub.tc+.delta.R.sub.tqs)n/(N-1)+i2.pi.(2- R.sub.tc+.delta.R.sub.tqs)/.lamda.) The frequency dependent term, B/c(2.DELTA.R.sub.tc+.delta.R.sub.tqs)n/(N-1), is not substantially affected since .delta.R.sub.tqs is less than 1/10 of the range resolution. In the time domain this quadrant's return is then, s.sub.qs=Aexp(i2.pi.(2R.sub.tc+.delta.R.sub.tqs)/.lamda.)). The phase term is altered by 2.pi..delta.R.sub.tqs/.lamda. and is the effect of interest in capturing a measurement of the angle offset of the target from the beam center. One determines the sums and differences: .DELTA..sub..alpha.=S.sub.q1+S.sub.q2-S.sub.q3-S.sub.q4, .DELTA..sub..beta.=S.sub.q1+S.sub.q3-S.sub.q1-S.sub.q4 and .SIGMA.=S.sub.q1+S.sub.q2+S.sub.q3+S.sub.q4. After a Fourier transformation of these signals to the time domain, the ratios .DELTA..sub..alpha./.SIGMA. and .DELTA..sub..beta./.SIGMA. then give the appropriate monopulse ratios for computing the horizontal (.alpha.) and vertical (.beta.) angle offsets of the target from the beam center. One can see this easily when the difference between the received time domain signals from two phase centers (e.g., 1 and 3) is divided by the sum of the two to form the monopulse input for the target range bin, the resulting form is, i tan(2.pi..delta.R.sub.t.alpha./.lamda.)=i tan(2.pi.d.sub..alpha./.lamda..delta..theta..sub.t) (i= -1) with d.sub..alpha. the distance from the array center to the quadrant's phase center. The symbol q which denoted the quadrant has been changed to .alpha. to denote the alpha channel which corresponds to the sums and differences of the different quadrants in the horizontal, .beta. would correspond to the sums and differences of the different quadrants in the vertical.

As used herein radar signature data is a radar pulse prior to any detection processing. I and Q data is a radar pulse prior to any detection processing, where I is the real component and Q is the imaginary component. B represents radar pulse bandwidth, and is the span of wave frequencies in the electromagnetic energy processed in the receive radar pulse. f.sub.c represents a center frequency of the transmitted radar pulse. PW is a pulse width of the radar pulse (duration in time). Monopulse is a term describing the collection and processing of data and the data itself to support the measurement of angular offsets of objects in radar data. RDC is Range Doppler Coupling effect in LFM radar data and is the apparent shift in range (time) of an object due to Doppler shift of radar carrier frequency given by object range rate multiplied by PW*f.sub.c/B. R, range, slant range.

Referring back to FIG. 1C, U is the direction cosine of the line of sight of the object from the radar with the transverse direction of the antenna. (x axis.). V is the direction cosine of the line of sight of the object from the radar with the vertical direction of the antenna. (y axis.). .DELTA.r.sub.bin size=c/(2B) (N.sub.s-1)/(N.sub.fft-1). Distance spanned between one range bin and another in the time domain radar pulse. (B is pulse bandwidth, c is speed of light, N.sub.s is the number of samples of the pulse in the frequency domain and N.sub.fft is the size of the Fourier Transform applied to the pulse). .SIGMA.=Sum of m phase centers channel data, .DELTA..sub..alpha.=Appropriate sum and difference of m phase centers' data is alpha difference channel pulse and .DELTA..sub..beta.=Appropriate sum and difference of m phase centers' data is beta difference channel pulse

As described herein a two part approach is used to generate chaff: an Offline Portion and a Real-Time Portion. In the Offline Portion part, the Final LUT table is generated that includes the radar LFM response of a particular instance of a chaff cloud (e.g., a single chaff puck or a dispenser of multiple pucks) centered in the radar beam and radar range window. This is done for the full operational bandwidth of the radar, for a sequence of radar pulses at a fixed PRF (Pulse Repetition Frequency) from the chaff deployment time until it is expected that, for the scenarios of use, the chaff cloud will exit the FOV (Field of View) or move beyond the detectable range for the radar. The Final LUT generation process accesses a chaff dipole electromagnetic Scattering LUT. This Scattering LUT would typically be generated with Method of Moments electromagnetic scattering code and stores the scattering amplitude and phase for an individual chaff piece as a function of frequency and angle of incidence of the radar beam.

In the Offline Portion, radar RCS (Radar Cross Section) prediction code to generate a LUT (called herein a "Scattering LUT") of the amplitude and phase of the monostatic reflected electromagnetic energy from the radar for a single piece of chaff as a function of frequency and aspect angle is used generate the Final LUT. The aspect angle is the angle between the incident radar electromagnetic wave and the vector normal to the long axis (axis of symmetry assuming a wire dipole, this can be easily generalized to asymmetric objects) of the chaff piece. The amplitude and phase are generated for a sampling of frequencies between the minimum and maximum operating frequency of the radar, for a sampling of chaff object orientations (aspect angles) relative to the radar line of sight. The Scattering LUT can be reused for multiple chaff puck types and scenarios. A different Scattering LUT is required for each chaff type (Different chaff piece types) to be simulated.

In the second part, referred to herein as the "Real Time Portion", the Final LUT (A single instance of radar data of a chaff puck, one could generate several Final LUTs and mix and match at random to introduce additional statistical variation of the chaff) is used in a real time or near real time radar simulation as the reference data to insert the chaff data into radar pulses as the pulses are placed by the real time tracking and scheduling algorithms of the radar software. This is accomplished by appropriate interpolation, time shifting and application of trajectory phase to be described in more detail later. This approach retains a very high fidelity representation of chaff data and reduces the real time computations from being done thousands or tens of thousands of times to a few times per radar pulse and allows for Monte Carlo runs to be performed with different target trajectories and statistical variation of the chaff cloud trajectories using the radar real time tracking, discrimination and pulse scheduling algorithms.

In the Real Time portion, the Final LUT is used to generate multiple chaff pucks. In particular, using the set (pulses) of deramped LFM spectra of a single chaff puck (The Final LUT), radar pulses with chaff from multiple chaff pucks along various trajectory offsets from the baseline complex trajectory can be produced with limited computation load. This process can be applied to multiple macro-objects of the same type (see processing block 230 of FIG. 2). For example, if the macro-object is a chaff dispenser that holds chaff pucks, then the RCS signatures of multiple chaff dispensers and their chaff may be determined as a way to reduce computational load further.

This is accomplished by flying (Numerical integration of the pucks' position and velocity using a gravity model) the chaff puck trajectories to the time of the radar pulse. The appropriate pulses from the Final LUT are selected based on time from the beginning of dispensing of chaff from each puck. The difference between the radar pulse time and the puck dispensing start time is used to select pulses from just before and just after the radar pulse time. Then the appropriate frequency samples from these selected Final LUT pulses are selected for the current pulse's center frequency and bandwidth. This data is used to interpolate to the time of the radar pulse under construction and interpolated in frequency to ensure a match in radar range bin size with the radar pulse's range bin size. There is interplay between the radar pulse's pulse width (The time duration of a pulse) and the Analog to Digital converter rate that produces the frequency sampling that is the equivalent of what is in the Final LUT. The radar will generate a specific number of frequency samples for a given pulse width and then zero pad (Zero padding is the appending of zeros to a data vector until the data vector has length that is a power of two in preparation of applying a FFT.) the data up to a power of two in preparation for application of a Fast Fourier Transform (FFT) to complete pulse compression. This combination of bandwidth, number of frequency samples for the pulse width and FFT size determine the range bin size in the Final compressed pulse. The frequency data is resampled from the Final LUT such that after zero padding and FFT the range bin size is the same as the radar pulse being added to the simulated chaff. The Final LUT may have substantially fewer or substantially greater number of frequency samples than the radar pulse. (e.g., equivalent to coming from a longer or shorter pulse width pulse, but as long as B (bandwidth) and range bin size are matched the final data combination is valid). In particular, the single puck radar signatures will in general be sampled at a different frequency sampling interval than the radar pulse and have a different number of frequency samples within the pulse's bandwidth. This, combined with the use of the FFT which is restricted to a power of two for the number of points, means that the range bin size of the single puck signature will in general not match the range bin size of the radar pulse. Given a radar pulse of bandwidth B with a number of samples, N.sub.s, and a size of the FFT, N.sub.fft>N.sub.s, the Final LUT has a number of samples N.sup.L.sub.s, then a number samples, N'.sup.L.sub.s, and FFT size, N'.sup.L.sub.fft, is needed such that .DELTA.r.sup.pulse.sub.bin size=.DELTA.r'.sup.L.sub.bin size. This gives the condition that N.sub.s/N.sub.fft=N'.sup.L.sub.s/N'.sup.L.sub.fft. This is a search process for N'.sup.L.sub.s and N'.sup.L.sub.fft, (With the constraint that B is fixed and N'.sup.L.sub.fft is a power of two), because N.sup.L.sub.s can be more than two times greater or more than two times smaller than N.

This data is then time shifted (e.g., shifted in relative range within the pulse) based on the difference between the radar scheduled pulse nominal range and the particular chaff puck centroid trajectory (Remember the Final LUT had the chaff centered in the radar range window as if the radar knew exactly where to place the range window). The trajectory phase of the puck centroid is applied to the data at this point (Since each chaff piece has a phase reference to the chaff puck centroid in the LUT, this combined with the chaff puck centroid trajectory phase produces the appropriate trajectory phase for each chaff piece. Recall the pulse model (exp(i2.pi.(B.tau..sub.roffn(N-1)+f.sub.c.tau..sub.roff))) for each chaff piece in the Final LUT, this is multiplied by exp(i2.lamda.f.sub.c.tau..sub.rt). .tau..sub.rt=2R(t)/c, R(t) is the chaff puck slant range from the radar as a function of time. The time shifted pulse for each puck is then transformed to the time domain from the frequency domain, (In a process matching the radar's process for windowing and ordering of FFT and data shifts which are typically somewhat unique for each radar), and summed with the radar pulse under construction with appropriate scaling for losses. (atmospheric, beam shape etc.)

This process can be applied to produce monopulse data to provide high fidelity radar angle measurements. A full implementation for a phase comparison monopulse system would produce four Final LUT's for each of four phase center offsets (Assuming four is the appropriate number for the modeled radar) by doing the Final LUT generation process for each phase center, but with the radar reference location shifted for each to that particular monopulse channel phase center. Then the real time process is performed for each of the four monopulse channel phase center offsets and appropriate sum and difference channel data is produced for each radar pulse. An initial implementation includes the real time portion described herein, but not the four Final LUT's. The four Final LUT's for highest fidelity monopulse data is more closely tied to the trajectories that the data was created for. The approach is based on recognition that for each monopulse channel the outgoing EM wave's radar slant range, R, is considered to be the same for all channels, which is the equivalent of a single transmitter located at the mid point of the monopulse channel locations. Then for each channel the receive slant range is slightly different, R+.delta.R.sub.q. Then when the Acceleration Method range shift is computed for the q.sup.th channel, 2R+.delta.R.sub.q is used for computation of the time (range) shift of the chaff pulse and for the trajectory phase addition.

Referring to FIG. 2, one example of determining RCS signatures of objects is a process 200. The approach for generating the Final LUT to be used in the real time processing is as follows. A sequence of radar pulses of the form of a deramped LFM waveform frequency spectrum is determined for an object (Chaff Puck) (210). This LFM is the equivalent of a radar LFM waveform that has the full operational bandwidth of the radar and operational center frequency of the radar. In the time domain the chaff cloud centroid is centered in the radar pulse range window. Centering here means the chaff cloud centroid trajectory is the radar range window center reference for each simulated pulse, the final sequence of simulated pulses is exactly the aligned and phase corrected data one would expect if the radar had placed a sequence of wideband LFM pulses on the chaff puck centroid and applied a trajectory phase correction based on the puck trajectory, but it is stored in the frequency domain.

The Final LUT is made up of radar pulses in the frequency domain. The m.sup.th radar pulse's n.sup.th frequency sample, aligned and phase corrected to target CG/chaff puck CG is given by: s.sub.f(n,m)=.tau..sub.kS.sub.k(m,n)exp(i2.pi.(B.tau.(m).sub.roff.sub.--.- sub.kn/(N-1)+f.sub.c.tau.(m).sub.roff.sub.--.sub.k)) Sum over k is sum over chaff pieces' (scatterers') radar slant range offset from chaff puck CG at m.sup.th time, n is pulse frequency sample index, n={-N/2 to N/2}. S.sub.k(m,n) is the k.sup.th chaff piece scattering amplitude and phase from the Scattering LUT for chaff piece orientation at m.sup.th time at n.sup.th frequency. A 6 degree of freedom trajectory (6-DOF, CG trajectory+rotational motion of piece) is flown for each chaff piece. The radar slant range offset of the chaff piece is computed from the chaff puck trajectory and the difference in radar slant range (Line of sight) between the chaff puck and each chaff piece is computed. (.tau.(m).sub.roff.sub.--.sup.k=2.DELTA.R(t)/c) This is the round trip time offset of an electro-magnetic wave along the radar slant range of the k.sup.th chaff piece from the chaff puck for time indexed by m. The spectrum is represented in a lookup table where the rows are indexed by time since beginning of chaff dispensing for the puck and the columns are indexed by radar operating frequency. For each chaff piece in the puck a trajectory and chaff piece orientation is maintained or computed as a function of time. At each pulse time, the Scattering LUT is referenced for the current orientation of the each chaff piece to get its scattering amplitude and phase. Then this amplitude and phase is used in combination with a LFM relative range and phase offset mathematical model above, (s.sub.f(n,m)=.SIGMA..sub.kS.sub.k(m,n)exp(i2.pi.(B.tau.(m).sub.roff.sub.- --.sub.kn/(N-1)+f.sub.c.tau.(m).sub.roff.sub.--.sub.k))) to locate the chaff piece response in the radar pulse's range window with the appropriate relative phase. All of the chaff piece's LFM radar pulse responses are summed under the assumption of linear superposition in this implementation. The PRF (Pulse Repetition Frequency) of the Final LUT should be high enough that the chaff cloud expansion rate is still within the unambiguous Doppler window if one expects to use Doppler processing. The radar pulse range window should be large enough to encompass the whole chaff cloud at the end of the data interval for true fidelity also. This one LUT is then useful for multiple Monte Carlo runs of the same baseline target complex trajectory and radar location. For example, threat object and chaff puck trajectory offsets can be varied around this baseline trajectory. This limitation of same baseline target complex trajectory and radar location may be relaxed somewhat also.

Using the spectrum of one object RCS, signatures of other objects of the same type may be determined (220). For example, if the spectrum is of one chaff puck, then the spectrum of one chaff puck may be applied to multiple chaff pucks. The process of applying to cover multiple objects of the same type is applied to macro-objects (230). For example, if the macro-object is a chaff dispenser that holds chaff pucks, then the RCS signatures of multiple chaff dispensers and their chaff may be determined.

Referring to FIG. 3, an example of a process to perform processing block 220 with respect to chaff, and in particular, determining RCS values for k pucks of chaff, is a process 300. A frequency response of a dipole (e.g., cylindrical) may be generated using low frequency RCS prediction tools (e.g., using method of moments, this is the Scattering LUT). The spectrum of one puck or dispenser may be generated from the Scattering LUT. This spectrum, a sequence of radar signature data pulses, is the Final LUT used in the real time process. This final LUT is accessed to generate radar pulses with multiple pucks inserted in real time thereby accelerating the high fidelity simulation capability timeline. Thus, simulating chaff signatures for multiple puck/dispenser activities may be achieved in real time by applying linear phase shifts, (e.g., according to scenario trajectories and radar beam scheduling), to pre-computed single puck radar signatures from the Final LUT and then coherently summing them with simulated real time radar pulses.

Radar data is obtained (302). For example for a time t, simulated I & Q radar pulse data for other targets within the radar beam is obtained; Radar beam placement data, R; U; V; Rd; Radar beam parameters, bandwidth, B; center frequency, f.sub.c; pulse width, PW; scaling, S.sub.c; and the pulse's range bin size, .DELTA.r.sub.bin size.

Puck data is obtained (304). For example, for each puck a trajectory and dispense time information is obtained. For example, for a k-th puck, the initial puck dispense time for the k-th puck, t.sub.0,k, is obtained and a slant range, R.sub.k is obtained. In one example, k.sup.th chaff puck trajectory is propagated to current radar pulse time. The slant range R.sub.k, from the radar antenna center to the chaff puck is determined. The slant range, R.sub.p, to the radar pulse's nominal reference point as placed by the radar pulse scheduler is obtained. Then the slant range offset .DELTA.R.sub.k=R.sub.k-R.sub.p, of the chaff puck relative to the radar pulse reference is determined.

Referring to FIG. 4, in one example, R.sub.k is measured from an origin 402 of a radar (Antenna center) to the k-th puck 404 along a direction of the radar beam 408. A slant range difference, .DELTA.R.sub.k, is the difference of the slant range, R.sub.k, of the chaff puck and the reference range of a radar pulse beam 408.

Referring back to FIG. 3, for each puck, processing blocks 310-360 are performed. A puck pulse is selected for time, t (310). For example, a frequency domain simulated pulse is obtained from a look-up table, the Final LUT, representing a spectrum of a puck shown in FIG. 5. The frequency domain simulated pulse to be added into the radar pulse is interpolated to the current time offset of the radar pulse by using the Final LUT pulses from just before and just after the current pulse time (320). For example, an interpolated time, .DELTA.t.sub.k, is equal to time, t, minus the dispense time of the k-th puck, t.sub.0,k. The dispense time corresponds to the time the chaff is first dispensed from the puck.

The frequencies are selected (330). For example, frequency bins are selected to match the current radar pulse's center frequency, f.sub.c, and bandwidth, B. For example, f.sub.c-B/2 to f.sub.c+B/2. A resampling is performed to frequency sampling of the radar to match range bin size, .DELTA.r.sub.bin size (340).

A range shift, trajectory phase and scaling for each radar phase center is applied (350). (Here the chaff puck data is shifted from the center of the pulse, as stored in the Final LUT, to the offset due to the radar's pulse placement, adding in the trajectory phase component of the chaff puck to the pulse from the Final LUT data and then scaling for beam shape, atmospheric and other losses and unit conversion from meters for chaff signal amplitude as stored in the Final LUT to the radar internal signal processing units.)

The radar pulse round trip time, .tau..sub.rt=(R.sub.k+R.sub.mrtk)/c and object time offset from the radar pulse's reference time, .tau..sub.roff=(.DELTA.R.sub.k+.DELTA.R.sub.mk)/c. These are the outbound slant range from the radar antenna center, R.sub.k, and the returning pulse slant range from the m.sup.th antenna phase center to the object R.sub.mk. The radar LFM signal model is exp(i2.pi.(B.tau..sub.roffn/(N.sub.s-1)+f.sub.c.tau..sub.rt)). For processing convenience the radar pulse's range bin size .DELTA.r.sub.bin size and c=.lamda..sub.cf.sub.c are used to change the dependence of the radar LFM signal model from .tau..sub.rt and .tau..sub.roff to .DELTA.R.sub.fmk and (R.sub.k+R.sub.mrtk). For example the fractional portion and integer portion, .DELTA.R.sub.fmk and .DELTA.R.sub.ik, of (.DELTA.R.sub.k+.DELTA.R.sub.mk)/2 modulo .DELTA.r.sub.bin size are computed where |.DELTA.R.sub.fmk|<.DELTA.r.sub.bin size and .DELTA.R.sub.ik=M .DELTA.r.sub.bin size, .DELTA.r.sub.bin size is a range bin size for the radar waveform and M is an integer such that (.DELTA.R.sub.k+.DELTA.R.sub.mk)/2=.DELTA.R.sub.fmk+.DELTA.R.sub.i,k. This split into fractional and integer portions is done so that the data is not wrapped around in the time domain after applying the Discreet Fourier Transform. The fractional component is required for narrow bandwidth pulses which can have range bin sizes of greater than 20 meters. It is at this point that the appropriate RDC adjustments are made, if necessary, based on the radar's handling of it. One should compute one value of .DELTA.R.sub.ik to be used for all four phase centers since the integer shift is most efficiently applied after the phase center data is combined and transformed to the time domain. The equation for the fractionally shifted puck data is then; S.sub.cw.sub.nS'.sub.nexp(i2.pi..DELTA.R.sub.fmk/.DELTA.r.sub.bin size n/(N.sub.s-1)+i2.pi.(R.sub.k+R.sub.mrtk).lamda..sub.c), with S'.sub.n the n.sup.th frequency sample of the resampled (frequency) and interpolated (time) Final LUT data selected for the k.sup.th puck and the current pulse. S.sub.c is the scaling for losses and unit conversion and w.sub.n the window, applied in preparation for Fourier transformation, for sidelobe reduction in the time domain. 2.lamda.(R.sub.k+R.sub.mrtk)/.lamda..sub.c is the trajectory phase term for the chaff puck.

A summation is performed (360). For example, sums and differences are performed to form .SIGMA., .DELTA..sub..alpha., .DELTA..sub..beta. data, this is the sum signal and monopulse information, so angular information can be measured on targets.

In this implementation these sums and differences are done in the frequency domain prior to Fourier transformation to reduce processing load. Then, for the current puck's data, they're transformed to the time domain, and the integer shift is applied according to .DELTA.R.sub.ik/.DELTA.r.sub.bin size. This shifted time domain puck data is added to the radar pulse.

It is determined if another puck is available (370). If another puck is available, processing blocks 310 to 370 are performed again.

Referring to FIGS. 6A to 6C, show sequences of radar pulses in the time domain. FIG. 6A shows aligned chaff cloud data as stored in the Final LUT. FIG. 6B shows the chaff cloud data shifted based on a radar's reference pulse alignment as done in the Real Time Process. FIG. 6C shows the chaff cloud with wrapping at the end of the data if the shift done in the Real Time Process were not broken into integer and fractional parts.

The following pseudo code may be used to determine RCS for chaff. References in the pseudo code are made to process 300.

TABLE-US-00001 Algorithm Pseudo Code Offline Construct Scattering LUT Final LUT Real Time/Near Real Time Portion (Example implementation of algorithm) While Radar Makes Pulse Requests (302) For each chaff puck (Index k) Propagate k.sup.th chaff puck trajectory to pulse time (304) Compute R.sub.k, the slant range from the radar antenna center to the chaff puck Compute .DELTA.R.sub.k = R.sub.k - R.sub.p, R.sub.p is the slant range to the pulse's nominal reference point as placed by the radar pulse scheduler. Interpolate to pulse time and frequency bins from Final LUT, S'.sub.n (310, 320, 330, 340) (Frequency Index n) (Could be M.sub.phase_.sub.center of Final LUT's and S'.sub.mn 's if highest fidelity desired) (End 310, 320, 330, 340) Z.sup.k.sub.mn is m.sup.th phase center frequency domain chaff puck radar pulse data For each radar antenna phase center (Index m) (350) Compute .DELTA.R.sub.mk and R.sub.mrtk, k is index of current puck, R.sub.mrtk is the slant range from the m.sup.th phase center to the chaff puck, and .DELTA.R.sub.mk = R.sub.mrtk - R.sub.p Compute .DELTA.R.sub.mfk and .DELTA.R.sub.ik, The fractional portion and integer portion of (.DELTA.R.sub.k + .DELTA.R.sub.mk)/2 modulo .DELTA.r.sub.bin size .tau..sub.rt = (R.sub.k + R.sub.mrtk)/c and .tau..sub.roff = (.DELTA.R.sub.k + .DELTA.R.sub.mk)/c, .DELTA.r.sub.bin size = c/(2B) (N.sub.s - 1)/(N.sub.s - 1) One should compute one value of .DELTA.R.sub.ik to be used for all four phase centers since the integer shift is most efficiently applied after the phase center data is combined. z.sup.k.sub.mn = S.sub.cS'.sub.nexp(i2.pi. .DELTA.R.sub.mfk/.DELTA.r.sub.bin size n/(N.sub.s - 1) + i2.pi.(R.sub.k + R.sub.mrtk)/.lamda.) n = {-N.sub.s/2 to N.sub.s/2} to minimize additional phase variation that can bias Doppler processing results S.sub.c is scaling for beam shape, atmospheric and other losses and scaling. (i = - 1) End for each radar phase center (End 350) Frequency domain fractional shift pulse data (360) f.DELTA..sup.k.sub..alpha.n = z.sup.k.sub.1n + z.sup.k.sub.2n - z.sup.k.sub.3n - z.sup.k.sub.4n, Alpha channel for k.sup.th puck f.DELTA..sup.k.sub..beta.n = z.sup.k.sub.1n + z.sup.k.sub.3n - z.sup.k.sub.2n - z.sup.k.sub.4n, Beta channel for k.sup.th puck f.SIGMA..sup.k.sub.n = z.sup.k.sub.1n + z.sup.k.sub.2n + z.sup.k.sub.3n + z.sup.k.sub.4n, Sum channel for k.sup.th puck, (End 360) Transform f.DELTA..sup.k.sub..alpha.n, f.DELTA..sup.k.sub..beta.n, f.pi..sup.k to time domain .DELTA..sup.k.sub..alpha.t = transform(f.DELTA..sup.k.sub..alpha.n) .DELTA..sup.k.sub..beta.t = transform(f.DELTA..sup.k.sub..alpha.n) .pi..sup.k.sub.t = transform(f.pi..sup.k.sub.n) Perform integer shift corresponding to .DELTA.R.sub.ik/.DELTA.r.sub.bin size .DELTA..sup.k.sub..alpha.s = integer shift of .DELTA..sup.k.sub..alpha.t .DELTA..sup.k.sub..beta.s = integer shift of .DELTA..sup.k.sub..beta.t .pi..sup.k.sub.s = integer shift of .pi..sup.k.sub.t Sum with current radar pulse under construction (380) Z.sub..DELTA..alpha. = Z.sub..DELTA..alpha. + .DELTA..sup.k.sub..alpha.s Z.sub..DELTA..beta. = Z.sub..DELTA..beta. + .DELTA..sup.k.sub..beta.s Z.sub..pi. = Z.sub..pi. + .pi..sup.k.sub.s (End 380) End for each chaff puck Pulse data is now constructed and ready for detection processing .pi. = Sum of m phase centers' Z.sub..pi. is sum channel pulse .DELTA..sub..alpha. = Appropriate sum and difference of m phase centers' Z.sub..DELTA..alpha. is alpha difference channel pulse .DELTA..sub..beta. = Appropriate sum and difference of m phase centers' Z.sub..DELTA..beta. is beta difference channel pulse End While Radar Makes Pulse Requests (End 302)

Referring to FIG. 7, a computer 700 includes a processor 722 and a volatile memory 724, a non-volatile memory 726 (e.g., a hard disk). The non-volatile memory 726 stores computer instructions 734, an operating system 736 and data 738 including a one puck spectrum lookup table 750 (e.g., Final Look-Up Table), radar data 760 and puck data 770. In one example, the radar data includes I & Q, time, R, U, V, B, f.sub.c, PW Scaling, and bin size. In one example, the puck data 770 includes puck dispense time and puck trajectory. In one example, the computer instructions 734 are executed by the processor 722 out of volatile memory 724 to perform all or part of the process 300.

Process 300 is not limited to use with the hardware and software of FIG. 7; they may find applicability in any computing or processing environment and with any type of machine or set of machines that is capable of running a computer program. Process 300 may be implemented in hardware, software, or a combination of the two. Process 300 may be implemented in computer programs executed on programmable computers/machines that each includes a processor, a storage medium or other article of manufacture that is readable by the processor (including volatile and non-volatile memory and/or storage elements), at least one input device, and one or more output devices. Program code may be applied to data entered using an input device to perform process 300 and to generate output information.

The system may be implemented, at least in part, via a computer program product, (e.g., in a machine-readable storage device), for execution by, or to control the operation of, data processing apparatus (e.g., a programmable processor, a computer, or multiple computers)). Each such program may be implemented in a high level procedural or object-oriented programming language to communicate with a computer system. However, the programs may be implemented in assembly or machine language. The language may be a compiled or an interpreted language and it may be deployed in any form, including as a stand-alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program may be deployed to be executed on one computer or on multiple computers at one site or distributed across multiple sites and interconnected by a communication network. A computer program may be stored on a storage medium or device (e.g., CD-ROM, hard disk, or magnetic diskette) that is readable by a general or special purpose programmable computer for configuring and operating the computer when the storage medium or device is read by the computer to perform process 300. Process 300 may also be implemented as a machine-readable storage medium, configured with a computer program, where upon execution, instructions in the computer program cause the computer to operate in accordance with process 300.

The processes described herein are not limited to the specific embodiments described. For example, the process 300 is not limited to the specific processing order of FIG. 3. Rather, any of the processing blocks of FIG. 3 may be re-ordered, combined or removed, performed in parallel or in serial, as necessary, to achieve the results set forth above.

The processing blocks in FIG. 3 associated with implementing the system may be performed by one or more programmable processors executing one or more computer programs to perform the functions of the system. All or part of the system may be implemented as, special purpose logic circuitry (e.g., an FPGA (field programmable gate array) and/or an ASIC (application-specific integrated circuit)).

Elements of different embodiments described herein may be combined to form other embodiments not specifically set forth above. Other embodiments not specifically described herein are also within the scope of the following claims.

* * * * *