没有合适的资源?快使用搜索试试~ 我知道了~
波束赋形
资源详情
资源评论
资源推荐
Beamforming:
A
Versatile
Approach to Spatial Filtering
Barry
D.
Van
A
beamformer
is
a processor used in conjunction with an
array of sensors to provide a versatile form of spatial filtering.
The sensor array collects spatial samples of propagating wave
fields, which are processed by the beamformer. The objective
is
to estimate the signal arriving from a desired direction in
the presence of noise and interfering signals.
A
beamformer
performs spatial filtering to separate signals that have over-
lapping frequency content but originate from different spatial
locations. This paper provides an overview of beamforming
from a signal processing perspective, with an emphasis on re-
cent research. Data independent, statistically optimum, adap-
tive, and partially adaptive beamforming are discussed.
1.
INTRODUCTION
The term beamforming derives from the fact that early
spatial filters were designed to form pencil beams (see
polar plot in Fig.
1
.I)
in order to receive a signal radiating
from
a
specific location and attenuate signals from other
locations. ”Forming beams” seems to indicate radiation of
energy; however, beamforming
is
applicable to either
radiation or reception of energy. In this paper we discuss
formation of beams for reception.
Systems designed to receive spatially propagating sig-
nals often encounter the presence of interference signals.
If the desired signal and interferers occupy the same tem-
poral frequency band, then temporal filtering cannot be
used to separate signal from interference. However, the
desired and interfering signals usually originate from dif-
ferent spatial locations. This spatial separation can be ex-
ploited to separate signal from interference using
a
spatial
filter at the receiver. Implementing
a
temporal filter requires
processing of data collected over
a
temporal aperture.
Similarly, implementing
a
spatial filter requires processing
of data collected over
a
spatial aperture.
Several applications that employ spatial filtering of data
are listed in Table
1
.1.
Fig.
1
.I
illustrates
a
microwave com-
munications antenna that employs
a
continuous spatial
aperture to accomplish spatial filtering with a single an-
tenna. Fig.
1.2
depicts a low frequency-towed sonar array
in which the spatial aperture
is
obtained through
a
dis-
crete spatial sampling by an array
of
sensors. When the
spatial sampling
is
discrete, the processor that performs
the spatial filtering
is
termed
a
beamformer. Typically
a
4
IEEE
ASSP MAGAZINE APRIL
1988
Veen and Kevin
M.
Buckley
’igure
1
.I
A
continuous spatial aperture provides one
nechanism for spatial filtering. Illustrated is a parabolic
nicrowave antenna system. The antenna dish provides
;he spatial aperture over which energy is collected. This
3nergy
is
reflected to the antenna feed. The dish and feed
iperate as a spatial integrator. The energy from a far field
source located directly in front
of
the antenna arrives at
;he feed temporarily aligned [i.e., all source-to-feed path
engths are equal1 and
is
coherently summed. In general,
?nergy from other sources will arrive at the feed via
iariable length paths, and add incoherently.
A
polar plot
if
a typical antenna beampattern (i.e., power gain
vs.,
in
;his case, azimuth angle1
is
shown for a selected fre-
iuency and for the elevation angle at which the antenna
is
minted.
0740-746718810400-0004$01
.OOGI9881
EEE
TABLE
1.1
ARRAYS
AN0 BEAMFORMERS PROVIDE AN EFFECTIVE AND VERSATILE MEANS
OF
SPATIAL
FILTERING.
THIS
TABLE LISTS A NUMBER OF APPLICATIONS OF SPATIAL FILTERING,
GIVES
EXAMPLES OF ARRAYS AN0 BEAMFORMERS, AND PROVIDES A FEW KEY REFERENCES.
Application Description References
RADAR
phased-array RADAR; air traffic Brookner [19851; Haykin 119851;
SONAR source localization and classification Knight et al. [19811; Owsley
Communications
directional transmission and Mayhan 119761; Compton
119781; Adams et al. 119801
control; synthetic aperture RADAR
Munson et al. (19831
[I9851
reception; sector broadcast
in
satellite communications
Imaging ultrasonic; optical; tomographic Macovski [19831; Pratt I19781;
Geophysical earth crust mapping;
oil
exploration Justice [I9851
Exploration
Astrophysical
high
resolution imaging of the Readhead [19821; Yen [I985
Exploration universe
Biomedical fetal heart monitoring; tissue
Kak [I9851
Widrow et al. 119751; Gee et al.
[19841; Peterson
et
al. [I9871
hyperthermia; hearing aids
beamformer linearly combines the spatially sampled time
series from each sensor to obtain a scalar output time
series in the same manner that an FIR filter linearly
combines temporally sampled data. Two principal advan-
tages of spatial sampling with an array of sensors are dis-
cussed below.
Spatial discrimination capability depends on the size of
the spatial aperture; as the aperture increases, discrimi-
nation improves. The absolute aperture size
is
not impor-
tant, rather its size in wavelengths
is
the critical parameter.
A
single physical antenna (continuous spatial aperture)
capable of providing the requisite discrimination
is
often
practical for high frequency signals since the wavelength
is
short. However, when low frequency signals are of in-
terest, an array of sensors can often synthesize a much
larger spatial aperture than that practical with
a
single
physical antenna.
A
second very significant advantage of using an array of
sensors, relevant at any wavelength,
is
the spatial filtering
versatility offered by discrete sampling. In many applica-
tion areas it
is
necessary to change the spatial filtering
function in real time to maintain effective suppression of
interfering signals. This change
is
easily implemented in
a
discretely sampled system by changing the way in which
the beamformer linearly combines the sensor data.
Changing the spatial filtering function of a continuous ap-
erture antenna
is
impractical.
The purpose of this paper
is
to describe beamforming
from
a
signal processing perspective, provide an overview
of beamformer design, and briefly discuss perform-
ance and implementation issues with an emphasis on re-
cent research. The paper begins with a section devoted to
defining basic terminology, notation, and concepts. Suc-
ceeding sections cover data-independent, statistically op-
timum, adaptive, and partially adaptive beamforming. We
then provide
a
brief discussion of implementation issues
and conclude with a summary.
Throughout the paper we use familiar methods and
techniques from FIR filtering to provide insight into vari-
ous aspects of spatial filtering with a beamformer. How-
ever, in some ways beamforming differs significantly from
FIR filtering. For example, in beamforming a source of
energy has several parameters that can be of interest:
range, azimuth and elevation angles, polarization, and
temporal frequency content. Different signals are often
mutually correlated as a result of multipath propagation.
The spatial sampling
is
often nonuniform and multi-
dimensional. Uncertainty must often be included in char-
acterization of individual sensor response and location,
motivating development of robust beamforming tech-
niques. These differences indicate that beamforming rep-
resents a more general problem than FIR filtering and as
a
result, more general design procedures and processing
structures are common.
Rather than making
a
futile attempt at attributing devel-
opments due to many different researchers in beam-
forming, we refer the reader to the following references:
books-J.
W.
R.
Griffiths et al., ed. [19731, Hudson [19811,
Monzingo and Miller [19801, Haykin, ed. [19851, Compton
[19881; special issues
-/€€E
Transactions on Antennas
and
Propagation [19761, [19861, journal
of
Ocean Engineering
[19871; tutorial
-
Gabriel [1976]; and bibliography- Marr
[19861. Papers devoted to beamforming are often found in
the
IEEE
Transactions on: Antennas
and
Propagation,
5
APRIL
1988
IEEE
ASSP MAGAZINE
Acoustics, Speech,
and
Signal Processing, Aerospace and
Electronic Systems, and in the lournal
of
the Acoustical
Society
of
America. There
is
a vast body of literature on
various aspects of beamforming and we can only refer-
ence a subset in support of our discussions. We often
refer to the
FIR
filtering literature in our discussions of
beamforming, since their histories are both parallel and
overlapping
.
II.
BASIC TERMINOLOGY AND CONCEPTS
In
this section we introduce terminology and concepts
employed throughout the paper. We begin with a subsec-
tion that defines the beamforming operation and discusses
spatial filtering. The next subsection, entitled “Second
Order Statistics,” develops representations for the covari-
ance of the data received at the array and discusses dis-
6
IEEE
ASSP MAGAZINE
APRIL
1988
tinctions between narrowband and broadband beam-
forming. The final subsection defines various types of
beamformers.
A.
Beam forming and Spatial Filtering
Fig.
2.1
depicts two beamformers. The first, which
samples the propagating wave field in space,
is
typically
used for processing narrowband signals. The output at
time k, y(k),
is
given by a linear combination of the data at
the
J
sensors at time k:
I
y(k)
=
C
w?xi(k)
(2.1)
L=l
where
*
represents complex conjugate. It
is
conventional
to multiply the data by conjugates of the weights to sim-
plify notation. We assume throughout that the data and
weights are complex since in many applications
a
quadra-
ture receiver
is
used at each sensor to generate in phase
and quadrature
(I
and
Q)
data. Each sensor
is
assumed
to
have any necessary receiver electronics and an
ND
con-
verter if beamforming
is
performed digitally.
The second beamformer in Fig. 2.1 samples the propa-
gating wave field in both space and time and
is
often used
when signals of significant frequency extent (broadband)
are of interest. The output in this case can be expressed as
(2.2)
where
K
-
1
is
the number of delays in each of the
J
sensor
channels. If the signal at each sensor
is
viewed as an input,
then a beamformer represents a multi-input single out-
put system.
It
is
convenient to develop notation which permits us to
treat both beamformers in Fig.
2.1
simultaneously. Note
that (2.1) and (2.2) can be written as
y(k)
=
wHx(k).
(2.3)
by appropriately defining a weight vector
w
and data vec-
tor x(k). We use lower and upper case boldface to denote
vector and matrix quantities respectively, and let super-
script H represent Hermitian (complex conjugate) trans-
pose. Vectors are assumed to be column vectors. Assume
that
w
and x(k) are
N
dimensional; this implies that
N
=
KJ
when referring to (2.2) and
N
=
J
when referring to (2.1).
Except for Section V on adaptive algorithms, we will drop
the time index and assume that
its
presence
is
understood
throughout the remainder of the paper. Thus (2.3)
is
writ-
ten as y
=
wHx.
Many of the techniques described in this
paper are applicable to continuous time as well
as
discrete
time beamforming.
The frequency response of an
FIR
filter with tap weights
w;,
0
5
p
5
J
and
a
tap delay of T seconds
is
given by
1
K-1
y(k)
=
c
c
WtpXl(k
-
P)
I=1
p=o
I
p=l
,.(@)
=
Wp+e-~~T(~-l)
(2.4a)
Alternatively
r(w)
=
wHd(w)
(2.4b)
where
wH
=
[w? w;.
. .
WT] and
d(w)
=
[I eJwTe12wT
...
e~il-l)wT
H
1
.
r(o) represents the response of the filter to
a
complex sinusoid of frequency
w
and
d(w)
is
a vector de-
scribing the phase of the complex sinusoid at each tap in
the
FIR
filter relative to the tap associated with wl.
Similarly, beamformer response
is
defined as the ampli-
tude and phase presented to a complex plane wave as a
function of location and frequency. Location
is
in general
a three dimensional quantity, but often we are only con-
cerned with one or two dimensional direction of arrival
(DOA). Throughout the remainder of the paper we do not
consider range. Fig. 2.2 illustrates the manner in which an
array of sensors samples a spatially propagating signal.
Assume that the signal
is
a
complex plane wave with
DOA
8
and frequency
w.
For convenience let the phase
be zero at the first sensor. This implies xl(k)
=
elwk and
xdk)
=
eJw[k-AI'o)l,
2
5
I
5
J.
Al(8)
represents the time delay
due to propagation from the first to the
Ith
sensor. Substi-
tution into (2.2) results in the beamformer output
I
Sensor
#I
-
Reference
I
Y--l
y(k)
=
2
'2'
Wtpe-lWIAl(o)+Pl
-
-
eIwkr(8,w) (2.5)
I=1
p=o
where
A,(@
=
0.
r(0,
w)
is
the beamformer response and
can be expressed in vector form as
r(8,
w)
=
wHd(8,
0).
(2.6)
The elements of
d(8,w)
correspond to the complex
exponentials e-lw[Al(o)+P'
.
In general it can be expressed
as
d(8,
w)
=
[I
eJWT2(o)
elWT3(8)
.
.
.
elwTN(@I]H
(2.7)
where the
7,(8),
2
5
i
5
N, are the time delays due to
propagation and any tap delays from the zero phase refer-
ence to the point at which the
ith
weight
is
applied. We
refer to
d(0,w)
as the array response vector. It
is
also
known as the steering vector or direction vector. Non-
ideal sensor characteristics can be incorporated into
d(8,w)
by multiplying each phase shift by
a
function
ado,
w),
which describes the associated sensor response
as
a function of frequency and direction.
APRIL
1988
IEEE
ASSP MAGAZINE
7
The beampattern
is
defined
as
the magnitude squared
of r(0,w). Note that each weight in w affects both the
temporal and spatial response of the beamformer. His-
torically, use of FIR filters has been viewed
as
providing
frequency dependent weights in each channel. This inter-
pretation
is
accurate but somewhat incomplete since the
coefficients in each filter also influence the spatial filtering
characteristics of the beamformer.
As
a multi-input single
output system, the spatial and temporal filtering that oc-
curs
is
a result of mutual interaction between spatial and
temporal sampling.
The correspondence between FIR filtering and beam-
forming
is
closest when the beamformer operates at
a
single
temporal frequency
w,)
and the array geometry
is
linear
and equi-spaced as illustrated in Fig.
2.3.
Letting the sen-
sor spacing be d, propagation velocity be c, and
0
repre-
sent
DOA
relative to broadside (perpendicular to the array)
we have
T,(@)
=
(i
-
1)
(d/c) sin
0.
In this case we identify
the relationship between temporal frequency
w
in
d(w)
(FIR filter) and direction
0
in
d(O,wO)
(beamformer)
as
w
=
wo(d/c) sin
0.
Thus, temporal frequency in an
FIR
fil-
m
s
m
m
I
Figure
2.3
The analogy between an equi-spaced omni-
directional narrowband line array and a single-channel
FIR
Filter
is
illustrated in this figure.
ter corresponds to the sine of direction in
a
narrowband,
linear equi-spaced beamformer. Complete interchange of
beamforming and
FIR
filtering methods
is
possible for this
special case provided the mapping between frequency
and direction
is
accounted for.
The vector notation introduced in (2.6) suggests
a
vector
space interpretation of beamforming. This point of view
is
useful both in beamformer design and analysis. We use it
here in consideration of spatial sampling and array geome-
try. The weight vector
w
and the array response vectors
d(0,
w)
are vectors in an
N
dimensional vector space. The
angles between w and
d(0,
w)
determine the response
r(0,
w).
For example, if for some
(0,
w)
the angle between
wand
d(0,
w)
is
90" (i.e., if w
is
orthogonal to
d(0,
w)),
then
the response
is
zero.
If
the angle
is
close to
o",
then the
response magnitude will be relatively large. The ability to
discriminate between sources at different locations and/or
frequencies, say
(01,
wl)
and
(02,
w2),
is
determined by the
angle between their array response vectors,
d(0,
w,)
and
The general effects
of
spatial sampling are similar to
temporal sampling. Spatial aliasing corresponds to an am-
biguity in source locations. The implication
is
that sources
at different locations have the same array response vector,
e.g., for narrowband sources
d(Ol,w0)
=
d(02,wo).
This
can occur if the sensors are spaced too far apart. If the
sensors are too close together, spatial discrimination suf-
fers as a result of the smaller than necessary aperture;
array response vectors are not well dispersed in the
N
dimensional vector space. Another type of ambiguity oc-
curs with broadband signals when
a
source at one location
and frequency cannot be distinguished from a source at a
different location and frequency, i.e.,
d(&,
w,)
=
d(02,
wJ.
For example, this occurs in
a
linear equi-spaced array
whenever
w1
sin
O1
=
w2
sin
02.
(The addition of temporal
samples at one sample prevents this particular ambiguity.)
A
primary focus of this paper
is
on designing response
via weight selection; however, (2.6) indicates that response
is
also a function of array geometry (and sensor character-
istics
if
the ideal omnidirectional sensor model
is
invalid).
In contrast with single channel filtering where
AID
con-
verters provide a uniform sampling in time, there
is
no
compelling reason to space sensors regularly. Sensor loca-
tions provide additional degrees of freedom in designing
a
desired response and can be selected
so
that over the
range of
(0,~)
of interest the array response vectors are
unambiguous and well dispersed in the N dimensional
vector space. Utilization of these degrees of freedom can
become very complicated due to the multidimensional
nature of spatial sampling and the nonlinear relationship
between r(0,
w)
and sensor locations. References discus-
sing array geometry design for response synthesis include
Unz [19561, Harrington [19611, lshimaru [19621, Lo [19631,
and Skolnik et al. [19641.
d(02,
~2)
(COX [19731).
B.
Second Order Statistics
Evaluation of beamformer performance usually involves
power
or
variance,
so
the second order statistics of the
剩余20页未读,继续阅读
qq_45884246
- 粉丝: 0
- 资源: 9
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
- 基于matlab实现用有限元法计算电磁场的Matlab工具 .rar
- 基于matlab实现有限元算法 计算电磁场问题 边界条件包括第一类边界和第二类边界.rar
- 基于matlab实现用于计算不同车重下的电动汽车动力性和经济性.rar
- 基于matlab实现遗传算法求解多车场车辆路径问题 有多组算例可以用.rar
- 浏览器.apk
- 基于matlab实现是一个matlab中的power system 中搭建的一个模型
- 基于JSP毕业设计-教学管理系统(源代码+论文).zip
- 基于JSP毕业设计-家政管理系统-毕业设计.zip
- 基于Python实现淘宝商品评论采集(含逆向)源代码
- 基于matlab实现多目标进化算法NSGAⅡ&Matlab讲解.rar
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功
评论0