CN113050053B - Method and system for acquiring coherent parameters of distributed coherent radar on moving platform - Google Patents

Method and system for acquiring coherent parameters of distributed coherent radar on moving platform Download PDF

Info

Publication number
CN113050053B
CN113050053B CN202110289756.2A CN202110289756A CN113050053B CN 113050053 B CN113050053 B CN 113050053B CN 202110289756 A CN202110289756 A CN 202110289756A CN 113050053 B CN113050053 B CN 113050053B
Authority
CN
China
Prior art keywords
coherent
model
moving platform
radar
transmission
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.)
Active
Application number
CN202110289756.2A
Other languages
Chinese (zh)
Other versions
CN113050053A (en
Inventor
杨琪
王元昊
曾旸
王宏强
邓彬
罗成高
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202110289756.2A priority Critical patent/CN113050053B/en
Publication of CN113050053A publication Critical patent/CN113050053A/en
Application granted granted Critical
Publication of CN113050053B publication Critical patent/CN113050053B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The application relates to a method and a system for acquiring phase parameters of a dynamic platform distributed phase-coherent radar. The method comprises the following steps: modeling a dynamic platform distributed coherent radar detection scene to obtain a dynamic platform emission signal model, a dynamic platform emission coherent model and a dynamic platform receiving coherent model, and then establishing a dynamic platform coherent calculation model according to the dynamic platform emission coherent model and the dynamic platform receiving coherent model; according to the coupling relation between the emission phase parameters, a state vector of the emission phase parameters is established, and a Singer model is adopted to model the state vector, so that a state equation of the emission phase parameters is obtained; according to the state equation, an observation equation is determined, and according to the observation equation, a Kalman filtering equation is determined, so that the coherent radar coherent parameters are calculated. The method can obtain the coherent radar coherent parameters of the mobile platform.

Description

动平台分布式相参雷达相参参数获取方法和系统Method and system for acquiring coherent parameters of distributed coherent radar on moving platform

技术领域Technical Field

本申请涉及信号处理技术领域,特别是涉及一种动平台分布式相参雷达相 参参数获取方法和系统。The present application relates to the field of signal processing technology, and in particular to a method and system for acquiring coherent parameters of a distributed coherent radar on a moving platform.

背景技术Background Art

低慢小无人机的广泛应用给社会各领域带来便捷的同时,也严重威胁着低空 领域的安全.近年来,由无人机造成的事故日益增多,研究针对低空领域无人机目 标的有效探测技术迫在眉睫。由于低慢小无人机目标反射区域少,单部雷达的 接收回波往往会存在信噪比不足的情况,检测困难。分布式相参雷达的核心思 想是使得多路发射信号在目标处同时叠加,使得回波信噪比大于单部雷达,从 而提升对小目标的探测能力。分布式相参雷达通过估计各个单元雷达之间的收 发时延差和相位差,这一参数称为相参参数(Coherentparameters,CPs),来使各 雷达之间信号相参。但是,由于目标无人机的可机动性,如果采取固定平台的 分布式相参雷达的方式,目标很容易飞入雷达的视野范围之外。而基于动平台 的分布式相参雷达有着固定平台式分布式相参雷达无可比拟的优势,以无人机 平台为例,其机动能力强,地形限制小,探测能力得以提升。但是,平台的机 动会带来诸多挑战,相参参数估计的滞后性就是动平台分布式相参雷达走向实 用化不可避免的一个问题,因为当目标和平台存在相对运动时,上一时刻认知 到的相参参数无法直接用于当前时刻的相参发射,需要研究相应的预测方法。The widespread use of low-speed, slow, and small drones has brought convenience to various fields of society, but it also seriously threatens the safety of low-altitude areas. In recent years, accidents caused by drones have increased day by day, and it is urgent to study effective detection technology for drone targets in low-altitude areas. Due to the small reflection area of low-speed, slow, and small drone targets, the received echo of a single radar often has insufficient signal-to-noise ratio, making detection difficult. The core idea of distributed coherent radar is to make multiple transmission signals superimposed at the target at the same time, so that the echo signal-to-noise ratio is greater than that of a single radar, thereby improving the detection capability of small targets. Distributed coherent radar makes the signals between radars coherent by estimating the transmission and reception delay difference and phase difference between each unit radar. This parameter is called coherent parameters (CPs). However, due to the maneuverability of the target drone, if a distributed coherent radar with a fixed platform is adopted, the target can easily fly out of the radar's field of view. The distributed coherent radar based on the moving platform has incomparable advantages over the fixed platform distributed coherent radar. Taking the UAV platform as an example, it has strong maneuverability, small terrain restrictions, and improved detection capabilities. However, the mobility of the platform will bring many challenges. The lag of coherent parameter estimation is an inevitable problem for the practical application of the distributed coherent radar of the moving platform. Because when the target and the platform are in relative motion, the coherent parameters recognized at the previous moment cannot be directly used for the coherent emission at the current moment, and the corresponding prediction method needs to be studied.

发明内容Summary of the invention

基于此,有必要针对上述技术问题,提供一种能够动平台分布式相参雷达 无法实现的动平台分布式相参雷达相参参数获取方法和系统。Based on this, it is necessary to provide a method and system for acquiring coherent parameters of a moving platform distributed coherent radar that cannot be achieved by a moving platform distributed coherent radar in response to the above-mentioned technical problems.

一种动平台分布式相参雷达相参参数获取方法,所述方法包括:A method for acquiring coherent parameters of a distributed coherent radar on a moving platform, the method comprising:

对动平台分布式相参雷达探测场景进行建模,得到动平台发射信号模型、 动平台发射相参模型、动平台接收相参模型;所述动平台分布式相参雷达探测 场景由多个独立的运动平台组成,所述运动平台之间相互独立;所述运动平台 之间通过无线链路传输数据;The moving platform distributed coherent radar detection scene is modeled to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model; the moving platform distributed coherent radar detection scene is composed of multiple independent moving platforms, and the moving platforms are independent of each other; the moving platforms transmit data through wireless links;

根据所述动平台发射相参模型和动平台接收相参模型,建立运动平台相参 计算模型;Establishing a motion platform coherence calculation model according to the motion platform transmitting coherence model and the motion platform receiving coherence model;

根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量,采用 Singer模型对所述状态向量进行建模,得到发射相参参数的状态方程;According to the coupling relationship between the emission coherent parameters, a state vector of the emission coherent parameters is established, and the state vector is modeled by using the Singer model to obtain a state equation of the emission coherent parameters;

根据所述状态方程,确定观测方程,以及根据所述观测方程,确定卡尔曼 滤波方程;Determine an observation equation according to the state equation, and determine a Kalman filter equation according to the observation equation;

根据所述卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参数序 列,根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参 雷达相参参数。According to the Kalman filter equation, the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter is determined, and according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained.

在其中一个实施例中,还包括:对动平台分布式相参雷达探测场景进行建 模,得到第m个运动平台发射信号的动平台发射信号模型为:In one embodiment, the method further includes: modeling a moving platform distributed coherent radar detection scene, and obtaining a moving platform transmission signal model of the m-th moving platform transmission signal as follows:

其中,Tp为发射脉宽,u为调频斜率,rect(t)为矩形函数,为载波,sm(t)=exp(j2π(m-1)Δft);Where, Tp is the transmit pulse width, u is the frequency modulation slope, rect(t) is the rectangular function, is the carrier, s m (t) = exp(j2π(m-1)Δft);

以m表示发射信号时运动平台的序号,l表示接收信号时运动平台的序号。 则第m个运动平台发射的信号为:m represents the number of the moving platform when transmitting the signal, and l represents the number of the moving platform when receiving the signal. Then the signal transmitted by the mth moving platform is:

式中κm表示雷达m较参考时钟的同步误差,为雷达m的初始相位。Where κ m represents the synchronization error of radar m compared to the reference clock, is the initial phase of radar m.

在其中一个实施例中,还包括:第m个运动平台发射的信号到达目标出的 信号表示为:In one embodiment, the signal emitted by the m-th moving platform and reaching the target is expressed as:

其中τm代表第m个运动平台发射的信号到达的时延,κm为该雷达较参考时 钟的同步误差,为该雷达较参考相位的同步误差;Where τ m represents the arrival delay of the signal emitted by the mth moving platform, κ m is the synchronization error of the radar relative to the reference clock, is the synchronization error of the radar relative to the reference phase;

则到达目标处总的信号为:The total signal reaching the target is:

设置雷达1为参考雷达,则调整后的各发射信号可以表示为:Set radar 1 as the reference radar, then the adjusted transmission signals can be expressed as:

其中为发射相参参数,得到动平台发射相参模型为:in and is the launch coherent parameter, and the launch coherent model of the moving platform is obtained as follows:

在其中一个实施例中,还包括:若第l个运动平台接收到目标反射回波表示 为:In one embodiment, the method further includes: if the lth moving platform receives the target reflected echo, it is expressed as:

其中,p(t)是目标处的回波信号;Where p(t) is the echo signal at the target;

则所有雷达接收到目标回波叠加为:Then the target echoes received by all radars are superimposed as:

设置雷达1为参考雷达,则调整后的各接收信号表示为:Set radar 1 as the reference radar, and the adjusted received signals are expressed as:

其中,为接收相参参数,得到动平台接收相参模型为:in, and To receive the coherent parameters, the moving platform receiving coherent model is obtained as:

在其中一个实施例中,还包括:根据所述动平台发射相参模型和动平台接 收相参模型,建立运动平台相参计算模型为:In one embodiment, the method further includes: establishing a moving platform coherent calculation model according to the moving platform transmitting coherent model and the moving platform receiving coherent model:

其中, in,

式中rl(n)是通过雷达直接量测得到的。Where r l (n) is directly measured by radar.

在其中一个实施例中,还包括:根据发射相参参数之间的耦合关系,建立 发射相参参数的状态向量为:In one embodiment, the method further includes: establishing a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters:

其中,R[n]表示状态向量,rl(n)表示发射相参参数;Where R[n] represents the state vector, r l (n) represents the transmission coherent parameter;

采用Singer模型对所述状态向量进行建模为:The state vector is modeled using the Singer model:

其中α是机动相关时间常数的倒数,即机动频率,是机动目标的加速度 方差;where α is the inverse of the maneuver-related time constant, i.e., the maneuver frequency, is the acceleration variance of the maneuvering target;

建立状态方程为:The state equation is established as:

R[n+1]=Φ(T,α)R[n]+u[n]R[n+1]=Φ(T,α)R[n]+u[n]

驱动噪声协方差为:The driving noise covariance is:

在其中一个实施例中,还包括:根据所述状态方程,确定观测方程为:In one embodiment, the method further includes: determining an observation equation according to the state equation as:

z[n]=HX[n]+v[n]z[n]=HX[n]+v[n]

其中,H=[1 0 0],v[n]=σ2,σ2为雷达量测噪声;Where, H = [1 0 0], v[n] = σ 2 , σ 2 is the radar measurement noise;

根据所述观测方程,确定卡尔曼滤波方程为:According to the observation equation, the Kalman filter equation is determined as:

P[n|n-1]=ΦP[n|n]ΦT+Q[n]P[n|n-1]=ΦP[n|n]Φ T +Q[n]

K[n]=P[n|n-1]HT(HP[n|n-1]HT+R)-1 K[n]=P[n|n-1]H T (HP[n|n-1]H T +R) -1

P[n|n]=(I-K[n]H)P[n|n-1]P[n|n]=(I-K[n]H)P[n|n-1]

在其中一个实施例中,还包括:据所述卡尔曼滤波方程,确定发射相参参 数对应的预测发射相参参数序列为:In one of the embodiments, the method further includes: determining, according to the Kalman filter equation, a predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter as:

根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参 雷达相参参数为:According to the predicted emission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained as follows:

其中, in,

一种动平台分布式相参雷达相参参数获取系统,所述系统包括:A coherent parameter acquisition system for a distributed coherent radar on a moving platform, the system comprising:

场景建模模块,用于对动平台分布式相参雷达探测场景进行建模,得到动 平台发射信号模型、动平台发射相参模型、动平台接收相参模型;所述动平台 分布式相参雷达探测场景由多个独立的运动平台组成,所述运动平台之间相互 独立;所述运动平台之间通过无线链路传输数据;A scene modeling module is used to model the moving platform distributed coherent radar detection scene to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model; the moving platform distributed coherent radar detection scene is composed of a plurality of independent moving platforms, and the moving platforms are independent of each other; the moving platforms transmit data through wireless links;

卡尔曼滤波模块,用于根据所述动平台发射相参模型和动平台接收相参模 型,建立运动平台相参计算模型;根据发射相参参数之间的耦合关系,建立发 射相参参数的状态向量,采用Singer模型对所述状态向量进行建模,得到发射 相参参数的状态方程;根据所述状态方程,确定观测方程,以及根据所述观测 方程,确定卡尔曼滤波方程;A Kalman filter module is used to establish a motion platform coherent calculation model according to the motion platform transmission coherent model and the motion platform receiving coherent model; establish a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters, use the Singer model to model the state vector, and obtain the state equation of the transmission coherent parameters; determine the observation equation according to the state equation, and determine the Kalman filter equation according to the observation equation;

相参参数计算模块,用于根据所述卡尔曼滤波方程,确定发射相参参数对 应的预测发射相参参数序列,根据所述预测发射相参参数序列和所述运动平台 相参计算模型,得到相参雷达相参参数。The coherent parameter calculation module is used to determine the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter according to the Kalman filter equation, and obtain the coherent radar coherent parameters according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model.

上述动平台分布式相参雷达相参参数获取方法和系统,利用发射相参参数 和不同平台与目标间距的耦合关系,通过量测目标与平台间距并获得该距离变 化序列,设计卡尔曼滤波器对距离变化序列进行一步预测,根据预测距离值估 计出下一时刻的发射相参参数,为动平台分布式相参雷达相参参数获取及系统 设计提供了研究支撑。The above-mentioned moving platform distributed coherent radar coherent parameter acquisition method and system utilize the coupling relationship between the transmission coherent parameters and the distances between different platforms and targets, measure the distance between the target and the platform and obtain the distance change sequence, design a Kalman filter to make a one-step prediction of the distance change sequence, and estimate the transmission coherent parameters at the next moment according to the predicted distance value, which provides research support for the coherent parameter acquisition and system design of the moving platform distributed coherent radar.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为一个实施例中动平台分布式相参雷达相参参数获取方法的流程示意 图;FIG1 is a schematic diagram of a process of acquiring coherent parameters of a distributed coherent radar on a moving platform in one embodiment;

图2为一个实施例中运动平台探测目标示意图;FIG2 is a schematic diagram of a moving platform detecting a target in one embodiment;

图3为一个实施例中相参参数估计流程示意图;FIG3 is a schematic diagram of a coherent parameter estimation process in one embodiment;

图4为一个实施例中动平台分布式相参雷达仿真场景图;FIG4 is a diagram of a simulation scene of a distributed coherent radar on a moving platform in one embodiment;

图5为一个实施例中各雷达量测距离误差图(基于singer模型);FIG5 is a diagram of the distance error measured by each radar in one embodiment (based on the Singer model);

图6为一个实施例中各雷达量测距离误差图(无singer模型)FIG. 6 is a diagram showing the distance error of each radar measurement in one embodiment (without the Singer model)

图7为一个实施例中雷达2相参参数误差变化图;FIG. 7 is a diagram showing a variation of a coherent parameter error of a radar 2 in one embodiment;

图8为一个实施例中雷达3相参参数误差变化图;FIG8 is a diagram showing a variation of a coherent parameter error of radar 3 in one embodiment;

图9为一个实施例中动平台分布式相参雷达相参参数获取系统的结构框图。FIG9 is a structural block diagram of a moving platform distributed coherent radar coherent parameter acquisition system in one embodiment.

具体实施方式DETAILED DESCRIPTION

为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实 施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅 用以解释本申请,并不用于限定本申请。In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not intended to limit the present application.

在一个实施例中,如图1所示,提供了一种动平台分布式相参雷达相参参 数获取方法,包括以下步骤:In one embodiment, as shown in FIG1 , a method for acquiring coherent parameters of a distributed coherent radar on a moving platform is provided, comprising the following steps:

步骤102,对动平台分布式相参雷达探测场景进行建模,得到动平台发射信 号模型、动平台发射相参模型、动平台接收相参模型。Step 102, modeling the moving platform distributed coherent radar detection scene to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model.

动平台分布式相参雷达可由若干运动平台组成,各个运动平台之间彼此相 互独立,可以通过无线链路传输数据,因此可以有多种工作方式供选择。在 MIMO模式下,各个运动平台发射正交信号,同时各个运动平台进行信号的接 收。因为发射的是正交信号,所以在接收运动平台端可以通过滤波等方式将不 同运动平台发射的回波进行分离,然后选取某一通道的信号作为参考信号,其 余通道的信号以参考通道信号为准,在时延和相位上与参考信号对齐,此时完 成接收相参。通过MIMO模式得到精确的相参参数后,系统可切换为全相参工 作模式,全相参模式下系统控制不同运动平台发射信号的时间和初相,使得所有运动平台发射的信号在同一时刻以同一相位到达目标处,从而在目标处形成 一次叠加,达到发射相参,在此基础之上再对叠加后的回波进行处理。以下分 别讨论动平台分布式相参系统的两种工作模式。系统工作在MIMO模式时,需 要将不同运动平台发射信号的回波进行分离,所以需要采用正交信号,此处采 用正交频分信号发射波形。The mobile platform distributed coherent radar can be composed of several mobile platforms. Each mobile platform is independent of each other and can transmit data through a wireless link. Therefore, there are multiple working modes to choose from. In the MIMO mode, each mobile platform transmits an orthogonal signal and receives the signal at the same time. Because the transmitted signal is an orthogonal signal, the echoes transmitted by different mobile platforms can be separated by filtering and other methods at the receiving mobile platform end, and then the signal of a certain channel is selected as the reference signal. The signals of the remaining channels are based on the reference channel signal and aligned with the reference signal in terms of delay and phase. At this time, the receiving coherence is completed. After obtaining accurate coherence parameters through the MIMO mode, the system can switch to the full coherence working mode. In the full coherence mode, the system controls the time and initial phase of the signals transmitted by different mobile platforms, so that the signals transmitted by all mobile platforms arrive at the target at the same time and with the same phase, thereby forming a superposition at the target to achieve transmission coherence. On this basis, the superimposed echoes are processed. The following discusses the two working modes of the mobile platform distributed coherence system. When the system works in MIMO mode, the echoes of signals transmitted by different moving platforms need to be separated, so an orthogonal signal is required. Here, an orthogonal frequency division signal is used to transmit the waveform.

步骤104,根据动平台发射相参模型和动平台接收相参模型,建立运动平台 相参计算模型。Step 104: Establish a motion platform coherent calculation model based on the motion platform transmitting coherent model and the motion platform receiving coherent model.

因为运动平台和目标相对运动产生的距离走动不可忽略,也会导致相参参 数随之发生变化,这一变化主要由路径差的变化所导致。发射相参的目的是要 在目标处发生正向干涉,但是由于平台和目标的相对位置的实时变化,且相参 参数认知又始终发生在相参合成之前,这就意味着仅靠一次观测下认知到的发 射相参参数无法校正发射去相参,即所认知知识具有滞后性(时间和相位同步 误差考虑不变),简单来说,能直接量测到的是。探测场景如图2所示,在平台 (目标)运动场景下,有必要对发射相参参数进行合理的预测,以补偿认知知 识滞后所带来的去相参。简单来说,n时刻的发射相参参数需要根据n时刻的状态来确定,但量测得到的最新状态还是处于n-1时刻。Because the distance movement caused by the relative motion of the moving platform and the target cannot be ignored, the coherence parameters will also change accordingly. This change is mainly caused by the change of the path difference. The purpose of transmission coherence is to cause forward interference at the target. However, due to the real-time change of the relative position of the platform and the target, and the coherence parameter cognition always occurs before the coherence synthesis, this means that the transmission coherence parameters recognized under only one observation cannot correct the transmission decoherence, that is, the recognized knowledge has a lag (time and phase synchronization errors are considered unchanged). In simple terms, what can be directly measured is. The detection scene is shown in Figure 2. In the platform (target) motion scene, it is necessary to make a reasonable prediction of the transmission coherence parameters to compensate for the decoherence caused by the lag of cognitive knowledge. In simple terms, the transmission coherence parameters at time n need to be determined based on the state at time n, but the latest state measured is still at time n-1.

步骤106,根据发射相参参数之间的耦合关系,建立发射相参参数的状态向 量,采用Singer模型对所述状态向量进行建模,得到发射相参参数的状态方程。Step 106: Establish a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters, use the Singer model to model the state vector, and obtain the state equation of the transmission coherent parameters.

步骤108,根据状态方程,确定观测方程,以及根据观测方程,确定卡尔曼 滤波方程。Step 108, determining the observation equation based on the state equation, and determining the Kalman filter equation based on the observation equation.

步骤110,根据卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参 数序列,根据预测发射相参参数序列和运动平台相参计算模型,得到相参雷达 相参参数。Step 110, according to the Kalman filter equation, determine the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter, and obtain the coherent radar coherent parameters according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model.

上述动平台分布式相参雷达相参参数获取方法中,利用发射相参参数和不 同平台与目标间距的耦合关系,通过量测目标与平台间距并获得该距离变化序 列,设计卡尔曼滤波器对距离变化序列进行一步预测,根据预测距离值估计出 下一时刻的发射相参参数,为动平台分布式相参雷达相参参数获取及系统设计 提供了研究支撑。In the above moving platform distributed coherent radar coherent parameter acquisition method, the coupling relationship between the transmission coherent parameters and the distances between different platforms and targets is utilized. By measuring the distance between the target and the platform and obtaining the distance change sequence, a Kalman filter is designed to make a one-step prediction of the distance change sequence, and the transmission coherent parameters at the next moment are estimated according to the predicted distance value, which provides research support for the coherent parameter acquisition and system design of the moving platform distributed coherent radar.

在其中一个实施例中,对动平台分布式相参雷达探测场景进行建模,得到 第m个运动平台发射信号的动平台发射信号模型为:In one embodiment, a moving platform distributed coherent radar detection scenario is modeled, and the moving platform transmission signal model of the m-th moving platform transmission signal is obtained as follows:

其中,Tp为发射脉宽,u为调频斜率,rect(t)为矩形函数,为载波,sm(t)=exp(j2π(m-1)Δft);Where, Tp is the transmit pulse width, u is the frequency modulation slope, rect(t) is the rectangular function, is the carrier, s m (t) = exp(j2π(m-1)Δft);

以m表示发射信号时运动平台的序号,l表示接收信号时运动平台的序号。 则第m个运动平台发射的信号为:m represents the number of the moving platform when transmitting the signal, and l represents the number of the moving platform when receiving the signal. Then the signal transmitted by the mth moving platform is:

式中κm表示雷达m较参考时钟的同步误差,为雷达m的初始相位。Where κ m represents the synchronization error of radar m compared to the reference clock, is the initial phase of radar m.

在其中一个实施例中,第m个运动平台发射的信号到达目标出的信号表示 为:In one embodiment, the signal emitted by the m-th moving platform and reaching the target is expressed as:

其中τm代表第m个运动平台发射的信号到达的时延,κm为该雷达较参考时 钟的同步误差,为该雷达较参考相位的同步误差;Where τ m represents the arrival delay of the signal emitted by the mth moving platform, κ m is the synchronization error of the radar relative to the reference clock, is the synchronization error of the radar relative to the reference phase;

则到达目标处总的信号为:The total signal reaching the target is:

设置雷达1为参考雷达,则调整后的各发射信号可以表示为:Set radar 1 as the reference radar, then the adjusted transmission signals can be expressed as:

其中为发射相参参数,得到动平台发射相参模型为:in and is the launch coherent parameter, and the launch coherent model of the moving platform is obtained as follows:

在其中一个实施例中,若第l个运动平台接收到目标反射回波表示为:In one embodiment, if the lth moving platform receives the target reflected echo, it is expressed as:

其中,p(t)是目标处的回波信号;Where p(t) is the echo signal at the target;

则所有雷达接收到目标回波叠加为:Then the target echoes received by all radars are superimposed as:

设置雷达1为参考雷达,则调整后的各接收信号表示为:Set radar 1 as the reference radar, and the adjusted received signals are expressed as:

其中,为接收相参参数,得到动平台接收相参模型为:in, and To receive the coherent parameters, the moving platform receiving coherent model is obtained as:

显然,在目标处,各个平台接收到的信号并不能同向叠加,为了实现电磁 波能量的相干叠加,也就是接收相参,需要调节各个雷达的接收延时和接收相 位。Obviously, at the target, the signals received by each platform cannot be superimposed in the same direction. In order to achieve the coherent superposition of electromagnetic wave energy, that is, receiving coherence, it is necessary to adjust the receiving delay and receiving phase of each radar.

在其中一个实施例中,假设目标和平台在单脉冲内可近似认为静止不动, 由式动平台发射相参模型和动平台接收相参模型可知,相参参数并不能直接量 测得到,为了估计相参参数,首先需要各个雷达发射正交信号,接收端通过匹 配滤波分离出自发自收回波和它发自收回波,记分离出的回波为In one embodiment, it is assumed that the target and the platform can be approximately considered to be stationary within a single pulse. According to the moving platform transmission coherence model and the moving platform reception coherence model, the coherence parameters cannot be directly measured. In order to estimate the coherence parameters, each radar needs to transmit orthogonal signals first, and the receiving end separates the self-transmitted and self-received echoes and the self-transmitted and self-received echoes through matched filtering. The separated echoes are recorded as

其中τlm=τlmml代表第m个运动平台发射的信号经过目标反射后到达 第l个运动平台的时延,αlm为该路径上的散射体响应,为两雷达的相 位同步误差。Where τ lm = τ l + τ m + κ m - κ l represents the time delay of the signal emitted by the mth moving platform reaching the lth moving platform after being reflected by the target, α lm is the response of the scatterer on the path, is the phase synchronization error between the two radars.

相参参数估计流程示意图如图3所示,通过相参参数的定义可得:The schematic diagram of the coherent parameter estimation process is shown in Figure 3. Through the definition of the coherent parameters, we can get:

当系统发射相同波形时,接收端依旧可以获得接收相参参数,但是无法再 次获得发射相参参数,因此必须建立发射相参参数和接收相参参数的转换关系, 比较动平台发射相参模型和动平台接收相参模型可得:When the system transmits the same waveform, the receiving end can still obtain the receiving coherent parameters, but cannot obtain the transmitting coherent parameters again. Therefore, the conversion relationship between the transmitting coherent parameters and the receiving coherent parameters must be established. By comparing the moving platform transmitting coherent model and the moving platform receiving coherent model, we can get:

在其中一个实施例中,根据动平台发射相参模型和动平台接收相参模型, 建立运动平台相参计算模型为:In one embodiment, according to the moving platform transmitting coherent model and the moving platform receiving coherent model, a moving platform coherent calculation model is established as follows:

其中, in,

式中rl(n)是通过雷达直接量测得到的。Where r l (n) is directly measured by radar.

在其中一个实施例中,已知发射相参参数与rl(n),l=1,2,…,N存在线性耦合关系,即如果能根据rl(n)预测出rl(n+1)便可以预测出n+1时刻的发射相参参数。In one embodiment, it is known that the transmission coherent parameters have a linear coupling relationship with r l (n), l=1, 2, ..., N, that is, if r l (n+1) can be predicted based on r l (n), the transmission coherent parameters at time n+1 can be predicted.

根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量为:According to the coupling relationship between the emission coherent parameters, the state vector of the emission coherent parameters is established as:

其中,R[n]表示状态向量,rl(n)表示发射相参参数;Where R[n] represents the state vector, r l (n) represents the transmission coherent parameter;

采用Singer模型对所述状态向量进行建模为:The state vector is modeled using the Singer model:

其中α是机动相关时间常数的倒数,即机动频率,是机动目标的加速度 方差;where α is the inverse of the maneuver-related time constant, i.e., the maneuver frequency, is the acceleration variance of the maneuvering target;

建立状态方程为:The state equation is established as:

R[n+1]=Φ(T,α)R[n]+u[n]R[n+1]=Φ(T,α)R[n]+u[n]

驱动噪声协方差为:The driving noise covariance is:

在其中一个实施例中,根据所述状态方程,确定观测方程为:In one embodiment, according to the state equation, the observation equation is determined as:

z[n]=HX[n]+v[n]z[n]=HX[n]+v[n]

其中,H=[100],v[n]=σ2,σ2为雷达量测噪声;Where, H = [100], v[n] = σ 2 , σ 2 is the radar measurement noise;

根据观测方程,确定卡尔曼滤波方程为:According to the observation equation, the Kalman filter equation is determined as:

P[n|n-1]=ΦP[n|n]ΦT+Q[n]P[n|n-1]=ΦP[n|n]Φ T +Q[n]

K[n]=P[n|n-1]HT(HP[n|n-1]HT+R)-1 K[n]=P[n|n-1]H T (HP[n|n-1]H T +R) -1

P[n|n]=(I-K[n]H)P[n|n-1]P[n|n]=(I-K[n]H)P[n|n-1]

在其中一个实施例中,根据卡尔曼滤波方程,确定发射相参参数对应的预 测发射相参参数序列为:In one embodiment, according to the Kalman filter equation, the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter is determined as:

根据预测发射相参参数序列和运动平台相参计算模型,得到相参雷达相参 参数为:According to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained as follows:

其中, in,

通过上述实施例,本发明的有益效果如下:Through the above embodiments, the beneficial effects of the present invention are as follows:

1.本发明方法是将卡尔曼预测滤波器应用于动平台分布式雷达相参合成, 解决了该情形下相参参数认知迟滞的问题,克服了基于平台运动补偿的方法, 需要已知目标的运动速度的不足。1. The method of the present invention applies the Kalman prediction filter to the coherent synthesis of the distributed radar of the moving platform, solves the problem of the coherent parameter recognition hysteresis in this case, and overcomes the deficiency of the method based on platform motion compensation that the moving speed of the target needs to be known.

2.本发明的运动建模是基于Singer模型,它考虑了目标所有机动的可能性, 能够应用于多种类型的机动,应用场景较为广泛。2. The motion modeling of the present invention is based on the Singer model, which takes into account the possibility of all maneuvers of the target and can be applied to various types of maneuvers, with a wide range of application scenarios.

3.本发明的方法具有实现简单、稳定性和普适性好等特点,利用卡尔曼滤 波和离线计算增益的特点,提高系统实时性。3. The method of the present invention has the characteristics of simple implementation, good stability and universality, and utilizes the characteristics of Kalman filtering and offline calculation gain to improve the real-time performance of the system.

以下以具体实施例进行说明。The following is a description of the specific embodiments.

本发明已经经过模拟验证。设计仿真参数为:The present invention has been verified by simulation. The design simulation parameters are:

仿真场景如图4所示:目标轨迹为有段高机动转弯,以此来检测本算法跟 踪高机动目标的性能。每部雷达之间间距5m。The simulation scenario is shown in Figure 4: The target trajectory has a high-maneuverability turn, which is used to test the performance of the algorithm in tracking high-maneuverability targets. The distance between each radar is 5m.

在开始阶段,各雷达发射正交信号,同时建立卡尔曼滤波方程,对距离进 行滤波估计,此时因信号能量弱,对于距离的量测误差较大,从图5中前20s 的距离误差值可以看出偏差还是比较大的,在状态切换时刻,各雷达发射相同 信号,此时信噪比得到提高,且利用基于singer模型的卡尔曼预测滤波器可以 较好的跟踪目标,在目标机动段最大距离误差也才4m,而图6给出了基于基本 卡尔曼预测滤波器的跟踪结果,可以看到在机动段目标量测偏差明显增加。证 明了本方法的有效性。图7和图8分别是雷达2和雷达3的发射相参参数预测 值和真实值的对比,可以看到也是在状态切换时刻之后预测值和真实值十分贴近。At the beginning, each radar transmits an orthogonal signal, and establishes a Kalman filter equation to filter and estimate the distance. At this time, due to the weak signal energy, the measurement error of the distance is large. From the distance error value of the first 20 seconds in Figure 5, it can be seen that the deviation is still relatively large. At the state switching moment, each radar transmits the same signal. At this time, the signal-to-noise ratio is improved, and the Kalman prediction filter based on the Singer model can better track the target. The maximum distance error in the target maneuvering stage is only 4m. Figure 6 shows the tracking result based on the basic Kalman prediction filter. It can be seen that the target measurement deviation increases significantly in the maneuvering stage. This proves the effectiveness of this method. Figures 7 and 8 are the comparisons of the predicted values and true values of the transmission coherent parameters of radar 2 and radar 3, respectively. It can be seen that the predicted values and true values are very close after the state switching moment.

应该理解的是,虽然图1的流程图中的各个步骤按照箭头的指示依次显示, 但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的 说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执 行。而且,图1中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些 子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行, 这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或 者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。It should be understood that, although the steps in the flowchart of FIG. 1 are shown in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear description in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, at least a part of the steps in FIG. 1 may include multiple sub-steps or multiple stages, and these sub-steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these sub-steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a part of the sub-steps or stages of other steps.

在一个实施例中,如图9所示,提供了一种动平台分布式相参雷达相参参 数获取系统,包括:场景建模模块902、卡尔曼滤波模块904和相参参数计算模 块906,其中:In one embodiment, as shown in FIG9 , a coherent parameter acquisition system for a distributed coherent radar on a moving platform is provided, comprising: a scene modeling module 902, a Kalman filter module 904 and a coherent parameter calculation module 906, wherein:

场景建模模块902,用于对动平台分布式相参雷达探测场景进行建模,得到 动平台发射信号模型、动平台发射相参模型、动平台接收相参模型;所述动平 台分布式相参雷达探测场景由多个独立的运动平台组成,所述运动平台之间相 互独立;所述运动平台之间通过无线链路传输数据;The scene modeling module 902 is used to model the moving platform distributed coherent radar detection scene to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model; the moving platform distributed coherent radar detection scene is composed of a plurality of independent moving platforms, and the moving platforms are independent of each other; the moving platforms transmit data via wireless links;

卡尔曼滤波模块904,用于根据所述动平台发射相参模型和动平台接收相参 模型,建立运动平台相参计算模型;根据发射相参参数之间的耦合关系,建立 发射相参参数的状态向量,采用Singer模型对所述状态向量进行建模,得到发 射相参参数的状态方程;根据所述状态方程,确定观测方程,以及根据所述观 测方程,确定卡尔曼滤波方程;The Kalman filter module 904 is used to establish a motion platform coherent calculation model according to the motion platform transmission coherent model and the motion platform reception coherent model; establish a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters, and use the Singer model to model the state vector to obtain a state equation of the transmission coherent parameters; determine an observation equation according to the state equation, and determine a Kalman filter equation according to the observation equation;

相参参数计算模块906,用于根据所述卡尔曼滤波方程,确定发射相参参数 对应的预测发射相参参数序列,根据所述预测发射相参参数序列和所述运动平 台相参计算模型,得到相参雷达相参参数。The coherent parameter calculation module 906 is used to determine the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter according to the Kalman filter equation, and obtain the coherent radar coherent parameters according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model.

在其中一个实施例中,场景建模模块902还用于对动平台分布式相参雷达 探测场景进行建模,得到第m个运动平台发射信号的动平台发射信号模型为:In one embodiment, the scene modeling module 902 is also used to model the moving platform distributed coherent radar detection scene, and the moving platform transmission signal model of the m-th moving platform transmission signal is obtained as:

其中,Tp为发射脉宽,u为调频斜率,rect(t)为矩形函数,为载波,sm(t)=exp(j2π(m-1)Δft);Where, Tp is the transmit pulse width, u is the frequency modulation slope, rect(t) is the rectangular function, is the carrier, s m (t) = exp(j2π(m-1)Δft);

以m表示发射信号时运动平台的序号,l表示接收信号时运动平台的序号。 则第m个运动平台发射的信号为:m represents the number of the moving platform when transmitting the signal, and l represents the number of the moving platform when receiving the signal. Then the signal transmitted by the mth moving platform is:

式中κm表示雷达m较参考时钟的同步误差,为雷达m的初始相位。Where κ m represents the synchronization error of radar m compared to the reference clock, is the initial phase of radar m.

在其中一个实施例中,场景建模模块902还用于第m个运动平台发射的信 号到达目标出的信号表示为:In one embodiment, the scene modeling module 902 is also used for the signal emitted by the m-th moving platform to reach the target. The signal is expressed as:

其中τm代表第m个运动平台发射的信号到达的时延,κm为该雷达较参考时 钟的同步误差,为该雷达较参考相位的同步误差;Where τ m represents the arrival delay of the signal emitted by the mth moving platform, κ m is the synchronization error of the radar relative to the reference clock, is the synchronization error of the radar relative to the reference phase;

则到达目标处总的信号为:The total signal reaching the target is:

设置雷达1为参考雷达,则调整后的各发射信号可以表示为:Set radar 1 as the reference radar, then the adjusted transmission signals can be expressed as:

其中为发射相参参数,得到动平台发射相参模型为:in and is the launch coherent parameter, and the launch coherent model of the moving platform is obtained as follows:

在其中一个实施例中,场景建模模块902还用于若第l个运动平台接收到目 标反射回波表示为:In one embodiment, the scene modeling module 902 is further configured to: if the l-th moving platform receives the target reflected echo, it is expressed as:

其中,p(t)是目标处的回波信号;Where p(t) is the echo signal at the target;

则所有雷达接收到目标回波叠加为:Then the target echoes received by all radars are superimposed as:

设置雷达1为参考雷达,则调整后的各接收信号表示为:Set radar 1 as the reference radar, and the adjusted received signals are expressed as:

其中,为接收相参参数,得到动平台接收相参模型为:in, and To receive the coherent parameters, the moving platform receiving coherent model is obtained as:

在其中一个实施例中,卡尔曼滤波模块904还用于根据所述动平台发射相 参模型和动平台接收相参模型,建立运动平台相参计算模型为:In one embodiment, the Kalman filter module 904 is further used to establish a motion platform coherence calculation model according to the motion platform transmission coherence model and the motion platform reception coherence model:

其中, in,

式中rl(n)是通过雷达直接量测得到的。Where r l (n) is directly measured by radar.

在其中一个实施例中,卡尔曼滤波模块904还用于根据发射相参参数之间 的耦合关系,建立发射相参参数的状态向量为:In one embodiment, the Kalman filter module 904 is further configured to establish a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters:

其中,R[n]表示状态向量,rl(n)表示发射相参参数;Where R[n] represents the state vector, r l (n) represents the transmission coherent parameter;

采用Singer模型对所述状态向量进行建模为:The state vector is modeled using the Singer model:

其中α是机动相关时间常数的倒数,即机动频率,是机动目标的加速度 方差;where α is the inverse of the maneuver-related time constant, i.e., the maneuver frequency, is the acceleration variance of the maneuvering target;

建立状态方程为:The state equation is established as:

R[n+1]=Φ(T,α)R[n]+u[n]R[n+1]=Φ(T,α)R[n]+u[n]

驱动噪声协方差为:The driving noise covariance is:

在其中一个实施例中,卡尔曼滤波模块904还用于根据所述状态方程,确 定观测方程为:In one embodiment, the Kalman filter module 904 is further configured to determine the observation equation according to the state equation:

z[n]=HX[n]+v[n]z[n]=HX[n]+v[n]

其中,H=[100],v[n]=σ2,σ2为雷达量测噪声;Where, H = [100], v[n] = σ 2 , σ 2 is the radar measurement noise;

根据所述观测方程,确定卡尔曼滤波方程为:According to the observation equation, the Kalman filter equation is determined as:

P[n|n-1]=ΦP[n|n]ΦT+Q[n]P[n|n-1]=ΦP[n|n]Φ T +Q[n]

K[n]=P[n|n-1]HT(HP[n|n-1]HT+R)-1 K[n]=P[n|n-1]H T (HP[n|n-1]H T +R) -1

P[n|n]=(I-K[n]H)P[n|n-1]P[n|n]=(I-K[n]H)P[n|n-1]

在其中一个实施例中,相参参数计算模块906还用于根据所述卡尔曼滤波 方程,确定发射相参参数对应的预测发射相参参数序列为:In one embodiment, the coherent parameter calculation module 906 is further used to determine the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter according to the Kalman filter equation:

根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参 雷达相参参数为:According to the predicted emission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained as follows:

其中, in,

关于动平台分布式相参雷达相参参数获取系统的具体限定可以参见上文中 对于动平台分布式相参雷达相参参数获取方法的限定,在此不再赘述。上述动 平台分布式相参雷达相参参数获取系统中的各个模块可全部或部分通过软件、 硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中 的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理 器调用执行以上各个模块对应的操作。The specific definition of the moving platform distributed coherent radar coherent parameter acquisition system can be found in the definition of the moving platform distributed coherent radar coherent parameter acquisition method in the above text, which will not be repeated here. Each module in the above moving platform distributed coherent radar coherent parameter acquisition system can be implemented in whole or in part by software, hardware and a combination thereof. The above modules can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory in the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程, 是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于 一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述 各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、 存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。 非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程 ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可 包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限, RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步 DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM (ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus) 直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器 总线动态RAM(RDRAM)等。A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and/or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述 实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特 征的组合不存在矛盾,都应当认为是本说明书记载的范围。The technical features of the above embodiments can be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细, 但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的 普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改 进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权 利要求为准。The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for ordinary technicians in this field, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims (9)

1.一种动平台分布式相参雷达相参参数获取方法,其特征在于,所述方法包括:1. A method for acquiring coherent parameters of a distributed coherent radar on a moving platform, characterized in that the method comprises: 对动平台分布式相参雷达探测场景进行建模,得到动平台发射信号模型、动平台发射相参模型、动平台接收相参模型;所述动平台分布式相参雷达探测场景由多个独立的运动平台组成,所述运动平台之间相互独立;所述运动平台之间通过无线链路传输数据;Modeling a moving platform distributed coherent radar detection scenario to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model; the moving platform distributed coherent radar detection scenario is composed of a plurality of independent moving platforms, and the moving platforms are independent of each other; the moving platforms transmit data via wireless links; 根据所述动平台发射相参模型和动平台接收相参模型,建立运动平台相参计算模型;Establishing a motion platform coherent calculation model according to the motion platform transmitting coherent model and the motion platform receiving coherent model; 根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量,采用Singer模型对所述状态向量进行建模,得到发射相参参数的状态方程;According to the coupling relationship between the emission coherent parameters, a state vector of the emission coherent parameters is established, and the state vector is modeled by using the Singer model to obtain a state equation of the emission coherent parameters; 根据所述状态方程,确定观测方程,以及根据所述观测方程,确定卡尔曼滤波方程;Determine an observation equation based on the state equation, and determine a Kalman filter equation based on the observation equation; 根据所述卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参数序列,根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参雷达相参参数。According to the Kalman filter equation, a predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter is determined, and according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained. 2.根据权利要求1所述的方法,其特征在于,所述对动平台分布式相参雷达探测场景进行建模,得到动平台发射信号模型,包括:2. The method according to claim 1 is characterized in that the step of modeling the moving platform distributed coherent radar detection scene to obtain the moving platform transmission signal model comprises: 对动平台分布式相参雷达探测场景进行建模,得到第m个运动平台发射信号的动平台发射信号模型为:The moving platform distributed coherent radar detection scenario is modeled, and the moving platform transmission signal model of the mth moving platform transmission signal is obtained as follows: 其中,Tp为发射脉宽,u为调频斜率,rect(t)为矩形函数,为载波,sm(t)=exp(j2π(m-1)Δft);Where, Tp is the transmit pulse width, u is the frequency modulation slope, rect(t) is the rectangular function, is the carrier, s m (t) = exp(j2π(m-1)Δft); 以m表示发射信号时运动平台的序号,l表示接收信号时运动平台的序号;则第m个运动平台发射的信号为:m represents the serial number of the motion platform when transmitting the signal, and l represents the serial number of the motion platform when receiving the signal; the signal transmitted by the mth motion platform is: 式中κm表示雷达m较参考时钟的同步误差,为雷达m的初始相位。Where κ m represents the synchronization error of radar m compared to the reference clock, is the initial phase of radar m. 3.根据权利要求2所述的方法,其特征在于,所述对动平台分布式相参雷达探测场景进行建模,得到动平台发射相参模型,包括:3. The method according to claim 2 is characterized in that the step of modeling the moving platform distributed coherent radar detection scene to obtain the moving platform emission coherent model comprises: 第m个运动平台发射的信号到达目标出的信号表示为:The signal emitted by the mth moving platform and reaching the target is expressed as: 其中τm代表第m个运动平台发射的信号到达的时延,κm为该雷达较参考时钟的同步误差,为该雷达较参考相位的同步误差;Where τ m represents the arrival delay of the signal emitted by the mth moving platform, κ m is the synchronization error of the radar relative to the reference clock, is the synchronization error of the radar relative to the reference phase; 则到达目标处总的信号为:The total signal reaching the target is: 设置雷达1为参考雷达,则调整后的各发射信号可以表示为:Set radar 1 as the reference radar, then the adjusted transmission signals can be expressed as: 其中为发射相参参数,得到动平台发射相参模型为:in and is the launch coherent parameter, and the launch coherent model of the moving platform is obtained as follows: . 4.根据权利要求3所述的方法,其特征在于,所述对动平台分布式相参雷达探测场景进行建模,得到动平台接收相参模型,包括:4. The method according to claim 3 is characterized in that the step of modeling the moving platform distributed coherent radar detection scene to obtain the moving platform receiving coherent model comprises: 若第l个运动平台接收到目标反射回波表示为:If the lth moving platform receives the target reflected echo, it can be expressed as: 其中,p(t)是目标处的回波信号;Where p(t) is the echo signal at the target; 则所有雷达接收到目标回波叠加为:Then the target echoes received by all radars are superimposed as: 设置雷达1为参考雷达,则调整后的各接收信号表示为:Set radar 1 as the reference radar, and the adjusted received signals are expressed as: 其中,为接收相参参数,得到动平台接收相参模型为:in, and To receive the coherent parameters, the moving platform receiving coherent model is obtained as: . 5.根据权利要求1至4任一项所述的方法,其特征在于,根据所述动平台发射相参模型和动平台接收相参模型,建立运动平台相参计算模型,包括:5. The method according to any one of claims 1 to 4, characterized in that, according to the moving platform transmitting coherent model and the moving platform receiving coherent model, a moving platform coherent calculation model is established, comprising: 根据所述动平台发射相参模型和动平台接收相参模型,建立运动平台相参计算模型为:According to the moving platform transmitting coherent model and the moving platform receiving coherent model, the moving platform coherent calculation model is established as follows: 其中, in, 式中rl(n)是通过雷达直接量测得到的。Where r l (n) is directly measured by radar. 6.根据权利要求1至4任一项所述的方法,其特征在于,所述根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量,采用Singer模型对所述状态向量进行建模,得到发射相参参数的状态方程,包括:6. The method according to any one of claims 1 to 4, characterized in that the state vector of the transmission coherent parameters is established according to the coupling relationship between the transmission coherent parameters, and the state vector is modeled by using the Singer model to obtain the state equation of the transmission coherent parameters, comprising: 根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量为:According to the coupling relationship between the emission coherent parameters, the state vector of the emission coherent parameters is established as: 其中,R[n]表示状态向量,rl(n)表示发射相参参数;Where R[n] represents the state vector, r l (n) represents the transmission coherent parameter; 采用Singer模型对所述状态向量进行建模为:The state vector is modeled using the Singer model: 其中α是机动相关时间常数的倒数,即机动频率,是机动目标的加速度方差;where α is the inverse of the maneuver-related time constant, i.e., the maneuver frequency, is the acceleration variance of the maneuvering target; 建立状态方程为:The state equation is established as: R[n+1]=Φ(T,α)R[n]+u[n]R[n+1]=Φ(T,α)R[n]+u[n] 驱动噪声协方差为:The driving noise covariance is: . 7.根据权利要求6所述的方法,其特征在于,根据所述状态方程,确定观测方程,以及根据所述观测方程,确定卡尔曼滤波方程,包括:7. The method according to claim 6, characterized in that determining the observation equation according to the state equation, and determining the Kalman filter equation according to the observation equation, comprises: 根据所述状态方程,确定观测方程为:According to the state equation, the observation equation is determined as: z[n]=HX[n]+v[n]z[n]=HX[n]+v[n] 其中,H=[100],v[n]=σ2,σ2为雷达量测噪声;Where, H = [100], v[n] = σ 2 , σ 2 is the radar measurement noise; 根据所述观测方程,确定卡尔曼滤波方程为:According to the observation equation, the Kalman filter equation is determined as: P[n|n-1]=ΦP[n|n]ΦT+Q[n]P[n|n-1]=ΦP[n|n]Φ T +Q[n] K[n]=P[n|n-1]HT(HP[n|n-1]HT+R)-1 K[n]=P[n|n-1]H T (HP[n|n-1]H T +R) -1 P[n|n]=(I-K[n]H)P[n|n-1]。P[n|n]=(I-K[n]H)P[n|n-1]. 8.根据权利要求7所述的方法,其特征在于,根据所述卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参数序列,根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参雷达相参参数,包括:8. The method according to claim 7 is characterized in that, according to the Kalman filter equation, a predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter is determined, and according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained, comprising: 根据所述卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参数序列为:According to the Kalman filter equation, the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter is determined as: 根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参雷达相参参数为:According to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model, the coherent radar coherent parameters are obtained as follows: 其中, in, 9.一种动平台分布式相参雷达相参参数获取系统,其特征在于,所述系统包括:9. A coherent parameter acquisition system for a distributed coherent radar on a moving platform, characterized in that the system comprises: 场景建模模块,用于对动平台分布式相参雷达探测场景进行建模,得到动平台发射信号模型、动平台发射相参模型、动平台接收相参模型;所述动平台分布式相参雷达探测场景由多个独立的运动平台组成,所述运动平台之间相互独立;所述运动平台之间通过无线链路传输数据;A scene modeling module is used to model a moving platform distributed coherent radar detection scene to obtain a moving platform transmission signal model, a moving platform transmission coherent model, and a moving platform reception coherent model; the moving platform distributed coherent radar detection scene is composed of a plurality of independent moving platforms, and the moving platforms are independent of each other; the moving platforms transmit data via wireless links; 卡尔曼滤波模块,用于根据所述动平台发射相参模型和动平台接收相参模型,建立运动平台相参计算模型;根据发射相参参数之间的耦合关系,建立发射相参参数的状态向量,采用Singer模型对所述状态向量进行建模,得到发射相参参数的状态方程;根据所述状态方程,确定观测方程,以及根据所述观测方程,确定卡尔曼滤波方程;A Kalman filter module is used to establish a motion platform coherent calculation model according to the motion platform transmission coherent model and the motion platform reception coherent model; establish a state vector of the transmission coherent parameters according to the coupling relationship between the transmission coherent parameters, and use the Singer model to model the state vector to obtain a state equation of the transmission coherent parameters; determine an observation equation according to the state equation, and determine a Kalman filter equation according to the observation equation; 相参参数计算模块,用于根据所述卡尔曼滤波方程,确定发射相参参数对应的预测发射相参参数序列,根据所述预测发射相参参数序列和所述运动平台相参计算模型,得到相参雷达相参参数。The coherent parameter calculation module is used to determine the predicted transmission coherent parameter sequence corresponding to the transmission coherent parameter according to the Kalman filter equation, and obtain the coherent radar coherent parameters according to the predicted transmission coherent parameter sequence and the motion platform coherent calculation model.
CN202110289756.2A 2021-03-17 2021-03-17 Method and system for acquiring coherent parameters of distributed coherent radar on moving platform Active CN113050053B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110289756.2A CN113050053B (en) 2021-03-17 2021-03-17 Method and system for acquiring coherent parameters of distributed coherent radar on moving platform

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110289756.2A CN113050053B (en) 2021-03-17 2021-03-17 Method and system for acquiring coherent parameters of distributed coherent radar on moving platform

Publications (2)

Publication Number Publication Date
CN113050053A CN113050053A (en) 2021-06-29
CN113050053B true CN113050053B (en) 2023-07-28

Family

ID=76513470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110289756.2A Active CN113050053B (en) 2021-03-17 2021-03-17 Method and system for acquiring coherent parameters of distributed coherent radar on moving platform

Country Status (1)

Country Link
CN (1) CN113050053B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114371461B (en) * 2022-03-21 2022-06-17 中国人民解放军国防科技大学 A method for obtaining coherent parameters of distributed coherent radar based on coherent parameter separation
CN115453508B (en) * 2022-09-21 2024-05-14 哈尔滨工业大学 Single-channel target detection and positioning system and method for maneuvering platform
CN115436905B (en) * 2022-09-30 2023-07-18 中国科学院空天信息创新研究院 Baseline estimation method, system and electronic equipment based on passive receiving system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102207548A (en) * 2010-03-31 2011-10-05 中国科学院电子学研究所 MIMO SAR imaging method by employing minimum mean square error estimation
CN109188387A (en) * 2018-08-31 2019-01-11 西安电子科技大学 Distributed coherent radar target component estimation method based on Interpolation compensation
CN110907910A (en) * 2019-11-27 2020-03-24 中国船舶重工集团公司第七二四研究所 Distributed coherent radar moving target echo coherent synthesis method
EP3712652A1 (en) * 2019-03-18 2020-09-23 NXP USA, Inc. Distributed aperture automotive radar system
CN112130139A (en) * 2020-08-21 2020-12-25 西安空间无线电技术研究所 Distributed full-coherent sparse linear array radar system optimization array arrangement method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415480A (en) * 2019-08-20 2021-02-26 是德科技股份有限公司 Multiple-input multiple-output (MIMO) target simulation system and method for testing millimeter wave radar sensor

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102207548A (en) * 2010-03-31 2011-10-05 中国科学院电子学研究所 MIMO SAR imaging method by employing minimum mean square error estimation
CN109188387A (en) * 2018-08-31 2019-01-11 西安电子科技大学 Distributed coherent radar target component estimation method based on Interpolation compensation
EP3712652A1 (en) * 2019-03-18 2020-09-23 NXP USA, Inc. Distributed aperture automotive radar system
CN110907910A (en) * 2019-11-27 2020-03-24 中国船舶重工集团公司第七二四研究所 Distributed coherent radar moving target echo coherent synthesis method
CN112130139A (en) * 2020-08-21 2020-12-25 西安空间无线电技术研究所 Distributed full-coherent sparse linear array radar system optimization array arrangement method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Using platform motion for improved spatial filtering in distributed antenna arrays;CHATTERJEE P et al.;2018 IEEE Radio and Wireless Symposium. Anaheim;全文 *
动平台分布式相参雷达系统分析;卢佳欣 等;信号处理;第35卷(第5期);全文 *
基于多特显点的无人机分布式相参雷达相位同步误差估计方法;刘晓瑜;吴建新;王彤;陈金铭;;系统工程与电子技术(04);全文 *

Also Published As

Publication number Publication date
CN113050053A (en) 2021-06-29

Similar Documents

Publication Publication Date Title
CN113050053B (en) Method and system for acquiring coherent parameters of distributed coherent radar on moving platform
CN111308437B (en) Entropy-solving and speed-ambiguity-solving method for millimeter wave MIMO traffic radar
CN110412558B (en) Method for resolving speed ambiguity of vehicle-mounted FMCW radar based on TDM MIMO
CN109100714B (en) Low-slow small target tracking method based on polar coordinate system
Rypkema et al. One-way travel-time inverted ultra-short baseline localization for low-cost autonomous underwater vehicles
CN114371461B (en) A method for obtaining coherent parameters of distributed coherent radar based on coherent parameter separation
US7804445B1 (en) Method and apparatus for determination of range and direction for a multiple tone phased array radar in a multipath environment
US11506745B2 (en) Vehicular self-positioning
CN112835005B (en) Micro Doppler feature extraction method based on super-resolution target tracking
CN105068049A (en) Split antenna MIMO radar Cramer-Rao bound calculation method
CN112147600B (en) Multi-base radar transmission parameter optimization method facing radio frequency stealth and target tracking
Ferri et al. A non-myopic, receding horizon control strategy for an AUV to track an underwater target in a bistatic sonar scenario
JP3595220B2 (en) Synthetic aperture radar device and target scattering point detection method
CN106291534A (en) A kind of track confirmation method of improvement
CN114286307A (en) Channel state information parameter estimation method based on matrix bundle
Allig et al. Heterogeneous track-to-track fusion using equivalent measurement and unscented transform
Wang et al. Robust TDOA/FDOA estimation from emitter signals for hybrid localization using UAVs
CN117388804B (en) Phased array radar system and control method thereof
Cho et al. Deep complex-valued network for ego-velocity estimation with millimeter-wave radar
Ferri et al. Results from COLLAB13 sea trial on tracking underwater targets with AUVs in bistatic sonar scenarios
CN117647797A (en) Method and device for correcting short-distance speed of radar, radar and storage medium
CN111432335B (en) Control method of international trade port system
CN112204422B (en) Speed measuring method, multiple input multiple output radar and movable platform
KR101837845B1 (en) System and method for obtaining information of underwater target
Han et al. Multi-vehicle velocity estimation using IEEE 802.11 ad waveform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant