本发明专利技术公开了一种基于杂波环境下的发射与接收联合优化自适应滤波方法,主要解决现有技术运算量较大,且未将接收机和发射机联合优化利用的问题。其实现过程是:1)发射一组相位编码信号;2)计算接收回波;3)对接收回波进行匹配滤波;4)将匹配滤波得到的散射中心幅度估值反馈至发射机,设置发射信号旁瓣权值;5)优化相位编码信号,使加权积分旁瓣电平最小;6)发射优化的相位编码信号并计算接收回波;7)对接收回波进行自适应滤波处理,得到散射中心的幅度估值。本发明专利技术对散射中心的幅度估计性能优于现有滤波方法,运算量低,并初步实现了认知雷达从接收机到发射机的闭环反馈,可用于雷达杂波抑制和目标检测。
【技术实现步骤摘要】
本专利技术属于雷达
,具体涉及一种发射与接收联合优化的自适应滤波方法,可用于雷达对杂波环境中感兴趣距离单元散射中心的幅度估计,实现从接收机到发射机的闭环反馈,提高估计精度并减小运算量。
技术介绍
随着现代电子科技技术的迅猛发展,雷达面临的工作环境日益复杂,现代雷达不仅要面对传统“四抗”问题,即抗电子干扰、抗反辐射导弹、反隐身、反低空突防,还要融合多种工作模式,处理多目标、多任务的综合能力。认知雷达作为一种可以根据目标和外部环境特性智能地选择发射信号、工作方式和资源分配的智能化雷达一经提出便受到了广泛的关注,它改变了传统自适应雷达单向的信息处理方式,实现了从接收到发射的闭环处理,被认为是未来雷达发展的重要方向。对环境特征的准确感知是认知雷达实现智能化的前提之一。雷达对外部环境的认知包括对地形特征,杂波和未知目标的统计特征等不同信息的认知。通常情况下,在雷达观测时间内,同一个距离单元散射中心的幅度是起伏变化的,因此雷达需要对多组接收回波进行滤波处理,以得到环境中感兴趣距离单元散射中心的幅度估值,实现对杂波和未知目标散射幅度统计信息的认知。估计散射中心幅度最简单直接的方法是匹配滤波,但这种方法仅适用于高斯白噪声下单个距离单元散射中心的幅度估计。对于多个距离单元散射中心幅度的估计问题,由于距离旁瓣的影响,匹配滤波会带来较大的估计误差,所以雷达接收机应对接收回波进行自适应滤波° T.Yardibi 等人在“Source Localization and Sensing:ANonparametricIterative Adaptive Approach Based on Weighted Least Squares, IEEETransactionsonAerospace and Electronic Systems, Vol.46,N0.1,January2010,,文章中基于加权最小均方代价函数提出了迭代的自适应滤波算法IAA,通过对邻近距离单元的去耦合实现自适应滤波,但该算法运算量较大,同时忽略了接收窗外邻近距离单元散射中心的影响。针对接收窗外散射中心的影响问题,S.D.Blunt等人在“Adaptive PulseCompression via MMSEEstimation, IEEE Transactions on Aerospace and ElectronicSystems, Vol.42, N0.2, April2006”中提出了一种递归最小均方误差算法RMMSE,通过递归的自适应信号处理方式实现各个距离单元散射中心的幅度估计。但同样存在运算量较大的问题,且容易受到设置参数的影响,易使滤波过程发散从而导致散射中心的幅度估计误差增大。
技术实现思路
本专利技术的目的在于针对上述已有技术的不足,在考虑接收窗外邻近距离单元散射中心影响的条件下,结合认知雷达的闭环反馈特点,提出一种,以减小滤波处理的运算量,提高散射中心幅度的估计精度。实现本专利技术目的的技术思路是:利用各个距离单元散射中心的幅度统计估值,联合优化发射信号和接收滤波器,以得到雷达观测环境中多个距离单元散射中心的幅度估值,其实现步骤包括如下:I)设置雷达发射信号模式由多个间隔i组成,i = 1,2,...,每个间隔i内包含M个相同的相位编码信号,编码方式采用P3码,用离散化向量s表示,编码长度为N ;2)计算雷达观测环境中第I个距离单元散射中心对应第m个脉冲信号的接收回波向量Y1,m,I = -N+1,…,0,…,L+N-2,L为雷达接收窗长度,m = 1,2,…,M, M为脉冲个数;3)对接收回波向量Y1, 进行匹配滤波,得到第I个距离单元散射中心的幅度估值权利要求1.一种,包括如下步骤: O设置雷达发射信号模式由多个间隔i组成,i = 1,2,…,每个间隔i内包含M个相同的相位编码信号,编码方式采用P3码,用离散化向量s表示,编码长度为N ; 2)计算雷达观测环境中第I个距离单元散射中心对应第m个脉冲信号的接收回波向量Y1,m, 1= -N+1,…,0,…,L+N-2,L为雷达接收窗长度,m = 1,2,…,M, M为脉冲个数; 3)对接收回波向量Y1,进行匹配滤波,得到第I个距离单元散射中心的幅度估值[h1|,其中,|.|表示模值; 4)将各个距离单元散射中心的幅度估值|h1|反馈至发射机,根据该值设置发射信号旁瓣的权值Y k:当 d = 1 时,Yk = 1, k = d,..., N-1,当 l<d<N 时2.根据权利要求1所述的方法,其中所述步骤(I)中相位编码信号向量S,表示如下: s = T, 其中,τ为转置操作,向量中第η个元素为3.根据权利要求1所述的方法,其中所述步骤(2)中计算雷达观测环境中第I个距离单元散射中心对应第m个脉冲信号的接收回波向量yi,m,按如下公式计算:4.根据权利要求1所述的方法,其中所述步骤(3)中对接收回波向量yi,m进行匹配滤波,得到第I个距离单元散射中心的幅度估值Pl,按如下公式进行:5.根据权利要求1所述的方法,其中所述步骤(5)中寻求优化的相位编码信号向量s',使加权积分旁瓣电平WISL最小,按如下公式进行: 其中,I.I2表示模值的平方,Yk为设定的权值,rk为信号旁瓣:6.根据权利要求1所述的方法,其中所述步骤(7b)中利用q-Ι次迭代得到的第I个距离单元散射中心对应第m个脉冲信号的幅度估值/ 广,求得第I个距离单元散射中心幅度均方估值£,按如下公式进行:7.根据权利要求1所述的方法,其中所述步骤(7c)中用第I个距离单元散射中心幅度均方估值£',求得第 个距离单元的滤波器向量Wl,按如下公式进行:W1 - RjiSr, 其中,—1表示矩阵求逆,R1为维数为NXN的相关矩阵,其表示式为:8.根据权利要求1所述的方法,其中所述步骤(7d)中用滤波器向量W1计算第q次迭代中第I个距离单元散射中心对应第m个脉冲信号的幅度估值,按如下公式进行:9.根据权利要求1所述的方法,其中所述步骤(7e)中计算第q-Ι次迭代和第q次迭代的幅度估值误差δ,按如下公式进行: 1],1V1 /=O L m=l_ 其中,I.I表示模值,m = 1,2,…,M,M为脉冲个数。全文摘要本专利技术公开了一种,主要解决现有技术运算量较大,且未将接收机和发射机联合优化利用的问题。其实现过程是1)发射一组相位编码信号;2)计算接收回波;3)对接收回波进行匹配滤波;4)将匹配滤波得到的散射中心幅度估值反馈至发射机,设置发射信号旁瓣权值;5)优化相位编码信号,使加权积分旁瓣电平最小;6)发射优化的相位编码信号并计算接收回波;7)对接收回波进行自适应滤波处理,得到散射中心的幅度估值。本专利技术对散射中心的幅度估计性能优于现有滤波方法,运算量低,并初步实现了认知雷达从接收机到发射机的闭环反馈,可用于雷达杂波抑制和目标检测。文档编号G01S7/28GK103116154SQ20131003057公开日2013年5月22日 申请日期2013年1月25日 优先权日2013年1月25日专利技术者刘峥, 吴旭姿, 刘韵佛 申请人:西安电子科技大学本文档来自技高网...
【技术保护点】
一种基于杂波环境下的发射与接收联合优化自适应滤波方法,包括如下步骤:?1)设置雷达发射信号模式由多个间隔i组成,i=1,2,…,每个间隔i内包含M个相同的相位编码信号,编码方式采用P3码,用离散化向量s表示,编码长度为N;?2)计算雷达观测环境中第l个距离单元散射中心对应第m个脉冲信号的接收回波向量yl,m,l=?N+1,…,0,…,L+N?2,L为雷达接收窗长度,m=1,2,…,M,M为脉冲个数;?3)对接收回波向量yl,m进行匹配滤波,得到第l个距离单元散射中心的幅度估值?其中,|·|表示模值;?4)将各个距离单元散射中心的幅度估值反馈至发射机,根据该值设置发射信号旁瓣的权值γk:?当d=1时,γk=1,k=d,…,N?1,?当1<d<N时,其中,d为旁瓣数;?5)寻求优化的相位编码信号向量s′,使加权积分旁瓣电平WISL最小;?6)将优化的相位编码信号向量s′作为下一个间隔i+1发射的脉冲信号,并计算第l个距离单元对应第m个脉冲信号的接收回波向量y′l,m;?7)对接收回波向量y′l,m进行自适应滤波:?(7a)对接收回波向量y′l,m进行匹配滤波,得到迭代初始值(7b)设迭代次数,q=1,2,…,在第q次迭代中,利用q?1次迭代得到的第l个距离单元散射中心对应第m个脉冲信号的幅度估值求得第l个距离单元散射中心幅度均方估值E[·]表示期望,|·|2表示模值的平方;?(7c)用第l个距离单元散射中心幅度均方估值求得第l个距离单元的?滤波器向量wl;?(7d)用滤波器向量wl计算第q次迭代中第l个距离单元散射中心对应第m个脉冲信号的幅度估值(7e)计算第q?1次迭代和第q次迭代的幅度估值误差δ,若幅度估值误差δ小于设定的阈值ρ=10?6,则停止迭代,得到第l个距离单元散射中心的幅度估值:?其中,表示最后一次迭代得到的第l个距离单元散射中心对应第m个脉冲信号的幅度估值;否则,重复步骤(7b)—(7e)继续迭代,直到幅度估值误差δ小于设定的阈值。?FDA00002781415200011.jpg,FDA00002781415200012.jpg,FDA00002781415200013.jpg,FDA00002781415200014.jpg,FDA00002781415200015.jpg,FDA00002781415200016.jpg,FDA00002781415200017.jpg,FDA00002781415200021.jpg,FDA00002781415200022.jpg,FDA00002781415200023.jpg...
【技术特征摘要】
【专利技术属性】
技术研发人员:刘峥,吴旭姿,刘韵佛,
申请(专利权)人:西安电子科技大学,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。