CN114563002B - 一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 - Google Patents
一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 Download PDFInfo
- Publication number
- CN114563002B CN114563002B CN202210274958.4A CN202210274958A CN114563002B CN 114563002 B CN114563002 B CN 114563002B CN 202210274958 A CN202210274958 A CN 202210274958A CN 114563002 B CN114563002 B CN 114563002B
- Authority
- CN
- China
- Prior art keywords
- data
- processing
- flexible gyroscope
- vehicle
- signal
- 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
Links
- 238000012545 processing Methods 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims description 19
- 238000005259 measurement Methods 0.000 claims abstract description 25
- 238000012935 Averaging Methods 0.000 claims abstract description 11
- 238000001914 filtration Methods 0.000 claims abstract description 10
- 238000003672 processing method Methods 0.000 claims abstract description 7
- 125000004122 cyclic group Chemical group 0.000 claims abstract description 4
- 238000010606 normalization Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 238000012360 testing method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/18—Stabilised platforms, e.g. by gyroscope
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Algebra (AREA)
- Life Sciences & Earth Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Computational Biology (AREA)
- Gyroscopes (AREA)
Abstract
本发明涉及一种在车载扰动环境下进行寻北测量的陀螺信号处理方法,包括下列步骤,S1:DSP响应2000Hz定时中断,对陀螺数据进行8次循环采样;S2:对8次数据进行由小到大排序,排序后取中间的4个数据进行平均;S3:对数据进行二阶低通滤波处理,滤波器截止频率166Hz;S4:对数据进行陷波处理,陷波器截止频率346Hz;S5:对数据进行滞后补偿处理;S6:对S5产生的数据按2000H解算周期,重复采集1000次,对1000次数据求和,也就是0.5s周期内的数据记为一个数据,存入数组;数组共存入30个数据,也就是15s周期数据;S7:对30个数据进行处理,先对数据进行由小到大排序,排序后取中间的10个数据进行平均;S8:对处理后的10个数据平均,进行数据归一化处理;数据代入惯性测量算法,解算测量值。
Description
技术领域
本发明涉及一种惯性测量方法,尤其是涉及一种适应车载扰动环境下寻北的挠性陀螺信号处理方法。
背景技术
挠性陀螺是惯性测量装置的重要元件之一,其用于测量载体的角速度,为惯性测量装置提供精确的角速度信息,挠性陀螺的性能指标直接影响惯性测量装置的精度。惯性测量装置应用到陆用车载环境下进行寻北测量时,要求载体处于停止状态,并且有良好的测试环境,不能有扰动。在实际的工程应用中,车载环境通常会受到人员走动、发动机工作、风力等的影响,处于扰动环境,对陀螺信号造成极大的影响,如果不进行处理,将无法实现惯性测量,得不到可用的测量数据。这一问题反映出陀螺信号的处理,不能仅满足试验室环境,更应该针对实际工程应用采取有效的方法。
发明内容
本发明提供了一种适应车载扰动环境下寻北的挠性陀螺信号处理方法,在于对工程应用中的车载扰动环境下的陀螺信号,通过对信号中的坏值剔除、综合滤波等过程,解决由人员走动、发动机工作、低频风扰等带来的影响,其技术方案如下所述:
一种适应车载扰动环境下寻北的挠性陀螺信号处理方法,包括以下步骤:
S1:对挠性陀螺产生的数据按2000Hz解算周期,将0.5s周期内的数据记为一个数据,存入数组,数组共存入30个数据,也就是15s周期数据;
S2:对30个数据由小到大排序,排序后取中间的10个数据,剔除最小的10个数和最大的10个数;
S3:对中间的10个数据进行平均,进行数据归一化处理;将归一化处理后的数据代入惯性测量算法,解算测量值。
进一步的,步骤S1中,挠性陀螺产生的数据包括以下处理:
S11:DSP响应2000Hz的定时中断,对挠性陀螺的数据进行8次循环采样;
S12:对8次数据进行由小到大排序,排序后取中间的4个数据进行平均,剔除最小的两个数和最大的两个数;
S13:对进行平均得到的数据进行二阶低通滤波处理,滤波器截止频率166Hz,用于处理数据中扰动带来的噪声;
S14:对二阶低通滤波处理得到的数据进行陷波处理,陷波器截止频率346Hz,用于处理数据中与陀螺干扰噪声相关部分;
S15:对陷波处理得到的数据进行滞后补偿处理。
进一步的,步骤S13中,采用二阶低通数字离散化后公式:
式中:xn为步骤S2得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
进一步的,步骤S14中,采用限波器数字离散化后公式:
式中:xn为步骤S3得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
进一步的,步骤S15中,采用滞后数字离散化后公式:
f3(x)=401xn-400xn-1 (3)
式中:xn为步骤S4得到的数据,xn-1为上一个周期解算数据。
本方法通过惯性测量装置进行处理,在惯性测量装置内部安装有挠性陀螺、第一处理电路板、第二信号处理板和电源板,所述电源板为第一处理电路板和第二信号处理板供电,第一处理电路板为挠性陀螺供电并进行前端信号处理,第二信号处理板采集第一处理电路板上的陀螺信号并进行信号处理。
所述挠性陀螺的信号经过前端第一处理电路板的处理,完成带通、解调、滤波环节,供第二信号处理板的DSP+FPGA采集处理。
所述适应车载扰动环境下寻北的挠性陀螺信号处理方法,能够有效的适应车载环境下的扰动,提高产品性能,提高测量精度,满足实际工程应用。
附图说明
图1是本发明配属的惯性测量装置构成图;
图2是本发明的电路原理图;
图3是本发明的信号处理流程图;
图4是本发明实施案例的测量数据。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明。但是应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要的混淆本发明的概念。以下结合附图进一步说明本发明的实施方式。
如图1和图2所示,所述适应车载扰动环境下寻北的挠性陀螺信号处理方法,配属到惯性测量装置2,在惯性测量装置2内部安装有挠性陀螺1、第一处理电路板3、第二信号处理板4和电源板5,所述电源板5为第一处理电路板3和第二信号处理板4供电,第一处理电路板3为挠性陀螺1供电并进行前端信号处理,第二信号处理板4采集第一处理电路板3上的陀螺信号并进行信号处理。
所述挠性陀螺1的信号经过前端第一处理电路板3的处理,完成带通、解调、滤波环节,供第二信号处理板4的DSP+FPGA采集处理,处理的信号一路控制施矩脉冲,完成恒流源施矩电流控制;另一路进行信号处理,实现惯性测量功能,具体惯性测量的处理如下所述。
如图3所示,所述适应车载扰动环境下寻北的挠性陀螺信号处理方法,在发出信号采集指令后,包括以下步骤:
S1:DSP响应2000Hz(0.0005s周期)定时中断,对挠性陀螺1的数据进行8次循环采样;其中,2000Hz定时器是解算周期的意思,在硬件图框中看不到,在图3中体现;
S2:对8次数据进行由小到大排序,排序后取中间的4个数据进行平均,最小的两个数和最大的两个数作为坏值进行剔除,主要处理数据中的高频尖峰干扰;之所以取8个数据排序,是因为在工程应用时有时间限制,取其他时间达不到预期的试验效果。
S3:对数据进行二阶低通滤波处理,滤波器截止频率166Hz,主要处理数据中扰动带来的噪声;
采用二阶低通数字离散化后公式:
式中:xn为步骤S2得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
S4:对数据进行陷波处理,陷波器截止频率346Hz,主要处理数据中与陀螺干扰噪声相关部分;
采用限波器数字离散化后公式:
式中:xn为步骤S3得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
S5:对数据进行滞后补偿处理;
采用滞后数字离散化后公式:
f3(x)=401xn-400xn-1 (3)
式中:xn为步骤S4得到的数据,xn-1为上一个周期解算数据。
S6:对步骤S5产生的数据按2000Hz解算周期,重复采集1000次,对1000次数据求和,也就是0.5s周期内的数据记为一个数据,存入数组;数组共存入30个数据,也就是15s周期数据;
其中,本步骤中取1000次是工程测试最优的结果,其他次数也可以。取30个数据,15s周期同样因为在工程应用时有时间限制,同时也方便下一步的排序计算。
S7:对30个数据进行处理,先对数据进行由小到大排序,排序后取中间的10个数据进行平均,最小的10个数和最大的10个数作为坏值进行剔除,主要处理数据中风扰带来的0.1Hz~5Hz低频干扰;
S8:对处理后的10个数据平均,进行数据归一化处理;将归一化处理后的数据代入惯性测量算法,解算测量值。
本实施例采用图1提供的惯性测量装置在车辆上开展实际测试,测试结果的数据如图4所示。经过验证可以看出,采用本发明的陀螺信号处理方法能够有效的适应车载环境下的扰动,提高产品性能。
以上显示和描述了本发明专利的基本原理、主要特征和本发明专利的优点,本行业的技术人员应该了解,本发明专利不受上述实施例的限制,上述实施例和说明书中描述的只是说明本专利的原理,在不脱离本发明专利精神和范围的前提下,本发明专利还会有各种变化和改进,这些变化和改进都落入要求保护的本发明专利范围内。本发明专利要求保护范围由所附的权利要求书及其等效物界定。
Claims (6)
1.一种适应车载扰动环境下寻北的挠性陀螺信号处理方法,包括以下步骤:
S1:对挠性陀螺产生的数据按2000Hz解算周期,将0.5s周期内的数据记为一个数据,存入数组,数组共存入30个数据,也就是15s周期数据;挠性陀螺产生的数据包括以下处理:
S11:DSP响应2000Hz的定时中断,对挠性陀螺的数据进行8次循环采样;
S12:对8次数据进行由小到大排序,排序后取中间的4个数据进行平均,剔除最小的两个数和最大的两个数;
S13:对进行平均得到的数据进行二阶低通滤波处理,滤波器截止频率166Hz,用于处理数据中扰动带来的噪声;
S14:对二阶低通滤波处理得到的数据进行陷波处理,陷波器截止频率346Hz,用于处理数据中与陀螺干扰噪声相关部分;
S15:对陷波处理得到的数据进行滞后补偿处理;
S2:对30个数据由小到大排序,排序后取中间的10个数据,剔除最小的10个数和最大的10个数;
S3:对中间的10个数据进行平均,进行数据归一化处理;将归一化处理后的数据代入惯性测量算法,解算测量值。
2.根据权利要求1所述的适应车载扰动环境下寻北的挠性陀螺信号处理方法,其特征在于:步骤S13中,采用二阶低通数字离散化后公式:
式中:xn为步骤S2得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
3.根据权利要求1所述的适应车载扰动环境下寻北的挠性陀螺信号处理方法,其特征在于:步骤S14中,采用限波器数字离散化后公式:
式中:xn为步骤S3得到的数据,xn-1、yn-1为上一个周期解算数据,xn-2、yn-2为上两个周期解算数据。
4.根据权利要求1所述的适应车载扰动环境下寻北的挠性陀螺信号处理方法,其特征在于:步骤S15中,采用滞后数字离散化后公式:
f3(x)=401xn-400xn-1 (3)
式中:xn为步骤S4得到的数据,xn-1为上一个周期解算数据。
5.根据权利要求1所述的适应车载扰动环境下寻北的挠性陀螺信号处理方法,其特征在于:本方法通过惯性测量装置进行处理,在惯性测量装置内部安装有挠性陀螺、第一处理电路板、第二信号处理板和电源板,所述电源板为第一处理电路板和第二信号处理板供电,第一处理电路板为挠性陀螺供电并进行前端信号处理,第二信号处理板采集第一处理电路板上的陀螺信号并进行信号处理。
6.根据权利要求5所述的适应车载扰动环境下寻北的挠性陀螺信号处理方法,其特征在于:所述挠性陀螺的信号经过前端第一处理电路板的处理,完成带通、解调、滤波环节,供第二信号处理板的DSP+FPGA采集处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210274958.4A CN114563002B (zh) | 2022-03-21 | 2022-03-21 | 一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210274958.4A CN114563002B (zh) | 2022-03-21 | 2022-03-21 | 一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114563002A CN114563002A (zh) | 2022-05-31 |
CN114563002B true CN114563002B (zh) | 2024-04-09 |
Family
ID=81719315
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210274958.4A Active CN114563002B (zh) | 2022-03-21 | 2022-03-21 | 一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114563002B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0077491A1 (de) * | 1981-10-19 | 1983-04-27 | LITEF Litton Techn. Werke der HELLIGE GmbH | Verfahren zur raschen und direkten Azimut-Winkelbestimmung mit Hilfe eines Kreisels |
EP0250608A1 (de) * | 1986-06-18 | 1988-01-07 | LITEF GmbH | Verfahren und Einrichtung zur Azimut-Winkelbestimmung mit Hilfe eines Strapdown-Kreisels |
CN103674062A (zh) * | 2013-12-05 | 2014-03-26 | 广东电网公司电力科学研究院 | 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法 |
CN106595628A (zh) * | 2016-11-18 | 2017-04-26 | 天津津航技术物理研究所 | 数字光纤陀螺信号的野值实时剔除方法 |
CN109084761A (zh) * | 2018-08-07 | 2018-12-25 | 合肥正阳光电科技有限责任公司 | 一种基于惯性寻北装置的抗冲击寻北解算方法 |
CN109766787A (zh) * | 2018-12-21 | 2019-05-17 | 哈尔滨工程大学 | 一种单轴光纤陀螺抗扰动寻北计算方法 |
WO2020253854A1 (zh) * | 2019-06-21 | 2020-12-24 | 台州知通科技有限公司 | 移动机器人姿态角解算方法 |
KR102285366B1 (ko) * | 2021-06-08 | 2021-08-04 | 주식회사 신보 | 자이로 1n 노이즈 예측 필터링 장치 및 그 방법 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10982959B2 (en) * | 2016-09-06 | 2021-04-20 | The Charles Stark Draper Laboratory, Inc. | Fused sensor ensemble for navigation and calibration process therefor |
-
2022
- 2022-03-21 CN CN202210274958.4A patent/CN114563002B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0077491A1 (de) * | 1981-10-19 | 1983-04-27 | LITEF Litton Techn. Werke der HELLIGE GmbH | Verfahren zur raschen und direkten Azimut-Winkelbestimmung mit Hilfe eines Kreisels |
EP0250608A1 (de) * | 1986-06-18 | 1988-01-07 | LITEF GmbH | Verfahren und Einrichtung zur Azimut-Winkelbestimmung mit Hilfe eines Strapdown-Kreisels |
CN103674062A (zh) * | 2013-12-05 | 2014-03-26 | 广东电网公司电力科学研究院 | 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法 |
CN106595628A (zh) * | 2016-11-18 | 2017-04-26 | 天津津航技术物理研究所 | 数字光纤陀螺信号的野值实时剔除方法 |
CN109084761A (zh) * | 2018-08-07 | 2018-12-25 | 合肥正阳光电科技有限责任公司 | 一种基于惯性寻北装置的抗冲击寻北解算方法 |
CN109766787A (zh) * | 2018-12-21 | 2019-05-17 | 哈尔滨工程大学 | 一种单轴光纤陀螺抗扰动寻北计算方法 |
WO2020253854A1 (zh) * | 2019-06-21 | 2020-12-24 | 台州知通科技有限公司 | 移动机器人姿态角解算方法 |
KR102285366B1 (ko) * | 2021-06-08 | 2021-08-04 | 주식회사 신보 | 자이로 1n 노이즈 예측 필터링 장치 및 그 방법 |
Non-Patent Citations (4)
Title |
---|
一种实用的激光陀螺寻北仪数据处理技术;钱硕楠 等;微计算机信息;第23卷(第01期);全文 * |
叶西宁 等.智能仪表原理与设计技术.华东理工大学出版社,2021,(第版次:2021年2月第3版,印次:2021年2月第1次版),235-238. * |
李杏春 等.8098单片机原理及实用接口技术.北京航空航天大学出版社,1996,(第版次:1996年5月第1版,印次:1996年5月第1次印刷版),423. * |
王化祥.现代传感技术及应用.天津大学出版社,2016,(第版次:2016年6月第1版,印次:2016年6月第一次版),200-204. * |
Also Published As
Publication number | Publication date |
---|---|
CN114563002A (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107422379B (zh) | 基于局部自适应凸化方法的多尺度地震全波形反演方法 | |
CN108764073B (zh) | 一种结合频谱能量形态拟合的加速度滤噪和积分方法 | |
CN102353500B (zh) | 一种用于动平衡测量的不平衡信号提取方法 | |
CN109142863B (zh) | 一种电力系统测频方法及系统 | |
CN113820003B (zh) | 一种适用于桥梁振动监测的加速度实时重构动态位移方法 | |
CN112326245A (zh) | 一种基于变分希尔伯特黄变换的滚动轴承故障诊断方法 | |
CN104713575A (zh) | 一种闭环光纤陀螺频率特性的测试方法 | |
CN111624671A (zh) | 旋转加速度计重力梯度仪重力梯度解调相位角确定方法及装置 | |
CN114699042A (zh) | 一种频率检测方法、装置、电子设备及存储介质 | |
CN114563002B (zh) | 一种适应车载扰动环境下寻北的挠性陀螺信号处理方法 | |
Henry et al. | Prism signal processing for sensor condition monitoring | |
CN102589551B (zh) | 一种基于小波变换的船用光纤陀螺信号实时滤波方法 | |
CN110231117A (zh) | 一种基于s变换的拉索基频特征辨识方法 | |
CN103940447B (zh) | 一种基于自适应数字滤波器的系泊状态初始对准方法 | |
US8412478B2 (en) | Device for determining an error induced by a high-pass filter and associated error correction method | |
CN108181486A (zh) | 加速度信号的处理方法及装置 | |
CN105865318A (zh) | 应用于lvdt位移传感器的非相敏数字解调系统及方法 | |
CN107643080B (zh) | 一种惯性测量单元信号的采集处理方法 | |
CN110596425B (zh) | 一种无人机mems加速度传感器噪声消除方法 | |
Sabat et al. | Characterization of fiber optics gyro and noise compensation using discrete wavelet transform | |
CN104570118A (zh) | 一种基于双因素的自动识别与去除工业干扰的方法 | |
CN111693685B (zh) | 血栓弹力信号提取方法、存储介质及计算机设备 | |
Zhang et al. | Research on the signal separation method based on multi-sensor cross-correlation fusion algorithm | |
CN117590431B (zh) | 一种抑制扫频干扰方法及装置 | |
CN116340730B (zh) | 一种火箭弹弹性振动频率在线辨识与抑制方法 |
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 |