CN107607937A - Radar target distance-finding method based on time reversal - Google Patents

Radar target distance-finding method based on time reversal Download PDF

Info

Publication number
CN107607937A
CN107607937A CN201710720014.4A CN201710720014A CN107607937A CN 107607937 A CN107607937 A CN 107607937A CN 201710720014 A CN201710720014 A CN 201710720014A CN 107607937 A CN107607937 A CN 107607937A
Authority
CN
China
Prior art keywords
msub
mrow
signal
distance
time
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.)
Granted
Application number
CN201710720014.4A
Other languages
Chinese (zh)
Other versions
CN107607937B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201710720014.4A priority Critical patent/CN107607937B/en
Publication of CN107607937A publication Critical patent/CN107607937A/en
Application granted granted Critical
Publication of CN107607937B publication Critical patent/CN107607937B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于时间反演的雷达目标测距方法,解决了多径环境下雷达目标准确测距的问题。实现步骤是:产生雷达发射信号得到回波数据;回波数据进行时间反演,得到时反信号;设置距离搜索范围,计算不同距离点对应的虚拟信道模型,将时反信号发射到对应的虚拟信道得到时反回波信号;构造对比信号,计算各时反回波信号和对比信号相关系数;在所有相关系数中找到最大相关系数对应的距离即为目标距离。本发明对回波数据时间反演并发射到虚拟信道中,将直达波和多径波在时间上聚焦,利用多径信息使真实距离对应的时反回波和对比信号相关性较高,实现对目标的准确测距,成本较低,适用于各类延迟多径情况,用于海面低空雷达点目标距离测量。

The invention discloses a radar target distance measurement method based on time inversion, which solves the problem of accurate radar target distance measurement in a multipath environment. The implementation steps are: generate the radar transmission signal to obtain the echo data; perform time inversion on the echo data to obtain the time-inversion signal; set the distance search range, calculate the virtual channel model corresponding to different distance points, and transmit the time-inversion signal to the corresponding virtual channel model. The channel obtains the time-inverse echo signal; constructs the comparison signal, and calculates the correlation coefficient between the time-inverse echo signal and the comparison signal; finds the distance corresponding to the maximum correlation coefficient among all correlation coefficients, which is the target distance. The present invention reverses the time of the echo data and transmits it to the virtual channel, focuses the direct wave and the multipath wave in time, uses the multipath information to make the time-inverse echo corresponding to the real distance highly correlated with the comparison signal, and realizes Accurate ranging of targets, low cost, suitable for various delay multipath situations, used for distance measurement of low-altitude radar points on the sea surface.

Description

基于时间反演的雷达目标测距方法Radar Target Ranging Method Based on Time Reversal

技术领域technical field

本发明属于雷达技术领域,特别涉及目标测距,具体是一种基于时间反演的目标测距方法,可用于对海面低空点目标的测距。The invention belongs to the technical field of radar, and in particular relates to target distance measurement, in particular to a method for target distance measurement based on time inversion, which can be used for distance measurement of low-altitude point targets on the sea surface.

背景技术Background technique

雷达在海面低空环境测量目标距离时会受到多径效应的影响,即雷达收到的回波信号包括来自目标的直达波信号及来自海面反射的多径波信号。此时若使用传统的匹配滤波方法测距,雷达将同时得到多条路径的距离信息,受到多径信号的干扰,雷达将会误测出多个假目标,从而无法准确测出目标的真实距离信息。When the radar measures the distance of the target in the low-altitude environment of the sea surface, it will be affected by the multipath effect, that is, the echo signal received by the radar includes the direct wave signal from the target and the multipath wave signal reflected from the sea surface. At this time, if the traditional matched filter method is used for distance measurement, the radar will obtain the distance information of multiple paths at the same time. Due to the interference of multipath signals, the radar will mistakenly detect multiple false targets, so that the real distance of the target cannot be accurately measured. information.

针对多径效应对雷达在海面低空环境目标距离测量的影响,现有的方法主要从空域和时域两个方面抑制多径信号来获取目标的真实距离信息。在空域方面,通过改变天线放置策略,使直达波信号和多径信号几乎同时到达天线,或设计地平面天线,遮挡来自天线下方的多径信号。空域方面的两种方法需要通过设置天线实现,实现成本较高。在时域方面,使用窄相关技术,令接收机使用更宽的预相关带宽和更窄的超前滞后相关器间隔,但其使用范围有一定限制,只适合中长延迟多径的抑制。In view of the influence of multipath effect on radar target range measurement in low-altitude sea environment, the existing methods mainly suppress multipath signals from two aspects of airspace and time domain to obtain the real distance information of the target. In terms of air space, by changing the antenna placement strategy, the direct wave signal and the multipath signal arrive at the antenna almost at the same time, or the ground plane antenna is designed to block the multipath signal from below the antenna. The two methods in the airspace need to be implemented by setting up antennas, and the implementation cost is relatively high. In the time domain, the use of narrow correlation technology allows the receiver to use a wider pre-correlation bandwidth and a narrower lead-lag correlator interval, but its application range is limited, and it is only suitable for the suppression of medium and long delay multipath.

受到多径效应的影响,现有技术中的匹配滤波方法会误测出假目标,而通过空域和时域的方法对多径进行抑制,存在着成本较高、适用范围小等问题。Affected by the multipath effect, the matched filtering method in the prior art will detect false targets by mistake, but suppressing multipath through the method of space domain and time domain has problems such as high cost and limited scope of application.

发明内容Contents of the invention

本发明的目的是针对上述方法的不足,提出一种更加准确的基于时间反演的海面低空点目标的测距方法,利用多径信息对目标进行测距。The object of the present invention is to address the deficiencies of the above methods, and propose a more accurate time-reversal-based ranging method for low-altitude point targets on the sea surface, using multipath information to measure the target.

本发明是一种基于时间反演的雷达目标测距方法,其特征在于,利用多径信息对目标进行测距,包括有如下步骤:The present invention is a radar target ranging method based on time inversion, which is characterized in that the multipath information is used to measure the target, including the following steps:

(1)雷达产生发射信号st(t),得到雷达回波数据sr(t)。(1) The radar generates the emission signal s t (t), and obtains the radar echo data s r (t).

(2)将回波数据sr(t)进行时间反演处理,得到时反信号sr *(-t),其中*为共轭符号。(2) Perform time inversion processing on the echo data s r (t) to obtain the time inversion signal s r * (-t), where * is the conjugate symbol.

(3)计算得到雷达时反回波信号STR(3) Calculate the radar time back echo signal S TR ,

(4)根据发射信号st(t)和时反回波信号STR,计算皮尔逊相关系数。(4) Calculate the Pearson correlation coefficient according to the transmitted signal s t (t) and the time-inverse echo signal S TR .

4a)根据发射信号st(t)构造对比信号st *(-t)。4a) Construct a comparison signal s t * (-t) from the transmitted signal s t (t).

4b)计算得到时反回波信号STR中每个回波与对比信号st *(-t)的皮尔逊相关系数ρ。4b) Calculate the Pearson correlation coefficient ρ of each echo in the time-inverse echo signal S TR and the comparison signal s t * (-t).

(5)在全部皮尔逊相关系数ρ中找到最大的相关系数,其对应的距离即为目标的距离,完成目标测距。(5) Find the largest correlation coefficient among all Pearson correlation coefficients ρ, and its corresponding distance is the distance of the target, and complete the target ranging.

本发明与现有技术相比,具有以下优点:Compared with the prior art, the present invention has the following advantages:

一、在多径环境下,传统的匹配滤波方法将受到多径信号的干扰,无法准确测出目标的距离。本发明将回波信号时间反演并重新发送到虚拟信道中,利用时反回波信号的相关性性质,可以更准确地测出目标距离信息。1. In a multipath environment, the traditional matched filter method will be interfered by multipath signals and cannot accurately measure the distance of the target. The invention reverses the time of the echo signal and resends it to the virtual channel, and uses the correlation property of the time-inverse echo signal to measure the target distance information more accurately.

二、为减弱多径效应对测距的影响,传统方法采取抑制多径效应的方法。本发明将时反信号发射到虚拟信道中,在真实目标距离处,时反回波中直达波信号和多径信号将在时间上产生聚焦效果。本发明利用多径信息,测出目标的真实距离。本发明方法成本较低,且适用于各种延迟多径的情况。Second, in order to weaken the impact of multipath effects on ranging, the traditional method adopts the method of suppressing multipath effects. The invention transmits the time-reverse signal into the virtual channel, and at the distance of the real target, the direct wave signal and the multi-path signal in the time-reverse echo will produce a focusing effect in time. The invention uses multipath information to measure the real distance of the target. The method of the invention has low cost and is applicable to various delay multipath situations.

附图说明Description of drawings

图1是本发明的实现流程图。Fig. 1 is the realization flowchart of the present invention.

图2是在设定的目标和环境参数下,使用匹配滤波方法得到的测距结果图。Figure 2 is a graph of ranging results obtained by using the matched filtering method under the set target and environmental parameters.

图3是在设定的目标和环境参数下,用本发明方法得到的测距结果图。Fig. 3 is a distance measurement result diagram obtained by using the method of the present invention under the set target and environment parameters.

具体实施方式Detailed ways

下面结合附图,对本发明详细说明:Below in conjunction with accompanying drawing, the present invention is described in detail:

实施例1Example 1

雷达在海面低空环境中进行目标测距时,传统方法对回波进行匹配滤波处理,受到多径效应的影响,雷达对距离的检测会被多径波干扰,误测出多个假目标,无法得到真实目标的距离。针对此现状,本发明展开了研究,针对性地提出了一种成本较低,能适用于各种延迟多径的海面低空点目标的测距方法。When the radar performs target ranging in the low-altitude environment on the sea surface, the traditional method performs matched filter processing on the echo, which is affected by the multipath effect. Get the distance to the real target. Aiming at this current situation, the present invention conducts research, and specifically proposes a low-cost ranging method applicable to low-altitude point targets on the sea surface with various delays and multipaths.

本发明是一种基于时间反演的雷达目标测距方法,参照图1,利用多径信息对目标进行测距,包括有如下步骤:The present invention is a radar target ranging method based on time inversion. With reference to Fig. 1, the multipath information is used to measure the target, including the following steps:

(1)雷达产生发射信号st(t),得到雷达回波数据sr(t)。(1) The radar generates the emission signal s t (t), and obtains the radar echo data s r (t).

(2)将回波数据sr(t)进行时间反演处理,得到时反信号sr *(-t),其中*为共轭符号。(2) Perform time inversion processing on the echo data s r (t) to obtain the time inversion signal s r * (-t), where * is the conjugate symbol.

(3)建立雷达到不同距离处的虚拟信道模型并把时反信号发射到这些虚拟信道中,进行距离搜索,计算得到雷达时反回波信号STR(3) Establish the virtual channel model of the radar to different distances and transmit the time-reverse signal to these virtual channels, carry out the distance search, and calculate the radar time-reverse echo signal S TR .

(4)根据发射信号st(t)和时反回波信号STR,计算皮尔逊相关系数。(4) Calculate the Pearson correlation coefficient according to the transmitted signal s t (t) and the time-inverse echo signal S TR .

4a)根据发射信号st(t)构造对比信号st *(-t),对发射信号进行时间上的共轭反转。4a) Construct the comparison signal s t * (-t) according to the transmitted signal s t (t), and perform time conjugate inversion on the transmitted signal.

4b)计算得到时反回波信号STR中每个回波与对比信号st *(-t)的皮尔逊相关系数ρ。4b) Calculate the Pearson correlation coefficient ρ of each echo in the time-inverse echo signal S TR and the comparison signal s t * (-t).

(5)在所有皮尔逊相关系数ρ中找到最大的相关系数,其对应的距离即为目标的距离,完成目标测距。(5) Find the largest correlation coefficient among all Pearson correlation coefficients ρ, and its corresponding distance is the distance of the target, and complete the target ranging.

在海面低空环境下,若受到多径信号的影响,传统的匹配滤波方法将测出多个假目标,无法准确测出目标的真实距离。本发明将回波信号时间反演并得到时反回波信号,通过计算时反回波信号和对比信号的皮尔逊相关系数,找出最大的相关系数对应的距离,可以准确地测出目标的真实距离。In the low-altitude sea environment, if affected by multipath signals, the traditional matched filtering method will detect multiple false targets and cannot accurately measure the real distance of the target. The present invention reverses the time of the echo signal and obtains the time-inverse echo signal, and by calculating the Pearson correlation coefficient of the time-inverse echo signal and the comparison signal, finds out the distance corresponding to the largest correlation coefficient, and can accurately measure the distance of the target. real distance.

实施例2Example 2

基于时间反演的雷达目标测距方法同实施例1,本发明步骤(3)中所述计算得到雷达时反回波信号STR,具体包括:The radar target ranging method based on time inversion is the same as embodiment 1, and the reverse echo signal S TR is calculated when the radar is obtained as described in the step (3) of the present invention, specifically comprising:

3a)选取距离范围R,在距离范围R内按从小到大的顺序选取不同距离作为搜索点。3a) Select the distance range R, and select different distances in the distance range R in ascending order as search points.

3b)根据环境参数和雷达参数,建立雷达到这些搜索点处的直达波和多径波虚拟信道模型H。3b) According to the environmental parameters and radar parameters, establish the direct wave and multipath wave virtual channel model H from the radar to these search points.

3c)将时反信号sr *(-t)发射到直达波和多径波虚拟信道H中,获得对应的时反回波信号STR,用于计算和对比信号的皮尔逊相关系数。3c) Transmit the time-reverse signal s r * (-t) into the direct wave and multipath wave virtual channel H, and obtain the corresponding time-reverse echo signal S TR , which is used to calculate and compare the Pearson correlation coefficient of the signals.

其中,搜索距离范围R=[R1 … Rk … Rn],R1为起始的搜索距离,Rk为目标的真实距离,Rn为最终的搜索距离,满足R1<Rk<Rn,虚拟信道模型表示为H=[H1 … Hk … Hn],其中H1,Hk,Hn分别为R1,Rk,Rn对应的虚拟信道模型,时反回波信号STR=[sTR1(t) … sTRk(t) …sTRn(t)],其中sTR1(t),sTRk(t),sTRn(t)分别为R1,Rk,Rn对应的时反回波信号。Among them, the search distance range R=[R 1 ... R k ... R n ], R 1 is the initial search distance, R k is the real distance of the target, R n is the final search distance, satisfying R 1 <R k < R n , the virtual channel model is expressed as H=[H 1 ... H k ... H n ], where H 1 , H k , and H n are the virtual channel models corresponding to R 1 , R k , and R n respectively. Signal S TR =[s TR1 (t) … s TRk (t) …s TRn (t)], where s TR1 (t), s TRk (t), s TRn (t) are R 1 , R k , The time-reverse echo signal corresponding to R n .

通常距离范围R是针对目标真实距离Rk选取的,可以先对Rk进行粗略预估从而设置大致搜索范围。Usually the distance range R is selected for the real distance R k of the target, and R k can be roughly estimated first to set a rough search range.

实施例3Example 3

基于时间反演的雷达目标测距方法同实施例1-2,步骤3c)将时反信号sr *(-t)发射到直达波和多径波虚拟信道H中,获得对应的时反回波信号STR,包括有如下步骤:The radar target ranging method based on time inversion is the same as in embodiment 1-2, step 3c) transmit the time-reversal signal s r * (-t) into the direct wave and multipath wave virtual channel H, and obtain the corresponding time-return signal The wave signal S TR includes the following steps:

根据直达波的信道传输响应函数hk(t)及多径波的信道传输响应函数hk′(t),当搜索到目标真实距离Rk位置时,得到此时的时反回波信号为:According to the channel transmission response function h k (t) of the direct wave and the channel transmission response function h k ′(t) of the multipath wave, when the real distance R k of the target is found, the time-inverse echo signal at this time is :

其中st *(-t)为对比信号,*为共轭符号,为卷积运算,Rk对应的虚拟信道模型为Hk=[hk(t)+hk′(t)],ηk为多径波信道复电磁散射系数,τk为直达波信道延时,τk′为多径波信道延时,N(t)为高斯白噪声信号。当雷达参数、环境参数及目标位置确定时,以上参数均为可确定的。Where s t * (-t) is the contrast signal, * is the conjugate symbol, is a convolution operation, the virtual channel model corresponding to R k is H k =[h k (t)+h k ′(t)], η k is the complex electromagnetic scattering coefficient of the multipath wave channel, and τ k is the direct wave channel delay , τ k ′ is the multipath wave channel delay, and N(t) is Gaussian white noise signal. When the radar parameters, environment parameters and target position are determined, all the above parameters can be determined.

为了减弱多径效应对测距的影响,可以通过改变天线放置或使用窄相关技术抑制多径效应。本发明将回波信号进行时间反演得到时反信号,将时反信号重新发射到虚拟信道中进行距离搜索。若搜索到了目标真实距离的位置,时反回波中的直达波和多径信号将在时间上产生聚焦效果,此时的时反回波和对比信号的相关系数产生极大值。本发明利用了多径信息,促使目标真实距离处的时反回波和对比信号产生了较高的相关性,从而实现目标测距。且本发明方法适用于各种延迟多径的情况,成本较低。In order to weaken the impact of multipath effect on ranging, the multipath effect can be suppressed by changing the antenna placement or using narrow correlation technology. The invention performs time inversion on the echo signal to obtain the time inversion signal, and retransmits the time inversion signal into the virtual channel for distance search. If the position of the real distance of the target is searched, the direct wave and multipath signals in the time-reverse echo will produce a focusing effect in time, and the correlation coefficient between the time-reverse echo and the contrast signal will produce a maximum value at this time. The invention utilizes the multipath information to promote a higher correlation between the time-reverse echo and the comparison signal at the real distance of the target, thereby realizing the target distance measurement. Moreover, the method of the present invention is applicable to various delay multipath situations, and the cost is low.

实施例4Example 4

基于时间反演的雷达目标测距方法同实施例1-3The radar target ranging method based on time inversion is the same as that in Embodiment 1-3

为实现在多径环境下保证测距准确度的目的,本发明采用的技术方案包括有如下步骤:In order to realize the purpose of ensuring ranging accuracy in a multipath environment, the technical solution adopted in the present invention includes the following steps:

(1)产生雷达发射信号st(t),得到雷达回波数据sr(t)。(1) Generate the radar emission signal s t (t), and obtain the radar echo data s r (t).

(2)将回波数据sr(t)进行时间反演处理,得到时反信号sr *(-t),其中*为共轭符号。时间反演操作将回波数据在时间进行了共轭反转,这一操作使得时反回波中的直达波和多径信号在时间上产生聚焦效果。(2) Perform time inversion processing on the echo data s r (t) to obtain the time inversion signal s r * (-t), where * is the conjugate symbol. The time-reversal operation performs conjugate inversion of the echo data in time, and this operation makes the direct wave and multipath signals in the time-reversal echo produce a focusing effect in time.

(3)计算时反信号回波STR(3) Calculate the time-inverse signal echo S TR :

3a)选取距离范围R,在距离范围R内按顺序取每一个距离作为搜索点。3a) Select a distance range R, and take each distance within the distance range R in sequence as a search point.

3b)根据环境参数和雷达参数,建立雷达到这些搜索点处的直达波和多径波信道模型H。3b) According to the environmental parameters and radar parameters, establish the direct wave and multipath wave channel model H from the radar to these search points.

3c)将时反信号sr *(-t)发射到直达波和多径波虚拟信道H中,获得对应的时反信号回波STR3c) Transmit the time-reversed signal s r * (-t) into the virtual channel H of the direct wave and the multipath wave, and obtain the corresponding time-reversed signal echo S TR .

其中,距离范围R=[R1,…,Rk,…,Rn],R1为初始的搜索距离,Rk为目标的真实距离,Rn为最终的搜索距离,虚拟信道模型H=[H1 … Hk … Hn],H1,Hk,Hn分别为R1,Rk,Rn对应的虚拟信道模型,时反信号回波STR=[sTR1(t) … sTRk(t) … sTRn(t)],sTR1(t),sTRk(t),sTRn(t)分别为R1,Rk,Rn对应的时反信号回波。Among them, the distance range R=[R 1 ,...,R k ,...,R n ], R 1 is the initial search distance, R k is the real distance of the target, R n is the final search distance, and the virtual channel model H= [H 1 ... H k ... H n ], H 1 , H k , and H n are the virtual channel models corresponding to R 1 , R k , and R n respectively, and the time-reverse signal echo S TR = [s TR1 (t) ... s TRk (t) … s TRn (t)], s TR1 (t), s TRk (t), s TRn (t) are the time-reverse signal echoes corresponding to R 1 , R k , and R n respectively.

(4)根据发射信号st(t)和时反信号回波STR,测出目标的距离:(4) Measure the distance of the target according to the transmitted signal s t (t) and the time-reverse signal echo S TR :

4a)根据发射信号st(t)构造对比信号st *(-t),即对发射信号进行时域的共轭反转。4a) Construct the comparison signal s t * (-t) according to the transmitted signal s t (t), that is, perform time-domain conjugate inversion on the transmitted signal.

4b)计算时反回波信号STR中每个回波与对比信号st *(-t)的皮尔逊相关系数ρ。4b) Calculating the Pearson correlation coefficient ρ of each echo in the counter-echo signal S TR and the comparison signal s t * (-t).

任意两个信号s1(t)和s2(t)的皮尔逊相关系数为ρ满足以下关系:The Pearson correlation coefficient of any two signals s 1 (t) and s 2 (t) is ρ and satisfies the following relationship:

其中cov[s1(t),s2(t)]为信号s1(t)和s2(t)的协方差,var[s1(t)]为信号s1(t)的方差,var[s2(t)]为信号s2(t)的方差。Where cov[s 1 (t), s 2 (t)] is the covariance of signal s 1 (t) and s 2 (t), var[s 1 (t)] is the variance of signal s 1 (t), var[s 2 (t)] is the variance of the signal s 2 (t).

由于本发明对回波数据进行了时间反演操作,时反回波信号STR中真实目标距离对应的信号和对比信号相关系数最大,在所有相关系数ρ中找出最大值对应的距离即为目标的真实距离。Since the present invention performs time inversion operation on the echo data, the correlation coefficient between the signal corresponding to the real target distance and the comparison signal in the time-inverse echo signal S TR is the largest, and the distance corresponding to the maximum value is found among all correlation coefficients ρ as The true distance to the target.

下面给出一个更加详尽的例子,对本发明进一步说明:A more detailed example is given below to further illustrate the present invention:

实施例5Example 5

基于时间反演的雷达目标测距方法同实施例1-4,参照图1,本发明的具体实现步骤如下:The radar target ranging method based on time inversion is the same as embodiment 1-4, with reference to Fig. 1, the concrete realization steps of the present invention are as follows:

步骤1:产生雷达发射信号st(t),得到雷达回波数据sr(t):Step 1: Generate the radar emission signal s t (t), and obtain the radar echo data s r (t):

其中为卷积符号,hk(t)为直达波信道传输响应函数,hk′(t)为多径波信道传输响应函数,n(t)为接收机噪声信号,可看作高斯白噪声信号。τk为直达波信道延时,ηk为多径波信道的复电磁散射系数,τk′为多径波信道延时。in is the convolution symbol, h k (t) is the direct wave channel transmission response function, h k ′(t) is the multipath wave channel transmission response function, n(t) is the receiver noise signal, which can be regarded as Gaussian white noise signal . τ k is the direct wave channel delay, η k is the complex electromagnetic scattering coefficient of the multipath wave channel, and τ k ′ is the multipath wave channel delay.

步骤2:将雷达回波数据sr(t)进行时间反演处理,得到时反信号sr *(-t):Step 2: Perform time inversion processing on the radar echo data s r (t) to obtain the time inversion signal s r * (-t):

sr *(-t)=st *(-t-τk)+ηk *st *(-t-τk′)+n*(-t)s r * (-t)=s t * (-t-τ k )+η k * s t * (-t-τ k ′)+n * (-t)

其中*为共轭符号。where * is the conjugate symbol.

步骤3:计算时反回波信号STRStep 3: Calculate the time-inverse echo signal S TR :

3a)选取距离范围R,在距离范围R内按从小到大的顺序取不同距离作为搜索点。3a) Select the distance range R, and take different distances in the distance range R in order of small to large as search points.

其中,距离范围R=[R1,…,Rk,…,Rn],满足关系R1<Rk<Rn,R1为初始的搜索距离,Rk为目标的真实距离,Rn为最终的搜索距离。Among them, the distance range R=[R 1 ,...,R k ,...,R n ], satisfying the relationship R 1 <R k <R n , R 1 is the initial search distance, R k is the real distance of the target, R n is the final search distance.

3b)根据环境参数和雷达参数,建立雷达到这些搜索点处的直达波和多径波虚拟信道模型H:3b) According to the environmental parameters and radar parameters, establish the direct wave and multipath wave virtual channel model H from the radar to these search points:

H=[H1 … Hk … Hn]H=[H 1 ... H k ... H n ]

Hk=[hk(t)+hk′(t)]H k =[h k (t)+h k '(t)]

其中H1,Hk,Hn分别为R1,Rk,Rn对应的信道模型,hk(t)为直达波的信道传输响应函数,hk′(t)为多径波的信道传输响应函数。where H 1 , H k , and H n are the channel models corresponding to R 1 , R k , and R n respectively, h k (t) is the channel transmission response function of the direct wave, and h k ′(t) is the channel of the multipath wave Transfer response function.

3c)将时反信号sr *(-t)发射到直达波和多径波虚拟信道H中,获得对应的时反回波信号STR3c) Transmit the time-reverse signal s r * (-t) into the direct wave and multipath wave virtual channel H, and obtain the corresponding time-reverse echo signal S TR :

STR=[sTR1(t) … sTRk(t) … sTRn(t)]S TR =[s TR1 (t) … s TRk (t) … s TRn (t)]

其中,sTR1(t),sTRk(t),sTRn(t)分别为R1,Rk,Rn对应的时反回波信号,st *(-t)为对比信号,*为共轭符号,为卷积运算,ηk为多径波信道复电磁散射系数,τk为直达波信道延时,τk′为多径波信道延时,N(t)为噪声信号。Among them, s TR1 (t), s TRk (t), s TRn (t) are the time-inverse echo signals corresponding to R 1 , R k , and R n respectively, s t * (-t) is the comparison signal, and * is conjugate symbol, is the convolution operation, η k is the multipath wave channel complex electromagnetic scattering coefficient, τ k is the direct wave channel delay, τ k ′ is the multipath wave channel delay, N(t) is the noise signal.

步骤4:根据发射信号st(t)和时反回波信号STR,计算皮尔逊相关系数。Step 4: Calculate the Pearson correlation coefficient according to the transmitted signal s t (t) and the time-inverse echo signal S TR .

4a)根据发射信号st(t)构造对比信号st *(-t)。4a) Construct a comparison signal s t * (-t) from the transmitted signal s t (t).

4b)计算时反回波信号STR中每个回波与对比信号st *(-t)的皮尔逊相关系数ρ。4b) Calculating the Pearson correlation coefficient ρ of each echo in the counter-echo signal S TR and the comparison signal s t * (-t).

步骤5:由于对回波信号进行了时间反演并重新发射到不同的虚拟信道中,目标真实的距离对应的时反回波和对比信号的皮尔逊相关系数最大,在所有的皮尔逊相关系数ρ中找出最大的相关系数,其对应的距离即为目标的真实距离,从而实现了对目标的准确测距。Step 5: Since the echo signal is time-reversed and re-transmitted to different virtual channels, the Pearson correlation coefficient of the time-inverse echo and the comparison signal corresponding to the real distance of the target is the largest, among all the Pearson correlation coefficients Find the largest correlation coefficient in ρ, and its corresponding distance is the real distance of the target, thus realizing the accurate ranging of the target.

在多径环境下,传统的匹配滤波测距方法无法准确测出距离,会误测出假目标,而本发明通过对回波信号进行时间反演并重新发送到虚拟信道中,利用时反回波的相关性性质准确测出目标距离信息。为减弱多径效应对测距的影响,传统方法采取抑制多径效应的方法,而本发明则通过利用多径信息进行测距,且成本较低,适用于各种延迟多径的情况。In a multi-path environment, the traditional matched filter ranging method cannot accurately measure the distance, and false targets will be detected by mistake. However, the present invention performs time inversion on the echo signal and resends it to the virtual channel, and uses time inversion The correlation property of the wave accurately measures the target distance information. In order to weaken the influence of multipath effect on distance measurement, the traditional method adopts the method of suppressing multipath effect, but the present invention uses multipath information for distance measurement, and the cost is low, and it is applicable to various delayed multipath situations.

下面通过仿真实验和数据对本发明的技术效果再做说明。The technical effects of the present invention will be further described through simulation experiments and data.

实施例6Example 6

基于时间反演的雷达目标测距方法同实施例1-5,本发明对多径环境下的目标测距性能通过以下仿真进一步验证。The radar target ranging method based on time inversion is the same as in Embodiments 1-5, and the present invention further verifies the target ranging performance in a multipath environment through the following simulation.

(1)实验场景:(1) Experimental scene:

在相同的环境参数、雷达参数情况下,分别采取传统的匹配滤波测距方法和本发明基于时间反演的雷达目标测距方法进行仿真实验。Under the same environmental parameters and radar parameters, the traditional matched filter ranging method and the radar target ranging method based on time inversion of the present invention are respectively used for simulation experiments.

(2)实验内容和实验结果分析:(2) Experimental content and analysis of experimental results:

实验条件:在计算机中模拟对海面低空点目标的测距,设置环境参数、雷达参数如下表。Experimental conditions: Simulate the distance measurement of low-altitude point targets on the sea surface in the computer, and set the environmental parameters and radar parameters as shown in the table below.

在上述参数情况下,只考虑一条多径对单目标测距的影响,发射信号为水平极化方式,采用传统的匹配滤波方法进行测距,得到测距结果图如图2所示,其中横坐标为距离值,纵坐标为归一化匹配滤波幅值。从图2中可以看出,受到多径效应的影响,匹配滤波测距结果图中出现2个尖峰,这两个尖峰的幅值相差不大,无法准确对目标测距,雷达将检测到假目标。In the case of the above parameters, only the influence of one multipath on single target ranging is considered, the transmitted signal is in the horizontal polarization mode, and the traditional matched filter method is used for ranging, and the ranging result diagram is shown in Figure 2, where the horizontal The coordinates are distance values, and the ordinate is the normalized matched filter amplitude. It can be seen from Figure 2 that due to the influence of the multipath effect, two peaks appear in the matched filter ranging result graph. Target.

在同样仿真条件下,采用本发明基于时间反演的雷达目标测距方法进行测距,得到测距结果图如图3所示,其中横坐标为距离值,纵坐标为归一化相关系数。从图3中可以看出,找出相关系数最大值对应的距离值即为目标距离,图中相关系数峰值点对应的距离为30000m,测出目标距离为30000m,与目标真实距离相同,从而实现了对目标距离的准确测量,测距结果误差极小。Under the same simulation conditions, the radar target ranging method based on time inversion of the present invention is used for ranging, and the ranging result diagram is obtained as shown in Figure 3, wherein the abscissa is the distance value, and the ordinate is the normalized correlation coefficient. It can be seen from Figure 3 that the distance value corresponding to the maximum value of the correlation coefficient is the target distance. The distance corresponding to the peak point of the correlation coefficient in the figure is 30000m, and the measured target distance is 30000m, which is the same as the real distance of the target. In order to accurately measure the target distance, the error of the ranging result is extremely small.

采用传统的匹配滤波方法进行目标测距,由于多径信号的干扰,测距结果中出现多个幅值相近的尖峰,雷达无法准确测出目标的距离。采用本发明方法进行测距,测距结果中有一个尖峰的相关系数远比其他位置处的相关系数大,只需找出尖峰对应的距离,即为目标的真实距离,不受假目标的影响。且发明方法无需对天线进行设置,成本较低,适用于不同延迟多径的情况。Using the traditional matched filter method for target ranging, due to the interference of multipath signals, multiple peaks with similar amplitudes appear in the ranging results, and the radar cannot accurately measure the distance of the target. Using the method of the present invention for distance measurement, the correlation coefficient of a peak in the distance measurement result is much larger than the correlation coefficients at other positions. It is only necessary to find the distance corresponding to the peak, which is the real distance of the target, and is not affected by false targets. . Moreover, the inventive method does not need to set up the antenna, has low cost, and is applicable to situations of different delay multipaths.

综上所述,本发明公开的一种基于时间反演的雷达目标测距方法,解决了在多径环境下目标准确测距的问题。其实现步骤是:1.产生雷达发射信号并得到回波数据;2.将回波数据进行时间反演,得到时反信号;3.设置距离搜索范围,计算不同距离点对应的虚拟信道模型,将时反信号发射到这些信道中,得到时反回波信号;4.构造对比信号,计算各时反回波信号和对比信号的皮尔逊相关系数;5.在所有相关系数中找到最大相关系数对应的距离即为目标的距离。本发明通过对回波数据进行时间反演并发射到虚拟信道中,将直达波和多径波信号在时间上实现聚焦,利用多径信息使真实距离对应的时反回波和对比信号的相关性较高,从而实现了对目标的准确测距,且成本较低,适用于各类延迟多径的情况,可用于海面低空点目标距离的测量。To sum up, the present invention discloses a radar target ranging method based on time inversion, which solves the problem of accurate target ranging in a multipath environment. The implementation steps are: 1. Generate the radar transmission signal and obtain the echo data; 2. Perform time inversion on the echo data to obtain the time inversion signal; 3. Set the distance search range and calculate the virtual channel model corresponding to different distance points, Transmit the time-reverse signal into these channels to obtain the time-reverse echo signal; 4. Construct the comparison signal, and calculate the Pearson correlation coefficient of each time-reverse echo signal and the comparison signal; 5. Find the maximum correlation coefficient among all correlation coefficients The corresponding distance is the target distance. The present invention performs time inversion on the echo data and transmits them to the virtual channel to focus the direct wave and multipath wave signals in time, and uses the multipath information to correlate the time-inverse echo corresponding to the real distance with the comparison signal High reliability, so as to achieve accurate ranging of the target, and the cost is low, suitable for various delay multipath situations, and can be used to measure the distance of the target at low altitude points on the sea surface.

以上描述仅是本发明的一个具体实例,并未构成对本发明的任何限制,显然对于本领域的专业人员来说,在了解了本发明内容和原理后,都可能在不背离本发明原理、结构的情况下,进行形式和细节上的各种修改和改变,但是这些基于本发明思想的修正和改变仍在本发明的权利要求保护范围之内。The above description is only a specific example of the present invention, and does not constitute any limitation to the present invention. Obviously, for those skilled in the art, after understanding the contents and principles of the present invention, it is possible without departing from the principles and structures of the present invention. Various modifications and changes in form and details are made under the circumstances of the present invention, but these modifications and changes based on the idea of the present invention are still within the protection scope of the claims of the present invention.

Claims (3)

1. a kind of radar target distance-finding method based on time reversal, it is characterised in that surveyed using multi-path information to target Away from, including have the following steps:
(1) radar produces transmission signal st(t) radar return data s, is obtainedr(t);
(2) by echo data sr(t) time reversal processing, inverted signal s when obtaining are carried outr *(- t), wherein * are conjugate of symbol;
(3) back echo signal S when radar is calculatedTR
(4) according to transmission signal st(t) and when back echo signal STR, obtain Pearson correlation coefficient:
4a) according to transmission signal st(t) contrast signal s is constructedt *(-t);
Back echo signal S when 4b) being calculatedTRIn each echo and contrast signal st *The Pearson correlation coefficient ρ of (- t);
(5) coefficient correlation of maximum is found in whole Pearson correlation coefficient ρ, its corresponding distance is the distance of target, Complete object ranging.
2. the radar target distance-finding method according to claim 1 based on time reversal, being calculated described in step 3 Back echo signal S during radarTR, specifically include:
3a) selected distance scope R, each distance is taken to be used as Searching point in order in distance range R;
3b) according to ambient parameter and radar parameter, direct wave and multipath ripple pseudo channel of the radar at these Searching points are established Model H;
3c) by when inverted signal sr *(- t) is transmitted into direct wave and multipath ripple pseudo channel H, back echo signal when corresponding to acquisition STR
Wherein, distance range R=[R1,…,Rk,…,Rn], R1For initial detection range, RkFor the actual distance of target, RnFor Final detection range, channel model H=[H1 … Hk … Hn], wherein H1,Hk,HnRespectively R1,Rk,RnIt is corresponding virtual Channel model, when back echo signal STR=[sTR1(t) … sTRk(t) … sTRn(t)], wherein sTR1(t),sTRk(t),sTRn (t) it is respectively R1,Rk,RnBack echo signal when corresponding.
3. the radar target distance-finding method according to claim 2 based on time reversal, wherein step 3c) by when inverted signal sr *(- t) is transmitted into direct wave and multipath ripple pseudo channel H, back echo signal S when corresponding to acquisitionTR, enter as follows OK:
According to the transmission receptance function h of direct wavek(t) and multipath ripple transmission receptance function h 'k(t), when searching RkDuring position, the when back echo signal obtained now is:
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>s</mi> <mrow> <mi>T</mi> <mi>R</mi> <mi>k</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <msub> <mi>s</mi> <mi>r</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CircleTimes;</mo> <msub> <mi>H</mi> <mi>k</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msup> <msub> <mi>s</mi> <mi>r</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;CircleTimes;</mo> <mo>&amp;lsqb;</mo> <msub> <mi>h</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>h</mi> <mi>k</mi> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <msup> <msub> <mi>&amp;eta;</mi> <mi>k</mi> </msub> <mo>*</mo> </msup> <msub> <mi>&amp;eta;</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <msup> <msub> <mi>s</mi> <mi>t</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>+</mo> <msup> <msub> <mi>&amp;eta;</mi> <mi>k</mi> </msub> <mo>*</mo> </msup> <msup> <msub> <mi>s</mi> <mi>t</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>t</mi> <mo>+</mo> <msub> <mi>&amp;tau;</mi> <mi>k</mi> </msub> <mo>-</mo> <msubsup> <mi>&amp;tau;</mi> <mi>k</mi> <mo>&amp;prime;</mo> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;eta;</mi> <mi>k</mi> </msub> <msup> <msub> <mi>s</mi> <mi>t</mi> </msub> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mo>-</mo> <mi>t</mi> <mo>-</mo> <msub> <mi>&amp;tau;</mi> <mi>k</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;tau;</mi> <mi>k</mi> <mo>&amp;prime;</mo> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <mi>N</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
Wherein st *(- t) is contrast signal, and * is conjugate of symbol,For convolution algorithm, RkCorresponding channel model multipath channel models Hk=[hk(t)+h′k(t)], ηkFor multipath ripple channel electromagnetic scattering coefficient, τkFor direct wave channel time delay, τ 'kBelieve for multipath ripple Road is delayed, and N (t) is noise signal.
CN201710720014.4A 2017-08-21 2017-08-21 Radar target ranging method based on time reversal Active CN107607937B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710720014.4A CN107607937B (en) 2017-08-21 2017-08-21 Radar target ranging method based on time reversal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710720014.4A CN107607937B (en) 2017-08-21 2017-08-21 Radar target ranging method based on time reversal

Publications (2)

Publication Number Publication Date
CN107607937A true CN107607937A (en) 2018-01-19
CN107607937B CN107607937B (en) 2020-07-28

Family

ID=61065613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710720014.4A Active CN107607937B (en) 2017-08-21 2017-08-21 Radar target ranging method based on time reversal

Country Status (1)

Country Link
CN (1) CN107607937B (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031231A (en) * 2018-08-03 2018-12-18 西安电子科技大学 Low-Altitude Target Elevation time reversal is concerned with angle estimating method
CN110109092A (en) * 2019-04-24 2019-08-09 西安电子科技大学 Radar velocity measurement method under multi-path environment based on time reversal
CN110515052A (en) * 2019-06-18 2019-11-29 西安电子科技大学 A time-reversal-based method for target detection with non-equally spaced sampling in the ultra-wideband frequency domain
CN111273294A (en) * 2020-02-19 2020-06-12 西北工业大学 A time-reversal-based underwater acoustic pulse echo ranging method
CN111273296A (en) * 2020-02-27 2020-06-12 浙江工业大学 An Iterative Deconvolution-Time-Inverse Target Detection and Distance Estimation Method
CN113009467A (en) * 2021-03-09 2021-06-22 森思泰克河北科技有限公司 Radar blind area target detection tracking method and device and terminal equipment
CN113009448A (en) * 2021-03-09 2021-06-22 森思泰克河北科技有限公司 Method, device, equipment and storage medium for detecting multipath target
TWI873509B (en) * 2022-01-15 2025-02-21 大陸商華為技術有限公司 Method, apparatus and readable storage medium for distance measurement signal transmission in uwb

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120127020A1 (en) * 2010-11-18 2012-05-24 Eung Gi Paek Distributed time-reversal mirror array
CN104656078A (en) * 2013-11-25 2015-05-27 北京环境特性研究所 A Method of Sea Ship Based on Ray Tracing
CN104833961A (en) * 2015-05-06 2015-08-12 电子科技大学 Time reversal electromagnetic cancellation for active stealth
CN106226754A (en) * 2016-08-22 2016-12-14 西安电子科技大学 The low elevation angle based on time reversal Wave arrival direction estimating method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120127020A1 (en) * 2010-11-18 2012-05-24 Eung Gi Paek Distributed time-reversal mirror array
CN104656078A (en) * 2013-11-25 2015-05-27 北京环境特性研究所 A Method of Sea Ship Based on Ray Tracing
CN104833961A (en) * 2015-05-06 2015-08-12 电子科技大学 Time reversal electromagnetic cancellation for active stealth
CN106226754A (en) * 2016-08-22 2016-12-14 西安电子科技大学 The low elevation angle based on time reversal Wave arrival direction estimating method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
X. L. LU ET AL.: "Ground-penetrating radar finite-difference reverse time migration from irregular surface by flattening surface topography", 《2016 16TH INTERNATIONAL CONFERENCE ON GROUND PENETRATING RADAR (GPR)》 *
刘泽龙 等: "时间反转技术在LFMCW雷达目标检测中的应用", 《太赫兹科学与电子信息学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031231A (en) * 2018-08-03 2018-12-18 西安电子科技大学 Low-Altitude Target Elevation time reversal is concerned with angle estimating method
CN109031231B (en) * 2018-08-03 2023-02-10 西安电子科技大学 Coherence Angle Estimation Method Based on Time Retrieval of Radar Low Altitude Targets
CN110109092A (en) * 2019-04-24 2019-08-09 西安电子科技大学 Radar velocity measurement method under multi-path environment based on time reversal
CN110515052A (en) * 2019-06-18 2019-11-29 西安电子科技大学 A time-reversal-based method for target detection with non-equally spaced sampling in the ultra-wideband frequency domain
CN111273294A (en) * 2020-02-19 2020-06-12 西北工业大学 A time-reversal-based underwater acoustic pulse echo ranging method
CN111273296A (en) * 2020-02-27 2020-06-12 浙江工业大学 An Iterative Deconvolution-Time-Inverse Target Detection and Distance Estimation Method
CN113009467A (en) * 2021-03-09 2021-06-22 森思泰克河北科技有限公司 Radar blind area target detection tracking method and device and terminal equipment
CN113009448A (en) * 2021-03-09 2021-06-22 森思泰克河北科技有限公司 Method, device, equipment and storage medium for detecting multipath target
CN113009448B (en) * 2021-03-09 2022-12-06 森思泰克河北科技有限公司 Method, device, equipment and storage medium for detecting multipath target
TWI873509B (en) * 2022-01-15 2025-02-21 大陸商華為技術有限公司 Method, apparatus and readable storage medium for distance measurement signal transmission in uwb

Also Published As

Publication number Publication date
CN107607937B (en) 2020-07-28

Similar Documents

Publication Publication Date Title
CN107607937B (en) Radar target ranging method based on time reversal
CN107942322A (en) Multipath utilizes radar time reversal object distance estimation method
CN109188470B (en) GNSS deception jamming detection method based on convolutional neural network
CN109031231B (en) Coherence Angle Estimation Method Based on Time Retrieval of Radar Low Altitude Targets
CN102156279B (en) Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output)
CN105807267B (en) A kind of MIMO radar extends mesh object detection method
CN103698759B (en) A kind of external radiation source radar system based on SFN and signal processing method thereof
US7468689B2 (en) System and method for determining position of radar apparatus based on reflected signals
CN109581317B (en) Corner target positioning method based on echo peak matching
CN113960558B (en) Non-line-of-sight target positioning method and system based on multiple-input multiple-output radar
CN109239683B (en) Clarmero boundary analysis method of broadband passive MIMO radar
CN109490845B (en) Method for suppressing main lobe suppressing interference by multi-station radar
CN104062638A (en) Multi-target detecting method used for through-wall radar
CN104678386B (en) Method for detecting target by utilizing relevant power of GNSS(global navigation satellite system) sea surface reflection signal
CN104155651A (en) Probability data correlation method for polarizing radar target tracking
CN104730498A (en) Target detection method based on Keystone and weighting rotating FFT
CN110109092B (en) Radar speed measurement method based on time reversal in multipath environment
CN107918115A (en) The radar target localization method utilized based on multipath
CN103197297A (en) Radar moving target detection method based on cognitive framework
CN111107626A (en) DOA positioning method based on time reversal
CN103135109B (en) Ultra wide band radar imaging method based on multipath signals
CN103487796B (en) A kind of method utilizing underwater acoustic channel Statistically invariant feature to realize passive ranging
CN104243064B (en) A kind of spread spectrum TOA distance-finding method under non line of sight multi-path dense scene
CN105891800B (en) Single-shot three receives antenna radar anti-multipath scattering detection method
CN104360333A (en) Phase-coherent accumulation detecting method capable of revising first-order and second-order range migration at the same time

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