CN114035436A - Backstepping control method based on saturation adaptive law, storage medium and equipment - Google Patents
Backstepping control method based on saturation adaptive law, storage medium and equipment Download PDFInfo
- Publication number
- CN114035436A CN114035436A CN202111405311.2A CN202111405311A CN114035436A CN 114035436 A CN114035436 A CN 114035436A CN 202111405311 A CN202111405311 A CN 202111405311A CN 114035436 A CN114035436 A CN 114035436A
- Authority
- CN
- China
- Prior art keywords
- function
- saturation
- adaptive law
- derivative
- control
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 47
- 230000006870 function Effects 0.000 claims abstract description 69
- 239000000126 substance Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 238000012938 design process Methods 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 239000000758 substrate Substances 0.000 claims 1
- 238000010276 construction Methods 0.000 abstract 1
- 239000000047 product Substances 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000012467 final product Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
A backstepping control method, a storage medium and equipment based on a saturation adaptive law belong to the technical field of nonlinear system control. The method aims to solve the problems that the existing self-adaptive backstepping control method cannot process unknown nonlinear functions in a system and the existing saturation control is unsmooth. Aiming at a controlled object, a two-dimensional nonlinear system state space model is established, and two state variables exist in the two-dimensional nonlinear system; then according to the state variable x of the system1And constructing an error variable z from the target signal1According to the state variable x2And a virtual control function alpha to be designed1Construction of the error variable z2Designing a Lyapunov function, solving a first derivative of the Lyapunov function with respect to time, and designing a virtual control function alpha based on the first derivative of the Lyapunov function1And controlling the input u, and finally designing to obtain a saturation adaptive law. The method is mainly used for controlling the nonlinear system.
Description
Technical Field
The invention belongs to the technical field of nonlinear system control, and particularly relates to a backstepping control method based on a saturation adaptive law, a storage medium and equipment.
Background
The basic idea of the method is to estimate the unknown constant in the system by using the adaptive parameter, and then to complete the design of the control law by using the adaptive parameter and the backstepping control idea. At present, the adaptive backstepping control method has been widely applied to industrial automation systems, automobile control systems, robot control systems, and the like, and for the design details of the adaptive backstepping control method, reference may be made to chinese patent CN108869420B, chinese patent CN106787940B, and chinese patent CN 105573125A.
It should be noted that in the conventional adaptive back-stepping control method, the adaptive parameters are designed to only handle the uncertainty of the system parameters, i.e.: the adaptive parameters can only estimate the unknown constants in the system. When the system contains unknown nonlinear functions, the traditional self-adaptive backstepping control method is not applicable any more. For an actual nonlinear system, the system uncertainty is often caused by some unknown nonlinear functions, such as viscous friction nonlinearity, actuator dead zone/saturation nonlinearity, and the like. The traditional self-adaptive backstepping control method cannot process unknown nonlinear functions of the system, so that the method is greatly limited in practical application. Therefore, how to design an adaptive backstepping control strategy capable of handling the unknown nonlinear function of the system is a key problem.
Disclosure of Invention
The invention aims to solve the problems that the existing self-adaptive backstepping control method cannot process unknown nonlinear functions in a system and the existing saturation control is unsmooth.
A backstepping control method based on a saturation adaptive law is characterized in that a two-dimensional nonlinear system state space model is established for a controlled object, and the controlled object is controlled by the saturation adaptive law;
the design process of the saturation adaptive law comprises the following steps:
the method comprises the following steps of firstly, establishing a two-dimensional nonlinear system state space model aiming at a controlled object, wherein the established state space model of the two-dimensional nonlinear system is in the following specific form:
y=x1
wherein x is1,x2Represents the state variable of the system and is,representing the system state variable x1,x2The first derivative of (a), y is system output, f is a nonlinear smooth function unknown to the system, d is unknown disturbance outside the system, u represents a system control input signal, and the control aim is to design the system control input u so that the system output y tracks a given target signal;
step two, according to the state variable x of the system1And a target signal ydAnd the state variable x of the system2And a virtual control function alpha to be designed1Separately constructing the error variable z1And z2;
Step three, utilizing the error variable z obtained in the step two1And z2Designing a Lyapunov function V;
step four, solving the first derivative of the Lyapunov function V in the step three to obtain the first derivative
Step five, according to the first derivative of the Lyapunov functionDesigning a virtual control function alpha1And based on a virtual control function alpha1Design control input u:
wherein the content of the first and second substances,representing a virtual control function alpha1The first derivative of (a) is,for estimation of an unknown non-linear function f, k2Is a constant greater than zero;the function is a saturation function;
Step six, based on the first derivative of the Lyapunov functionAnd a virtual control function alpha1And finally designing a control input u to obtain a saturation adaptive law.
Furthermore, the controlled object is a motor system, and the corresponding system state variable x1,x2The rotational angle and the rotational speed of the motor system,the rotating speed and the angular acceleration of the motor system are shown, and the system output y is the rotating angle of the motor system.
Further, a solution of a state space model (1) of the two-dimensional nonlinear system exists and is unique; unknown non-linear smoothing function f andits first order second derivative satisfiesWhereinIs a constant greater than zero; target signal ydAnd its first and second derivatives are bounded.
Further, the error variable z in step two1And z2The following were used:
wherein alpha is1Representing the virtual control function to be designed.
Further, the first derivative of the lyapunov function V with respect to time in step three is:
wherein the content of the first and second substances,representing a virtual control function alpha1The first derivative of (a).
Further, the virtual control function α1As follows
Wherein k is1Is a constant greater than zero.
Further, the saturation adaptive law in the step six is designed as follows:
Wherein, mu1,μ2,μ3Is a constant greater than zero, η is a constant greater than zero, and s is a filter state variable.
A storage medium having stored therein at least one instruction that is loaded and executed by a processor to implement a method of backstepping control based on a saturated adaptive law.
An apparatus comprising a processor and a memory, the memory having stored therein at least one instruction that is loaded and executed by the processor to implement the method for saturation adaptive law based backstepping control.
The invention has the beneficial effects that:
the invention provides a backstepping control method based on a saturation self-adaptive law. Compared with the traditional self-adaptive backstepping control method which cannot process the unknown nonlinear function of the system, the method can directly estimate and process the unknown nonlinear function of the system; meanwhile, the saturation function designed by the invention can be conducted in a second order, so that the control effect of the system is smoother.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a system tracking performance curve under the method of the present invention;
FIG. 3 is a system tracking error curve under the method of the present invention;
FIG. 4 shows a system state x according to the method of the present invention2A response curve;
FIG. 5 is a diagram illustrating a system adaptive estimation performance curve according to the present invention;
FIG. 6 is a response curve of the system control input u according to the method of the present invention.
Detailed Description
The first embodiment is as follows: the present embodiment is described in connection with figure 1,
the embodiment is a backstepping control method based on a saturation adaptive law, aiming at a controlled object, establishing a two-dimensional nonlinear system state space model, and controlling the controlled object by using the saturation adaptive law;
the design process of the saturation adaptive law comprises the following steps:
step one, aiming at a controlled object, establishing a two-dimensional nonlinear system state space model, wherein two state variables x exist in the two-dimensional nonlinear system1And x2A control input and a given target signal yd;
Step two, according to the state variable x of the system1And a target signal ydAnd the state variable x of the system2And a virtual control function alpha to be designed1Separately constructing the error variable z1And z2;
Step three, utilizing the error variable z obtained in the step two1And z2Designing a Lyapunov function V;
step four, solving the first derivative of the Lyapunov function V in the step three to obtain the first derivative
Step five, according to the first derivative of the Lyapunov functionDesigning a virtual control function alpha1And a control input u;
step sixFirst derivative based on Lyapunov functionAnd a virtual control function alpha1And finally designing a control input u to obtain a saturation adaptive law.
The second embodiment is as follows:
the embodiment is a backstepping control method based on a saturation adaptive law, and in the step one, the specific form of the established state space model of the two-dimensional nonlinear system is as follows:
wherein x is1,x2Represents the state variable of the system and is,representing the system state variable x1,x2Is the first derivative of (a), y is the system output, f is a nonlinear smooth function unknown to the system, d is an unknown disturbance external to the system, u represents the system control input signal, the control objective is to design the system control input u so that the system output y tracks a given target signal yd。
The controlled object of the invention can be a motor system, and can also be other systems or objects, and when the controlled object is the motor system, the corresponding system state variable x1,x2The rotational angle and the rotational speed of the motor system,the rotating speed and the angular acceleration of the motor system are shown, and the system output y is the rotating angle of the motor system.
Other steps and parameters are the same as in the first embodiment.
The third concrete implementation mode:
the embodiment is a backstepping control method based on a saturation self-adaption law, and the solution of a state space model (1) of the two-dimensional nonlinear system exists and is unique; unknown non-linear smoothnessFunction f and its first and second derivatives satisfyWhereinIs a constant greater than zero; target signal ydAnd its first and second derivatives are bounded.
Other steps and parameters are the same as in the second embodiment.
The fourth concrete implementation mode:
the embodiment is a backstepping control method based on a saturation adaptive law, and an error variable z in the step two1And z2The following were used:
wherein alpha is1Representing the virtual control function to be designed.
Other steps and parameters are the same as in one of the first to third embodiments.
The fifth concrete implementation mode:
the present embodiment is a back-stepping control method based on the saturation adaptive law, which uses the error variable z set in the second step1And z2The designed Lyapunov function V is as follows:
other steps and parameters are the same as in embodiment four.
The sixth specific implementation mode:
the embodiment is a reverse step control method based on a saturation adaptive law, and the first derivative of the lyapunov function V in the third step to the time is:
wherein the content of the first and second substances,representing a virtual control function alpha1The first derivative of (a).
Other steps and parameters are the same as those in the fifth embodiment.
The seventh embodiment:
the embodiment is a backstepping control method based on a saturation adaptive law, and a virtual control function alpha1And the control input u is:
wherein the content of the first and second substances,representing a virtual control function alpha1The first derivative of (a) is,for estimation of an unknown non-linear function f, k1、k2Respectively being a constant greater than zero
The saturation function of the invention is different from the traditional saturation function, the traditional saturation function is continuous but can not be conducted, and the saturation function designed by the invention can be conducted in a second order, so that the system control effect is smoother.
Other steps and parameters are the same as in one of the first to sixth embodiments.
The specific implementation mode is eight:
the embodiment is a backstepping control method based on a saturation adaptive law, and the saturation adaptive law in the sixth step is designed as follows:
where η is a constant greater than zero and s is a filter state variable.
The traditional adaptive rate only estimates the system unknown constant, and the adaptive law designed by the invention can directly estimate the system unknown function.
The other steps and parameters are the same as in the seventh embodiment.
It will be demonstrated below that the back-stepping controller (6) based on the saturated adaptation law enables the system tracking error to converge to a smaller neighborhood near the origin. The demonstration process is as follows:
substituting the formula (6) into the formula (4) to obtain the final product
is obtainable from the formula (10)
The formula indicates the error variable z1And z2Are bounded.
To z2Derived by derivation
From f, z1And z2Are all bounded and availableIs bounded, so there is a constant for the filter (9)So that
Where ξ is the filter estimation error.
Substituting the formulas (12) and (13) into the formula (8) to obtain the final product
If it isSelecting Lyapunov functionCan obtain the productThe first derivative with respect to time is
Wherein a suitable c is selected5So thatCan obtain the productThe constant c is chosen such thatIt can be obtained that t is more than or equal to all t1,This is true.
Thus, it is possible to obtain t ≧ t for all t ≧ t1Is obtained by
By substituting formula (19) for formula (20)
From the formula (21)
From the formulae (19) and (22)
Equation (23) shows that as time goes to infinity, the system tracking error converges to a small neighborhood around the origin, and the result is verified.
The specific implementation method nine:
the embodiment is a storage medium, wherein at least one instruction is stored in the storage medium, and the at least one instruction is loaded and executed by a processor to realize a backstepping control method based on a saturation adaptive law.
The detailed implementation mode is ten:
the embodiment is an apparatus, which includes a processor and a memory, where the memory stores at least one instruction, and the at least one instruction is loaded and executed by the processor to implement a method for controlling backstepping based on a saturation adaptive law.
Example one
For the system (1), the initial value is taken as x1(0)=0,x2(0) 0, function f 1.25sin (x)1x2) +0.25cos (0.5t), system target signal ydSint. Taking the controller parameter as k1=3,k2=0.5,c=100,η=50,μ1=100,μ2=80,μ 35. The system sampling interval time was 0.002 seconds.
FIG. 2 is a system tracking performance curve under the method of the present invention; FIG. 3 is a system tracking error curve under the method of the present invention; FIG. 4 shows a system state x according to the method of the present invention2A response curve; FIG. 5 shows a system for implementing the method of the present inventionSelf-adaptively estimating a performance curve; FIG. 6 is a response curve of the system control input u according to the method of the present invention;
and (4) conclusion: from fig. 2, the back-stepping controller based on the saturation adaptive law designed by equation (6) can make the system tracking error converge to a smaller neighborhood near the origin.
The above-described calculation examples of the present invention are merely to explain the calculation model and the calculation flow of the present invention in detail, and are not intended to limit the embodiments of the present invention. It will be apparent to those skilled in the art that other variations and modifications of the present invention can be made based on the above description, and it is not intended to be exhaustive or to limit the invention to the precise form disclosed, and all such modifications and variations are possible and contemplated as falling within the scope of the invention.
Claims (10)
1. A backstepping control method based on a saturation adaptive law is characterized in that a two-dimensional nonlinear system state space model is established for a controlled object, and the controlled object is controlled by the saturation adaptive law;
the design process of the saturation adaptive law comprises the following steps:
the method comprises the following steps of firstly, establishing a two-dimensional nonlinear system state space model aiming at a controlled object, wherein the established state space model of the two-dimensional nonlinear system is in the following specific form:
y=x1
wherein x is1,x2Represents the state variable of the system and is,representsSystem state variable x1,x2The first derivative of (a), y is system output, f is a nonlinear smooth function unknown to the system, d is unknown disturbance outside the system, u represents a system control input signal, and the control aim is to design the system control input u so that the system output y tracks a given target signal;
step two, according to the state variable x of the system1And a target signal ydAnd the state variable x of the system2And a virtual control function alpha to be designed1Separately constructing the error variable z1And z2;
Step three, utilizing the error variable z obtained in the step two1And z2Designing a Lyapunov function V;
step four, solving the first derivative of the Lyapunov function V in the step three to obtain the first derivative
Step five, according to the first derivative of the Lyapunov functionDesigning a virtual control function alpha1And based on a virtual control function alpha1Design control input u:
wherein the content of the first and second substances,representing a virtual control function alpha1The first derivative of (a) is,for estimation of an unknown non-linear function f, k2Is a constant greater than zero;the function is a saturation function;
2. The saturation adaptive law-based backstepping control method according to claim 1, wherein the controlled object is a motor system, and the corresponding system state variable x1,x2The rotational angle and the rotational speed of the motor system,the rotating speed and the angular acceleration of the motor system are shown, and the system output y is the rotating angle of the motor system.
3. A method for backstepping control based on saturation adaptive law according to claim 1 or 2, characterized in that the solution of the state space model (1) of the two-dimensional non-linear system exists and is unique; the unknown nonlinear smooth function f and the first-order second-order derivative thereof satisfyWhereinIs a constant greater than zero; target signal ydAnd its first and second derivatives are bounded.
9. A storage medium having stored therein at least one instruction that is loaded and executed by a processor to implement a method of saturation adaptive law based backstepping control as claimed in any one of claims 1 to 8.
10. An apparatus comprising a processor and a memory, the memory having stored therein at least one instruction that is loaded and executed by the processor to implement a method of saturation adaptive law based backstepping control as claimed in any one of claims 1 to 8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111405311.2A CN114035436B (en) | 2021-11-24 | 2021-11-24 | Backstepping control method, storage medium and device based on saturation adaptive law |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111405311.2A CN114035436B (en) | 2021-11-24 | 2021-11-24 | Backstepping control method, storage medium and device based on saturation adaptive law |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114035436A true CN114035436A (en) | 2022-02-11 |
CN114035436B CN114035436B (en) | 2024-04-02 |
Family
ID=80138687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111405311.2A Active CN114035436B (en) | 2021-11-24 | 2021-11-24 | Backstepping control method, storage medium and device based on saturation adaptive law |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114035436B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009228605A (en) * | 2008-03-24 | 2009-10-08 | Daihatsu Motor Co Ltd | Method for controlling engine speed and device for controlling engine speed |
WO2019165801A1 (en) * | 2018-03-02 | 2019-09-06 | 曾喆昭 | Disturbance perception control method |
GB201911738D0 (en) * | 2019-07-12 | 2019-10-02 | Huaiyin Inst Technology | Adaptive backstepping optimal control method of fractional-order chaotic electromechanical transducer system |
CN110647104A (en) * | 2019-09-02 | 2020-01-03 | 华南理工大学 | Flexible riser backstepping boundary control method based on boundary disturbance observer |
CN113093553A (en) * | 2021-04-13 | 2021-07-09 | 哈尔滨工业大学 | Adaptive backstepping control method based on instruction filtering disturbance estimation |
-
2021
- 2021-11-24 CN CN202111405311.2A patent/CN114035436B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009228605A (en) * | 2008-03-24 | 2009-10-08 | Daihatsu Motor Co Ltd | Method for controlling engine speed and device for controlling engine speed |
WO2019165801A1 (en) * | 2018-03-02 | 2019-09-06 | 曾喆昭 | Disturbance perception control method |
GB201911738D0 (en) * | 2019-07-12 | 2019-10-02 | Huaiyin Inst Technology | Adaptive backstepping optimal control method of fractional-order chaotic electromechanical transducer system |
CN110647104A (en) * | 2019-09-02 | 2020-01-03 | 华南理工大学 | Flexible riser backstepping boundary control method based on boundary disturbance observer |
CN113093553A (en) * | 2021-04-13 | 2021-07-09 | 哈尔滨工业大学 | Adaptive backstepping control method based on instruction filtering disturbance estimation |
Non-Patent Citations (1)
Title |
---|
XIAOLONG ZHENG等: "Command Filter and Universal Approximator Based Backstepping Control Design for Strict-Feedback Nonlinear Systems With Uncertainty", 《IEEE TRANSACTIONS ON AUTOMATIC CONTROL》, vol. 65, no. 3, 27 February 2020 (2020-02-27) * |
Also Published As
Publication number | Publication date |
---|---|
CN114035436B (en) | 2024-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107121932B (en) | Motor servo system error symbol integral robust self-adaptive control method | |
CN110673472B (en) | Adaptive robust control method based on neural network compensation dead zone inversion error | |
CN107870570B (en) | Terminal sliding mode mechanical arm track tracking method based on fractional order power approach law | |
CN112817231B (en) | High-precision tracking control method for mechanical arm with high robustness | |
CN104333280B (en) | Robustness adaptive control (RAC) method of direct driving motor system | |
CN108228975B (en) | Motor servo system parameter identification method and backlash elimination control method | |
CN107544256A (en) | Underwater robot sliding-mode control based on adaptive Backstepping | |
CN107561935A (en) | Motor position servo system friciton compensation control method based on multilayer neural network | |
CN106774379B (en) | Intelligent supercoiled strong robust attitude control method | |
Hung et al. | Adaptive control for nonlinearly parameterized uncertainties in robot manipulators | |
CN113589689A (en) | Sliding mode controller design method based on multi-parameter adaptive neural network | |
CN111965976B (en) | Robot joint sliding mode control method and system based on neural network observer | |
CN114815618B (en) | Adaptive neural network tracking control method based on dynamic gain | |
Hu et al. | Precision motion control of a small launching platform with disturbance compensation using neural networks | |
CN111736472A (en) | Motor self-adaptive preset performance asymptotic control method based on RISE | |
CN109108964B (en) | Spatial mechanical arm coordination control method based on adaptive dynamic programming Nash game | |
CN112904726B (en) | Neural network backstepping control method based on error reconstruction weight updating | |
Wu et al. | Characteristic modeling and control of servo systems with backlash and friction | |
CN109995278B (en) | Motor servo system self-adjustment control method considering input limitation | |
CN114035436A (en) | Backstepping control method based on saturation adaptive law, storage medium and equipment | |
Li et al. | Repetitive learning control of nonlinear continuous-time systems using quasi-sliding mode | |
CN115720061A (en) | Fuzzy self-adaptive backstepping control method of electromechanical servo system based on finite time | |
CN114859725A (en) | Self-adaptive event trigger control method and system for nonlinear system | |
Zhang et al. | Adaptive backstepping control of dual-motor driving servo systems with friction | |
CN111856941B (en) | Adaptive terminal dynamic sliding mode control method based on active disturbance rejection |
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 |