CN110987846A - 一种基于iPLS-PA算法的硝酸盐浓度预测方法 - Google Patents

一种基于iPLS-PA算法的硝酸盐浓度预测方法 Download PDF

Info

Publication number
CN110987846A
CN110987846A CN201911182738.3A CN201911182738A CN110987846A CN 110987846 A CN110987846 A CN 110987846A CN 201911182738 A CN201911182738 A CN 201911182738A CN 110987846 A CN110987846 A CN 110987846A
Authority
CN
China
Prior art keywords
prediction
interval
nitrate
modeling
nitrate concentration
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
CN201911182738.3A
Other languages
English (en)
Other versions
CN110987846B (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.)
XiAn Institute of Optics and Precision Mechanics of CAS
Original Assignee
XiAn Institute of Optics and Precision Mechanics of CAS
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 XiAn Institute of Optics and Precision Mechanics of CAS filed Critical XiAn Institute of Optics and Precision Mechanics of CAS
Priority to CN201911182738.3A priority Critical patent/CN110987846B/zh
Publication of CN110987846A publication Critical patent/CN110987846A/zh
Application granted granted Critical
Publication of CN110987846B publication Critical patent/CN110987846B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/33Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明属于硝酸盐浓度预测方法,具体涉及一种基于iPLS‑PA算法的硝酸盐浓度预测方法。解决传统的硝酸盐测定方法存在的价格昂贵、操作复杂且分析时间长、需要消耗试剂、存在对水体的二次污染等问题。该方法利用每一个子区间的峰面积与对应的硝酸盐浓度进行交叉验证建模,选择出最优的峰面积计算区间以获得更高的硝酸盐浓度预测精度。

Description

一种基于iPLS-PA算法的硝酸盐浓度预测方法
技术领域
本发明属于硝酸盐浓度预测方法,具体涉及一种基于iPLS-PA算法的硝酸盐浓度预测方法。
背景技术
水体中硝酸盐浓度过高会引起水体富营养化、危害人体健康等问题。传统的硝酸盐测定方法包括离子色谱法、镉柱还原法、离子电极法等,但普遍存在一些缺点,如大多数方法价格昂贵、操作复杂且分析时间长、需要消耗试剂、存在对水体的二次污染等。
发明内容
本发明的目的是提供一种基于iPLS-PA算法的硝酸盐浓度预测方法,解决传统的硝酸盐测定方法存在的价格昂贵、操作复杂且分析时间长、需要消耗试剂、存在对水体的二次污染等问题。
紫外可见光谱技术具有测量方便快捷、无需引入其他试剂等优点,因此可以用于硝酸盐浓度的测量。该方法利用测量物质吸收紫外可见光谱辐射的原理,代入预测模型对新样本的浓度进行预测,灵敏度高、准确度好。
在预测模型的建模分析时,若采用全谱建立多元校正模型,通常存在计算量大、模型不够稳健、且模型的预测精度往往也不能达到最优等问题,因此,通过特定方法筛选特征波长或波长区间有可能得到效果更好的浓度预测模型。
针对建模波长区间选择的问题,Norgaard等在2000年提出了间隔偏最小二乘(Interval PLS,iPLS)方法。iPLS是基于偏最小二乘法(Partial least squares,PLS)内部形成的交互式扩展,可在全谱区域的等距子区间内开发局部PLS模型。
它的主要作用是提供不同光谱区间中相关信息的整体情况,从而关注重要的光谱区域并消除其他区域的干扰。iPLS方法提供的信息突出了PLS算法对噪声变量的敏感性。iPLS模型是在等宽度的光谱子区间上进行计算的,并且比较了这些局部模型和全局(全光谱)模型的预测性能。
比较主要基于验证参数交叉验证的均方根误差(RMSECV),同时也评估其他参数,如决定系数(R2),斜率和偏移等,以保证对模型全面的描述。
该算法的基本思路是:将测得吸光度谱线区间均匀的划分为n个子区间,在每个子区间内分别使用PLS进行交叉验证建模,选取最优模型对应的子区间为最优建模区间。
峰面积(Peak area,PA)指测得光谱曲线中吸收峰的面积,峰面积与样品的浓度成正比。由于吸收光谱的峰面积受样品和仪器的影响比吸收峰峰高小,因此通常使用峰面积进行定量计算的结果比使用峰高更准确。
基于以上分析,本发明提出了一种基于iPLS-PA算法的硝酸盐浓度预测方法,可以快速、准确、稳定的预测水中硝酸盐浓度,包括以下步骤:
步骤一、获得紫外可见光谱曲线;
步骤二、预处理紫外可见光谱曲线;
对全光谱曲线进行预处理,去除由仪器产生的噪声及液体中固体沉淀物质颗粒对光谱曲线产生的影响;
步骤三、划分子区间;
将预处理后的全光谱曲线在全光谱区间等间距划分为I个子区间;
步骤四、计算每个子区间对应的峰面积;
计算不同浓度下每一个子区间对应的峰面积记为Si,j,其中,i=1,2,3,…,I,j=1,2,3,…,J,J表示水样样本数;
步骤五、针对每个子区间分别建模;
在每个子区间内分别使用峰面积与对应的硝酸盐浓度进行线性回归交叉验证建模;
步骤六、选出预测性能高的若干个子区间;
计算并比较每个子区间模型的预测性能,选择预测性能大与预设值的若干个子区间,将选择的若干个子区间进行合并;
步骤七、确定合并后区间的预测精度;
在该合并后的区间内使用峰面积与对应的硝酸盐浓度进行线性回归交叉验证建模,确定预测精度;
步骤八、获得最优建模区间;
判断合并后区间的预测精度是否达到要求或达到循环次数,若是,进入步骤九;若否,以合并后子区间为中心,单向、双向扩充、或消减波长变量,作为新的光谱区间;将该光谱区间等间距划分为I个子区间;重复步骤四至步骤六,直至预测精度达到要求或达到循环次数,选出最优建模区间;
步骤九、获得预测模型;
采用线性回归的方法计算所有样本在最优建模区间内峰面积与对应硝酸盐浓度间的回归系数,得到的回归系数即可作为硝酸盐的浓度预测模型;
步骤十、获得硝酸盐浓度;
将新样本的最优建模区间内的峰面积带入预测模型,输出预测结果,获得新样本的硝酸盐浓度。
进一步地,步骤五中采用蒙特卡洛交叉验证方法进行线性回归交叉验证建模。
进一步地,步骤六中预测性能通过均方根误差表征。
进一步地,步骤七中预测精度通过决定系数R2,斜率和偏移表征。
进一步地,步骤四中采用定积分的元素法计算峰面积:
Figure BDA0002291704480000031
其中f(x)为峰面积处的光谱曲线。
本发明的有益效果是:
1、本发明利用每一个子区间的峰面积与对应的硝酸盐浓度进行交叉验证建模,选择出最优的峰面积计算区间以获得更高的硝酸盐浓度预测精度。
首先,使用吸收光谱的峰面积的方法相比于使用吸收峰峰高的方法,其受到样品自身和仪器的影响较小。
其次,借鉴并改进了iPLS方法中区间优化选择的思路,选出了更适合峰面积计算的区间,减小了水中其他物质对硝酸盐浓度预测的干扰,从而进一步提高了硝酸盐浓度预测的精度。
2、紫外可见光谱技术测量方便快捷、无需引入其他试剂,因此不会对水体产生二次污染。
3、本发明iPLS-PA算法既能保证对海水样本中硝酸盐浓度预测准确性,同时能有效减少参与运算的数据量。在使用iPLS-PA方法后,针对小麦岛海水样本这类型仅有硝酸盐浓度差异的样本,其预测精度相比单波长法得到了提升,与PLS方法的预测精度相近,但参与建模的波长点数仅为60,是PLS建模时的1/9。
4、本发明iPLS-PA算法既能提高受浊度色度干扰的样本中硝酸盐浓度预测准确性,同时能有效减少参与运算的数据量。在使用iPLS-PA方法后,针对不同色度浊度的硝酸盐混合溶液样本,其预测精度高于单波长法及PLS建模方法,同时仅使用了全部450个变量中的15个。
附图说明
图1为本发明方法流程图;
图2为本发明采用的定积分元素法计算峰面积的示意图;
图3为实施例中小麦岛海水样本吸光度光谱图;
图4为实施例中混合溶液样本的光谱曲线。
具体实施方式
以下结合附图及具体实施例对本发明做进一步地描述。
常规的基于峰面积的方法是通过对不同浓度水样进行测量计算,得出测量的吸光度曲线峰面积与硝酸盐浓度的回归曲线,即可根据相应关系计算出未知水样的硝酸盐浓度。峰面积的计算方式对预测的精度有着直接影响。
本发明提出一种可以快速、准确、稳定的预测水中硝酸盐浓度的方法,基于间隔偏最小二乘法结合峰面积(iPLS-PA)的建模预测方法,用以优化计算峰面积时所选择的波长区间,以获得最适用于硝酸盐浓度建模的峰面积,进而使得硝酸盐浓度的预测精度提高。
具体步骤如图1所示:
(1)获取紫外可见光谱曲线,对全光谱曲线进行预处理,去除由仪器产生的噪声及液体中固体沉淀物质颗粒对光谱曲线产生的影响;
(2)划分子区间,并计算每个子区间的面积;
将全光谱区间等间距划分为I个子区间,计算不同浓度下每一个子区间内的峰面积记为Si,j,其中,i=1,2,3,…,I,j=1,2,3,…,J,J表示水样样本数;
(3)每个子区间分别回归建模;
将每一个子区间内不同浓度样本的峰面积并与对应的硝酸盐浓度yj进行交叉验证回归建模;在这里使用了重复随机子采样验证也称作蒙特卡洛交叉验证,将数据集随机分成训练集和测试集进行若干次建模验证,平均所有验证的结果。
(4)计算出每个子区间的RMSECV,以此确定预测效果最好的若干个波长子区间;将该若干个波长子区间进行合并,并在该合并后的区间内使用峰面积与对应的硝酸盐浓度进行线性回归交叉验证建模,确定预测精度;此处预测精度是通过R2及RMSECV表现的。
(5)判断预测精度是否达到要求,或是否达到循环次数,若是,进入(6),若否,再以选出的这些区间为中心单向或双向扩充(或消减)波长变量,将扩充后的建模区间再划分为若干个子区间进行第(2)步子区间的筛选,重复此过程,组合预测精度高的建模波长子区间,从而得到最佳的波长区间。
(6)获得预测模型;
采用线性回归的方法计算所有样本在最优建模区间内峰面积与对应硝酸盐浓度间的回归系数,得到的回归系数即可作为硝酸盐的浓度预测模型;
(7)、获得硝酸盐浓度;
将新样本的最优建模区间内的峰面积带入预测模型,输出预测结果,获得新样本的硝酸盐浓度。
计算峰面积时,峰高乘半高峰宽法和三角形法一方面受峰的对称性限制,另一方面与峰高乘峰宽的平均法和分割加和法类似,都属于近似估算的方法,其计算误差较大,不适用于高精度的建模。
蒙特卡洛法需要大量的复杂抽样,模拟时要求样本数量足够大,才能使随机点的分布与真实面积的比例接近。
因此本实施例采用了定积分的元素法,该方法简单直接,且不受波形对称性的影响。
定积分的元素法假设所求峰面积处的光谱曲线可以用y=f(x)表示,x∈[l1,l2],则峰面积可以表示为以[l1,l2]为积分区间的定积分:
Figure BDA0002291704480000061
如图2,峰面积可看作一个曲边梯形,其面积的微分为:
dA(x)=f(x)dx …(2),
以dx为宽的曲边梯形面积为:
Figure BDA0002291704480000062
在点x处,高为f(x)、宽为dx的矩形面积为:f(x)dx。则有:
ΔA≈f(x)dx …(4)
f(x)dx便是曲边梯形的面积元素,[l1,l2]区间内的峰面积则是以面积元素f(x)dx为底的被积表达式,即[l1,l2]区间内的定积分。
以下使用本发明预测方法分别在使用小麦岛海水配制的不同浓度硝酸盐溶液样本(光谱曲线如图3),以及加入不同浊度色度的硝酸盐溶液样本(光谱曲线如图4)上进行实验。同时对比了使用单波长法、PLS的预测结果,如表1、表2。
表1不同方法预测小麦岛海水样本中硝酸盐浓度结果对比
Figure BDA0002291704480000063
表2不同方法预测混合溶液样本中硝酸盐浓度结果对比
Figure BDA0002291704480000064
Figure BDA0002291704480000071
可以看出,在使用iPLS-PA方法后,针对小麦岛海水样本这类型仅有硝酸盐浓度差异的样本,其预测精度相比单波长法得到了提升,与PLS方法的预测精度相近,但参与建模的波长点数仅为60,是PLS建模时的1/9。
在使用iPLS-PA方法后,针对不同色度浊度的硝酸盐混合溶液样本,其预测精度高于单波长法及PLS建模方法,同时仅使用了全部450个变量中的15个。

Claims (5)

1.一种基于iPLS-PA算法的硝酸盐浓度预测方法,其特征在于,包括以下步骤:
步骤一、获得紫外可见光谱曲线;
步骤二、预处理紫外可见光谱曲线;
对全光谱曲线进行预处理,去除由仪器产生的噪声及液体中固体沉淀物质颗粒对光谱曲线产生的影响;
步骤三、划分子区间;
将预处理后的全光谱曲线在全光谱区间等间距划分为I个子区间;
步骤四、计算每个子区间对应的峰面积;
计算不同浓度下每一个子区间对应的峰面积记为Si,j,其中,i=1,2,3,…,I,j=1,2,3,…,J,J表示水样样本数;
步骤五、针对每个子区间分别建模;
在每个子区间内分别使用峰面积与对应的硝酸盐浓度进行线性回归交叉验证建模;
步骤六、选出预测性能高的若干个子区间;
计算并比较每个子区间模型的预测性能,选择预测性能大与预设值的若干个子区间,将选择的若干个子区间进行合并;
步骤七、确定合并后区间的预测精度;
在该合并后的区间内使用峰面积与对应的硝酸盐浓度进行线性回归交叉验证建模,确定预测精度;
步骤八、获得最优建模区间;
判断合并后区间的预测精度是否达到要求或达到循环次数,若是,进入步骤九;若否,以合并后子区间为中心,单向、双向扩充、或消减波长变量,作为新的光谱区间;将该光谱区间等间距划分为I个子区间;重复步骤四至步骤六,直至预测精度达到要求或达到循环次数,选出最优建模区间;
步骤九、获得预测模型;
采用线性回归的方法计算所有样本在最优建模区间内峰面积与对应硝酸盐浓度间的回归系数,得到的回归系数即可作为硝酸盐的浓度预测模型;
步骤十、获得硝酸盐浓度;
将新样本的最优建模区间内的峰面积带入预测模型,输出预测结果,获得新样本的硝酸盐浓度。
2.根据权利要求1所述的基于iPLS-PA算法的硝酸盐浓度预测方法,其特征在于:步骤五中采用蒙特卡洛交叉验证方法进行线性回归交叉验证建模。
3.根据权利要求2所述的基于iPLS-PA算法的硝酸盐浓度预测方法,其特征在于:步骤六中预测性能通过均方根误差表征。
4.根据权利要求3所述的基于iPLS-PA算法的硝酸盐浓度预测方法,其特征在于:步骤七中预测精度通过决定系数R2,斜率和偏移表征。
5.根据权利要求4所述的基于iPLS-PA算法的硝酸盐浓度预测方法,其特征在于:步骤四中采用定积分的元素法计算峰面积:
Figure FDA0002291704470000021
其中f(x)为峰面积处的光谱曲线。
CN201911182738.3A 2019-11-27 2019-11-27 一种基于iPLS-PA算法的硝酸盐浓度预测方法 Active CN110987846B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911182738.3A CN110987846B (zh) 2019-11-27 2019-11-27 一种基于iPLS-PA算法的硝酸盐浓度预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911182738.3A CN110987846B (zh) 2019-11-27 2019-11-27 一种基于iPLS-PA算法的硝酸盐浓度预测方法

Publications (2)

Publication Number Publication Date
CN110987846A true CN110987846A (zh) 2020-04-10
CN110987846B CN110987846B (zh) 2021-04-20

Family

ID=70087285

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911182738.3A Active CN110987846B (zh) 2019-11-27 2019-11-27 一种基于iPLS-PA算法的硝酸盐浓度预测方法

Country Status (1)

Country Link
CN (1) CN110987846B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948155A (zh) * 2020-07-30 2020-11-17 中国科学院西安光学精密机械研究所 精细全光谱结合改进gs-svr的复杂水体硝酸盐定量分析方法
CN112014344A (zh) * 2020-08-21 2020-12-01 浙江全世科技有限公司 一种污水在线监测方法
CN112504983A (zh) * 2020-12-11 2021-03-16 中国科学院西安光学精密机械研究所 一种基于浊度色度补偿的硝酸盐浓度预测方法
CN113588591A (zh) * 2021-08-11 2021-11-02 江门市华讯方舟科技有限公司 一种用于蛋氨酸含量快速检验的方法
CN116008208A (zh) * 2023-03-27 2023-04-25 山东省科学院海洋仪器仪表研究所 一种海水硝酸盐浓度特征光谱波段的选择方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108107019A (zh) * 2017-12-15 2018-06-01 暨南大学 一种基于近红外光谱法快速检测玉米中杂色曲菌素含量的方法
CN109142249A (zh) * 2018-08-31 2019-01-04 青岛卓建海洋装备科技有限公司 一种海水硝酸盐浓度测试方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108107019A (zh) * 2017-12-15 2018-06-01 暨南大学 一种基于近红外光谱法快速检测玉米中杂色曲菌素含量的方法
CN109142249A (zh) * 2018-08-31 2019-01-04 青岛卓建海洋装备科技有限公司 一种海水硝酸盐浓度测试方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
李丹 等: "一种基于紫外光谱法的海水硝酸盐在线监测系统", 《光谱学与光谱分析》 *
杨鹏程 等: "间隔偏最小二乘紫外光谱法海水硝酸盐最佳建模波长区间选取", 《海洋环境科学》 *
汪晨: "紫外-可见光谱技术在连缸染色浓度分析中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》 *
胡莉军 等: "基于长光程吸收池的痕量有害气体定量检测研究", 《科学技术与工程》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111948155A (zh) * 2020-07-30 2020-11-17 中国科学院西安光学精密机械研究所 精细全光谱结合改进gs-svr的复杂水体硝酸盐定量分析方法
CN111948155B (zh) * 2020-07-30 2022-05-10 中国科学院西安光学精密机械研究所 精细全光谱结合改进gs-svr的复杂水体硝酸盐定量分析方法
CN112014344A (zh) * 2020-08-21 2020-12-01 浙江全世科技有限公司 一种污水在线监测方法
CN112014344B (zh) * 2020-08-21 2022-11-22 浙江全世科技有限公司 一种污水在线监测方法
CN112504983A (zh) * 2020-12-11 2021-03-16 中国科学院西安光学精密机械研究所 一种基于浊度色度补偿的硝酸盐浓度预测方法
CN113588591A (zh) * 2021-08-11 2021-11-02 江门市华讯方舟科技有限公司 一种用于蛋氨酸含量快速检验的方法
CN113588591B (zh) * 2021-08-11 2023-05-26 江门市华讯方舟科技有限公司 一种用于蛋氨酸含量快速检验的方法
CN116008208A (zh) * 2023-03-27 2023-04-25 山东省科学院海洋仪器仪表研究所 一种海水硝酸盐浓度特征光谱波段的选择方法

Also Published As

Publication number Publication date
CN110987846B (zh) 2021-04-20

Similar Documents

Publication Publication Date Title
CN110987846B (zh) 一种基于iPLS-PA算法的硝酸盐浓度预测方法
JP3245157B2 (ja) スペクトルデータの測定および補正
US7570357B2 (en) Visible/near-infrared spectrometry and its device
CN103983595B (zh) 一种基于紫外-可见光谱处理的水质浊度解算方法
CN104062256B (zh) 一种基于近红外光谱的软测量方法
CN111965140B (zh) 基于特征峰的波长点重组方法
CN106918567B (zh) 一种测量痕量金属离子浓度的方法和装置
CN101413885A (zh) 一种快速定量蜂蜜品质的近红外光谱方法
CN101825567A (zh) 一种近红外光谱和拉曼光谱波长的筛选方法
CN109669023A (zh) 一种基于多传感器融合的土壤属性预测方法
CN105486655A (zh) 基于红外光谱智能鉴定模型的土壤有机质快速检测方法
CN109085136B (zh) 近红外漫反射光谱测定水泥生料氧化物成分含量的方法
CN109540836A (zh) 基于bp人工神经网络的近红外光谱糖度检测方法及系统
CN103592255A (zh) 一种基于近红外光谱技术的阿胶化皮液中总蛋白含量的软测量方法
CN106950192A (zh) 一种基于近红外光谱分析技术的植物蛋白饮料中主要成分含量快速检测的方法
CN111855608A (zh) 基于融合特征波长选择算法的苹果酸度近红外无损检测方法
CN105784672A (zh) 一种基于双树复小波算法的毒品检测仪标准化方法
CN111999258B (zh) 一种面向光谱基线校正的加权建模局部优化方法
CN107505283A (zh) 近红外光谱测定次生盐渍土壤中硝酸根离子含量的方法
CN111141809B (zh) 一种基于非接触式电导信号的土壤养分离子含量检测方法
CN111896497A (zh) 一种基于预测值的光谱数据修正方法
CN115165770B (zh) 基于宽光谱及bpnn的水体cod与浊度同时检测方法
CN110174362B (zh) 一种检测中性糖和酸性糖含量的方法
CN109459408A (zh) 一种基于稀疏回归lar算法的近红外定量分析方法
CN105651735A (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