CN103852648A - 获取空间电磁强度数据的方法 - Google Patents
获取空间电磁强度数据的方法 Download PDFInfo
- Publication number
- CN103852648A CN103852648A CN201410103112.XA CN201410103112A CN103852648A CN 103852648 A CN103852648 A CN 103852648A CN 201410103112 A CN201410103112 A CN 201410103112A CN 103852648 A CN103852648 A CN 103852648A
- Authority
- CN
- China
- Prior art keywords
- signal
- radar
- frequency
- communication
- transmitter
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种获取空间电磁强度数据的方法,主要解决现有技术运算速度缓慢,造成电子设备无法正常工作或目标干扰大的问题,其实现步骤是:1)确定雷达发射机、通信发射机的数量和参数,并分别计算其采样频率和采样点数;2)分别对雷达信号和通信信号进行采样,并对采样信号进行傅里叶变换得到相应信号频谱;3)设置观测点;4)通过对信号频谱运算,得到观测点雷达信号叠加频谱和通信信号叠加频谱;5)根据叠加频谱计算信号总功率,进而得到电磁强度;6)对整个观测区域重复步骤3至步骤5,获得观测区域的电磁强度数据。本发明大幅度减小了运算量,提高了电子设备反应速度,可用于分析电磁环境,避免各种电磁设备相互干扰。
Description
技术领域
本发明属于空间电磁态势感知技术,具体涉及一种获取空间电磁强度数据的方法。可用于分析电磁环境,避免各种电磁设备相互干扰。
背景技术
随着信息化的发展,现实环境中的电磁设备越来越多,造成电磁环境异常复杂,如何得到电磁环境中的功率,电磁强度数据对研究电磁环境具有重要价值:一是有助于避免各电磁设备间的相互干扰;二是可引导民航飞机在运动过程中避免强电磁干扰,以使自身安全起飞和降落。
传统获取空间电磁强度数据的方法是:在观测点处先将所有经过该点的时域信号叠加,再对叠加信号进行频谱变换获得空间电磁强度数据。该方法需要在每个观测点处均进行一次运算量较大的频谱变换计算,在有大量观测点的情况下,运算速度十分缓慢,将造成电子设备无法正常工作或飞机无法及时躲避干扰等情况。
发明内容
本发明的目的在于针对上述已有技术的不足,提出一种获取空间电磁强度数据的方法,以大幅度减小运算量,提高电子设备反应速度。
为实现上述目的,本发明的实现步骤包括如下:
1)在空间放置M台雷达发射机和N台通信发射机,其中,M≥1,N≥1,雷达发射机的功率为:Prm>0,雷达发射机的天线增益为:Grm>0,通信发射机的功率为:Pcn>0,通信发射机的天线增益为:Gcn>0,雷达发射机的发射信号为srm(t),通信发射机的发射信号为scn(t),其中,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N;
2)计算雷达发射信号的采样频率fr和采样点数Kr,通信发射信号的采样频率fc和采样点数Kc;
3)计算雷达发射信号频谱Srm(f)和通信发射信号频谱Scn(f):
对各雷达发射机信号均以采样频率fr进行采样,并对采样信号进行傅里叶变换得到雷达发射信号频谱Srm(f);
对各通信发射机信号均以采样频率fc进行采样,并对采样信号进行傅里叶变换得到通信发射信号频谱Scn(f);
4)设置一个观测点,计算观测点处各雷达信号的幅值Arm和各通信信号的幅值Acn:
其中,σ表示空间观测点的截面积,Rrm和Rcn分别表示雷达发射机和通信发射机到该观测点的距离;
5)计算观测点处的雷达信号叠加频谱Sr(f)和通信信号叠加频谱Sc(f);
其中,τrm=Rrm/c表示观测点与雷达发射机间的时间延迟,c表示光速;
其中,τcn=Rcn/c表示观测点与通信发射机间的时间延迟;
6)根据观测点处的雷达信号叠加频谱Sr(f)和通信信号的叠加频谱Sc(f),计算观测点处信号总功率P:
6a)根据雷达信号叠加频谱Sr(f)每个频点上的谱线幅值Hr(k),计算观测点处雷达信号功率Pr;
6b)根据通信信号叠加频谱Sc(f)每个频点上的谱线幅值Hc(k),计算观测点处通信信号功率Pc;
6c)将观测点处雷达信号功率Pr与通信信号功率Pc相加,得到观测点处信号总功率P;
8)对整个观测区域的其它观测点重复步骤4)至步骤7),即可获得观测区域的空间电磁强度数据。
本发明与现有技术相比具有以下优点:
1)本发明采用全采样方法对中心频率较低的通信信号进行采样,而对中心频率较高的雷达信号,仅仅取其带宽内的信号进行采样,避免了对每个发射信号进行一次全采样,因而很大程度上减轻了系统的负担。
2)本发明事先对所有发射信号做频域变换得到发射信号的频谱,在每个观测点处,只需要将发射信号频谱进行一系列运算量较小的乘法和加法运算即可获得叠加频谱,避免了在每个观测点处均进行一次运算量较大的频域变换,大大提高了电子设备的运算速率,节省了运算时间,使电子装备更具灵活性。
附图说明
图1是本发明的实现流程图;
图2是观测点数与运行时间的关系图;
图3是观测区域电磁强度平面分布图;
图4是观测区域电磁强度三维分布图。
具体实施方式
参照图1,本发明的具体步骤如下:
步骤1:确定雷达发射机和通信发射机的数量及相关参数。
在空间放置M台雷达发射机和N台通信发射机,其中,M≥1,N≥1,雷达发射机的功率为:Prm>0,雷达发射机的天线增益为:Grm>0,通信发射机的功率为:Pcn>0,通信发射机的天线增益为:Gcn>0,雷达发射机的发射信号为srm(t),通信发射机的发射信号为scn(t),其中,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N。
步骤2:计算雷达发射信号的采样频率fr和采样点数Kr。
2a)计算雷达发射信号采样频率fr:
2a1)根据各雷达发射信号的中心频率frm和带宽Brm,计算各雷达发射信号的起始频率frmL和终止频率frmH:
frmL=frm-Brm/2
frmH=frm+Brm/2
2a2)取所有雷达发射信号的起始频率中的最小值,终止频率中的最大值,得到雷达发射信号的最低频率frL和最高频率frH:
frL=min{frmL}, frH=max{frmH};
2a3)根据雷达发射信号的最低频率frL和最高频率frH,计算中间变量o:
2a4)根据雷达发射信号的最高频率frH和中间变量o,计算雷达发射信号的采样频率fr:
fr=2frH/(o+1);
2b)计算雷达发射信号采样点数Kr:
为了保证得到的离散信号中至少包含一个完整的脉冲周期,故将雷达发射信号的采样点数Kr按下式计算:
Kr=1.2frmax{trm},
其中,trm为各雷达发射信号的脉冲重复周期;为了能够进行快速傅里叶变换,雷达发射信号的采样点数Kr必须是2的幂次方,如果Kr不是2的幂次方,则将Kr增大到2的幂次方。
步骤3:计算通信发射机信号的采样频率fc和采样点数Kc。
3a)根据各通信发射信号的中心频率fcn和带宽Bcn,计算各通信发射信号的终止频率fcnH:
fcnH=fcn+Bcn/2;
3b)取所有通信发射信号终止频率中的最大值,得到通信发射信号的最高频率fcH:
fcH=max{fcnH};
3c)通信发射信号采样点数Kc:
根据信号频谱叠加的要求,雷达信号频谱和通信信号频谱需要具有相同的采样率,即fc/Kc=fr/Kr;根据奈奎斯特采样定理,通信发射信号采样频率fc应大于等于信号最高频率fcH的2倍,即fc≥2×fcH;
由以上两式得到不等式Kc×(fr/Kr)≥2fcH,解此不等式得到通信发射信号采样点数Kc的取值范围:
为了能够进行快速傅里叶变换,Kc必须是2的幂次方,因此,在Kc的取值范围内,Kc取2的幂次方的最小整数;
3d)根据通信发射信号的采样点数Kc,计算通信发射信号的采样频率fc:
fc=Kc×(fr/Kr)。
步骤4:计算雷达发射信号频谱Srm(f)和通信发射信号频谱Scn(f)。
4a)对各雷达发射机信号均以采样频率fr进行采样,并对采样信号进行傅里叶变换得到雷达发射信号频谱Srm(f);
4b)对各通信发射机信号均以采样频率fc进行采样,并对采样信号进行傅里叶变换得到通信发射信号频谱Scn(f)。
步骤5:计算观测点处雷达信号的幅值Arm和观测点处通信信号的幅值Acn。
5a)在观测区域确定一个观测点,其位置为(x,y,z),其中x,y,z分别表示经度、纬度和高度,分别计算观测点到各雷达发射机的距离Rrm,观测点到各通信发射机的距离Rcn:
其中,(xrm,yrm,zrm)表示第m个雷达发射机的位置坐标,(xcn,ycn,zcn)表示第n个通信发射机的位置坐标;
5b)计算观测点处雷达信号的幅值Arm和观测点处通信信号的幅值Acn:
其中,σ表示空间观测点的截面积,σ取值为1。
步骤6:计算观测点处的雷达信号叠加频谱Sr(f)和通信信号叠加频谱Sc(f)。
根据傅里叶变换的基本性质:
其中,τrm=Rrm/c表示观测点与各雷达发射机之间的时间延迟,c表示光速,
其中,τcn=Rcn/c表示观测点到各通信发射机间的时间延迟;
6b)计算观测点处所有雷达信号的叠加频谱Sr(f)和所有通信信号的叠加频谱Sc(f):
步骤7:计算观测点处信号总功率P。
7a)根据雷达信号叠加频谱Sr(f)每个频点上的谱线幅值Hr(k),计算观测点处雷达信号功率Pr:
7b)根据通信信号叠加频谱Sc(f)每个频点上的谱线幅值Hc(k),计算观测点处通信信号功率Pc:
7c)将观测点处雷达信号功率Pr与通信信号功率Pc相加,得到观测点处信号总功率P:
P=Pr+Pc。
步骤8:获得观测点处的场强E。
8a)根据观测点处信号总功率P和观测点处的面积σ,得出功率谱密度D
D=P/σ;
8b)根据功率谱密度D,由式得出观测点处场强。
步骤9:对整个观测区域的其它观测点重复步骤4)至步骤8)即可获得观测区域的空间电磁强度数据。
本发明的效果可以通过以下仿真结果进一步说明:
1.仿真条件:
空间有四部不同的装备,即雷达装备、通信装备、雷达干扰装备、通信干扰装备各一个,位置坐标分别为(115°,15°,100m),(124°,15°,300m),(124°,24°,50m),(115°,24°,100m)。它们的发射信号的中心频率、带宽、天线增益均相同。其中发射信号中心频率为2500MHz,带宽为10MHz,天线增益为25dB,发射功率分别为50kw,1kw,150kw,20kw。
观测区域经度范围:115°~124°,纬度范围:15°~24°,高度:10000m。
采用本发明在MatlabR2010a软件平台进行仿真实验。
2.仿真内容与结果:
仿真1,采用本发明和现有的时域叠加求场强法,分别获取各观测点的电磁强度数据,并统计软件运行时间,绘制观测点数与运行时间的关系如图2所示。图2中横坐标表示观测点数,纵坐标表示运行时间。
从图2可以看出,本发明与现有时域叠加求场强法相比,大幅度降低了运算时间,且随着观测点数的增加,运行时间降低幅度越大。
仿真2,采用本发明和现有的时域叠加求场强法,分别获取观测区域的各点的电磁强度数据,根据观测点位置坐标和对应该点的电磁强度数据,分别绘制电磁强度平面分布如图3和电磁强度三维分布如图4。
由图3、4可见,本发明和现有的时域叠加求场强法得到的磁强度分布是相同的,证明本发明获得的电磁强度数据是正确的。
综上,本发明在保证观测点电磁强度计算准确的同时,大幅度提高运算速度,不仅使电子装备更具灵活性,还有助于更迅速的感知电磁环境,做出合理决策。
Claims (5)
1.一种获取空间电磁强度数据的方法,包括以下步骤:
1)在空间放置M台雷达发射机和N台通信发射机,其中,M≥1,N≥1,雷达发射机的功率为:Prm>0,雷达发射机的天线增益为:Grm>0,通信发射机的功率为:Pcn>0,通信发射机的天线增益为:Gcn>0,雷达发射机的发射信号为srm(t),通信发射机的发射信号为scn(t),其中,m表示雷达发射机序号,m=1,2,...,M,n表示通信发射机序号,n=1,2,...,N;
2)计算雷达发射信号的采样频率fr和采样点数Kr,通信发射信号的采样频率fc和采样点数Kc;
3)计算雷达发射信号频谱Srm(f)和通信发射信号频谱Scn(f):
对各雷达发射机信号均以采样频率fr进行采样,并对采样信号进行傅里叶变换得到雷达发射信号频谱Srm(f);
对各通信发射机信号均以采样频率fc进行采样,并对采样信号进行傅里叶变换得到通信发射信号频谱Scn(f);
4)设置一个观测点,计算观测点处各雷达信号的幅值Arm和各通信信号的幅值Acn:
其中,σ表示空间观测点的截面积,Rrm和Rcn分别表示雷达发射机和通信发射机到该观测点的距离;
5)计算观测点处的雷达信号叠加频谱Sr(f)和通信信号叠加频谱Sc(f);
其中,τrm=Rrm/c表示观测点与雷达发射机间的时间延迟,c表示光速;
其中,τcn=Rcn/c表示观测点与通信发射机间的时间延迟;
6)根据观测点处的雷达信号叠加频谱Sr(f)和通信信号的叠加频谱Sc(f),计算观测点处信号总功率P:
6a)根据雷达信号叠加频谱Sr(f)每个频点上的谱线幅值Hr(k),计算观测点处雷达信号功率Pr;
6b)根据通信信号叠加频谱Sc(f)每个频点上的谱线幅值Hc(k),计算观测点处通信信号功率Pc;
6c)将观测点处雷达信号功率Pr与通信信号功率Pc相加,得到观测点处信号总功率P;
8)对整个观测区域的其它观测点重复步骤4)至步骤7),即可获得观测区域的空间电磁强度数据。
2.根据权利要求1所述方法,其中所述步骤2)中计算所有雷达发射信号的采样频率fr和采样点数Kr,按如下步骤进行:
2a)计算雷达发射信号的采样频率fr:
2a1)根据各雷达发射信号的中心频率frm和带宽Brm,计算各雷达发射信号的起始频率frmL和终止频率frmH:
2a2)取所有雷达发射信号的起始频率中的最小值,终止频率中的最大值,得到雷达发射信号的最低频率frL和最高频率frH:
frL=min{frmL}, frH=max{frmH};
2a3)根据雷达发射信号的最低频率frL和最高频率frH,计算中间变量o:
2a4)根据雷达发射信号的最高频率frH和中间变量o,计算雷达发射信号的采样频率fr:
fr=2frH/(o+1);
2b)根据雷达发射信号的采样频率fr,得出雷达发射信号的采样点数Kr:
Kr=1.2frmax{trm},
其中,trm表示各雷达发射信号的脉冲重复周期,如果采样点数Kr不是2的幂次方,则将Kr增大到2的幂次方。
3.根据权利要求1所述方法,其中所述步骤2)中计算通信发射信号的采样频率fc和采样点数Kc,按如下步骤进行:
2.1)根据各通信发射信号的中心频率fcn和带宽Bcn,计算各通信发射信号的终止频率fcnH:
fcnH=fcn+Bcn/2,
2.2)取所有通信发射信号终止频率中的最大值,得到通信发射信号的最高频率fcH:
fcH=max{fcnH}
2.3)根据通信发射信号的最高频率fcH,雷达发射信号的采样频率fr和采样点数Kr,得到通信发射信号采样点数Kc的取值范围:
在此范围内,Kc取2的幂次方的最小整数;
2.4)根据通信发射信号的采样点数Kc,计算通信发射信号的采样频率fc:
fc=Kc×(fr/Kr)。
4.根据权利要求1所述方法,其中步骤6a)所述的计算观测点处雷达信号功率Pr,按如下公式计算:
其中,Hr(k)表示雷达信号叠加频谱Sr(f)每个频点上的谱线幅值,Kr表示雷达发射信号的采样点数。
5.根据权利要求1所述方法,其中步骤6b)所述的计算观测点处通信信号功率Pc,按如下公式计算:
其中,Hc(k)表示通信信号叠加频谱Sc(f)每个频点上的谱线幅值,Kc表示通信发射信号采样点数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410103112.XA CN103852648B (zh) | 2014-03-20 | 2014-03-20 | 获取空间电磁强度数据的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410103112.XA CN103852648B (zh) | 2014-03-20 | 2014-03-20 | 获取空间电磁强度数据的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103852648A true CN103852648A (zh) | 2014-06-11 |
CN103852648B CN103852648B (zh) | 2016-07-06 |
Family
ID=50860549
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410103112.XA Expired - Fee Related CN103852648B (zh) | 2014-03-20 | 2014-03-20 | 获取空间电磁强度数据的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103852648B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954586A (zh) * | 2016-07-15 | 2016-09-21 | 西安电子科技大学 | 基于变采样率的空间电磁信号频谱获取方法 |
CN105974283A (zh) * | 2016-05-18 | 2016-09-28 | 东北电力大学 | 基于小波包生存指数奇异熵的电缆局放特征提取方法 |
CN106291126A (zh) * | 2016-07-15 | 2017-01-04 | 西安电子科技大学 | 基于变采样率的复杂电磁环境时域信号模拟方法 |
CN109521407A (zh) * | 2018-12-20 | 2019-03-26 | 陕西黄河集团有限公司 | 雷达发射分系统工作带宽测试方法及系统 |
CN112711032A (zh) * | 2020-12-07 | 2021-04-27 | 中国人民解放军海军航空大学 | 一种基于图数据和gcn的雷达目标检测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5122748A (en) * | 1990-08-13 | 1992-06-16 | The Trustees Of Columbia University In The City Of New York | Method and apparatus for spatial localization of magnetic resonance signals |
CN1618023A (zh) * | 2002-01-18 | 2005-05-18 | 加拿大工业部 | 用于测量复杂电磁场的天线阵列 |
US20050206384A1 (en) * | 2004-03-22 | 2005-09-22 | Hideta Habara | Radiofrequency wave transmission circuit using arbitray waveform generator and NMR apparatus using the same |
CN103018759A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种基于场强值的同频多辐射源定位方法 |
-
2014
- 2014-03-20 CN CN201410103112.XA patent/CN103852648B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5122748A (en) * | 1990-08-13 | 1992-06-16 | The Trustees Of Columbia University In The City Of New York | Method and apparatus for spatial localization of magnetic resonance signals |
CN1618023A (zh) * | 2002-01-18 | 2005-05-18 | 加拿大工业部 | 用于测量复杂电磁场的天线阵列 |
US20050206384A1 (en) * | 2004-03-22 | 2005-09-22 | Hideta Habara | Radiofrequency wave transmission circuit using arbitray waveform generator and NMR apparatus using the same |
CN103018759A (zh) * | 2012-11-23 | 2013-04-03 | 北京航空航天大学 | 一种基于场强值的同频多辐射源定位方法 |
Non-Patent Citations (1)
Title |
---|
蒲鋆: "空间电磁辐射场强的建模与仿真", 《中国优秀硕士学位论文全文数据库·信息科技辑》, no. 07, 15 July 2011 (2011-07-15), pages 136 - 834 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105974283A (zh) * | 2016-05-18 | 2016-09-28 | 东北电力大学 | 基于小波包生存指数奇异熵的电缆局放特征提取方法 |
CN105954586A (zh) * | 2016-07-15 | 2016-09-21 | 西安电子科技大学 | 基于变采样率的空间电磁信号频谱获取方法 |
CN106291126A (zh) * | 2016-07-15 | 2017-01-04 | 西安电子科技大学 | 基于变采样率的复杂电磁环境时域信号模拟方法 |
CN106291126B (zh) * | 2016-07-15 | 2018-11-20 | 西安电子科技大学 | 基于变采样率的复杂电磁环境时域信号模拟方法 |
CN109521407A (zh) * | 2018-12-20 | 2019-03-26 | 陕西黄河集团有限公司 | 雷达发射分系统工作带宽测试方法及系统 |
CN112711032A (zh) * | 2020-12-07 | 2021-04-27 | 中国人民解放军海军航空大学 | 一种基于图数据和gcn的雷达目标检测方法及系统 |
CN112711032B (zh) * | 2020-12-07 | 2022-05-13 | 中国人民解放军海军航空大学 | 一种基于图数据和gcn的雷达目标检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103852648B (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103852648A (zh) | 获取空间电磁强度数据的方法 | |
CN103869311B (zh) | 实波束扫描雷达超分辨成像方法 | |
CN103023586B (zh) | 一种天波超视距雷达电离层信道仿真方法 | |
CN104614713A (zh) | 一种适合于艇载雷达系统的雷达回波信号模拟器 | |
CN102866388B (zh) | 一种空时自适应处理中的自适应权值迭代计算方法 | |
CN104199022B (zh) | 一种基于目标模态估计的临近空间高超声速目标跟踪方法 | |
CN103869298B (zh) | 一种分布式mimo天波超视距雷达海杂波仿真方法 | |
CN110174643A (zh) | 一种无需噪声功率信息的基于到达时间差的定位方法 | |
CN103048651B (zh) | 多参数模拟气象雷达回波发生装置及发生方法 | |
CN104035095A (zh) | 基于空时最优处理器的低空风切变风速估计方法 | |
CN103777186B (zh) | 基于波束分解和局部照射下运动目标近场雷达回波特性计算方法 | |
KR102013205B1 (ko) | 레이더 신호 처리 모의 장치 및 방법 | |
JP6346082B2 (ja) | パルス圧縮レーダ装置及びそのレーダ信号処理方法 | |
CN106291489B (zh) | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 | |
CN107607943A (zh) | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 | |
CN106405533A (zh) | 基于约束加权最小二乘的雷达目标联合同步与定位方法 | |
CN106597425A (zh) | 一种基于机器学习的雷达目标定位方法 | |
KR101828332B1 (ko) | 차량용 레이더 시뮬레이터 | |
CN113949439A (zh) | 一种面向无人机空对空通信的几何随机信道建模方法 | |
CN103267964A (zh) | 一种基于低秩矩阵恢复的弹载导引头雷达σδ-stap方法 | |
CN103699810A (zh) | 一种粗糙面微波段双向反射分布函数的建模方法 | |
CN103163509B (zh) | 一种基于电磁散射的高频近似法的模拟合成孔径雷达 | |
CN104101869B (zh) | 一种极坐标下的地波雷达运动目标仿真建模方法 | |
CN105372637A (zh) | 基于运动车辆目标微动特性的对雷达的窄带欺骗干扰方法 | |
CN105548959B (zh) | 一种基于稀疏重建的多传感器多目标的定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160706 Termination date: 20210320 |