CN106053945B - 一种短时s函数史密斯圆图分析装置及方法 - Google Patents
一种短时s函数史密斯圆图分析装置及方法 Download PDFInfo
- Publication number
- CN106053945B CN106053945B CN201610520360.3A CN201610520360A CN106053945B CN 106053945 B CN106053945 B CN 106053945B CN 201610520360 A CN201610520360 A CN 201610520360A CN 106053945 B CN106053945 B CN 106053945B
- Authority
- CN
- China
- Prior art keywords
- function
- measured piece
- term
- short
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/16—Measuring impedance of element or network through which a current is passing from another source, e.g. cable, power line
Abstract
本发明公开了一种短时S函数史密斯圆图分析装置及方法，属于测试领域，本发明将网络分析仪中的固态史密斯圆图分析分解为相关电距离的瞬态史密斯圆图，用以分析电缆、连接器、夹具等器部件连阻抗的连续性变化，可以便捷、直观、准确的进行连续阻抗测试，快速直观的给出连续阻抗变换，并补偿了常规方法的测试误差。
Description
技术领域
本发明属于测试领域，具体涉及一种短时S函数史密斯圆图分析装置及方法。
背景技术
目前矢量网络分析仪中，使用“黑箱理论”，将电缆、连接器、夹具等具有一定电长度的器部件统一当做被测件测试一组S参数，对连续阻抗的变换一般采用频时分析方法，观察脉冲响应函数和阶跃响应函数，不能直观的表达阻抗的连续变化，且响应幅度有误差。
发明内容
针对现有技术中存在的上述技术问题，本发明提出了一种短时S函数史密斯圆图分析装置及方法，设计合理，克服了现有技术的不足，可以便捷、直观、准确的进行连续阻抗测试。
为了实现上述目的，本发明采用如下技术方案：
一种短时S函数史密斯圆图分析装置，包括信号发生器、耦合器、接收机和负载；
所述信号发生器，被配置为用于产生激励信号；
所述接收机，包括第一I/Q接收机和第二I/Q接收机；
所述第一I/Q接收机，被配置为用于测试信号发生器的发射信号；
所述第二I/Q接收机，被配置为用于测试被测件的反射信号；
所述耦合器，被配置为用于从信号发生器的激励信号耦合部分信号送入第一I/Q接收机，以及从被测件的反射信号耦合部分信号送入第二I/Q接收机；
所述负载，被配置为用于吸收被测件的输出信号；
信号发生器输出激励信号，其中一小部分激励信号经过耦合器的耦合进入第一I/Q接收机，大部分激励信号经过被测件的输入端进入被测件，其中产生的反射信号经过被测件的输入端输出，经过耦合器的耦合进入第二I/Q接收机，产生的输出信号经过被测件的输出端加载至负载被负载吸收。
此外，本发明还提到一种短时S函数史密斯圆图分析方法，该方法采用如上所述的一种短时S函数史密斯圆图分析装置，包括如下步骤：
步骤1：测试被测件的反射S函数，具体包括如下步骤：
步骤1.1：信号发生器发射信号激励被测件；
步骤1.2：通过第一I/Q接收机测试信号发生器的发射信号的矢量，读数记为R(f0+nΔf)；
步骤1.3：通过第二I/Q接收机测试被测件的反射信号的矢量，读数记为A(f0+nΔf)；
步骤1.4：被测件的反射S函数表示为S′11(f0+nΔf)＝A(f0+nΔf)/R(f0+nΔf)；
其中，f0表示测试的初始频率，Δf表示频率变化间隔；
步骤2：根据公式(1)对被测件的反射S函数进行包括时间轴在内的修正，修正结果表示为
S11(f0+nΔf)＝FFT([IFFT(S′11(f0+nΔf))*δ(t-t0)]·u(t)) (1)；
其中，FFT(·)为快速傅里叶变换计算；IFFT(·)为逆傅里叶变换计算；*为卷积计算；δ(t-t0)为单位脉冲函数；t0是测试空置，即信号发生器不连接被测件和负载时最大波峰对应的时延值；u(t)为单位阶跃函数；
步骤3：采用频时转换技术并根据公式(2)求解脉冲响应Pluse(kΔTs)
其中，P＝A0 exp(jθ0)，A0为取样点的半径长度；θ0为取样点的相角；W0为螺旋线的伸展曲率；为相邻两取样点间的等分角，为变换的中间变量；N为取样点数；ΔTs＝L/[C(N-1)]为采样的时间间隔；L为被测件的物理长度；C为被测件的微波传播速度；
步骤4：根据公式(3)对脉冲响应Pluse(kΔTs)进行修正，得到修正数据Pluse′(kΔTs)
Pluse′(kΔTs)＝Pluse(kΔTs)·B(kΔTs) (3)；
其中，B(kΔTs)＝σΔTs/C，B(kΔTS)为脉冲响应的修正系数，σ为被测件单位长度的损耗；
步骤5：设置窗函数及滑动门函数，采用短时S函数分析显示动态史密斯圆图。
优选地，在步骤5中，具体包括
步骤5.1：设置窗函数Window及滑动门函数Gateq(n,ΔN,Δt)，Window选用汉宁窗或贝塞尔窗，滑动门函数Gateq(n,ΔN,Δt)的取值如下：
其中，n＝0,1,2,...,N-1，ΔN为门的长度，Δt为分析的步长，t0为拟分析的初始时刻，q＝0,1,2,...为自然数，Gateq表示第q个滑动门；
步骤5.2：令q＝0，分析被测件在t0C处的S函数；
步骤5.3：令短时脉冲函数Pluseshort＝Pluse′(kΔTs)·Window·Gateq(n,ΔN,Δt)；
步骤5.4：根据步骤5.3中的脉冲函数Pluseshort计算被测件在(t0+qΔt)C处的短时S函数
短时S函数＝FFT(Pluseshort)；其中，FFT为傅里叶变换；
步骤5.5：根据步骤5.4中的短时S函数计算并显示被测件在(t0+qΔt)C处的短时史密斯圆图，史密斯圆图为短时S函数的极坐标表示；
步骤5.6：判断Δt与Δtterm-ΔTsΔN/2的大小关系；其中，Δtterm为分析终止长度的对应时刻，设置为Δtterm＝L/C；
若：判断结果是Δt≤Δtterm-ΔTsΔN/2，则令q＝q+1，然后执行步骤5.3；
或判断结果是Δt＞Δtterm-ΔTsΔN/2，则分析结束。
本发明所带来的有益技术效果：
本发明提出了一种短时S函数史密斯圆图分析装置及方法，与现有技术相比，本发明将网络分析仪中的固态史密斯圆图分析分解为相关电距离的瞬态史密斯圆图，用以分析电缆、连接器、夹具等器部件连阻抗的连续性变化，可以便捷、直观、准确的进行连续阻抗测试，快速直观的给出连续阻抗变换，并补偿了常规方法的测试误差。
附图说明
图1为本发明一种短时S函数史密斯圆图分析装置的硬件原理图。
图2为本发明一种短时S函数史密斯圆图分析方法的流程框图。
图3为本发明双端口非线性网络的示意图。
图4为采用短时S函数分析显示动态史密斯圆的流程框图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
实施例1：
如图1所示的一种短时S函数史密斯圆图分析装置，包括信号发生器、耦合器、接收机和负载；
所述信号发生器，被配置为用于产生激励信号；
所述接收机，包括第一I/Q接收机和第二I/Q接收机；
所述第一I/Q接收机，被配置为用于测试信号发生器的发射信号；
所述第二I/Q接收机，被配置为用于测试被测件的反射信号；
所述耦合器，被配置为用于从信号发生器的激励信号耦合部分信号送入第一I/Q接收机，以及从被测件的反射信号耦合部分信号送入第二I/Q接收机；
所述负载，被配置为用于吸收被测件的输出信号；
信号发生器输出激励信号，其中一小部分激励信号经过耦合器的耦合进入第一I/Q接收机，大部分激励信号经过被测件的输入端进入被测件，其中产生的反射信号经过被测件的输入端输出，经过耦合器的耦合进入第二I/Q接收机，产生的输出信号经过被测件的输出端加载至负载被负载吸收。
实施例2：
在上述实施例的基础上，本发明还提到一种短时S函数史密斯圆图分析方法(如图2所示)，包括如下步骤：
步骤1：测试被测件的反射S函数
此步骤为借用网络参数分析理论，双端口非线性网络的示意图如图3所示，校准后，使用网络分析仪测试被测件的S参数。S参数应表达为对应频率下的幅度与相位响应格式，采用对数幅度，线性相位，保存S参数。这里推荐使用S11。
具体包括如下步骤：
步骤1.1：信号发生器发射信号激励被测件；
步骤1.2：通过第一I/Q接收机测试信号发生器的发射信号的矢量，读数记为R(f0+nΔf)；
步骤1.3：通过第二I/Q接收机测试被测件的反射信号的矢量，读数记为A(f0+nΔf)；
步骤1.4：被测件的反射S函数表示为S′11(f0+nΔf)＝A(f0+nΔf)/R(f0+nΔf)；
其中，f0表示测试的初始频率，Δf表示频率变化间隔；
步骤2：根据公式(1)对被测件的反射S函数进行包括时间轴在内的修正，修正结果表示为
S11(f0+nΔf)＝FFT([IFFT(S′11(f0+nΔf))*δ(t-t0)]·u(t)) (1)；
其中，FFT(·)为快速傅里叶变换计算；IFFT(·)为逆傅里叶变换计算；*为卷积计算；δ(t-t0)为单位脉冲函数；t0是测试空置，即信号发生器不连接被测件和负载时最大波峰对应的时延值；u(t)为单位阶跃函数；
步骤3：采用频时转换技术并根据公式(2)求解脉冲响应Pluse(kΔTs)
使用时频转换技术求得微波器部件的脉冲响应。因本发明目的在于进行短时动态分析，因此选用CHIRPZ变换代替常用的IFFT变换。
其中，P＝A0exp(jθ0)，A0为取样点的半径长度；θ0为取样点的相角；W0为螺旋线的伸展曲率；为相邻两取样点间的等分角，为变换的中间变量；N为取样点数；ΔTs＝L/[C(N-1)]为采样的时间间隔；L为被测件的物理长度；C为被测件的微波传播速度。
步骤4：脉冲响应修正。
得到的滑动脉冲响应需要进行修正。
根据公式(3)对脉冲响应Pluse(kΔTs)进行修正，得到修正数据Pluse′(kΔTs)
Pluse′(kΔTs)＝Pluse(kΔTs)·B(kΔTs) (3)；
其中，B(kΔTs)＝σΔTs/C，B(kΔTS)为脉冲响应的修正系数，σ为被测件单位长度的损耗；
步骤5：设置窗函数及滑动门函数，采用短时S函数分析显示动态史密斯圆，具体包括如下步骤(如图4所示)：
步骤5.1：设置窗函数Window及滑动门函数Gateq(n,ΔN,Δt)，Window选用汉宁窗或贝塞尔窗，滑动门函数Gateq(n,ΔN,Δt)的取值如下：
其中，n＝0,1,2,...,N-1，ΔN为门的长度，Δt为分析的步长，t0为拟分析的初始时刻， q＝0,1,2,...为自然数，Gateq表示第q个滑动门；
步骤5.2：令q＝0，分析被测件在t0C处的S函数；
步骤5.3：令短时脉冲函数Pluseshort＝Pluse′(kΔTs)·Window·Gateq(n,ΔN,Δt)；
步骤5.4：根据步骤5.3中的脉冲函数Pluseshort计算被测件在(t0+qΔt)C处的短时S函数
短时S函数＝FFT(Pluseshort)；其中，FFT为傅里叶变换；
步骤5.5：根据步骤5.4中的短时S函数计算并显示被测件在(t0+qΔt)C处的短时史密斯圆图，史密斯圆图为短时S函数的极坐标表示；
步骤5.6：判断Δt与Δtterm-ΔTsΔN/2的大小关系；其中，Δtterm为分析终止长度的对应时刻，设置为Δtterm＝L/C；
若：判断结果是Δt≤Δtterm-ΔTsΔN/2，则令q＝q+1，然后执行步骤5.3；
或判断结果是Δt＞Δtterm-ΔTsΔN/2，则分析结束。
本发明一种短时S函数史密斯圆图分析方法，将网络分析仪中的固态史密斯圆图分析分解为相关电距离的瞬态史密斯圆图，用以分析电缆、连接器、夹具等器部件连阻抗的连续性变化，可以便捷、直观、准确的进行连续阻抗测试，快速直观的给出连续阻抗变换，并补偿了常规方法的测试误差。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。
Claims (3)
1.一种短时S函数史密斯圆图分析装置，其特征在于：包括信号发生器、耦合器、接收机和负载；
所述信号发生器，被配置为用于产生激励信号；
所述接收机，包括第一I/Q接收机和第二I/Q接收机；
所述第一I/Q接收机，被配置为用于测试信号发生器的发射信号；
所述第二I/Q接收机，被配置为用于测试被测件的反射信号；
所述耦合器，被配置为用于从信号发生器的激励信号耦合部分信号送入第一I/Q接收机，以及从被测件的反射信号耦合部分信号送入第二I/Q接收机；
所述负载，被配置为用于吸收被测件的输出信号；
信号发生器输出激励信号，其中一小部分激励信号经过耦合器的耦合进入第一I/Q接收机，大部分激励信号经过被测件的输入端进入被测件，其中产生的反射信号经过被测件的输入端输出，经过耦合器的耦合进入第二I/Q接收机，产生的输出信号经过被测件的输出端加载至负载被负载吸收。
2.一种短时S函数史密斯圆图分析方法，其特征在于：采用如权利要求1所述的一种短时S函数史密斯圆图分析装置，包括如下步骤：
步骤1：测试被测件的反射S函数，具体包括如下步骤：
步骤1.1：信号发生器发射信号激励被测件；
步骤1.2：通过第一I/Q接收机测试信号发生器的发射信号的矢量，读数记为R(f0+nΔf)；
步骤1.3：通过第二I/Q接收机测试被测件的反射信号的矢量，读数记为A(f0+nΔf)；
步骤1.4：被测件的反射S函数表示为S′11(f0+nΔf)＝A(f0+nΔf)/R(f0+nΔf)；
其中，f0表示测试的初始频率，Δf表示频率变化间隔；
步骤2：根据公式(1)对被测件的反射S函数进行包括时间轴在内的修正，修正结果表示为
S11(f0+nΔf)＝FFT([IFFT(S′11(f0+nΔf))*δ(t-t0)]·u(t)) (1)；
其中，FFT(·)为快速傅里叶变换计算；IFFT(·)为逆傅里叶变换计算；*为卷积计算；δ(t-t0)为单位脉冲函数；t0是测试空置，即信号发生器不连接被测件和负载时最大波峰对应的时延值；u(t)为单位阶跃函数；
步骤3：采用频时转换技术并根据公式(2)求解脉冲响应Pluse(kΔTs)
其中， P＝A0exp(jθ0)，A0为取样点的半径长度；θ0为取样点的相角；W0为螺旋线的伸展曲率；为相邻两取样点间的等分角，为变换的中间变量；N为取样点数；ΔTs＝L/[C(N-1)]为采样的时间间隔；L为被测件的物理长度；C为被测件的微波传播速度；
步骤4：根据公式(3)对脉冲响应Pluse(kΔTs)进行修正，得到修正数据Pluse′(kΔTs)
Pluse′(kΔTs)＝Pluse(kΔTs)·B(kΔTs)(3)；
其中，B(kΔTs)＝σΔTs/C，B(kΔTS)为脉冲响应的修正系数，σ为被测件单位长度的损耗；
步骤5：设置窗函数及滑动门函数，采用短时S函数分析显示动态史密斯圆图。
3.根据权利要求2所述的短时S函数史密斯圆图分析方法，其特征在于：在步骤5中，具体包括
步骤5.1：设置窗函数Window及滑动门函数Gateq(n,ΔN,Δt)，Window选用汉宁窗或贝塞尔窗，滑动门函数Gateq(n,ΔN,Δt)的取值如下：
其中，n＝0,1,2,...,N-1，ΔN为门的长度，Δt为分析的步长，t0为拟分析的初始时刻，q＝0,1,2,...为自然数，Gateq表示第q个滑动门；
步骤5.2：令q＝0，分析被测件在t0C处的S函数；
步骤5.3：令短时脉冲函数Pluseshort＝Pluse′(kΔTs)·Window·Gateq(n,ΔN,Δt)；
步骤5.4：根据步骤5.3中的脉冲函数Pluseshort计算被测件在(t0+qΔt)C处的短时S函数短时S函数＝FFT(Pluseshort)；其中，FFT为傅里叶变换；
步骤5.5：根据步骤5.4中的短时S函数计算并显示被测件在(t0+qΔt)C处的短时史密斯圆图，史密斯圆图为短时S函数的极坐标表示；
步骤5.6：判断Δt与Δtterm-ΔTsΔN/2的大小关系；其中，Δtterm为分析终止长度的对应时刻，设置为Δtterm＝L/C；
若：判断结果是Δt≤Δtterm-ΔTsΔN/2，则令q＝q+1，然后执行步骤5.3；
或判断结果是Δt＞Δtterm-ΔTsΔN/2，则分析结束。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610520360.3A CN106053945B (zh) | 2016-07-05 | 2016-07-05 | 一种短时s函数史密斯圆图分析装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610520360.3A CN106053945B (zh) | 2016-07-05 | 2016-07-05 | 一种短时s函数史密斯圆图分析装置及方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106053945A CN106053945A (zh) | 2016-10-26 |
| CN106053945B true CN106053945B (zh) | 2019-01-22 |
Family
ID=57201858
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610520360.3A Active CN106053945B (zh) | 2016-07-05 | 2016-07-05 | 一种短时s函数史密斯圆图分析装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106053945B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107203495B (zh) * | 2017-06-09 | 2020-12-22 | 中国电子科技集团公司第四十一研究所 | 一种s参数圆形插值方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2514902A1 (fr) * | 1981-10-20 | 1983-04-22 | Cables De Lyon Geoffroy Delore | Procede de mesure de l'impedance d'une prise de terre et dispositif en faisant application |
| US5274333A (en) * | 1992-03-17 | 1993-12-28 | La Corporation De L'ecole Polytechnique | Frequency balanced six-port reflectometer with a variable test port impedance |
| CN102540106A (zh) * | 2012-01-30 | 2012-07-04 | 江苏省计量科学研究院 | 新能源逆变系统直流侧传导emi噪声测量装置及测校方法 |
| CN103063926A (zh) * | 2012-12-19 | 2013-04-24 | 中国电力科学研究院 | 一种低压电力线窄带载波接入阻抗的测试方法及测试装置 |
| CN105137196A (zh) * | 2015-08-25 | 2015-12-09 | 浪潮电子信息产业股份有限公司 | 一种简易的sas、sata信号走线阻抗测试方法 |
-
2016
- 2016-07-05 CN CN201610520360.3A patent/CN106053945B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2514902A1 (fr) * | 1981-10-20 | 1983-04-22 | Cables De Lyon Geoffroy Delore | Procede de mesure de l'impedance d'une prise de terre et dispositif en faisant application |
| US5274333A (en) * | 1992-03-17 | 1993-12-28 | La Corporation De L'ecole Polytechnique | Frequency balanced six-port reflectometer with a variable test port impedance |
| CN102540106A (zh) * | 2012-01-30 | 2012-07-04 | 江苏省计量科学研究院 | 新能源逆变系统直流侧传导emi噪声测量装置及测校方法 |
| CN103063926A (zh) * | 2012-12-19 | 2013-04-24 | 中国电力科学研究院 | 一种低压电力线窄带载波接入阻抗的测试方法及测试装置 |
| CN105137196A (zh) * | 2015-08-25 | 2015-12-09 | 浪潮电子信息产业股份有限公司 | 一种简易的sas、sata信号走线阻抗测试方法 |
Non-Patent Citations (3)
| Title |
|---|
| 同轴电缆特性阻抗测试方法及应用研究;周军霞 等;《机电工程技术》;20101231;第39卷(第06期);第35-38页 |
| 基于史密斯圆图分析射频放大器的设计;罗建国 等;《武汉大学学报（理学版）》;20050228;第51卷(第1期);第123-128页 |
| 应用史密斯圆图提取慢波微带线特征阻抗方法;董峦;《电子测试》;20110531(第5期);第94-96页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106053945A (zh) | 2016-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107390100B (zh) | 一种基于时间反演的电力电缆局部放电定位方法 | |
| Aiello et al. | A chirp-z transform-based synchronizer for power system measurements | |
| CN103941087B (zh) | 欠采样速率下的高频余弦信号的频率测量方法及其装置 | |
| CN109782063A (zh) | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 | |
| CN101871975A (zh) | 电缆转移阻抗时域测试系统及测试方法 | |
| CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
| US9823281B2 (en) | Time frame measuring method with calibration in the frequency range | |
| CN108318852A (zh) | 一种智能电能表方波影响试验方法 | |
| CN104270208B (zh) | 一种远端射频单元rru驻波比检测的方法及装置 | |
| CN103760425A (zh) | 一种基于时域准同步的介损角快速测量方法及其测量装置 | |
| CN104406509A (zh) | 基于hht电力电缆长度测量方法 | |
| CN105572480A (zh) | 原位测试双导体形式线缆的宽频带传输线参数的方法 | |
| CN102998528B (zh) | 动态正弦畸变信号中频率分量的提取方法及其实现装置 | |
| CN105353330A (zh) | 一种基于虚拟仪器技术的电子式电流互感器在线校验系统及其算法 | |
| CN106053945B (zh) | 一种短时s函数史密斯圆图分析装置及方法 | |
| CN111624405A (zh) | 一种基于小波去噪的长线屏蔽线缆转移阻抗测量系统及方法 | |
| Hashmi et al. | Determining wave propagation characteristics of MV XLPE power cable using time domain reflectometry technique | |
| CN106483445B (zh) | 宽带电路相位非线性失真的一种内建测量方法及装置 | |
| CN103439580B (zh) | 大尺度电气设备阻抗宽频特性的时域测量方法 | |
| Analyzers | THE UNIVERSITY OF BRITISH COLUMBIA | |
| Rodrigues et al. | Low-cost embedded measurement system for power quality frequency monitoring | |
| CN115128404A (zh) | 一种非接触式电缆故障定位方法 | |
| CN105158562A (zh) | 一种用于整机雷电间接效应试验的时频域数据处理方法 | |
| CN106324540A (zh) | 一种数字示波器频率响应的快速测量方法及数字示波器 | |
| CN104849575B (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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |