发明名称 基于时间反演的低仰角波达方向估计方法
摘要 本发明公开了一种基于时间反演的低仰角波达方向估计方法,其思路为:确定天线阵列A,所述天线阵列A包含M个阵元,并且所述天线阵列A发送前向探测信号,所述前向探测信号经过目标散射后传播至天线阵列A,得到天线阵列A中第n个阵元接收经过目标散射传播回来的前向回波时域信号Y<sub>n</sub>(t),并计算天线阵列A中第n个阵元的时间反演接收信号<img file="DDA0001086532580000011.GIF" wi="146" he="70" />进而计算天线阵列A中第n个阵元的采样频率点ω<sub>q</sub>处时间反演频域接收信号的协方差矩阵<img file="DDA0001086532580000012.GIF" wi="195" he="71" />对<img file="DDA0001086532580000013.GIF" wi="169" he="71" />进行宽带聚焦处理,依次计算第n个阵元的宽带聚焦频率ω<sub>0</sub>处的聚焦协方差矩阵<img file="DDA0001086532580000014.GIF" wi="169" he="63" />和第n个阵元的宽带聚焦频率ω<sub>0</sub>处的聚焦平滑矩阵<img file="DDA0001086532580000015.GIF" wi="268" he="78" />并计算天线阵列A中第n个阵元的空间谱函数S(θ),进而计算天线阵列A中第n个阵元的空间谱函数S(θ)的H个极大值,每一个极大值对应的角度为波达方向。
申请公布号 CN106226754A 申请公布日期 2016.12.14
申请号 CN201610702654.8 申请日期 2016.08.22
申请人 西安电子科技大学 发明人 陈伯孝;曾小路;杨明磊
分类号 G01S7/41(2006.01)I;G01S7/539(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 西安睿通知识产权代理事务所(特殊普通合伙) 61218 代理人 惠文轩
主权项 一种基于时间反演的低仰角波达方向估计方法,其特征在于,包括以下步骤:步骤1,确定天线阵列A,所述天线阵列A包含M个阵元,并且所述天线阵列A发送前向探测信号,所述前向探测信号经过目标散射后传播至天线阵列A,得到天线阵列A中第n个阵元接收经过目标散射传播回来的前向回波时域信号Y<sub>n</sub>(t);t∈{1,2,…,S},S表示前向回波时域信号包含的时域快拍数,n∈{1,2,..,M},M表示天线阵列A包含的阵元个数;所述第n个阵元为天线阵列A中的任意一个阵元;步骤2,对天线阵列A中第n个阵元接收经过目标散射传播回来的前向回波时域信号Y<sub>n</sub>(t)进行时间反演操作,得到天线阵列A中第n个阵元发射的时间反演探测信号z<sub>n</sub>(t);t∈{1,2,…,S},S表示前向回波时域信号包含的时域快拍数;步骤3,根据天线阵列A中第n个阵元发射的时间反演探测信号z<sub>n</sub>(t),计算得到天线阵列A中第n个阵元的时间反演接收信号<img file="FDA0001086532550000011.GIF" wi="147" he="63" />步骤4,根据天线阵列A中第n个阵元的时间反演接收信号<img file="FDA0001086532550000012.GIF" wi="150" he="64" />计算得到第n个阵元中采样频率点ω<sub>q</sub>处时间反演频域接收信号的协方差矩阵<img file="FDA0001086532550000013.GIF" wi="195" he="71" />q∈{1,…,Q},ω<sub>q</sub>表示第q个采样点的采样频率,Q表示天线阵列A中第n个阵元的时间反演频域接收信号<img file="FDA0001086532550000014.GIF" wi="150" he="62" />进行带内采样的采样点个数;步骤5,对第n个阵元中采样频率点ω<sub>q</sub>处时间反演频域接收信号的协方差矩阵<img file="FDA0001086532550000015.GIF" wi="171" he="70" />进行宽带聚焦处理,依次计算得到第n个阵元的宽带聚焦频率ω<sub>0</sub>和第n个阵元中宽带聚焦频率ω<sub>0</sub>处的聚焦协方差矩阵<img file="FDA0001086532550000016.GIF" wi="195" he="63" />步骤6,对第n个阵元中宽带聚焦频率ω<sub>0</sub>处的聚焦协方差矩阵<img file="FDA0001086532550000017.GIF" wi="171" he="63" />进行解相干处理,计算得到第n个阵元中宽带聚焦频率ω<sub>0</sub>处的聚焦平滑矩阵<img file="FDA0001086532550000018.GIF" wi="267" he="78" />步骤7,根据第n个阵元中宽带聚焦频率ω<sub>0</sub>处的聚焦平滑矩阵<img file="FDA0001086532550000019.GIF" wi="268" he="79" />计算得到天线阵列A中第n个阵元的空间谱函数S(θ),进而计算天线阵列A中第n个阵元的空间谱函数S(θ)的H个极大值,每一个极大值对应的角度为波达方向;其中,H=P',P'表示对第n个阵元中宽带聚焦频率ω<sub>0</sub>处的聚焦协方差矩阵<img file="FDA0001086532550000021.GIF" wi="170" he="71" />划分的子阵个数。
地址 710071 陕西省西安市太白南路2号