课程设计(论文)任务书
专 业
电子信息工程
班 级
0802503
学 生
080250303
指导教师
刘帅、王军、周洪娟
题 目
自适应波束形成算法仿真
子 题
设计时间
2011 年 11 月 28 日 至 2011 年 12 月 9 日 共 2 周
设计要求
设计(论文)的任务和基本要求,包括设计任务、查阅文献、方案设
计、说明书(计算、图纸、撰写内容及规范等)、工作量等内容。
1. 设计任务:实现自适应波束形成算法,并考察算法性能。
2. 参考文献:《现代信号处理》,《雷达原理》, Matlab 相关书籍等
3. 方案设计:
1) 由均匀线阵形式,确定阵列的导向矢量;
2) 由阵列导向矢量,对接收信号进行建模仿真;
3) 根据约束条件,实现自适应波束形成算法;
4) 考察算法性能与信噪比,采样率,观测时间等参数的关系。
4. 说明书:设计报告内容要求包括:信号模型的建立、分析过程;建模
源程序;仿真结果分析。
本专业综合设计要求在两周内完成。
指导教师签字: 系(教研室)主任签字:
年 月 日
前言
近十年,自适应波束形成算法在通信雷达、声纳、生物医学工程等科技领域
中到了极为广泛的应用。在实际应用中,如果信号源、天线阵列出现误差,传统
的自适应波束形成算法性能将会下降。但是对于稳健的自适应波束形成算法。环
境及天线阵列的误差和不确定性是必须要考虑的关键问题。这里对稳健的自适应
波束形成算法的研究现状与发展动态进行了较为详细的评述。
传统的算法在设计波束形成器时,都是假定在其训练数据中不含有期望信号。
这时的波束形成器对于阵列响应误差和有限次快拍数据的稳健性非常好。但
是在许多情况下,干扰和噪声信号的观测数据不可避免地要被期望信号所污染。
而即使在理想情况下,这时精确的知道期望信号的方向向量,在有限次快拍训练
数据中含有的 响应之间的不匹配。自适应阵列信号处理对于这类误差的反应非
常敏感。因为这时期望信号分量会被当作干扰抑制掉。阵列校正误差及观测方向
的偏差等现象在实际中经常出现。
引起自适应波束形成算法性能降低的原因除了阵列响应误差外,由于信道的
不稳定、干扰和天线的移动、天线的转动等所造成的波束形成器训练数据的不稳
定性也是一个主要原因。首先,即使在训练数据 中不包含期望信号,这种不稳
定性也会限制自适应波束形成算法的性能;进而,当在训练数据中含有期望信号
时,波束形成器性能会进一步严重降低;最后,如果出现快速移动的干扰。波束
形成器的性能将会崩亏。因为这时自适应权向量的收敛速度跟不上干扰的变化,
不能有效的对干扰形成零陷。
前言....................................................................................................................................................2
一.设计原理....................................................................................................................................3
二.设计过程....................................................................................................................................6
(1)使用 MVDR 算法............................................................................................................6
(2)EC 算法............................................................................................................................9
(3)GUI 界面设计................................................................................................................12
三.设计总结..................................................................................................................................12
四.参考文献..................................................................................................................................13
一.设计原理
传统算法,波束形成器的输出可表示为 ( 1 ) ,这里复向量
x( k )是在观测阵元上第 k 次快拍的采样值,它是 Mx l 维的复向量;W 是波束形
成器 Mx l 维的权向量设阵元数目为 M,( - ) n 代表共轭转置,其 中向量 x ( k )
可表示为( k ) = s ( k ) + i ( k ) + n ( k )。这里 S ( k ) ,i ( k ) ,n ( k )分别代表期
望信号、干扰和噪声,它们都是统计独立的分量。当电波从点辐射源以球面波向
外传播,只要离辐射源足够远,在接收的局部区域,球面波就可以近似为平面波,
通信信号的传播一般都满足这一远场条件。这时有 其中 是信号
的方向向量。
则最佳权向量可以通过使信干噪比( S I N R ) 的最大化来获得:
wRw
wRw
ni
H
s
H
�
�SINR
(2)
这里 , ,它们分别是 MxM
维的信号协方差矩阵、干扰加噪声协方差矩阵,其中 E { .}代表数学期望。为
了求得式( 2 )中 S I N R 的最大值,可以通过在约束条件 的约束下,求
满足
mi n (3 )
的权向量,也就是使输出信噪比的功率最小。
当 时,式( 3 )的约束条件可简化为
这 种 方 法 称 为 最 小 方 差 无 畸 变 响 应 ( MVDR : minimum variance
distortionless response ) 波束形成。 式( 3 ) 的结果为
( 4 )
式中, 在实际计算中,式(4)中的协方差矩阵 用有限次
快拍数据得到的估计值 代替
(5)
式中 X= 【x ( 1 ), x ( 2 ),…,x ( N )】是 Mx N 维矩阵,N 为快拍数。
稳健的自适应波束形成算法, 目前最流行的算法,是使用对角加载
( diagonal loading ) 技术,它在小快拍数据和存在阵列信号响应误差的情况下。
仍然具有较好的稳健性。它的核心思想是:通过在目标函数式( 3 ) 中增加一个
� � � �
kxw
H
�ky
� � � �
ss
aksks �
s
a
� � � �
� �
kSkSER
H
sss
�
� � � �
� �
� � � �
� �
� �
H
ni
knkiknkiER ���
�
1�wRw
s
H
wRw
ni
H
�
� �
1�
s
Rrank
1�
s
H
aw
sniopt
aRw
1�
�
�
�
� �
1
1
�
�
�
�
sni
H
s
aRa
�
ni
R
�
�
R
� � � �
HH
N
n
XX
N
nxnx
N
R
11
1
���
�
�
二次补偿项来调整权向量。可以表示为
s . t .
( 6 )
式中 为加载因子。式( 6 ) 的解为
( 7 )
式中 I 为单位矩阵。
对角加载技术能够提高自适应波束形成算法的性能,但是这种方法的
缺陷是对于
参数 的选择不可靠。而且也不易获得。另外一个比较流行的稳健自适应
波束形成算 法是基于特征空间的波束形成算法。这种方
法仅适用于点信号源。其核心思想是利用方 向向量在采样信号与干扰信号子空
间上的投影来代替信号的方向向量。
式( 5 ) 的特征分解为
( 8 )
式中 M× ( L + 1 ) 维矩阵 包含 中的 L + 1 个信号与干扰特征向
量, 为其 特征值构成的 ( L + I ) x ( L + 1 )维对角矩阵;M× ( M— L— 1 )维
矩阵 包含 中的 M—L 一 1 个噪声特征向量, 为其特征值构成的对角矩阵。
假设干扰信号数为 L ,则基于特征空间的波束形成算法的权向量
可表示为
( 9 )
式 , 中, 是在信号与干扰子
空间上 的正交投影。
如果环境条件匹配, 基于特征空间的波束形成算法的性能会大幅度的
提高。但是这种方法的致命缺陷是要求数据必须是基于低秩的固定模型且干
扰源数必须精确知道,而且这种方法只适用于高信噪比的情 况下,当信噪
比较低时,由于可能会发生子空间交换 而导致投影矩阵失效。
针对上述算法中的缺陷,已经提出了许多改进的算法。为了解决权向量
wwwRw
HH
�
�
�
min
1�wRw
S
H
�
sDL
aIRw
1�
�
�
�
�
�
�
�
��
�
�
HH
GTGEER
�������
���
�
E
�
R
�
�
�
G
�
R
�
T
vRw
eig
1�
�
�
s
E
aPv
�
�
HHH
E
EEEEEEP
���
�
���
�
�
�
�
�
�
�
�
�
�
�
1