CN113346901B - 一种交替并行采样系统变温条件下修正系数自动适配方法 - Google Patents
一种交替并行采样系统变温条件下修正系数自动适配方法 Download PDFInfo
- Publication number
- CN113346901B CN113346901B CN202110590735.4A CN202110590735A CN113346901B CN 113346901 B CN113346901 B CN 113346901B CN 202110590735 A CN202110590735 A CN 202110590735A CN 113346901 B CN113346901 B CN 113346901B
- Authority
- CN
- China
- Prior art keywords
- frequency
- temperature
- correction
- calibration
- under
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005070 sampling Methods 0.000 title claims abstract description 32
- 230000006978 adaptation Effects 0.000 title claims abstract description 17
- 230000003044 adaptive effect Effects 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims abstract description 7
- 230000007613 environmental effect Effects 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 13
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 238000011088 calibration curve Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/06—Continuously compensating for, or preventing, undesired influence of physical parameters
- H03M1/08—Continuously compensating for, or preventing, undesired influence of physical parameters of noise
- H03M1/089—Continuously compensating for, or preventing, undesired influence of physical parameters of noise of temperature variations
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M1/00—Analogue/digital conversion; Digital/analogue conversion
- H03M1/10—Calibration or testing
- H03M1/1009—Calibration
- H03M1/1014—Calibration at one point of the transfer characteristic, i.e. by adjusting a single reference value, e.g. bias or gain error
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
- Analogue/Digital Conversion (AREA)
Abstract
本发明公开了波形数字化领域的一种交替并行采样系统变温条件下修正系数自动适配方法,包括步骤:步骤1:预先获取在不同温度下,频带范围内的每一个单频点信号对应的输入信号频率下的通道失配误差标定结果;步骤2:根据当前环境温度从预先获取的失配误差标定结果中线性插值,得到当前温度下的失配误差系数;步骤3:通过硬件逻辑计算当前温度下的适配误差系数对应的修正系数,对当前温度下的测试结果实时修正。本发明使用经过温度适配后的修正参数对TIADC系统进行修正可使系统在变温下的宽带动态性能基本保持不变,从而使TIADC可以适应变温的工作环境。此外,标定过程通过硬件自动实现,有效地提高了标定效率,为TIADC系统的多通道应用提供了便利。
Description
技术领域
本发明涉及波形数字化领域,具体是一种交替并行采样系统变温条件下修正系数自动适配方法。
背景技术
波形数字化技术是高速信号测量领域的就热点,其中一种方法是通过高速模拟-数字转换器(Analog-to-Digital Convertor,ADC)直接对波形进行实时采样。如图1所示,在波形数字化系统中,通过并行交替采样(TIADC)技术可使系统采样率突破ADC采样率极限,通过将M个采样率为Fs的ADC的采样时刻以1/(MFs)等间隔排布,即可得到等效采样率为MFs的波形数字化系统。
而在TIADC系统中,不同ADC之间存在不一致性;这种不一致性会导致系统的动态性能,如信噪比、有效位性能等相比于单个ADC的性能会产生下降。由ADC不一致性导致的误差称为失配误差,具体包括增益误差、相位误差和偏置误差三种。
在实际应用中,TIADC系统必须进行失配误差的修正,一般地可采用自适应修正方法与基于前景校准的数字修正方法。前者能够跟随失配误差的变化而更正修正系数,但要求输入信号是广义平稳的,但在实际应用中信号往往在时间和幅度上都是随机的,短时间尺度内不能满足自适应修正方法的要求。后者可以在失配误差参数不变的情况下实现失配误差修正,但是,在实际应用中,经常会遇到失配误差系数变化的情况。一种较为常见的情况就是变温导致系统频率响应变化引起的失配误差系数变化。当失配误差系数变化时,如仍采用原有的修正参数,会导致系统的动态性能明显下降。例如,对于一个20Gsps,12bit,由8个并行采样通道组成的TIADC系统,其中4个通道的相位误差参数温漂约为0.015ps/℃。如图3所示,若采用常温(25℃)下的修正参数对0℃下的采样结果修正,其有效位性能相对于常温下经失配误差修正后的系统性能相差约1bit。
发明内容
本发明的目的在于提供一种交替并行采样系统变温条件下修正系数自动适配方法,通过失配误差系数的温度线性插值、基于FPGA的失配误差系数自动标定以及完美重构修正系数的自动计算实现变温条件下修正参数的自动适配,进而实现失配误差修正。
为实现上述目的,本发明提供如下技术方案:
一种交替并行采样系统变温条件下修正系数自动适配方法,包括步骤:
步骤1:预先获取在不同温度下,频带范围内的每一个单频点信号对应的输入信号频率下的通道失配误差标定结果;
步骤2:根据当前环境温度从预先获取的失配误差标定结果中线性插值,得到当前温度下的失配误差系数;
步骤3:通过硬件逻辑计算当前温度下的适配误差系数对应的修正系数,对当前温度下的测试结果实时修正。
在一些实施例中,步骤1具体包括:
步骤1.1:设计自动标定硬件逻辑,存储某一单频点信号下的ADC的采样数据,顺次对ADC各通道的采样数据进行四参数拟合,将拟合结果转化为相同形式;
步骤1.2:设定ADC某一通道为基准,基准的增益误差为1,相位误差为0,根据拟合结果得到ADC各通道相对于基准的失配误差标定结果,将失配误差标定结果写入存储,完成当前频点下的标定;
步骤1.3:改变输入信号频率,重复步骤1.1-1.2,完成设定的频带范围内所有频点下的ADC各通道的失配误差标定;
步骤1.4:改变温度,重复步骤1.3,完成在不同温度下ADC各通道的失配误差标定。
在一些实施例中,步骤3中修正系数的计算通过完美重构法实现。
在一些实施例中,步骤3具体包括:
步骤3.1:启动修正系数计算逻辑,读入当前频率、当前温度对应的适配误差系数;
步骤3.2:计算完美重构修正所需的频点的角频率,对每个角频率下的完美重构矩阵的各元素进行求解计算,获得当前角频率下的修正滤波器的频域数值解;
步骤3.3:对频域数值解进行IFFT变换,得到对应的修正滤波器的时域数值解;
步骤3.4:加窗截断处理,得到修正滤波器的系数,可实现对当前温度下的测试结果实时修正。
在一些实施例中,步骤3.2中,频点的数量N=2P。
有益效果:本发明采用温度插值可得到当前系统所处的环境温度下对应的失配误差系数,并根据此系数计算获得合适的修正参数。经实验,本发明使用经过温度适配后的修正参数对TIADC系统进行修正可使系统在变温下的宽带动态性能基本保持不变,从而使TIADC可以适应变温的工作环境。此外,标定过程通过硬件自动实现,有效地提高了标定效率,为TIADC系统的多通道应用提供了便利。
附图说明
图1为TIADC系统的采样示意图;
图2为TIADC系统采用完美重构算法的采样修正过程图;
图3为某一TIADC系统在不同温度下的修正效果图;
图4为本发明自动标定逻辑流程示意图;
图5为本发明某一TIADC系统在不同温度下的增益误差系数标定曲线;
图6是本发明某一TIADC系统在不同温度下的相位误差系数标定曲线;
图7是本发明通过硬件语言的高层次综合编程得到的失配误差计算逻辑流程示意图;
图8为本发明的实施方式中修正参数经温度自动适配后的修正结果对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种交替并行采样系统变温条件下修正系数自动适配方法,主要分为自动标定过程、变温失配误差插值以及修正系数自动计算过程,分别对应步骤1、2、3。
步骤1:预先获取在不同温度下,频带范围内的每一个单频点信号对应的输入信号频率下的通道失配误差标定结果。
具体地,启动自动标定逻辑,通过信号源输入频率为fi的单频点信号并将输入频率通过寄存器写入逻辑。在该频率下将各通道的若干个连续的采样点存储在EEPROM中。通过硬件逻辑进行四参数拟合得到当前频率下的通道失配误差标定结果。改变输入信号频率,使其覆盖所关心的信号频率范围,完成当前温度下的标定。
改变温度,重复上述操作,可得到不同温度下的标定结果。将标定温度和对应的标定频率以及失配误差标定结果写入EEPROM中,完成失配误差自动标定过程。
步骤2:在实际的应用环境中,根据温度传感器获取当前的环境温度,也即系统的工作温度。根据当前环境温度从EEPROM内存储的失配误差标定结果中线性插值,得到当前温度下的失配误差系数。
步骤3:启动失配误差计算流程,通过硬件逻辑计算当前温度下的适配误差系数对应的修正系数,对当前温度下的测试结果实时修正。本步骤中修正系数的计算通过完美重构法实现,采用完美重构法既能实现宽带修正同时又适合硬件的求解,在修正时能够维持系统的宽带动态性能不变,这是其他计算方法所不能实现的。
实施例,以一个由8通道(M=8)构成的20Gsps 12位的TIADC系统为例,对上述方法步骤进行具体说明。
步骤1.1:设计自动标定硬件逻辑,其状态机如4所示。将标定频率ω下的ADC的采样数据存储在RAM中,启动逻辑状态机,按顺序将ADC的1-8个通道的采样数据进行四参数拟合。具体地,首先将第1个ADC子通道的1024个采样数据读入启动四参数拟合过程,采样序列可写作向量x:
根据IEEE数字波形记录标准中的正弦波四参数拟合过程,通过迭代求解将采样序列x拟合为同一形式,即正交形式的拟合结果:
xm_fit(t)=Acos(ωt)+Bsin(ωt)+DC
通过迭代和最小二乘求解超定方程组过程依次可得到ADC各通道如上式形式的拟合结果。四参数拟合具体的求解方法为现有技术,本申请仅将其求解过程采用硬件化实现,在此不对其具体求解过程进行说明。
步骤1.2:根据得到的四参数拟合结果求得各通道之间的失配误差。以Chn0通道为基准,设增益误差g0=1,相位误差Δt0=0。将各通道四参数拟合结果转化为如下形式:
xm_fit(t)=Gmcos(ωi(t+tm))+DCm_i
据拟合结果得到ADC各通道(Chn1、Chn2……Chn7)相对于基准的失配误差,将标定结果写入EEPROM,完成当前频点下的标定。
步骤1.3:通过改变信号源输出频率,即输入信号频率,重复步骤1.1-1.2,实现当前温度下的所设定的信号频率区间内所有特征频点的标定。
步骤1.4:改变系统的工作温度,重复步骤1.3,完成在不同温度下ADC各通道的失配误差标定。
在本实例中,以-10℃、10℃、25℃三个标定温度点为例,得到其中第1通道相对第0通道(基准)的增益误差和相位误差随温度变化的曲线如图5和图6所示。由于偏置误差为常量,可通过计算系统实际测量时的噪声谱平均值得到,因此无需单独存储标定结果。
步骤2:当系统需要计算修正系数时,通过启动信号启动图7所示的状态机,首先从EEPROM中读入各温度下的各特征频点的失配误差标定结果,再根据当前系统所处的环境温度,线性插值得到当前温度下的失配误差系数。
步骤3.1:启动修正系数计算逻辑,通过求解完美重构方程组进行修正系数的自动计算过程,读入当前频率、当前温度对应的适配误差系数;系统通过完美重构法进行采样修正过程如图2所示。
步骤3.2:计算完美重构修正所需的频点的角频率,在每个角频率下对完美重构矩阵的各元素进行计算。频点的数量N需要足够的充分,使计算结果在一定的误差范围内满足之后的修正精度要求,一般地,至少需要满足N=2P个。
每个角频率下可得到一个M×M的完美重构矩阵,以第一Nyquist区为例,通过消元法求解如下方程组:
即可得到当前角频率下的修正滤波器频域数值解Fm[jω]。上述方程中,各参数意义如下:
M:ADC通道个数;
第m通道的增益误差随频率变化的函数;
第m通道的增益误差随频率变化的函数;
Ts:系统等最终等效的采样周期;
c:为一固定增益系数,可自行设置;
d:为一固定延迟系数,可自行设置,一般为滤波器阶数的一半。
其中各频率下的增益误差g和相位误差Δt的值由上述温度插值得到的失配误差标定序列再进行频率插值得到。目标函数中c=1,采用80阶滤波器修正时对应的延迟d=40。
步骤3.3:将上述得到的滤波器频域解进行2P个点的IFFT变换,得到对应的修正滤波器的时域解;
步骤3.4:对信号时域加窗和截断处理,得到80阶修正滤波器的系数,即可用于TIADC系统在当前温度下的测试结果实时修正。
图8所示为在0℃下,修正参数经过温度适配后得到的系统宽带动态性能测试结果,表明其有效位基本与常温下保持一致。可见本发明使用经过温度适配后的修正参数对TIADC系统进行修正可使系统在变温下的宽带动态性能基本保持不变,从而使TIADC可以适应变温的工作环境。此外,标定过程通过硬件自动实现,有效的提高了标定效率,为TIADC系统的多通道应用提供了便利。
虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
故以上所述仅为本申请的较佳实施例,并非用来限定本申请的实施范围;即凡依本申请的权利要求范围所做的各种等同变换,均为本申请权利要求的保护范围。
Claims (3)
1.一种交替并行采样系统变温条件下修正系数自动适配方法,其特征在于,包括步骤:
步骤1:预先获取在不同温度下,频带范围内的每一个单频点信号对应的输入信号频率下的通道失配误差标定结果;
步骤2:根据当前环境温度从预先获取的失配误差标定结果中线性插值,得到当前温度下的失配误差系数;
步骤3:通过硬件逻辑计算当前温度下的适配误差系数对应的修正系数,对当前温度下的测试结果实时修正;修正系数的计算通过完美重构法实现;
步骤3具体包括:
步骤3.1:启动修正系数计算逻辑,读入当前频率、当前温度对应的适配误差系数;
步骤3.2:计算完美重构修正所需的频点的角频率,对每个角频率下的完美重构矩阵的各元素进行求解计算,获得当前角频率下的修正滤波器的频域数值解;
步骤3.3:对频域数值解进行IFFT变换,得到对应的修正滤波器的时域数值解;
步骤3.4:加窗截断处理,得到修正滤波器的系数,可实现对当前温度下的测试结果实时修正。
2.根据权利要求1所述的一种交替并行采样系统变温条件下修正系数自动适配方法,其特征在于,步骤1具体包括:
步骤1.1:设计自动标定硬件逻辑,存储某一单频点信号下的ADC的采样数据,顺次对ADC各通道的采样数据进行四参数拟合,将拟合结果转化为相同形式;
步骤1.2:设定ADC某一通道为基准,基准的增益误差为1,相位误差为0,根据拟合结果得到ADC各通道相对于基准的失配误差标定结果,将失配误差标定结果写入存储,完成当前频点下的标定;
步骤1.3:改变输入信号频率,重复步骤1.1-1.2,完成设定的频带范围内所有频点下的ADC各通道的失配误差标定;
步骤1.4:改变温度,重复步骤1.3,完成在不同温度下ADC各通道的失配误差标定。
3.根据权利要求1所述的一种交替并行采样系统变温条件下修正系数自动适配方法,其特征在于,步骤3.2中,频点的数量N=2P。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110590735.4A CN113346901B (zh) | 2021-05-28 | 2021-05-28 | 一种交替并行采样系统变温条件下修正系数自动适配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110590735.4A CN113346901B (zh) | 2021-05-28 | 2021-05-28 | 一种交替并行采样系统变温条件下修正系数自动适配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113346901A CN113346901A (zh) | 2021-09-03 |
CN113346901B true CN113346901B (zh) | 2024-03-29 |
Family
ID=77471910
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110590735.4A Active CN113346901B (zh) | 2021-05-28 | 2021-05-28 | 一种交替并行采样系统变温条件下修正系数自动适配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113346901B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117526943B (zh) * | 2024-01-08 | 2024-03-29 | 成都能通科技股份有限公司 | 一种基于fpga的高速adc性能测试系统及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7782235B1 (en) * | 2007-04-30 | 2010-08-24 | V Corp Technologies, Inc. | Adaptive mismatch compensators and methods for mismatch compensation |
CN102857225A (zh) * | 2012-09-13 | 2013-01-02 | 电子科技大学 | 一种多通道高速并行交替采样系统的失配误差校准方法 |
CN108494402A (zh) * | 2018-03-14 | 2018-09-04 | 东南大学 | 一种基于正弦拟合的tiadc系统误差估计和补偿方法 |
CN110266311A (zh) * | 2019-05-29 | 2019-09-20 | 中国科学技术大学 | 一种tiadc系统失配误差校准方法、装置、设备及介质 |
CN112751564A (zh) * | 2019-10-31 | 2021-05-04 | 深圳市中兴微电子技术有限公司 | 采样时钟相位失配误差估计方法及装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7839323B2 (en) * | 2008-12-29 | 2010-11-23 | Intersil Americas, Inc. | Error estimation and correction in a two-channel time-interleaved analog-to-digital converter |
US8144040B2 (en) * | 2009-07-01 | 2012-03-27 | Intersil Americas, Inc. | Randomization of sample window in calibration of time-interleaved analog to digital converter |
US20170117914A1 (en) * | 2015-10-23 | 2017-04-27 | Industry-Academic Cooperation Foundation, Chosun University | Method and apparatus for providing digital background calibration for mismatches in m-channel time-interleved adcs (ti-adcs) |
-
2021
- 2021-05-28 CN CN202110590735.4A patent/CN113346901B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7782235B1 (en) * | 2007-04-30 | 2010-08-24 | V Corp Technologies, Inc. | Adaptive mismatch compensators and methods for mismatch compensation |
CN102857225A (zh) * | 2012-09-13 | 2013-01-02 | 电子科技大学 | 一种多通道高速并行交替采样系统的失配误差校准方法 |
CN108494402A (zh) * | 2018-03-14 | 2018-09-04 | 东南大学 | 一种基于正弦拟合的tiadc系统误差估计和补偿方法 |
CN110266311A (zh) * | 2019-05-29 | 2019-09-20 | 中国科学技术大学 | 一种tiadc系统失配误差校准方法、装置、设备及介质 |
CN112751564A (zh) * | 2019-10-31 | 2021-05-04 | 深圳市中兴微电子技术有限公司 | 采样时钟相位失配误差估计方法及装置 |
Non-Patent Citations (4)
Title |
---|
An 8-Gs/s 12-Bit TIADC System With Real-Time Broadband Mismatch Error Correction;Lei Zhao;IEEE Transactions on Nuclear Science;第65卷(第12期);2892 - 2900 * |
一种并行交替采样中时基非均匀信号自适应重构方法;潘卉青;田书林;叶;曾浩;;电子测量与仪器学报(第01期);全文 * |
基于循环自相关的TIADC通道失配校正;印茂伟;唐斌;向利;;微计算机信息(第16期);全文 * |
高速TIADC系统误差校准技术研究;安珊;杨波;王春阳;;仪器仪表用户(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113346901A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU2003253088B2 (en) | Improvements relating to time-interleaved samplers | |
CN108494402B (zh) | 一种基于正弦拟合的tiadc系统误差估计和补偿方法 | |
US7330140B2 (en) | Interleaved analog to digital converter with compensation for parameter mismatch among individual converters | |
JP4076553B2 (ja) | 校正装置及び線形補正器校正方法 | |
US8290032B2 (en) | Distortion identification apparatus, test system, recording medium and distortion identification method | |
CN102163972A (zh) | 对adc中非线性的修正 | |
CN113346901B (zh) | 一种交替并行采样系统变温条件下修正系数自动适配方法 | |
CN113258930B (zh) | 一种数字示波器及时间交织模数转换器的校正方法 | |
JPH1168567A (ja) | 集積回路用a/dコンバータの較正方法及び装置 | |
CN111044962B (zh) | 一种基于fft的带宽平坦度校准方法 | |
CN113114243B (zh) | 一种tiadc系统失配误差校正方法及系统 | |
CN115776299A (zh) | 一种低复杂度的tiadc时间失配误差校准方法 | |
CN113063978B (zh) | 一种数字示波器及采样时刻失配的校正方法 | |
CN110716167B (zh) | 任意波形发生器幅度扫频自适应分辨率校准方法及系统 | |
Kasher et al. | Memory-efficient SFDR-optimized post-correction of analog-to-digital converters via frequency-selective look-up tables | |
US7443326B2 (en) | Digital input signals constructor providing analog representation thereof | |
Lundin et al. | A criterion for optimizing bit-reduced post-correction of AD converters | |
CN112505413B (zh) | 一种时频分析方法和系统 | |
EP1678868A2 (en) | Analog signals sampler providing digital representation thereof | |
CN115236634A (zh) | 一种激光雷达中激光线性调频信号的校正方法及装置 | |
CN106788807B (zh) | 一种无线信道测试方法及装置 | |
Qian et al. | The fast time-interleaved ADC channel mismatch calibration method based on equivalent sampling and dynamic reference point | |
EP0945987B1 (en) | Process for the linearization of the transfer characteristic of an analog-to-digital converter and a respective circuit | |
CN113691258B (zh) | 一种基于遗传算法的adc动态性能参数计算方法 | |
CN113804981B (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 |