发明名称 一种实波束雷达迭代最小均方误差角超分辨方法
摘要 本发明公开了一种实波束雷达迭代最小均方误差角超分辨方法,克服了传统实波束雷达方位向角分辨率低的问题。通过建立方位回波信号的最小均方误差目标函数并求解,利用方位信号功率和协方差矩阵的计算关系构建迭代表达式,该方法通过迭代计算使对方位目标的估计逐步逼近真实目标位置,实现实波束雷达的方位角超分辨。
申请公布号 CN104166129A 申请公布日期 2014.11.26
申请号 CN201410415712.X 申请日期 2014.08.22
申请人 电子科技大学 发明人 张寅;黄钰林;邓敏;王园园;武俊杰;杨建宇
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 成都宏顺专利代理事务所(普通合伙) 51227 代理人 周永宏
主权项 一种实波速雷达迭代最小均方误差角超分辨方法,其特征在于,包括: S1:成像系统二维回波生成,具体包括以下步骤: S11:对回波信号进行相干解调: <img file="FDA0000557980040000011.GIF" wi="1633" he="145" />其中,S(t,τ)表示回波相干解调信号,rect(·)表示距离时间窗,α(·)表示方位时间窗,t表示距离向时间变量,τ表示方位向时间变量,λ表示发射信号波长,R<sub>0</sub>表示运动平台与目标之间的距离,K<sub>a</sub>表示发射信号的时间调频斜率,c表示光速;L表示发射信号个数,θ<sub>l</sub>表示第l个发射信号时天线的指向位置,σ<sub>l</sub>表示幅度参数; S12:扫描雷达成像区域的方位时间向量T<sub>a</sub>为: T<sub>a</sub>=[‑PRI·N<sub>a</sub>/2,‑PRI·(N<sub>a</sub>/2‑1),…,PRI·(N<sub>a</sub>/2‑1)]; 其中,PRI表示表示脉冲重复间隔,N<sub>a</sub>表示回波信号方位向采样点数; S13:扫描雷达成像区域的距离时间向量T<sub>r</sub>为: T<sub>r</sub>=[‑1/f<sub>r</sub>·N<sub>r</sub>/2,‑1/f<sub>r</sub>·(N<sub>r</sub>/2‑1),…,1/f<sub>r</sub>·(N<sub>r</sub>/2‑1)]; 其中,f<sub>r</sub>表示距离向采样率,N<sub>r</sub>表示回波信号距离向采样点数; S2:回波数据距离向脉冲压缩,具体包括以下分步骤: S21:对回波信号进行距离向脉冲压缩处理,通过距离向FFT得到距离向频域、方位向时域的回波信号S(f<sub>r</sub>,τ): <img file="FDA0000557980040000012.GIF" wi="1534" he="182" />S22:构造距离向匹配滤波函数H(f<sub>r</sub>): <img file="FDA0000557980040000013.GIF" wi="511" he="159" />S23:将H(f<sub>r</sub>)与回波信号相乘得到距离压缩后的距离向频域、方位向时域 的回波信号S<sub>1</sub>(f<sub>r</sub>,τ): <img file="FDA0000557980040000021.GIF" wi="1258" he="182" />S24:对S<sub>1</sub>(f<sub>r</sub>,τ)进行距离向IFFT变换得到二维时域信号S<sub>2</sub>(t,τ): <img file="FDA0000557980040000022.GIF" wi="1244" he="173" />S3:实波束雷达方位向信号建模,得到方位向回波信号向量y为: y=H(θ)f+n; 其中,H(θ)表示方向矩阵,f表示方位向离散目标的幅度信息,n表示附加噪声向量; 所述方位向回波信号向量y=[y<sub>1</sub>,…,y<sub>K</sub>]; 所述方位向离散目标的幅度信息f=[f<sub>1</sub>,…,f<sub>N</sub>]; S4:构造最小均方估计,具体包括以下分步骤: S41:针对方位回波信号,构造K×N维的加权矩阵w,求最优加权矩阵: minE{|f‑w<sup>H</sup>y|<sup>2</sup>}; 其中,E{·}表示均值计算,(·)<sup>H</sup>表示矩阵或者向量的共轭转置,|·|<sup>2</sup>表示向量各元素的平方和; S42:得到目标函数关于加权矩阵的最优解: w=(H(θ)PH(θ)<sup>H</sup>+R<sub>n</sub>)<sup>‑1</sup>H(θ)P 其中,H(θ)PH(θ)<sup>H</sup>表示信号的协方差矩阵,<img file="FDA0000557980040000024.GIF" wi="646" he="386" />表示信号功率谱矩阵,R<sub>n</sub>为噪声协方差矩阵;S42:求解出关于方位目标的估计值<img file="FDA0000557980040000031.GIF" wi="64" he="80" /><img file="FDA0000557980040000032.GIF" wi="978" he="157" />S5:迭代最小均方估计,通过构建迭代表达式以实现方位目标的精确估计; f<sub>σ</sub><sub>+1</sub>=(H(θ)P<sub>σ</sub>H(θ)<sup>H</sup>+R<sub>n</sub>)<sup>‑1</sup>H(θ)P<sub>k</sub>y; 其中σ表示迭代次数,<img file="FDA0000557980040000033.GIF" wi="736" he="396" />S6:判断迭代是否收敛; S7:实波束角超分辨成像。 
地址 611731 四川省成都市高新区(西区)西源大道2006号
您可能感兴趣的专利