发布时间:2023-06-12 09:00
2021年2月27日发(作者:融合英文)
信息与通信工程学院
阵列信号处理实验报告
<
p>
(自适应波束形成
Matlab
仿真)<
p>
学
号:
XXXXXX
专
业:
XXXXXX
学生姓名:
XXX
任课教师:
XXX
2015
年
X
月
题目
:自适应波束形成
Matlab
仿真
1.
算法简述:
自适应波束形成,
源于自适应天线的一个概念。
接收端的信号处理,
可以通
过将各阵元输出进行加权求和,将天线阵列波束“导向”到一个方向上,对期
望
信号得到最大输出功率的导向位置即给出波达方向估计。
<
/p>
波束形成算法是在一定准则下综合个输入信息来计算最优权值的数学方法,
线性约束最小方差准则(
LCMV
)是最重要、最常用
的方法之一。
LCMV
是对有
用信号形
式和来向完全已知,
在某种约束条件下使阵列输出的方差最小。
该准则
属于广义约束,缺点是需要知道期望分量的波达方向。准则的代价函数为
J
(
w
)
?
w
H
Rw
,约束条件为
w
H
a
(
?
)
?
f
;最佳解为
w
H
R
?
1
c
[
c
H
R
<
br>?
1
c
]
?
1
f
。
2.
波束形成原理
以一维
M
元等距离线阵为例,如图
1
所示,设空间信号为窄带信号,每个
通道用一个附加权值系数来调整该通道的幅度和相位
。
x
1
(<
/p>
n
)
x
1
(
n
)<
p>
x
1
(
n
)
……..
w
1
(
?
)
w
1
(
?
)
w
1
(
?
)
y
p>
(
n
)
图
1
波束形成算法结构图
这时阵列的输出可以表示为:
y
p>
(
t
)
?
w
i
(
?
)
x
i
(
<
br>t
)
?
i
*
?
1
M
如果采用矢量来表示各阵元输出及加权系数,即
x
(
t
)
<
br>?
[
x
1
(
t
)
x
2
(
t
)
x
p>
M
(
t
)]
T
w
(<
p>
?
)
?
[
w
1
(
?
)
w
2
(
?
)
w
M
(<
/p>
?
)]
T
那么,阵列的输出也可以用矢量表示为
y
(
t
)
?
w
H
(
?
p>
)
x
(
t
)
为了在某一方向
?
上补偿各阵元之间的时延以形成一个主瓣,常规波束形成
器在期
望方向上的加权矢量可以构成为
w
(
?
)
?
[
p>
1
e
?
j
??
e
?
j<
p>
(
M
?
1)
??
]
T
p>
观察此加权矢量,
发现若空间只有一个来自方向
?
的信号,
其方向矢量
?
(
?
)
的表示形式与此
权值矢量相同。则有
y
(<
p>
t
)
?
w
H
(
?
)
x
(
t
)
?
?
H
(
?<
/p>
)
x
(
t
)
这时常规波束形成器的输出功率可以表示为
<
br>P
CBF
(
?
<
br>)
?
E
[
y
(
t
)
2
]
?
w
H
p>
(
?
)
Rw
(
?
)
?<
p>
?
H
(
?
)
R
?
(
?
)
式中矩阵
<
br>R
为阵列输出
x
(
t
)
的协方差矩阵。
3.
实验内容与结果:
实验使用均匀线阵
,
阵元间距为信号波长的一半,
输入信号为
1
个
BPSK
信
<
br>号,
2
个非相干的单频干扰,设置载波频率
10MHz
、采样频率
50MHz
<
br>、快拍数
300
、信噪比
-25
dB
、信干比
-90dB
、信号方位角
0
、干扰方位角
?
40
和
50
,分析
阵元数分别在
3
、
6
、
9
和
12<
p>
时波束图的变化。实验结果见图
1
。
图
1
不同阵元数情况下的波束图
4.
仿真分析
由上图可知,在已知实验条
件下,均匀线阵阵元数不小于信号干扰总数时,
应用
LCMV<
/p>
方法得到的波束图在信噪比很低的情况下可以有效抑制非相干干
扰
;对于信号方向虽然波束主瓣可以较准确地指向信号来向,但阵元数为
4
时,
主瓣波束宽度较大,
旁瓣增益较大;
当阵元数增加,
主瓣波束逐渐变窄并且能精
准地指向
信号来向,同时旁瓣数目增加、增益逐渐变小。因此,增加阵元的数目
可以有效地提高波
束形成精度,
同时主瓣变窄,
旁瓣变小,
使得输出信噪比有效
提高。
5.
程序
clear all;
close all;
clc;
%---------------------------------
%
信号参数设置
Signal_No = 1;
Interference_No =
2;%
小于五
S_No =
Signal_No + Interference_No;
%
信号总数
sensor_No = 12;
%
阵元数
azimuth = [0*pi/180 -40*pi/180
50*pi/180 30*pi/180 20*pi/180];
%
方位角
Fs = 50e6;
Fc = 10e6;
F1 = [8e6 9e6
10e6 11e6];
%
干扰功率
RB = 1e6;
M = 2;
wavelength = 3e8/Fc;
d =
wavelength/2;
K = 300;
%
快拍数
Ps = 0;
%dB
SNR = -15;
SIR = -90;
%
二进制