CN102221689A - 超宽带合成孔径雷达中的多通道校正方法 - Google Patents

超宽带合成孔径雷达中的多通道校正方法 Download PDF

Info

Publication number
CN102221689A
CN102221689A CN 201110145827 CN201110145827A CN102221689A CN 102221689 A CN102221689 A CN 102221689A CN 201110145827 CN201110145827 CN 201110145827 CN 201110145827 A CN201110145827 A CN 201110145827A CN 102221689 A CN102221689 A CN 102221689A
Authority
CN
China
Prior art keywords
channel
receiving cable
expression
correcting
correction
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
CN 201110145827
Other languages
English (en)
Other versions
CN102221689B (zh
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 CN201110145827A priority Critical patent/CN102221689B/zh
Publication of CN102221689A publication Critical patent/CN102221689A/zh
Application granted granted Critical
Publication of CN102221689B publication Critical patent/CN102221689B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明提供一种超宽带合成孔径雷达中的多通道校正方法,应用于小时宽带宽积线性调频信号体制下的多通道UWB SAR/GMTI系统。技术方案包括以下步骤:第一步,利用谱修正技术构造匹配滤波器系统函数;第二步,对理想接收通道进行修正以获得参考通道;第三步,为每一个接收通道构造通道均衡器;第四步,完成通道不一致性的校正。本发明结合谱修正技术,在每一个接收通道中插入一个通道均衡器后,不仅校正了频率特性失配,完成了通道间的不一致性校正;而且有效提高了小时宽带宽积线性调频信号的脉冲压缩性能。

Description

超宽带合成孔径雷达中的多通道校正方法
技术领域
本发明属于雷达信号处理技术领域,涉及多通道UWB SAR(Ultra Wide-bandSynthetic Aperture Radar,超宽带合成孔径雷达)的通道校正方法,尤其是小时宽带宽积线性调频信号体制下的多通道UWB SAR通道校正方法。
背景技术
近十几年的研究和战争实践表明,将UWB SAR和GMTI(Ground MovingTarget Indication,地面运动目标指示)相结合更能满足现代战场侦察的需要。它能更好地捕捉实际场景,提高情报的收集能力和情报的准确性。其中,杂波抑制是运动目标检测的首要问题,而通道校正是影响杂波抑制性能的一个重要因素。因此通道校正是多通道UWB SAR/GMTI雷达系统的一个重要课题。
在多通道UWB SAR/GMTI雷达系统中,每个通道都包括阵元和馈电线路、射频放大与变频、中频处理、同相/正交支路和模拟数字变换部件,任一环节有误差都会引起各通道间的不一致性。此外,由于模拟器件及其所构成的电路不可避免地存在特性上的差异,并且在工作中模拟电路的状态也会不断的变化,所以对于众多的接收通道要求其工作特性在任何时刻完全一致是不可能的。因此为保证雷达良好的检测性能,除了严格的硬件设计和加工外,利用一定的信号处理方法可以进一步减小幅相误差,提高各通道间的一致性。
各通道间频率响应的不一致性称为通道不一致。为保证所有通道的频率响应达到一致或匹配,需要在每一个通道中插入一个通道均衡器,实现与频率有关的幅度和相位的补偿,校正频率特性失配,完成通道间的不一致性校正。
当前,对于多通道雷达系统,实现通道校正有三种方案。这些方案的基本原理和技术特点如下:
一、选择一个带内最为平坦的、畸变最小的通道作为参考通道,其余各通道与其作比较,以取得一致的频率特性。但这样处理需要一个比较选择过程,增大了运算量。
二、简单地以任意一个通道作为参考通道,在其他各个通道中插入通道均衡器,使其他所有通道都与参考通道相一致。这样N个通道需要N-1个均衡器。但这样处理仅仅保证了各个通道频率特性相一致,由于参考通道选择的“任意性”,无法保证每一个通道都具有良好的频率特性。
三、以理想接收通道作为参考通道,所有通道中全部插入均衡器,使所有通道的频率特性尽可能接近理想状态。然而,当线性调频信号的时宽带宽积小于100(此时称为小时宽带宽积线性调频信号),理想回波信号的幅度谱不具有近似矩形的特性,通带内有较严重的菲涅尔波纹,所以选取理想接收通道作为参考通道进行通道校正不能获得良好的脉冲压缩性能,不利于UWB SAR/GMTI后端系统的处理。
因此,如何对小时宽带宽积线性调频信号体制下的多通道UWB SAR/GMTI雷达系统进行通道校正是一个亟待解决的关键问题。
发明内容
本发明的目的是,提供一种小时宽带宽积线性调频信号体制下的多通道校正方法,以经过修正的理想接收通道作为参考通道,在完成多通道间不一致性校正的同时,改善回波信号的脉冲压缩性能。
本发明的技术方案是,一种超宽带合成孔径雷达中的多通道校正方法,应用于小时宽带宽积线性调频信号体制下的多通道UWB SAR/GMTI系统,其特征在于,包括以下步骤:
第一步,利用谱修正技术构造匹配滤波器系统函数。
由于时域压缩产生的距离旁瓣与频域上的边缘跃变以及带内菲涅尔波纹紧紧相关,不同时宽带宽积的线性调频信号具有不同的压缩性能。为抑制线性调频信号频谱的边缘跃变,降低距离旁瓣电平,在距离压缩时需要进行加窗处理。但小时宽带宽积信号频谱的菲涅尔波纹较为严重,影响了加窗的效果,所以与大时宽带宽积信号相比,小时宽带宽积信号脉冲压缩输出的峰值旁瓣比将会下降。因此,鉴于小时宽带宽积信号频谱的菲涅耳波纹对于距离旁瓣的突出影响,一种直观的方案是在距离压缩时,利用谱修正技术将小时宽带宽积线性调频信号频谱修正为矩形频谱。
线性调频信号频谱的修正通过修改匹配滤波器的系统函数来实现。本发明中构造的匹配滤波器系统函数M(f)为
M ( f ) = U * ( f ) I ( f ) I * ( f ) | U ( f ) | 2 - - - ( 1 )
上式中f表示频率,U(f)表示理想接收通道输出信号的频谱,I(f)=rect(f/B)表示频率宽度为B的矩形窗函数,B表示接收通道带宽,(g)*表示共轭操作,|g|2表示模值的平方。
第二步,对理想接收通道进行修正以获得参考通道。
对理想接收通道进行修正,得到修正后的参考通道Href(f)为
Href(f)=Ro(f)M(f)            (2)
上式中Ro(f)表示理想接收通道的系统函数,由雷达系统信号参数决定。
第三步,为每一个接收通道构造通道均衡器。
为补偿接收通道间的不一致性,需要在每一个接收通道中插入一个通道均衡器Em(f)
Em(f)=Href(f)/Rm(f),m=1,2,L,L    (3)
上式中Rm(f)表示第m个接收通道的系统函数,L表示多通道UWB SAR/GMTI雷达系统的通道数目。
第四步,完成通道不一致性的校正。
校正后的接收通道Hm(f)为
Hm(f)=Rm(f)Em(f),m=1,2,L,L    (4)
由公式(3)和公式(4)可知,Hm(f)=Href(f),从而所有接收通道的频率响应达到了一致,完成了多通道间不一致性的校正。
采用本发明可以取得以下技术效果:
本发明结合谱修正技术,在每一个接收通道中插入一个通道均衡器后,不仅校正了频率特性失配,完成了通道间的不一致性校正;而且有效提高了小时宽带宽积线性调频信号的脉冲压缩性能。
附图说明
图1是本发明的通道校正方法流程图;
图2是通道校正前三通道实测数据的相对幅相误差;
图3是采用本发明通道校正方法后三通道实测数据的相对幅相误差;
图4是三通道实测数据回波信号脉冲压缩输出响应波形对比图。
具体实施方式
图1是本发明的通道校正方法流程图,具体实施方案分为四步:
第一步,利用谱修正技术构造匹配滤波器系统函数M(f);
第二步,对理想接收通道进行修正以获得参考通道Href(f);
第三步,为每一个接收通道构造通道均衡器Em(f),m=1,2,L,L;
第四步,完成通道不一致性的校正。
图2至图4是在实验室利用实测数据进行实验的结果。使用的实际雷达系统的通道数为3。
图2是通道校正前三通道实测数据的相对幅相误差。失配通道相对幅相误差因子Dmn(f)定义为两个失配通道系统函数的比值,即
Figure BDA0000065477810000041
上式中dmn(f)、
Figure BDA0000065477810000042
分别表示通道校正前的相对幅度误差和相对相位误差,j为虚数单位。图2中,CH1、2表示通道1和通道2间的相对幅相误差,CH1、3表示通道1和通道3间的相对幅相误差,CH2、3表示通道2和通道3间的相对幅相误差。(a)表示通道校正前三通道间相对幅度误差dmn(f),横坐标为频率值坐标,纵坐标为幅度值坐标;(b)表示通道校正前三通道间相对相位误差
Figure BDA0000065477810000051
横坐标为频率值坐标,纵坐标为角度值坐标。从图中可以看出,通道校正前三通道间具有较大的幅度不一致性和相对固定的相位差。
图3是采用本发明通道校正方法后三通道实测数据的相对幅相误差。校正后通道相对幅相误差因子Gmn(f)定义为两个校正后通道系统函数的比值,即
G mn ( f ) = H m ( f ) / H n ( f ) = g mn ( f ) e j ψ mn ( f ) , m = 1,2 , L , Ln = 1,2 , L , L - - - ( 6 )
其中,gmn(ω)、ψmn(ω)分别表示通道校正后的相对幅度误差和相对相位误差,j为虚数单位。图3中,CH1、2表示通道1和通道2间的相对幅相误差,CH1、3表示通道1和通道3间的相对幅相误差,CH2、3表示通道2和通道3间的相对幅相误差。(a)表示通道校正后三通道间相对幅度误差gmn(ω),横坐标为频率值坐标,纵坐标为幅度值坐标;(b)表示通道校正后三通道间相对相位误差ψmn(ω),横坐标为频率值坐标,纵坐标为角度值坐标。对比图2可以看出,采用本发明通道校正方法后,三通道间的幅度不一致性几乎被消除,相位差也得到大大降低。
图4是三通道实测数据回波信号脉冲压缩输出响应波形对比图。三幅图中,横坐标均为时间坐标,纵坐标均为归一化的幅度值坐标;CH1、CH2、CH3分别表示通道1、通道2和通道3回波信号脉冲压缩输出响应波形。(a)为通道校正前三通道回波信号脉冲压缩输出响应波形;(b)为采用常规通道校正方法(以理想接收通道作为参考通道)后的三通道回波信号脉冲压缩输出响应波形;(c)为采用本发明通道校正方法后的三通道回波信号脉冲压缩输出响应波形。从(a)中可以看出波形主瓣宽度为28ns,比理想值略高,峰值旁瓣比为20dB,这说明通道校正前三通道回波信号脉压性能较差;波形的不重合反映了三通道间的不一致性。(b)中,波形较好的重合反映了三通道间较好的一致性,但较高的旁瓣反映了较差的脉压性能,这说明在小时宽带宽积线性调频信号体制下,常规的通道校正方法仅能校正通道间的不一致性,而较差的脉压性能不利于UWBSAR/GMTI后端系统的处理。(c)体现了本发明通道校正方法的优势所在,从图中可以看出波形旁瓣基本消除,峰值旁瓣比低于40dB,主瓣宽度为27ns,这说明通道校正后三通道回波信号脉压性能得到较大的改善;波形较好的重合,反映了三通道间较好的一致性。实测数据的应用结果进一步证明了本发明通道校正方法的可行性和有效性。

Claims (1)

1.一种超宽带合成孔径雷达中的多通道校正方法,应用于小时宽带宽积线性调频信号体制下的多通道UWB SAR/GMTI(Ultra Wide-band SyntheticAperture Radar/Ground Moving Target Indication,超宽带合成孔径雷达/地面运动目标指示)系统,其特征在于,包括以下步骤:
第一步,利用谱修正技术构造匹配滤波器系统函数;
构造的匹配滤波器系统函数M(f)为
Figure FDA0000065477800000011
上式中f表示频率,U(f)表示理想接收通道输出信号的频谱,I(f)=rect(f/B)表示频率宽度为B的矩形窗函数,B表示接收通道带宽,(g)*表示共轭操作,|g|2表示模值的平方;
第二步,对理想接收通道进行修正以获得参考通道;
对理想接收通道进行修正,得到修正后的参考通道Href(f)为
Href(f)=Ro(f)M(f)            (2)
上式中Ro(f)表示理想接收通道的系统函数,由雷达系统信号参数决定;
第三步,为每一个接收通道构造通道均衡器;
为每一个接收通道构造的通道均衡器Em(f)为:
Em(f)=Href(f)/Rm(f),m=1,2,L,L    (3)
上式中Rm(f)表示第m个接收通道的系统函数,L表示多通道UWB SAR/GMTI系统的通道数目;
第四步,完成通道不一致性的校正;
校正后的接收通道Hm(f)为
Hm(f)=Rm(f)Em(f),m=1,2,L,L    (4)。 
CN201110145827A 2011-06-01 2011-06-01 超宽带合成孔径雷达中的多通道校正方法 Active CN102221689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110145827A CN102221689B (zh) 2011-06-01 2011-06-01 超宽带合成孔径雷达中的多通道校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110145827A CN102221689B (zh) 2011-06-01 2011-06-01 超宽带合成孔径雷达中的多通道校正方法

Publications (2)

Publication Number Publication Date
CN102221689A true CN102221689A (zh) 2011-10-19
CN102221689B CN102221689B (zh) 2012-09-05

Family

ID=44778291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110145827A Active CN102221689B (zh) 2011-06-01 2011-06-01 超宽带合成孔径雷达中的多通道校正方法

Country Status (1)

Country Link
CN (1) CN102221689B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944885A (zh) * 2012-11-19 2013-02-27 中国人民解放军国防科学技术大学 一种卫星导航信号模拟多通道一致性方法及装置
CN103176169A (zh) * 2011-12-21 2013-06-26 中国科学院电子学研究所 利用多级均衡滤波器改善宽带有源定标器幅相特性的方法
CN103701537A (zh) * 2013-12-17 2014-04-02 电子科技大学 一种宽带接收通道综合校正法
CN104914416A (zh) * 2015-04-16 2015-09-16 南京航空航天大学 一种uhf波段sar gmti动目标筛选方法
CN106324598A (zh) * 2016-08-02 2017-01-11 北京无线电测量研究所 一种多通道子带信号的均衡、同步方法及其系统
CN106772281A (zh) * 2016-11-03 2017-05-31 上海机电工程研究所 基于滑窗滤波的多通道合成孔径雷达系统通道校准方法
WO2017173446A3 (en) * 2016-04-01 2017-11-09 Texas Instruments Incorporated Dynamic iq mismatch correction in fmcw radar
CN108267723A (zh) * 2017-01-04 2018-07-10 南京智慧魔方电子科技有限公司 陆基短距k波段雷达的离线在线式幅相误差校正方法
CN108318883A (zh) * 2018-01-16 2018-07-24 中国科学院声学研究所 一种基于全频谱换能器自然加权的距离向计算方法
CN108777671A (zh) * 2018-06-04 2018-11-09 中国电子科技集团公司第三十八研究所 一种超宽带正交解调接收机的补偿方法及装置
CN109343421A (zh) * 2018-11-15 2019-02-15 北京遥感设备研究所 一种线性调频信号的幅相补偿系统及方法
CN109387818A (zh) * 2018-10-10 2019-02-26 电子科技大学 一种基于回波数据的双基地sar通道均衡方法
CN109444883A (zh) * 2018-11-15 2019-03-08 北京遥感设备研究所 一种共口径三通道动目标检测系统
CN112578345A (zh) * 2020-11-20 2021-03-30 福瑞泰克智能系统有限公司 一种雷达遮挡检测方法、装置、设备及存储介质
US10989793B2 (en) * 2017-12-01 2021-04-27 Electromagnetic Systems, Inc. Adaptive correction of radar channel-to-channel time-dependent errors
CN113567939A (zh) * 2021-05-17 2021-10-29 陕西航天技术应用研究院有限公司 一种合成孔径雷达系统的预失真补偿方法、模块及收发系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6518914B1 (en) * 2000-11-02 2003-02-11 Totalförsvarets Forskningsinstitut Synthetic aperture radar system capable of detecting moving targets
CN201555948U (zh) * 2009-11-27 2010-08-18 中国电子科技集团公司第三十八研究所 子带合成实现超大带宽信号产生与接收系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6518914B1 (en) * 2000-11-02 2003-02-11 Totalförsvarets Forskningsinstitut Synthetic aperture radar system capable of detecting moving targets
CN201555948U (zh) * 2009-11-27 2010-08-18 中国电子科技集团公司第三十八研究所 子带合成实现超大带宽信号产生与接收系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《信号处理》 20090831 刘向阳 等 机载多通道GMTI雷达自适应通道均衡方法研究 619-622 1 第25卷, 第8A期 *
《信号处理》 20100331 张月 等 宽带数字阵列雷达通道均衡方法的设计与实现 453-457 1 第26卷, 第3期 *
《国防科技大学学报》 20100430 周红 等 基于子孔径图像的低频SAR多通道均衡方法 78-83 1 第32卷, 第4期 *

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103176169A (zh) * 2011-12-21 2013-06-26 中国科学院电子学研究所 利用多级均衡滤波器改善宽带有源定标器幅相特性的方法
CN102944885A (zh) * 2012-11-19 2013-02-27 中国人民解放军国防科学技术大学 一种卫星导航信号模拟多通道一致性方法及装置
CN102944885B (zh) * 2012-11-19 2014-07-09 中国人民解放军国防科学技术大学 一种卫星导航信号模拟多通道一致性方法及装置
CN103701537A (zh) * 2013-12-17 2014-04-02 电子科技大学 一种宽带接收通道综合校正法
CN103701537B (zh) * 2013-12-17 2016-03-30 电子科技大学 一种宽带接收通道综合校正法
CN104914416A (zh) * 2015-04-16 2015-09-16 南京航空航天大学 一种uhf波段sar gmti动目标筛选方法
CN104914416B (zh) * 2015-04-16 2018-01-23 南京航空航天大学 一种uhf波段sar gmti动目标筛选方法
WO2017173446A3 (en) * 2016-04-01 2017-11-09 Texas Instruments Incorporated Dynamic iq mismatch correction in fmcw radar
CN109154645B (zh) * 2016-04-01 2023-05-23 德州仪器公司 Fmcw雷达中的动态iq失配校正
US11460540B2 (en) 2016-04-01 2022-10-04 Texas Instruments Incorporated Dynamic IQ mismatch correction in FMCW radar
US10145937B2 (en) 2016-04-01 2018-12-04 Texas Instruments Incorporated Dynamic IQ mismatch correction in FMCW radar
CN109154645A (zh) * 2016-04-01 2019-01-04 德州仪器公司 Fmcw雷达中的动态iq失配校正
CN106324598A (zh) * 2016-08-02 2017-01-11 北京无线电测量研究所 一种多通道子带信号的均衡、同步方法及其系统
CN106772281A (zh) * 2016-11-03 2017-05-31 上海机电工程研究所 基于滑窗滤波的多通道合成孔径雷达系统通道校准方法
CN108267723A (zh) * 2017-01-04 2018-07-10 南京智慧魔方电子科技有限公司 陆基短距k波段雷达的离线在线式幅相误差校正方法
US20210215792A1 (en) * 2017-12-01 2021-07-15 Electromagnetic Systems, Inc. Adaptive Correction of Radar Channel-to-Channel Time-Dependent Errors
US10989793B2 (en) * 2017-12-01 2021-04-27 Electromagnetic Systems, Inc. Adaptive correction of radar channel-to-channel time-dependent errors
CN108318883B (zh) * 2018-01-16 2020-04-03 中国科学院声学研究所 一种基于全频谱换能器自然加权的距离向计算方法
CN108318883A (zh) * 2018-01-16 2018-07-24 中国科学院声学研究所 一种基于全频谱换能器自然加权的距离向计算方法
CN108777671B (zh) * 2018-06-04 2020-11-10 中国电子科技集团公司第三十八研究所 一种超宽带正交解调接收机的补偿方法及装置
CN108777671A (zh) * 2018-06-04 2018-11-09 中国电子科技集团公司第三十八研究所 一种超宽带正交解调接收机的补偿方法及装置
CN109387818A (zh) * 2018-10-10 2019-02-26 电子科技大学 一种基于回波数据的双基地sar通道均衡方法
CN109444883A (zh) * 2018-11-15 2019-03-08 北京遥感设备研究所 一种共口径三通道动目标检测系统
CN109343421A (zh) * 2018-11-15 2019-02-15 北京遥感设备研究所 一种线性调频信号的幅相补偿系统及方法
CN112578345A (zh) * 2020-11-20 2021-03-30 福瑞泰克智能系统有限公司 一种雷达遮挡检测方法、装置、设备及存储介质
CN113567939A (zh) * 2021-05-17 2021-10-29 陕西航天技术应用研究院有限公司 一种合成孔径雷达系统的预失真补偿方法、模块及收发系统

Also Published As

Publication number Publication date
CN102221689B (zh) 2012-09-05

Similar Documents

Publication Publication Date Title
CN102221689B (zh) 超宽带合成孔径雷达中的多通道校正方法
CN103176172B (zh) 一种机载干涉sar基于同步内定标信号的相位测量补偿方法
CN101149429B (zh) 适用于均匀圆阵的阵列互耦校正与信源测向方法
CN103616679B (zh) 基于差波束调制和波形分析的pd雷达测距测角方法
CN109471095B (zh) 一种基于快速迭代插值的fmcw雷达距离估计方法
CN103954938B (zh) 一种sar回波信号的多子带接收合成方法
CN102890271B (zh) 一种外辐射源雷达阵列天线幅相一致性校正方法
CN104267386B (zh) 一种外辐射源雷达旋转阵列的无源通道校正方法及系统
CN106338731B (zh) 一种多子带并发的mimo-sar雷达成像方法
CN102520419A (zh) 一种gnss海洋反射信号接收系统及接收方法
CN104391282A (zh) 通过频谱校正来改善成像质量的方法
CN103323822A (zh) 一种估计通道误差的方法及装置
CN103226190A (zh) 一种基于真延时的宽带相控阵波束形成方法
CN104779989A (zh) 一种宽带阵列校正滤波器系数计算方法
CN103760534B (zh) 一种星载sar数据的电离层色散效应校正方法
CN105591629A (zh) 自动频率调谐阻抗匹配的匹配频率的获取方法和装置
Vizitiu et al. The synthesis of some NLFM laws using the stationary phase principle
CN106597392B (zh) 一种sar系统相位误差补偿方法
CN105656834A (zh) 一种新型宽带接收机iq通道失配的数字校正方法
CN109946643A (zh) 基于music求解的非圆信号波达方向角估计方法
CN103701537B (zh) 一种宽带接收通道综合校正法
CN110308426B (zh) 基于数字滤波的单频连续波导航雷达零频泄露抑制方法
CN109061599B (zh) 一种基于循环平稳和对称先验知识的stap方法
CN106154257A (zh) 基于FFT与apFFT的精密测量雷达二次测频方法
CN108680907A (zh) 一种基于观测矩阵的压缩感知mimo雷达干扰抑制方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant