CN112859587B - PID target tracking control method based on additional integrated module - Google Patents

PID target tracking control method based on additional integrated module Download PDF

Info

Publication number
CN112859587B
CN112859587B CN202110044826.8A CN202110044826A CN112859587B CN 112859587 B CN112859587 B CN 112859587B CN 202110044826 A CN202110044826 A CN 202110044826A CN 112859587 B CN112859587 B CN 112859587B
Authority
CN
China
Prior art keywords
tracking
pid
integrated module
control method
target
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
CN202110044826.8A
Other languages
Chinese (zh)
Other versions
CN112859587A (en
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.)
Institute of Optics and Electronics of CAS
Original Assignee
Institute of Optics and Electronics 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 Institute of Optics and Electronics of CAS filed Critical Institute of Optics and Electronics of CAS
Priority to CN202110044826.8A priority Critical patent/CN112859587B/en
Publication of CN112859587A publication Critical patent/CN112859587A/en
Application granted granted Critical
Publication of CN112859587B publication Critical patent/CN112859587B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P.I., P.I.D.
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Abstract

The invention provides a PID target tracking control method based on an additional integrated module. In the field of target tracking control, since a target always moves at a low intermediate frequency, the tracking accuracy of the low intermediate frequency is a problem of general attention of the system. The method comprises the steps of firstly identifying a controlled object on a photoelectric tracking experiment platform, then designing a PID controller based on an additional integrated module, solving the problem of parameter adjustment in the controller by adopting a multi-target optimization algorithm, and finally completing the tracking of a target after implementing the control method provided by the invention. The performance index representing the tracking accuracy of the photoelectric tracking system is known to be error suppression capability, and after comparing a PID controller based on an additional integrated module with a PID control method without adding the integrated module, the error suppression capability of the control method provided by the invention at low and medium frequencies is obviously superior to that of a PID, thereby achieving the effect of improving the tracking accuracy.

Description

PID target tracking control method based on additional integrated module
Technical Field
The invention belongs to the field of tracking control of a photoelectric system, and particularly relates to a PID (proportion integration differentiation) controller based on an additional integrated module, which is used for target tracking. The additional integrated module is added into the traditional PID controller to control the voltage signal of the motor, so that the error suppression capability of the system is effectively improved, and the tracking precision of the photoelectric tracking system is fully improved.
Background
The photoelectric tracking system consists of a main control system, a servo system, an image system, a signal transmission and equal division system, and the servo system drives the motor to rotate according to the output signal of the controller and the voltage after receiving the deviation signal of the image sensor, so that the tracking error is reduced. The design of the controller in the servo system is therefore an important factor affecting the tracking accuracy. In most cases, the system is always designed as a type-I system, that is, the open-loop transfer function of the system only contains one integral element, and at the moment, the system has a certain low-frequency error suppression capability. However, the tracking capabilities of this type of system are somewhat inadequate when performing fast object tracking tasks. Therefore, the document PID-I controller of charged coupled tracking loop for fast-steering mirror [ J ]. Optical Engineering,2011,50(4):043002. the PID-I controller is proposed to enable the system to have two integration links, thereby improving the non-static-error degree of the system; in addition, the adjustment parameters of the controller are considered from an amplitude margin and a phase margin, namely the amplitude margin is larger than or equal to 6dB, and the phase margin is larger than or equal to 45 degrees, so that the error suppression effect of 10dB improvement compared with the first-type system under the sampling frequency of 2KHz is obtained. However, the document does not consider the problem that the stability is reduced due to the phase lag because the system type is improved. Based on this, the document "A comprehensive performance improvement control method by fractional order control" [ J ]. IEEE Photonics Journal,2018,10(5):1-11. the controller is subjected to fractional order exploration, namely the type of the system is between one type and two types, but not integral multiple of the number of integral links; the controller has the advantage of higher error rejection than the one-type system, but has less phase lag than the two-type system, so the method adopts a compromise method to improve the tracking accuracy of the system. From the basic theory of signal and information processing, it can be seen that any complex signal is compounded by a plurality of sine waves with different frequencies and different sizes. From the control theory, a higher-type control loop has the capability of suppressing higher-order tracking errors, so that a high-type system always belongs to a research hotspot in order to reduce the tracking errors of the system as much as possible. Meanwhile, the disadvantage of the high-type loop is an unstable phenomenon caused by phase lag, and parameter adjustment work of the controller is still important in order to better apply the advantages of the control algorithm.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: in order to improve the tracking precision of a tracking system, the invention provides a PID target tracking control method based on an additional integrated module. The method adds an integrated module on the basis of the original PID controller, and aims to improve the type of a control loop to three types, thereby improving the tracking accuracy of a tracking system. Meanwhile, the problem of parameter adjustment of the controller based on the multi-target optimization algorithm NSGA-II is also solved.
The technical scheme adopted by the invention for solving the technical problems is as follows: a PID target tracking control method based on an additional integrated module comprises the following specific steps:
step (1), installing an image sensor (CCD) on a tracking platform, and testing the position characteristics of a controlled object by a Dynamic Signal Analyzer (DSA), thereby obtaining a position object model G p (s);
Step (2) at the acquisition positionControlled object model G p On the basis of(s), a Proportional-Integral-Derivative (PID) controller C is designed P (s), wherein the proportional link exists to reduce the tracking error of the system, and the integral link exists to improve the no-static-error degree of the system; the differential link exists to improve the dynamic performance of the system;
step (3) in the designed PID controller C P Adding an additional integrated module C on the basis of(s) a (s) forming a PID controller C based on the additional integrated module m (s) wherein the module C is integrated a (s) the purpose is to upgrade the type of system to three, thus upgrading the system to be free of slack;
step (4) aiming at PID controller C based on additional integrated module m The undetermined parameters of(s) are obtained by optimizing a multi-target optimizing Algorithm, namely a Non-poor ordering Genetic Algorithm (NSGA-II), wherein the performance indexes are as follows: system bandwidth, system error rejection ratio and stability margin interval;
step (5), designing a PID controller C based on an additional integrated module m (s), the tracking method is applied to an actual system to generate a voltage signal to control a motor, so that the target is tracked; after comparing with the PID control method without adding an integrated module, the effect of improving the tracking precision is achieved.
By utilizing the five steps, the target can be tracked through the PID controller based on the additional integrated module. Firstly, the image sensor detects a system tracking error signal, secondly, the error signal is input into the position controller to obtain a voltage signal, and the motor is controlled, so that the tracking platform rotates, the tracking error is gradually reduced, and the tracking task is completed. Finally, by comparison, the tracking accuracy of the control method provided by the invention at low and medium frequency (target motion frequency band) is obviously superior to that of PID, and the effect of improving the tracking accuracy is achieved.
Compared with the prior art, the invention has the following advantages:
(1) compared with the traditional one-type loop control method, the tracking controller provided by the invention has higher tracking precision and has the advantage of tracking a faster target;
(2) compared with the common control method for improving the tracking precision, such as sensor signal fusion, an error observer and the like, the tracking controller provided by the invention does not need to introduce an additional control loop or sensor;
(3) compared with the traditional method for adjusting the parameters of the controller by using the numerical values, the method provided by the invention fully considers the principle of conservation of system performance, considers the problems from the perspective of multi-performance indexes, and is more comprehensive;
(4) the invention has low cost and small calculation amount, fully utilizes the conditions of the system and does not add extra cost, space, load and power consumption to the system.
Drawings
FIG. 1 is a control block diagram of the present invention;
FIG. 2 is a physical schematic of the platform drive motor and load;
FIG. 3 is a diagram of a system position model fit;
FIG. 4 is a closed loop Bode diagram of the system position;
the system of fig. 5 uses pure PID versus Bode plot for error rejection based on additional integration module PID.
Detailed Description
Embodiments of the present invention are described below with reference to the drawings. The following examples are only illustrative of the present invention, and the scope of the present invention shall include the full contents defined by the claims; and the claims of the present invention can be realized in their entirety by those skilled in the art from the following examples.
The invention relates to a PID target tracking control method based on an additional integrated module, which specifically comprises the following steps:
(1) the control block diagram of the present invention is shown in fig. 1. The block diagram uses the output signal y to perform a tracking function on the input signal r. Wherein G(s) is a position controlled object measured by the image sensor, C(s) is a controller for outputting a voltage signal u, e is a tracking error detected by the image sensor, F(s) is an open-loop transfer function, and T(s) is a closed-loop transfer function.
Taking the tracking process of the photoelectric tracking experiment platform driven by the voice coil motor in a single direction as an example, the model is simplified and analyzed, as shown in fig. 2. The circuit balance principle can be used as follows:
Figure BDA0002896762020000031
wherein U is a Is the armature voltage of the voice coil motor, K b Is the motor back emf coefficient, θ a Being the deflection angle of the platform, R a Is the motor armature loop resistance, I a Is the motor armature loop current, L a Is the motor armature loop inductance, and t is time.
The motor moment balance equation can be obtained:
Figure BDA0002896762020000041
wherein C is m Is the force coefficient of the voice coil motor, I a For the motor armature loop current, d is the distance from the elastic support to the center of the voice coil motor, J is the load moment of inertia of the platform, and theta a Angle of deflection of the platform, f m Is the mechanical damping coefficient of the platform, K is the elastic support elastic coefficient, and t is the time.
The circuit balance principle and the motor moment balance equation are combined to eliminate I a (t) and the driving input voltage U is obtained by Laplace change a And an output deflection angle theta a The transfer function of (c):
Figure BDA0002896762020000042
where s is the laplacian operator.
Therefore, the transfer model of the platform deflection angle and the input voltage is presented as an approximate third-order filter link. Since the numerator is constant, at a low frequency band, the input voltage and the output angle can be regarded as a proportional characteristic, which also indicates that the input voltage corresponds to the output angle position in a physical dimension. According to the polynomial theory, a third-order polynomial of a denominator always has a real number, so that a transfer function of a platform can be decomposed into series connection of a proportional link, an inertia link and a second-order oscillation link.
Figure BDA0002896762020000043
Wherein K is the coefficient of inertia element, T e Is an inertia coefficient of a first-order inertia link,
Figure BDA0002896762020000044
the second-order oscillation link natural frequency is obtained, xi is the second-order oscillation link damping ratio, and s is a Laplace operator.
(2) Fig. 3 shows a graph of the position transfer function Bode of the system obtained from the system sweeping the measured values. The solid line represents the measured curve, and the dotted line is the object obtained after identification according to the typical link.
After the object identification of fig. 3, the position loop PID controller is designed by using the zero-pole cancellation principle:
Figure BDA0002896762020000045
wherein the proportional link parameters
Figure BDA0002896762020000046
Differential link parameter
Figure BDA0002896762020000047
The undetermined parameter is an integral link parameter k i And T c
Controller molecule part
Figure BDA0002896762020000048
The second-order oscillation link in the denominator of the controlled position object is used for zero-pole cancellation and the inertia link
Figure BDA0002896762020000051
Is a low-pass filter, in order to ensure controlSystem ware C p Causality of(s) and at the same time for filtering out high frequency noise.
(3) Designing PID controller C based on additional integrated module m (s) driving the motor.
Integrated module C a (s) and PID controller C based on additional integrated module m (s) are respectively represented as:
Figure BDA0002896762020000052
Figure BDA0002896762020000053
wherein k, a and b respectively represent the proportional element and the coefficients of two first-order differential elements in the integrated module, and k c =kk i . Due to T c Is a low-pass filter which is always set near one-half sampling frequency point, the sampling frequency is 1 KHz, and T is set in the experiment c Set to 0.00032. The undetermined parameter is therefore k c ,a,b。
(4) Method for solving undetermined parameter k through NSGA-II multi-target optimization algorithm c ,a,b。
Wherein, the performance indexes are respectively:
Figure BDA0002896762020000054
wherein, J 1 ,J 2 ,J 3 Respectively representing the bandwidth, the error rejection ratio and the performance index of a stability margin interval, omega B Represents the frequency point, omega, corresponding to the amplitude value of-3 dB in the closed loop Bode diagram of the system 0 And ω c Respectively representing the initial frequency and the open-loop cut-off frequency, | E of the system m (j ω) | is the absolute value of the amplitude of the system frequency domain error suppression function in the Bode graph amplitude-frequency curve, and s ═ j ω; omega b And omega a Respectively represents two frequency points corresponding to the phase margin of 45 degrees in the open loop Bode graph phase-frequency curve of the system, and omega b >ω a . Selecting a set of parameters with the highest error suppression ratio, namely k, in the final optimization result c =5990,a=1.2,b=0.03。
(5) The closed loop characteristics were measured by applying the controller to the actual system, and a Bode plot of the closed loop position of the system is shown in fig. 4, which shows the tracking bandwidth around 5 hz.
(6) Pure PID controller C is adopted for comparison p (s) and PID controller C based on additional integrated module m Error suppression capability of(s), fig. 5 shows two systematic error suppression characteristics compared to Bode plots.
By comparison, it can be seen that the tracking error suppression curve is improved below 3.3 hz (low-if) by using the method of the present invention. Specifically, there is a boost effect of around 10dB at 1 hz. Because the target always moves at a low intermediate frequency, the method can effectively improve the tracking accuracy of the tracking system and has the potential of tracking a faster target.
And (4) completing the control method for realizing target tracking through the PID controller based on the additional integrated module by utilizing the six steps. The tracking accuracy of the tracking system is improved under the condition that a control loop is not added and an additional sensor is not introduced.
Parts of the invention not described in detail are well known in the art.

Claims (4)

1. A PID target tracking control method based on an additional integrated module is characterized in that: the method comprises the following steps:
step (1), installing an image sensor on a tracking platform, and carrying out CCD: Charge-Coupled Device, measured by dynamic signal analyzer, DSA: a Dynamic Signal Analyzer for testing the position characteristics of the controlled object so as to obtain a position controlled object model G p (s);
Step (2) of obtaining the controlled object model G at the position p On the basis of(s), a Proportional-Integral-Derivative (PID) controller C is designed P (s), wherein the proportional element is used for reducing the tracking error of the system, and the integral element is used for improving the no-static-error of the systemDegree; the differential link exists to improve the dynamic performance of the system;
step (3) in the designed PID controller C P Adding an additional integrated module C on the basis of(s) a (s) forming a PID controller C based on the additional integrated module m (s) wherein an integrated module C is added a (s) the purpose is to upgrade the type of system to three, thus upgrading the system to be free of slack;
for additional integrated modules C in step (3) a (s) and PID controller C based on additional integrated module m (s), the expression is as follows:
Figure FDA0003651181590000011
Figure FDA0003651181590000012
in the above formula, k, a and b are each C a (s) the coefficient of the proportional element, the two first order differential elements, and k c =kk i From a PID controller C based on an add-on integrated module m (s) expression, it can be seen that P (s) two integral operators, two first order differential links and a proportion link are introduced on the basis, wherein the integral link is used for improving the tracking precision of the tracking system, but 180-degree phase lag is introduced for the system, so that the two first order differential links are used for improving the phase of the system, and the phenomenon of system instability caused by too low open loop stability margin is avoided; the proportion link exists for finally adjusting the open loop gain of the open loop system;
step (4) aiming at PID controller C based on additional integrated module m The undetermined parameter(s) is obtained by optimizing a multi-target optimizing Algorithm, namely a Non-poor ordering Genetic Algorithm NSGA-II, a Non-Dominated sequencing Genetic Algorithm II, wherein the performance indexes are as follows: system bandwidth, system error rejection ratio and stability margin interval;
step (5),By designing PID controller C based on additional integrated module m (s), the tracking method is applied to an actual system to generate a voltage signal to control a motor, so that the target is tracked; after comparing with the PID control method without adding an integrated module, the effect of improving the tracking precision is achieved.
2. The PID target tracking control method based on the additional integration module as claimed in claim 1, wherein: when DSA is adopted to measure the controlled object in the step (1), the input signal is a sine sweep voltage signal, and the output signal is a position signal of an image sensor, so that a position controlled object model G is obtained p (s); secondly, offline identification is carried out on the actually measured curve, the type of the curve belonging to a typical link is preliminarily judged according to the trend of the curve, and then the size of parameters in the typical link is finely adjusted.
3. The PID target tracking control method based on the additional integrated module according to claim 1, wherein: the PID controller expression in the step (2) is as follows:
Figure FDA0003651181590000021
in the above formula, k p ,k i ,k d The coefficients of the proportional, integral and differential links, T c Is to make C P (s) satisfy causal inertial element coefficients, s being the frequency domain "laplacian" operator.
4. The PID target tracking control method based on the additional integration module as claimed in claim 1, wherein: in the step (4), because the performance of the system is conserved and each performance index of the actual system is sometimes mutually exclusive, based on the bandwidth and the error rejection ratio, the stability margin interval is used as a measure index for solving the controller parameter by the NSGA-II, and the performance index expression is as follows:
Figure FDA0003651181590000022
wherein, J 1 ,J 2 ,J 3 Respectively representing the performance indexes of bandwidth, error rejection ratio and stability margin interval, omega B Represents the frequency point, omega, corresponding to the amplitude value of-3 dB in the closed loop Bode diagram of the system 0 And ω c Respectively representing the initial frequency and the open-loop cut-off frequency, | E of the system m (j ω) | is the absolute value of the amplitude of the system frequency domain error suppression function in the Bode graph amplitude-frequency curve, and s ═ j ω; omega b And omega a Respectively represents two frequency points corresponding to the phase margin of 45 degrees in the open loop Bode graph phase-frequency curve of the system, and omega ba
By utilizing the steps (1) to (5), the target can be tracked through a PID controller based on an additional integrated module, firstly, an image sensor detects a system tracking error signal, secondly, the error signal is input into the position controller to obtain a voltage signal, and a motor is controlled, so that a tracking platform rotates, the tracking error is gradually reduced, and a tracking task is completed; finally, by comparison, the tracking accuracy of the control method at low and medium frequencies, namely the target motion frequency band, is obviously superior to that of a PID (proportion integration differentiation), and the effect of improving the tracking accuracy is achieved.
CN202110044826.8A 2021-01-13 2021-01-13 PID target tracking control method based on additional integrated module Active CN112859587B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110044826.8A CN112859587B (en) 2021-01-13 2021-01-13 PID target tracking control method based on additional integrated module

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110044826.8A CN112859587B (en) 2021-01-13 2021-01-13 PID target tracking control method based on additional integrated module

Publications (2)

Publication Number Publication Date
CN112859587A CN112859587A (en) 2021-05-28
CN112859587B true CN112859587B (en) 2022-08-23

Family

ID=76003536

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110044826.8A Active CN112859587B (en) 2021-01-13 2021-01-13 PID target tracking control method based on additional integrated module

Country Status (1)

Country Link
CN (1) CN112859587B (en)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8072204B1 (en) * 2007-05-09 2011-12-06 Zilker Labs, Inc. Control system optimization via digital diode emulation
CN105045141A (en) * 2015-05-27 2015-11-11 中国科学院光电技术研究所 Analog control circuit capable of enlarging control bandwidth of fast steering mirror
CN105119471B (en) * 2015-08-01 2017-07-21 安徽铜冠铜箔有限公司 A kind of circuit and method for eliminating electrolytic copper foil surface high frequency switch power vortex
CN105425580B (en) * 2015-12-14 2018-02-06 北京理工大学 A kind of ITAE optimal Ns type system constituting method
US10009121B2 (en) * 2016-04-29 2018-06-26 Commscope Technologies Llc Testing methods using spread-spectrum excitation signals and related test apparatus
CN108681242B (en) * 2018-05-02 2021-08-13 中国科学院光电技术研究所 Fractional order control method for fast reflector in photoelectric system
CN108646572B (en) * 2018-07-16 2020-12-25 广西师范大学 Control method of three-axis pan-tilt servo motor based on combination of BP neural network and active disturbance rejection controller
CN109245534B (en) * 2018-10-23 2020-04-17 电子科技大学 Digital proportional-integral-differential compensation circuit based on double signal paths
CN109884882A (en) * 2019-02-25 2019-06-14 中国科学院光电技术研究所 A kind of photoelectric follow-up control method based on differential tracker
CN110262222B (en) * 2019-06-26 2021-07-13 东北大学 Optimal interval PID controller parameter optimization method for interval system

Also Published As

Publication number Publication date
CN112859587A (en) 2021-05-28

Similar Documents

Publication Publication Date Title
CN101546173B (en) Method and apparatus for controlling system
CN101454969B (en) Servo controller
CN110784149B (en) Mechanical resonance suppression method and system for alternating current servo system
US7706899B2 (en) Method and apparatus for adaptive cascade proportional-integral-derivative controller
CN104242769A (en) Permanent magnet synchronous motor speed composite control method based on continuous terminal slip form technology
Sira-Ramírez et al. On the equivalence between ADRC and Flat Filter based controllers: A frequency domain approach
CN111817638A (en) Phase advance linear active disturbance rejection controller of permanent magnet synchronous linear motor platform
CN111258220A (en) Flexible mechanical arm cascade control method and system based on disturbance observer
CN110209122B (en) Control method, device, medium and equipment of multi-axis motion platform
WO2018173654A1 (en) Method for designing filter of delay compensator, feedback control method using same, and motor control device
JP3266931B2 (en) Motor control device
CN112859587B (en) PID target tracking control method based on additional integrated module
CN112334845A (en) Feedback control method and feedback control device
US10558179B2 (en) Shaft torque control device
JP5017984B2 (en) Servo control device and speed tracking control method thereof
Garba et al. Crone controller based speed control of permanent magnet direct current motor
CN116805849A (en) Continuous set model prediction control method of permanent magnet synchronous motor
CN113765453B (en) Suspension control system of magnetic suspension switch reluctance motor with wide-narrow pole characteristics
JP2011078192A (en) Controller for motor
CN111796509A (en) Gyro self-stabilization control method
CN111208728A (en) Linear active disturbance rejection control method, device equipment and storage medium
CN113517835A (en) PMSM drive system field loss fault control method and permanent magnet synchronous motor
Casavola et al. H∞ and H2 simple controllers for robotic applications
CN116520684B (en) Controller optimization method based on active disturbance rejection control and Youla parameterization
CN116610035B (en) Completely decoupled active disturbance rejection control system and parameter setting method thereof

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