CN112305525A - Underwater acoustic ranging signal processing method based on spread spectrum technology - Google Patents
Underwater acoustic ranging signal processing method based on spread spectrum technology Download PDFInfo
- Publication number
- CN112305525A CN112305525A CN202011095675.0A CN202011095675A CN112305525A CN 112305525 A CN112305525 A CN 112305525A CN 202011095675 A CN202011095675 A CN 202011095675A CN 112305525 A CN112305525 A CN 112305525A
- Authority
- CN
- China
- Prior art keywords
- signal
- underwater acoustic
- underwater
- data
- spread spectrum
- 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
- 238000005516 engineering process Methods 0.000 title claims abstract description 12
- 238000001228 spectrum Methods 0.000 title claims abstract description 12
- 238000003672 processing method Methods 0.000 title claims abstract description 10
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000005070 sampling Methods 0.000 claims abstract description 7
- 230000008054 signal transmission Effects 0.000 claims abstract description 6
- 230000003321 amplification Effects 0.000 claims abstract description 4
- 230000003750 conditioning effect Effects 0.000 claims abstract description 4
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 4
- 230000021615 conjugation Effects 0.000 claims abstract 2
- 238000000605 extraction Methods 0.000 claims abstract 2
- 238000004364 calculation method Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000004308 accommodation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S11/00—Systems for determining distance or velocity not using reflection or reradiation
- G01S11/14—Systems for determining distance or velocity not using reflection or reradiation using ultrasonic, sonic, or infrasonic waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域technical field
本发明属于水下定位导航领域,尤其涉及一种用于水声测距的信号处理方法。The invention belongs to the field of underwater positioning and navigation, and in particular relates to a signal processing method for underwater acoustic ranging.
背景技术Background technique
水下机器人可以在深水中代替人类完成对海洋资源的探测和开发,为生产节省了大量的成本。同时,在保护国家安全和人民财产方面体现重要价值。水下定位导航技术是水下机器人技术的难点和重点,在实际应用中,由于水下航行器的体积小和重量轻、易受海洋水底环境的影响,导致水下定位导航成为较为棘手的问题。目前,水下机器人中应用的导航系统一般根据其工作方式分为两类:一类是不需要借助外部信息的自主导航系统,如惯性导航系统;另一类是借助于外部信息的导航系统,如GPS/北斗/伽利略等卫星无线电导航系统、水声定位导航系统等。惯性导航系统需要借助于惯性器件的输出完成导航。但由于惯性器件的自身漂移问题会产生误差,其在不修正的情况下定位误差随时间积累,因此在水下长航程任务中无法单独使用。GPS等卫星无线电导航系统定位精度虽然比较高,但电磁波在水中传播时容易被吸收,也无法在长航程任务中使用。水声定位导航系统事先在工作的海域布置参考信标,通过测量水声信号从发射机到接收机的传播时长来测距,可以弥补惯性导航系统和卫星无线电系统在水下机器人长航程任务的不足。但传统的水声定位导航系统缺点也非常明显,由于其问答工作模式导致其容纳水下机器人的数目有限,并且定位服务更新速率太慢。因此,可以考虑借助于扩频技术,研究信号的处理算法,解决以上问题。Underwater robots can replace humans in deep water to complete the detection and development of marine resources, saving a lot of costs for production. At the same time, it reflects important value in protecting national security and people's property. Underwater positioning and navigation technology is the difficulty and focus of underwater robot technology. In practical applications, due to the small size and light weight of underwater vehicles, which are easily affected by the underwater environment, underwater positioning and navigation has become a more difficult problem. . At present, the navigation systems used in underwater robots are generally divided into two categories according to their working methods: one is an autonomous navigation system that does not require external information, such as an inertial navigation system; the other is a navigation system that relies on external information. Such as GPS/Beidou/Galileo and other satellite radio navigation systems, underwater acoustic positioning and navigation systems, etc. The inertial navigation system needs to complete the navigation with the help of the output of the inertial device. However, due to the drift problem of the inertial device itself, errors will occur, and the positioning error will accumulate over time without correction, so it cannot be used alone in long-range underwater missions. Although the positioning accuracy of satellite radio navigation systems such as GPS is relatively high, electromagnetic waves are easily absorbed when propagating in water and cannot be used in long-range missions. The underwater acoustic positioning and navigation system arranges reference beacons in the working sea area in advance, and measures the distance by measuring the propagation time of the underwater acoustic signal from the transmitter to the receiver, which can make up for the inertial navigation system and the satellite radio system in the long-range mission of the underwater robot. insufficient. However, the shortcomings of the traditional underwater acoustic positioning and navigation system are also very obvious. Due to its question-and-answer mode, the number of underwater robots it can accommodate is limited, and the update rate of the positioning service is too slow. Therefore, it can be considered to study the signal processing algorithm with the help of spread spectrum technology to solve the above problems.
本发明基于扩频技术,不同的信标发射机向外发射定位服务信号,水下机器人接收机接收水声信号进行自我位置解算,避免了传统水声定位系统问答式作业方式造成的系统容纳水下机器人数目有限的问题。通过同一码族不同扩频码来区分不同的发射机和发射信号时刻,从实现更新速度较快的测距服务。为了保证以上功能的实现,需要信号发射机和信号接收机的时钟同步,信号发射时刻为时钟的整秒。本发明为水下机器人协同作业提供了理论支撑和工程应用保障。The invention is based on the spread spectrum technology, different beacon transmitters transmit positioning service signals to the outside, and the underwater robot receiver receives the underwater acoustic signals for self-position calculation, avoiding the system accommodation caused by the question-and-answer operation mode of the traditional underwater acoustic positioning system. The problem of limited number of underwater robots. Different transmitters and transmission signal timings are distinguished by different spreading codes of the same code family, so as to realize the ranging service with faster update speed. In order to ensure the realization of the above functions, the clocks of the signal transmitter and the signal receiver need to be synchronized, and the signal transmission time is an entire second of the clock. The invention provides theoretical support and engineering application guarantee for the collaborative operation of underwater robots.
发明内容SUMMARY OF THE INVENTION
本发明涉及一种基于扩频技术的水声测距信号处理方法,包括以下步骤:The present invention relates to an underwater acoustic ranging signal processing method based on spread spectrum technology, comprising the following steps:
A1:对水听器接收到的数据进行放大滤波等信号调理后,进而通过模拟数字采集电路进行数字采集形成原始数据C1;A1: After the data received by the hydrophone is subjected to signal conditioning such as amplification and filtering, the original data C1 is formed by digital acquisition through the analog-digital acquisition circuit;
A2:对原始数据C1进行正交基带解调、低通滤波、降采样等信号处理过程形成初步处理数据C2;A2: Perform signal processing processes such as quadrature baseband demodulation, low-pass filtering, and downsampling on the original data C1 to form preliminary processed data C2;
A3:根据已有信号参数与生成原始参考数据C3,进行正交基带解调、低通滤波、降采样和共轭处理形成初步参考数据C4;A3: According to the existing signal parameters and the generated original reference data C3, perform quadrature baseband demodulation, low-pass filtering, down-sampling and conjugate processing to form preliminary reference data C4;
A4:将初步处理数据C2和初步参考数据C4进行卷积运算后取模,得到检测信号C5;A4: Perform a convolution operation on the preliminary processed data C2 and the preliminary reference data C4 and take the modulo to obtain a detection signal C5;
A5:将检测信号C5与阈值进行比较,大于阈值后判断信号传播时间,进而求出发射机与接收机距离。A5: Compare the detection signal C5 with the threshold value, and judge the signal propagation time when it is greater than the threshold value, and then obtain the distance between the transmitter and the receiver.
为了准确估计出信号传播时长,需要保证信标信号发射机和水下机器人信号接收机的时钟同步。In order to accurately estimate the signal propagation time, it is necessary to ensure that the clocks of the beacon signal transmitter and the underwater robot signal receiver are synchronized.
信号发射时刻为时钟的整秒开始,即整数秒零毫秒零微秒零纳秒。The signal transmission time is the start of the whole second of the clock, that is, the whole number of seconds, zero milliseconds, zero microseconds, and zero nanoseconds.
通过同一码族不同扩频码来区分不同发射机和发射信号时刻。如果发射机数目为N,标志发射信号整秒时刻数目记做M,同一码族共有D组不同的伪随机码,则三者关系需要满足:Different transmitters and transmission signal moments are distinguished by different spreading codes of the same code family. If the number of transmitters is N, the number of signals transmitted in a whole second is marked as M, and there are D groups of different pseudorandom codes in the same code family, the relationship between the three needs to be satisfied:
N×M≤DN×M≤D
为了减小由于遍历不同扩频码产生的计算量,需要尽量减少发射信号整秒时刻数目M。如果水声测距最大作用距离为L,水中声速为c,水声定位服务更新频率为H秒,则M需要满足:In order to reduce the amount of calculation caused by traversing different spreading codes, it is necessary to minimize the number of times M in a whole second of the transmitted signal. If the maximum working distance of underwater acoustic ranging is L, the speed of sound in water is c, and the update frequency of underwater acoustic positioning service is H seconds, then M needs to satisfy:
其中,为向上取整。in, to round up.
本发明的有益效果与现有技术相比,本发明的有益效果是:Compared with the prior art, the beneficial effects of the present invention are:
针对多个水下机器人协同作业需求定位导航服务时,可以采用信标发射机发射信号,水下机器人只接收声学信号的服务模式。由此避免了传统水声定位系统问答模式导致的系统容纳水下机器人的数目有限的情况。此外,本发明采用不同的伪随机码信号标致信号的发射起始时刻,使得水下机器人较高的定位服务更新速率需求成为可能。本发明为多水下机器人协同工作提供了理论基础和工程应用方法。When positioning and navigation services are required for the collaborative operation of multiple underwater robots, beacon transmitters can be used to transmit signals, and the underwater robots only receive acoustic signals. In this way, the situation that the system accommodates the limited number of underwater robots caused by the question-and-answer mode of the traditional underwater acoustic positioning system is avoided. In addition, the present invention adopts different pseudo-random code signals to mark the start time of transmitting the signal, which makes it possible for the underwater robot to have a higher update rate requirement of the positioning service. The invention provides a theoretical basis and an engineering application method for the cooperative work of multiple underwater robots.
附图说明Description of drawings
图1是信号发生器产生的信号仿真图。Figure 1 is a simulation diagram of a signal generated by a signal generator.
图2是模拟水听器接收到的测距信号。Figure 2 is the ranging signal received by the simulated hydrophone.
图3是通过低通滤波后信号仿真图。Figure 3 is a signal simulation diagram after low-pass filtering.
图4是经过处理最后得到检测信号示意图。FIG. 4 is a schematic diagram of a detection signal finally obtained after processing.
具体实施方式Detailed ways
本发明涉及一种基于扩频技术的水声测距信号处理方法,包括以下步骤:The present invention relates to an underwater acoustic ranging signal processing method based on spread spectrum technology, comprising the following steps:
A1:对水听器接收到的数据进行放大滤波等信号调理后,进而通过模拟数字采集电路进行数字采集形成原始数据C1;A1: After the data received by the hydrophone is subjected to signal conditioning such as amplification and filtering, the original data C1 is formed by digital acquisition through the analog-digital acquisition circuit;
A2:对原始数据C1进行正交基带解调、低通滤波、降采样等信号处理过程形成初步处理数据C2;A2: Perform signal processing processes such as quadrature baseband demodulation, low-pass filtering, and downsampling on the original data C1 to form preliminary processed data C2;
A3:根据已有信号参数与生成原始参考数据C3,进行正交基带解调、低通滤波、降采样和共轭处理形成初步参考数据C4;A3: According to the existing signal parameters and the generated original reference data C3, perform quadrature baseband demodulation, low-pass filtering, down-sampling and conjugate processing to form preliminary reference data C4;
A4:将初步处理数据C2和初步参考数据C4进行卷积运算后取模,得到检测信号C5;A4: Perform a convolution operation on the preliminary processed data C2 and the preliminary reference data C4 and take the modulo to obtain a detection signal C5;
A5:将检测信号C5与阈值进行比较,大于阈值后判断信号传播时间,进而求出发射机与接收机距离。A5: Compare the detection signal C5 with the threshold value, and judge the signal propagation time when it is greater than the threshold value, and then obtain the distance between the transmitter and the receiver.
进一步的,为了准确估计出信号传播时长,需要保证信标信号发射机和水下机器人信号接收机的时钟同步。Further, in order to accurately estimate the signal propagation time, it is necessary to ensure that the clocks of the beacon signal transmitter and the underwater robot signal receiver are synchronized.
进一步的,信号发射时刻为时钟的整秒开始,即整数秒零毫秒零微秒零纳秒。Further, the signal transmission time is the start of an entire second of the clock, that is, an integer second, zero milliseconds, zero microseconds, and zero nanoseconds.
进一步的,通过同一码族不同扩频码来区分不同发射机和发射信号时刻。如果发射机数目为N,标志发射信号整秒时刻数目记做M,同一码族共有D组不同的伪随机码,则三者关系需要满足:Further, different transmitters and transmission signal times are distinguished by different spreading codes of the same code family. If the number of transmitters is N, the number of signals transmitted in a whole second is marked as M, and there are D groups of different pseudorandom codes in the same code family, the relationship between the three needs to be satisfied:
N×M≤DN×M≤D
进一步的,为了减小由于遍历不同扩频码产生的计算量,需要尽量减少发射信号整秒时刻数目M。如果水声测距最大作用距离为L,水中声速为c,水声定位服务更新频率为H秒,则M需要满足:Further, in order to reduce the amount of calculation caused by traversing different spreading codes, it is necessary to reduce the number M of the whole second of the transmitted signal as much as possible. If the maximum working distance of underwater acoustic ranging is L, the speed of sound in water is c, and the update frequency of underwater acoustic positioning service is H seconds, then M needs to satisfy:
其中,为向上取整。in, to round up.
结合仿真数据来详述上述解算步骤。The above calculation steps are described in detail in combination with the simulation data.
假设系统的工作深度要求为6000m;Assume that the working depth requirement of the system is 6000m;
信号中心频率定为11kHz;The signal center frequency is set to 11kHz;
传播损失为86dB;The propagation loss is 86dB;
噪声源级为85dB;The noise source level is 85dB;
换能器的声源级为171dB;The sound source level of the transducer is 171dB;
检测阈值为20dB;The detection threshold is 20dB;
扩频码采用7级Gold码,码速率为2.2kHz;The spread spectrum code adopts 7-level Gold code, and the code rate is 2.2kHz;
信号采样率为70kHz;The signal sampling rate is 70kHz;
解调频率为8kHz;The demodulation frequency is 8kHz;
低通滤波的通带为0~6kHz。The passband of the low-pass filter is 0 to 6 kHz.
仿真得到的结果如图1~图4所示。从图4可以看出,按照图1产生的发射信号在信噪比为0dB的情况下,本发明提供的信号处理方法可以实现水下声学测距信号的正确处理。The simulation results are shown in Figures 1 to 4. As can be seen from FIG. 4 , when the signal-to-noise ratio of the transmitted signal generated according to FIG. 1 is 0 dB, the signal processing method provided by the present invention can realize the correct processing of the underwater acoustic ranging signal.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011095675.0A CN112305525B (en) | 2020-10-14 | 2020-10-14 | A method of underwater acoustic ranging signal processing based on spread spectrum technology |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011095675.0A CN112305525B (en) | 2020-10-14 | 2020-10-14 | A method of underwater acoustic ranging signal processing based on spread spectrum technology |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112305525A true CN112305525A (en) | 2021-02-02 |
CN112305525B CN112305525B (en) | 2023-09-22 |
Family
ID=74488035
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011095675.0A Active CN112305525B (en) | 2020-10-14 | 2020-10-14 | A method of underwater acoustic ranging signal processing based on spread spectrum technology |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112305525B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08114665A (en) * | 1994-10-14 | 1996-05-07 | Japan Radio Co Ltd | Sonobui position detection method |
CN103988094A (en) * | 2011-08-16 | 2014-08-13 | 欧洲太空署 | Navigation system using spreading codes based on pseudorandom noise sequences |
CN109407098A (en) * | 2018-08-30 | 2019-03-01 | 中国船舶重工集团公司第七〇五研究所 | A kind of implementation method for the parallel processing algorithm handling orthogonal pulses string signal |
CN111427030A (en) * | 2020-06-12 | 2020-07-17 | 飞马滨(北京)科技研发有限公司 | Underwater acoustic signal modulation and processing method for positioning underwater robot |
CN111580038A (en) * | 2020-04-02 | 2020-08-25 | 浙江工业大学 | An FPGA-based acoustic underwater beacon signal processing system |
-
2020
- 2020-10-14 CN CN202011095675.0A patent/CN112305525B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08114665A (en) * | 1994-10-14 | 1996-05-07 | Japan Radio Co Ltd | Sonobui position detection method |
CN103988094A (en) * | 2011-08-16 | 2014-08-13 | 欧洲太空署 | Navigation system using spreading codes based on pseudorandom noise sequences |
CN109407098A (en) * | 2018-08-30 | 2019-03-01 | 中国船舶重工集团公司第七〇五研究所 | A kind of implementation method for the parallel processing algorithm handling orthogonal pulses string signal |
CN111580038A (en) * | 2020-04-02 | 2020-08-25 | 浙江工业大学 | An FPGA-based acoustic underwater beacon signal processing system |
CN111427030A (en) * | 2020-06-12 | 2020-07-17 | 飞马滨(北京)科技研发有限公司 | Underwater acoustic signal modulation and processing method for positioning underwater robot |
Also Published As
Publication number | Publication date |
---|---|
CN112305525B (en) | 2023-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20190250244A1 (en) | Underwater source node positioning method | |
CN103823205B (en) | Underwater locating navigation system and method | |
CN103823229B (en) | A DGPS buoy-based underwater positioning and navigation system and method | |
CN110366095B (en) | An underwater node distributed positioning system and method | |
CN105137394B (en) | Ultra-short baseline acoustic positioning system and its location algorithm based on Maximum-likelihood estimation | |
CN104272132B (en) | Determining position of underwater node | |
CN101377542B (en) | Method and device for estimating GNSS reflected signal parameter | |
Han et al. | Error control and adjustment method for underwater wireless sensor network localization | |
CN105974363B (en) | A kind of underwater enhancing localization method based on sounding error correlation | |
CN109814069B (en) | Underwater mobile node passive positioning method and system based on single positioning beacon | |
CN103323815B (en) | A kind of under-water acoustic locating method based on the equivalent velocity of sound | |
CN107505597A (en) | A kind of ultra short baseline locating system and method based on depth information | |
CN111366897B (en) | Signal design and processing method for high-precision underwater sound positioning system | |
CN104180804A (en) | Single reference node underwater vehicle integrated navigation method based on underwater information network | |
CN111580080A (en) | A method and system for effective wave height inversion based on shipborne GNSS-R | |
CN102183770A (en) | GPS (Global Positioning System) pseudo-random code tracking loop for preventing multi-path interference and method for preventing multi-path interference thereof | |
CN104197939A (en) | Multi-reference-point under-water vehicle combination navigation method based on underwater information network | |
US20150061931A1 (en) | Methods and systems for pseudo-random coding in a wide area positioning system | |
CN106842209A (en) | A kind of multirobot cooperates with underwater terrain matching navigation system and method | |
CN102183741A (en) | Method for realizing high frame rate unambiguous positioning in long baseline asynchronous acoustic positioning system | |
Li et al. | Long baseline acoustic localization based on track-before-detect in complex underwater environments | |
CN108882167B (en) | A Navigation and Positioning Method Using Signal Strength in Underwater Acoustic Communication Network | |
CN112305525B (en) | A method of underwater acoustic ranging signal processing based on spread spectrum technology | |
CN111427030B (en) | Underwater acoustic signal modulation and processing method for positioning underwater robot | |
Kottege et al. | Underwater acoustic localization for small submersibles |
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 |