Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Diss / (Springer Series in Information Sciences 25) S. Haykin, J. Litva, T. J. Shepherd (auth.), Professor Simon Haykin, Dr. John Litva, Dr. Terence J. Shepherd (eds.)-Radar Array Processing-Springer-Verlag

.pdf
Скачиваний:
68
Добавлен:
27.03.2016
Размер:
14.79 Mб
Скачать

304 B.D. Steinberg

Signal

I

I R, Ri

L

Beamformer

 

Array output

Fig. 7.7. Array model and signal processor configuration for self-synchronization. (After [7.6])

The spatial correlation algorithm was developed to handle these cases [7.6]. Figure 7.7 sketches a periodic array of N elements steered to 80 the angie from broadside. The phase-shift fjJn in the nth channel is nkduo where k = 21t/Ais the wavenumber, d is the element spacing and Uo = sin 80 , There is a phase error ofjJn in the nth channel, resulting both from element position error and electrical mistuning. The left channel is considered the reference. Its phase error is arbitrarily assigned the value of O. In each channel except the reference channel is a phase-error correcting weight. The sum is the array output.

Information for determining the correct values of the weights is obtained from the products of signals in adjacent channels. The first product is eoet; the second product is el e! ; and so on. Each product is integrated, or averaged over J range bins:

~

 

1

J

.et .

(7.1)

Ro

.

1 = -

L eo

 

J

j= 1

,J ,J

 

is an estimate of the correlation coefficient of the radiation field between the zeroth and first antenna elements modified by the phase error OfjJ. Its expectation is

 

(7.2)

Therefore,

 

Ro.l :::: Ro, 1 exp( -jofjJd

 

R1,2 :::: Rl,2 exP[j(OfjJl - OfjJ2)J

 

Rn•n +l :::: Rn.n+ l exp[j(ofjJn - ofjJn+dJ·

(7.3)

306 B.D. Steinberg

2687

Q)

Cl c: co a:

..,,;...... ..'.:-

""'. .

...•

.'.,

..

 

J-..... ,

. .

 

 

• •

' • ":':rl

• , .....

r~

 

.'AII.

 

••• .. ~o o· '.'

a

 

b

2627'--_______~

-16

Milliradians

4

 

 

Fig. 7.9. Comparison of the dominant scatterer algorithm <a) and the spatial correlation algorithm

(b) when no range bin contained a suitably prominent reflector. Observe that the DSA fails completely while the SCA indicates that the target is an airplane. [7.10]

O~---T-----~----~-----'

(I)

 

 

 

 

 

 

-10

 

 

 

 

-20

 

 

 

 

dB

 

 

 

 

-30

 

 

 

 

-40

 

 

 

 

-50~----~----~--__~____-J

-w

o

 

-w

W

20

 

Milliradians

 

Milliradians

 

Fig. 7.10a, b. One-dimensional images of a comer reflector. Adaptive beamforming is based upon echoes from farmland. The dominant scatterer algorithm (a) fails to self-calibrate the array. The spatial correlation algorithm (b) succeeds; its image exhibits a main lobe and sidelobe fall-otT. [7.10]

[7.8J. In this case, the dominant scatterer algorithm (left) failed, while the image of the spatial correlation algorithm indicates that the target is an airplane, permits its length and wing span to be measured, and provides other clues as to its class. Unlike the earlier comparison, diversity combining was not employed in Fig. 7.9. The comparison is made between two primal images. Figure 7.10 compares the two algorithms when ground clutter echoes drove the system. This figure compares the one-dimensional images of a comer reflector when the dominant scatterer algorithm (left) and the spatial correlation algorithm formed weight vectors based upon measurements of echoes from farmland. The superiority of the spatial correlation algorithm is evident in this case; the central region of the image is prominent and distinct, and the sidelobe artifacts falloff rapidly.

7. The Radio Camera

307

7.6 Number of Elements

Array thinning is essential in the radio camera. The resolution of its array is determined by its length, and not the number of elements. The latter determines its SNR gain, and its sidelobe or grating lobe properties. The SNR gain of a receiving array is directly proportional to N. The effect upon sidelobe level and the existence of grating lobes is more complicated, and is our primary concern. Grating lobes are formed when thinning is periodic, i.e., the inter-element space d remains constant, but is considerably larger than the "filled" spacing A/2. The radiation pattern

(7.5)

consists of the patternfo(u) of the filled array repeated at intervals of Ajd. There are two conditions under which (7.5) is an acceptable pattern for

imaging purposes. The first exists when the transmitter beamwidth is smaller than the lobe spacing A./d. This occurs when the transmitting antenna is larger than d, in which case, only a single lobe of(7.5) is illuminated by the transmitter. Because no scatterers outside the central lobe are excited, no lobe ambiguity results, nor does undesired back-scattered energy from the central lobe enter the system. This procedure was used in taking data for the image of Fig. 7.3. The transmitting antenna was a 1.2 m diameter dish; the receiving array sampled the radiation field every 25 cm, which was 8 wavelengths.

The second condition under which grating lobes are acceptable is when the target is smaller than A/d, and is free of a clutter background. In this case, the transmitter may illuminate a sector containing several ofthe receiving lobes, but because no scatterers are outside the central lobe, no ambiguity results. The data for Fig. 7.5 were obtained in this manner. The transmitting antenna, again, was a 1.2 m dish; the receiving samples were spaced at about 0.75 m, or 24 wavelengths.

Periodic thinning and grating lobe suppression, as previously described, preserves satisfactory sidelobe characteristics in the neighborhood of the main lobe. When circumstances do not permit periodic thinning, an aperiodic procedure such as random thinning must be employed. Aperiodic thinning destroys grating lobes, but not the energy in the grating lobes. The result is undesirably high sidelobes. The average sidelobe power level is N- 1, and the peak sidelobe is typically about 10 dB higher [7.1]. When random thinning is employed, the system must be designed to tolerate the high sidelobe level.

One favorable aspect of random thinning is that the resolution and sidelobe properties of the radiation pattern are almost totally decoupled. The angular resolution remains A./L, where L is the length over which the elements are randomly deployed. The average sidelobe level is N - 1 independent of the length, while the peak sidelobe grows very slowly with length; the relationship is logarithmic. Thus, once the penalty for using a random array is paid, and N is

308 B.D. Steinberg

determined by the tolerable sidelobe level, the aperture size can be made almost arbitrarily large. In essence, the higher the resolution, the higher the thinning factor. N is determined by the tolerable sidelobe level which, in tum, is governed by the required dynamic range in the image. For an imaging problem in which a 20 dB dynamic range is desired, about 1000 elements are required.

7.7 Number of Bits per Sample

Savings can also be affected by reducing b, the number of bits per sample. Figure 7.11 is a repeat of Fig. 7.3, but with a different data quantization procedure. For Fig. 7.3, the I and Q samples were each quantized to 8 bits. For Fig. 7.11, they were each quantized to 1 bit, which is a data reduction of 8: 1. Furthermore, the Fourier kernel of the image processor was similarly quantized to 2 bits, thereby eliminating the need for all complex multiplication in the signal processor. Comparison of the images show them to be nearly identical, and equally good. Their resolutions are the same, and the pixel intensities and locations are almost identical. The reader is referred to [7.9] for the theory of data compression in microwave imaging, and examples of several combinations of quantization of amplitude and phase.

7.8 Data Truncation

Array thinning and bit compression offer about three orders of magnitude reduction in data handling for the example given in the introduction. While it is a large reduction, it is still insufficient, for the reduced NBb product still exceeds

,..

 

 

.,. c

 

 

\'

 

-4490

 

". ,.

I. ~. !~.

 

\..

 

 

 

 

,,·tll~

...

,.:~

 

-4450

 

 

.••~I

J~,

 

"

 

 

 

.,

 

 

 

-4400

 

 

 

\"

"~ ,

·t· "

.' ,4350

I

I

 

,

,. ","w:

 

I

 

 

I

 

40

0

10

 

20

 

 

30

 

Azimuth (radians)

Fig. 7.11. Example of bit compression and simplification of the radio camera signal processor. Image of the housing development of Fig. 7.3. Input data to the phased array are hard limited and quantized to two bits of phase. The kernel of the signal processor is also quantized to two bits of phase. (From [7.9])

7. The Radio Camera

309

1011 bps. The input data to the high-resolution signal processor must be truncated to reduce the average rate to an acceptable level. In discussing the problem, [7.9] recommends a low or conventional resolution screening process to discover those targets or volume elements in space for which high-resolution imagery is desirable. Only those targets are processed by the high-resolution signal processor. The system consists of a scanning radar operating in conventional mode. Input data are also passed to a temporary buffer store with highresolution capability and limited capacity. The input data reside in store only long enough for the main receiver to decide whether or not a high-resolution image is called for. Most of the data are discarded. The undiscarded data are read out at a modest rate acceptable to the signal processor. By these three procedures, the total average data rate may be kept within bounds.

7.9 Conclusions

The radio camera is shown to produce very high-resolution images of targets and scenes from microwave radar data. The huge aperture that is necessary to achieve the high-resolution is self-calibrated by a process called adaptive beamforming (ABF). ABF deduces the required corrections for the phase errors in the array, directly from the back-scattered radiation field.

One ABF procedure utilizes the known spherical wavefront of the reradiation from a point reflector. It is called the dominant scatterer algorithm. The second procedure operates upon the spatial correlation function of the reradiation field from a quasi-homogeneous distribution of scatterers. It is called the spatial correlation algorithm.

Data rate reduction is essential for practical utilization of the radio camera techniques. Array thinning, bit compression, and data truncation are discussed for this purpose.

References

7.1B.D. Steinberg: Principles ofAperture and Array System Design: Including Random and Adaptive Arrays (Wiley, New York 1976)

7.2B.D. Steinberg: Microwave Imaging with Large Antenna Arrays: Radio Camera Principles and Techniques (Wiley, New York 1983)

7.3A.V. Oppenheim, J.S. Lim: The Importance of Phase in Signals. Proc. IEEE 69, 529-541 (1981)

7.4B.D. Steinberg, J. Tsao, D. Carlson, T. Seeleman: Imaging experiments. Valley Forge Research Center Quarterly Progress Report No. 45, University of Pennsylvania (1984) pp. 1-13

7.5B.D. Steinberg: Microwave imaging of aircraft. Proc. IEEE 76, 1578-1592 (1988)

7.6E.H. Attia, B.D. Steinberg: Self-Cohering Large Antenna Arrays Using the Spatial Correlation Properties of Radar Clutter. IEEE Trans. AP-37, 30-38 (1989)

7.7E.H. Attia: Phase synchronizing large antenna arrays using the spatial correlation properties of radar clutter. Ph.D. Dissertation, University of Pennsylvania (1984)

310 B.D. Steinberg

7.8B. Kang, B.D. Steinberg: Self-calibration of phased arrays using the modified Muller-Buffing- ton theorem and the spatial correlation algorithm, Valley Forge Research Center Quarterly Progress Report No. 54, University of Pennsylvania (1988) pp. 32-45

7.9B.D. Steinberg: A theory of the effect of hard limiting and other distortions upon the quality of microwave images. IEEE Trans. ASSP-3S, 1462-1472 (1987)

7.10B.D. Steinberg, H.M.S.: Microwave Imaging Techniques (Wiley, New York 1991)

Subject Index

a posteriori residual see residual a priori residual see residual

Accuracy

227

 

 

 

 

 

 

 

 

 

 

Adaptive

 

 

 

 

 

 

 

 

 

 

 

 

-

angular response

 

62

 

 

 

 

 

 

-

antenna

 

155, 157, 158, 159, 164,210, 213

-

array

153-155, 187

 

 

 

 

 

 

-

beamforrner

 

156

 

 

 

 

 

 

 

 

-

beamforming (ABF)

1, 153, 157, 160, 166,

 

168, 173,

175,

185,210,213,215,218,225,

 

227,228

 

 

 

 

 

 

 

 

 

 

 

-

beam space

296, 301, 302, 309

 

 

 

-

cancellation

 

173

 

 

 

 

 

 

 

 

-

classical

251 ; 254

 

 

 

 

 

 

 

-

closed-loop superresolution algorithm

82

-

conventional

 

249

 

 

 

 

 

 

 

-

nulling

227

 

 

 

 

 

 

 

 

 

 

Akaike

Information

Criterion

(AIC)

15,

59,

 

67,80

 

 

 

 

 

 

 

 

 

 

 

 

Alternating projection (AP) method

77, 126

Angular frequency

86

 

 

 

 

 

 

Antenna array see array

 

 

 

 

 

 

Antenna elements

297

 

 

 

 

 

 

-

array of I

 

 

 

 

 

 

 

 

 

 

 

Anti-phase multipath see multipath

 

 

 

AR see autoregressive

 

 

 

 

 

 

Arithmetic precision

173

 

 

 

 

 

ARMA see autoregressive moving average

 

Array (processor, network)

 

 

 

 

 

-

Gentleman-Kung

 

216, 226,228

 

 

 

-

linear

226

 

 

 

 

 

 

 

 

 

 

-

rectangular

227

 

 

 

 

 

 

 

 

-

square

226

 

 

 

 

 

 

 

 

 

 

-

systolic see systolic

 

 

 

 

 

 

 

-

trapezoidal

183, 233, 234

 

 

 

 

-

triangular

(-systolic)

156, 166,

167,

168,

 

169,

173,

179,

180,

181,

184,

185,

187,

189,

 

195,

196,

197, 200,

202,

203,

204,

205,

206,

 

207,

209,

210,

211,

212,

215,

221,

225,

226,

 

230, 233, 240

 

 

 

 

 

 

 

 

 

 

-

wavefront

155, 155-157, 210, 211, 213

 

Array

 

 

 

 

 

 

 

 

 

 

 

 

 

-

calibration errors

 

43

 

 

 

 

 

 

-

covariance

106

-

filled

307

 

 

-

(linear) antenna

157, 159, 174, 210, 249

-

manifold

105

 

-

phased

 

195, 299

-

planar

 

249

 

-

response

 

105

 

-

simple

 

250

 

-

super

250

 

-

thinning

 

307, 308, 309

-

unambiguous

106

Asymptotic equivalence 121, 122, 125, 125 Asymptotic robustness 118 Autoregressive (AR) process 58, 88 Autoregressive moving average (ARMA)

 

model

60

 

Auxiliary

 

-

channels

159, 161, 194, 221, 227

-

elements

349

-

input

188

-

signal

158, 160, 161

Averaged likelihood ratio 19

-

test

II

 

Back-substitution

164,168, 169,173,179,187,

201, 203

 

 

 

 

 

 

 

 

Banded matrix see matrix

 

 

 

 

Beam pattern

154, 158, 159

 

 

 

Beam steering

218

 

 

 

 

 

Beamforrner

(beamforrning) 154,

218,

220,

222,225

 

 

 

 

 

 

 

 

Bessel function, modified, of zero order

9

Binary phase coding

91

 

 

 

 

Bit compression

 

308, 309

 

 

 

 

Blocking matrix

 

197, 198, 199,200

 

 

Boundary cells

 

155,

166,167,168,

171, 178,

184,

185,

186,

187,

188,

204,

207,

208,

209,

218,221,234,235,238,241,262

Boundary nodes, (proct:ssors) see boundary cells

312

Subject Index

 

 

 

 

 

Broad-band

 

157,222

 

 

 

 

 

-

adaptive beamfonning

220-222

 

 

Burg algorithm

88

 

 

 

 

 

 

Burg-Levinson method

59

 

 

 

Calibration signal

85

 

 

 

 

 

Canonical

157, 160, 191, 220, 236

 

 

-

combiner

 

204

 

 

 

 

 

 

-

configuration (form)

158-160, 189

 

-

least-squares

157

 

 

 

 

 

-

problem

157-163, 192, 216, 219, 227

 

-

processor

 

189, 190, 191, 192, 194, 195, 196,

 

197,198,200,233,235

 

 

 

 

Capon-Pisarenko-type methods

62

 

 

Cell (node, processor)

 

 

 

 

 

-

final processing

262

 

 

 

 

-

frozen

184

 

 

 

 

 

 

 

Centre vector

223

 

 

 

 

 

 

Cholesky square root

202

 

 

 

 

Clock skew

 

156, 168, 181, 204, 206, 209, 211,

 

224

 

 

 

 

 

 

 

 

 

 

Closed loop

 

154, 157, 162

 

 

 

Clutter

303,306, 307

 

 

 

 

 

Coherent

 

 

 

 

 

 

 

 

 

 

-

signals

118

 

 

 

 

 

 

 

-

targets

57

 

 

 

 

 

 

 

 

Coherent radar detection

6-16, 133-137

 

-

signal and noise model for 6-7

 

 

-

of targets with known directions

7-13

 

-

of targets with unknown directions 13-16

Complex amplitude

73, 93

 

 

 

-

envelope

103

 

 

 

 

 

 

 

-

estimation

94

 

 

 

 

 

 

 

Complex arithmetic

227

 

 

 

 

Complex gradient

161

 

 

 

 

Computational complexity

130

 

 

 

Condition number

162, 164,260

 

 

Conditional density

111

 

 

 

 

Consistent estimate

115, 116, 122, 136

 

Constant false alarm rate (CFAR)

7

 

Constrained least-squares see least-squares

 

Constraint

154,

157,

159,

160,

191,

193,

196,

 

197, 198, 201, 204, 213, 226, 233, 234, 235

-

array

215

 

 

 

 

 

 

 

 

-

cell 208, 209

 

 

 

 

 

 

 

-

column 207, 209

 

 

 

 

 

-

look direction

173

 

 

 

 

 

-

matrix

193, 196, 200, 233, 235

 

 

-

multiple

193-196, 199, 226

 

 

 

-

post-processor

160,207,226,235

 

-pre-processor 173, 189-200,201,211,225, 232.233

-vector 190, 192, 193, 207, 209, 233, 235

Control bit

205, 206, 207

 

Conventional beamfonning 52

Conventional resolution limit

49

Convergence

154, 162,225

 

Comer reflector

297,298,306

 

Correlated

 

 

 

-

noise

76

 

 

 

-

targets

57, 75, 85

 

Correlation

 

 

 

-

coefficient

304, 305

 

-

function

90

 

 

Covariance filter

229

 

Covariance matrix 161, 162,

163, 173, 190,

 

201,202,232

 

 

-

error

228

 

 

 

-

estimation

56

 

 

Cramer-Rao (lower) bound (CRLB) 77, 91, 113, 114, 117, 125

Cross-correlation vector 161 Cube-edge-Iength (CEL) test 68

Data matrix 163, 185, 187, 190, 193,200,201, 202,204,216,227,230,234,235

-weighted 236 Data

-data skew see clock skew

-

truncation

308, 309

-

vector

195

 

Decoupled beams 75

Degrees-of-freedom 249

DESAS

81

 

Desired signal

154,158,159,173,174

Detection

1, 73, 141

-coherent radar see coherent radar detection

-noncoherent radar see noncoherent

-passive radar see passive

-radar 3

Deterministic nulling 50

Difference beam weighting with array

antennas

54

 

 

Diffraction-limited image

299

Direct-residual extraction

175-180,227

-

a posteriori

178-179

 

-

a priori

179-180

 

Direction of arrival (DOA)

104

DML see maximum likelihood

Doppler

 

 

 

-

filter bank

51, 86

 

-

frequency

51

 

Dual (dual-cell) 278

 

Efficient estimator 113, 115

Eigenbeams

250