发布时间 : 星期六 文章阵列信号处理—music、Capon更新完毕开始阅读6e191969a4e9856a561252d380eb6294dd882225
宽带信号中的三种二维平面阵DOA估计
宽带信号中的三种二维平面阵DOA估计
一． 背景
目前关于阵列窄带信号的高分辨算法已比较成熟，但是随着信号处理技术的发展，信号环境日趋复杂，信号形式多样，信号密度日渐增大，窄带阵列探测系统的确定逐渐显示出来。
由于宽带信号具有目标回波携带的信息量大，有利于目标探测、参量估计和目标特征提取等特点，在有源探测系统中越来越多地得到应用。而在无源探测系统中，利用目标辐射的宽带连续谱进行目标检测是有效发现目标的一种重要手段。
ISM方法把宽带信号在频域分解为J个窄带分量，然后在每一个子带上直接进行窄带处理。因为信号为调频信号，所以信号在时域的分段实际上就是频域的分段。将信号分解为窄带信号后，我们就可以利用窄带算法进行处理，最后将各个结果进行加权综合，即可得到最终的结果。 二维DOA估计是阵列信号处理中的重要内容，通过二维DOA估计可以得到信号源在平面中的角度信息。一般采用L型、面阵和平行阵或矢量传感器实现二维参数的估计，多数有效的二维DOA估计算法是在一维DOA估计的基础上，直接针对空间二维谱提出的，如二维MUSIC算法以及二维CAPON算法等。这两种算法可以产生渐进无偏估计，但要在二维参数空间搜索谱峰，计算量相当大。而采用二维ROOT MUSIC算法可以减小计算量，但是需要付出精度下降的代价。 本次报告将结合宽带信号和二维DOA估计算法，进行相关的算法介绍和仿真。
二． 算法介绍
1. 接收信号模型：
图 1 平面阵列示意图
如图1所示，设平面阵元数为M×N，信源数为K。信源的波达方向为(?1,?1),?,(?k,?k)，
第i个阵元与参考阵元之间的波程差为：
??2?(xicos?sin??yisin?sin??zicos?)/?
设子阵1沿x轴的方向矩阵为Ax，而子阵2的每个阵元相对于参考阵元的波程差就等于子阵1的阵元的波程差加上2?dsin?sin?/?，所以接收信号为
?AyD1(Ax)???AD(A)y2x?S?N X?????????AyDM?1(Ax)??协方差矩阵为
R?XXH?EsDsEsH?EnDnEnH
其中，Ds代表由最大的K个特征值构成的一个K×K对角阵，Dn代表由MN-k个较小的特征值构成的对角矩阵, Es和En分别代表由Ds和Dn对应的特征值构成的特征矢量。
沿x轴的方向矩阵可以表示为：
11????j2?dcos?`1sin?1/???j2?dcos?`1sin?1/?e?e? Ax????????j2?(N?1)dcos?sin?/???j2?(N?1)dcos?`1sin?1/?`11?e???e?沿y轴的方向矩阵可以表示为：
11????j2?dsin?`1sin?1/???j2?dsin?`1sin?1/?e?e? Ay????????j2?(N?1)dsin?sin?/???j2?(N?1)dsin?`1sin?1/?`11?e?e???2. 二维MUSIC算法
2.1 原理介绍
阵列协方差矩阵通过奇异分解，可以划分为噪声子空间和信号子空间，即
R?Us?sUsH?UN?NUNH
因为方向矩阵A中的各个列向量与噪声子空间正交，所以当方向矩阵中的角度为波达
方向时两者相乘的值会很小，根据这个性质，得到该阵列空间谱函数为
PMUSIC(?)?1??ay(?，?)?ax(?，?)??ENEN??ay(?，?)?ax(?，?)??HH
通过变化角度，找到的波峰位置就是估计的信源的二维角度。
2.2 算法流程
2.3 算法仿真
快拍数L=100，目标数K=3，8×8的方阵，假设源信号的仰角为10°，25，方位角为35°，45°，信噪比为20dB，宽带信号为基带频率为80Hz，带宽为40Hz的信号。并将该信号在时域上均分为5段。
二维宽带MUSIC第1段X: 36Y: 11Z: 11X: 46Y: 26Z: 0.87630.80.60.41005001008060402000.20
图 2 二维MUSIC第一段信号