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

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

Info

Publication number
CN113489552A
CN113489552A CN202110556552.0A CN202110556552A CN113489552A CN 113489552 A CN113489552 A CN 113489552A CN 202110556552 A CN202110556552 A CN 202110556552A CN 113489552 A CN113489552 A CN 113489552A
Authority
CN
China
Prior art keywords
frequency
time
spectrum
signal
hopping signal
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
CN202110556552.0A
Other languages
English (en)
Other versions
CN113489552B (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.获取观测信号的短时傅里叶变换
Figure DDA0003077454830000013
并计算时频谱Px(m,n)；S3.将时频谱Px(m,n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure DDA0003077454830000012
S4.在原时频谱Px(m,n)上减去对应频率的平均功率谱
Figure DDA0003077454830000011
得到时频对消谱Psub(m,n)；S5.对所得到的时频对消谱Psub(m,n)进行归一化处理，并计算每个频点时频谱的方差d(m)和局部方差T。本方案将短时傅里叶变换、时频对消方法与方差特点结合起来对跳频信号进行检测，具有有效性和实用性等优点。

Description

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

Claims (10)

1.一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，该方法包括以下步骤：
S1.根据跳频信号模型和定频信号模型确定接收机截获的观测信号是否存在跳频信号或者定频信号，若存在，则执行后续步骤；
S2.获取观测信号的短时傅里叶变换STFTx(m，n)，并计算时频谱Px(m，n)；
S3.将时频谱Px(m，n)的不同频率分量在所有时间点的功率进行累加，得到只与频率相关的平均功率谱
Figure FDA0003077454800000011
S4.在原时频谱Px(m，n)上减去对应频率的平均功率谱
Figure FDA0003077454800000012
得到时频对消谱Psub(m，n)；
S5.对所得到的时频对消谱Psub(m，n)进行归一化处理，并计算每个频点时频谱的方差d(m)和局部方差T。
2.根据权利要求1所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，所述的跳频信号模型为公式(1)
Figure FDA0003077454800000013
其中，T为总观测时间，0≤t≤T，K为观测时间内的总跳数，TH为跳周期，a(t)为基带复包络，As为跳频信号的幅度，fk和
Figure FDA0003077454800000014
分别表示第k跳的载波频率和相位，
Figure FDA0003077454800000015
3.根据权利要求2所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，所述定频信号模型为公式(2)
Figure FDA0003077454800000016
其中，AJ、fJ和
Figure FDA0003077454800000017
分别表示为定频信号的幅度、载波频率和相位。
4.根据权利要求3所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S1中，通过以下方式确定所述的观测信号：
对于给定的连续信号x(t)，其短时傅里叶变换表示为公式(3)：
Figure FDA0003077454800000021
其中h(t)是窗函数，h*(t)是它的共轭形式；
令仅存在噪声的假设为H0，噪声中存在跳频信号和定频信号的假设为H1，得到接收机截获的观测信号为公式(4)：
Figure FDA0003077454800000022
则观测信号的短时傅里叶变换表示为公式(5)：
Figure FDA0003077454800000023
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 FDA0003077454800000024
Figure FDA0003077454800000031
7.根据权利要求6所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S3中，通过以下公式(7)得到时频对消谱Psub(m，n)：
Figure FDA0003077454800000032
8.根据权利要求7所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S4中，通过以下公式(8)对所得到的时频对消谱进行归一化处理，以避免噪声方差对信号检测结果产生影响：
Figure FDA0003077454800000033
9.根据权利要求8所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S5中，通过以下公式(9)计算每个频点时频谱的方差d(m)：
Figure FDA0003077454800000034
10.根据权利要求9所述的一种基于时频谱矩阵局部方差的跳频信号检测方法，其特征在于，步骤S6中，通过以下公式(10)计算时频谱的局部方差T：
Figure FDA0003077454800000035
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 true CN113489552A (zh) 2021-10-08
CN113489552B 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)

Cited By (1)

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

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
冯维婷等: ""一种改进的跳频信号时频分析方法"", 《西安邮电大学学报》 *

Cited By (1)

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

Also Published As

Publication number Publication date
CN113489552B (zh) 2022-10-18

Similar Documents

Publication Publication Date Title
Bhargavi et al. Performance comparison of energy, matched-filter and cyclostationarity-based spectrum sensing
JP5350377B2 (ja) 異なった周波数に同調される複数の並列ｓｒ増幅器を有している超再生（ｓｒ）装置
Manesh et al. Performance evaluation of spectrum sensing techniques for cognitive radio systems
Ziafat et al. Spectrum sensing techniques for cognitive radio networks: Performance analysis
Kumar et al. Improved wavelet transform based edge detection for wide band spectrum sensing in Cognitive Radio
CN113489552B (zh) 一种基于时频谱矩阵局部方差的跳频信号检测方法
Swetha et al. Analysis of spectrum sensing based on energy detection method in cognitive radio networks
Jindal et al. Wavelet based spectrum sensing in a multipath Rayleigh fading channel
US20090063101A1 (en) Apparatus and method for detecting and identifying raman signatures
CN106154236B (zh) 一种基于cemd抑制高频地波雷达射频干扰的方法
Belgibaev et al. Software-Defined Radio Ionosonde for Diagnostics of Wideband HF Channels with the Use of USRP Platform
CN113595588B (zh) 一种基于时频谱熵的跳频信号感知方法
Ahn et al. Snapping shrimp noise detection and mitigation for underwater acoustic orthogonal frequency division multiple communication using multilayer frequency
CN108900211A (zh) 一种采用相关接收机模板设计抑制超宽带脉冲无线电干扰的方法
Imani et al. Using weighted multilevel wavelet decomposition for wideband spectrum sensing in cognitive radios
Hosseini et al. Enhanced FCME thresholding for wavelet‐based cognitive UWB over fading channels
Sansoy et al. Cyclostationary feature based detection using window method in SIMO cognitive radio system
Yeste-Ojeda et al. Cyclostationarity-based signal separation in interceptors based on a single sensor
Arjun et al. Performance Analysis of Wavelet based Spectrum Sensing and Conventional Spectrum Sensing in Fading Environment for Cognitive Radios
Avila et al. Improved method of sensing the spectrum holes in cognitive radio
Bektas et al. Energy based spectrum sensing using wavelet transform for fading channels
Hassan et al. Adaptive spectrum sensing of wireless microphones with noise uncertainty
Gautier et al. Wideband frequency domain detection using Teager-Kaiser energy operator
Kumar et al. Discrete wavelet transform based hybrid edge detection technique for wideband spectrum sensing in CRNs
Khader Enhanced performance of fh detection system using adaptive threshold level

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