1.一种冲击噪声环境下相干信号波达方向动态跟踪方法，其特征在于：
（1）获取信号采样数据：
T
X(t)＝AS(t)+N(t)，X(t)＝[x1(t),x2(t),...,xL(t)] 为快拍数为t时的L×1维接收数据矢量，A为阵列流型矩阵，S(t)为Μ×1维窄带信号矢量，N(t)为L×1维加性噪声矢量，噪声类型为复冲击噪声，A＝[a(θ1),a(θ2),...,a(θM)]， 为导向矢量，θi为第i个信源的入射角度，vi＝exp(-jπsin(θi)),i＝1,2,...,M；
（2）对接收到的信号数据矢量进行类归一化预处理：
z(t)为经过类归一化预处理的接收数据矢量，参数p的值由噪
声的特征指数α确定；
（3）将均匀线阵划分为若干个子阵：
将阵元数为L的均匀线阵划分为K个子阵，每个子阵的阵元数为N，N=L-K+1，则第k个子阵的前向空间平滑的接收数据矢量为：
zfk(t)＝[zk(t),zk+1(t),...,zk+N-1(t)]T
后向空间平滑的接收数据矢量为：
zbk(t)＝[zL-k+1(t),zL-k(t),...,zL-k-N+2(t)]H
其中，k=1,2,…,K；
（4）将经过前后向空间平滑的接收数据矢量带入PAST算法，得到信号子空间：
参数yfk(t)＝WH(t-1)zfk(t),ybk(t)＝WH(t-1)zbk(t)，W(t)为第t个快拍时的权矩阵，更新参数：
β为遗忘系数，0＜β＜1，更新参
数：
由先验信息得估计误差为efk(t)＝
zfk(t)-W(t-1)yfk(t)，
ebk( t) ＝ z bk(t )- W( t- 1) ybk( t) ， 权 矩 阵 更 新 为
权矩阵W(t)即为快拍数为t时对应
的信号特征矢量，由W(t)张成的空间为信号子空间，信号子空间US(t)＝orth[W(t)]；
（5）对信号子空间的数据使用MUSIC算法进行处理，经谱峰搜索获得目标信号波达方向：
信号子空间的MUSIC谱函数： 其中，I为单位
阵，对MUSIC谱峰值进行搜索，得到入射信号的波达方向角度；
（6）令t=t+1，将下一块拍数的数据矢量带入步骤（1）中，得到下一块拍数对应的信号波达方向角度，直到达到最大快拍数，实现冲击噪声环境下相干信号源的DOA动态跟踪。