[发明专利]一种冲击噪声背景下小快拍数信号波达方向估计方法有效
申请号: | 201410384625.2 | 申请日: | 2014-08-07 |
公开(公告)号: | CN104155629A | 公开(公告)日: | 2014-11-19 |
发明(设计)人: | 刘磊;刁鸣;高洪元 | 申请(专利权)人: | 哈尔滨工程大学 |
主分类号: | G01S3/00 | 分类号: | G01S3/00 |
代理公司: | 无 | 代理人: | 无 |
地址: | 150001 黑龙江省哈尔滨市南岗区*** | 国省代码: | 黑龙江;23 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明涉及一种冲击噪声环境下小快拍数信号波达方向估计方法。本发明包括:采集信号采样数据;对接收到的信号数据矢量进行去冲击预处理;对去冲击后的矩阵的协方差矩阵;将经过重采样构建后的协方差矩阵分解子空间得到其噪声子空间;对B次重采样的协方差矩阵分别进行求解,得到由B个入射角度估计值构成的矢量矩阵,利用置信空间对矩阵进行校正后求均值即可得到最终波达方向角度的估计值本发明经过对接收数据的预处理,能对冲击噪声进行有效地抑制;采用block-bootstrap方法对数据矩阵进行重采样构建,在小快拍数下能对信号进行有效地扩展;该方法稳健性高,适用于冲击噪声背景下的小快拍数DOA估计,具有良好的测向效果。 | ||
搜索关键词: | 一种 冲击 噪声 背景 下小快拍数 信号 方向 估计 方法 | ||
【主权项】:
一种冲击噪声背景下小快拍数信号波达方向估计方法,其特征在于:(1)采集信号采样数据:X(t)=AS(t)+N(t),其中,X(t)=[x1(t),x2(t),...,xL(t)]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)对接收到的信号数据矢量进行去冲击预处理:判断接收阵列数据矢量矩阵绝对值最大值Ximax与门限值之间的大小,其中门限值为为接收阵列数据矢量矩阵绝对值的均值,如果Ximax大于等于门限值,则用替换矩阵中最大值,如小于门限值则执行下一步骤;重复上述比较过程,直至矩阵中的值均小于门限值;Ximax≥aX‾i,Ximax=bX‾iXimax<aX‾i,Ximax=Ximax,]]>其中,b=0.6(2+α),α为冲击噪声特征系数,α∈(0,2],GSNR是广义信噪比;(3)对去冲击后的矩阵的协方差矩阵RX:RX=E[XXH],以协方差矩阵中每一列为单位进行可放回的随机均匀重采样B次后重构协方差矩阵RX*,B=200~1000,B∝1/N,N为快拍数;(4)将经过重采样构建后的协方差矩阵RX*分解子空间得到其噪声子空间UN,f(z)=eiHp(z),i=M+1,...,L]]>式中eiH是数据协方差矩阵RX*中小特征值对应的L‑M个特征矢量即代表着噪声子空间,其中p(z)=[1 z…zL‑1]Τ求解多项式的根:zi=exp(jwi)(1≤i≤M)得到入射角度:θi=arcsin(λ2πdarg{z^i}),i=1,...,M;]]>(5)对B次重采样的协方差矩阵分别进行求解,得到由B个入射角度估计值构成的矢量矩阵,利用置信空间对矩阵进行校正后求均值即可得到最终波达方向角度的估计值
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于哈尔滨工程大学,未经哈尔滨工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201410384625.2/,转载请声明来源钻瓜专利网。