CN109828570B - 一种自适应边界层水面无人艇控制导引方法 - Google Patents
一种自适应边界层水面无人艇控制导引方法 Download PDFInfo
- Publication number
- CN109828570B CN109828570B CN201910120131.6A CN201910120131A CN109828570B CN 109828570 B CN109828570 B CN 109828570B CN 201910120131 A CN201910120131 A CN 201910120131A CN 109828570 B CN109828570 B CN 109828570B
- Authority
- CN
- China
- Prior art keywords
- speed
- tracking
- target point
- boundary layer
- current
- 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
Images
Abstract
本发明属于控制领域,具体涉及一种自适应边界层水面无人艇控制导引方法。包括神经网络离线训练,初始化算法控制参数,获取规划点及传感器信息,判断规划点是圆弧还是直线,根据当前水面无人艇至目标点距离(若为圆,则为已完成的跟踪角度),通过安全阈值观测器判断是否达到目标点,如果到达目标点,则将上一个目标点删除,再跳至步骤5,否则输出期望航向和期望航速。本发明通过自适应边界层在LOS算法中引入航速的考虑,提高其跟踪响应速度,利用双曲正切修正器优化LOS算法的跟踪控制精度。同时本算法采取模块化设计,不论是自适应边界层中的水面无人艇制动长度神经网络,速度优化层还是基于向量场改进的圆弧导引策略,都能够运用在其他算法中。
Description
技术领域
本发明属于控制领域,具体涉及一种自适应边界层水面无人艇控制导引方法。
背景技术
水面无人艇将传统船舶技术与无人技术相结合,具有无人自主、机动性强、隐身性能好、成本较低等优点。水面无人艇大多在恶劣海况下执行任务,因此高精度,强适应能力的跟踪控制方法是实现自主的基本。如何实现水面无人艇高精度控制,强适应能力,同时兼顾最优速度已成为现在的研究热点。视线法(Line-of-Sight,LOS)作为工程中常见的跟踪控制方法,在已有的文献中,许多学者针对该方法进行了改进从而提高其自适应能力。但是这些方法大多基于艇体数学模型,忽略速度耦合,同时敏感度过高,故其鲁棒性和普适性较差。
挪威科技大学的Fossen I等首先将LOS跟踪算法引入到欠驱动水面无人艇中,将横向距离,纵向距离,艏向角度这三个自由度的跟踪问题转化为艏向角度,前进速度的二自由度跟踪问题,其缺点是收敛半径为常数,故不具备自适应能力,且为只能实现局部收敛。
哈尔滨工程大学的廖煜雷等人,通过改进LOS算法,在视线角中引入积分项,增强其抗扰动能力,同时通过线性反馈设计自适应的收敛半径使得算法能够实现全局收敛。但其算法参数过多,且敏感度高,调节困难。
海军工程大学的陈宵等人,将在LOS中的积分项中引入一种漂角观测器,以减小定积分增益选择不当带来的积分饱和以及超调等不利影响,在此基础上,在前视距离中引入比例增益使其变为时间变量,从而使得无人艇操纵更加灵活。通过级联系统理论证明了该闭环系统为全局k-指数稳定,但同样具有敏感度高的问题,同时忽略了速度耦合的影响。
武汉理工大学的陈俊等人,考虑到速度耦合的影响,提出基于LOS的变速循迹控制,通过状态反馈和反步法分别计算艇体加速度,再利用最小二乘法对其进行估计,使得水面无人艇能够根据横向偏差以及变化率调整速度,更快速地达到预设路径,但该算法基于模型设计,且只考虑直线跟踪,普适性较差。
发明内容
本发明的目的在于提供一种基于无模型的、高精度的、高鲁棒性且具备一定自适应能力的自适应边界层水面无人艇控制导引方法。
一种自适应边界层水面无人艇控制导引方法,具体包括如下步骤:
步骤1、收集对当前速度制动长度有影响的参数,包括航速、海风的风速与风向、海流的流速和流向、波浪的波高与波长及波向,制动长度五个体系指标,并对这五个体系指标进行识别与处理以及无量纲化;
步骤2、采用水面无人艇制动长度神经网络对所述的五个体系指标进行训练,然后对其进行泛化验证;
步骤3、初始化双曲正切修正器,需要初始化的控制参数为tanh系数,增益系数,收敛半径;
其中,δ表示横侧偏差,ktanh表示tanh系数,ξ表示增益系数,表示收敛半径,ktanh在横侧偏差越接近0时越影响收敛速度,ktanh越大,收敛速度越慢;ξ抑制横侧偏差的振荡,同时也会降低收敛速度,影响收敛速度的快慢,过大会增加横侧偏差的超调。
步骤4、通过传感器获取水面无人艇当前位姿势,包括航速V,航向角χ,经度Lat,纬度Lon,以及海洋环境信息,通过规划系统获取规划路径p(p1,p2,......,pn);
步骤5、根据当前给出的规划点,判断跟踪直线还是圆,如果是直线跳至步骤七,如果是圆跳至步骤八;
步骤6、根据当前目标点p1和上一个目标点p0,求解期望艏向:
步骤7、根据当前目标点p1即圆心,跟踪方向ζ以及跟踪半径r,求解期望艏向:
其中,kc为控制参数,kc越大,越能抑制超调,同时抑制振荡,当过大会令横侧偏差存在静态误差;χorbit为舰体和圆心的连线与正北坐标系的夹角,ζ为跟踪方向,若ζ=1则顺时针跟踪,反之亦然;
步骤8、根据航速V,以及传感器得到的环境信息,通过水面无人艇制动神经网络预测其制动距离,得到当前规划路径的边界层,然后结合当前边界层和制动距离得到期望速度:
其中,umax为最大航速,u*为最优航速;
步骤9、根据当前水面无人艇至目标点距离,对比安全阈值ε,判断是否达到目标点,若为圆,则为已完成的跟踪角度;如果到达目标点则删除上一个目标点,并跳至步骤5;否则,输出期望航速和航向。
所述一种自适应边界层水面无人艇控制导引方法,步骤2中神经网络模型具备一个输入层,两个全连接层,一个输出层,返回一个单一且连续的数值,其中两层全连接层的神经元数均为64个,激活函数采用ReLU,采用RMS进行优化,学习律设定为0.001,损失函数采用交叉熵,具体过程如下所示:
式中,a为当前神经网络输出值,y为目标值,x为样本,m为该批量样本数量;
式中,g为反向梯度,θ为该神经元权重,▽θ表示对权重进行求偏导;
式中,r为累计平方梯度,初始值为0;为衰减速率,ρ给定常值0.9;
式中,ε为学习律,给定常值0.001,δ为常数10-10。
θ←θ+Δθ
所述一种自适应边界层水面无人艇控制导引方法,步骤6中χp,δ由下式确定:
χp=a tan2(p0,p1)δ=(nt-nk)sinχp+(et-ek)cosχp
其中,atan2表示四象限反正切函数,χp的方向为p0指向p1,(nt,yt)为水面无人艇经纬度转坐标,为目标点经纬度转xy坐标。
所述一种自适应边界层水面无人艇控制导引方法,步骤8中所述边界层的含义为当艇体在远离跟踪路线的时候,加快速度使其尽快在路径周围,当艇体在路线周围时,采用最优的速度去跟踪路径,加速横侧偏差的收敛,当艇体稳定在路线上时,采用最大速度尽快完成本次路径跟踪,
u2 Target=u*,DtoTarget<Ldecelerate
其中,Ldecelerate为当前航速V的制动距离,DtoTarget为水面无人艇至目标点距离,最终的期望速度优先级为u1 Target<u2 Target。
本发明的有益效果在于:
本发明通过自适应边界层在LOS算法中引入航速的考虑,提高其跟踪响应速度,利用双曲正切修正器优化LOS算法的跟踪控制精度。同时本算法采取模块化设计,不论是自适应边界层中的水面无人艇制动长度神经网络,速度优化层还是基于向量场改进的圆弧导引策略,都能够运用在其他算法中。
附图说明
图1为自适应边界层的欠驱动水面无人艇控制导引方法;
图2为自适应边界层流程图;
图3为安全阈值观测器流程图;
图4为基于自适应边界层水面无人舰控制导引方法概念图。
具体实施方式
下面结合附图对本发明作进一步描述。
如附图1所示,为自适应边界层的欠驱动水面无人艇控制导引方法;如附图2所示,为自适应边界层流程图;如附图3所示,为安全阈值观测器流程图;如附图4所示,为基于自适应边界层水面无人舰控制导引方法概念图。
一种自适应边界层水面无人艇控制导引方法,具体包括如下步骤:
步骤1、收集对当前速度制动长度有影响的参数,包括航速、海风的风速与风向、海流的流速和流向、波浪的波高与波长及波向,制动长度五个体系指标,并对这五个体系指标进行识别与处理以及无量纲化;
步骤2、采用水面无人艇制动长度神经网络对所述的五个体系指标进行训练,然后对其进行泛化验证;
步骤3、初始化双曲正切修正器,需要初始化的控制参数为tanh系数,增益系数,收敛半径;
其中,δ表示横侧偏差,ktanh表示tanh系数,ξ表示增益系数,表示收敛半径,ktanh在横侧偏差越接近0时越影响收敛速度,ktanh越大,收敛速度越慢;ξ抑制横侧偏差的振荡,同时也会降低收敛速度,影响收敛速度的快慢,过大会增加横侧偏差的超调。
步骤4、通过传感器获取水面无人艇当前位姿势,包括航速V,航向角χ,经度Lat,纬度Lon,以及海洋环境信息,通过规划系统获取规划路径p(p1,p2,......,pn);
步骤5、根据当前给出的规划点,判断跟踪直线还是圆,如果是直线跳至步骤七,如果是圆跳至步骤八;
步骤6、根据当前目标点p1和上一个目标点p0,求解期望艏向:
步骤7、根据当前目标点p1即圆心,跟踪方向ζ以及跟踪半径r,求解期望艏向:
其中,kc为控制参数,kc越大,越能抑制超调,同时抑制振荡,当过大会令横侧偏差存在静态误差;χorbit为舰体和圆心的连线与正北坐标系的夹角,ζ为跟踪方向,若ζ=1则顺时针跟踪,反之亦然;
步骤8、根据航速V,以及传感器得到的环境信息,通过水面无人艇制动神经网络预测其制动距离,得到当前规划路径的边界层,然后结合当前边界层和制动距离得到期望速度:
其中,umax为最大航速,u*为最优航速;
步骤9、根据当前水面无人艇至目标点距离,对比安全阈值ε,判断是否达到目标点,若为圆,则为已完成的跟踪角度;如果到达目标点则删除上一个目标点,并跳至步骤5;否则,输出期望航速和航向。
所述一种自适应边界层水面无人艇控制导引方法,步骤2中神经网络模型具备一个输入层,两个全连接层,一个输出层,返回一个单一且连续的数值,其中两层全连接层的神经元数均为64个,激活函数采用ReLU,采用RMS进行优化,学习律设定为0.001,损失函数采用交叉熵,具体过程如下所示:
式中,a为当前神经网络输出值,y为目标值,x为样本,m为该批量样本数量;
式中,g为反向梯度,θ为该神经元权重,▽θ表示对权重进行求偏导;
式中,r为累计平方梯度,初始值为0;为衰减速率,ρ给定常值0.9;
式中,ε为学习律,给定常值0.001,δ为常数10-10。
θ←θ+Δθ
所述一种自适应边界层水面无人艇控制导引方法,步骤6中χp,δ由下式确定:
χp=a tan2(p0,p1)δ=(nt-nk)Sinχp+(et-ek)cosχp
其中,atan2表示四象限反正切函数,χp的方向为p0指向p1,(nt,yt)为水面无人艇经纬度转坐标,为目标点经纬度转xy坐标。
所述一种自适应边界层水面无人艇控制导引方法,步骤8中所述边界层的含义为当艇体在远离跟踪路线的时候,加快速度使其尽快在路径周围,当艇体在路线周围时,采用最优的速度去跟踪路径,加速横侧偏差的收敛,当艇体稳定在路线上时,采用最大速度尽快完成本次路径跟踪,
u2 Target=u*,DtoTarget<Ldecelerate
Claims (4)
1.一种自适应边界层水面无人艇控制导引方法,其特征在于,具体包括如下步骤:
步骤1、收集对当前速度制动长度有影响的参数,包括航速、海风的风速与风向、海流的流速和流向、波浪的波高与波长及波向,制动长度五个体系指标,并对这五个体系指标进行识别与处理以及无量纲化;
步骤2、采用水面无人艇制动长度神经网络对所述的五个体系指标进行训练,然后对其进行泛化验证;
步骤3、初始化双曲正切修正器Rk和χmax,需要初始化的控制参数为tanh系数,增益系数,收敛半径;
其中,δ表示横侧偏差,ktanh表示tanh系数,ξ表示增益系数,表示收敛半径,横侧偏差越接近0,ktanh越影响收敛速度;ktanh越大,收敛速度越慢;ξ抑制横侧偏差的振荡,同时也会降低收敛速度,影响收敛速度的快慢,过大会增加横侧偏差的超调,Rk代表直线跟踪的修正器,χmax代表圆弧跟踪的修正器;
步骤4、通过传感器获取水面无人艇当前位姿势,包括航速V,航向角χ,经度Lat,纬度Lon,以及海洋环境信息,通过规划系统获取规划路径p={p0,p1,p2,......pn};
步骤5、根据当前给出的规划点,判断跟踪直线还是圆,如果是直线则跳至步骤6,如果是圆则跳至步骤7;
步骤6、根据当前目标点p1和上一个目标点p0,求解期望艏向
步骤7、根据当前目标点p1即圆心,跟踪方向ζ以及跟踪半径r,求解期望艏向
其中,kc为控制参数,kc越大,越能抑制超调,同时抑制振荡,当过大会令横侧偏差存在静态误差;χorbit为舰体和圆心的连线与正北坐标系的夹角,ζ为跟踪方向,若ζ=1则顺时针跟踪,反之亦然;
步骤8、根据航速V,以及传感器得到的环境信息,通过水面无人艇制动神经网络预测其制动距离,得到当前规划路径的边界层,然后结合当前边界层和制动距离得到期望速度
其中,umax为最大航速,u*为最优航速;
步骤9、根据当前水面无人艇至目标点距离,对比安全阈值ε,判断是否达到目标点,若为圆,则为已完成的跟踪角度;如果到达目标点则删除上一个目标点,并跳至步骤5;否则,输出期望航速和航向。
2.根据权利要求1所述一种自适应边界层水面无人艇控制导引方法,其特征在于,步骤2中所述神经网络具备一个输入层,两个全连接层,一个输出层,返回一个单一且连续的数值,其中两层全连接层的神经元数均为64个,激活函数采用ReLU,采用RMS进行优化,学习律设定为0.001,损失函数采用交叉熵,具体过程如下所示:
式中,a为当前神经网络输出值,y为目标值,x为样本,m为批量样本数量;
ω←ρω+(1-ρ)g⊙g
式中,ω为累计平方梯度,初始值为0;为衰减速率,ρ给定常值0.9;
式中,s为学习律,给定常值0.001,μ为常数10-10,
θ←θ+Δθ。
3.根据权利要求1所述一种自适应边界层水面无人艇控制导引方法,其特征在于,步骤6中所述χp,δ由下式确定
χp=atan2(p0,p1) δ=(nt-nk)sinχp+(et-ek)cosχp
其中,atan2表示四象限反正切函数,χp的方向为p0指向p1,(nt,et)为水面无人艇经纬度转xy坐标,(nk,ek)为目标点经纬度转xy坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910120131.6A CN109828570B (zh) | 2019-02-18 | 2019-02-18 | 一种自适应边界层水面无人艇控制导引方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910120131.6A CN109828570B (zh) | 2019-02-18 | 2019-02-18 | 一种自适应边界层水面无人艇控制导引方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109828570A CN109828570A (zh) | 2019-05-31 |
CN109828570B true CN109828570B (zh) | 2021-10-01 |
Family
ID=66862197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910120131.6A Active CN109828570B (zh) | 2019-02-18 | 2019-02-18 | 一种自适应边界层水面无人艇控制导引方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109828570B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442135B (zh) * | 2019-08-06 | 2020-08-21 | 南京赛沃夫海洋科技有限公司 | 一种基于改进遗传算法的无人艇路径规划方法及系统 |
CN110609553B (zh) * | 2019-09-16 | 2022-07-15 | 哈尔滨工程大学 | 一种用于铺管船舶圆弧路径的los导引控制方法 |
CN110631587B (zh) * | 2019-09-26 | 2021-05-25 | 浙江圣奥家具制造有限公司 | 一种自动归位椅复合路径规划方法 |
CN110716574B (zh) * | 2019-09-29 | 2023-05-02 | 哈尔滨工程大学 | 一种基于深度q网络的uuv实时避碰规划方法 |
CN111045332A (zh) * | 2019-12-27 | 2020-04-21 | 哈尔滨工程大学 | 一种无人艇路径跟踪导引策略和扰动补偿方法 |
CN111506086A (zh) * | 2020-05-22 | 2020-08-07 | 中国石油大学(华东) | 一种改进los引导律结合模糊pid的无人船路径跟踪控制方法 |
CN112393731B (zh) * | 2020-10-10 | 2023-04-25 | 上海钛米机器人股份有限公司 | 跟踪路径的方法、装置、电子设备、及存储介质 |
CN112346465B (zh) * | 2020-11-27 | 2022-09-02 | 哈尔滨工程大学 | 基于ialos导引律的欠驱无人船自适应模糊控制方法 |
CN113110458B (zh) * | 2021-04-19 | 2023-09-01 | 大连海事大学 | 一种无人船虚拟目标跟踪控制系统 |
CN112947494B (zh) * | 2021-04-22 | 2022-10-04 | 中国船舶重工集团公司第七0七研究所九江分部 | 一种基于模糊pid的船舶自动靠泊控制方法 |
CN113625725B (zh) * | 2021-09-02 | 2024-05-07 | 中国舰船研究设计中心 | 一种水面无人艇路径跟踪控制方法 |
CN114063622B (zh) * | 2022-01-06 | 2022-05-17 | 陕西欧卡电子智能科技有限公司 | 无人船自主停泊定位方法、装置及相关组件 |
CN114489091B (zh) * | 2022-01-21 | 2023-01-13 | 湖南大学无锡智能控制研究院 | 一种自主水下机器人路径跟踪的制导律控制方法、设备及系统 |
CN115993831B (zh) * | 2023-03-23 | 2023-06-09 | 安徽大学 | 基于深度强化学习的机器人无目标网络的路径规划方法 |
CN117647981A (zh) * | 2023-11-23 | 2024-03-05 | 中船(北京)智能装备科技有限公司 | 一种无人艇智能航行控制方法、装置及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012034119A2 (en) * | 2010-09-11 | 2012-03-15 | Wayne State University | Guidance and control system for under-actuated marine surface ships and other automomous platforms |
CN104765370A (zh) * | 2015-04-15 | 2015-07-08 | 哈尔滨工业大学 | 有环境干扰下考虑侧滑角的uuv航迹视线导引法 |
CN108267955A (zh) * | 2018-01-17 | 2018-07-10 | 哈尔滨工程大学 | 面向无人艇自主靠泊的运动控制方法 |
CN108415423A (zh) * | 2018-02-02 | 2018-08-17 | 哈尔滨工程大学 | 一种高抗扰自适应路径跟随方法及系统 |
CN108459602A (zh) * | 2018-02-28 | 2018-08-28 | 哈尔滨工程大学 | 多障碍复杂环境下欠驱动无人艇的自主靠泊方法 |
CN109283842A (zh) * | 2018-08-02 | 2019-01-29 | 哈尔滨工程大学 | 一种无人艇航迹跟踪智能学习控制方法 |
-
2019
- 2019-02-18 CN CN201910120131.6A patent/CN109828570B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012034119A2 (en) * | 2010-09-11 | 2012-03-15 | Wayne State University | Guidance and control system for under-actuated marine surface ships and other automomous platforms |
CN104765370A (zh) * | 2015-04-15 | 2015-07-08 | 哈尔滨工业大学 | 有环境干扰下考虑侧滑角的uuv航迹视线导引法 |
CN108267955A (zh) * | 2018-01-17 | 2018-07-10 | 哈尔滨工程大学 | 面向无人艇自主靠泊的运动控制方法 |
CN108415423A (zh) * | 2018-02-02 | 2018-08-17 | 哈尔滨工程大学 | 一种高抗扰自适应路径跟随方法及系统 |
CN108459602A (zh) * | 2018-02-28 | 2018-08-28 | 哈尔滨工程大学 | 多障碍复杂环境下欠驱动无人艇的自主靠泊方法 |
CN109283842A (zh) * | 2018-08-02 | 2019-01-29 | 哈尔滨工程大学 | 一种无人艇航迹跟踪智能学习控制方法 |
Non-Patent Citations (2)
Title |
---|
基于切换视线法的欠驱动无人艇鲁棒自适应路径跟踪控制;曾江峰 等;《兵工学报》;20181231;第39卷(第12期);第2427-2437页全文 * |
基于非对称模型的欠驱动USV路径跟踪控制;董早鹏 等;《中国造船》;20160331;第57卷(第1期);第116-126页全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109828570A (zh) | 2019-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109828570B (zh) | 一种自适应边界层水面无人艇控制导引方法 | |
CN107168335B (zh) | 一种考虑混合多目标避障的水面无人艇路径跟踪制导方法 | |
CN106950974B (zh) | 欠驱动自主水下航行器的对三维路径进行理解及跟踪控制方法 | |
Wang et al. | A paradigm for path following control of a ribbon-fin propelled biomimetic underwater vehicle | |
CN108052009B (zh) | 基于滤波反步法的水面目标救援跟踪观测控制器设计方法 | |
CN112130566B (zh) | 一种基于模糊逻辑和滑模控制策略的无人艇、无人机混合编队控制方法及其控制系统 | |
Sawada et al. | Path following algorithm application to automatic berthing control | |
CN108415423B (zh) | 一种高抗扰自适应路径跟随方法及系统 | |
Amundsen et al. | Autonomous ROV inspections of aquaculture net pens using DVL | |
Kong et al. | Dynamics modeling and motion control of an new unmanned underwater vehicle | |
CN106444838A (zh) | 一种自主水下航行器的精确航路跟踪控制方法 | |
CN112650231B (zh) | 一种实现避碰和避障的欠驱船编队控制系统 | |
KR20230011310A (ko) | 선박의 자동 유도 방법, 선박의 자동 유도 프로그램, 선박의 자동 유도 시스템 및 선박 | |
Zhao et al. | Broken lines path following algorithm for a water-jet propulsion USV with disturbance uncertainties | |
CN113093804B (zh) | 基于反演滑模控制的无人艇编队控制方法和控制系统 | |
Woo et al. | Vector field based guidance method for docking of an unmanned surface vehicle | |
Grenestedt et al. | LORCA: A high performance USV with applications to surveillance and monitoring | |
CN111538242A (zh) | 一种不确定性估计和饱和补偿的动力定位t-s模糊抗饱和控制方法 | |
Nguyen et al. | Robust Adaptive Depth Control of hybrid underwater glider in vertical plane | |
CN113296505B (zh) | 一种基于速变los的无人船多模式路径跟踪控制方法 | |
CN115014355A (zh) | 一种双体无人船的定点返航调控方法及装置 | |
Miller | Model predictive ship trajectory tracking system based on line of sight method | |
Jiang et al. | Design of novel S-plane controller of autonomous underwater vehicle established on sliding mode control | |
Liu et al. | Active disturbance rejection path following control of USV based on fuzzy optimization | |
CN114035574A (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 |