一种基于稀疏阵列的DOA跟踪方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN113777555A(43)申请公布日2021.12.10(21)申请号202110989967.7(22)申请日2021.08.26(71)申请人南京航空航天大学地址210016江苏省南京市秦淮区御道街29号(72)发明人张淇婷李建峰朱珂慧汤亚为(74)专利代理机构北京华沛德权律师事务所11302代理人马苗苗(51)Int.Cl.G01S3/14(2006.01)G06F17/14(2006.01)G06F17/16(2006.01)权利要求书1页说明书4页附图3页(54)发明名称一种基于稀疏阵列的DOA跟踪方法(57)摘要本发明公开了一种基于稀疏阵列的DOA跟踪方法，隶属于波达方向(DirectionofArrival,DOA)跟踪领域，在由一定的形式排列的无源阵元组成的波达方向跟踪系统中应用。具体包括：首先，使用稀疏阵列接收信号，得到每一个时刻的测量信息，对接收到的信号进行虚拟化处理；然后利用DFT算法获得信源起始时刻的DOA估计值；接着通过一阶泰勒展开构造总体最小二乘(TLS)来获取相邻时刻偏移量；最后在已估计出的角度基础上补偿相邻两个时刻的偏移量实现DOA跟踪。该方法可获得很高的估计精度，在性能和复杂度之间取得了良好折中，具有重要的应用价值。CN113777555ACN113777555A权利要求书1/1页1.一种基于偏移量补偿的DOA跟踪方法，其特征在于，包括以下步骤：步骤1)，稀疏阵列接收信号；步骤2)，在初始时刻利用DFT算法进行DOA估计，得到初始的角度值；步骤3)，通过一阶泰勒展开求解相邻时刻的偏移量；步骤4)，当前角度与所得偏移量进行叠加，得到下一时刻的角度值；步骤5)，通过估计出每一时刻的角度，实现了DOA跟踪。2.根据权利要求1所述的基于偏移量补偿的DOA跟踪方法，其特征在于，步骤2)中所述的对初始时刻的接收信号进行DOA估计，估计算法为DFT：求yini＝F·z1的K个峰值得到DOA估计值其中，K表示信源数，F是一个归一化矩阵，z1是连续虚拟阵列接收信号。3.根据权利要求1所述的基于偏移量补偿的DOA跟踪方法，其特征在于，步骤3)中所述的在处进行一阶泰勒展开求解t时刻与t+1时刻之间的偏移量：其中ξk为偏移量的值，t＝0,1,...,Ns，Ns为信源运动时间。4.根据权利要求1所述的基于偏移量补偿的DOA跟踪方法，其特征在于，步骤4)中所述的通过补偿偏移量获取下一时刻的DOA估计值：其中ξk为偏移量的值。2CN113777555A说明书1/4页一种基于稀疏阵列的DOA跟踪方法技术领域[0001]本发明属于波达方向(DirectionofArrival,DOA)跟踪领域，涉及一种基于稀疏阵列的DOA跟踪方法。背景技术[0002]稀疏阵列相比传统阵列能够获得更大的阵列孔径及更高的自由度，从而在谱估计精度和分辨率上更有优势。以二级嵌套阵列为例，它的构造方式可以表述为：由两个均匀线阵组成，第一级均匀线阵的阵元数为N1，阵元间隔为d，第二级均匀线阵的阵元数为N2，阵元间隔为(N1+1)d。该阵列结构突破了传统天线阵阵元间距半波长的限制，使得天线孔径得到扩展，能够获得角度估计性能的提升。[0003]现有的DOA估计方法假设信号在一定的观测时间内是不动的，然而，在动态环境中信号源是移动的，因此传统的DOA估计算法不能直接运用到DOA跟踪中。[0004]子空间方法是经典的跟踪算法，它通过更新信号子空间实现DOA跟踪。PAST及其变体PASTd算法都是典型的子空间跟踪算法。PAST通过最小化接收数据与当前子空间估计构建的投影近似之间的加权最小二乘(WLS)误差，递归地找到信号子空间，再借助DOA估计算法实现跟踪。虽然PAST算法在计算复杂度上有优势，但性能较为一般，尤其是在低信噪比时性能较差。发明内容[0005]本发明提出一种基于稀疏阵列的DOA跟踪方法，该方法首先利用离散傅里叶变换(DFT)算法估计出初始时刻的DOA值，然后通过补偿相邻时刻的偏移量估计出每一个时刻的角度，实现DOA跟踪。[0006]为解决上述问题，本发明采用以下技术方案：[0007]一种基于偏移量补偿的DOA跟踪方法，包括以下步骤：[0008](1)阵列接收信号；[0009](2)对初始时刻，即t＝0时刻的接收信号进行DOA估计，估计算法为DFT：求yini＝F·z1的K个峰值得到DOA估计值[0010][0011]其中，K表示信源数，F是一个归一化矩阵，z1是连续虚拟阵列接收信号。[0012](3)在处进行泰勒一阶展开求解第t时刻与第t+1时刻之间的偏移量：[0013][0014]其中ξk为偏移量的值，t＝0,1,...,Ns，Ns为信源运动时间。[0015](4)对角