CN107729694A - 一种基于神经网络的多参数电磁场建模仿真方法 - Google Patents
一种基于神经网络的多参数电磁场建模仿真方法 Download PDFInfo
- Publication number
- CN107729694A CN107729694A CN201711145631.2A CN201711145631A CN107729694A CN 107729694 A CN107729694 A CN 107729694A CN 201711145631 A CN201711145631 A CN 201711145631A CN 107729694 A CN107729694 A CN 107729694A
- Authority
- CN
- China
- Prior art keywords
- neutral net
- training
- transmission function
- branch
- parameter
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Data Mining & Analysis (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Computational Linguistics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
本发明提供了一种基于神经网络的多参数电磁场建模仿真方法,属于电磁场建模仿真技术领域。相比于现有技术的方法只能输出一个参数,本发明所述方法共有三个相互平行独立地分支用于输出三种不同的参数。同时,数据分类技术和支持向量机同样也被包含在每个分支中,用于将输入进行分类。因此当几何变量被输入后,本方法可以快速从各分支分别输出S参数、增益和方向图。
Description
技术领域
本发明属于电磁场建模仿真技术领域,具体涉及一种基于神经网络的多参数电磁场建模仿真方法。
背景技术
近些年来,基于神经网络的计算机辅助设计方法被认为是一种高效的电磁场仿真建模技术。一旦完成训练过程,训练成熟的神经网络能替代耗时的电磁场仿真而明显的加快电磁场仿真建模和设计。将神经网络和传递函数构建为神经网络-传递函数模型被认为是一种先进的设计方法。
文献“A new training approach for parametric modeling of microwavepassive components using combined neural networks and transfer functions”公开了一种解决传递函数不连续性问题的神经网络-传递函数的方法;文献“High-dimensionalneural-network technique and applications to microwave filter modeling”公开了一种分解高维器件并对每一个低维部分进行建模仿真的神经网络-传递函数的方法;文献“Parametric modeling of EM behavior of microwave components using combinedneural networks and pole-residue-based transfer functions”公开了一种基于新型传递函数的神经网络-传递函数的方法。
但目前所提出的传递函数-神经网络模型只有一个输出,即只能输出器件的S参数,对于建模仿真复杂天线,需同时考虑工作频率,带宽,增益和方向图等等,仅仅一个S参数是不够的。对于传递函数-神经网络模型内部的阶数改变问题,通常采用使阶数最大化或者动态调整阶数的方法;但是阶数最大化得到的神经网络仿真结果的精度往往不是很高,而动态调整阶数的方法在计算中又太过复杂,难以实现。
发明内容
本发明的目的是提供一种基于神经网络的多参数电磁场建模仿真方法,以解决现有的神经网络-传递函数的方法不能为复杂天线提供多参数建模仿真的问题,并且利用支持向量机从另一个方向解决传递函数-神经网络模型内部的阶数改变问题。
本发明所提出的技术问题是这样解决的:
一种基于神经网络的多参数电磁场建模仿真方法,该方法结构上包含三个分支,其中第一分支用于仿真S参数,第二分支用于仿真增益,第三分支用于仿真方向图,具体包括以下步骤:
步骤1:神经网络-传递函数建模仿真方法的训练:
本方法采用的传递函数为:
其中,pi和ri分别为传递函数的极点和留数,i=1,…,Q,Q为传递函数的阶数,s在第一分支和第二分支为频率,在第三分支为角度。
每个分支的训练具体过程如下:
步骤1-1.将几何变量输入到电磁仿真软件,通过全波仿真得到对应于几何变量的天线的S参数、增益或方向图O={O1,…,OW},其中,W是频率采样点的个数。
步骤1-2.根据步骤1-1中得到的仿真结果,通过使用矢量拟合法,得到神经网络训练数集中每组数据所对应传递函数的极点和留数。由于传递函数的阶数不同会影响神经网络的训练精度,因此为了准确地映射几何变量和训练数集,根据传递函数的不同的阶数,每个分支都包含有若干个神经网络。
同样的,根据传递函数的阶数,原始的训练数集被分为若干类Ck(k=1,2,…,K),其中K是各个分支的类别的个数;对于具有相同的传递函数阶数的训练样本被分为一类,每一类的阶数可以用Qk(k=1,2,…,K)表示;每一类只用于训练一个神经网络。
步骤1-3.通过神经网络的训练,得到传递函数的极点、留数和几何变量间的映射关系,进而得到传递函数的具体形式;带入自变量s,得到神经网络的输出O'={O'1,…,O'W}。
步骤1-4.通过调整神经网络内部的权重和阈值,对不同的几何变量,使O和O'之间的误差最小化。
步骤1-5.将几何变量输入到支持向量机,支持向量机的输出为每个训练样本的阶数Q'={Q'1,…,Q'K1},通过调整支持向量机内部的权重和阈值使Q'和传递函数实际的阶数{Q1,…,QK}之间的误差最小。
步骤2.神经网络-传递函数建模仿真方法的测试:
在测试过程中,将几何变量输入至训练好的支持向量机中,从而得到所对应的传递函数的阶数;然后,将几何变量输入至对应的神经网络得到仿真结果。
本发明的有益效果是:输入几何参数后,本发明可以同时输出天线的S参数,增益及方向图,从而解决了以往神经网络-传递函数模型无法仿真结构复杂天线的问题,另外由于支持向量机的引入,将阶数改变的问题以一种简单快速的方式得到了解决。
附图说明
图1为本发明所述方法训练步骤的流程图;
图2为本发明所述方法测试步骤的流程图;
图3为任一分支的训练步骤流程图;
图4为支持向量机的训练步骤流程图;
图5为电磁带隙结构图,其中(a)俯视图,(b)侧视图,(c)仰视图,(d)单元;
图6为馈电天线的结构,其中(a)俯视图,(b)侧视图;
图7为法布里帕罗天线的结构图;
图8为支持向量机S参数的分类结果;
图9为数据在训练数据范围内本发明所提出的模型和HFSS的比较,其中(a)样本1和样本2的S11参数,(b)样本1和样本2的增益,(c)样本1和样本2在10GHz的方向图;
图10为数据在训练数据范围外本发明所提出的模型和HFSS的比较,其中(a)样本1和样本2的S11参数,(b)样本1和样本2的增益,(c)样本1和样本2在10GHz的方向图;
图11为天线1和天线2的优化结果,其中(a)S11参数,(b)增益,(c)在10GHz处的方向图。
具体实施方式
下面结合附图和实施例对本发明进行进一步说明。
一种基于神经网络的多参数电磁场建模仿真方法,该方法结构上包含三个分支,其中第一分支用于仿真S参数,第二分支用于仿真增益,第三分支用于仿真方向图,具体包括以下步骤:
步骤1:神经网络-传递函数建模仿真方法的训练,其流程图如图1所示:
本方法采用的传递函数为:
其中,pi和ri分别为传递函数的极点和留数,i=1,…,Q,Q为传递函数的阶数,s在第一分支和第二分支为频率,在第三分支为角度。
每个分支的训练具体过程如下:
步骤1-1.将几何变量输入到电磁仿真软件,通过全波仿真得到对应于几何变量的天线的S参数、增益或方向图O={O1,…,OW},其中,W是频率采样点的个数。
步骤1-2.根据步骤1-1中得到的仿真结果,通过使用矢量拟合法,得到神经网络训练数集中每组数据所对应传递函数的极点和留数。由于传递函数的阶数不同会影响神经网络的训练精度,因此为了准确地映射几何变量和训练数集,根据传递函数的不同的阶数,每个分支都包含有若干个神经网络。
同样的,根据传递函数的阶数,原始的训练数集被分为若干类Ck(k=1,2,…,K),其中K是各个分支的类别的个数;对于具有相同的传递函数阶数的训练样本被分为一类,每一类的阶数可以用Qk(k=1,2,…,K)表示;每一类只用于训练一个神经网络。
步骤1-3.通过神经网络的训练,得到传递函数的极点、留数和几何变量间的映射关系,进而得到传递函数的具体形式;带入自变量s,得到神经网络的输出O'={O'1,…,O'W},本步骤流程图如图3所示。
步骤1-4.通过调整神经网络内部的权重和阈值,对不同的几何变量,使O和O'之间的误差最小化。
步骤1-5.将几何变量输入到支持向量机,支持向量机的输出为每个训练样本的阶数Q'={Q'1,…,Q'K1},通过调整支持向量机内部的权重和阈值使Q'和传递函数实际的阶数{Q1,…,QK}之间的误差最小,本步骤流程图如图4所示。
步骤2.神经网络-传递函数建模仿真方法的测试,其流程图如图2所示:
在测试过程中,将几何变量输入至训练好的支持向量机中,从而得到所对应的传递函数的阶数;然后,将几何变量输入至对应的神经网络得到仿真结果。
本实施例采用如图5所示的法布里帕罗天线作为算例进行有效性的验证。金属贴片和正方形孔径分别被焊接在厚度为T电介质常数为εr的介质基板上,其中贴片单元lb的尺寸为8mm,lp和la分别为正方形孔径的内外尺寸。
本实施例所采用的法布里帕罗天线的馈电结构如图6所示,寄生贴片被设计在RogersRT/duroid 5880介质上(εr=2.2andtanδ=0.0009)。馈电天线被放置接地板的上方和馈电线通过接地板上的开槽结合在一起。该馈电天线的结构参数如下:wp=9.3mm,w1=1.2mm,w2=2.3mm,ws=2.3mm,L1=9.5mm,Ls=8.2mm,Lstub=3mm,hair=2.5mm,and h1=h2=0.787mm。
法布里帕罗的整体结构如图7所示,馈电天线被放置在腔体中间,hc为电磁带隙层和地板层间的间隙。三个在法布里帕罗天线扮演重要角色的几何变量x=[lp la hc]T被设为输入。
如上表所示,几何变量有两种形式的训练数据:形式1的数据相对较窄,形式2的数据相对较宽。在两种形式中,训练和测试数据分别根据八层(64个训练数据)和六层(36个测试数据)的实验设计方法确定。通过电磁场仿真得到训练和测试数据的时间分别为21.33小时和12小时。同时,本实施例用Hecht–Nelson方法确定了神经网络中的隐含层的节点数,即输入层的节点数为n,隐含层的节点数为2n+1。HFSS 15.0软件被用于电磁场全波仿真,并产生建模所用的训练和测试数据。本例所有的计算都在包含16GB RAM的Intel i7-48702.50GHz平台下完成。
(1)S参数
在第一分支中,对于形式1的数据S参数传递函数阶数范围为8到10,形式2的传递函数的阶数为8到12。如下表所示,训练样本根据传递函数的阶数被分为若干类用于神经网络的训练。
同时,几何变量及其对应的传递函数的阶数分别被设为支持向量机的输入和输出。对于36个测试样本,分类的结果如图8所示。训练后的支持向量机的分类精度为97.22%。建模完成后,第一分支的形式1和形式2的训练误差分别为0.399%和0.424%,同时平均测试误差分别为0.615%和0.672%。
(2)增益
在第二分支中,增益所对应的传递函数的参数被设定为神经网络的输出。建模之后,形式1和形式2增益的平均训练误差分别为0.873%和0.857%,同时测试误差分别为0.954%和0.971%。
(3)方向图
在第三分支中,神经网络被用于映射集合变脸和在10GHz处提取的方向图的传递函数的参数。和以上两个分支不同的是,第三分支将角度作为额外的输入变量,范围为[-2π,2π]。经过建模过程,对于形式1和形式2的数据平均训练误差分别为2.645%和2.448%,平均测试误差分别为2.964%和2.912%。由于不平滑的曲线,所以分支3的训练误差比起第一分支和第二分支要大一些。本发明所提出的神经网络-传递函数模型的建模结果下表所示:
(4)综合测试
图9展示了对两个不同的法布里帕罗天线几何测试样本,本发明所提出的神经网络-传递函数模型和HFSS仿真结果。两个几何变量的样本x1=[5.62 5.44 14.63]T和x2=[5.87 5.51 15.36]T都在训练数据分范围内。可以看出本发明所提出的仿真模型可以实现较高的仿真的精度,即便这些数据未出现在训练样本中。
另外两个在训练数据之外的几何样本x′1=[5.3 4.9 14.2]T和x'2=[6.5 4.916.1]T被选为测试本仿真模型。从图10看出,尽管这些数据在训练数据之外,但是本仿真模型依然达到了较高的仿真精度。
(5)优化应用
当模型训练成熟后,它可替代费时的电磁场仿真应用于优化中。作为优化的实例,两个单独的法布里帕罗天线被优化以达到两个不同的设计规格。
天线1的优化目标:
1)在8.75-11.25GHz范围内,|S11|≤-10dB;
2)相对3dB增益带宽达到32%;
3)主瓣增益Gmax≥12.5dB。
天线2的优化目标:
1)在10-11GHz范围内,|S11|≤-10dB;
2)相对3dB增益带宽达到21%;
3)主瓣增益Gmax≥14dB。
通过NSGA-II优化算法不断调用训练好的仿真模型即可得到最后的优化结果。变量的初始值设为xinitial=[15 5.5 6]T。优化过程仅用了60秒,即可实现每个天线最后的优化结果,如图11所示。从帕累托前沿选择出了两个达到优化目标的几何尺寸:xopt1=[14.748 5.189 5.904]T和xopt2=[14.733 6.011 6.401]T。如下表所示,和直接调用电磁场仿真优化相比,本发明节约了大量了优化时间。
Claims (2)
1.一种基于神经网络的多参数电磁场建模仿真方法,其特征在于,该方法结构上包含三个分支,其中第一分支用于仿真S参数,第二分支用于仿真增益,第三分支用于仿真方向图,具体包括以下步骤:
步骤1:神经网络-传递函数建模仿真方法的训练:
本方法采用的传递函数为:
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>Q</mi>
</munderover>
<mfrac>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mrow>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中,pi和ri分别为传递函数的极点和留数,i=1,…,Q,Q为传递函数的阶数,s在第一分支和第二分支为频率,在第三分支为角度;
每个分支的训练具体过程如下:
步骤1-1.将几何变量输入到电磁仿真软件,通过全波仿真得到对应于几何变量的天线的S参数、增益或方向图O={O1,...,OW},其中,W是频率采样点的个数;
步骤1-2.根据步骤1-1中得到的仿真结果,通过使用矢量拟合法,得到神经网络训练数集中每组数据所对应传递函数的极点和留数;根据传递函数的阶数,原始的训练数集被分为若干类Ck(k=1,2,…,K),其中K是各个分支的类别的个数;对于具有相同的传递函数阶数的训练样本被分为一类,每一类的阶数可以用Qk(k=1,2,…,K)表示;每一类只用于训练一个神经网络;
步骤1-3.通过神经网络的训练,得到传递函数的极点、留数和几何变量间的映射关系,进而得到传递函数的具体形式;带入自变量s,得到神经网络的输出O'={O'1,...,O'W};
步骤1-4.通过调整神经网络内部的权重和阈值,对不同的几何变量,使O和O'之间的误差最小化;
步骤2.神经网络-传递函数建模仿真方法的测试:
在测试过程中,将几何变量输入至对应的神经网络得到仿真结果。
2.根据权利要求1所述的基于神经网络的多参数电磁场建模仿真方法,其特征在于,步骤1-4后还包括如下步骤:步骤1-5.将几何变量输入到支持向量机,支持向量机的输出为每个训练样本的阶数Q'={Q'1,...,Q'K1},通过调整支持向量机内部的权重和阈值使Q'和传递函数实际的阶数{Q1,...,QK}之间的误差最小;
在步骤2前,还包括如下步骤:
将几何变量输入至训练好的支持向量机中,从而得到所对应的传递函数的阶数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711145631.2A CN107729694B (zh) | 2017-11-17 | 2017-11-17 | 一种基于神经网络的多参数电磁场建模仿真方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711145631.2A CN107729694B (zh) | 2017-11-17 | 2017-11-17 | 一种基于神经网络的多参数电磁场建模仿真方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107729694A true CN107729694A (zh) | 2018-02-23 |
CN107729694B CN107729694B (zh) | 2020-09-25 |
Family
ID=61216240
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711145631.2A Active CN107729694B (zh) | 2017-11-17 | 2017-11-17 | 一种基于神经网络的多参数电磁场建模仿真方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107729694B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5056039A (en) * | 1988-07-14 | 1991-10-08 | University Of Alabama In Huntsville | Holographic interconnect system |
CN102393884A (zh) * | 2011-10-12 | 2012-03-28 | 杭州电子科技大学 | 基于bp神经网络的热连轧电磁感应加热温度预测方法 |
CN102879704A (zh) * | 2012-10-08 | 2013-01-16 | 南京航空航天大学 | 一种高压输电线高阻接地故障多参数数据融合检测方法 |
CN104535865A (zh) * | 2014-12-30 | 2015-04-22 | 西安工程大学 | 基于多参数的电力变压器运行故障综合诊断方法 |
CN105787558A (zh) * | 2016-04-11 | 2016-07-20 | 江苏科技大学 | 基于ads的知识神经网络微带滤波器设计方法 |
CN106355245A (zh) * | 2016-09-12 | 2017-01-25 | 哈尔滨工业大学 | 一种基于神经网络算法的阵列天线方向图综合方法 |
CN106355249A (zh) * | 2016-08-29 | 2017-01-25 | 国网江苏省电力公司电力科学研究院 | 一种基于bp神经网络的继电保护设备状态识别方法 |
-
2017
- 2017-11-17 CN CN201711145631.2A patent/CN107729694B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5056039A (en) * | 1988-07-14 | 1991-10-08 | University Of Alabama In Huntsville | Holographic interconnect system |
CN102393884A (zh) * | 2011-10-12 | 2012-03-28 | 杭州电子科技大学 | 基于bp神经网络的热连轧电磁感应加热温度预测方法 |
CN102879704A (zh) * | 2012-10-08 | 2013-01-16 | 南京航空航天大学 | 一种高压输电线高阻接地故障多参数数据融合检测方法 |
CN104535865A (zh) * | 2014-12-30 | 2015-04-22 | 西安工程大学 | 基于多参数的电力变压器运行故障综合诊断方法 |
CN105787558A (zh) * | 2016-04-11 | 2016-07-20 | 江苏科技大学 | 基于ads的知识神经网络微带滤波器设计方法 |
CN106355249A (zh) * | 2016-08-29 | 2017-01-25 | 国网江苏省电力公司电力科学研究院 | 一种基于bp神经网络的继电保护设备状态识别方法 |
CN106355245A (zh) * | 2016-09-12 | 2017-01-25 | 哈尔滨工业大学 | 一种基于神经网络算法的阵列天线方向图综合方法 |
Non-Patent Citations (3)
Title |
---|
LI-YE XIAO 等: "Artificial neural network with data mining techniques for antenna design", 《2017 IEEE INTERNATIONAL SYMPOSIUM ON ANTENNAS AND PROPAGATION》 * |
XIAOLEI DING 等: "Neural-network approaches to electromagnetic-based modeling of passive components and their applications to high-frequency and high-speed nonlinear circuit optimization", 《IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 》 * |
林培群 等: "多分支BP网络模型及其在车型分类中的应用", 《微计算机信息》 * |
Also Published As
Publication number | Publication date |
---|---|
CN107729694B (zh) | 2020-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104992001B (zh) | 大规模mimo阵列天线远场辐射场的精确快速计算方法 | |
CN104199536B (zh) | 一种基于bp神经网络的fpga动态功耗估算方法 | |
CN104112051B (zh) | 基于单元形函数的反射面天线机电集成设计方法 | |
CN104485513B (zh) | 宽频带天线阵列实时综合方法 | |
CN102708229A (zh) | 复杂分层媒质结构的矩阵分解结合新奇异值分解方法 | |
CN106777620A (zh) | 一种用于功率晶体管的神经网络空间映射建模方法 | |
Zhao et al. | Displacement discontinuity analysis of a nonlinear interfacial crack in three-dimensional transversely isotropic magneto-electro-elastic bi-materials | |
CN107992696A (zh) | 一种复杂色散媒质中的改进的指数时间积分构造方法 | |
CN107748806A (zh) | 微波热解污泥电磁场强度及其分布规律的计算和模拟方法 | |
CN106372348B (zh) | 一种线性系统中基于误差控制的矢量拟合模型降阶方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN107729694A (zh) | 一种基于神经网络的多参数电磁场建模仿真方法 | |
CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN106815380A (zh) | 一种提取寄生电阻的方法及系统 | |
CN107305536A (zh) | 混合阶时域不连续伽略金方法 | |
CN111898320A (zh) | 基于神经网络的射频功放散射参数提取方法及装置 | |
CN104778293B (zh) | 非均匀介质目标电磁散射的体积分Nystrom分析方法 | |
CN106815379A (zh) | 一种提取寄生电容的方法及系统 | |
CN105205299A (zh) | 电大目标电磁散射特性快速降维分析方法 | |
Adedayo et al. | Optimetric analysis of 1x4 array of circular microwave patch antennas for mammographic applications using adaptive gradient descent algorithm | |
CN114491992A (zh) | 一种基于等效偶极矩和物理光学法的高效电磁散射方法 | |
Qi et al. | Numerical dispersion compensation for FDTD via deep learning | |
Kaya et al. | ANFIS models for synthesis of micro‐coplanar stripline and asymmetric coplanar stripline with an infinitely wide strip | |
CN106156431B (zh) | 基于nurbs曲面建模的导体目标电磁散射仿真方法 | |
CN105095154A (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 |