CN104166134A - 一种实波束前视扫描雷达目标二维定位方法 - Google Patents
一种实波束前视扫描雷达目标二维定位方法 Download PDFInfo
- Publication number
- CN104166134A CN104166134A CN201410422691.4A CN201410422691A CN104166134A CN 104166134 A CN104166134 A CN 104166134A CN 201410422691 A CN201410422691 A CN 201410422691A CN 104166134 A CN104166134 A CN 104166134A
- Authority
- CN
- China
- Prior art keywords
- target
- distance
- orientation
- centerdot
- sigma
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种实波束前视扫描雷达目标二维定位方法,包括以下步骤:S1:成像***参数初始化,计算成像区域任意点目标与运动平台的距离,设置实波束扫描雷达点目标仿真参数;S2:进行距离向匹配滤波;S3:进行距离向运动补偿处理;S4:对扫描雷达方位向回波信号进行建模;S5:构造加权最小二乘目标函数;S6:进行目标方位定位。本发明利用发射信号和天线方向图参数信息,将目标幅度估计问题转化为目标函数关于加权向量wn的最优解的问题,通过求解目标函数的最优解求解出目标方位维的位置信息;有效的解决了运动平台前视作用区域中目标方位维定位精度低的问题,实现了运动平台前视区域目标距离维和方位维的二维精确定位。
Description
技术领域
本发明属于雷达技术领域,特别涉及一种实波束前视扫描雷达目标二维定位方法。
背景技术
雷达由于不受恶劣天气因素的影响,并能够全天时工作,因此在军事侦察、海洋及水文观测、陆/海追踪与救援等军用和民用领域发挥了不可或缺的作用。对实现运动平台前视区域目标的精确定位,对敌我目标侦察与探测、目标追踪与打击、海面目标搜救等功能的实现具有重要的意义。
运动平台前视区域目标的距离维高精度定位可以通过发射大带宽的线性调频信号和使用脉冲压缩技术处理实现。然而对于前视成像区域的方位维,由于平台与成像区域内目标相对运动产生的多普勒频率梯度几乎为零,使得现有的合成孔径等成像技术很难实现方位维目标的准确定位只能通过扫描成像的方式获得方位维目标的低精度定位结果。由于该方式目标位置的定位和和幅度的估计精度低,严重影响了运动平台的侦察、监视、定位和识别能力。
针对运动平台前视区域二维目标定位问题,特别是如何提高方位维目标的定位能力,一般采用两种方法。其一如文献:Blair W D,Brandt-Pearce M.Monopulse DOA estimation oftwo unresolved Rayleigh targets[J].Aerospace and Electronic Systems,IEEE Transactions on,2001,37(2):452-469.所采用单脉冲技术进行方位维处理。该技术基于单脉冲测角原理,主要适用于单个强点目标的定位,虽然对特定条件下的两点目标有效,但对于存在多散射中心的复杂目标环境下,目标定位会存在严重偏差,甚至会产生虚假目标等现象;其二如文献:Mahafza B R,Knight D L,Audeh N F.Forward-looking SAR imaging using a linear array withtransverse motion[C]//Southeastcon'93,Proceedings.,IEEE.IEEE,1993:4 p.的方法,该文章提出一种线性阵列合成孔径雷达前视成像方法,利用线性阵列与前视区域内目标之间的相对运动产生多普勒带宽,再利用匹配滤波技术实现目标方位定位。但该方法需要尽量长的线阵以增加孔径尺寸,同时,由于前视区域内目标的多普勒带宽很小,能够获得的目标定位精度依然有限。
发明内容
本发明的目的在于克服现有技术的不足,提供一种通过求解目标函数的最优解求解出目标方位维的位置信息,有效的解决了运动平台前视作用区域中目标方位维定位精度低的问题,实现了运动平台前视区域目标距离维和方位维的二维精确定位的实波束前视扫描雷达目标二维定位方法。
本发明的目的是通过以下技术方案来实现的:一种实波束前视扫描雷达目标二维定位方法,包括以下步骤:
S1:成像***参数初始化,计算成像区域任意点目标与运动平台的距离,设置实波束扫描雷达点目标仿真参数;
S2:进行距离向匹配滤波;
S3:进行距离向运动补偿处理;
S4:对扫描雷达方位向回波信号进行建模;
S5:构造加权最小二乘目标函数;
S6:进行目标方位定位。
进一步地,所述的步骤S1中计算成像区域任意点目标与运动平台的距离的具体方法为:运动平台零时刻位置记为(0,0,h),运动平台沿y轴运动,运动速度为V,目标相对平台的方位角记为q,雷达天线下视角记为雷达天线扫描速度记为ω;则t时刻时运动平台与场景中目标距离通过***参数表示为:
其中,R0为运动平台与目标之间的初始距离;
设置实波束扫描雷达点目标仿真参数的具体方法为:假设在扫描区域中同一距离R处的不同方位采样位置上都幅值,令产生这些运动幅值的目标的位置参数为θ=(θ1,θ2,...θN),幅度参数为σ=(σ1,σ2,...,σN),雷达发射信号为线性调频信号,扫描雷达作用区域的回波经过相干解调记为S(t,τ):
其中,τ为距离向时间变量,rect(·)和a(·)分别代表距离时间窗和方位时间窗,K是发射信号的时间调频斜率,c为光速,R(τ)代表运动平台与成像区域内各目标之间的距离变化;
扫描雷达成像区域的方位时间向量记为:
Ta=[-PRI·Na/2,-PRI·(Na/2-1),…,PRI·(Na/2-1)]
距离时间向量为:
Tr=[-1/fr·Nr/2,-1/fr·(Nr/2-1),…,1fr·(Nr/2-1)]
其中fr为距离向采样率,PRI为发射信号脉冲重复间隔,Na为方位向采样点数,Nr为距离向采样点数。。
进一步地,所述的步骤S2中距离向匹配滤波具体包括以下子步骤:
S21:对经过相干解调的回波S(t,τ)进行距离向脉冲压缩处理,获得距离维目标高分辨率,并通过距离向FFT得到距离向频域、方位向时域的回波信号S(fr,τ):
S22:构造距离向匹配滤波函数H(fr):
S23:将H(fr)与回波信号S(fr,τ)相乘得到距离压缩后的距离向频域、方位向时域的回波信号S1(fr,τ):
进一步地,所述的步骤S3中进行距离向运动补偿处理具体包括以下子步骤:
S31:将步骤S1中的进行泰勒展开;
S32:忽略展开后的距离关系表达式中的二次项,同时由于与较小,所以cosθ≈1,因此,使得R(x,y,t)≈R0-Vt;
S33:构造距离徙动因子H(fr,t):
S34:将H(fr,t)与S1(fr,τ)相乘消除雷达平台运动造成的距离徙动,并进行距离向IFFT变换得到高距离定位精度和低方位定位精度的二维时域信号S2(t,τ):
进一步地,所述的步骤S4中建模的具体方法为:对于各距离单元,方位扫描成像的回波模型及处理方式是相同的,因此任意选取任一距离单元的回波数据进行信号建模,方位向回波信号向量y表示为:
y=A(θ)x+n
其中,为方向矩阵,由各个方位采样点对应的方向向量组成,a(n)=[a1,…,aN]∈RL×1为天线方向图序列,N为一个波束宽度的采样点数,x=[x1,...,xN]表示方位向离散目标的幅度信息,M为方位向采样点数,y=[y1,...,yM]为方位向接收的回波信号,n为附加噪声向量。
进一步地,所述的步骤S5中构造加权最小二乘目标函数的具体方法为:对该距离单元的第n个目标,构造M×1维的加权向量wn,令并建立求该目标幅值的目标函数的最小二乘解:
其中,K为扫描雷达扫过目标场景的扫描次数,xn为第n个目标的幅值,展开目标函数得到:
其中, 为信号的协方差矩阵,令 则上式化简为:
进一步地,所述的步骤S6中进行目标方位定位具体包括以下子步骤:
S61:求等式 关于xn的导数并令其导数为零,得到关于xn幅度最优估计函数:
S62:将等式的目标幅度最优估计函数代入到式 中,并令 得到重构后的目标函数:
S63:求目标函数关J1(w)关于加权向量wn的最优解,其计算方法为:求目标函数关J1(w)关于wn的导数并令其为零,得到关于wn的最优解:
S64:将wn的计算结果代入到等式 中,得到关于xn的幅度为:
S65:利用步骤S63和S64的方法计算出该距离单元的所有目标幅度,并定位目标的角度,实现目标方位维的精确定位,再将算法应用到整个扫描雷达作用区域中,逐距离单元对整个面目标场景处理,实现成像区域内目标的二维精确定位。
进一步地,所述的方位向采样点数M的计算方法为:
其中,PRF为脉冲重复频率,ω为扫描速度,Φ为扫描范围。
本发明的有益效果是:本发明中采用实波束雷达以匀速扫描模式工作,雷达发射线性调频信号,利用发射信号和天线方向图参数信息,首先通过对回波信号的距离维匹配滤波处理获得目标距离维的精度定位,再进行运动补偿消除平台运动对回波信号的影响,实现目标距离维的高精度定位,随后建立实波束前视扫描雷达的信号回波模型,并基于加权最小二乘准则,构造了关于求解方位目标位置的目标函数,将目标幅度估计问题转化为目标函数关于加权向量wn的最优解的问题,通过求解目标函数的最优解求解出目标方位维的位置信息;有效的解决了运动平台前视作用区域中目标方位维定位精度低的问题,最终实现了运动平台前视区域目标距离维和方位维的二维精确定位,与传统二维目标定位方法相比,目标方位维定位精度更高,幅度的复原更准确。本发明可以应用于动目标追踪、精确制导等领域。
附图说明
图1为本发明的定位方法流程图;
图2为本发明实施例采用的实波束扫描雷达成像***结构图;
图3为本发明的实施实施例采用的仿真目标场景布置图;
图4为本发明的根据***参数生成的二维回波信号;
图5为本发明实施例的距离向脉压后的数据加入SNR=20dB高斯白噪声图形;
图6为本发明实施例的距离向脉压后数据加入SNR=20dB高斯白噪声沿方位向的剖面图;
图7为本发明的实施例图3中9个点目标进行二维目标定位处理的结果;
图8为本发明对应图6的处理结果沿方位向的剖面图。
具体实施方式
下面结合附图和具体实施例进一步说明本发明的技术方案,但本发明所保护的内容不局限于以下所述。
如图1所示,一种实波束前视扫描雷达目标二维定位方法,包括以下步骤:
S1:成像***参数初始化,计算成像区域任意点目标与运动平台的距离,设置实波束扫描雷达点目标仿真参数;
S2:进行距离向匹配滤波;
S3:进行距离向运动补偿处理;
S4:对扫描雷达方位向回波信号进行建模;
S5:构造加权最小二乘目标函数;
S6:进行目标方位定位。
进一步地,所述的步骤S1中计算成像区域任意点目标与运动平台的距离的具体方法为:运动平台零时刻位置记为(0,0,h),其中,0、0和h分别为接收站的x轴、y轴和z轴坐标;运动平台沿y轴运动,运动速度为V,目标相对平台的方位角记为θ,雷达天线下视角记为雷达天线扫描速度记为ω;则t时刻时运动平台与场景中目标距离通过***参数表示为:
其中,R0为运动平台与目标之间的初始距离;
设置实波束扫描雷达点目标仿真参数的具体方法为:假设在扫描区域中同一距离R处的不同方位采样位置上都幅值,令产生这些运动幅值的目标的位置参数为θ=(θ1,θ2,...θN),幅度参数为σ=(σ1,σ2,...,σN),雷达发射信号为线性调频信号,扫描雷达作用区域的回波经过相干解调记为S(t,τ):
其中,τ为距离向时间变量,rect(·)和a(·)分别代表距离时间窗和方位时间窗,K是发射信号的时间调频斜率,c为光速,R(τ)代表运动平台与成像区域内各目标之间的距离变化;
扫描雷达成像区域的方位时间向量记为:
Ta=[-PRI·Na/2,-PRI·(Na/2-1),…,PRI·(Na/2-1)]
距离时间向量为:
Tr=[-1/fr·Nr/2,-1/fr·(Nr/2-1),…,1/fr·(Nr/2-1)]
其中fr为距离向采样率,PRI为发射信号脉冲重复间隔,Na为方位向采样点数,Nr为距离向采样点数。
进一步地,所述的步骤S2中距离向匹配滤波具体包括以下子步骤:
S21:对经过相干解调的回波S(t,τ)进行距离向脉冲压缩处理,获得距离维目标高分辨率,并通过距离向FFT得到距离向频域、方位向时域的回波信号S(fr,τ):
S22:构造距离向匹配滤波函数H(fr):
S23:将H(fr)与回波信号S(fr,τ)相乘得到距离压缩后的距离向频域、方位向时域的回波信号S1(fr,τ):
进一步地,所述的步骤S3中进行距离向运动补偿处理具体包括以下子步骤:
S31:将步骤S1中的进行泰勒展开;
S32:忽略展开后的距离关系表达式中的二次项,同时由于与较小,所以cosθ≈1,因此,使得R(x,y,t)≈R0-Vt;
S33:构造距离徙动因子H(fr,t):
S34:将H(fr,t)与S1(fr,τ)相乘消除雷达平台运动造成的距离徙动,并进行距离向
IFFT变换得到高距离定位精度和低方位定位精度的二维时域信号S2(t,τ):
进一步地,所述的步骤S4中建模的具体方法为:对于各距离单元,方位扫描成像的回波模型及处理方式是相同的,因此任意选取任一距离单元的回波数据进行信号建模,方位向回波信号向量y表示为:
y=A(θ)x+n
其中,为方向矩阵,由各个方位采样点对应的方向向量组成,a(n)=[a1,…,aN]∈RL×1为天线方向图序列,N为一个波束宽度的采样点数,x=[x1,...,xN]表示方位向离散目标的幅度信息,M为方位向采样点数,y=[y1,...,yM]为方位向接收的回波信号,n为附加噪声向量。
进一步地,所述的步骤S5中构造加权最小二乘目标函数的具体方法为:对该距离单元的第n个目标,构造M×1维的加权向量wn,令并建立求该目标幅值的目标函数的最小二乘解:
其中,K为扫描雷达扫过目标场景的扫描次数,xn为第n个目标的幅值,展开目标函数得到:
其中, 为信号的协方差矩阵,令 则上式化简为:
进一步地,所述的步骤S6中进行目标方位定位具体包括以下子步骤:
S61:求等式 关于xn的导数并令其导数为零,得到关于xn幅度最优估计函数:
S62:将等式的目标幅度最优估计函数代入到式 中,并令 得到重构后的目标函数:
S63:通过上述操作,将目标幅度估计问题转化为求目标函数J1(w)关于加权向量wn的最优解问题,求目标函数关J1(w)关于加权向量wn的最优解的计算方法为:求目标函数关J1(w)关于wn的导数并令其为零,得到关于wn的最优解:
S64:将wn的计算结果代入到等式 中,得到关于xn的幅度为:
S65:利用步骤S63和S64的方法计算出该距离单元的所有目标幅度,并定位目标的角度,实现目标方位维的精确定位,再将算法应用到整个扫描雷达作用区域中,逐距离单元对整个面目标场景处理,实现成像区域内目标的二维精确定位。
进一步地,所述的方位向采样点数M的计算方法为:
其中,PRF为脉冲重复频率,ω为扫描速度,Φ为扫描范围。
下面结合具体实施例对本发明的技术方案进行进一步说明,主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2010上验证正确性。
步骤一:对成像区域任意点目标,计算目标与运动平台的距离,设置实波束扫描雷达点目标仿真参数,本发明采用如图2所示的实波束扫描雷达成像***结构图,对应的雷达成像***参数如表一所示。
表一
参数 | 符号 | 数值 |
载频 | fc | 10GHz |
带宽 | B | 20MHz |
发射信号时宽 | T | 50μs |
平台速度 | υ | 150m/s |
发射信号带宽 | B | 40MHz |
平台高度 | H | 5Km |
脉冲采样频率 | PRF | 1000Hz |
天线扫描速度 | ω | 30°/s |
天线波束宽度 | θ | 3° |
扫描范围 | Φ | -8°~8° |
本实施例采用的成像场景如图3所示,途中圆点为布置与地面上的3×3的点目标,沿y轴正方向,幅度依次为1、0.9、0.8,点目标沿方位向位置分别为-4°、2°和3.5°;沿x轴方向间隔为500m,雷达平台初始时刻位置坐标为(0,0,5km),xoy平面内目标散射函数记为f(x,y),t时刻xoy平面内的点(x,y)与雷达平台d的距离记为R(x,y,t)。
步骤二:根据步骤一设置的成像***参数和成像场景产生回波矩阵S(t,τ)并进行距离向FFT得到S(fr,τ),再根据发射信号调频斜率K和距离向参考时间τ,在频域构造距离向脉压参考函数,将S(fr,τ)与脉压参考函数进行最大自相关运算,完成距离向脉冲压缩,脉压后的距离向频域方位向时域的二维回波数据表示为S1(fr,τ),生成的回波信号如图4所示。
步骤三:根据前视区域内目标的斜距距离历史R(x,y,t)的泰勒级数展开结果R(x,y,t)≈R0-Vt,对数据S1(fr,τ)进行尺度变换和消除雷达平台运动造成的距离徙动,并进行距离向IFFT变换得到二维时域信号S2(t,τ)。为了模拟存在噪声的实际情况,在数据S2(t,τ)中加入SNR=20dB的高斯白噪声,相应的结果如图5所示,沿方位向剖面如图6所示。
步骤四:利用***设置的扫描速度、脉冲重复时间、天线波束宽度等***参数等生成方向向量a(θk)和方向矩阵A。
步骤五:根据生成的回波信号,利用公式 和 分别计算协方差矩阵R和向量g,并构造维数为M×1的加权向量wn。
步骤六:针对该距离单元第n个目标,首先,将步骤五计算出的协方差矩阵R和向量g代入到表达式中计算矩阵并将计算结果和构造的方向向量h(θn)代入到等式 中计算出求该目标幅度的最优加权向量wn,将wn的计算结果代入等式 中计算 计算出该目标的幅值和定位该目标的位置,并利用该方法逐点计算出方位向各目标的幅度和位置。
最后用步骤三到步骤五的方法逐距离单元对整个扫描雷达前视作用区域的距离向进行处理,得到整个运动平台扫描雷达前视成像区域内目标二维定位结果,成像结果如图7、图8所示,图示中,三个目标的方位定位位置分别为-4°、2.045°和3.459°,定位误差分别为0°、0.045°和0.041°。从图中可以看出,本发明可以实现运动平台扫描雷达前视区域进行二维目标定位处理,可以显著提高实波束扫描雷达目标二维定位精度,定位误差低,处理结果对于目标的方位维定位、相邻目标的分辨等具有良好的改善效果。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (8)
1.一种实波束前视扫描雷达目标二维定位方法,其特征在于:包括以下步骤:
S1:成像***参数初始化,计算成像区域任意点目标与运动平台的距离,设置实波束扫描雷达点目标仿真参数;
S2:进行距离向匹配滤波;
S3:进行距离向运动补偿处理;
S4:对扫描雷达方位向回波信号进行建模;
S5:构造加权最小二乘目标函数;
S6:进行目标方位定位。
2.根据权利要求1所述的定位方法,其特征在于:所述的步骤S1中计算成像区域任意点目标与运动平台的距离的具体方法为:运动平台零时刻位置记为(0,0,h),运动平台沿y轴运动,运动速度为V,目标相对平台的方位角记为q,雷达天线下视角记为雷达天线扫描速度记为ω;则t时刻时运动平台与场景中目标距离通过***参数表示为:
其中,R0为运动平台与目标之间的初始距离;
设置实波束扫描雷达点目标仿真参数的具体方法为:假设在扫描区域中同一距离R处的不同方位采样位置上都幅值,令产生这些运动幅值的目标的位置参数为θ=(θ1,θ2,...θN),幅度参数为σ=(σ1,σ2,...,σN),雷达发射信号为线性调频信号,扫描雷达作用区域的回波经过相干解调记为S(t,τ):
其中,τ为距离向时间变量,rect(·)和a(·)分别代表距离时间窗和方位时间窗,K是发射信号的时间调频斜率,c为光速,R(τ)代表运动平台与成像区域内各目标之间的距离变化;
扫描雷达成像区域的方位时间向量记为:
Ta=[-PRI·Na/2,-PRI·(Na/2-1),…,PRI·(Na/2-1)]
距离时间向量为:
Tr=[-1/fr·Nr/2,-1/fr·(Nr/2-1),…,1/fr·(Nr/2-1)]
其中fr为距离向采样率,PRI为发射信号脉冲重复间隔,Na为方位向采样点数,Nr为距离向采样点数。
3.根据权利要求2所述的定位方法,其特征在于:所述的步骤S2中距离向匹配滤波具体包括以下子步骤:
S21:对经过相干解调的回波S(t,τ)进行距离向脉冲压缩处理,获得距离维目标高分辨率,并通过距离向FFT得到距离向频域、方位向时域的回波信号S(fr,τ):
S22:构造距离向匹配滤波函数H(fr):
S23:将H(fr)与回波信号S(fr,τ)相乘得到距离压缩后的距离向频域、方位向时域的回波信号S1(fr,τ):
4.根据权利要求3所述的定位方法,其特征在于:所述的步骤S3中进行距离向运动补偿处理具体包括以下子步骤:
S31:将步骤S1中的进行泰勒展开;
S32:忽略展开后的距离关系表达式中的二次项,同时由于与较小,所以cosθ≈1,因此,使得R(x,y,t)≈R0-Vt;
S33:构造距离徙动因子H(fr,t):
S34:将H(fr,t)与S1(fr,τ)相乘消除雷达平台运动造成的距离徙动,并进行距离向IFFT变换得到高距离定位精度和低方位定位精度的二维时域信号S2(t,τ):
5.根据权利要求4所述的定位方法,其特征在于:所述的步骤S4中建模的具体方法为:对于各距离单元,方位扫描成像的回波模型及处理方式是相同的,因此任意选取任一距离单元的回波数据进行信号建模,方位向回波信号向量y表示为:
y=A(θ)x+n
其中,为方向矩阵,由各个方位采样点对应的方向向量组成,a(n)=[a1,…,aN]∈RL×1为天线方向图序列,N为一个波束宽度的采样点数,x=[x1,...,xN]表示方位向离散目标的幅度信息,M为方位向采样点数,y=[y1,...,yM]为方位向接收的回波信号,n为附加噪声向量。
6.根据权利要求5所述的定位方法,其特征在于:所述的步骤S5中构造加权最小二乘目标函数的具体方法为:对该距离单元的第n个目标,构造M×1维的加权向量wn,令并建立求该目标幅值的目标函数的最小二乘解:
其中,K为扫描雷达扫过目标场景的扫描次数,xn为第n个目标的幅值,展开目标函数得到:
其中, 为信号的协方差矩阵,令 则上式化简为:
7.根据权利要求6所述的定位方法,其特征在于:所述的步骤S6中进行目标方位定位具体包括以下子步骤:
S61:求等式 关于xn的导数并令其导数为零,得到关于xn幅度最优估计函数:
S62:将等式的目标幅度最优估计函数代入到式 中,并令 得到重构后的目标函数:
S63:求目标函数关J1(w)关于加权向量wn的最优解,其计算方法为:求目标函数关J1(w)关于wn的导数并令其为零,得到关于wn的最优解:
S64:将wn的计算结果代入到等式 中,得到关于xn的幅度为:
S65:利用步骤S63和S64的方法计算出该距离单元的所有目标幅度,并定位目标的角度,实现目标方位维的精确定位,再将算法应用到整个扫描雷达作用区域中,逐距离单元对整个面目标场景处理,实现成像区域内目标的二维精确定位。
8.根据权利要求5所述的定位方法,其特征在于:所述的方位向采样点数M的计算方法为:
其中,PRF为脉冲重复频率,ω为扫描速度,Φ为扫描范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410422691.4A CN104166134A (zh) | 2014-08-25 | 2014-08-25 | 一种实波束前视扫描雷达目标二维定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410422691.4A CN104166134A (zh) | 2014-08-25 | 2014-08-25 | 一种实波束前视扫描雷达目标二维定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104166134A true CN104166134A (zh) | 2014-11-26 |
Family
ID=51910034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410422691.4A Pending CN104166134A (zh) | 2014-08-25 | 2014-08-25 | 一种实波束前视扫描雷达目标二维定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104166134A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977582A (zh) * | 2015-06-10 | 2015-10-14 | 电子科技大学 | 一种实现扫描雷达方位超分辨成像的解卷积方法 |
CN106680817A (zh) * | 2016-12-26 | 2017-05-17 | 电子科技大学 | 一种实现前视雷达高分辨成像的方法 |
CN108226891A (zh) * | 2018-01-26 | 2018-06-29 | 中国电子科技集团公司第三十八研究所 | 一种扫描雷达回波计算方法 |
CN109765554A (zh) * | 2018-11-14 | 2019-05-17 | 北京遥感设备研究所 | 一种雷达前视成像***及方法 |
CN110402548A (zh) * | 2017-03-20 | 2019-11-01 | 华为技术有限公司 | 基于用户设备定位精度配置参考信号波束的设备 |
RU2741333C1 (ru) * | 2019-10-28 | 2021-01-25 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Санкт-Петербургский государственный университет телекоммуникаций им. проф. М.А. Бонч-Бруевича" | Способ определения местоположения работающей РЛС пассивным многолучевым пеленгатором |
CN117234217A (zh) * | 2023-11-13 | 2023-12-15 | 华中科技大学 | 基于三维时空域的水面无人艇轨迹跟踪制导方法及*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412305A (zh) * | 2013-07-15 | 2013-11-27 | 电子科技大学 | 扫描雷达超分辨成像方法 |
CN103869311A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 实波束扫描雷达超分辨成像方法 |
-
2014
- 2014-08-25 CN CN201410422691.4A patent/CN104166134A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103412305A (zh) * | 2013-07-15 | 2013-11-27 | 电子科技大学 | 扫描雷达超分辨成像方法 |
CN103869311A (zh) * | 2014-03-18 | 2014-06-18 | 电子科技大学 | 实波束扫描雷达超分辨成像方法 |
Non-Patent Citations (1)
Title |
---|
YIN ZHANG 等: "WEIGHTED LEAST SQUARES METHOD FOR FORWARD-LOOKING IMAGING OF SCANNING RADAR", 《GEOSCIENCE AND REMOTE SENSING SYMPOSIUM (IGARSS), 2014 IEEE INTERNATIONAL》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104977582B (zh) * | 2015-06-10 | 2018-09-04 | 电子科技大学 | 一种实现扫描雷达方位超分辨成像的解卷积方法 |
CN104977582A (zh) * | 2015-06-10 | 2015-10-14 | 电子科技大学 | 一种实现扫描雷达方位超分辨成像的解卷积方法 |
CN106680817B (zh) * | 2016-12-26 | 2020-09-15 | 电子科技大学 | 一种实现前视雷达高分辨成像的方法 |
CN106680817A (zh) * | 2016-12-26 | 2017-05-17 | 电子科技大学 | 一种实现前视雷达高分辨成像的方法 |
CN110402548A (zh) * | 2017-03-20 | 2019-11-01 | 华为技术有限公司 | 基于用户设备定位精度配置参考信号波束的设备 |
CN110402548B (zh) * | 2017-03-20 | 2021-05-04 | 华为技术有限公司 | 基于用户设备定位精度配置参考信号波束的设备 |
US11082104B2 (en) | 2017-03-20 | 2021-08-03 | Huawei Technologies Co., Ltd. | Apparatus for configuring reference signal beams based on accuracy of user equipment localization |
CN108226891A (zh) * | 2018-01-26 | 2018-06-29 | 中国电子科技集团公司第三十八研究所 | 一种扫描雷达回波计算方法 |
CN108226891B (zh) * | 2018-01-26 | 2021-09-03 | 中国电子科技集团公司第三十八研究所 | 一种扫描雷达回波计算方法 |
CN109765554A (zh) * | 2018-11-14 | 2019-05-17 | 北京遥感设备研究所 | 一种雷达前视成像***及方法 |
RU2741333C1 (ru) * | 2019-10-28 | 2021-01-25 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Санкт-Петербургский государственный университет телекоммуникаций им. проф. М.А. Бонч-Бруевича" | Способ определения местоположения работающей РЛС пассивным многолучевым пеленгатором |
CN117234217A (zh) * | 2023-11-13 | 2023-12-15 | 华中科技大学 | 基于三维时空域的水面无人艇轨迹跟踪制导方法及*** |
CN117234217B (zh) * | 2023-11-13 | 2024-02-02 | 华中科技大学 | 基于三维时空域的水面无人艇轨迹跟踪制导方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN104166134A (zh) | 一种实波束前视扫描雷达目标二维定位方法 | |
Chen et al. | Radon-fractional ambiguity function-based detection method of low-observable maneuvering target | |
CN102147469B (zh) | 一种双基地前视合成孔径雷达成像方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN103207387B (zh) | 一种机载相控阵pd雷达杂波的快速模拟方法 | |
CN104833972B (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN103487803A (zh) | 迭代压缩模式下机载扫描雷达成像方法 | |
CN102004250B (zh) | 基于频域展开的星机联合双基地合成孔径雷达成像方法 | |
CN104166129A (zh) | 一种实波束雷达迭代最小均方误差角超分辨方法 | |
CN107831480A (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN104749570A (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN102707269B (zh) | 一种机载雷达距离走动校正方法 | |
CN104698453B (zh) | 基于合成孔径天线阵列的雷达信号被动定位方法 | |
CN104035095A (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN104515971A (zh) | 宽带多目标机载单站无源定位方法 | |
CN100561250C (zh) | 星载双基地雷达的杂波基带模拟信号产生方法 | |
CN103487802A (zh) | 扫描雷达角超分辨成像方法 | |
CN105487074A (zh) | 一种双基地合成孔径雷达数值距离多普勒成像方法 | |
CN103308913A (zh) | 一种高速飞行器载双天线前视sar解模糊算法 | |
CN103529437A (zh) | 系留气球载相控阵雷达在多目标下分辨空地目标的方法 | |
CN105093224A (zh) | 一种大斜视合成孔径雷达成像处理方法 | |
CN112098999B (zh) | 一种高动态雷达导引头掠海目标电磁信号建模方法 | |
CN105137424A (zh) | 一种杂波背景下实波束扫描雷达角超分辨方法 | |
CN105137425A (zh) | 基于卷积反演原理的扫描雷达前视角超分辨方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141126 |