CN113703464A - 2D simulated fish heading control method - Google Patents

2D simulated fish heading control method Download PDF

Info

Publication number
CN113703464A
CN113703464A CN202110885848.7A CN202110885848A CN113703464A CN 113703464 A CN113703464 A CN 113703464A CN 202110885848 A CN202110885848 A CN 202110885848A CN 113703464 A CN113703464 A CN 113703464A
Authority
CN
China
Prior art keywords
point
coordinate point
fish
degrees
tangent
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
CN202110885848.7A
Other languages
Chinese (zh)
Other versions
CN113703464B (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.)
Army Engineering University of PLA
Original Assignee
Army Engineering University of PLA
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 Army Engineering University of PLA filed Critical Army Engineering University of PLA
Priority to CN202110885848.7A priority Critical patent/CN113703464B/en
Publication of CN113703464A publication Critical patent/CN113703464A/en
Application granted granted Critical
Publication of CN113703464B publication Critical patent/CN113703464B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/04Control of altitude or depth
    • G05D1/06Rate of change of altitude or depth
    • G05D1/0692Rate of change of altitude or depth specially adapted for under-water vehicles

Abstract

The invention discloses a 2D simulation fish heading control method, which comprises the following steps: calculating an angle deviation threshold value caused by a water wave random interference function according to the ideal ball hitting point, the fish head and the coordinate point of the target position; selecting a heading strategy according to the angle deviation threshold, and performing heading action according to the heading strategy when the simulated fish moves to the target position; and judging whether the ball pushing action is finished or not according to the distance between the center of the ball and the target position, if so, finishing, and if not, circularly executing the steps. The invention can construct the entity environment by means of analog simulation, and can obtain the research and development target with controllable cost.

Description

2D simulated fish heading control method
Technical Field
The invention relates to a top ball control method of a 2D (two-dimensional) simulated fish, belonging to the technical field of simulated fish.
Background
With the continuous deepening of the aquatic strategy in China, the aquatic bionic technology gets more and more attention, the aquatic bionic robot realizes the functions of underwater target detection and identification, underwater navigation (positioning), water surface battlefield simulation, energy collection and the like on the basis of a cybernetics and a system theory, and practical application products of the aquatic bionic robot are more and more extensive in the field of underwater science and technology application. However, the underwater bionic robot is expensive and has high research and development cost; in order to solve the problems, the application provides a 2D simulated fish heading control method.
Disclosure of Invention
The invention aims to overcome the defects in the prior art and provide a 2D simulated fish heading control method, which is used for constructing a physical environment by means of simulation and can obtain a research and development target with controllable cost.
In order to achieve the purpose, the invention is realized by adopting the following technical scheme:
the invention provides a 2D (two-dimensional) simulation fish heading control method, which comprises the following steps:
calculating an angle deviation threshold value caused by a water wave random interference function according to the ideal ball hitting point, the fish head and the coordinate point of the target position;
selecting a heading strategy according to the angle deviation threshold, and performing heading action according to the heading strategy when the simulated fish moves to the target position;
judging whether the ball pushing action is finished or not according to the distance between the center of the ball and the target position, if so, finishing, and if not, circularly executing the steps;
and the coordinate point of the ideal hitting point is obtained according to the fish head of the simulated fish, the sphere center of the sphere and the initialization coordinate point of the target position.
Preferably, the obtaining of the coordinate point of the ideal hitting point includes:
making a straight line based on the fish head coordinate point and the spherical center coordinate point;
if the coordinate point of the target position is on one side of the straight line, two symmetrical tangent lines are made to the collision circle of the sphere through the coordinate point of the fish head, two tangent points are obtained, and the tangent point on the other side of the straight line is taken as an ideal hitting point;
if the coordinate point of the target position is on the straight line, two intersection points of the straight line and a collision circle of the sphere are obtained, and the intersection point far away from the target point is taken as an ideal hitting point;
wherein the collision circle of the sphere is a central section of the sphere appearing on 2D.
Preferably, the calculating the angle deviation threshold caused by the water wave random interference function includes:
making a triangle according to the ideal hitting point, the fish head and the coordinate point of the target position, and solving a first included angle value with the ideal hitting point as a vertex according to the cosine law;
making a triangle according to the ideal hitting point, the fish head offset and the coordinate point of the target position under the water wave random interference function, and solving a second included angle value taking the ideal hitting point as a vertex according to the cosine law;
and obtaining an angle deviation threshold value according to the absolute value of the difference value of the angle value of the first included angle and the angle value of the second included angle.
Preferably, the selecting a heading strategy according to the angle deviation threshold includes:
if the angle deviation threshold is less than or equal to 30 degrees, adopting a linear ball jacking algorithm as a ball jacking strategy;
if the angle deviation threshold value is larger than 30 degrees and smaller than or equal to 180 degrees, an arc cut-in dome ball algorithm is adopted as the top ball strategy.
Preferably, the straight line heading algorithm includes:
when the angle value of the first included angle is less than 5 degrees, the direction of the simulated fish is unchanged, and the speed is 14 grades of the preset grade;
when the angle value of the first included angle is more than or equal to 5 degrees and less than 20 degrees, the direction of the simulated fish turns at a preset gear 6-5 or 8-9, and the speed is at a preset level of 12;
when the angle value of the first included angle is more than or equal to 20 degrees and less than 40 degrees, the direction of the simulated fish turns at a preset gear 4-3 or 10-11, and the speed is at a preset level of 8;
when the angle value of the first included angle is more than or equal to 40 degrees and less than 60 degrees, the direction of the simulated fish turns at a preset gear 2-1 or 12-13, and the speed is at a preset level 2;
when the angle value of the first included angle is more than or equal to 60 degrees and less than or equal to 180 degrees, the direction of the simulated fish turns at a preset gear 0 or 14, and the speed is at a preset grade 0 level;
the preset speed level is 0-14 levels, and the preset speed levels are 15 levels and are sequentially increased; the preset gears in the direction are 0-14 gears and 15 gears in total, 0-6 gears are left-turning gears and are sequentially decreased in number, 7 gears are unchanged in direction, and 8-14 gears are right-turning gears and are sequentially increased in number.
Preferably, the arc plunge dome algorithm comprises:
simulating a coordinate point of the tangential point towards which the fish is oriented;
wherein the obtaining of the coordinate point of the tangent point comprises:
making a tangent line to the incision circle of the sphere through the coordinate point of the fish head and obtaining a tangent point;
making a triangle according to the coordinate point of the fish head, the circle center of the incision circle and the tangent point, and solving the tangent length according to the pythagorean theorem;
calculating to obtain an arc deflection included angle according to the radius of the cut-in circle and the tangent length;
acquiring a coordinate point of the rotation center of the simulated fish, and calculating a direction vector by combining the coordinate point of the fish head;
calculating to obtain a rotation vector according to the direction vector and the arc deflection included angle;
calculating to obtain a coordinate point of the tangent point according to the rotation vector and the tangent length;
wherein the incised circle of the sphere is a central section of the sphere appearing in 2D.
Compared with the prior art, the invention has the following beneficial effects:
according to the method for controlling the top ball of the 2D simulated fish, the position of the simulated fish can be effectively and reasonably adjusted in each cycle period, the intelligent top ball which is in opposition to a water wave random interference function is realized, the speed and the stability of the top ball are improved, and the time for the top ball is shortened under the condition of deviation; the entity environment is constructed by means of simulation, and a research and development target with controllable cost can be obtained.
Drawings
FIG. 1 is a flow chart of a method for controlling a 2D simulated fish heading according to an embodiment of the invention;
FIG. 2 is a first schematic diagram of obtaining an ideal hitting point according to the embodiment of the present invention;
FIG. 3 is a schematic diagram of obtaining an ideal hitting point according to the embodiment of the present invention;
FIG. 4 is a schematic diagram of a third exemplary embodiment of an ideal hitting point acquisition;
fig. 5 is a schematic diagram of obtaining an angle deviation threshold according to an embodiment of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
The first embodiment is as follows:
the embodiment performs simulation based on URWPGSim2D platform:
as shown in fig. 1, the present embodiment provides a method for controlling a heading of a 2D simulated fish, including the following steps:
step 1, obtaining a coordinate point of an ideal hitting point T according to an initialization coordinate point of a fish head F of the simulated fish, a sphere center C of the sphere and a target position S; the method specifically comprises the following steps:
step 1.1, making a straight line based on a coordinate point of a fish head F and a coordinate point of a sphere center C;
step 1.2, if the coordinate point of the target position S is on one side of the straight line, making two symmetrical tangent lines to the collision circle of the sphere through the coordinate point of the fish head F, obtaining two tangent points, and taking the tangent point on the other side of the straight line as an ideal batting point T; as shown in fig. 2;
step 1.3, if the coordinate point of the target position S is on a straight line, acquiring two intersection points of a collision circle of the straight line and a ball body, and taking the intersection point far away from the target point as an ideal hitting point T; as shown in fig. 3-4;
wherein the collision circle of the sphere is a central section of the sphere appearing on 2D.
Step 2, calculating an angle deviation threshold value caused by a water wave random interference function according to the ideal hitting point, the fish head and the coordinate point of the target position; the method specifically comprises the following steps:
step 2.1, as shown in fig. 5, making a triangle according to coordinate points of an ideal hitting point T, a fish head F and a target position S, and solving a first included angle value & lt FTS with the ideal hitting point T as a vertex according to a cosine law;
Figure BDA0003194080820000051
acquiring an angle value of the FTS, wherein the FTS is more than or equal to 0 degree and less than or equal to 180 degrees;
step 2.2, according to a water wave random interference function, making a triangle according to an ideal hitting point T, a fish head offset F 'and a coordinate point of a target position S, and solving a second included angle value & lt F' TS with the ideal hitting point T as a vertex according to a cosine theorem;
Figure BDA0003194080820000052
acquiring an angle value of an angle F 'TS, wherein the angle value of the angle F' TS is more than or equal to 0 degree and less than or equal to 180 degrees;
and 2.3, obtaining an angle deviation threshold value xzdangle according to the difference value of the angle value of the first included angle and the angle value of the second included angle, wherein xzdangle is | < F' TS- < FTS |.
Step 3, selecting a heading strategy according to the angle deviation threshold, and performing heading action according to the heading strategy when the simulated fish moves to the target position;
the angle deviation threshold value selection heading strategy specifically comprises the following steps:
if the angle deviation threshold is less than or equal to 30 degrees, adopting a linear ball jacking algorithm as a ball jacking strategy;
if the angle deviation threshold value is larger than 30 degrees and smaller than or equal to 180 degrees, an arc cut-in dome ball algorithm is adopted as the top ball strategy.
The straight line heading algorithm comprises the following steps:
when the angle value of the first included angle is less than 5 degrees, the direction of the simulated fish is unchanged, and the speed is 14 grades of the preset grade;
when the angle value of the first included angle is more than or equal to 5 degrees and less than 20 degrees, the direction of the simulated fish turns at a preset gear 6-5 or 8-9, and the speed is at a preset level of 12;
when the angle value of the first included angle is more than or equal to 20 degrees and less than 40 degrees, the direction of the simulated fish turns at a preset gear 4-3 or 10-11, and the speed is at a preset level of 8;
when the angle value of the first included angle is more than or equal to 40 degrees and less than 60 degrees, the direction of the simulated fish turns at a preset gear 2-1 or 12-13, and the speed is at a preset level 2;
when the angle value of the first included angle is more than or equal to 60 degrees and less than or equal to 180 degrees, the direction of the simulated fish turns at a preset gear 0 or 14, and the speed is at a preset grade 0 level;
the preset speed level is 0-14 levels, and 15 levels are added in sequence; the preset gears in the direction are 0-14 gears, 15 gears in total, 0-6 gears in left rotation are sequentially decreased, 7 gears in the direction are unchanged, and 8-14 gears in right rotation are sequentially increased.
The arc plunge dome ball algorithm comprises:
simulating a coordinate point of the tangential point towards which the fish is oriented;
wherein, the acquisition of the coordinate point of the tangent point comprises:
making a tangent line to the incision circle of the sphere through the coordinate point of the fish head and obtaining a tangent point;
making a triangle according to the coordinate point of the fish head F', the circle center C of the incision circle and the tangent point, and solving the tangent length according to the pythagorean theorem
Figure BDA0003194080820000071
A modulo length as a rotation vector;
calculating the included angle of deflection of the arc line according to the radius r of the incised circle, namely CT, and the length of the incised line
Figure BDA0003194080820000072
Coordinate points (C.a, C.b) of the rotation center of the simulated fish are obtained, and a direction vector U.x which is (C.a-F '. x)/d and U.y which is (C.b-F'. y)/d are calculated by combining the coordinate points (F '. x, F', y) of the fish head;
calculating a rotation vector Q.x ═ U.x ═ cos (|) to U.y | _ sin (| _ angle |), Q.y ═ U.x | _ sin (| _ angle |) -U.y | _ cos (|);
calculating a coordinate point Q.x (Q.x d + F '. x) and Q.y (Q.y d + F'. y) of the tangent point according to the rotation vector and the tangent length;
wherein, the incised circle of the sphere is a central section of the sphere appearing on 2D.
And 4, judging whether the ball pushing action is finished or not according to the distance between the sphere center of the sphere and the target position (generally, the distance is less than or equal to 15 pixels, and one pixel is a coordinate point), if so, finishing, and if not, circularly executing the steps 2-4.
According to the method for controlling the heading of the 2D simulated fish, the position of the simulated fish is effectively and reasonably adjusted in each cycle period, intelligent heading against an interference function of a bottom platform URWPGSim2D is achieved, heading speed and heading stability are improved, and heading time is shortened under the condition of deviation.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.

Claims (6)

1. A2D artificial fish heading control method is characterized by comprising the following steps:
calculating an angle deviation threshold value caused by a water wave random interference function according to the ideal ball hitting point, the fish head and the coordinate point of the target position;
selecting a heading strategy according to the angle deviation threshold, and performing heading action according to the heading strategy when the simulated fish moves to the target position;
judging whether the ball pushing action is finished or not according to the distance between the center of the ball and the target position, if so, finishing, and if not, circularly executing the steps;
and the coordinate point of the ideal hitting point is obtained according to the fish head of the simulated fish, the sphere center of the sphere and the initialization coordinate point of the target position.
2. The method as claimed in claim 1, wherein the obtaining of the coordinate point of the ideal hitting point comprises:
making a straight line based on the fish head coordinate point and the spherical center coordinate point;
if the coordinate point of the target position is on one side of the straight line, two symmetrical tangent lines are made to the collision circle of the sphere through the coordinate point of the fish head, two tangent points are obtained, and the tangent point on the other side of the straight line is taken as an ideal hitting point;
if the coordinate point of the target position is on the straight line, two intersection points of the straight line and a collision circle of the sphere are obtained, and the intersection point far away from the target point is taken as an ideal hitting point;
wherein the collision circle of the sphere is a central section of the sphere appearing on 2D.
3. The method as claimed in claim 1, wherein the calculating the threshold of the angle deviation caused by the random interference function of the water wave comprises:
making a triangle according to the ideal hitting point, the fish head and the coordinate point of the target position, and solving a first included angle value with the ideal hitting point as a vertex according to the cosine law;
making a triangle according to the ideal hitting point, the fish head offset and the coordinate point of the target position under the water wave random interference function, and solving a second included angle value taking the ideal hitting point as a vertex according to the cosine law;
and obtaining an angle deviation threshold value according to the absolute value of the difference value of the angle value of the first included angle and the angle value of the second included angle.
4. The method as claimed in claim 3, wherein the selecting the heading strategy according to the angle deviation threshold comprises:
if the angle deviation threshold is less than or equal to 30 degrees, adopting a linear ball jacking algorithm as a ball jacking strategy;
if the angle deviation threshold value is larger than 30 degrees and smaller than or equal to 180 degrees, an arc cut-in dome ball algorithm is adopted as the top ball strategy.
5. The method as claimed in claim 4, wherein the linear heading algorithm comprises:
when the angle value of the first included angle is less than 5 degrees, the direction of the simulated fish is unchanged, and the speed is 14 grades of the preset grade;
when the angle value of the first included angle is more than or equal to 5 degrees and less than 20 degrees, the direction of the simulated fish turns at a preset gear 6-5 or 8-9, and the speed is at a preset level of 12;
when the angle value of the first included angle is more than or equal to 20 degrees and less than 40 degrees, the direction of the simulated fish turns at a preset gear 4-3 or 10-11, and the speed is at a preset level of 8;
when the angle value of the first included angle is more than or equal to 40 degrees and less than 60 degrees, the direction of the simulated fish turns at a preset gear 2-1 or 12-13, and the speed is at a preset level 2;
when the angle value of the first included angle is more than or equal to 60 degrees and less than or equal to 180 degrees, the direction of the simulated fish turns at a preset gear 0 or 14, and the speed is at a preset grade 0 level;
the preset speed level is 0-14 levels, and the preset speed levels are 15 levels and are sequentially increased; the preset gears in the direction are 0-14 gears and 15 gears in total, 0-6 gears are left-turning gears and are sequentially decreased in number, 7 gears are unchanged in direction, and 8-14 gears are right-turning gears and are sequentially increased in number.
6. The method as claimed in claim 4, wherein the arc cut dome algorithm comprises:
simulating a coordinate point of the tangential point towards which the fish is oriented;
wherein the obtaining of the coordinate point of the tangent point comprises:
making a tangent line to the incision circle of the sphere through the coordinate point of the fish head and obtaining a tangent point;
making a triangle according to the coordinate point of the fish head, the circle center of the incision circle and the tangent point, and solving the tangent length according to the pythagorean theorem;
calculating to obtain an arc deflection included angle according to the radius of the cut-in circle and the tangent length;
acquiring a coordinate point of the rotation center of the simulated fish, and calculating a direction vector by combining the coordinate point of the fish head;
calculating to obtain a rotation vector according to the direction vector and the arc deflection included angle;
calculating to obtain a coordinate point of the tangent point according to the rotation vector and the tangent length;
wherein the incised circle of the sphere is a central section of the sphere appearing in 2D.
CN202110885848.7A 2021-08-03 2021-08-03 Ball jacking control method for 2D simulated fish Active CN113703464B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110885848.7A CN113703464B (en) 2021-08-03 2021-08-03 Ball jacking control method for 2D simulated fish

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110885848.7A CN113703464B (en) 2021-08-03 2021-08-03 Ball jacking control method for 2D simulated fish

Publications (2)

Publication Number Publication Date
CN113703464A true CN113703464A (en) 2021-11-26
CN113703464B CN113703464B (en) 2023-11-10

Family

ID=78651415

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110885848.7A Active CN113703464B (en) 2021-08-03 2021-08-03 Ball jacking control method for 2D simulated fish

Country Status (1)

Country Link
CN (1) CN113703464B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106485024A (en) * 2016-10-25 2017-03-08 广东工业大学 A kind of emulation machine fish control method based on extreme learning machine and device
CN107065563A (en) * 2016-11-30 2017-08-18 北京理工大学 A kind of monomer machine fish Behavior- Based control strategy based on three-dimensional fuzzy PTP algorithms
CN109291060A (en) * 2018-09-28 2019-02-01 北京信息科技大学 Machine fish heading control method based on region division and interim point auxiliary
KR20190084177A (en) * 2017-12-22 2019-07-16 (주)아이로 (Aqua game system using robot fish
CN111572742A (en) * 2020-05-25 2020-08-25 中国科学院自动化研究所 Underwater robot fish-imitating propelling mechanism
CN113009923A (en) * 2021-03-05 2021-06-22 武汉工程大学 Composite obstacle avoidance method and device for double-joint robot fish

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106485024A (en) * 2016-10-25 2017-03-08 广东工业大学 A kind of emulation machine fish control method based on extreme learning machine and device
CN107065563A (en) * 2016-11-30 2017-08-18 北京理工大学 A kind of monomer machine fish Behavior- Based control strategy based on three-dimensional fuzzy PTP algorithms
KR20190084177A (en) * 2017-12-22 2019-07-16 (주)아이로 (Aqua game system using robot fish
CN109291060A (en) * 2018-09-28 2019-02-01 北京信息科技大学 Machine fish heading control method based on region division and interim point auxiliary
CN111572742A (en) * 2020-05-25 2020-08-25 中国科学院自动化研究所 Underwater robot fish-imitating propelling mechanism
CN113009923A (en) * 2021-03-05 2021-06-22 武汉工程大学 Composite obstacle avoidance method and device for double-joint robot fish

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陶金: "基于动作决策的机器鱼顶球算法", 《兵工自动化》, vol. 29, no. 11, pages 70 - 73 *

Also Published As

Publication number Publication date
CN113703464B (en) 2023-11-10

Similar Documents

Publication Publication Date Title
CN109582027B (en) Improved particle swarm optimization algorithm-based USV cluster collision avoidance planning method
CN106934833B (en) One kind stacking material pick device at random and method
CN101887586B (en) Self-adaptive angular-point detection method based on image contour sharpness
TWI402506B (en) Method and system for motion tracking
CN102915044B (en) A kind of robot head eye coordinate motion control method based on bionic principle
CN102385748B (en) Image registration method
Röfer Evolutionary gait-optimization using a fitness function based on proprioception
CN102930245B (en) Method and system for tracking vehicles
WO2005104654A3 (en) Calculating three-dimensional (3d) voronoi diagrams
CN109190704A (en) The method and robot of detection of obstacles
CN107633518A (en) A kind of product design detection method based on Kinect
CN106989758A (en) Neighbouring underwater robot based on artificial lateral-line system closely cognitive method
CN102012501B (en) Three-dimensional visualization correcting method for radar detecting range influenced by terrain
CN108562872A (en) A kind of method of exceptional value detection when ultra-short baseline under-water acoustic locating calibration
CN113703464A (en) 2D simulated fish heading control method
CN108801262B (en) Method for planning and correcting route of automatic navigation controller of ship
CN107545596B (en) Method for extracting optimal cutting plane of point cloud model
CN110516857A (en) A kind of dead fish fishing paths planning method, server and dead fish finishing device
CN106204724A (en) A kind of 3 D laser scanning measures the method for single wood apparatus derivatorius
Wang et al. Study on edge detection of LIDAR point cloud
CN103559723B (en) A kind of human body tracing method based on self-adaptive kernel function and mean shift
CN103871068A (en) High-precision calibration method based on genetic algorithm
CN105005993A (en) Heterogeneous projection based fast and exact matching method for three-dimensional terrain
Zhang et al. Residual Reinforcement Learning for Motion Control of a Bionic Exploration Robot-RoboDact
CN116975504B (en) Rapid calculation method for satellite reconnaissance coverage area target

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