CN103166638B - 测试装置以及模拟至数字转换器的测试方法 - Google Patents

测试装置以及模拟至数字转换器的测试方法 Download PDF

Info

Publication number
CN103166638B
CN103166638B CN201210019061.3A CN201210019061A CN103166638B CN 103166638 B CN103166638 B CN 103166638B CN 201210019061 A CN201210019061 A CN 201210019061A CN 103166638 B CN103166638 B CN 103166638B
Authority
CN
China
Prior art keywords
analog
digital converter
value
test device
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
Application number
CN201210019061.3A
Other languages
English (en)
Other versions
CN103166638A (zh
Inventor
张宏盛
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.)
Silicon Motion Inc
Original Assignee
Silicon Motion Inc
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 Silicon Motion Inc filed Critical Silicon Motion Inc
Publication of CN103166638A publication Critical patent/CN103166638A/zh
Application granted granted Critical
Publication of CN103166638B publication Critical patent/CN103166638B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/10Calibration or testing
    • H03M1/1071Measuring or testing
    • H03M1/1095Measuring or testing for ac performance, i.e. dynamic testing
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/12Analogue/digital converters

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Tests Of Electronic Circuits (AREA)
  • Analogue/Digital Conversion (AREA)

Abstract

本发明提供一种测试装置。于一实施例中,该测试装置接收一模拟至数字转换器输出的多个位信号,包括多个频率计数器、一比较模块、以及一效能分析模块。所述频率计数器分别计算所述位信号的位值变换的次数,以得到多个转换次数。该比较模块分别比较所述转换次数与对应于所述转换次数的多个转换次数理想值以得到多个误差次数。该效能分析模块依据所述误差次数估计该模拟至数字转换器的一效能值。

Description

测试装置以及模拟至数字转换器的测试方法
【技术领域】
本发明系有关于模拟至数字转换器,特别是有关于模拟至数字转换器的测试。
【背景技术】
模拟至数字转换器是数字装置常用的组件,用以将模拟信号转换为数字信号。由于模拟至数字转换器通常用以将模拟输入信号转换为数字信号以作为数字处理器的输入,若模拟至数字转换器的效能不佳,将会于数字处理器的输入信号引起错误,导致数字处理器的输出亦发生错误。因此,模拟至数字转换器的效能是十分重要的,而必须以严格的测试以淘汰具较差效能的模拟至数字转换器。
一般而言,模拟至数字转换器的测试系于晶圆产出阶段或芯片封装后进行。一般系以测试装置对模拟至数字转换器进行测试工作。图1为模拟至数字转换器190的测试装置100的区块图。模拟至数字转换器190接受一模拟输入信号后将模拟输入信号转换为多个输出比特流bit 0、bit 1、…、bit 9。测试装置100接收输出比特流后依据输出比特流评估模拟至数字转换器190的效能。于一实施例中,测试装置100包括十进制转换器110、长条图绘制器120、以及效能分析模块130。
十进制转换器110将模拟至数字转换器190的二进制输出比特流bit 0、bit 1、…、bit 9转换为一连串的十进制值。长条图绘制器120再把一连串的十进制值绘制成长条图(histogram)。效能分析模块130再依据长条图分析积分非线性值(integral nonlinearity,INL)以及微分非线性值(differentialnonlinearity,DNL),以决定模拟至数字转换器190的效能。图2A显示模拟至数字转换器190的模拟输入信号,而图2B显示长条图绘制器120依据十进制值所绘制的长条图。
然而,图1的测试装置100具有较高的成本与较复杂的电路设计。首先,为了储存十进制转换器110所产生的多个十进制值,长条图绘制器120必须具有大容量的内存,而大容量的内存会提高测试装置100的硬件成本。其次,长条图绘制器120与效能分析模块130均需要高效能的复杂计算,因而需要复杂的电路,亦提高了测试装置100的硬件成本。为了减少测试装置的硬件成本,需要用更简单的方法以计算模拟至数字转换器的效能。因此,需要一种模拟至数字转换器的新型态的测试装置。
【发明内容】
有鉴于此,本发明的目的在于提供一种测试装置,以解决现有技术存在的问题。于一实施例中,该测试装置接收一模拟至数字转换器输出的多个位信号,包括多个频率计数器、一比较模块、以及一效能分析模块。所述频率计数器分别计算所述位信号的位值变换的次数,以得到多个转换次数。该比较模块分别比较所述转换次数与对应于所述转换次数的多个转换次数理想值以得到多个误差次数。该效能分析模块依据所述误差次数估计该模拟至数字转换器之一效能值。
本发明更提供一种模拟至数字转换器的测试方法。首先,接收一模拟至数字转换器输出的多个位信号。接着,以多个频率计数器(frequency counter)分别计算所述位信号的位值变换的次数,以得到多个转换次数。接着,以一比较模块分别比较所述转换次数与对应于所述转换次数的多个转换次数理想值以得到多个误差次数。接着,以一效能分析模块依据所述误差次数估计该模拟至数字转换器之一效能值。
为了让本发明的上述和其它目的、特征、和优点能更明显易懂,下文特举数较佳实施例,并配合所附图示,作详细说明如下:
【附图说明】
图1为模拟至数字转换器的测试装置的区块图;
图2A显示模拟至数字转换器的模拟输入信号;
图2B显示图1的测试装置的长条图绘制器依据十进制值所绘制的长条图;
图3依据本发明的模拟至数字转换器的测试装置的区块图;
图4为依据本发明测试模拟至数字转换器的方法的流程图;
图5为依据本发明的模拟至数字转换器输出的多个位信号与频率计数器产生的转换次数间的关系之一实施例的示意图;
图6A为理想的模拟至数字转换器所产生的的多个位信号的理想值的示意图;
图6B为实际的模拟至数字转换器所产生的的多个位信号的实际值的示意图;
图7为依据本发明的各位信号的转换次数的理想值与实际值的关系之一实施例的示意图。
【主要组件符号说明】
(图1)
190~模拟至数字转换器;
100~测试装置;
110~十进制转换器;
120~长条图绘制器;
130~效能分析模块;
(图3)
390~模拟至数字转换器;
300~测试装置;
310,317,318,319~低通滤波器;
320,327,328,329~频率计数器;
330~比较模块;
340~效能分析模块。
【具体实施方式】
图3为依据本发明的模拟至数字转换器390的测试装置300的区块图。模拟至数字转换器390接收一模拟输入信号,并对模拟输入信号进行模拟至数字转换而产生多个位信号S0、…、S7、S8、S9。位信号S0、…、S7、S8、S9为模拟输入信号经转换所得的位值bit 0、…、bit 7、bit 8、bit 9。于一实施例中,模拟至数字转换器390所接收的模拟输入信号为三角波输入信号(trianglewave)或斜面波输入信号(ramp wave)。
测试装置300接收模拟至数字转换器390输出的多个位信号S0、…、S7、S8、S9,并依据所述位信号S0、…、S7、S8、S9估计反映模拟至数字转换器390的效能之一效能值。于一实施例中,测试装置300包括多个低通滤波器310~319、频率计数器320~329、比较模块330、以及效能分析模块340。低通滤波器310~319过滤模拟至数字转换器390输出的多个位信号S0、…、S7、S8、S9,以产生多个过滤位信号S0’、…、S7’、S8’、S9’。于一实施例中,低通滤波器310~319分别将位信号S0、…、S7、S8、S9的每连续三个位值平均,以得到过滤位信号S0’、…、S7’、S8’、S9’。于另一实施例中,低通滤波器310~319分别将位信号S0、…、S7、S8、S9的每连续四个位值平均,以得到过滤位信号S0’、…、S7’、S8’、S9’。
频率计数器320~329分别计算所述过滤位信号S0’、…、S7’、S8’、S9’的位值变换的次数,以得到多个转换次数F0、…、F7、F8、F9。于一实施例中,频率计数器320~329分别计算过滤位信号S0’、…、S7’、S8’、S9’由位值1变换至位值0的第一总次数f1、…、f7、f8、f9,再计算位信号S0’、…、S7’、S8’、S9’由位值0变换至位值1的第二总次数f1’、…、f7’、f8’、f’9,并分别将所述第一总次数与所述第二总次数相加(f1+f1’)、…、(f7+f7’)、(f8+f8’)、(f9+f9’)以得到转换次数F0、…、F7、F8、F9
比较模块330接收频率计数器320~329产生的转换次数F0、…、F7、F8、F9。比较模块330接着分别比较转换次数F0、…、F7、F8、F9与分别对应于所述转换次数的多个转换次数理想值F0’、…、F7’、F8’、F9’以得到多个误差次数E0、…、E7、E8、E9。于一实施例中,比较模块330将转换次数F0、…、F7、F8、F9分别与对应的转换次数理想值F0’、…、F7’、F8’、F9’相减以得到误差次数E0、…、E7、E8、E9。接着,效能分析模块340依据所述误差次数E0、…、E7、E8、E9估计模拟至数字转换器390之一效能值。于一实施例中,效能分析模块340依据多个权重加总所述误差次数E0、…、E7、E8、E9,以得到模拟至数字转换器390的效能估计值。当误差次数E0、…、E7、E8、E9愈高,模拟至数字转换器390的效能估计值愈高,表示模拟至数字转换器390的效能愈低。
由于测试装置300直接以多个频率计数器320~329计算模拟至数字转换器390所输出的位信号的转换次数,再依据转换次数估计模拟至数字转换器390的效能,因此测试装置300不似现有技术的测试装置100般需具备大容量的内存以储存模拟至数字转换器390所输出的位信号值。因此,测试装置300的硬件成本较现有技术的测试装置100减少许多。其次,测试装置300的频率计数器320~329、比较模块330、效能分析模块340的计算都相当的简单,因此具有简单的电路设计,再次减少了测试装置300的硬件成本。因此,本发明的测试装置300可以藉较低的硬件成本实现模拟至数字转换器390的效能评估。
图4为依据本发明测试模拟至数字转换器的方法400的流程图。首先,输入一三角波至模拟至数字转换器390(步骤402)。接着,测试装置300接收模拟至数字转换器390输出的多个位信号S0、…、S9(步骤404)。接着,测试装置300以低通滤波器310、…、319分别过滤所述位信号S0、…、S9以得到多个过滤位信号S0’、…、S9’(步骤406)。接着,测试装置300以频率计数器320、…、329分别计算所述过滤位信号S0’、…、S9’的位值转换的次数以得到多个转换次数F0、…、F9(步骤408)。接着,比较模块330将所述转换次数F0、…、F9与多个理想转换次数比较以得到多个误差次数E0、…、E9(步骤410)。最后,效能分析模块340依据所述误差次数E0、…、E9计算模拟至数字转换器390之一效能值后输出(步骤412)。于一实施例中,若最低的三个有效位以外的任何其它所述误差次数为正,则效能分析模块340判定模拟至数字转换器390的效能值为不及格。
图5为依据本发明的模拟至数字转换器输出的多个位信号与频率计数器产生的转换次数间的关系之一实施例的示意图。假设模拟至数字转换器390共输出10个位信号S0、…、S9,分别对应图5中的位0、…、位9,其中位0为最低有效位(least significant bit,LSB),而位9为最高有效位(mostsignificant bit,MSB)。假设测试装置300自模拟至数字转换器390依序接收到分别对应于16个时点的样本。例如对应于时序10的(位9、…、位0)的样本为(0,0,0,0,0,0,1,0,0,1),而对应于时序16的(位9、…、位0)的样本为(0,0,0,0,0,0,1,1,1,1)。图1的现有技术的测试装置100会将对应于时序10的位信号样本转换为十进制值9,并将对应于时序16的位信号样本转换为十进制值15,再依据十进制值计算模拟至数字转换器的效能估计值,因此需要储存每个时点的十进制值,而增加了内存的硬件成本。
相反的,依据本发明的测试装置300仅运用频率计数器320、…、329分别计数字元9、…、位0的位值由0变动为1或由1变动为0的转换次数,而分别得到对应于位9、…、位2、位1、位0的转换次数0、…、3、7、15,因此不需储存每个时点的位信号的值,而省去了内存的硬件成本。举例来说,位1的值于时序3时由0转为1、于时序5时由1转为0、于时序7时由0转为1、于时序9时由1转为0、于时序11时由0转为1、于时序13时由1转为0、于时序15时由0转为1,因此共计有转换次数7。位2的值于时序5时由0转为1、于时序9时由1转为0、于时序13时由0转为1,因此共计有转换次数3。
理想的模拟至数字转换器不会于转换过程产生噪声(noise),藉此产生的位信号为理想值。图6A为理想的模拟至数字转换器所产生的的多个位信号S0、…、S9的理想值的示意图。图6A中显示了时序1~时序14的各位信号的位值。最低有效位0的值每隔一时序即变换一次位值(由0变为1或由1变为0),而最高有效位9的值于时序1~时序14均未改变。各时序所累加产生的十进制值逐次累加一。由于实际的模拟至数字转换器会于模拟至数字转换过程产生噪声,因此产生的位信号为带有误差的实际值。图6B为实际的模拟至数字转换器所产生的的多个位信号S0、…、S9的实际值的示意图。由于噪声所产生的误差,使得最低有效位0的值并非每隔一时序即变换一次位值,而各时序所累加产生的十进制值亦未逐次累加一。
图6B中的实际的模拟至数字转换器所产生的的位信号S0~S9因噪声所带有的误差会导致图3的频率计数器320~329计数得到错误的转换次数F0~F9。因此,必须以低通滤波器310~319对位信号S0~S9进行过滤以产生过滤位信号S0’~S9’。图7为依据本发明的各位信号的转换次数的理想值与实际值的关系之一实施例的示意图。理想状态下的模拟至数字转换器390输出的多个位信号S0、S1、S2、…、S9的转换次数的理想值分别为4091、2044、1020、…、4,如表中第1列所示。经过低通滤波器将位信号的每三个位值平均一次所得的过滤位信号S0’、S1’、S2’、…、S9’的转换次数的理想值分别为4085、2044、1020、…、4,如表中第2列所示。经过低通滤波器将位信号的每四个位值平均一次所得的过滤位信号S0’、S1’、S2’、…、S9’的转换次数的理想值分别为5、4、1020、…、4,如表中第3列所示。实际的模拟至数字转换器390输出的多个位信号S0、S1、S2、…、S9的转换次数的实际值分别为2393、2471、1808、…、12,如表中第4列所示。经过低通滤波器将位信号的每四个位值平均一次所得的过滤位信号S0’、S1’、S2’、…、S9’的转换次数的实际值分别为583、565、914、…、4,如表中第5列所示。由图5中可见,由于低通滤波器310、…、319会对位信号S0、S1、S2、…、S9进行过滤及平均,因此对应于最低有效位S0及次低有效位S1的过滤后位信号S0’、S1’的转换次数的实际测试值与理想值相差很大,而无法用以作为模拟至数字转换器390的效能估计值的计算依据。因此,效能分析模块340仅依据未对应最低有效位S0及次低有效位S1的误差次数E9、…、E3以估计模拟至数字转换器390的效能值。详细言的,图3的比较模块330与效能分析模块340用以判断模拟至数字转换器390的效能的判断基准如下。首先由比较模块330比较转换次数实际值F0~F9与相对应的转换次数理想值F0’~F9’以得到误差次数E0~E9。接着由效能分析模块340依据较高有效位的误差次数计算模拟至数字转换器390的效能值。以图7中的LPF4的转换次数的实际值与理想值比较,可以发现仅最低有效位0、1、2的实际值与理想值有误差。因此,若除了对应于最低有效位信号S0、S1、S2的误差次数E0、E1、E2外,效能分析模块340发现其它的误差次数E3~E9任一为正,则效能分析模块340判定模拟至数字转换器390的效能为不合格。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明,任何熟习此项技术者,在不脱离本发明的精神和范围内,当可作些许的更动与润饰,因此本发明的保护范围当视后附的申请专利范围所界定者为准。

Claims (9)

1.一种测试装置,接收一模拟至数字转换器输出的多个位信号,包括:
多个频率计数器(frequency counter),分别计算所述位信号的位值变换的次数,以得到多个转换次数;
一比较模块,分别比较所述转换次数与对应于所述转换次数的多个转换次数理想值以得到多个误差次数;以及
一效能分析模块,依据所述误差次数估计该模拟至数字转换器之一效能值。
2.根据权利要求1所述的测试装置,其特征在于,该模拟至数字转换器接收一三角波输入信号(triangle wave)或一斜面波输入信号(ramp wave),并对该三角波输入信号或该斜面波输入信号进行模拟至数字转换,以产生所述位信号。
3.根据权利要求1所述的测试装置,其特征在于,该模拟至数字转换器所输出的所述位信号于同一时序的多个位值的组合为该模拟至数字转换器之一输入值的二进制值。
4.根据权利要求1所述的测试装置,其特征在于,所述频率计数器分别计算所述位信号由位值1变换至位值0的多个第一总次数以及所述位信号由位值0变换至位值1的多个第二总次数,并分别将所述第一总次数与所述第二总次数相加以得到所述转换次数。
5.根据权利要求1所述的测试装置,其特征在于,该测试装置更包括多个低通滤波器,分别耦接于该模拟至数字转换器与所述频率计数器其中之一的间,分别过滤所述位信号以得到多个过滤的位信号,再将所述过滤的位信号送至所述频率计数器作为所述频率计数器的输入。
6.根据权利要求5所述的测试装置,其特征在于,所述低通滤波器分别将所述位信号的每连续三个或四个位值平均,以得到所述过滤的位信号。
7.根据权利要求1所述的测试装置,其特征在于,若最低的三个有效位以外的任何其它所述误差次数为正,该效能分析模块判定该模拟至数字转换器的该效能值为不及格。
8.根据权利要求1所述的测试装置,其特征在于,该效能分析模块依据多个权重加总所述误差次数以得到该效能值;其中当所述误差次数愈高而使该效能值愈高,表示该模拟至数字转换器的效能愈低。
9.根据权利要求1所述的测试装置,其特征在于,该效能分析模块仅依据未对应最低有效位及次低有效位的所述误差次数以估计该效能值。
CN201210019061.3A 2011-12-15 2012-01-20 测试装置以及模拟至数字转换器的测试方法 Active CN103166638B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
TW100146406A TWI483555B (zh) 2011-12-15 2011-12-15 測試裝置以及類比至數位轉換器之測試方法
TW100146406 2011-12-15

Publications (2)

Publication Number Publication Date
CN103166638A CN103166638A (zh) 2013-06-19
CN103166638B true CN103166638B (zh) 2016-12-07

Family

ID=48589410

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210019061.3A Active CN103166638B (zh) 2011-12-15 2012-01-20 测试装置以及模拟至数字转换器的测试方法

Country Status (3)

Country Link
US (1) US8648740B2 (zh)
CN (1) CN103166638B (zh)
TW (1) TWI483555B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4352160A (en) * 1980-01-21 1982-09-28 The United States Of America As Represented By The Secretary Of The Air Force Statistical method of measuring the differential linearity of an analog/digital converter using a pseudo-random triangle wave stimulus
US7482958B2 (en) * 2006-10-27 2009-01-27 Oki Semiconductor Co., Ltd. OFDM receiver

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3683369A (en) * 1971-06-25 1972-08-08 Honeywell Inc Analog to digital converter
FR2636740B1 (fr) * 1988-09-16 1990-12-28 Sgs Thomson Microelectronics Detecteur d'enveloppe logarithmique de signal analogique
US5909186A (en) * 1997-07-01 1999-06-01 Vlsi Technology Gmbh Methods and apparatus for testing analog-to-digital and digital-to-analog device using digital testers
CN1163057C (zh) * 1999-03-31 2004-08-18 富士通将军股份有限公司 画质补正电路
JP3810318B2 (ja) * 2001-12-28 2006-08-16 有限会社ニューロソリューション アナログデジタル変換装置
DE50303186D1 (de) * 2003-02-05 2006-06-08 Alcatel Sa Verfahren zur Digital-Analog-Konvertierung und Digital-Analog-Konverter
JP4266350B2 (ja) * 2004-02-12 2009-05-20 株式会社ルネサステクノロジ テスト回路
JP4364175B2 (ja) * 2005-08-01 2009-11-11 株式会社東芝 乗算器及びこれを用いる無線通信装置
JP5270173B2 (ja) * 2006-02-01 2013-08-21 富士通株式会社 半導体装置及びノイズ計測方法
JP2008131583A (ja) * 2006-11-24 2008-06-05 Toshiba Corp 並列型a/d変換器
US7564386B2 (en) * 2007-11-16 2009-07-21 Agilent Technologies, Inc. Pre-processing data samples from parallelized data converters

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4352160A (en) * 1980-01-21 1982-09-28 The United States Of America As Represented By The Secretary Of The Air Force Statistical method of measuring the differential linearity of an analog/digital converter using a pseudo-random triangle wave stimulus
US7482958B2 (en) * 2006-10-27 2009-01-27 Oki Semiconductor Co., Ltd. OFDM receiver

Also Published As

Publication number Publication date
US8648740B2 (en) 2014-02-11
CN103166638A (zh) 2013-06-19
TW201325097A (zh) 2013-06-16
US20130154861A1 (en) 2013-06-20
TWI483555B (zh) 2015-05-01

Similar Documents

Publication Publication Date Title
CN110543860B (zh) 基于tjm迁移学习的机械故障诊断方法及系统
CN101615010B (zh) 基于fpga的多路数据采集系统
CN101226133B (zh) 一种血细胞脉冲信号的分类识别方法
CN107064628A (zh) 高精度频率测量系统及方法
CN106556342B (zh) 一种基于fpga的光栅细分装置及方法
CN101710304A (zh) 一种软件过程实施质量评测方法
CN110673532B (zh) 一种脉冲信号处理和采集装置
CN104657169A (zh) 稼动率计算方法及其系统和嵌入式装置
CN107315991A (zh) 一种基于小波阈值去噪的ifra频响曲线去噪方法
CN103529379B (zh) 用于高精度adc测试中低分辨率信号源的选取方法
CN101424919A (zh) 半导体制造系统的重调度决策系统
CN103475369A (zh) 基于信号源误差一次性校准识别的高精度adc测试方法
CN111626360A (zh) 用于检测锅炉故障类型的方法、装置、设备和存储介质
JP2024524666A (ja) 衝撃振動センサー信号のピックアップ方法
CN103166638B (zh) 测试装置以及模拟至数字转换器的测试方法
CN105182219A (zh) 一种基于汉明纠错码支持向量机的功率变换器故障分类方法
CN113189513B (zh) 一种基于纹波的冗余电源均流状态识别方法
CN101840451B (zh) 一种集成电路工艺参数模型的优化方法
CN106802293B (zh) 波形峰值检测方法及装置
WO2024139035A1 (zh) 模数转换器的校正方法
CN103095249B (zh) 一种中值滤波电路及方法
CN105718240B (zh) 一种基于总和增量调制的任意路输入信号的比特流加法器
CN111610428A (zh) 一种基于响应混叠性度量小波包分解算法的参数优化方法
CN115130748A (zh) 一种建筑电气系统故障诊断方法
CN101672664B (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
C14 Grant of patent or utility model
GR01 Patent grant