CN113489552B - 一种基于时频谱矩阵局部方差的跳频信号检测方法 - Google Patents

一种基于时频谱矩阵局部方差的跳频信号检测方法 Download PDF

Info

Publication number
CN113489552B
CN113489552B CN202110556552.0A CN202110556552A CN113489552B CN 113489552 B CN113489552 B CN 113489552B CN 202110556552 A CN202110556552 A CN 202110556552A CN 113489552 B CN113489552 B CN 113489552B
Authority
CN
China
Prior art keywords
frequency
time
spectrum
signal
frequency hopping
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
CN202110556552.0A
Other languages
English (en)
Other versions
CN113489552A (zh
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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202110556552.0A priority Critical patent/CN113489552B/zh
Publication of CN113489552A publication Critical patent/CN113489552A/zh
Application granted granted Critical
Publication of CN113489552B publication Critical patent/CN113489552B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B17/00—Monitoring; Testing
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B1/00—Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69—Spread spectrum techniques
    • H04B1/713—Spread spectrum techniques using frequency hopping
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Noise Elimination (AREA)

Abstract

本发明提供了一种基于时频谱矩阵局部方差的跳频信号检测方法，该方法包括以下步骤：S1.根据跳频信号模型和定频信号模型确定接收机截获的观测信号是否存在跳频信号或者定频信号，若存在，则执行后续步骤；S2.获取观测信号的短时傅里叶变换STFTx(m,n)，并计算时频谱Px(m,n)；S3.将时频谱Px(m,n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure DEST_PATH_FDA0003832071300000011
S4.在原时频谱Px(m,n)上减去对应频率的平均功率谱
Figure DEST_PATH_FDA0003832071300000012
得到时频对消谱Psub1(m,n)；S5.对所得到的时频对消谱Psub1(m,n)进行归一化处理，并计算每个频点时频谱的方差d(m)和局部方差T。本方案将短时傅里叶变换、时频对消方法与方差特点结合起来对跳频信号进行检测，具有有效性和实用性等优点。

Description

一种基于时频谱矩阵局部方差的跳频信号检测方法
技术领域
本发明属于无线通信中认知无线电领域，尤其是涉及一种基于时频谱矩阵局部方差的跳频信号检测方法。
背景技术
跳频是最常用的扩频方式之一，其工作原理是指收发双方传输信号的载波频率按照预定规律进行离散变化的通信方式，也就是说，通信中使用的载波频率受伪随机变化码的控制而随机跳变。从通信技术的实现方式来说，"跳频"是一种用码序列进行多频频移键控的通信方式，也是一种码控载频跳变的通信系统。
跳频通信具有截获率低、抗干扰能力强和组网方便等特点，目前已成为军事通信的主导技术，所以跳频信号的检测是通信侦察的主要任务之一。现有的跳频信号检测方法大多是从跳频信号的频域出发，而从跳频信号的频域出发不能充分利用跳频信号的时频谱特点，检测性能也有待提高。也有一些关于功率谱对消方法的研究，但是检测概率不能在低信噪比下达到很好的效果，很难满足实际需求。
发明内容
本发明的目的是针对上述问题，提供一种基于时频谱矩阵局部方差的跳频信号检测方法。
为达到上述目的，本发明采用了下列技术方案：
一种基于时频谱矩阵局部方差的跳频信号检测方法，该方法包括以下步骤：
S1.根据跳频信号模型和定频信号模型确定接收机截获的观测信号是否存在跳频信号或者定频信号，若存在，则执行后续步骤；
S2.获取观测信号的短时傅里叶变换STFTx(m,n)，并计算时频谱Px(m,n)；
S3.将时频谱Px(m,n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure GDA0003824632690000021
S4.在原时频谱Px(m,n)上减去对应频率的平均功率谱
Figure GDA0003824632690000022
得到时频对消谱Psub1(m，n)；
S5.对所得到的时频对消谱Psub1(m，n)进行归一化处理，并计算每个频点时频谱的方差d(m)和局部方差T。
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S1中，所述的跳频信号模型为公式(1)
Figure GDA0003824632690000023
其中，T为总观测时间，0≤t≤T，K为观测时间内的总跳数，TH为跳周期，a(t)为基带复包络，As为跳频信号的幅度，fk和
Figure GDA0003824632690000024
分别表示第k跳的载波频率和相位，
Figure GDA0003824632690000025
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S1中，所述定频信号模型为公式(2)
Figure GDA0003824632690000026
其中，AJ、fJ和
Figure GDA0003824632690000027
分别表示为定频信号的幅度、载波频率和相位。
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S1中，通过以下方式确定所述的观测信号：
对于给定的连续信号x(t)，其短时傅里叶变换表示为公式(3)：
Figure GDA0003824632690000028
其中h(t)是窗函数，h*(t)是它的共轭形式；
令仅存在噪声的假设为H0，噪声中存在跳频信号和定频信号的假设为H1，得到接收机截获的观测信号为公式(4)：
Figure GDA0003824632690000031
则观测信号的短时傅里叶变换表示为公式(5)：
Figure GDA0003824632690000032
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，在步骤S2中，通过以下方式获取所述的短时傅里叶变换STFTx(m,n)：
对观测信号x(t)进行采样，窗长为P，每次滑动一定的长度，计算得到其短时傅里叶变换STFTx(m，n)；
通过以下方式计算所述的时频谱Px(m,n)：
对观测信号的短时傅里叶变换进行取模值平方计算，即令时频谱为Px(m，n)＝|STFTx(m，n)|2，其中，频率维度m＝1，2，…，M，M＝P/2，时间维度n＝1，2，…，N。
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S3中，通过以下公式(6)获得只与频率相关的平均功率谱
Figure GDA0003824632690000033
Figure GDA0003824632690000034
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S4中，通过以下公式(7)得到时频对消谱Psub1(m，n)：
Figure GDA0003824632690000035
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S5中，通过以下公式(8)对所得到的时频对消谱进行归一化处理，以避免噪声方差对信号检测结果产生影响：
Figure GDA0003824632690000036
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S5中，通过以下公式(9)计算每个频点时频谱的方差d(m)：
Figure GDA0003824632690000041
在上述的基于时频谱矩阵局部方差的跳频信号检测方法中，步骤S5中，通过以下公式(10)计算时频谱的局部方差T：
Figure GDA0003824632690000042
本发明的优点在于：
1、本文提出了一种基于时频谱矩阵局部方差的跳频信号检测方法，将短时傅里叶变换、时频对消方法与方差特点结合起来，具有有效性和实用性等优点。
2、提出改进的时频对消算法，减少了定频信号对跳频信号检测结果的影响。
3、充分利用噪声、定频信号和跳频信号的对消时频谱局部方差的差异性，使得跳频信号在低信噪比下仍能获得较高的识别率。
附图说明
图1为噪声对消时频谱局部方差的差异性示图；
图2为噪声+跳频信号对消时频谱局部方差的差异性示图；
图3为噪声+定频信号对消时频谱局部方差的差异性示图；
图4为噪声时频谱；
图5为跳频信号时频谱。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细的说明。
如图1和图3所示，本方案提出一种基于时频谱矩阵局部方差的跳频信号检测方法，该方法包括以下步骤：
S1.根据跳频信号模型和定频信号模型确定接收机截获的观测信号是否存在跳频信号或者定频信号，若存在，则执行后续步骤；
其中跳频信号模型为公式(1)：
Figure GDA0003824632690000051
其中，T为总观测时间，0≤t≤T，K为观测时间内的总跳数，TH为跳周期，a(t)为基带复包络，AS为跳频信号的幅度，fk和
Figure GDA0003824632690000057
分别表示第k跳的载波频率和相位，
Figure GDA0003824632690000052
定频信号模型为公式(2)：
Figure GDA0003824632690000053
其中，AJ、fJ和
Figure GDA0003824632690000054
分别表示为定频信号的幅度、载波频率和相位。
对于给定的连续信号x(t)，其短时傅里叶变换的一般表达式为公式(3)：
Figure GDA0003824632690000055
其中h(t)是窗函数，h*(t)是它的共轭形式；
令仅存在噪声的假设为H0，噪声中存在跳频信号和定频信号的假设为H1，得到接收机截获的观测信号为公式(4)：
Figure GDA0003824632690000056
则观测信号的短时傅里叶变换表示为公式(5)，将观测信号分为了噪声和信号，将噪声、跳频和定频信号分别用各自的短时傅里叶代替：
Figure GDA0003824632690000061
S2.获取观测信号的短时傅里叶变换STFTx(m,n)，并计算时频谱Px(m,n)；
具体通过以下方式获取短时傅里叶变换STFTx(m,n)：
对观测信号x(t)进行采样，窗长为P，每次滑动一定的长度，计算得到其短时傅里叶变换STFTx(m，n)，
通过以下方式计算时频谱Px(m,n)：
对观测信号的短时傅里叶变换进行取模值平方计算，即令时频谱为Px(m，n)＝|STFTx(m，n)|2，其中，频率维度m＝1，2，…，M，M＝P/2，时间维度n＝1，2，…，N。
通过该步骤S2获取观测信号的短时傅里叶变换，最终目的是得到信号和噪声的时频谱，通过比较时频谱的差异，进行后续操作，如图4和图5所示，噪声的时频谱在整个M×N维平面上均匀分布；跳频信号在不同时间段、不同频点上时频谱不为0；因此可得跳频信号时频谱的局部方差最大，本方案正是基于该原理进行的跳频信号检测。当有定频干扰信号存在时，它在某个频率上存在时频谱不为0，其他为0。所以首先需要降低定频信号对跳频信号检测的影响，本方案具体采用步骤S3，S4消除定频信号时频谱分量。
S3.将时频谱Px(m,n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure GDA0003824632690000062
具体如公式(6)
Figure GDA0003824632690000063
S4.在原时频谱Px(m,n)上减去对应频率的平均功率谱
Figure GDA0003824632690000064
得到时频对消谱Psub1(m，n)具体如公式(7)，
Figure GDA0003824632690000065
由于噪声与跳频信号的差异，又因为定频信号在某个频率上存在时频谱不为0，其他为0，所以通过S4的时频对消方法，消去定频信号的影响，最终达到检测出跳频信号的作用。
S5.对所得到的时频对消谱Psub1(m，n)进行归一化处理，并计算每个频点时频谱的方差d(m)和局部方差T。具体通过公式(8)对时频对消谱进行归一化处理的能够避免噪声方差对信号检测结果产生影响。
Figure GDA0003824632690000071
具体通过以下公式(9)计算每个频点时频谱的方差d(m)：
Figure GDA0003824632690000072
通过公式(10)计算时频谱的局部方差T：
Figure GDA0003824632690000073
通过步骤S5去除噪声方差，目的是为了消除噪声方差对信号检测的影响。最终得到的局部方差T，是为了可以体现出噪声与信号的差异性，利用差异性，最终得到跳频信号的检测规则：
Figure GDA0003824632690000074
其中λ为判决门限，需根据纽曼皮尔逊准则确定，即给定虚警概率，通过蒙特卡洛仿真得到门限。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
尽管本文较多地使用了跳频信号、定频信号、观测信号、短时傅里叶变换、时频谱、平均功率谱、时频对消谱等术语，但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质；把它们解释成任何一种附加的限制都是与本发明精神相违背的。

Claims (10)

1.一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，该方法包括以下步骤：
S1.根据跳频信号模型和定频信号模型确定接收机截获的观测信号是否存在跳频信号或者定频信号，若存在，则执行后续步骤；
S2.获取观测信号的短时傅里叶变换STFTx(m,n)，并计算时频谱Px(m,n)；
S3.将时频谱Px(m,n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure FDA0003832071300000011
S4.在原时频谱Px(m,n)上减去对应频率的平均功率谱
Figure FDA0003832071300000012
得到时频对消谱Psub1(m，n)；
S5.对所得到的时频对消谱Psub1(m，n)进行归一化处理得到Psub(m，n)，并计算每个频点时频谱的方差d(m)和局部方差T，然后根据跳频信号的检测规则识别跳频信号：
Figure FDA0003832071300000013
其中λ为判决门限，T为局部方差，H0表示检测结果为仅存在噪声，H1表示检测结果存在跳频信号。
2.根据权利要求1所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，所述的跳频信号模型为公式(1)
Figure FDA0003832071300000014
其中，T为总观测时间，0≤t≤T，K为观测时间内的总跳数，TH为跳周期，a(t)为基带复包络，As为跳频信号的幅度，fk和
Figure FDA0003832071300000015
分别表示第k跳的载波频率和相位，
Figure FDA0003832071300000021
3.根据权利要求2所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，所述定频信号模型为公式(2)
Figure FDA0003832071300000022
其中，AJ、fJ和
Figure FDA0003832071300000023
分别表示为定频信号的幅度、载波频率和相位。
4.根据权利要求3所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，通过以下方式确定所述的观测信号：
对于给定的连续信号x(t)，其短时傅里叶变换表示为公式(3)：
Figure FDA0003832071300000024
其中h(t)是窗函数，h*(t)是它的共轭形式；
令仅存在噪声的假设为H0，噪声中存在跳频信号和定频信号的假设为H1，得到接收机截获的观测信号为公式(4)：
Figure FDA0003832071300000025
则观测信号的短时傅里叶变换表示为公式(5)：
Figure FDA0003832071300000026
5.根据权利要求4所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，在步骤S2中，通过以下方式获取所述的短时傅里叶变换STFTx(m,n)：
对观测信号x(t)进行采样，窗长为P，每次滑动一定的长度，计算得到其短时傅里叶变换STFTx(m，n)；
通过以下方式计算所述的时频谱Px(m,n)：
对观测信号的短时傅里叶变换进行取模值平方计算，即令时频谱为Px(m，n)＝|STFTx(m，n)|2，其中，频率维度m＝1，2，…，M，M＝P/2，时间维度n＝1，2，…，N。
6.根据权利要求5所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S3中，通过以下公式(6)获得只与频率相关的平均功率谱
Figure FDA0003832071300000031
Figure FDA0003832071300000032
7.根据权利要求6所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S4中，通过以下公式(7)得到时频对消谱Psub1(m，n)：
Figure FDA0003832071300000033
8.根据权利要求7所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S5中，通过以下公式(8)对所得到的时频对消谱进行归一化处理，以避免噪声方差对信号检测结果产生影响：
Figure FDA0003832071300000034
9.根据权利要求8所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S5中，通过以下公式(9)计算每个频点时频谱的方差d(m)：
Figure FDA0003832071300000035
10.根据权利要求9所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S5中，通过以下公式(10)计算时频谱的局部方差T：
Figure FDA0003832071300000036
CN202110556552.0A 2021-05-21 2021-05-21 一种基于时频谱矩阵局部方差的跳频信号检测方法 Active CN113489552B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110556552.0A CN113489552B (zh) 2021-05-21 2021-05-21 一种基于时频谱矩阵局部方差的跳频信号检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110556552.0A CN113489552B (zh) 2021-05-21 2021-05-21 一种基于时频谱矩阵局部方差的跳频信号检测方法

Publications (2)

Publication Number Publication Date
CN113489552A CN113489552A (zh) 2021-10-08
CN113489552B true CN113489552B (zh) 2022-10-18

Family

ID=77933419

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110556552.0A Active CN113489552B (zh) 2021-05-21 2021-05-21 一种基于时频谱矩阵局部方差的跳频信号检测方法

Country Status (1)

Country Link
CN (1) CN113489552B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113820006B (zh) * 2021-11-22 2022-03-11 北京思安利鑫科技有限公司 一种弱信噪比单频正弦信号的参数估计方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001244851A (ja) * 2000-02-25 2001-09-07 Nippon Telegr & Teleph Corp <Ntt> 周波数ホッピング通信方式用送受信装置及び周波数ホッピング通信方法
CN107407703A (zh) * 2016-01-12 2017-11-28 阿澳尼亚股份有限公司 具有提高的实时带宽的频谱分析仪
CN109802736A (zh) * 2019-03-13 2019-05-24 西安电子科技大学 一种频谱感知的方法、装置及计算机存储介质
CN110336587A (zh) * 2019-07-16 2019-10-15 电子科技大学 一种多跳频信号侦察中获取组合时频分布的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001244851A (ja) * 2000-02-25 2001-09-07 Nippon Telegr & Teleph Corp <Ntt> 周波数ホッピング通信方式用送受信装置及び周波数ホッピング通信方法
CN107407703A (zh) * 2016-01-12 2017-11-28 阿澳尼亚股份有限公司 具有提高的实时带宽的频谱分析仪
CN109802736A (zh) * 2019-03-13 2019-05-24 西安电子科技大学 一种频谱感知的方法、装置及计算机存储介质
CN110336587A (zh) * 2019-07-16 2019-10-15 电子科技大学 一种多跳频信号侦察中获取组合时频分布的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"一种改进的跳频信号时频分析方法";冯维婷等;《西安邮电大学学报》;20180531;第23卷(第3期);全文 *

Also Published As

Publication number Publication date
CN113489552A (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
Bhargavi et al. Performance comparison of energy, matched-filter and cyclostationarity-based spectrum sensing
Manesh et al. Performance evaluation of spectrum sensing techniques for cognitive radio systems
CN102904604B (zh) 一种窄带干扰抑制方法和装置
JP2008206152A (ja) マルチキャリア受信信号における狭帯域干渉を検出するための受信装置
WO2008055094A2 (en) Method and apparatus for spur cancellation in an orthogonal frequency division multiplexing communication system
Suseela et al. Non-cooperative spectrum sensing techniques in cognitive radio-a survey
Khalaf et al. Blind spectrum detector for cognitive radio using compressed sensing and symmetry property of the second order cyclic autocorrelation
Kumar et al. Improved wavelet transform based edge detection for wide band spectrum sensing in Cognitive Radio
CN113489552B (zh) 一种基于时频谱矩阵局部方差的跳频信号检测方法
US7899650B2 (en) Apparatus and method for detecting and identifying raman signatures
Ren et al. Wideband spectrum sensing based on coprime sampling
CN113595588B (zh) 一种基于时频谱熵的跳频信号感知方法
Molchanov et al. A novel bispectrum-based frequency encoding technique developed for noisy and multipath fading radio links
Honghan et al. A synchronization acquisition algorithm based on the frequency hopping pulses combining
Yeste-Ojeda et al. Cyclostationarity-based signal separation in interceptors based on a single sensor
Deepa et al. Cyclostationary-based architectures for spectrum sensing in IEEE 802.22 WRAN
Sansoy et al. Cyclostationary feature based detection using window method in SIMO cognitive radio system
CN111262754B (zh) 一种基于畸变功率去除的窄带干扰检测方法及其通信装置
Arjun et al. Performance Analysis of Wavelet based Spectrum Sensing and Conventional Spectrum Sensing in Fading Environment for Cognitive Radios
Bektas et al. Energy based spectrum sensing using wavelet transform for fading channels
Kobayashi et al. Improved MB cognitive radio spectrum sensing using wavelet spectrum filtering
Avila et al. Improved method of sensing the spectrum holes in cognitive radio
Hosseini et al. Enhanced FCME thresholding for wavelet‐based cognitive UWB over fading channels
CN113904703B (zh) 连续频谱随机信号载波通信方法
KR101203718B1 (ko) 통신 장치 및 방법

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