跳转到主要内容年代pan>
下载文章年代pan>
伊朗库姆技术大学年代pan>
检查更新一个><一个rticle class="widget-listing people-also-looked-at side-article-related hidden">
人也看了 一个rt我cle>
在哪里<我nl我ne-formula id="inf3">
J米我><米o>=米o><米年代问rt><米row><米o>
量米o><米n>1米n>米row>米年代问rt>
。
beamformer输出是由
Z米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米年代up>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
y米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
,米o><米年代pace width="5em">
(米o><米n>3米n><米o年代tretchy="false">
)米o>米在h>
在哪里<我nl我ne-formula id="inf4">
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
波束形成的权向量的长度锰 和<年代up>H年代up>表示共轭转置。frequency-invariant获得期望的方向(φ <年代ub>d年代ub>,θ <年代ub>d年代ub>)。因此,我们对常见的无失真约束:
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
一个米我><米fencedopen="(" close=")">
f米我><米o>,米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
,米o><米年代ub><米row><米i>
θ米我>米row><米row>
d米text>米row>
=米o><米n>1米n><米o>。米o><米年代pace width="5em">
(米o><米n>4米n><米o年代tretchy="false">
)米o>米在h>
beampattern表达的反应beamformer平面波和定义
B米我><米fencedopen="(" close=")">
f米我><米o>,米o><米我>φ米我><米o>,米o><米我>θ米我><米o>,米o><米我米在hv一个r我一个nt="bold">
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米年代up>
一个米我>米row><米row>
H米我>米row>米年代up>
f米我><米o>,米o><米我>φ米我><米o>,米o><米我>θ米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
。米o><米年代pace width="5em">
(米o><米n>5米n><米o年代tretchy="false">
)米o>米在h>
为了提高可读性,我们有时使用更简洁的符号指导向量,<年代trong>一个年代trong>,波束形成的权向量,<年代trong>h年代trong>beampattern,<我nl我ne-formula id="inf5">
B米我><米row><米o年代tretchy="false">
(米o><米row><米我米在hvariant="bold">
h米我>米row><米o年代tretchy="false">
)米o>米row>米在h>
通过省略变量f ,φ ,θ 。
2.2性能的措施 我们现在开发宽带性能措施量化绕线和方向性频带。窄带和宽带输入信噪比(信噪比)
我米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米fr一个c><米row>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
(米o><米n>6米n><米o年代tretchy="false">
)米o>米在h>
和
我米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我><米o>=米o><米fr一个c><米row>
∫米o>米row><米row>
f米我>米row>米年代ub>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
d米我><米我>f米我>米row><米row>
∫米o>米row><米row>
f米我>米row>米年代ub>
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
,米o><米年代pace width="5em">
(米o><米n>7米n><米o年代tretchy="false">
)米o>米在h>
分别在哪里<我nl我ne-formula id="inf6">
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
的方差<我nl我ne-formula id="inf7">
X米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
和<我nl我ne-formula id="inf8">
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
加性噪声的方差在参考传感器。从<一个href="#e1">(1)一个>和<一个href="#e3">(3)一个>窄带和宽带输出信噪比
o米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米fr一个c><米row>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
一个米我><米fencedopen="(" close=")">
f米我><米o>,米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
,米o><米年代ub><米row><米i>
θ米我>米row><米row>
d米text>米row>
2米n>米row>米年代up>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o>米row><米row>
4米n>米row>米fenced>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
(米o><米n>8米n><米o年代tretchy="false">
)米o>米在h>
和
o米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我>米td><米td columnalign="left">
=米o><米fr一个c><米row>
∫米o>米row><米row>
f米我>米row>米年代ub>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
一个米我><米fencedopen="(" close=")">
f米我><米o>,米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
,米o><米年代ub><米row><米i>
θ米我>米row><米row>
d米text>米row>
2米n>米row>米年代up>
d米我><米我>f米我>米row><米row>
∫米o>米row><米row>
f米我>米row>米年代ub>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
=米o>米row><米row>
4米n>米row>米fenced>
∫米o>米row><米row>
f米我>米row>米年代ub>
ϕ米我>米row><米row>
X米我>米row>米年代ub>
f米我>米row>米fenced>
d米我><米我>f米我>米row><米row>
∫米o>米row><米row>
f米我>米row>米年代ub>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
,米o>米td>米tr>
(米o><米n>9米n><米o年代tretchy="false">
)米o>米在h>
分别在哪里<我nl我ne-formula id="inf9">
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
的相关矩阵<我nl我ne-formula id="inf10">
问米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
。给出了窄带和宽带阵列收益
G米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米fr一个c><米row>
o米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
我米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
(米o><米n>10米n><米o年代tretchy="false">
)米o>米在h>
和
G米我><米o>=米o><米fr一个c><米row>
o米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我>米row><米row>
我米我><米我米在hv一个r我ant="normal">
年代米我><米我米在hv一个r我ant="normal">
N米我><米我米在hv一个r我ant="normal">
R米我>米row>米frac>
,米o><米年代pace width="5em">
(米o><米n>11米n><米o年代tretchy="false">
)米o>米在h>
分别。
当<我nl我ne-formula id="inf11">
问米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
球各向同性噪声,<我nl我ne-formula id="inf12">
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
是由
Φ米我>米row><米row>
问米我>米row>米年代ub>
f米我>米row>米fenced>
=米o><米年代ub><米row><米i>
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
Γ米我>米row><米row>
d米text>米row>
f米我>米row>米fenced>
,米o><米年代pace width="5em">
(米o><米n>12米n><米o年代tretchy="false">
)米o>米在h>
的元素<我nl我ne-formula id="inf13">
Γ米我>米row><米row>
d米text>米row>
f米我>米row>米fenced>
是由
Γ米我>米row><米row>
d米text>米row>
f米我>米row>米fenced>
r米我><米o>,米o><米年代up>
r米我>米row><米row>
′米o>米row>米年代up>
=米o><米我米在hv一个r我一个nt="normal">
年代米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
n米我><米我米在hv一个r我ant="normal">
c米我><米fencedopen="(" close=")">
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
d米我>米row><米row>
r米我><米o>,米o><米年代up>
r米我>米row><米row>
′米o>米row>米年代up>
;米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>r米我><米o>,米o><米年代up>
r米我>米row><米row>
′米o>米row>米年代up>
≤米o><米我>米米我><米我>N米我><米o>;米o><米年代pace width="5em">
(米o><米n>13米n><米o年代tretchy="false">
)米o>米在h>
在哪里d <年代ub>r ,r ′年代ub>之间的欧几里得距离相应的传感器。
窄带的指向性因子(DF)被定义为窄带阵列增益的球面各向同性噪声(<一个href="#B3">Benesty et al ., 2018一个>),是由
D米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米fr一个c><米row>
1米n>米row><米row>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Γ米我>米row><米row>
d米text>米row>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
。米o><米年代pace width="5em">
(米o><米n>14米n><米o年代tretchy="false">
)米o>米在h>
我们定义频段宽带方向性因子<我nl我ne-formula id="inf14">
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
随着宽带阵列增益的球面各向同性噪声。的场景<我nl我ne-formula id="inf15">
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
是frequency-invariant,即<我nl我ne-formula id="inf16">
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
≡米o><米年代ub><米row><米i>
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
的紧凑表达式宽带方向性因素是由:
D米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
=米o><米fr一个c><米row>
1米n>米row><米row>
1米n>米row><米row>
f米我>米row><米row>
H米text>米row>
−米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
∫米o>米row><米row>
f米我><米o>=米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
f米我>米row><米row>
H米text>米row>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
Γ米我>米row><米row>
d米text>米row>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
=米o>米row><米row>
14米n>米row>米fenced>
1米n>米row><米row>
1米n>米row><米row>
f米我>米row><米row>
H米text>米row>
−米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
∫米o>米row><米row>
f米我><米o>=米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
f米我>米row><米row>
H米text>米row>
1米n>米row><米row>
D米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
。米o>米td>米tr>
(米o><米n>15米n><米o年代tretchy="false">
)米o>米在h>
在这里,我们定义宽带指向性指数(DI)频带<我nl我ne-formula id="inf17">
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
作为
D米我><米我>我米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
=米o><米n>10米n><米o>米o><米年代ub><米row><米i>
日志米我>米row><米row>
10米n>米row>米年代ub>
D米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
。米o><米年代pace width="5em">
(米o><米n>16米n><米o年代tretchy="false">
)米o>米在h>
绕线是一种常用的性能度量,反映了物理缺陷(数组的鲁棒性<一个href="#B3">Benesty et al ., 2018一个>)。它被定义为数组空间白噪声增益,即,<我nl我ne-formula id="inf18">
Φ米我>米row><米row>
问米我>米row>米年代ub>
(米o><米row><米我>f米我>米row><米o年代tretchy="false">
)米o>米row><米o>=米o><米年代ub><米row><米i>
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
(米o><米row><米我>f米我>米row><米o年代tretchy="false">
)米o>米row><米年代ub>
我米我>米row><米row>
米米我><米我>N米我>米row>米年代ub>
,在那里<年代trong>我年代trong>锰年代ub>是锰 ×锰 单位矩阵。的场景<我nl我ne-formula id="inf19">
ϕ米我>米row><米row>
问米我>米row><米row>
0米n>米row>米年代ub>
f米我>米row>米fenced>
frequency-invariant,窄带和宽带绕线的紧凑表达式:
W米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
=米o><米fr一个c><米row>
1米n>米row><米row>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
(米o><米n>17米n><米o年代tretchy="false">
)米o>米在h>
和
W米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
=米o><米fr一个c><米row>
1米n>米row><米row>
1米n>米row><米row>
f米我>米row><米row>
H米text>米row>
−米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
∫米o>米row><米row>
f米我><米o>=米o><米年代ub><米row><米i>
f米我>米row><米row>
l米text>米row>
f米我>米row><米row>
H米text>米row>
1米n>米row><米row>
W米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
d米我><米我>f米我>米row>米frac>
,米o><米年代pace width="5em">
(米o><米n>18米n><米o年代tretchy="false">
)米o>米在h>
分别。
2.3问题公式化 我们假定传感器的数量,锰 是给定的。我们的目标是找到最优位置的传感器和最佳beamformer重量最大化DI,同时保持一个恒定的波束宽度和可接受的绕线。
让b <年代ub>X年代ub>和b <年代ub>Y年代ub>表示在XZ半功率波束宽度和YZ平面,分别。然后,我们的优化问题的定义是
米米我><米我米在hv一个r我ant="normal">
一个米我><米我米在hv一个r我ant="normal">
x米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
米米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
z米我><米我米在hv一个r我ant="normal">
e米我>米row><米row>
x米我>米row><米row>
米米我>米row>米年代ub>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>米米我><米o>≤米o><米我>米米我>米td>米tr>
y米我>米row><米row>
n米我>米row>米年代ub>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>n米我><米o>≤米o><米我>N米我>米td>米tr>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
D米我><米我>我米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
年代米我><米我米在hv一个r我ant="normal">
u米我><米我米在hv一个r我ant="normal">
b米我><米我米在hv一个r我ant="normal">
j米我><米我米在hv一个r我ant="normal">
e米我><米我米在hv一个r我ant="normal">
c米我><米我米在hv一个r我ant="normal">
t米我><米年代pace width="0.3333em" class="nbsp">
t米我><米我米在hv一个r我ant="normal">
o米我>米td>米tr>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
一个米我><米fencedopen="(" close=")">
f米我><米o>,米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
,米o><米年代ub><米row><米i>
θ米我>米row><米row>
d米text>米row>
=米o><米n>1米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
b米我>米row><米row>
X米text>米row>
f米我>米row>米fenced>
=米o><米年代ub><米row><米i>
θ米我>米row><米row>
X米text>米row>
,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
b米我>米row><米row>
Y米text>米row>
f米我>米row>米fenced>
=米o><米年代ub><米row><米i>
θ米我>米row><米row>
Y米text>米row>
,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
W米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
=米o><米我>ℵ米我><米o>,米o>米td>米tr>
(米o><米n>19米n><米o年代tretchy="false">
)米o>米在h>
在哪里θ <年代ub>X年代ub>和θ <年代ub>Y年代ub>在XZ所需的半功率波束宽度和YZ平面,分别和ℵ ≥1可接受的宽带绕线。
3平面阵列Constant-Beamwidth Beampattern设计 本节将展示如何开发一个常数波束宽度的平面阵列基于线性阵列的设计。在XZ平面波束宽度可能不同于YZ平面的波束宽度。这可以通过设计两个线性阵列,每一个都有不同的波束宽度。
一个线性阵列beampattern设计是这样的,当躺在线性数组x 设在beampattern保持无失真约束在侧向方向,即。,(φ <年代ub>d年代ub>,θ <年代ub>d年代ub>)=(90°、90°)。它的波束宽度,b <年代ub>φ 年代ub>以XY平原。我们的目标是设计一个beampattern的半功率波束宽度θ <年代ub>X年代ub>在XZ平面θ <年代ub>Y年代ub>在YZ平面。假设给出了两个线性阵列(在下一节中,我们建议一个方法设计线性阵列,然而,本节设计并不局限于这个方法)。第一个线性数组米 传感器固定波束宽度的θ <年代ub>X年代ub>。表示位置的<我nl我ne-formula id="inf20">
{米o><米row><米年代ub><米row>
x米我>米row><米row>
米米我>米row>米年代ub>
}米o>米row>米row>
米米我><米o>=米o><米n>1米n>米row><米row>
米米我>米row>米年代ubsup>
与频率相关的权重<年代trong>h年代trong>X年代ub>(f )。第二个线性数组N 传感器固定波束宽度的θ <年代ub>Y年代ub>。表示位置的<我nl我ne-formula id="inf21">
{米o><米row><米年代ub><米row>
y米我>米row><米row>
n米我>米row>米年代ub>
}米o>米row>米row>
n米我><米o>=米o><米n>1米n>米row><米row>
N米我>米row>米年代ubsup>
与频率相关的权重<年代trong>h年代trong>Y年代ub>(f )。为了简化分析,我们假设x <年代ub>米 年代ub>=0的米 这y <年代ub>n 年代ub>=0的n 。我们位置平面阵列传感器基于线性阵列传感器的位置。我们变换使用的权重线性阵列平面beamformer得到权重。转换确保XZ和YZ平面波束宽度保持不变。转换并不依赖于频率,只有物理位置的传感器。变换使constant-beamwidth使用平面阵列波束形成与许多传感器设计简单的线性数组。
首先,我们展示如何使用克罗内克积波束形成(<一个href="#B2">Benesty et al ., 2019一个>非均匀平面阵列)。然后我们提供波束宽度如何使用较少的传感器保持不变;该方法使绕线和指向性之间的权衡。
3.1克罗内克积波束形成 从两个给定的线性数组,我们构造一个米 ×N 网格的x 坐标是<我nl我ne-formula id="inf22">
{米o><米row><米年代ub><米row>
x米我>米row><米row>
米米我>米row>米年代ub>
}米o>米row>米row>
米米我><米o>=米o><米n>1米n>米row><米row>
米米我>米row>米年代ubsup>
和y 坐标是<我nl我ne-formula id="inf23">
{米o><米row><米年代ub><米row>
y米我>米row><米row>
n米我>米row>米年代ub>
}米o>米row>米row>
n米我><米o>=米o><米n>1米n>米row><米row>
N米我>米row>米年代ubsup>
——形成一个平面阵列共有锰 传感器。让<年代trong>一个年代trong>P年代ub>(f ,φ ,θ )表示指导向量的长度锰 在是谁的条目<一个href="#e2">(2)一个>。表示由<年代trong>一个年代trong>X年代ub>(f ,φ ,θ ),<年代trong>一个年代trong>Y年代ub>(f ,φ ,θ )指导向量(各自的长度米 和N )的子阵,躺在x 和y 给出了轴的条目
e米我>米row><米row>
−米o><米我米在hv一个r我一个nt="script">
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>因为米我><米row><米o>(米o><米row><米我>φ米我>米row><米o>
)米o>米row><米我>
罪米我><米row><米o>(米o><米我>θ米我><米o>)米o>米row>米row>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>米米我><米o>≤米o><米我>米米我><米o>,米o><米年代pace width="5em">
(米o><米n>20.。米n><米o年代tretchy="false">
)米o>米在h>
和
e米我>米row><米row>
−米o><米我米在hv一个r我一个nt="script">
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
y米我>米row><米row>
n米我>米row>米年代ub>
米o><米我>罪米我><米row><米o>(米o><米row><米我>φ米我>米row><米o>
)米o>米row><米我>
罪米我><米row><米o>(米o><米row><米我>θ米我>米row><米o>
)米o>米row>米row>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>n米我><米o>≤米o><米我>N米我><米o>,米o><米年代pace width="5em">
(米o><米n>21米n><米o年代tretchy="false">
)米o>米在h>
分别。请注意,
一个米我>米row><米row>
P米text>米row>
=米o><米年代ub><米row><米i mathvariant="bold">
一个米我>米row><米row>
Y米text>米row>
⊗米o><米年代ub><米row><米i mathvariant="bold">
一个米我>米row><米row>
X米text>米row>
,米o><米年代pace width="5em">
(米o><米n>22米n><米o年代tretchy="false">
)米o>米在h>
克罗内克积⊗表示。
给出了子阵的beampatterns
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
X米text>米row>
=米o><米年代up>
一个米我>米row><米row>
X米text>米row>
H米我>米row>米年代up>
h米我>米row><米row>
X米text>米row>
(米o><米n>23米n><米o年代tretchy="false">
)米o>米在h>
和
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
Y米text>米row>
=米o><米年代up>
一个米我>米row><米row>
Y米text>米row>
H米我>米row>米年代up>
h米我>米row><米row>
Y米text>米row>
。米o><米年代pace width="5em">
(米o><米n>24米n><米o年代tretchy="false">
)米o>米在h>
线性阵列,因为它存在沿一维线,每2 d切片(包含1 d线)的beampattern是相同的。子数组躺在x 设在(y 设在),每一个2 d切片(包含x 设在(y 设在)beampattern是相同的;例如,XZ (YZ)和XY平面。因此,由constant-beamwidth设计的线性数组:
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
X米text>米row>
=米o><米fencedopen="{" close="">
1米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>=米o><米n>90年年米n><米o>°米o><米年代pace width="0.3333em" class="nbsp">
o米我><米我米在hv一个r我ant="normal">
r米我><米年代pace width="0.3333em" class="nbsp">
θ米我><米o>=米o><米n>0米n><米o>°米o><米o>,米o>米td>米tr>
1米n>米row><米row>
2米n>米row>米年代问rt>
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>,米o><米我>θ米我>米row>米fenced>
=米o><米fencedopen="(" close=")">
0米n><米o>°米o><米o>,米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
,米o>米td>米tr>
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
Y米text>米row>
=米o><米fencedopen="{" close="">
1米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>=米o><米n>0米n><米o>°米o><米o>,米o>米td>米tr>
1米n>米row><米row>
2米n>米row>米年代问rt>
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>,米o><米我>θ米我>米row>米fenced>
=米o><米fencedopen="(" close=")">
90年米n><米o>°米o><米o>,米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
Y米text>米row>
。米o>米td>米tr>
(米o><米n>25米n><米o年代tretchy="false">
)米o>米在h>
让<年代trong>h年代trong>3年代ub>是一个二维窗口构造成的外积两个1 d窗口:
h米我>米row><米row>
3米n>米row>米年代ub>
米米我><米o>,米o><米我>n米我>米row>米fenced>
=米o><米年代ub><米row><米i>
h米我>米row><米row>
1米n>米row>米年代ub>
米米我>米row>米fenced>
h米我>米row><米row>
2米n>米row>米年代ub>
n米我>米row>米fenced>
;米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>米米我><米o>≤米o><米我>米米我><米年代pace width="0.17em">
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>n米我><米o>≤米o><米我>N米我><米o>;米o><米年代pace width="5em">
(米o><米n>26米n><米o年代tretchy="false">
)米o>米在h>
在哪里<年代trong>h年代trong>1年代ub>和<年代trong>h年代trong>2年代ub>1dw我ndow年代和各自的长度吗米 和N 。例如,每个1 d窗口可以变得闷闷不乐皇帝那么2 d窗口是一个2 d Kaiser窗。的建设,<年代trong>h年代trong>3年代ub>两个变量是可分的。因此,如果我们安排它的值放入一个向量的长度锰 ,得到的向量可以制定的克罗内克积两个1 d窗口。我们设计的平面阵列波束形成权向量,<年代trong>h年代trong>P年代ub>,因为
h米我>米row><米row>
P米text>米row>
=米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
Y米text>米row>
⊗米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
X米text>米row>
。米o><米年代pace width="5em">
(米o><米n>27米n><米o年代tretchy="false">
)米o>米在h>
以下的方法(<一个href="#B2">Benesty et al ., 2019一个>)(而非均匀间隔的矩形数组)概括,给出了平面阵列的beampattern
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
P米text>米row>
=米o><米年代up>
一个米我>米row><米row>
P米text>米row>
H米我>米row>米年代up>
f米我><米o>,米o><米我>φ米我><米o>,米o><米我>θ米我>米row>米fenced>
h米我>米row><米row>
P米text>米row>
f米我>米row>米fenced>
=米o><米年代up>
一个米我>米row><米row>
Y米text>米row>
⊗米o><米年代ub><米row><米i mathvariant="bold">
一个米我>米row><米row>
X米text>米row>
H米我>米row>米年代up>
h米我>米row><米row>
Y米text>米row>
⊗米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
X米text>米row>
=米o><米fencedopen="(" close=")">
一个米我>米row><米row>
X米text>米row>
H米我>米row>米年代up>
h米我>米row><米row>
X米text>米row>
一个米我>米row><米row>
Y米text>米row>
H米我>米row>米年代up>
h米我>米row><米row>
Y米text>米row>
=米o><米我米在hv一个r我一个nt="script">
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
X米text>米row>
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
Y米text>米row>
,米o>米td>米tr>
(米o><米n>28米n><米o年代tretchy="false">
)米o>米在h>
第三平等,我们使用了克罗内克积身份
u米我>米row><米row>
1米n>米row>米年代ub>
⊗米o><米年代ub><米row><米i mathvariant="bold">
u米我>米row><米row>
2米n>米row>米年代ub>
H米我>米row>米年代up>
v米我>米row><米row>
1米n>米row>米年代ub>
⊗米o><米年代ub><米row><米i mathvariant="bold">
v米我>米row><米row>
2米n>米row>米年代ub>
=米o><米fencedopen="(" close=")">
u米我>米row><米row>
2米n>米row><米row>
H米我>米row>米年代ubsup>
v米我>米row><米row>
2米n>米row>米年代ub>
u米我>米row><米row>
1米n>米row><米row>
H米我>米row>米年代ubsup>
v米我>米row><米row>
1米n>米row>米年代ub>
。米o><米年代pace width="5em">
(米o><米n>29日日米n><米o年代tretchy="false">
)米o>米在h>
XZ的波束宽度和YZ平面由替代实现<一个href="#e25">(25)一个>成<一个href="#e28">(28)一个>:
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
P米text>米row>
=米o><米fencedopen="{" close="">
1米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>,米o><米我>θ米我>米row>米fenced>
=米o><米fencedopen="(" close=")">
0米n><米o>°米o><米o>,米o><米n>0米n><米o>°米o>米row>米fenced>
,米o>米td>米tr>
1米n>米row><米row>
2米n>米row>米年代问rt>
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>,米o><米我>θ米我>米row>米fenced>
=米o><米fencedopen="(" close=")">
0米n><米o>°米o><米o>,米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
,米o>米td>米tr>
1米n>米row><米row>
2米n>米row>米年代问rt>
,米o><米年代pace width="0.3333em" class="nbsp">
φ米我><米o>,米o><米我>θ米我>米row>米fenced>
=米o><米fencedopen="(" close=")">
90年米n><米o>°米o><米o>,米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
Y米text>米row>
。米o>米td>米tr>
(米o><米n>30.。米n><米o年代tretchy="false">
)米o>米在h>
每架飞机的beampattern维护所需的波束宽度和满足无失真约束。这表明我们可以分解的2 d问题设计一个constant-beamwidth beamformer平面阵列的两个一维问题的设计constant-beamwidth beamformers与线性阵列。克罗内克积是有效的只有当传感器被放置在所有锰 网格的位置。接下来,我们开发一种方法对于实现constant-beamwidth beamformer即使只有一个子集的网格位置传感器。
克罗内克Beamformer的计算复杂度:
在本文中,我们考虑计算复杂度随着乘法所需数量的方法。假设两个线性阵列,平面阵列波束形成加权计算<一个href="#e27">(27)一个>,要求锰 乘法/频率。让F 表示频率垃圾箱的数量。因此,计算复杂度<我nl我ne-formula id="inf24">
O米我><米row><米o年代tretchy="false">
(米o><米row><米我>米米我><米我>N米我><米我>F米我>米row><米o年代tretchy="false">
)米o>米row>米在h>
。
3.2线性平面体重变换 本节提出了一种方法,不需要在每个网格位置传感器,可以绕线和指向性之间的权衡。从两个线性阵列,我们构造米 ×N 网格,但这一次选择年代 的可能锰 职位的年代 传感器。条件对传感器的位置稍后讨论。
beampattern给出的
B米我><米fencedopen="(" close=")">
f米我><米o>,米o><米我>φ米我><米o>,米o><米我>θ米我><米o>,米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
P米text>米row>
f米我>米row>米fenced>
=米o><米under><米row>
∑米o>米row><米row>
米米我>米row>米under>
∑米o>米row><米row>
n米我>米row>米under>
一个米我>米row><米row>
P米text>米row>
*米o>米row>米年代up>
米米我><米o>,米o><米我>n米我>米row>米fenced>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
,米o><米年代pace width="5em">
(米o><米n>31日日米n><米o年代tretchy="false">
)米o>米在h>
求和只是在哪里年代 值(只有指数对哪里有一个传感器)和*表示复共轭。在这里,h <年代ub>P年代ub>(米 ,n ]beamformer重量传感器定位在(x <年代ub>米 年代ub>,y <年代ub>n 年代ub>)。方括号内的其他表达式在本节中类似的定义。
让我们检查beampattern在XZ平面。在海拔<我nl我ne-formula id="inf25">
θ米我><米o>=米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
:
B米我><米fencedopen="(" close=")">
f米我><米o>,米o><米我>φ米我><米o>=米o><米n>0米n><米o>,米o><米我>θ米我><米o>=米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
,米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
P米text>米row>
f米我>米row>米fenced>
=米o><米under><米row>
∑米o>米row><米row>
米米我>米row>米under>
∑米o>米row><米row>
n米我>米row>米under>
e米我>米row><米row>
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>罪米我><米fencedopen="(" close=")">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
=米o><米under><米row>
∑米o>米row><米row>
米米我>米row>米under>
e米我>米row><米row>
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>罪米我><米fencedopen="(" close=")">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
∑米o>米row><米row>
n米我>米row>米under>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
。米o>米td>米tr>
(米o><米n>32米n><米o年代tretchy="false">
)米o>米在h>
我们渴望的绝对值<一个href="#e32">(32)一个>就等于<我nl我ne-formula id="inf26">
1米n><米o>/米o><米年代问rt><米row><米n>
2米n>米row>米年代问rt>
,这正是振幅的线性阵列beampattern<我nl我ne-formula id="inf27">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
(25)一个>:
B米我><米fencedopen="(" close=")">
h米我>米row><米row>
X米text>米row>
=米o><米fencedopen="|" close="|">
∑米o>米row><米row>
米米我>米row>米under>
e米我>米row><米row>
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>罪米我><米fencedopen="(" close=")">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
h米我>米row><米row>
X米text>米row>
米米我>米row>米fenced>
=米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米年代问rt>
。米o><米年代pace width="5em">
(米o><米n>33米n><米o年代tretchy="false">
)米o>米在h>
因此,我们需要:
∑米o>米row><米row>
米米我>米row>米under>
e米我>米row><米row>
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>罪米我><米fencedopen="(" close=")">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
∑米o>米row><米row>
n米我>米row>米under>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
=米o><米fencedopen="|" close="|">
∑米o>米row><米row>
米米我>米row>米under>
e米我>米row><米row>
J米我><米fr一个c><米row>
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
x米我>米row><米row>
米米我>米row>米年代ub>
米o><米我>罪米我><米fencedopen="(" close=")">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
X米text>米row>
h米我>米row><米row>
X米text>米row>
米米我>米row>米fenced>
。米o><米年代pace width="5em">
(米o><米n>34米n><米o年代tretchy="false">
)米o>米在h>
达到<一个href="#e34">(34)一个>它足以要求:
∑米o>米row><米row>
n米我>米row>米under>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
=米o><米年代ub><米row><米i>
h米我>米row><米row>
X米text>米row>
米米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>米米我><米o>≤米o><米我>米米我><米o>。米o><米年代pace width="5em">
(米o><米n>35米n><米o年代tretchy="false">
)米o>米在h>
同样,达到理想的波束宽度在YZ平面,我们要求:
∑米o>米row><米row>
米米我>米row>米under>
h米我>米row><米row>
P米text>米row>
米米我><米o>,米o><米我>n米我>米row>米fenced>
=米o><米年代ub><米row><米i>
h米我>米row><米row>
Y米text>米row>
n米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>n米我><米o>≤米o><米我>N米我><米o>。米o><米年代pace width="5em">
(米o><米n>36米n><米o年代tretchy="false">
)米o>米在h>
总而言之,<一个href="#e35">(35)一个>持有米 约束和<一个href="#e36">(36)一个>持有N 约束。因此,我们必须解决一个线性系统米 +N 方程,年代 传感器参数(权重,<年代trong>h年代trong>P年代ub>)。在矩阵表示法:
C米我><米年代ub><米row>
h米我>米row><米row>
P米text>米row>
=米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
l米text>米row>
,米o><米年代pace width="5em">
(米o><米n>37米n><米o年代tretchy="false">
)米o>米在h>
在哪里<年代trong>C年代trong>是一个(米 +N )×年代 矩阵,<我nl我ne-formula id="inf28">
h米我>米row><米row>
l米text>米row>
≐米o><米年代up>
h米我>米row><米row>
X米text>米row>
T米我>米row>米年代up>
,米o>米td><米tdcolumnalign="center">
h米我>米row><米row>
Y米text>米row>
T米我>米row>米年代up>
T米我>米row>米年代up>
,<年代up>T年代up>表示转置。我们只需要计算<年代trong>C年代trong>一次;注意,它不是frequency-all的函数的值<年代trong>C年代trong>是0或1。
一个解决方案<一个href="#e37">(37)一个>存在的等级<我nl我ne-formula id="inf29">
C米我><米o>,米o>米td><米tdcolumnalign="center">
h米我>米row><米row>
l米text>米row>
=排名 {<年代trong>C年代trong>}。当年代 =锰 ,然后排名 {<年代trong>C年代trong>}=米 +N −1(的行<年代trong>C年代trong>是线性相关的,因为第一次的总和米 行=最后的总和N 行)。中的元素的总和<年代trong>h年代trong>X年代ub>等于元素之和<年代trong>h年代trong>Y年代ub>和= 1来满足无失真约束。由此可见,<我nl我ne-formula id="inf30">
排名米o><米row><米o年代tretchy="false">
{米o><米row><米fenced open="(" close=")">
C米我><米o>,米o>米td><米tdcolumnalign="center">
h米我>米row><米row>
l米text>米row>
}米o>米row><米o>=米o><米我>米米我><米o>+米o><米我>N米我><米o>−米o><米n>1米n>米在h>我nline-formula>。因此,在考虑删除一个传感器时,我们必须确保没有对应的列,排名的<年代trong>C年代trong>仍然是米 +N −1。beamformer权重方法(后被发现<一个href="#B3">Benesty et al ., 2018一个>)。
以下优化问题最大化约束下的绕线<一个href="#e37">(37)一个>:
h米我>米row><米row>
P米text>米row><米row>
W米text>米row>
=米o><米under><米row>
argmin米我>米row><米row>
h米我>米row>米under>
h米我>米row><米row>
H米我>米row>米年代up>
h米我><米年代pace width="0.3333em" class="nbsp">
年代米我><米o>。米o><米我米在hv一个r我一个nt="normal">
t米我><米o>。米o><米年代pace width="0.3333em" class="nbsp">
C米我><米我米在hv一个r我ant="bold">
h米我><米o>=米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
l米text>米row>
。米o><米年代pace width="5em">
(米o><米n>38米n><米o年代tretchy="false">
)米o>米在h>
其给出的解决方案是:
h米我>米row><米row>
P米text>米row><米row>
W米text>米row>
=米o><米年代up>
C米我>米row><米row>
H米我>米row>米年代up>
C米我><米年代up>
C米我>米row><米row>
H米我>米row>米年代up>
−米o><米n>1米n>米row>米年代up>
h米我>米row><米row>
l米text>米row>
,米o><米年代pace width="5em">
(米o><米n>39米n><米o年代tretchy="false">
)米o>米在h>
这也是最小二乘解<一个href="#e37">(37)一个>。以下优化问题最大化约束下的方向性<一个href="#e37">(37)一个>:
h米我>米row><米row>
P米text>米row><米row>
D米text>米row>
=米o><米under><米row>
argmin米我>米row><米row>
h米我>米row>米under>
h米我>米row><米row>
H米我>米row>米年代up>
Γ米我><米我米在hv一个r我ant="bold">
h米我><米年代pace width="0.3333em" class="nbsp">
年代米我><米o>。米o><米我米在hv一个r我一个nt="normal">
t米我><米o>。米o><米年代pace width="0.3333em" class="nbsp">
C米我><米我米在hv一个r我ant="bold">
h米我><米o>=米o><米年代ub><米row><米i mathvariant="bold">
h米我>米row><米row>
l米text>米row>
,米o><米年代pace width="5em">
(米o><米n>40米n><米o年代tretchy="false">
)米o>米在h>
的元素<年代trong>Γ年代trong>是由
Γ米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
年代米我><米o>,米o><米年代up>
年代米我>米row><米row>
′米o>米row>米年代up>
=米o><米我米在hv一个r我一个nt="normal">
年代米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
n米我><米我米在hv一个r我ant="normal">
c米我><米fencedopen="(" close=")">
2米n><米我>π米我><米我>f米我>米row><米row>
v米我>米row>米frac>
d米我>米row><米row>
年代米我><米o>,米o><米年代up>
年代米我>米row><米row>
′米o>米row>米年代up>
;米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>年代米我><米o>,米o><米年代up>
年代米我>米row><米row>
′米o>米row>米年代up>
≤米o><米我>年代米我><米o>;米o><米年代pace width="5em">
(米o><米n>41米n><米o年代tretchy="false">
)米o>米在h>
在哪里d <年代ub>年代 ,年代 ′年代ub>之间的欧几里得距离传感器对应的列年代 和年代 ′的<年代trong>C年代trong>。其给出的解决方案是:
h米我>米row><米row>
P米text>米row><米row>
D米text>米row>
=米o><米年代up>
Γ米我>米row><米row>
−米o><米n>1米n>米row>米年代up>
C米我>米row><米row>
H米我>米row>米年代up>
C米我><米年代up>
Γ米我>米row><米row>
−米o><米n>1米n>米row>米年代up>
C米我>米row><米row>
H米我>米row>米年代up>
−米o><米n>1米n>米row>米年代up>
h米我>米row><米row>
l米text>米row>
。米o><米年代pace width="5em">
(米o><米n>42米n><米o年代tretchy="false">
)米o>米在h>
最大化之间的权衡绕线和指向性可以获得与给出的权重
h米我>米row><米row>
P米text>米row><米row>
α米我>米row>米年代ubsup>
=米o><米年代ub年代up>
Γ米我>米row><米row>
α米我>米row><米row>
−米o><米n>1米n>米row>米年代ub年代up>
C米我>米row><米row>
H米我>米row>米年代up>
C米我><米年代ub年代up>
Γ米我>米row><米row>
α米我>米row><米row>
−米o><米n>1米n>米row>米年代ub年代up>
C米我>米row><米row>
H米我>米row>米年代up>
−米o><米n>1米n>米row>米年代up>
h米我>米row><米row>
l米text>米row>
,米o><米年代pace width="5em">
(米o><米n>43米n><米o年代tretchy="false">
)米o>米在h>
在哪里<年代trong>Γ年代trong>α年代ub>=(1−α )<年代trong>Γ年代trong>+α <年代trong>我年代trong>年代年代ub>为0≤α ≤1。选择α 控制绕线和指向性之间的权衡:一个更大的α 增加了绕线时小α 增加了方向性。我们可以设置α 频率相关。
权衡Beamformers的计算复杂度:
我们假设两个线性阵列。首先,我们分析的计算复杂度最大绕线beamformer<一个href="#e39">(39)一个>。回想一下,<年代trong>C年代trong>不是一个频率的函数。因此,<我nl我ne-formula id="inf31">
C米我>米row><米row>
H米我>米row>米年代up>
(米o><米row><米我米在hvariant="bold">
C米我><米年代up>
C米我>米row><米row>
H米我>米row>米年代up>
)米o>米row>米row>
−米o><米n>1米n>米row>米年代up>
可以计算一次<我nl我ne-formula id="inf32">
O米我><米fencedopen="(" close=")">
年代米我><米fencedopen="(" close=")">
米米我>米row><米row>
2米n>米row>米年代up>
+米o><米年代up>
N米我>米row><米row>
2米n>米row>米年代up>
+米o><米年代up>
米米我>米row><米row>
3米n>米row>米年代up>
+米o><米年代up>
N米我>米row><米row>
3米n>米row>米年代up>
操作。作为年代 ≥米 +N −1一个可行的解决方案,更紧凑的形式的复杂性可以写:<我nl我ne-formula id="inf33">
O米我><米fencedopen="(" close=")">
年代米我><米fencedopen="(" close=")">
米米我>米row><米row>
2米n>米row>米年代up>
+米o><米年代up>
N米我>米row><米row>
2米n>米row>米年代up>
。随后,乘以<年代trong>h年代trong>l年代ub>需要<我nl我ne-formula id="inf34">
年代米我><米fencedopen="(" close=")">
米米我><米o>+米o><米我>N米我>米row>米fenced>
乘法/频率。因此,总计算复杂度<我nl我ne-formula id="inf35">
O米我><米fencedopen="(" close=")">
年代米我><米fencedopen="(" close=")">
米米我><米o>+米o><米我>N米我>米row>米fenced>
米米我><米o>+米o><米我>N米我><米o>+米o><米我>F米我>米row>米fenced>
。
接下来,我们分析的计算复杂度最大的方向性beamformer<一个href="#e42">(42)一个>和权衡beamformer<一个href="#e43">(43)一个>。矩阵乘法和反演需要<我nl我ne-formula id="inf36">
O米我><米fencedopen="(" close=")">
年代米我>米row><米row>
3米n>米row>米年代up>
计算每一个频率。因此,总计算复杂度<我nl我ne-formula id="inf37">
O米我><米fencedopen="(" close=")">
年代米我>米row><米row>
3米n>米row>米年代up>
F米我>米row>米fenced>
。
4线性阵列Constant-Beamwidth Beampattern设计 本节提出了一种方法来获得对称的非均匀线性阵列波束形成的权值,假设传感器位置。然后我们展示如何选择传感器的最佳位置。
考虑一个对称的线性阵列米 = 2l +1传感器。我们使用对称阵列,阵列相位中心频率不变(<一个href="#B47">沃德et al ., 2001一个>),因此,我们可以简单地概括窗口技术。用传感器的位置x <年代ub>l 年代ub>,−l ≤l ≤l 。由于数组的对称性,
x米我>米row><米row>
l米我>米row>米年代ub>
=米o><米n>0米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
l米我><米o>=米o><米n>0米n><米o>,米o>米td>米tr>
x米我>米row><米row>
−米o><米我>l米我>米row>米年代ub>
=米o><米o>−米o><米年代ub><米row><米i>
x米我>米row><米row>
l米我>米row>米年代ub>
,米o><米年代pace width="0.3333em" class="nbsp">
1米n><米o>≤米o><米我>l米我><米o>≤米o><米我>l米我><米o>。米o>米td>米tr>
(米o><米n>44米n><米o年代tretchy="false">
)米o>米在h>
4.1波束形成的权向量的设计 在本节中,我们假设传感器位置。我们的目标是设计beamformer权重最大的DI,同时保持一个常数的波束宽度φ <年代ub>B年代ub>。在形式上,
米米我><米我米在hv一个r我ant="normal">
一个米我><米我米在hv一个r我ant="normal">
x米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
米米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
z米我><米我米在hv一个r我ant="normal">
e米我>米row><米row>
h米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
D米我><米我>我米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
年代米我><米我米在hv一个r我ant="normal">
u米我><米我米在hv一个r我ant="normal">
b米我><米我米在hv一个r我ant="normal">
j米我><米我米在hv一个r我ant="normal">
e米我><米我米在hv一个r我ant="normal">
c米我><米我米在hv一个r我ant="normal">
t米我><米年代pace width="0.3333em" class="nbsp">
t米我><米我米在hv一个r我ant="normal">
o米我>米td>米tr>
h米我>米row><米row>
H米我>米row>米年代up>
f米我>米row>米fenced>
一个米我><米fencedopen="(" close=")">
f米我><米o>,米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
,米o><米年代ub><米row><米i>
θ米我>米row><米row>
d米text>米row>
=米o><米n>1米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
b米我>米row><米row>
φ米我>米row>米年代ub>
=米o><米年代ub><米row><米i>
φ米我>米row><米row>
B米text>米row>
,米o><米年代pace width="0.3333em" class="nbsp">
f米我>米row><米row>
l米text>米row>
≤米o><米我>f米我><米o>≤米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
。米o>米td>米tr>
(米o><米n>45米n><米o年代tretchy="false">
)米o>米在h>
constant-beamwidth波束形成一个简单的方法是限制beamformer权重等于一个窗口的水龙头。调整的窗口有一个参数,因此波束宽度。这个方法被引入长等人(<一个href="#B30">长et al ., 2019一个>这种)。该方法逼近期望的beampattern和计算复杂度之间的妥协。我们推广这个方法对于对称线性数组在使用凯撒窗口。现有的形式相比,我们的推广提高DI。它使维护所需的波束宽度在一个广泛的频率,甚至对这种结果所示。
Kaiser窗取决于一个参数,窗口的形状因子,β 。对于每一个频率,<我nl我ne-formula id="inf38">
β米我><米fencedopen="(" close=")">
f米我>米row>米fenced>
根据需要选择,因此波束宽度。由于逆窗口宽度和beamformer宽度之间的关系(<一个href="#B48">沃德et al ., 1995一个>),增加β 扩大beampattern。可以直观地理解这个逆关系<一个href="#e5">(5)一个>注意到beampattern是(不均匀)的离散傅里叶变换beamformer权重。如果重量是频率无关,beampattern缩小随着频率的增加。因此,更高的值β 在更高的频率用于扩大beampattern波束宽度保持不变。
我们不限制数组是统一的;因此,我们样品连续Kaiser窗在传感器位置。连续Kaiser窗口(<一个href="#B18">凯撒和谢弗,1980一个>)是由
w米我>米row><米row>
K米text>米row>
x米我>米row>米fenced>
=米o><米fencedopen="{" close="">
我米我>米row><米row>
0米n>米row>米年代ub>
β米我><米年代问rt><米row>
1米n><米o>−米o><米年代up>
x米我>米row><米row>
W米我><米o>/米o><米n>2米n>米row>米frac>
2米n>米row>米年代up>
我米我>米row><米row>
0米n>米row>米年代ub>
β米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
x米我>米row>米fenced>
≤米o><米我>W米我><米o>/米o><米n>2米n><米o>,米o>米td>米tr>
0米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
x米我>米row>米fenced>
>米o><米我>W米我><米o>/米o><米n>2米n><米o>,米o>米td>米tr>
(米o><米n>46米n><米o年代tretchy="false">
)米o>米在h>
在哪里<我nl我ne-formula id="inf39">
我米我>米row><米row>
0米n>米row>米年代ub>
x米我>米row>米fenced>
是第一类零阶贝塞尔函数修改,然后呢W 是窗口支持。Kaiser窗的长度R +1)用于数字信号处理的
w米我>米row><米row>
齿龈米text>米row><米row>
R米我>米row>米fenced>
l米我>米row>米fenced>
=米o><米年代ub><米row><米i>
w米我>米row><米row>
K米text>米row>
x米我><米o>=米o><米fr一个c><米row>
W米我>米row><米row>
R米我>米row>米frac>
l米我>米row>米fenced>
=米o><米fencedopen="{" close="">
我米我>米row><米row>
0米n>米row>米年代ub>
β米我><米年代问rt><米row>
1米n><米o>−米o><米年代up>
l米我>米row><米row>
R米我><米o>/米o><米n>2米n>米row>米frac>
2米n>米row>米年代up>
我米我>米row><米row>
0米n>米row>米年代ub>
β米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
l米我>米row>米fenced>
≤米o><米我>R米我><米o>/米o><米n>2米n><米o>,米o>米td>米tr>
0米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
l米我>米row>米fenced>
>米o><米我>R米我><米o>/米o><米n>2米n><米o>,米o>米td>米tr>
(米o><米n>47米n><米o年代tretchy="false">
)米o>米在h>
在哪里R + 1 =米 。我们建议取样连续Kaiser窗在传感器位置,也就是说,
w米我>米row><米row>
拉米text>米row>
l米我>米row>米fenced>
=米o><米年代ub><米row><米i>
w米我>米row><米row>
K米text>米row>
x米我><米o>=米o><米年代ub><米row><米i>
x米我>米row><米row>
l米我>米row>米年代ub>
=米o><米fencedopen="{" close="">
我米我>米row><米row>
0米n>米row>米年代ub>
β米我><米年代问rt><米row>
1米n><米o>−米o><米年代up>
x米我>米row><米row>
l米我>米row>米年代ub>
W米我><米o>/米o><米n>2米n>米row>米frac>
2米n>米row>米年代up>
我米我>米row><米row>
0米n>米row>米年代ub>
β米我>米row>米fenced>
,米o><米年代pace width="0.3333em" class="nbsp">
x米我>米row><米row>
l米我>米row>米年代ub>
≤米o><米我>W米我><米o>/米o><米n>2米n><米o>,米o>米td>米tr>
0米n><米年代pace width="1em">
,米o><米年代pace width="0.3333em" class="nbsp">
x米我>米row><米row>
l米我>米row>米年代ub>
>米o><米我>W米我><米o>/米o><米n>2米n><米o>。米o>米td>米tr>
(米o><米n>48米n><米o年代tretchy="false">
)米o>米在h>
通过选择窗口支持平等的W = 2x <年代ub>l 年代ub>(在完成(<一个href="#B30">长et al ., 2019一个>)),剩下的一个单一的参数调整,β 。我们增加了自由度的问题通过优化参数W 。我们考虑多个值W (每个频率),为每一个W ,调β 所以,波束宽度。然后,每个频率,我们选择W 达到最高DF。减少可能的值的数量W 我们限制窗口支持集<我nl我ne-formula id="inf40">
W米我>米row><米row>
我米我>米row>米年代ub>
我米我><米o>=米o><米n>1米n>米row><米row>
l米我>米row>米年代ubsup>
的大小l 在哪里W <年代ub>我 年代ub>= 2x <年代ub>我 年代ub>。在这里,W <年代ub>我 年代ub>对应于只有2我 +1中心传感器具有非零权重。我们称这些传感器为活跃传感器。在结果中,我们表明,不限制的值的性能改进W 这组是可以忽略不计。因此,我们限制的值W 这将减少找到权重的计算复杂度。请注意,W 每个频率优化,例如,W =W (f )。
反映,一些传感器密集,我们使用梯形集成技术(甲),即,我们把每一个重量乘以<我nl我ne-formula id="inf41">
Δ米我>米row><米row>
l米我>米row>米年代ub>
=米o><米fr一个c><米row>
x米我>米row><米row>
l米我><米o>+米o><米n>1米n>米row>米年代ub>
−米o><米年代ub><米row><米i>
x米我>米row><米row>
l米我><米o>−米o><米n>1米n>米row>米年代ub>
2米n>米row>米frac>
在(<一个href="#B47">沃德et al ., 2001一个>)。我们把最外层的传感器的重量乘以Δ<年代ub>l 年代ub>=Δ<年代ub>−l 年代ub>=x <年代ub>l 年代ub>−x <年代ub>l −1年代ub>。总而言之,
Δ米我>米row><米row>
l米我>米row>米年代ub>
=米o><米fencedopen="{" close="">
x米我>米row><米row>
l米我><米o>+米o><米n>1米n>米row>米年代ub>
−米o><米年代ub><米row><米i>
x米我>米row><米row>
l米我><米o>−米o><米n>1米n>米row>米年代ub>
2米n>米row>米frac>
,米o><米fencedopen="|" close="|">
l米我>米row>米fenced>
<米o><米我>l米我><米o>,米o>米td>米tr>
x米我>米row><米row>
l米我>米row>米年代ub>
−米o><米年代ub><米row><米i>
x米我>米row><米row>
l米我><米o>−米o><米n>1米n>米row>米年代ub>
,米o><米fencedopen="|" close="|">
l米我>米row>米fenced>
=米o><米我>l米我><米o>。米o>米td>米tr>
(米o><米n>49米n><米o年代tretchy="false">
)米o>米在h>
通过使用多在频率使用的所有传感器(W = 2x <年代ub>l 年代ub>中心),传感器,靠的很近没有影响。这使得维护所需的波束宽度在低频率。
最后,我们权重来实现正常化<一个href="#e4">(4)一个>。最后,beamformer重量是:
h米我><米fencedopen="[" close="]">
l米我>米row>米fenced>
=米o><米fr一个c><米row>
Δ米我>米row><米row>
l米我>米row>米年代ub>
w米我>米row><米row>
K米text>米row>
x米我>米row><米row>
l米我>米row>米年代ub>
∑米o>米row><米row>
l米我>米row><米row>
′米o>米row>米年代up>
=米o><米o>−米o><米我>l米我>米row><米row>
l米我>米row>米underover>
Δ米我>米row><米row>
l米我>米row><米row>
′米o>米row>米年代up>
w米我>米row><米row>
K米text>米row>
x米我>米row><米row>
l米我>米row><米row>
′米o>米row>米年代up>
,米o><米年代pace width="0.3333em" class="nbsp">
−米o><米我>l米我><米o>≤米o><米我>l米我><米o>≤米o><米我>l米我><米o>。米o><米年代pace width="5em">
(米o><米n>50米n><米o年代tretchy="false">
)米o>米在h>
每个频率的重量被发现后,迪可以评估。<一个href="#alg1">算法1一个>总结了过程。结果beamformer重量、<我nl我ne-formula id="inf42">
h米我>米row><米o年代tretchy="false">
̂米o>米over>米row>
(米o><米row><米我>f米我>米row><米o年代tretchy="false">
)米o>米row>米在h>
,满足约束<一个href="#e45">(45)一个>。因此,我们继续确定传感器的位置。
计算的复杂性<一个href="#alg1">算法1一个>:
计算的重量<一个href="#e50">(50)一个>需要米 计算产品的乘法<我nl我ne-formula id="inf43">
Δ米我>米row><米row>
l米我>米row>米年代ub>
w米我>米row><米row>
K米text>米row>
x米我>米row><米row>
l米我>米row>米年代ub>
,−l ≤l ≤l 。因为我们从最大的开始β 最初,波束宽度超出预期,因此的绝对值beampattern方位<我nl我ne-formula id="inf44">
φ米我><米o>=米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
+米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
φ米我>米row><米row>
B米text>米row>
大于<我nl我ne-formula id="inf45">
1米n><米o>/米o><米年代问rt><米row><米n>
2米n>米row>米年代问rt>
。因此,要检查是否b <年代ub>φ 年代ub>>φ <年代ub>B年代ub>,它是足以评估的绝对值beampattern方位<我nl我ne-formula id="inf46">
φ米我><米o>=米o><米年代ub><米row><米i>
φ米我>米row><米row>
d米text>米row>
+米o><米fr一个c><米row>
1米n>米row><米row>
2米n>米row>米frac>
φ米我>米row><米row>
B米text>米row>
并检查它是否大于<我nl我ne-formula id="inf47">
1米n><米o>/米o><米年代问rt><米row><米n>
2米n>米row>米年代问rt>
,要求米 乘法。让B 表示数量的窗口形状系数的选择。beampattern是评估F +B 次窗口支持选项,因此评价的<我nl我ne-formula id="inf48">
F米我><米o>+米o><米我>B米我>米row>米fenced>
l米我>米在h>我nline-formula>次了。总之,算法需要的第一部分<我nl我ne-formula id="inf49">
O米我><米fencedopen="(" close=")">
米米我>米row><米row>
2米n>米row>米年代up>
F米我><米o>+米o><米我>B米我>米row>米fenced>
计算。
在第二部分的算法,计算DF要求米 计算重量和乘法米 <年代up>2年代up>+米 为矩阵乘法,乘法/频率和窗口支持选项。因此,需要<我nl我ne-formula id="inf50">
O米我><米fencedopen="(" close=")">
米米我>米row><米row>
3米n>米row>米年代up>
F米我>米row>米fenced>
计算。
总的来说,<一个href="#alg1">算法1一个>需要<我nl我ne-formula id="inf51">
O米我><米fencedopen="(" close=")">
米米我>米row><米row>
2米n>米row>米年代up>
米米我><米我>F米我><米o>+米o><米我>B米我>米row>米fenced>
计算。
4.2对称线性阵列传感器定位 我们的目标是找到在哪里位置传感器,DI最大化而使用重量设计提出了在前面的小节。在形式上,
米米我><米我米在hv一个r我ant="normal">
一个米我><米我米在hv一个r我ant="normal">
x米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
米米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
z米我><米我米在hv一个r我ant="normal">
e米我>米row><米row>
x米我>米row><米row>
l米我>米row>米年代ub>
,米o><米年代pace width="0.3333em" class="nbsp">
−米o><米我>l米我><米o>≤米o><米我>l米我><米o>≤米o><米我>l米我>米row>米under>
D米我><米我>我米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
,米o><米年代pace width="5em">
(米o><米n>51米n><米o年代tretchy="false">
)米o>米在h>
权重达到使用在哪里<一个href="#alg1">算法1一个>;这些重量满足约束<一个href="#e45">(45)一个>。蛮力搜索所有可能的位置是不可行的。所以,我们首先提供了一个迭代算法来确定传感器的位置。由此产生的性能是一个上限DI(使用重量设计提出了在前一节中)。之后,我们提出一个简单的算法,实现竞争力的性能。
4.2.1迭代算法准备 传感器位置的迭代更新,以最大化DI使用Nelder-Mead单纯形算法(<一个href="#B21">Lagarias et al ., 1998一个>)。一组l 变量代表了传感器位置<我nl我ne-formula id="inf52">
x米我>米row><米row>
l米我>米row>米年代ub>
l米我><米o>=米o><米n>1米n>米row><米row>
l米我>米row>米年代ubsup>
。的<我nl我ne-formula id="inf53">
D米我><米我>我米我>米row><米row>
f米我>米row><米row>
l米text>米row>
,米o><米年代ub><米row><米i>
f米我>米row><米row>
H米text>米row>
评估在每个迭代中通过构造对称阵列<一个href="#e44">(44)一个>据,找到beamformer的权重<一个href="#alg1">算法1一个>,最后计算DI<一个href="#e16">(16)一个>。初始化算法多次与不同的起点,以确保全球最大的发现。因此,由此产生的DI是一个上限DI(使用重量设计提出了在前一节中)。
4.2.2价值算法 在本节中,我们概述的步骤来实现一个常数波束宽度最大的频率范围。之后,我们修改的步骤为代价提高DI减少频率范围。
初始化算法与对称线性阵列米 <年代ub>0年代ub>= 2l <年代ub>0年代ub>+ 1 (<米 )传感器。这些传感器是活跃的时,能保持理想的波束宽度的最低频率达到使用最广泛的窗口(β =0),表示这个频率<我nl我ne-formula id="inf54">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
。我们添加一对对称位置的传感器。表示由<我nl我ne-formula id="inf55">
d米我>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
+米o><米n>1米n>米row><米o年代tretchy="false">
)米o>米row>米row>
新的传感器之间的元件间的间距前面最外层的传感器。它是可取的<我nl我ne-formula id="inf56">
d米我>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
+米o><米n>1米n>米row><米o年代tretchy="false">
)米o>米row>米row>
大数组可以支持低频率。然而,如果太大,在更高频率(具体来说,波束宽度<我nl我ne-formula id="inf57">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
)太窄。所以,<我nl我ne-formula id="inf58">
d米我>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
+米o><米n>1米n>米row><米o年代tretchy="false">
)米o>米row>米row>
选择这样<我nl我ne-formula id="inf59">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
下面不活跃和一双新传感器<我nl我ne-formula id="inf60">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
一双新传感器是活跃的。因此,我们增加<我nl我ne-formula id="inf61">
d米我>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
+米o><米n>1米n>米row><米o年代tretchy="false">
)米o>米row>米row>
只要数组可以维护所需的波束宽度<我nl我ne-formula id="inf62">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
。因为<我nl我ne-formula id="inf63">
f米我>米row><米row>
最小值米text>米row><米row>
(米o><米row><米年代ub><米row>
l米我>米row><米row>
0米n>米row>米年代ub>
)米o>米row>米row>
是最高频率,新的传感器是活跃的,获得了理想的波束宽度和最窄窗口:β =β <年代ub>马克斯年代ub>(=10)。添加一对传感器的过程不断重复,直到数组组成米 传感器。
结果显示,旁瓣高值低β :减少DI。因此,我们限制最宽的窗口β =β <年代ub>最小值年代ub>>0。这一修改提高DI但减少了频率范围,因为beampattern不能当一样窄β =0。我们需要初始化米 <年代ub>0年代ub>=5传感器因为三不能达到足够窄波束宽度f <年代ub>H年代ub>=8kHz没有表现出较高的旁瓣。如何选择初始的一个选择米 <年代ub>0年代ub>传感器:齿龈米 <年代ub>0年代ub>的传感器元件间的间距最大化DF的决定f <年代ub>H年代ub>。详细的步骤<一个href="#alg2">算法2一个>。
该算法可用于确定有多少传感器需要保持一个恒定的波束宽度的频带(f <年代ub>l年代ub>,f <年代ub>H年代ub>]。一次f ≤f <年代ub>l年代ub>,我们可以添加新的双传感器完成。一旦找到位置<一个href="#alg2">算法2一个>,be一个米为米er权重了<一个href="#alg1">算法1一个>。计算的复杂性<一个href="#alg2">算法2一个> :
详细的复杂性分析<一个href="#alg1">算法1一个>,<我nl我ne-formula id="inf64">
O米我><米fencedopen="(" close=")">
米米我>米row>米fenced>
操作是否满足检查b <年代ub>φ 年代ub>>φ <年代ub>B年代ub>。让D 表示位置选项的数目。D 来源于最大容许阵列孔径的决议d <年代up>(l )年代up>在每个迭代中增加。beampattern是评估D 次部分找哪里的地方下一对传感器 ,和F 次部分找到最低频率,维护所需的波束宽度 。完全,beampattern评估D +F 次了。因此,<一个href="#alg2">算法2一个>需要<我nl我ne-formula id="inf65">
O米我><米fencedopen="(" close=")">
米米我><米fencedopen="(" close=")">
D米我><米o>+米o><米我>F米我>米row>米fenced>
计算。
5实验结果 在本节中,我们检查我们建议的方法的性能,并对现有方法进行比较。我们首先检查线性阵列和比较现有方法提出了迭代和价值算法。接下来,我们设计平面阵列beamformers使用张量积和权衡。
5.1线性阵列与先前的方法的比较 在本节中,我们检查我们建议的方法的性能和比较之前的方法使用米 =11传感器。所需的半功率波束宽度是任意选择的φ <年代ub>B年代ub>=15°。我们比较设计方法三个方法:儿子(的稀疏阵列设计<一个href="#B39">儿子,2021一个>),长期的齿龈设计等。<一个href="#B30">长et al ., 2019一个>),病房的对数间隔阵列设计等。<一个href="#B47">沃德et al ., 2001一个>)。数组显示在布局<一个href="#F1">图1一个>。
使用前一节中提出的迭代算法,由此产生的传感器位置的对称线性阵列(l ≥1)(3。8,7。9,14.3,29.2,74.8)厘米。生育儿子的数组(<一个href="#B39">儿子,2021一个>),我们将参考beampattern相等
年代米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
n米我><米我米在hv一个r我ant="normal">
c米我><米fencedopen="(" close=")">
φ米我><米o>−米o><米fr一个c><米row>
π米我>米row><米row>
2米n>米row>米frac>
0.2955米n>米row>米frac>
,米o><米年代pace width="5em">
(米o><米n>52米n><米o年代tretchy="false">
)米o>米在h>
15°的半功率波束宽度。我们还检查其他参考beampatterns,然而他们导致类似演出使用<一个href="#e52">(52)一个>。一个初始生成能量最高的传感器选择从齿龈2001传感器和一个元件间的间距为0.1厘米。剩下的2000个传感器被分成五组。从每一组,两个传感器接近该集团被选中的重心。最后,创建稀疏阵列的11所选择的传感器。接下来,权重被发现。mainlobe误差公差是0.006E <年代ub>1年代ub>为φ ∈(82.5°,97.5°)。旁瓣误差公差是0.028E <年代ub>2年代ub>为φ ∈(0°,82.5°)∪(97.5°、180°)。在这里,E <年代ub>1年代ub>=31日E <年代ub>2年代ub>=330是方位角度采样的数量(分辨率为0.5°)。确保波束宽度是根据需要,我们添加了两个约束条件:误差公差是0.0001φ ∈{82.5°,97.5°}。由此产生的稀疏阵列是一个对称的线性阵列与实际重量。长期使用的齿龈权重等。<一个href="#B30">长et al ., 2019一个>)复制使用<一个href="#alg1">算法1一个>,窗口支持仅限于只等于数组的支持,例如,W ≡2x <年代ub>l 年代ub>。的比较,我们使用了一个元件间的间距2.8厘米,因为它导致了最大<我nl我ne-formula id="inf66">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
而实现波束宽度至少15°的频率。复制数组的病房等。<一个href="#B47">沃德et al ., 2001一个>),我们将参考beampattern相等<一个href="#e52">(52)一个>。我们使用对数比3(用问 在(<一个href="#B47">沃德et al ., 2001一个>),因为它导致了DI最高。
最优窗口参数(β 和W 为我们的数组,β 长时间)计算使用<一个href="#alg1">算法1一个>并绘制在<一个href="#F2">图2一个>作为频率的函数。正如预期的那样,对于相同数量的活跃传感器的频率,β 增加随着频率的增加。这是必要的,以抵消的缩小beampattern随着频率的增加。增加β 凯撒窗口缩小,进而扩大了波束宽度由于逆窗口宽度和beamformer宽度之间的关系。我们还观察到有源传感器(2的数量我 +1)减少随着频率的增加。减少活动传感器缩小窗口的数量,产生的效果类似于增加β 。
beamformer权重(计算<一个href="#e50">(50)一个>为我们的数组<一个href="#e47">(47)一个>其次是归一化来实现<一个href="#e4">(4)一个>长时间)所示<一个href="#F3">图3一个>。提议和病房数组内传感器有效活性在低频率由于甲(甲不影响这种由于他们统一元件间的间距)。外传感器给出重量更大,因为他们更大的相邻传感器之间的距离<一个href="#e49">(49)一个>。这使得维护所需的波束宽度在低频有效扩大窗口。这些权重beampatterns构造的根据<一个href="#e5">(5)一个>和说明<一个href="#F3">图3一个>。绕线,波束宽度和DF绘制<一个href="#F4">图4一个>作为频率的函数。儿子的方法性能只是报道在频率绕线至少0分贝(f ≥350年赫兹)。数组保持理想的波束宽度620赫兹,而长期不能使所需的波束宽度低于3.77 kHz。病房的数组没有表现出一个常数波束宽度;病房需要更多的传感器来获得成功。儿子的数组保持理想的波束宽度到350 Hz绕线和DI但低于我们。儿子和我们的数组的绕线低于长和病房由于更少的活跃传感器在每个频率,如图所示<一个href="#F3">图3一一个>。
从DF, DI计算使用<一个href="#e16">(16)一个>。的<我nl我ne-formula id="inf70">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
提出,长,病房数组是7.6,5.9,和5.7 dB,分别;的<我nl我ne-formula id="inf71">
D米我><米我>我米我>米row><米row>
350年米n><米年代pace width="0.3333em" class="nbsp">
H米我><米我米在hv一个r我ant="normal">
z米我><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
拟议的数组是1.4 dB比儿子的。这说明我们提出的设计的优势。
5.2 Beamformer权重设计 在本节中,我们将讨论提出的属性权重的设计。我们报告时提出配置的性能退化的一个设计步骤省略。
乳头改善<我nl我ne-formula id="inf72">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
0.6 dB和启用维护所需的波束宽度较低频率(一个额外的410 Hz),证明该技术的附加值。多给予更大的重量外传感器由于其更大的相邻传感器之间的距离<一个href="#e49">(49)一个>。这使得维护所需的波束宽度在低频有效扩大窗口。
取样连续Kaiser窗显示元素的位置<一个href="#e48">(48)一个>导致温和改善0.05 dB的均匀采样窗口显示<一个href="#e47">(47)一个>。的<我nl我ne-formula id="inf73">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
改进在使用连续的值W 不限制W 集<我nl我ne-formula id="inf74">
2米n><米年代ub><米row><米i>
x米我>米row><米row>
我米我>米row>米年代ub>
我米我><米o>=米o><米n>1米n>米row><米row>
l米我>米row>米年代ubsup>
是微不足道的(0.01 dB),因此,它足以限制窗口支持这个集合。尽管如此,如果我们要限制窗口支持只等于数组的支持,例如,W ≡2x <年代ub>l 年代ub>,将窄波束宽度超过15°以上频率1.63 kHz。这表明设置传感器的一个子集的优点是积极有效地扩大beampattern在更高的频率。
最后,我们检查性能有限的配置是一个数组米 =11传感器齿龈。找到齿龈的最优配置,在前一节中提出的迭代算法优化的一个变量,代表了元件间的间距。由此产生的齿龈的元件间的间距是4.1厘米,<我nl我ne-formula id="inf75">
D米我><米我>我米我>米row><米row>
350年米n><米年代pace width="0.3333em" class="nbsp">
H米我><米我米在hv一个r我ant="normal">
z米我><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
7.2德国铁路股份公司依旧比儿子,长,病房。
5.3价值算法的性能 我们显示的属性提出价值算法确定传感器的位置。我们比较了迭代和价值算法。对于迭代算法,我们在前一节中使用相同的阵列设计。价值的算法,我们初始化的齿龈米 <年代ub>0年代ub>=5传感器元件间的间距为3.4厘米。我们使用这个空间,因为它最大化的DFf <年代ub>H年代ub>=8kHz。我们初始化米 <年代ub>0年代ub>=5传感器因为三不能达到理想的波束宽度8 kHz没有表现出较高的旁瓣。我们设置f <年代ub>H年代ub>= 8 kHz,米 = 11,β <年代ub>马克斯年代ub>=10,φ <年代ub>B年代ub>=15°。
我们构建的两个数组使用价值算法(详细<一个href="#alg2">算法2一个>)。一个数组,我们集β <年代ub>最小值年代ub>=0,和其他数组,我们选择β <年代ub>最小值年代ub>=1。36,最低频率所需的波束宽度是一样的,用迭代算法(620赫兹)。结果显示在传感器位置<一个href="#F5">图5一个>。职位是选择后,参数被发现使用的最佳窗口<一个href="#alg1">算法1一个>并绘制在<一个href="#F6">图6一个>。绕线,波束宽度和DF绘制<一个href="#F7">图7一个>作为频率的函数。
阵列设计β <年代ub>最小值年代ub>=0,我们看到,随着频率减少,另一个传感器被激活之前,尽可能宽的窗口(β =0)的阵列设计β <年代ub>最小值年代ub>=1。36,我们看到最优窗口参数,β ,有时是小于1.36。这是合理的因为β <年代ub>最小值年代ub>是仅用于选择传感器的位置。之后,<一个href="#alg1">算法1一个>可以使用任何的价值吗β 。
我们看到,设置β <年代ub>最小值年代ub>>0导致相比,改进的方向性β <年代ub>最小值年代ub>=0。正如所预料的那样,产生一个数组迭代算法比价值更大的DI算法。然而,<我nl我ne-formula id="inf79">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
获得的价值算法仅为0.2 dB小。这是有益的考虑价值的计算简单算法迭代算法相比,不需要优化工具箱价值算法。
5.4平面阵列 本节比较我们提出技术constant-beamwidth波束形成平面阵列对儿子的设计(<一个href="#B39">儿子,2021一个>使用99个传感器)。数组显示在布局<一个href="#F8">图8一个>。所需的半功率波束宽度是任意选择的θ <年代ub>X年代ub>=15°在XZ平面θ <年代ub>Y年代ub>=30.°在YZ平面。
构建平面阵列,我们首先设计两个线性阵列。为了找到的位置传感器,我们在前一节中使用价值算法:我们设计一个对称的线性数组米 =11传感器固定波束宽度的θ <年代ub>X年代ub>=15°和另一个对称的线性数组N =9传感器固定波束宽度的θ <年代ub>Y年代ub>=30.°。我们初始化的齿龈米 <年代ub>0年代ub>=5传感器元件间的间距为3.4厘米。对于第一个数组,我们集β <年代ub>最小值年代ub>=1。36,第二个数组,我们集β <年代ub>最小值年代ub>=1。00,最低频率beampatterns仍然保持理想的波束宽度都是620赫兹。第二个数组,因为需要较少的传感器所需的波束宽度更宽。我们使用价值算法的迭代算法由于其低计算复杂度而达到性能相提并论。beamformer重量为两个对称的线性数组计算使用<一个href="#alg1">算法1一个>。
基于线性阵列的传感器位置,我们构造一个对称平面阵列共有锰 =99传感器。对于这个提议布局,四个不同beamformer权重向量进行比较。克罗内克积的波束形成,给出了权重的克罗内克积线性阵列的权重<一个href="#e27">(27)一个>。beamformer权重最大的绕线、权衡和最大给出方向性<一个href="#e43">(43)一个>,我们组α 分别∈{1,0.5,0.01}。
生育儿子的数组(<一个href="#B39">儿子,2021一个>),我们将参考beampattern相等
年代米我><米我米在hv一个r我ant="normal">
我米我><米我米在hv一个r我ant="normal">
n米我><米我米在hv一个r我ant="normal">
c米我><米fencedopen="(" close=")">
θ米我>米row><米row>
−米o><米n>0.1477米n><米o>米o><米我>因为米我><米fencedopen="(" close=")">
2米n><米我>φ米我>米row>米fenced>
+米o><米n>0.4432米n>米row>米frac>
,米o><米年代pace width="5em">
(米o><米n>53米n><米o年代tretchy="false">
)米o>米在h>
有半功率波束宽度的XZ平面上15°和30°在YZ平面。从一个统一的平面阵列201×201传感器和一个元件间的间距为1厘米,最初30传感器生成能量最高的选择。剩下的传感器被分成23组。从每一组,三个传感器接近该集团被选中的重心。最后,稀疏阵列成立与99年选择传感器。接下来,权重被发现。mainlobe误差公差是0.006E <年代ub>1年代ub>为θ ∈(0°15°)。旁瓣误差公差是0.028E <年代ub>2年代ub>为θ ∈(15°、90°)。在这里,E <年代ub>1年代ub>=31日×181E <年代ub>2年代ub>=150×181的数量角度采样(方位角和仰角决议2°和0.5°,分别)。确保波束宽度是根据需要,我们添加了四个约束:误差公差是0.0001 (φ ,θ )∈{(0°、7.5°),(90°15°),(180°,7.5°),(270°15°)}。由此产生的稀疏阵列是不对称的,复值权重。
这些权重beampatterns构造的根据<一个href="#e5">(5)一个>说明了几个频率<一个href="#F9">图9一个>。XZ的波束宽度和YZ平面的提议和儿子的beamformers绘制<一个href="#F10">图10一个>作为频率的函数。拟议中的beamformers的波束宽度一样的线性阵列的波束宽度他们建造的。克罗内克beamformer,这是由于财产<一个href="#e28">(28)一个>和其他beamformers,这是因为约束<一个href="#e37">(37)一个>。因此,拟议的技术维护所需的波束宽度620 Hz。儿子的数组保持理想的波束宽度在更低的频率。我们看到在<一个href="#F10">图10一个>儿子的波束宽度窄比期望的频率。这是因为beampattern跌破<我nl我ne-formula id="inf80">
1米n><米o>/米o><米年代问rt><米row><米n>
2米n>米row>米年代问rt>
然后回到达到顶峰<我nl我ne-formula id="inf81">
1米n><米o>/米o><米年代问rt><米row><米n>
2米n>米row>米年代问rt>
在<我nl我ne-formula id="inf82">
1米n>米row><米row>
2米n>米row>米frac>
θ米我>米row><米row>
Y米text>米row>
作为约束。
绕线和DF绘制<一个href="#F11">图11一个>作为频率的函数。儿子的表现只是报道在频率绕线至少0分贝(f ≥890年赫兹)。对于大多数的频率,我们设计的绕线技术低于儿子的。这是预期的,因为儿子的目标是最大化绕线。我们也看到,大α 是该技术,绕线越大。
从DF, DI计算使用<一个href="#e16">(16)一个>。的<我nl我ne-formula id="inf85">
D米我><米我>我米我>米row><米row>
890年米n><米年代pace width="0.3333em" class="nbsp">
H米我><米我米在hv一个r我ant="normal">
z米我><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
拟议的技术高于通过儿子。这说明该技术的优势。我们也看到,较小的α 是该技术,DI越大。重要的是要注意,权衡beamformer绕线和迪克罗内克beamformer。克罗内克beamformer的性能较低,因为它没有解决的优化问题获得理想的波束宽度。
我们的技术的主要优势在儿子的计算复杂度。我们报告的运行时间英特尔<年代up>ⓒ年代up>™核心i7 - 8550 u与8 GB RAM和CPU在MATLAB中实现的代码。权衡beamformer, 3 s才找到两个线性阵列的传感器位置使用价值算法,17世纪计算权重<一个href="#alg1">算法1一个>(该决议β 0。01),33个年代计算beamformer权重(801频率)。总计:小于1分钟。另一方面,对于儿子,279年代才选择的位置传感器和286年代beamformer权重计算一个频率。为所有801个频率计算,儿子的技术比我们花了1000多倍的时间。
最后,我们检查时的性能限制了数组配置一个11×11平面阵列均匀间距。这种我们构造两个(每个都有米 =11传感器)的元件间的间距4.2厘米的传感器的位置。我们使用这个空间,因为它导致最高<我nl我ne-formula id="inf86">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
同时保持的宽波束宽度至少为所有频率。我们使用<一个href="#alg1">算法1一个>找到这种重量的:对于一个理想的波束宽度θ <年代ub>X年代ub>=15°和其他所需的波束宽度θ <年代ub>Y年代ub>=30.°。从这两个这种,我们构造一个平面阵列(4.2厘米)的元件间的间距的权重beamformer的权衡<一个href="#e43">(43)一个>与α =0。01最大化的方向性。它的<我nl我ne-formula id="inf87">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
11.2 dB-lower至少1.3 dB比我们提出非均匀间隔的平面阵列beamformers,尽管它有121个传感器。这类似于我们看到线性阵列:非均匀间隔的数组实现更大的DI。
5.5平面阵列与失踪的传感器 在本节中,我们模拟一个平面阵列传感器只在一些网格的位置。我们使用相同的平面阵列设计一节。然而,我们只传感器年代 =39的锰 =99网格的位置。此外,我们选择把传感器的方式很像一个明星。选定的网格位置显示在<一个href="#F8">图8一个>和称为明星数组。因为只有使用网格位置的一个子集,在克罗内克积beamformer是无效的。因此,我们只比较三种不同beamformers: beamformers给出的<一个href="#e43">(43)一个>与α ∈{1,0.5,0.01},最大绕线,权衡,分别和最大的方向性。beamformers的重量、beampatterns构造的根据<一个href="#e5">(5)一个>和几个频率<一个href="#F12">图12一个>。显示了波束宽度<一个href="#F10">图10一个>。尽管我们只放置在网格的一个子集位置传感器,beamformers仍然保持理想的波束宽度620 Hz。绕线和DF绘制<一个href="#F13">图13一个>作为频率的函数。正如所料,大α ,高是绕线和低DI。的<我nl我ne-formula id="inf88">
D米我><米我>我米我>米row><米row>
890年米n><米年代pace width="0.3333em" class="nbsp">
H米我><米我米在hv一个r我ant="normal">
z米我><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
的年代 =39传感器明星数组仍高于99年通过儿子的稀疏阵列传感器。这说明该技术的优势。
我们比较明星阵列的性能与均匀间隔7×7平面阵列。我们找到的位置传感器,构建两个这种(每个都有米 =7传感器元件间的间距为4.3厘米。我们使用这个空间,因为它导致最高<我nl我ne-formula id="inf91">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
同时保持至少一样宽的波束宽度所需的频率。我们使用<一个href="#alg1">算法1一个>找到这种重量的:对于一个理想的波束宽度θ <年代ub>X年代ub>=15°,和其他所需的波束宽度θ <年代ub>Y年代ub>=30.°。从这两个这种,我们构建一个统一的平面阵列(4.3厘米)的元件间的间距的权重beamformer的权衡<一个href="#e43">(43)一个>与α =0。01最大化的方向性。它的<我nl我ne-formula id="inf92">
D米我><米我>我米我>米row><米row>
0米n><米o>,米o><米年代pace width="0.3333em" class="nbsp">
8米n><米年代pace width="0.3333em" class="nbsp">
k米我><米我米在hv一个r我ant="normal">
H米我><米我米在hv一个r我ant="normal">
z米我>米row>米fenced>
8.7 dB-lower至少比我们的非均匀间隔的星阵beamformers 0.4 dB,尽管它有49个传感器。这又一次显示了使用非均匀间隔设计的价值。
6结论 我们提出了在XZ平面beamformers产生不同的波束宽度和YZ平面设计constant-beamwidth线性阵列。平面波束宽度等于所有频率的线性阵列波束宽度。第一种技术包括利用克罗内克积性质,不需要找到权重矩阵求逆。我们的第二个技术提供了一个封闭的解决方案,允许一个绕线和指向性之间的权衡。第二个方法是适用的,即使只有一个子集的传感器使用。这还需要进一步的研究来选择最优的子集传感器。
平面beamformer作为构建块,我们提出一个方法来找到最优传感器安置在一个对称的非均匀线性阵列和最佳beamformer DI最大化。设计包括取样连续Kaiser窗,不同窗口支持,使用甲。取样连续Kaiser窗在传感器位置适当考虑到非均匀对称阵列元件间的间距。不同的窗口支持和窗口形状系数使更多的控制窗口的有效宽度。使用多以低频率允许维护所需的波束宽度。实验结果表明,该技术产生高指向性和涉及比最先进的计算复杂度低。另外,由此产生的对称的非均匀阵列比uniformly-spaced数组。未来的工作将集中在推广该方法的三维数组和寻址的问题constant-beamwidth转向与非均匀阵列。
数据可用性声明 最初的贡献提出了研究中都包含在本文/辅料,可以针对相应的作者进一步询问。
作者的贡献 房颤导致方法、调查、正式的分析软件,验证和原创作品。集成电路导致概念化、方法、监督和writing-review和编辑。
资金 这项工作是Pazy研究基金会的支持。
的利益冲突 作者声明,这项研究是在没有进行任何商业或财务关系可能被视为一个潜在的利益冲突。
出版商的注意 本文表达的所有索赔仅代表作者,不一定代表的附属组织,或出版商、编辑和审稿人。任何产品,可以评估在这篇文章中,或声称,可能是由其制造商,不保证或认可的出版商。
确认 作者要感谢茎音频提供设备和技术指导。
引用
Aroudi,。,Br一个un,年代。(2021). “Dbnet: Doa-Driven Beamforming Network for End-To-End Reverberant Sound Source Separation,” in<年代pan class="conf-name">ICASSP 2021 - 2021年IEEE国际会议音响、演讲和信号处理(ICASSP年代pan>(<年代pan class="publisher-name">IEEE年代pan>),211 - 215。doi: 10.1109 / icassp39728.2021.9414187
CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=Dbnet:+Doa-Driven+Beamforming+Network+for+End-To-End+Reverberant+Sound+Source+Separation&btnG=">谷歌学术搜索一个>
Benesty, J。科恩,I。,Chen,J。(2019)。数组处理:克罗内克积波束形成 ,18岁。<年代pan class="publisher-name">施普林格年代pan>。
谷歌学术搜索一个>
Benesty, J。科恩,I。,Chen,J。(2018)。信号增强和阵列信号处理的基础 。<年代pan class="publisher-name">威利在线图书馆年代pan>。
谷歌学术搜索一个>
Buchris Y。阿玛,。,Benesty, J。,科恩,I。(2018)。非相干合成Frequency-Invariant稀疏阵列的波束形成。我EEE/ACM反式。音频,演讲,朗。过程。 27日,482 - 495。
谷歌学术搜索一个>
Chan)。,Chen,H。(2006)。 Uniform Concentric Circular Arrays with Frequency-Invariant Characteristics—Theory, Design, Adaptive Beamforming and Doa Estimation.我EEE反式。信号的过程。 55岁,165 - 177。
谷歌学术搜索一个>
科恩,I。,Benesty, J。,G一个nnot,年代。(2009).现代通信中语音处理:挑战和视角 3。<年代pan class="publisher-name">施普林格科学与商业媒体年代pan>。
谷歌学术搜索一个>
Gazor, S。,Gren我er,Y。(1995). Criteria for Positioning of Sensors for a Microphone Array.我EEE反式。演讲音频的过程。 3,294- 303。
谷歌学术搜索一个>
库恩,e . V。,P我tz,C。A., Matsuo, M. V., Bakri, K. J., Seara, R., and Benesty, J. (2021). A Kronecker Product Clms Algorithm for Adaptive Beamforming.数字信号。的过程。 111年,102968年。doi: 10.1016 / j.dsp.2021.102968
CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=A+Kronecker+Product+Clms+Algorithm+for+Adaptive+Beamforming&btnG=">谷歌学术搜索一个>
李。,杨,X。,N我ng, L., Long, T., and Sarkar, T. K. (2017). “Broadband Constant Beamwidth Beamforming for Suppressing Mainlobe and Sidelobe Interferences,” in<年代pan class="conf-name">2017年IEEE雷达会议(RadarConf)年代pan>(<年代pan class="publisher-name">IEEE年代pan>),1041 - 1045。doi: 10.1109 / radar.2017.7944358
CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=Broadband+Constant+Beamwidth+Beamforming+for+Suppressing+Mainlobe+and+Sidelobe+Interferences&btnG=">谷歌学术搜索一个>
刘,R。,Cu我,G。,lu,问。, Yu, X., Feng, L., and Zhu, J. (2021). “Constant Beamwidth Receiving Beamforming Based on Template Matching,” in<年代pan class="conf-name">2021年IEEE雷达会议(RadarConf21)年代pan>(<年代pan class="publisher-name">IEEE年代pan>),1 - 5。doi: 10.1109 / radarconf2147009.2021.9455192
CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=Constant+Beamwidth+Receiving+Beamforming+Based+on+Template+Matching&btnG=">谷歌学术搜索一个>
刘,W。,维斯,S。(2004)。“新类的宽带阵列频率不变的光束模式,”<年代pan class="conf-name">2004年IEEE国际会议音响、演讲和信号处理(IEEE)年代pan>。2 - 185.2
谷歌学术搜索一个>
刘,W。,维斯,S。(2010)。宽带波束形成:概念和技术 ,17岁。<年代pan class="publisher-name">约翰威利& Sons年代pan>。
谷歌学术搜索一个>
刘,Y。,Cheng, J., Xu, K. D., Yang, S., Liu, Q. H., and Guo, Y. J. (2018). Reducing the Number of Elements in the Synthesis of a Broadband Linear Array with Multiple Simultaneous Frequency-Invariant Beam Patterns.我EEE反式。天线Propagat。 66年,5838 - 5848。doi: 10.1109 / tap.2018.2862361
CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=Reducing+the+Number+of+Elements+in+the+Synthesis+of+a+Broadband+Linear+Array+with+Multiple+Simultaneous+Frequency-Invariant+Beam+Patterns&btnG=">谷歌学术搜索一个>
长,T。,科恩,I。,Berdugo, B., Yang, Y., and Chen, J. (2019). Window-based Constant Beamwidth Beamformer.传感器 19日,2091年。doi: 10.3390 / s19092091
《公共医学图书馆摘要》一个>|<一个href="https://doi.org/10.3390/s19092091">CrossRef全文一个>|<一个href="https://scholar.google.com/scholar?hl=en&as_sdt=0%2C5&q=Window-based+Constant+Beamwidth+Beamformer&btnG=">谷歌学术搜索一个>
Ramezanpour, P。,Aghababaie, M., Mosavi, M., and M de Andrés, D. (2022). Multi-stage Beamforming Using Dnns.伊朗j .电工实习。电子。Eng。 ,2266年。
谷歌学术搜索一个>
范树,h·l . (2004)。最佳阵处理:第四部分的检测、估计、调制理论 。<年代pan class="publisher-name">约翰威利& Sons年代pan>。
谷歌学术搜索一个>
王,X。,Benesty, J。、黄G。陈,J。,科恩,I。(2019)。克罗内克积设计与长方体beamformers麦克风阵列(Universitatsbibliothek der亚琛工业) 。
谷歌学术搜索一个>
王,X。,黄G。,Benesty, J。陈,J。,科恩,I。(2020)。到达时差估计基于克罗内克积分解。我EEE信号。的过程。列托人。 28日,51-55。
谷歌学术搜索一个>
燕,S。,米一个,Y。(2005)。“设计的冷杉Beamformer频率不变的模式通过共同优化空间和频率响应,”<年代pan class="conf-name">程序。(ICASSP 05)。IEEE国际会议音响、演讲和信号处理,2005年代pan>(<年代pan class="publisher-name">IEEE年代pan>)。4所示。4 - 789。
谷歌学术搜索一个>
燕,S。,年代un,H。,米一个,X。, Svensson, U. P., and Hou, C. (2010). Time-domain Implementation of Broadband Beamformer in Spherical Harmonics Domain.我EEE反式。音频,演讲,朗。过程。 19日,1221 - 1230。
谷歌学术搜索一个>
杨,W。,黄G。,Chen,J。,Benesty, J。科恩,I。,Kellermann, W. (2020). Robust Dereverberation with Kronecker Product Based Multichannel Linear Prediction.我EEE信号。的过程。列托人。 28日,101 - 105。
谷歌学术搜索一个>
赵,Y。,刘,W。,l一个ngley, R. (2009). “A Least Squares Approach to the Design of Frequency Invariant Beamformers,” in<年代pan class="conf-name">2009年欧洲17信号处理的会议年代pan>(<年代pan class="publisher-name">IEEE年代pan>),844 - 848。
谷歌学术搜索一个>
赵,Y。,刘,W。(2012)。健壮的固定频率不变的Beamformer设计主题Norm-Bounded错误。我EEE信号。的过程。列托人。 20.岁,169 - 172。
谷歌学术搜索一个>
脚注 1一个>年代up> https://github.com/Ariel12321/Constant_Beamwidth_Beamforming_Nonuniform一个>
关键词:年代pan>宽带beamformer constant-beamwidth beamformer、麦克风阵列,克罗内克积,指向性因子,白噪声增益,Kaiser窗,非均匀阵列
引用:年代pan>弗兰克·科恩和我(2022)Constant-Beamwidth克罗内克积与非均匀平面阵列波束形成。前面。Proc团体。 2:829日463。doi: 10.3389 / frsip.2022.829463
收到:年代pan>2021年12月05;<年代pan>接受:年代pan>2022年4月01; <年代pan>发表:年代pan>2022年5月11日。
人也看了
下载年代pan>