CN108733880A - 系统以及用于调整测量设备的方法、计算机程序和装置 - Google Patents

系统以及用于调整测量设备的方法、计算机程序和装置 Download PDF

Info

Publication number
CN108733880A
CN108733880A CN201810325458.2A CN201810325458A CN108733880A CN 108733880 A CN108733880 A CN 108733880A CN 201810325458 A CN201810325458 A CN 201810325458A CN 108733880 A CN108733880 A CN 108733880A
Authority
CN
China
Prior art keywords
operating point
adjusted
output parameter
affiliated
point
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
CN201810325458.2A
Other languages
English (en)
Other versions
CN108733880B (zh
Inventor
B.比肖夫
E.克洛彭布格
M.迈斯特
N-T.世维
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.)
Robert Bosch GmbH
Original Assignee
Robert Bosch GmbH
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 Robert Bosch GmbH filed Critical Robert Bosch GmbH
Publication of CN108733880A publication Critical patent/CN108733880A/zh
Application granted granted Critical
Publication of CN108733880B publication Critical patent/CN108733880B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M15/00Testing of engines
    • G01M15/04Testing internal-combustion engines

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Combustion & Propulsion (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Chemical & Material Sciences (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Complex Calculations (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

本发明涉及系统以及用于调整测量设备的方法、计算机程序和装置。一种用于调整测量设备(100)的方法,其中产生待试样品的下一个要由所述测量设备(100)调整的工作点(),其中从可能的所要调整的工作点()的集合(C)中选择所要调整的工作点(),其中所述选择根据已经调整出的工作点()和所属的被确定的输出参量()来进行。

Description

系统以及用于调整测量设备的方法、计算机程序和装置
技术领域
本发明涉及一种用于调整测量设备的方法、一种计算机程序和一种装置,所述计算机程序和装置被设立为实施该方法,而且本发明涉及一种包括该测量设备和该装置的系统。
背景技术
从DE 10 2006 201 559 A1中公知一种用于利用多个测量点来测量机动车的内燃机以便分别得到至少一个输出参量的值的方法。在这种情况下规定:从测量点的集合中选择一个测量点,其方式是将这些测量点根据技术系统的工作点来分组,而且针对这些工作点中的每个工作点,根据距被分配给那个工作点的初始测量点的升高的间距来从测量点的集合中依次选择这些测量点。
发明内容
本发明的优点
与此相应地,具有独立权利要求1的特征的方法具有如下优点:该方法提供了使用已经测量过的数据点从而找到输入空间的如下区域的可能性,对于所述区域来说,数据点还没有反映出足够的信息。借此,运行测量设备是特别高效的。
有利的扩展方案是从属权利要求的主题。
在第一方面,本发明涉及一种用于调整、尤其是运行测量设备的方法,其中产生待试样品的下一个要由测量设备调整的工作点,其中从可能的要调整的工作点的集合中选择所要调整的工作点,其中该选择根据已经调整出的工作点和所属的被确定的输出参量来进行。接着,尤其可以规定:测量设备调整待试样品的所选择的要调整的工作点并且确定所属的输出参量。所要调整的工作点与所属的输出参量构成的对各被称作一个数据点。已经调整出的工作点和所属的被确定的输出参量的集合也被称作训练数据。
每个新测量的点都可以立即被添加到训练数据,以便制成新的模型,用于描述在调整出的工作点与所属的被确定的输出参量之间的函数关系。现在包含在整个的、被更新的数据组中的信息被用于从已经调整出的工作点推断出新的工作点。基于该模型可以确定如下工作点,所述工作点根据可预先给定的标准来提供用于总模型的最高的信息增益,所述总模型描述了工作点与所属的输出参量之间的函数关系。在此,可预先给定的标准可以以基于数据的模型为基础。接着,该工作点作为下一个要调整的工作点来运转。以这种方式可以在线地制成试验计划,所述试验计划提供特别有参考性的训练数据,以便制成尽可能好的模型。
可以规定:所述可预先给定的标准将点压缩在用于测量的输入空间(即可能的所要调整的工作状态的空间)的输出参量表现得高度动态的那个区域内。以这种方式,测量过的数据点可以反映出输出参量的(也许出现的)局部变化的动态性。
因而,在另一方面可以规定,根据如下情况来选择所要调整的工作点:所确定的输出参量与已经调整出的工作点的函数关系通过根据可能的所要调整的工作点中的相应的工作点的优选地线性的函数被描述得多好。也就是说根据所要调整的工作点中的相应的工作点,确定用于描述所确定的输出参量与调整出的工作点之间的函数关系的线性函数,而且针对该线性函数来确定在实际上被确定的输出参量与借助于该函数针对该工作点所期望的输出参量之间的偏差。换句话说,因此可以规定:测量压缩在输入空间的输出参量表现得动态的那个区域内的工作点,其方式是将数据与模型、尤其是局部的线性模型进行比较。接着,将新的测量点刚好设置在数据和局部线性的模型没有良好地相配的那个地方。
在其它方面可以规定:通过根据可能的所要调整的工作点中的相应的工作点的线性函数来描述所确定的输出参量与已经调整出的工作点的函数关系的质量由不合理值来表征。
因而,在一个扩展方案中可以规定:从可能的所要调整的工作点的集合中选择如下那个可能的所要调整的工作点,对于该工作点来说,所确定的输出参量与已经调整出的工作点的函数关系通过根据所述可能的所要调整的工作点的优选地线性的函数被描述得最差。
还可以规定:优选地线性的函数与可能的所要调整的工作点中的相应的工作点的关系通过以所述相应的可能的所要调整的工作点为中心的核密度函数来表征。所述核密度函数尤其可以是高斯核密度函数。但是所述核密度函数也可以通过其它函数来给定,例如柯西核(Cauchy-Kern)、皮卡核(Picard-Kern)或者叶帕涅奇尼科夫核(Epanechnikov-Kern)。
为此,尤其是可以根据不合理值来确定核密度函数的核宽度(英文kernelwidth)。不合理值表征:输出参量通过如下那个线性函数被表征得多好,所述线性函数通过以工作点之一为中心的核密度函数来表征。如果模型不合理,那么不合理值升高,而且核宽度变得更小,使得更少的数据点影响局部模型。
所述优选地线性的函数尤其可以利用局部线性核回归来确定。
在另一方面,本发明也涉及一种计算机程序,所述计算机程序被设立为在所述计算机程序在计算机上运行时执行该方法。
附图说明
随后,本发明的实施方式参考随附的附图进一步予以阐述。在所述附图中:
图1示意性地示出了其中应用该方法的测量试验台;
图2示出了关于本发明的第一方面的流程图;
图3示出了关于本发明的第二方面的流程图。
具体实施方式
图1示例性地示出了试验台100,在所述试验台中可以应用按照本发明的方法。试验台100由操控装置200来操控,所述操控装置可以在构造上集成到试验台100中,但是也可以在构造上与该试验台分开地来布置。操控装置200将操控指令p传输给试验台100,所述试验台100借此被操控为使得所述试验台100在要由其测量的待试样品(未示出)中调整出可预先给定的工作点x。在这种情况下,工作点x表征待试样品的工作状态。待试样品可以是硬件,或者可以是软件,或者可以是混合物体、由硬件和软件构成的混合物体。
试验台100例如可以是发动机试验台,其中待试样品是内燃机。
在待试样品上,试验台100借助于例如传感器(同样未示出)确定输出参量值y,所述输出参量值表征待试样品的输出参量。试验台100将所确定的输出参量值y传输给操控装置200。如果在下文谈到多个工作点x或所属的输出参量值y中的一个,则这些参量x、y都配备下标i。测量参量x是d维参量,其中d是自然数,输出参量y例如通过一维数值来给定。
操控装置200将所确定的输出参量值y传输给工作点发生器300,所述工作点发生器根据一个或多个所确定的输出参量值y来确定如下那个工作点x,接着应该在所述工作点x测量待试样品的输出参量y。工作点发生器300可以与操控装置200和/或试验台100在构造上分开,而且例如可以在远离操控装置200和试验台100的服务器上实现。
图2和3图解说明了按照优选的实施方式的用于生成工作点x的方法的实施方式。该方法可以在工作点发生器300中和/或在操控装置200中运行。
在一个优选的实施方式中,操控装置200和工作点发生器300包括机器可读的存储介质201、301,在所述机器可读的存储介质上分布式地存储有计算机程序,所述计算机程序被设立为:在所述计算机程序在操控装置200和/或工作点发生器300上实施时,实施按照本发明的方法。
图2示出了按照本发明的第一方面的流程图。在步骤1000中,接收初始数据组,其中维度n等于初始维度n0。测量点 例如可以是离线地生成的,所属的被确定的输出参量例如是利用试验台100来确定的。在确定该初始数据组Dn之后,不仅是步骤1010而且可选地是步骤2001,在所述步骤2001的情况下将该初始数据组Dn提供给在图3中图解说明的可选的方法,用来估计方差σ。
在步骤1010中,接收候选点。借助于在图2中图解说明的方法,从这些候选点C中选择下一测量点。这些候选点C例如可以在这里所描述的方法的准备阶段借助于DoE(Design-of-Experiment(实验设计))方法来确定。
在步骤1010中,还按照如下公式来确定核宽度的初始值,其中每个j=1, …,m0:
在这种情况下,参数是可预先给定的常数,该常数优选被选择>3。方程(1)例如可以用数值方法来求解。
还将迭代计数器t初始化到值t=1。
接着是步骤1020,在所述步骤1020中,针对候选点中的每个候选点形成一个局部线性核回归模型(在下文简称:LLKR模型)。为此,分别对优化问题
优选地利用如下高斯核函数来进行求解
现在,利用公式(2)的解来确定如下预测值
接着是步骤1030,在所述步骤1030中,确定权重
和不合理值(Unplausibilität)
而且从中针对候选点中的每个候选点确定成本项
在这种情况下,针对每个数据组Dn,可以通过步骤1001从在图3中图解说明的方法中获得方差估计值,或者可以针对初始数据组Dn0从在图3中图解说明的方法中一次性地获得方差估计值,或者将方差估计值预先给定为固定值。
接着是步骤1040,在所述步骤1040中,按照公式
来确定优选的值
所属的候选点被选为下一测量数据点,而且作为第(n+1)个测量数据点被添加给测量数据点的集合。以这种方式形成集合。同时,将所选择的候选点从候选点C的列表中除去。
接着是步骤1060,在所述步骤1060中,接收属于测量数据点的输出参量。该输出参量有利地是由试验台100确定的,其方式是该试验台已经将待试样品置于工作点。该输出参量被添加给所确定的输出参量的集合并且这样形成
测量数据点的所述集合以及所确定的输出参量的集合通过步骤2002可选地被提供给在图3中图解说明的方法。
接着是步骤1070,在所述步骤1070中,首先借助于公式(4)来确定针对新添加的测量数据点的估计值。紧接着,针对每个和每个,借助于公式(6)重新计算不合理值
接着是步骤1080,在所述步骤1080中,利用如下公式来确定核宽度
其中每个。在这种情况下,参数可以与在公式(1)中的选择得相同,但是该参数也可以取其它值。对公式(10)的求解可以利用与对公式(9)的求解相同的方法来实现。
在步骤1090中检查是否满足收敛标准。这例如可以通过迭代计数器t超过可预先给定的迭代最高值tmax来给定。可替换地,收敛标准也可以通过询问在该实施例中通过公式(7)限定的成本项是否收敛来给定。如果满足收敛标准(输出“j”),那么该方法在步骤1100中结束。否则,变量n和t的值加1,而变量m的值减1,并且该方法重新开始步骤1020。
图3示出了本发明的优选的实施方式的另一方面。该方法的所述可选的部分对被用于计算不合理值的所估计的方差进行估计。
该方法开始于步骤2001,在所述步骤2001中,从在图2中图解说明的方法部分接收初始数据组Dn。同样接收n的初始值。迭代计数器t2被设置到值1。
在紧接着的步骤2000中,确定核宽度,所述核宽度在该方法部分中用来表示。针对每个,这利用公式
来实现。在这种情况下,参数可以如在公式(1)中那样取可预先给定的值,优选地>3。对公式(10)的求解可以利用与对公式(1)或(9)的求解相同的方法来实现。
接着是步骤2010,类似于步骤1020地,在所述步骤2010中形成LLKR模型。不过在该步骤中,不同于步骤1020地,不是围绕着候选点来形成LLKR模型,而是围绕着针对的测量数据点来形成LLKR模型。
为此,针对每个,首先形成核函数,所述核函数优选地具有如下高斯形式:
紧接着,将参数确定为
现在,利用公式(12)的解来确定预测值
在接下来的步骤2020中,确定所估计的方差。首先针对每个确定参量
而且从中确定
该值通过步骤1001被反馈给该方法的在图2中图解说明的方面。在接下来的步骤2030中,通过步骤2002来接收工作参量和输出参量的新的值。针对并且针对来确定参量,其方式是首先借助于公式(11)和(12)形成围绕着的LLKR模型,其中作为替代,取针对在步骤1040中借助于公式(4)确定的下标的核宽度的在步骤1010或1080中利用公式(1)或(9)形成的值的值,作为(还没有被确定的)核宽度的值。
现在,针对所有,以
来计算不合理值
在接下来的步骤2050中,通过对公式
的求解来确定核宽度。在这种情况下,参数可以被选择得与在公式(10)中的相同,但是该参数也可以取其它值。对公式(17)的求解可以利用与对公式(1)、(9)或(10)的求解相同的方法来实现。
现在,n的值以及迭代计数器t2的值加1,而该方法分岔回到步骤2010,在所述步骤2010中开始下一次迭代。
在一个优选的实施方式中,在操控装置200中在操控试验台100的情况下运行步骤1060,并且在工作点发生器300中运行所有其它步骤。其它划分当然是可设想的。
对于本领域技术人员来说易于理解的是,本发明可以以软件来实现,或者是以硬件来实现的,或者是以硬件和软件的混合形式来实现的。

Claims (15)

1.一种用于调整测量设备(100)的方法,其中产生待试样品的下一个要由所述测量设备(100)调整的工作点(),其中从可能的所要调整的工作点()的集合(C)中选择所要调整的工作点(),其中所述选择根据已经调整出的工作点()和所属的被确定的输出参量()来进行。
2.根据权利要求1所述的方法,其中在由所述测量设备(100)调整出所述待试样品的下一个要调整的工作点()以及检测所属的输出参量()之后,所述下一个要调整的工作点()被添加给已经调整出的工作点(),以便形成已经调整出的工作点()的扩展的集合,而且所检测到的所属的输出参量()被添加给所属的被确定的输出参量(),以便形成所属的被确定的输出参量()的扩展的集合,而且其中所述下一个要调整的工作点()从可能的所要调整的工作点()的集合(C)中被除去,以便形成可能的所要调整的工作点的减小的集合,其中根据已经调整出的工作点()的扩展的集合和所属的被确定的输出参量()的扩展的集合从可能的所要调整的工作点的减小的集合中选择所述待试样品的再下一个要调整的工作点。
3.根据权利要求1或2所述的方法,其中根据所确定的输出参量()与已经调整出的工作点()的函数关系通过根据所述可能的所要调整的工作点()中的相应的工作点的函数被描述得多好,选择所要调整的工作点()。
4.根据权利要求3所述的方法,其中通过根据所述可能的所要调整的工作点()中的相应的工作点的线性函数来描述所述所确定的输出参量()与所述已经调整出的工作点()的函数关系的质量由不合理值()来表征。
5.根据权利要求3或4所述的方法,其中从可能的所要调整的工作点()的集合(C)中选择如下那个可能的所要调整的工作点(),对于所述工作点来说,所述所确定的输出参量()与所述已经调整出的工作点()的函数关系通过根据所述可能的所要调整的工作点()的线性函数被描述得最差。
6.根据权利要求3至5之一所述的方法,其中线性函数与所述可能的所要调整的工作点()中的相应的工作点的关系通过以所述相应的可能的所要调整的工作点()为中心的核密度函数()来表征。
7.根据权利要求5和6所述的方法,其中根据所述不合理值()来确定核密度函数()的核宽度()。
8.根据权利要求5至7之一所述的方法,其中所述线性函数利用局部线性核回归来确定。
9.根据权利要求3或4所述的方法,其中线性函数借助于高斯过程来确定。
10.根据权利要求3至9之一和权利要求4所述的方法,其中质量()表征所述所确定的输出参量()与所述已经调整出的工作点()的函数关系通过根据所述可能的所要调整的工作点()中的相应的工作点的线性函数被描述得多好。
11.根据上述权利要求之一和权利要求4所述的方法,其中根据噪声值()来选择所述不合理值(),所述噪声值表征在调整出的工作点()和所属的被确定的输出参量()中的噪声。
12.一种装置,其被设立为实施根据上述权利要求之一所述的方法。
13.一种系统,其包括测量设备(100)和根据权利要求12所述的装置。
14.一种计算机程序,其被设立为实施根据权利要求1至11之一所述的方法。
15.一种机器可读的存储介质,在其上存储有根据权利要求14所述的计算机程序。
CN201810325458.2A 2017-04-13 2018-04-12 系统以及用于调整测量设备的方法、计算机程序和装置 Active CN108733880B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102017206395.5A DE102017206395A1 (de) 2017-04-13 2017-04-13 Verfahren, Computerprogramm und Vorrichtung zum Einstellen einer Messvorrichtung und System
DE102017206395.5 2017-04-13

Publications (2)

Publication Number Publication Date
CN108733880A true CN108733880A (zh) 2018-11-02
CN108733880B CN108733880B (zh) 2024-04-23

Family

ID=63679156

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810325458.2A Active CN108733880B (zh) 2017-04-13 2018-04-12 系统以及用于调整测量设备的方法、计算机程序和装置

Country Status (2)

Country Link
CN (1) CN108733880B (zh)
DE (1) DE102017206395A1 (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050263727A1 (en) * 2004-05-31 2005-12-01 Mitutoyo Corporation Surface scan measuring device, surface scan measuring method, surface scan measuring program and recording medium
WO2015110163A1 (de) * 2014-01-23 2015-07-30 Carl Zeiss Industrielle Messtechnik Gmbh Verfahren zur auswertung von koordinatenmessdaten und vorrichtung
JP2016024058A (ja) * 2014-07-22 2016-02-08 キヤノン株式会社 波面計測方法および装置、面形状計測方法および装置、並びに加工装置
DE102015200295A1 (de) * 2015-01-13 2016-07-14 Robert Bosch Gmbh Verfahren und Vorrichtung zum Ermitteln von Trainingsdatenpunkten zum Erstellen eines Modells für ein physikalisches System
DE102015203771A1 (de) * 2015-03-03 2016-09-08 Robert Bosch Gmbh Verfahren und Vorrichtung zum Auswerten eines datenbasierten Funktionsmodells
DE102015216953A1 (de) * 2015-09-04 2017-03-09 Robert Bosch Gmbh Verfahren und Vorrichtung zum Vermessen einer zu testenden Einheit

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102016201559A1 (de) 2016-02-02 2017-08-03 Robert Bosch Gmbh Verfahren und Vorrichtung zum Vermessen eines zu testenden Systems

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050263727A1 (en) * 2004-05-31 2005-12-01 Mitutoyo Corporation Surface scan measuring device, surface scan measuring method, surface scan measuring program and recording medium
WO2015110163A1 (de) * 2014-01-23 2015-07-30 Carl Zeiss Industrielle Messtechnik Gmbh Verfahren zur auswertung von koordinatenmessdaten und vorrichtung
EP3077763A1 (de) * 2014-01-23 2016-10-12 Carl Zeiss Industrielle Messtechnik GmbH Verfahren zur auswertung von koordinatenmessdaten und vorrichtung
JP2016024058A (ja) * 2014-07-22 2016-02-08 キヤノン株式会社 波面計測方法および装置、面形状計測方法および装置、並びに加工装置
DE102015200295A1 (de) * 2015-01-13 2016-07-14 Robert Bosch Gmbh Verfahren und Vorrichtung zum Ermitteln von Trainingsdatenpunkten zum Erstellen eines Modells für ein physikalisches System
DE102015203771A1 (de) * 2015-03-03 2016-09-08 Robert Bosch Gmbh Verfahren und Vorrichtung zum Auswerten eines datenbasierten Funktionsmodells
DE102015216953A1 (de) * 2015-09-04 2017-03-09 Robert Bosch Gmbh Verfahren und Vorrichtung zum Vermessen einer zu testenden Einheit

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
OTTO DYKSTRA,JR.: "The Augmentation of Experimental Data to Maximize", TECHNOMETMC, pages 682 - 685 *

Also Published As

Publication number Publication date
DE102017206395A1 (de) 2018-10-18
CN108733880B (zh) 2024-04-23

Similar Documents

Publication Publication Date Title
JP4898805B2 (ja) システムの構成パラメータ・セットを最適化するための方法、プログラムおよび装置
CN103577658B (zh) 智能终端硬件的评估方法及评估系统
JP6404909B2 (ja) 技術的な系の出力量のモデルを算出する方法
CN108304966A (zh) 一种短期月售电量预测方法
JP2005172830A5 (zh)
CN110262959B (zh) 底层服务压力测试方法、装置、电子设备及存储介质
CN110287332A (zh) 云环境下仿真模型选择方法与装置
CN108133223A (zh) 确定卷积神经网络cnn模型的装置和方法
CN108733880A (zh) 系统以及用于调整测量设备的方法、计算机程序和装置
CN110377941A (zh) 建立卫星温度场的罚盲似然克里金代理模型的方法
US20100223035A1 (en) Test Time Calculator
CN104079566B (zh) 终端、服务器和页面中结果属性值的计算方法
CN112214827A (zh) 基于多应力的轨道交通电子控制装置寿命评估方法及装置
Sim et al. Performance of new signal‐to‐noise ratio estimation for SEM images based on single image noise cross‐correlation
CN113704942B (zh) 供水管网未知节点压力估计方法、装置及计算机系统
CN106888133B (zh) 应用测速方法及装置
CN108804670A (zh) 数据推荐方法、装置、计算机设备和存储介质
Shen et al. Evaluation of multi-order derivatives by local radial basis function differential quadrature method
CN113032229A (zh) 一种java性能测试方法、系统及介质
CN107292045B (zh) 一种含阻尼结构有限元模型的复模态评估方法
JP2010510509A (ja) ディメンション計測システム用の信号処理方法及び装置
US20230333080A1 (en) Method and system of predicting fabric feature
CN110248372A (zh) 一种模拟器检测的方法、装置、存储介质及计算机设备
CN110826816B (zh) 船舶造修的能耗确定方法、装置、设备和存储介质
CN111595433A (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