CN113110511B - 一种基于广义模糊双曲模型的智能船舶航向控制方法 - Google Patents

一种基于广义模糊双曲模型的智能船舶航向控制方法 Download PDF

Info

Publication number
CN113110511B
CN113110511B CN202110546654.4A CN202110546654A CN113110511B CN 113110511 B CN113110511 B CN 113110511B CN 202110546654 A CN202110546654 A CN 202110546654A CN 113110511 B CN113110511 B CN 113110511B
Authority
CN
China
Prior art keywords
ship
fuzzy
control
input
course
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
CN202110546654.4A
Other languages
English (en)
Other versions
CN113110511A (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.)
Dalian Maritime University
Original Assignee
Dalian Maritime University
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 Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN202110546654.4A priority Critical patent/CN113110511B/zh
Publication of CN113110511A publication Critical patent/CN113110511A/zh
Application granted granted Critical
Publication of CN113110511B publication Critical patent/CN113110511B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/0206Control of position or course in two dimensions specially adapted to water vehicles

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于广义模糊双曲模型的智能船舶航向控制方法,包括:建立以船舶舵角为输入,船舶航向角及航向角速度为输出的非线性船舶航向控制模型;将船舶实际艏向角度与期望艏向角差值作为非线性船舶航向控制模型的输入信息,对该模型中的非线性函数进行逼近、获得智能模糊船舶航向控制模型,设计智能船舶航向控制的模糊虚拟控制函数;将存在饱和限制的输入信号u与控制信号v进行比较获得误差信号、将该误差信号进行反馈和辅助补偿处理获得辅助补偿信号、从而循环地对船舶舵角输入限制进行补偿;通过广义模糊双曲正切逼近方法建立模糊航向控制模型,采用饱和补偿方法分析辅助补偿信号,设计模糊自适应更新率,从而实现船舶航向控制过程。

Description

一种基于广义模糊双曲模型的智能船舶航向控制方法
技术领域
本发明涉及智能船舶控制领域,尤其涉及一种基于广义模糊双曲模型的智能船舶航向控制方法。
背景技术
自动舵是智能船舶最重要的操控系统之一,智能船舶运动具有大时滞、大惯性、强非线性等特点,在实际的船舶控制中包括建模误差、船舶载况变化引起的参数摄动和不确定的外部干扰如风、流等环境力作用均会对控制的精度造成影响甚至造成严重事故,但目前大多数基于状态反馈进行的船舶航向控制设计并不能妥善解决数学模型不确定性问题。针对这个问题模型中的不确定项来描述,目前主要结合神经网络或者模糊逻辑系统进行逼近处理。
然而在实际工程应用中神经网络与模糊逻辑系统需要获取较多设计参数,模型复杂程度较高,且航向角追踪误差存在受限和执行航向追踪任务精度有限的情况。因此,现有智能船舶航向控制成果中,考虑控制性能与控制成本的折中,使用成本较高不利于工程实现。
发明内容
根据现有技术存在的问题,本发明公开了一种基于广义模糊双曲模型的智能船舶航向控制方法,具体采用如下步骤:
考虑船舶航向控制的回转非线性和外界干扰,建立以船舶舵角为输入,船舶航向角及航向角速度为输出的非线性船舶航向控制模型;
将船舶实际艏向角度与期望艏向角差值作为非线性船舶航向控制模型的输入信息,对该模型中的非线性函数进行逼近、获得智能模糊船舶航向控制模型,设计智能船舶航向控制的模糊虚拟控制函数;
将存在饱和限制的输入信号u与控制信号v进行比较获得误差信号、将该误差信号进行反馈和辅助补偿处理获得辅助补偿信号、从而循环地对船舶舵角输入限制进行补偿;
通过广义模糊双曲正切逼近方法建立模糊航向控制模型,采用饱和补偿方法分析辅助补偿信号,设计模糊自适应更新率,从而基于模糊双曲正切的饱和船舶航向控制。
进一步的,所述非线性船舶航向控制模型为:
δ为船舶舵角,φ为船舶航向,其导数为船舶航向角变化速度,K为船舶的回转性指数,T为船舶的跟从性指数,由于船舶的回转性和直线稳定性之间是相互矛盾的,是/>的非线性函数,近似表示为:/> 为船舶非线性系数,是实值常数;定义x1=φ,/>u=δ,将公式(1)转化得到船舶航向角及航向角速度为输出的非线性船舶航向控制模型:
其中x=[x1 x2]T,f1(x1)=0,u是系统饱和极限的控制输入,y是系统的输出,d1=0,d2是未知有界的干扰满足|d2|≤d*,f2(x)是满足利普希兹条件的未知非线性项。
进一步的,所述获得智能模糊船舶航向控制模型时:
首先对输入变量进行线性变换,从而得到广义输入变量
x=[x1 (t),x2 (t),…,xn (t)]T (3)
其中xi =xz-dzjxi 的个数为wi为将xz线性变换的个数,dzj为xz线性变换点。
对船舶航向控制模型中的非线性项逼近采用如下描述:
其中为广义模糊双曲模型逼近误差,/>为最优参数向量,/>
结合上述广义模糊双曲正切逼近方法获得智能模糊船舶航向控制模型
采集船舶罗经设备中船舶实际艏向角与期望艏向角差值作为输入Z1=x1-yd,Z2=x22,设计智能船舶航向控制的模糊虚拟控制函数:
其中α2是智能船舶航向控制的模糊虚拟控制函数,yd是期望的船舶艏向角信号。
进一步的,对船舶舵角输入限制进行补偿时:
在船舶航向控制设计中,输入限制的船舶舵角u,-um≤u≤uM,其中,-um和uM分别表示根据船舶海上航行已知输入舵角u的最小值和最大值,饱和限制得到输入值表示为
其中,v表示船舶航向控制中的控制输入;
将存在饱和限制的输入信号u与控制信号v进行比较获得误差信号,采用辅助设计模块进行辅助补偿处理获得辅助补偿信号:
其中,C21>0,Δu=u-v,ε是一个正的设计参数,e是辅助设计模块引入的变量。
进一步的,通过广义模糊双曲正切逼近方法建立的模糊航向控制模型以及饱和补偿模块的辅助补偿信号,设计模糊自适应更新率:
其中误差变量Z2=x222是设计的虚拟控制率,ρi,ri为大于零的设计参数。
设计了基于模糊双曲正切的饱和船舶航向最终控制率:
其中K2>0为设计参数。
由于采用了上述技术方案,本发明提供的一种基于广义模糊双曲模型的智能船舶航向控制方法,该方法一方面针对船舶航向系统存在船舶舵角输入受限影响,从而引入辅助设计系统从而避免输入饱和对系统性能和稳定性的影响。另一方面本发明基于广义模糊双曲模型建立智能船舶自动控制算法,有效解决在实际的船舶控制中包括建模误差、船舶载况变化引起的参数摄动引起的未知非线性不确定项的问题。通过设计系统控制参数,可以以期望精度完成航向追踪任务。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为发明控制方法流程图
图2为本发明控制系统框图
图3为船舶航向历时曲线
图4为船舶航向跟踪误差历时曲线
图5为船舶舵角历时曲线
图6为船舶转首角速度历时曲线
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1所示的一种基于广义模糊双曲模型的智能船舶航向控制方法,具体包括如下步骤:
S1:通过工程的观点,采集归纳智能船舶相关信息,将船舶看成一个动态系统,以船舶控制量(舵角)作为系统输入,船舶运动状态(航向或首摇角速度)作为输出,考虑外界风浪流干扰和模型不确定性,建立了船舶舵角,船舶航向,艏摇角速度的非线性船舶运动数学模型;
式(1)中δ为舵角,φ为航向,K为船舶的回转性指数(每秒),T为船舶的跟从性指数(秒),均是船舶前进速度和船长的函数。另外,由于船舶的回转性和直线稳定性之间是相互矛盾的,是/>的非线性函数,可以近似表示为:/>为船舶非线性系数,是实值常数。
定义x1=φ,u=δ.可将公式(1)转化得到具有未知函数和外部扰动的非线性船舶航向控制系统模型
式(2)其中x=[x1 x2]Tf1(x1)=0,u是系统饱和极限的控制输入,y是系统的输出。d1=0,d2是未知有界的干扰满足|d2|≤d*。f2(x)是满足利普希兹条件的未知非线性项。
S2:利用广义模糊双曲模型对输入变量进行线性变换,进一步扩大模糊模型逼近范围,广义输入变量可以被表示为
x=[x1 (t),x2 (t),…,xn (t)]T (3)
公式(3)中xi =xz-dzjxi 的个数为wi为将xz线性变换的个数,dzj为xz线性变换点。通过扩展一系列模糊双曲基函数,可以表示GFHM,因此基于GFHM的非线性项可以表示如下:
应用广义模糊双曲模型对智能船舶自动控制系统中的非线性项f2(x)进行逼近,可以得到如下逼近模型:
其中为广义模糊双曲模型逼近误差,/>为最优参数向量,/>
将(4),(5)式结合,系统(2)可表示为如下模糊船舶航向控制系统:
S3:为了解决舵角输入的幅度和速率限制的输入饱和问题,引进一种辅助设计系统去消除输入饱和限制的影响
在船舶航向控制设计中,首先我们考虑船舶舵角u的输入限制如下:
-um≤u≤uM (7)
式(7)中,-um和uM分别表示已知输入舵角u的最小值和最大值(船舶舵角最大幅度一般为35°),
其中,v表示整个系统要设计的控制输入。
为了方便地考虑船舶航向系统饱和输入的影响,引入了一种辅助设计系统
式(9)中C21>0,Δu=u-v,ε是一个正的设计参数,e是辅助设计系统引入的变量。
S4:结合辅助设计系统提出了一种考虑输入饱和的基于广义模糊双曲模型的智能船舶自动控制系统的控制器设计方法。
智能船舶自动控制系统基于广义模糊双曲模型的自适应更新率为:
式中误差变量Z2=x222是设计的虚拟控制率,ρi,ri为大于零的设计参数。
设计了考虑输入饱和影响下智能船舶自动控制系统的最终控制率形式如下:
其中K2>0为设计参数
本实施例利用Matlab进行计算机仿真,以“育鲲”轮为例,验证本发明控制算法的有效性。根据实际情况选取跟踪信号模型为
其中,φm表征了船舶航向的理想系统性能,φr(t)是命令输入信号,其值从0度到30度,周期为500s
选择控制器设计参数K1=5,K2=0.01,C21=0.52,ε*=0.1
本实施例利用MATLAB进行计算机仿真研究,结果如图3-6所示,图1为船舶航向跟踪轨迹,实线为实际航向轨迹,虚线为参考轨迹。图2是跟踪误差。从图中可以看出,在跟踪过程中误差往往在一个以零为中心的小邻域内,说明跟踪效果较好;图3和图4给出了船舶输入舵角和转向角速度的时移曲线。控制系统状态均有大的超调量,并且趋于稳定,表明本发明满足要求。上图验证了本发明提出的考虑输入饱和的基于GFHM的模糊自适应控制算法的有效性。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (5)

1.一种基于广义模糊双曲模型的智能船舶航向控制方法,其特征在于包括:
考虑船舶航向控制的回转非线性和外界干扰,建立以船舶舵角为输入,船舶航向角及航向角速度为输出的非线性船舶航向控制模型;
将船舶实际艏向角度与期望艏向角差值作为非线性船舶航向控制模型的输入信息,对该模型中的非线性函数进行逼近,获得智能模糊船舶航向控制模型,设计智能船舶航向控制的模糊虚拟控制函数;
将存在饱和限制的输入信号u与控制信号v进行比较获得误差信号,将该误差信号进行反馈和辅助补偿处理获得辅助补偿信号,从而循环地对船舶舵角输入限制进行补偿;
通过广义模糊双曲正切逼近方法建立模糊航向控制模型,采用饱和补偿方法分析辅助补偿信号,设计模糊自适应更新率,从而实现基于模糊双曲正切的饱和船舶航向控制过程。
2.根据权利要求1所述的方法,其特征在于:所述非线性船舶航向控制模型为:
δ为船舶舵角,φ为船舶航向,其导数为船舶航向角变化速度,K为船舶的回转性指数,T为船舶的跟从性指数,由于船舶的回转性和直线稳定性之间是相互矛盾的,是/>的非线性函数,近似表示为:/>ai为船舶非线性系数,是实值常数,i=1,2,3......,定义x1=φ,/>u=δ,将公式(1)转化得到船舶航向角及航向角速度为输出的非线性船舶航向控制模型:
其中x=[x1 x2]T,f1(x1)=0,u是系统饱和极限的控制输入,y是系统的输出,d1=0,d2是未知有界的干扰满足|d2|≤d*,f2(x)是满足利普希兹条件的未知非线性项。
3.根据权利要求1所述的方法,其特征在于:所述获得智能模糊船舶航向控制模型时:
首先对输入变量进行线性变换,从而得到广义输入变量
x=[x1 (t),x2 (t),…,xn (t)]T
其中xi =xz-dzjxi 的个数为wi为将xz线性变换的个数,dzj为xz线性变换点;
对船舶航向控制模型中的非线性项逼近采用如下描述:
其中εi为广义模糊双曲模型逼近误差,/>为最优参数向量,/>
结合上述广义模糊双曲正切逼近方法获得智能模糊船舶航向控制模型
采集船舶罗经设备中船舶实际艏向角与期望艏向角差值作为输入Z1=x1-yd,Z2=x22,设计智能船舶航向控制的模糊虚拟控制函数:
其中α2是设计的智能船舶航向控制的模糊虚拟控制律,yd是期望的船舶艏向角信号。
4.根据权利要求1所述的方法,其特征在于:对船舶舵角输入限制进行补偿时:
在船舶航向控制设计中,输入限制的船舶舵角u,-um≤u≤uM,其中,-um和uM分别表示根据船舶海上航行已知输入舵角u的最小值和最大值,饱和限制得到输入值表示为
其中,v表示船舶航向控制中的控制输入;
将存在饱和限制的输入信号u与控制信号v进行比较获得误差信号,采用辅助设计模块进行辅助补偿处理获得辅助补偿信号:
其中,C21>0,Δu=u-v,ε是一个正的设计参数,e是辅助设计模块引入的变量。
5.根据权利要求1所述的方法,其特征在于:通过广义模糊双曲正切逼近方法建立的模糊航向控制模型以及饱和补偿模块的辅助补偿信号,设计模糊自适应更新率:
其中误差变量Z2=x222是设计的虚拟控制律,ρi,ri为大于零的设计参数;
设计了基于模糊双曲正切的饱和船舶航向最终控制律:
其中K2>0为设计参数。
CN202110546654.4A 2021-05-19 2021-05-19 一种基于广义模糊双曲模型的智能船舶航向控制方法 Active CN113110511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110546654.4A CN113110511B (zh) 2021-05-19 2021-05-19 一种基于广义模糊双曲模型的智能船舶航向控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110546654.4A CN113110511B (zh) 2021-05-19 2021-05-19 一种基于广义模糊双曲模型的智能船舶航向控制方法

Publications (2)

Publication Number Publication Date
CN113110511A CN113110511A (zh) 2021-07-13
CN113110511B true CN113110511B (zh) 2024-04-05

Family

ID=76722786

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110546654.4A Active CN113110511B (zh) 2021-05-19 2021-05-19 一种基于广义模糊双曲模型的智能船舶航向控制方法

Country Status (1)

Country Link
CN (1) CN113110511B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113536463B (zh) * 2021-07-20 2024-04-02 大连海事大学 一种基于改进梯度下降法的神经网络船舶整体模型逼近方法
CN113625727A (zh) * 2021-09-13 2021-11-09 大连海事大学 一种基于广义模糊双曲模型控制存在复杂噪声的船舶航向系统的方法
CN114595898B (zh) * 2022-03-22 2023-05-02 中国船舶科学研究中心 一种考虑航道环境的船舶智能操纵方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014004911A (ja) * 2012-06-25 2014-01-16 Tokyo Univ Of Marine Science & Technology 非線形の自己回帰モデルによる船舶の航路保持方法
CN108762091A (zh) * 2018-06-25 2018-11-06 东南大学 一种基于未知控制方向的自适应编队控制算法
CN109933074A (zh) * 2019-04-10 2019-06-25 大连海事大学 一种有领导者的多无人船群集运动控制器结构及设计方法
CN110308726A (zh) * 2019-07-11 2019-10-08 上海海事大学 基于非线性反步法的欠驱动船舶航向控制方法
CN110687799A (zh) * 2019-11-12 2020-01-14 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN111880546A (zh) * 2020-08-11 2020-11-03 大连海事大学 一种基于自适应动态规划算法的虚拟引导船舶自动靠泊控制方法
CN111897225A (zh) * 2020-08-14 2020-11-06 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN112782981A (zh) * 2020-12-30 2021-05-11 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈指定性能控制方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8498796B2 (en) * 2007-02-12 2013-07-30 Deere & Company Perception model for trajectory following autonomous and human augmented speed control

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2014004911A (ja) * 2012-06-25 2014-01-16 Tokyo Univ Of Marine Science & Technology 非線形の自己回帰モデルによる船舶の航路保持方法
CN108762091A (zh) * 2018-06-25 2018-11-06 东南大学 一种基于未知控制方向的自适应编队控制算法
CN109933074A (zh) * 2019-04-10 2019-06-25 大连海事大学 一种有领导者的多无人船群集运动控制器结构及设计方法
CN110308726A (zh) * 2019-07-11 2019-10-08 上海海事大学 基于非线性反步法的欠驱动船舶航向控制方法
CN110687799A (zh) * 2019-11-12 2020-01-14 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN111880546A (zh) * 2020-08-11 2020-11-03 大连海事大学 一种基于自适应动态规划算法的虚拟引导船舶自动靠泊控制方法
CN111897225A (zh) * 2020-08-14 2020-11-06 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN112782981A (zh) * 2020-12-30 2021-05-11 大连海事大学 一种智能船舶自动舵系统的模糊自适应输出反馈指定性能控制方法及系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Application-oriented Adaptive Fuzzy Design of Ship’s Autopilot with Rudder Dynamics;Cheng Liu 等;《IEEE Xplore》;第127-131页 *
基于线性自抗扰控制的船舶航迹积分滑模控制器;邱峰 等;《上海海事大学学报》;第38卷(第03期);第12-17页 *
基于非线性增益递归滑模的船舶轨迹跟踪动态面自适应控制;沈智鹏;《自动化学报》;第44卷(第10期);第1833-1840页 *
船舶航向非线性系统自适应模糊补偿控制;朱冬健 等;《上海交通大学学报》;第49卷(第02期);第250-254页 *

Also Published As

Publication number Publication date
CN113110511A (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN113110511B (zh) 一种基于广义模糊双曲模型的智能船舶航向控制方法
Yang et al. Adaptive NN backstepping control design for a 3-DOF helicopter: Theory and experiments
CN110687799B (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN111650929B (zh) 一种自适应滑模控制方法、系统及移动机器人控制器
CN111273549B (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈容错控制方法及系统
CN110647154B (zh) 基于模糊状态观测器的智能船舶自动舵系统航向轨迹跟踪设计方法
CN111948937B (zh) 多智能体系统的多梯度递归强化学习模糊控制方法及系统
CN111290387B (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈指定性能控制方法及系统
CN112612209A (zh) 基于指令滤波神经网络控制器的全驱动船舶轨迹跟踪控制方法及系统
CN111897225B (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈控制方法及系统
CN111221335A (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈有限时间控制方法及系统
CN110703605B (zh) 一种面向智能船舶自动舵系统的自适应模糊最优控制方法及系统
CN113848887B (zh) 一种基于mlp方法的欠驱动无人艇轨迹跟踪控制方法
CN113467231A (zh) 基于侧滑补偿ilos制导律的无人艇路径跟踪方法
CN116719319B (zh) 水面无人艇的跟踪控制器设计方法
CN112782981A (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈指定性能控制方法及系统
CN111930124A (zh) 一种智能船舶自动舵系统的模糊自适应输出反馈有限时间控制方法及系统
Zhang Stable adaptive fuzzy sliding mode control of interconnected systems
Chen et al. Cooperative tracking of vessel trajectories based on curved dynamic coordinates
CN109828586B (zh) 基于非线性不确定性的无人船舶鲁棒h∞航向控制方法
CN113110527B (zh) 一种自主水下航行器有限时间路径跟踪的级联控制方法
CN112363538B (zh) 一种速度信息不完备下的auv区域跟踪控制方法
Witkowska et al. Adaptive backstepping tracking control for an over–actuated DP marine vessel with inertia uncertainties
Ihle et al. Robust formation control of marine craft using lagrange multipliers
CN112904719B (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