1.圆形声矢量传感器阵列近场源多参数估计方法,其特征在于:所述声矢量传感器阵列由M个等间隔布置于半径为R的圆周上的阵元构成,以圆周的圆心为坐标原点且在坐标原点放置相同的声矢量传感器作为参考阵元,圆形声矢量传感器阵列的相邻阵元间的间隔为λmin/(8sin(π/M)),阵列的阵元是由声压传感器和x轴、y轴和z轴方向振速传感器构成的声矢量传感器,其中,λmin为入射信号的最小波长;
近场源多参数估计方法的步骤如下:阵列接收K个不同频率、互不相关窄带、随机平稳近场源信号,步骤一、利用圆形声矢量传感器阵列获取近场声波信号的接收数据;
该圆形声矢量传感器阵列接收信号的N次快拍数据构成直接采样数据Z1,接收信号延迟ΔT后的N次同步快拍数据构成延时采样数据Z2,由Z1和Z2这两组数据构成接收信号全数据其中 fs为奈奎斯特采样频率;
步骤二、利用接收信号全数据获取信号子空间和噪声子空间;
估计全数据相关矩阵 其中,A是全数据信号导向矢量矩阵,Rs为入射信号相关矩阵,σ2是高斯白噪声的功率,I是8M×8M的单位矩阵,根据子空间理论,对数据相关矩阵RZ进行特征分解获取信号子空间Us和噪声子空间UN,其中,EVD表示特征分解,λi是特征分解得到的第i个特征值,vi是特征值对应的第i个特征矢量,Us=[v1,…,vK]为K个大特征值对应的特征矢量构成的信号子空间,UN=[vK+1,…,v8M]为8M-K个小特征值对应的特征矢量构成的噪声子空间;
步骤三、估计信号导向矢量矩阵 和信号频率
将8M×K的信号子空间Us分成上下两块4M×K的矩阵U1和U2,利用时间旋转不变关系结构,由U1和U2通过矩阵运算得到ΨT=ΩT,其中 对矩阵Ψ进行特征分解,特征值构成矩阵 特征矢量构成矩阵 其中 是Ω的估计值,是T的估计值,从而得到信号导向矢量矩阵和信号频率的估计:其中,arg(·)表示取幅角, 表示矩阵 取的第k行第k列元素, Us=AT,U1=A1T,U2=A2T,A就是步骤二中的全数据阵列导向矢量矩阵,A1是直接采样数据阵列导向矢量矩阵,A2是延时采样数据阵列导向矢量矩阵,T是K×K的非奇异矩阵,是矩阵U1的伪逆矩阵,
步骤四、根据信号导向矢量的构成和阵列排布,获得第k个信号的x轴、y轴振速和声压对z轴振速的归一化矢量,并利用归一化矢量获得到达角和声源距离的粗略估计值;
根据声矢量传感器的结构和阵列的排布形式,将信号导向矢量估计值 的第k(1≤k≤K)列分成每个阵元对应的矩阵块, 表示 的第k列,表示第m个阵元的第k列, 的x轴、y轴方向振速和声压三个分量都与z轴方向振速分量相比得到归一化矢量 由M个阵元的归一化矢量平均得到第k个信号的归一化矢量 从而得到信号到达角粗略估计值 和声源距离粗略估计值第k个信号的归一化矢量为:
从而得到到达角和距离
的估计值:
其中,为第k个信号与坐标原点阵元之间的距离,λk为第k个信号的波长,ρ0是环境流体密度,c是声波传播速度, 向量Γk的第
1、2、4个元素,exp(·)为求指数运算,tan(·)和arctan(·)分别表示求正切和反正切运算;
步骤五、利用MUSIC算法在粗略值附近搜索信号到达角和距离的精确估计值;
利用圆形声矢量传感器阵列的结构形式,给出粗略估计值附近的小区域内的全数据阵列导向矢量 利用MUSIC谱峰搜索方法在粗略值附近搜索得到信号的精确估计值;
其中, 是信号到达阵元m和参考阵元的相位差构成的空域导向矢量, 为入射声源信号
到达阵元m和参考阵元的相位差,Un是步骤二得到的噪声子空间,单位能量信号的x轴、y轴和z轴方向的振速分量和声压标量为:θ,φ,r是搜索变量,
和 分别
是步骤四中的方位角、俯仰角和距离的粗略估计值,εθ、εφ和εr分别用来设置俯仰角、方位角和距离的搜索区间长度;
前述步骤中的k=1,...,K,m=1,...,M,n=1,...,M,j表示虚数单位。