Search Images Maps Play YouTube News Gmail Drive More »
Sign in
Screen reader users: click this link for accessible mode. Accessible mode has the same essential features but works better with your reader.

Patents

  1. Advanced Patent Search
Publication numberUS4713666 A
Publication typeGrant
Application numberUS 07/017,627
Publication dateDec 15, 1987
Filing dateFeb 24, 1987
Priority dateFeb 28, 1986
Fee statusLapsed
Also published asDE3764121D1, EP0237404A1, EP0237404B1
Publication number017627, 07017627, US 4713666 A, US 4713666A, US-A-4713666, US4713666 A, US4713666A
InventorsJean P. Poux
Original AssigneeThomson-Csf
Export CitationBiBTeX, EndNote, RefMan
External Links: USPTO, USPTO Assignment, Espacenet
Method of processing the sum and difference signals of a radar of the monopulse type for estimating the parasite phase introduced between these signals by the ultrahigh frequency formation circuits of the sum and difference channels
US 4713666 A
Abstract
The method of the invention, applied for example to estimating the parasite phase φa between the azimuth sum Σa and difference Δa channels, consists in calculating the expressions S=Σa Σe * and D=Δa Σe * for different successive measurements made by the radar, then in calculating the expressions [S.sub.(k+1) -S.sub.(k-1) ]D*(k) for different successive times k-1, k, k+1 and in averaging the result of this expression or different times k, it being understood that the radar antenna sweeps in azimuth and that the radar wave emitted is polarized circularly, the single pulse azimuth and elevation receivers receiving respectively only one of the circular polarization types: right hand or left hand.
Images(2)
Previous page
Next page
Claims(4)
What is claimed is:
1. A method of processing the sum Σ and difference Δ signals, in azimuth: Σa, Δa and in elevation: Σe, Δe, of a radar of the monopulse type, for estimating the parasite phase Ψ introduced between these signals by the ultrahigh frequency formation circuits of the sum and difference channels, comprising:
in calculating for each measurement made by the radar at least one case for the expressions S and D:
first case: ##EQU14## for estimating the azimuth parasite phase Ψa ; second case: ##EQU15## for estimating the elevational parasite phase Ψe ; third case: ##EQU16## for estimating the differential parasite phase δΨ=Ψae, respectively,
(where the symbol * designates the complex conjugate quantity), it being understood that the radar wave emitted is polarized circularly and that the azimuth
and elevation monopulse receivers receive respectively only one of the circular polarization types: right hand or left hand;
in calculating the expression
[S(k+1)-S(k-1)]D*(k)
it being understood that the radar antenna sweeps during time in azimuth in the first case and in the second case, in elevation in the second case and that k-1, k, k+1 designate three successive times such that the difference S(k+1) (θ)-s.sub.(k-1) (θ) is of the same sign as d.sub.(k) (θ), where S.sub.(k+1) (θ) and s.sub.(k-1) (θ) respectively designate the gain of the sum channel at angle θ, in azimuth in the first case and in the third case, in elevation in the second case, respectively at times k+1 and k-1, and where d.sub.(k) (θ) designates the gain of the difference channel at angle θ and at time k, in the azimuth in the first and third cases, in elevation in the second case;
in calculating the mean value of the result of the preceding expression for several successive times corresponding to several measurements successively made by the radar, which allows an expression of the form ejΨ. a where Ψ is the parasite phase sought and "a" a real number of given absolute value and sign.
2. Method as claimed in claim 1, wherein the differential parasite phase δΨ is estimated in a calibration phase prior to the phase of using the radar, and using an antenna sweep, and wherein the azimuth parasite phase Ψa is then deduced, during the phase of using the radar, from the differential parasite phase δΨ thus estimated during the calibration phase, and from the elevational parasite phase Ψe estimated during the phase of using the radar, without it being necessary to use antenna sweeping, by averaging the expression Σa Δa over several successive measurements made by the radar.
3. The method as claimed in claim 2, wherein, since only the product [S(k+1)-S(k-1)]D*(k) being of importance the expressions S and D are calculated by means of a single device used successively for calculating S and D.
4. The method as claimed in claim 1, wherein, since only the product [S(k+1)-S(k-1)]D*(k) being of importance the expressions S and D are calculated by means of a single device used successively for calculating S and D.
Description
BACKGROUND OF THE INVENTION

1. Field of the Invention

The present invention relates generally to the field of radar and concerns more particularly a method of processing the signals of a radar of the monopulse type, for estimating a parasite phase introduced between these signals by the ultrahigh frequency formation circuits of the sum and difference channels.

2. Description of the Prior Art

It is known that with a radar of the monopulse type the angular location of targets may be obtained by processing the different signals received simultaneously from these targets and corresponding to different beam directivities.

Thus, as shown very schematically in FIG. 1, by providing two sources respectively in the azimuth plane (1a, 2a) and in the elevational plane (1e, 2e), at reception there are available for each of these two planes two signals of which the sum Σ and the difference Δ may be formed by ultrahigh frequency means (3, 4), which is tantamount to having an antenna with a beam corresponding to the Σ channel and another beam corresponding to the Δ channel, and thus allows angular location of the targets to be obtained simply.

It is also known that, with the reradiation of a target being a sum of the waves reflected by each of the reflecting elements which form it, the signals S and D received respectively in the sum and difference channels, in azimuth or in elevation, have as complex representation: ##EQU1## with: s(θi): gain of the sum channel at the angle θi ;

d(θi): gain of the difference channel at the angle θi ;

ai : amplitude relative to the ith reflector;

φi : phase shift introduced by the ith reflector;

di : distance between the radar and the ith reflector;

Ψ: parasite phase introduced between the Σ and Δ channels by the ultrahigh frequency formation circuits of the sum and difference channels.

It is also known that, in the particular case of the elevational plane, the parasite phase Ψ between the sum and difference channels may be evaluated in the following way:

let ##EQU2## where D* designates the complex conjugate quantity of D, and where: ##EQU3##

If we calculate the mean value of P over different distances divisions of the radar, the phase ##EQU4## being considered as random, we have: ##EQU5## where E designates the expectation operator.

The term "a" is real, but is a priori of unknown sign because of the factors d(θi) which may be positive or negative. In fact, referring to FIG. 2 which shows the trend of the sum and difference diagrams, Σ and Δ, in the azimuth plane or in the elevational plane, the sum diagram has the maximum in the direction of the axis of the antenna, whereas the difference Δ diagram has on the contrary a minimum in this direction and is, on each side of this direction, either in phase (Δ+) or in phase opposition (Δ-) with Σ.

By averaging P, we can then a priori only estimate the phase Ψ to with in the sign.

Now, it so happens that in the elevational plane this ambiguity of sign may be relatively simply removed for the sign of d(θi) is always related to the rank of the range bin considered being negative for the closest range bins and positive for the furthest range bins. By averaging P over the different range bins, an estimation of the phase Ψ is obtained without ambiguity of sign after correction of the sign as a function of the considered range bin.

SUMMARY OF THE INVENTION

The present invention relates to a method of estimating the parasite phase Ψ for resolving these ambiguities of sign also in the azimuth plane, and whose principle is also applicable to the estimation of the parasite phase Ψ in the elevational plane, as well as to the estimation of the differential phase δφ between the azimuth and elevational planes, this latter possibility moreover forming, as will be seen further on, the basis of another method of estimating the parasite phase in the azimuth plane.

According to the invention, the method of processing the sum Σ and difference Δ signals, in azimuth: Σa, Δa, and in elevation: Σe, Δe, of a radar of the monopulse type, for estimating the parasite phase Ψ introduced between these signals by the ultrahigh frequency formation circuits of the sum and difference channels consists:

in calculating for each measurement made by the radar the expressions S and D, with ##EQU6## (where the symbol * designates the complex conjugate quantity), depending on whether it is desired to estimate the parasite phase Ψa in azimuth in the first case, or the parasite phase Ψe in elevation in the second case or the differential parasite phase δΨ=Ψae in the third case, it being understood that the radar wave emitted is polarized circularly and that the azimuth and elevation monopulse receivers only receive respectively one of the circular polarization types: right or left,

in calculating the expression

[S(k+1)-S(k-1)]D*(k)

it being understood that the radar antenna sweeps during time in azimuth in the first case and in the third case, in elevation in the second case and that k-1, k, k+1 designate three successive times such that the difference s.sub.(k+ 1) (θ)-s.sub.(k-1) (θ) is of the same sign as d.sub.(k) (θ), where s.sub.(k+1) (θ) and s.sub.(k-1) (θ) designate respectively the gain of the sum channel at the angle θ, in azimuth in the first case and in the third case, in elevation in the second case, respectively at times k+1 and k-1, and where d.sub.(k) (θ) designates the gain of the difference channel at the angle θ at time k, in azimuth in the first and in the third cases, in elevation in the second case;

in calculating the mean value of the result of the preceding expression for several successive times corresponding to several measurements successively made by the radar, which allows an expression to be obtained of the form ejΨ.a, where Ψ is the parasite phase sought and "a" a real number of given absolute value and sign.

BRIEF DESCRIPTION OF THE DRAWINGS

Other objects and features of the present invention will appear more clearly from the following description of embodiments made with reference to the accompanying drawings in which:

FIG. 1 shows schematically the formation of the sum and difference channels in azimuth and in elevation in a radar of monopulse type;

FIG. 2 shows the form of the sum and difference diagrams of a monopulse radar in the azimuth plane or in the elevational plane;

FIG. 3 illustrates the sweeping of the radar antenna, according to the invention;

FIG. 4 is a block diagram of the processing carried out in accordance with the invention; and

FIG. 5 is a diagram of the product detector used in the processing according to the invention.

DESCRIPTION OF THE PREFERRED EMBODIMENT

The invention is used in a radar system with circularly polarized coherent or non coherent transmission, capable of receiving and discriminating the right hand circular polarization and left hand circular polarization. For the sake of clarity in what follows, the azimuth channels will receive the right hand polarization whereas the elevational channels receive the left hand polarization; the reverse would also be possible.

Thus, referring to an absolute hypothetical phase the sum signals Σa and Σe, respectively in azimuth and in elevation, and the difference signals Δa and Δe, respectively in azimuth and in elevation, are written: ##EQU7## where the indices D and G designate the channels receiving the right hand and left hand polarizations, and the indices a and e the angular azimuth and elevational dimensions, and where φt designates the non coherent phase component which changes from pulse to pulse.

From the sum and difference signals Σa, Σe, Δa, Δe, the following expressions are calculated.

S=Σa Σe *

D=Δa Σe *

With a unit of time fixed sufficiently small for the phase shift between channels to be considered as constant, we have at a time k: ##EQU8##

It will be noted that calculation of the expressions S and D allows, among other things, to be free of the non coherent phase component φt.

We then calculate the expression:

[S(k+1)-S(k-1)]D*(k)

where k-1, k and k+1 designate three successive times.

We have: ##EQU9##

In the expression ##EQU10## the factors bn are assigned with a phase term φn which, since it may be considered as random, will disappear by averaging the expression: [S(k+1)-S(k-1)]D*(k) over several successive measurements of Σa, Δa, Σe, Δe made by the radar and corresponding to several successive times k.

The notion of measurements made by the radar relates to the notion of range bins and pulses. For a given range bin, the averaging is carried out over several successive pulses, and averaging is then carried out over several distance divisions. Because of the double polarization. the factors to which such a phase term has not been assigned are limited to the expression "a" for which we have the equality i=l and j=m. All the other cases (such as i=j or m=l) are in fact excluded for they would correspond to individual reflectors which would reflect the two types of polarization: right hand and left hand, that is to say which would be both of the dihedron type (even number of successive reflections) and of the trihedral type (uneven number of successive reflections).

The invention allows then any angular information between the signals received to be kept, while eliminating the non coherence of the transmission, by referencing them to the response to the same pulse of different reflectors. This is achieved through the combination of the calculation of expressions S and D and the double polarization.

During the preceding development, the important assumption was made that between times k-1 and k+1 the angles of elevation θe (in this case θG since we are considering in this example the case where the azimuth channels receive the right hand polarization whereas the elevational channels receive the left hand polarization) are invariable, which allows the common factorization in the expression "a". This assumption is reasonable in the case of a fixed radar or of an airborne radar if the time unit is sufficiently small so that advance of the platform does not change the geometry of the system.

The invention is moreover used in an azimuth sweep radar system, the rotational speed of the antenna (in azimuth) being fixed at a value such that the lobes Δ- and Δ+ at time k coincide respectively with the lobe Σ at time k-1 and with the lobe Σ at time k+1, as shown in FIG. 3.

Under these conditions, the product

[sa(k+1) (θDi)-sa(k-1) (θDi)]da(k) (θDi)

is positive for any elementary reflector "i". In fact, if we consider for example an elementary reflector which is at time k on lobe Δ+, we then have da(k) which is positive, sa(k+1) which is also positive and sa(k-1) which is very small. The product is therefore positive; it would be the same for an elementary reflector which at time k would be on lobe Δ-.

The azimuth speed of rotation of the antenna is not critical and could be chosen differently. The proposed choice seems however optimum in that it gives the largest values of the difference (sa(k+1) -Sa(k-1).

The factor "a" is therefore written as the sum of a large number of positive definites. On the contrary the factors bn are assigned with a phase term φn which may be considered as random.

Examination of several successive measurements of [S(k+1)-S(k-1)]D*(k), by averaging, allows the phase correction Ψ13 to be estimated whose variance will vary inversely with the number of measurements made, this number of measurements being fixed a priori or adapted to the desired accurancy by software.

Adaptation of the repetition period of the measurements made by the radar (coresponding to the notion of range bins and pulses) and of the repetition period of times k, that is to say the sweeping speed of the antenna, may be made in different ways, for example by suitably choosing the repetition period of the pulses or else by carrying out post integration over several pulses.

In FIG. 5 a diagram has been shown of the product detector used in the processing carried out in accordance with the invention for estimating the parasite phase between the sum and difference channels in azimuth.

The signals Σa, Δa and Σe are applied to the input of this product detector which calculates the expressions:

S=Σa Σe *

D=Δa Δe *

by means of two mixers 5 and 6 one of which receives on the one hand the signal Σe and on the other hand either the signal Σa, or the signal Δa, and the other of which receives on the one hand the signal Σa phase shifted by 90 by means of a phase shifter 7 and, on the other hand, either the signal Σa or the signal Δa. The signals Σa Δa and Σe are further transformed into medium frequency signals by means of frequency change stages 8 and 9.

It should be noted that S and D are not used at the same times, since only the product [S(k+1)-S(k-1)]D*(k) is important. It is then possible to use only a single product detector, addressed successively through a switch 10 by channels Σa and Σc, then Δa and Σc.

The magnitudes S and D are complex magnitudes which may be written in the form I+jQ where I and Q designate respectively their real part and their imaginary part; magnitudes I and Q are obtained at two outputs of the product detector.

The other steps of the method may be achieved using adapted software.

In the foregoing, the invention has been described for estimating the parasite phase between the sum and difference channels in azimuth. The principle of the invention is nevertheless applicable to estimating the parasite phase between the sum and difference channels in elevation.

In this case, the following expressions would be calculated:

S=Σe Σa *

D=Δe Σa *

The constraint concerning the sweep would then become a constraint concerning the elevational sweep, and the assumption that between times k-1, k and k+1 the azimuth angles are invariable would also have to be made.

As will now be seen the principle of the invention is also applicable to estimating the differential parasite phase Ψae, where Ψa and Ψe designate respectively the parasite phase in azimuth and in elevation. Estimation of this differential parasite phase may be interesting for the following reasons.

In the case more especially of a radar carried by a missile whose supersonic speed decreases during flight to subsonic values, the high temperature differences between the beginning and the end of the mission cause a variation of the azimuth and elevation parasite phases.

Now, the above mentioned antenna sweep constraints mean that estimation of a phase correcting factor is not possible during the whole period of use and processing of the radar.

Starting from the fact that the four receivers (two in azimuth and two in elevation) are physically close to each other and thermally interdependent, so that the variations of Ψa and Ψe will be similar, it is then sufficient to evaluate Ψa and Ψe during a so called calibration phase (previous to any operation in the search or tracking mode) and to calculate the difference Ψae which will remain constant during the whole flight time. In the tracking phase, it is still possible, as recalled in the introduction, to calculate Ψe by a method of evaluation which offers no antenna position or frequency constraint. It is then possible to deduce Ψa therefrom.

However, a successive evaluation of Ψa and Ψe during the calibration phase risks inducing an error in Ψae as for as the temperature may have changed between the two evaluations. Moreover, the calibration phase lasts a considerable time with respect to the total time of the mission of a short range missile.

The method of estimating the parasite phase of the invention allows the differential phase Ψae to be calculated directly without passing through the successive calculations of Ψa and Ψe, which removes the above mentioned risk of error.

This method will now be described applied to estimating the differential phase Ψae. As before, the signals received by the four angle error measurement channels are written: ##EQU11## where the indices D and G designate the channels receiving the right hand and left hand polarizations (still assuming that the azimuth channels receive the right hand polarization and the elevational channels the left hand polarization) and indices a and e the angular azimuth and elevation dimensions.

ai is the amplitude of the ith individual reflector;

φi is the phase shift introduced by the ith reflector;

di is the radar-ith reflector distance;

φt is the non coherent phase component which changes from pulse to pulse;

Ψ1, Ψ2, Ψ3, Ψ4 are the parasite phases introduced in the four reception channels.

According to the same method as before, a product detector calculates:

S=Σa Σe *

D=Δa Δe *

Then ##EQU12## where the phases φn may be considered as random phases and where: ##EQU13##

Now, the azimuth rotational speed of the radar has been chosen so that the product: [sa(k+1) (θDi)-sa(k-1) (θDi)]da(k) (θDi) is positive for all the elementary reflectors.

Moreover, seGj) is positive, a2 Di a2 Gj is positive and deGj) is negative for the nearest range bins and is positive for the furthest range bins

The ambiguity of sign may therefore be removed and corrected so that the factor "a" is the sum of a large number of positive real numbers

Examination of several successive measurements of [S(k+1)-S(k-1)]D*(k) allows by averaging to estimate the phase correction Ψ1234 that is to say Ψae.

The average is calculated for several range bins and several pulses so as to obtain the desired accuracy.

Patent Citations
Cited PatentFiling datePublication dateApplicantTitle
US3794998 *Apr 26, 1972Feb 26, 1974Raytheon CoMonopulse radar receiver with error correction
US4368468 *Dec 22, 1980Jan 11, 1983Westinghouse Electric Corp.Monopulse radio receiver compensation apparatus
US4568940 *Dec 22, 1981Feb 4, 1986Hughes Aircraft CompanyDual-mode radar receiver
US4646095 *Aug 16, 1985Feb 24, 1987Raytheon CompanyMethod of resolving closely spaced targets
US4656480 *Nov 15, 1984Apr 7, 1987U.S. Philips CorporationMonopulse FM-CW radar system with improved axis stability
EP0038734A1 *Apr 3, 1981Oct 28, 1981Thomson-CsfApparatus for normalizing the gradient of the off-boresight measurement by radar and air-earth radar comprising such an apparatus
Referenced by
Citing PatentFiling datePublication dateApplicantTitle
US4968969 *Feb 12, 1990Nov 6, 1990Thomson-CsfDevice and method for data transmission and/or acquisition using two cross polarizations of an electromagnetic wave and magnetic recording device
US5808578 *Dec 20, 1996Sep 15, 1998Barbella; Peter F.Guided missile calibration method
US5986605 *Apr 6, 1998Nov 16, 1999Raytheon CompanyMethod for improving monopulse processing of aperture segment outputs
US6686874 *Sep 23, 2002Feb 3, 2004Eads Deutschland GmbhProcess for calibrating radar signals at subapertures of an antenna of two-channel SAR/MTI radar system
US7812761 *Mar 13, 2008Oct 12, 2010Fujitsu LimitedDSSS radar, method implemented by radar and computer-readable storage medium
CN102323579A *Aug 12, 2011Jan 18, 2012西安天伟电子系统工程有限公司Height measurement method for continuous wave search radar
Classifications
U.S. Classification342/152, 342/188, 342/194
International ClassificationG01S13/44
Cooperative ClassificationG01S13/4436
European ClassificationG01S13/44E
Legal Events
DateCodeEventDescription
Feb 24, 1987ASAssignment
Owner name: THOMSON-CSF, 173, BL. HAUSSMANN 75008 PARIS FRANCE
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST.;ASSIGNOR:POUX, JEAN P.;REEL/FRAME:004672/0201
Effective date: 19870211
Owner name: THOMSON-CSF, FRANCE
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:POUX, JEAN P.;REEL/FRAME:004672/0201
Owner name: THOMSON-CSF, FRANCE
Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:POUX, JEAN P.;REEL/FRAME:004672/0201
Effective date: 19870211
May 6, 1991FPAYFee payment
Year of fee payment: 4
May 23, 1995FPAYFee payment
Year of fee payment: 8
Jul 6, 1999REMIMaintenance fee reminder mailed
Dec 12, 1999LAPSLapse for failure to pay maintenance fees
Feb 22, 2000FPExpired due to failure to pay maintenance fee
Effective date: 19991215