[发明专利]双基前视SAR动目标成像方法与动目标速度估计方法有效
申请号: | 201410136027.3 | 申请日: | 2014-04-08 |
公开(公告)号: | CN103885062A | 公开(公告)日: | 2014-06-25 |
发明(设计)人: | 李中余;武俊杰;孙稚超;杨海光;黄钰林;杨建宇 | 申请(专利权)人: | 电子科技大学 |
主分类号: | G01S13/90 | 分类号: | G01S13/90;G01S13/58 |
代理公司: | 成都宏顺专利代理事务所(普通合伙) 51227 | 代理人: | 周永宏 |
地址: | 611731 四川省成*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种基于失配压缩的双基前视SAR动目标成像方法与动目标速度估计方法。本发明的成像方法首先推导出双基前视SAR动目标多普勒质心和多普勒调频率并构造出动目标方位信号;接着利用存在速度估计误差的参考函数与动目标方位信号进行失配压缩处理,然后基于驻定相位原理得到失配压缩处理的动目标散焦结果图像;最后利用该散焦图像的位置及散焦展宽宽度,与多普勒质心误差、多普勒调频率误差的关系,完成对多普勒质心误差、多普勒调频率误差的估计;进而可以得到动目标的准确多普勒质心、多普勒调频率,最后利用估计出的动目标准确多普勒质心、多普勒调频率完成对动目标的成像处理。 | ||
搜索关键词: | 双基前视 sar 目标 成像 方法 速度 估计 | ||
【主权项】:
1.一种基于失配压缩的BFL-SAR动目标成像方法,具体包括如下步骤:步骤一:建立BFL-SAR成像几何模型,完成参数初始化;设P为成像区域中的动目标,假设其距离向和方位向的运动速度分别为vr和va;BFL-SAR发射站与动目标P的斜视距离为RT,发射站速度为VT,发射站飞行方向与波束中心夹角为θ;接收站与动目标P的斜视距离为RR,发射站速度为VR,接收站飞行方向与波束中心夹角为零度;步骤二:获取BFL-SAR动目标的多普勒质心和多普勒调频率;BFL-SAR模式下,动目标P的多普勒质心fdc为:f dc = V R - v r λ + V T cos θ - v r cos θ - v a sin θ λ ]]> 其中,λ为发射信号载波波长,动目标P的多普勒调频率fdr为:f dr = v a 2 λ R R + ( V T sin θ + v r sin θ - v a cos θ ) 2 λ R T ]]> 假设动目标方位向速度估计误差和距离向速度估计误差分别为Δva和Δvr,则存在速度估计误差情况下的多普勒质心f′dc和多普勒调频率f′dr分别为:f dc ′ = V R - ( v r + Δ v r ) λ + V T cos θ - ( v r + Δ v r ) cos θ - ( v a + Δ v a ) sin θ λ = f dc + Δ f dc ]]> 其中,多普勒中心估计误差Δfdc为:Δ f dc = - Δ v r λ - Δ v r cos θ + Δ v a sin θ λ ; ]]>f dr ′ = ( v a + Δ v a ) 2 λ R R + ( V T sin θ - ( v r + Δ v r ) sin θ - ( v a + Δ v a ) cos θ ) 2 λ R T = f dr + Δ f dr ]]> 其中,多普勒调频率估计误差Δfdr为:Δ f dr = 2 v a Δ v a + Δ v a 2 λ R R + 2 ( V T sin θ + v r sin θ - v a cos θ ) ( Δ v r sin θ - Δ v a cos θ ) + ( Δ v r sin θ - Δ v a cos θ ) 2 λ R T ]]> 步骤三:构造出动目标方位信号和存在速度估计误差的参考函数;由步骤二,可得BFL-SAR动目标方位信号S(t)为:S ( t ) = rect [ t T ] exp { j 2 π ( f dc t + 1 2 f dr t 2 ) } ]]> 其中,rect[·]为方位时间窗,T为方位时宽,t为方位向时间。存在速度估计误差的参考函数Sref(t)为:S ref ( t ) = rect [ t T p ] exp { j 2 π ( f dc ′ t + 1 2 f dr ′ t 2 ) } = rect [ t T p ] exp { j 2 π ( f dc t + 1 2 f dr t 2 ) } exp { j 2 π ( Δ f dc t + 1 2 Δ f dr t 2 ) } ]]> 其中,Tp为参考函数时宽。步骤四:利用参考函数与动目标方位信号进行失配压缩处理,可得:S cmp ( t ) = ∫ - ∞ ∞ S ( ζ ) · S ref * ( ζ - t ) dζ = ∫ - ∞ ∞ rect [ ζ T ] · rect [ ζ - t T p ] exp { - j 2 π [ 1 2 Δ f dr ζ 2 + ( - f dr t + Δ f dc - Δ f dr t ) ζ + 1 2 f dr t 2 - f dc t - Δ f dc t + 1 2 Δ f dr t 2 ] } dζ = exp { - j 2 π ( c - b 2 4 a ) } ∫ - ∞ ∞ rect [ ζ T ] · rect [ ζ - t T p ] exp { - 2 πa ( ζ + b 2 a ) 2 } dζ = exp { - 2 π ( c - b 2 4 a ) } · L ( t ) ]]> 其中,ζ为时间变量,(*)表示共轭,
b=-fdrt+Δfdc-Δfdrt,c = 1 2 f dr t 2 - f dc t - Δ f dc t + 1 2 Δ f dr t 2 , ]]> 且L(t)的表达式为:L ( t ) = ∫ - ∞ ∞ rect [ ζ T ] exp { - j 2 πa ( ζ + b 2 a ) 2 } dζ ]]> 利用驻定相位原理,完成L(t)的积分;首先获取L(t)驻定相位点ζk为:
再将该驻定相位点表达式带入上式L(t)中,可得:L ( t ) = rect [ t - Δ f dc / f dr ′ Δ f dr / f dr ′ · T ] ∫ ζ k - Δ ζ k + Δ exp { - j 2 πa ( ζ - ζ k ) 2 } dζ ]]> 其中,[ζk-Δ,ζk+Δ]表示预先选取的积分区间,令η = 2 a ( ζ - ζ k ) , ]]> 上式可简化为:L ( t ) = rect [ t - Δ f dc / f dr ′ Δ f dr / f dr ′ · T ] 1 2 a ∫ - 2 a Δ 2 a Δ exp { - j π 2 η 2 } dη ]]> 对上式取模操作,可得失配压缩后动目标的散焦图像结果为:| L ( t ) | = rect [ f dr ′ t - Δ f dc Δ f dr T ] 1 2 a = rect [ t - Δ f dc / f dr ′ Δ f dr f dr ′ T ] 1 2 a ]]> 其中,|·|表示取模操作。步骤五:利用散焦图像的位置及散焦展宽宽度,与多普勒质心误差、多普勒调频率误差的关系,完成对多普勒质心误差、多普勒调频率误差的估计;由步骤四可得失配压缩后动目标的散焦图像的位置tmid与散焦展宽宽度Tdcmp,又因为:tmid=Δfdc/′drT dcmp = Δ f dr f dr ′ T ]]> 则可得多普勒质心误差、多普勒调频率误差的估计为:Δ f dc = f dr ′ · t mid Δ f dr = f dr ′ · T dcmp / T ]]> 步骤六:利用多普勒质心误差、多普勒调频率误差的估计得到动目标的准确多普勒质心、多普勒调频率,构造动目标方位准确参考信号,完成对动目标的成像处理;由步骤二及步骤五,可得动目标的准确多普勒质心fdc、多普勒调频率为fdr:f dc = f dr ′ - Δ f dc f dr = f dr ′ - Δ f dr ]]> 最后利用得到的动目标的准确多普勒质心fdc和多普勒调频率fdr,构造动目标方位准确参考信号S(t),完成对动目标的成像处理,成像处理结果Simge(t)为:S imge ( t ) = ∫ - ∞ ∞ S ( ζ ) · S * ( ζ - t ) dζ = sin c ( t - t 0 ) , ]]> 其中,sinc(·)为辛格函数,t0为动目标的成像结果位置点。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于电子科技大学,未经电子科技大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201410136027.3/,转载请声明来源钻瓜专利网。