CN105049128B - 一种音频播放中嵌入多载波声波通信的方法 - Google Patents
一种音频播放中嵌入多载波声波通信的方法 Download PDFInfo
- Publication number
- CN105049128B CN105049128B CN201510329624.2A CN201510329624A CN105049128B CN 105049128 B CN105049128 B CN 105049128B CN 201510329624 A CN201510329624 A CN 201510329624A CN 105049128 B CN105049128 B CN 105049128B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msup
- audio signal
- original audio
- msub
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000006854 communication Effects 0.000 title claims abstract description 28
- 238000004891 communication Methods 0.000 title claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims abstract description 70
- 230000005236 sound signal Effects 0.000 claims abstract description 54
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 230000003595 spectral effect Effects 0.000 claims abstract description 10
- 238000005311 autocorrelation function Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Transmission Systems Not Characterized By The Medium Used For Transmission (AREA)
Abstract
本发明涉及一种音频播放中嵌入多载波声波通信的方法,将发送端滤波器组发送各路子载波的方法简化为相乘的IFFT和PPN矩阵;然后提取一段原始音频信号的功率谱包络信息,将原始音频信号的功率谱和经过IFFT后的各路子载波的频谱进行相乘从而对各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号并将该声波信号向接收端滤波器组进行发射;接收端滤波器组接收到声波信号后进行解调,获得原始数据。本发明提供的方法能够掩盖FBMC声波通信***通信过程中产生的噪声。
Description
技术领域
本发明涉及无线通信技术领域,更具体地,涉及一种音频播放中嵌入多载波声波通信的方法。
背景技术
近年来,无线通信技术飞速发展,滤波器组多载波(FBMC,FilterBankMulti-carrier)由于其相对OFDM技术有更高频率利用率,在感知无线电和多接入通信***中得到广泛研究和应用。
在FBMC声波通信***中,由于原始数据经过发送端编码调制以后,通过麦克风产生的是一定频率范围内的声波,在人耳听到的是让人难以接受的噪声,容易对周围环境和用户体验造成不良影响。因此,如何掩盖这种噪声成为必须要解决的问题。
发明内容
本发明为解决以上现有技术的缺陷,提供了一种音频播放中嵌入多载波声波通信的方法,用以掩盖FBMC声波通信***通信过程中产生的噪声。
为实现以上发明目的,采用的技术方案是:
一种音频播放中嵌入多载波声波通信的方法,将发送端滤波器组发送各路子载波的方法简化为相乘的IFFT和PPN矩阵;然后提取一段原始音频信号的功率谱包络信息,将原始音频信号的功率谱和经过IFFT后的各路子载波的频谱进行相乘从而对各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号并将该声波信号向接收端滤波器组进行发射;接收端滤波器组接收到声波信号后进行解调,获得原始数据。
上述方案中,本发明提供的方法通过对多载波声波通信的各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号,该声波信号能够掩盖FBMC声波通信***通信过程中产生的噪声。
优选地,将发送端滤波器组简化为相乘的IFFT和PPN矩阵的具体过程如下:
发送端滤波器组包含的滤波器的频率响应H(f)表示为:其中Z=ei2πf,hi表示脉冲响应系数,L表示脉冲响应系数hi的长度;根据扩展FFT方法,滤波器频率响应的Z变换H(Z)表示如下:
令滤波器的脉冲响应系数hi的长度为两个因子的乘积即L=kM,则Z变换表示为:
Hp(ZM)包含了频率响应的相位偏移,设发送端滤波器组第一个滤波器频率响应为B0(f),那么发送端滤波器组第k个滤波器的频率响应由B0(f)经过k/M个单位频偏得到,如下所示:
相对应的Z变换为:
用多相结构表示:
令W=e-j2π/M,得到如下所示的各路子载波的频谱矩阵:
为各路子载波的频谱矩阵,为IFFT,为多相结构PPN矩阵。
优选地,对音频信号进行线性预测分析从而提取原始音频信号的功率谱,线性预测分析过程为:计算原始音频信号的自相关函数,然后结合计算的自相关函数使用列文森-杜宾算法求得原始音频信号的AR系数a=[a(1),……,a(Na)]和增益常数σ2,即可通过下式获取原始音频信号的功率谱:
其中,E(ejΩ)表示原始音频信号的功率谱,A(ejΩ)表示AR系数a的频谱,Na为线性预测分析的阶数。
优选地,提取获得一段原始音频信号的功率谱后,使用低通滤波器去除原始音频信号功率谱中的低频部分,修改后的各路子载波的频谱与PPN矩阵相乘后再与去除低频部分的原始音频信号功率谱相加。
与现有技术相比,本发明的有益效果是:
本发明提供的方法通过对多载波声波通信的各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号,该声波信号能够掩盖FBMC声波通信***通信过程中产生的噪声。
附图说明
图1为IFFT-PPN方法的实施过程示意图。
图2为OQAM调制的示意图。
图3为多载波声波通信方法的实施过程示意图。
图4为子载波功率控制的示意图。
图5为提取音频信号功率谱的原理示意图。
图6为子载波功率控制方法的原理示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
以下结合附图和实施例对本发明做进一步的阐述。
实施例1
本发明提供的多载波声波通信的方法具体如下:
将发送端滤波器组发送各路子载波的方法简化为相乘的IFFT和PPN矩阵;然后提取一段原始音频信号的功率谱包络信息,将原始音频信号的功率谱和经过IFFT后的各路子载波的频谱进行相乘从而对各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号并将该声波信号向接收端滤波器组进行发射;接收端滤波器组接收到声波信号后进行解调,获得原始数据。
本实施例中,通过以下操作将发送端滤波器组简化为相乘的IFFT和PPN矩阵:
在时域上,原型滤波器可以定义为脉冲响应系数hi和输入输出序列的关系,如下式所示:
其中L为脉冲响应系数的长度,采样频率假设为归一化,频率响应可以表示为:
数字处理中,特别是数字滤波器一般使用Z变换来表示频率响应,令Z=ei2πf,那么上式可表示为:
根据扩展FFT方法,滤波器频率响应的Z变换H(Z)表示如下:
令滤波器的脉冲响应系数hi的长度为两个因子的乘积即L=kM,则Z变换表示为:
Hp(ZM)包含了频率响应的相位偏移,设发送端滤波器组第一个滤波器频率响应为B0(f),那么发送端滤波器组第k个滤波器的频率响应由B0(f)经过k/M个单位频偏得到,如下所示:
相对应的Z变换为:
用多相结构表示:
令W=e-j2π/M,得到如下所示的各路子载波的频谱矩阵:
为各路子载波的频谱矩阵,为IFFT,可以由IFFT实现,为多相结构PPN矩阵,这是滤波器组中所有滤波器的共有部分。
步骤S2中,对音频信号进行线性预测分析从而提取音频信号的功率谱,具体过程如图5所示,其中线形较细的曲线代表音频信号的功率谱曲线,线形较粗的曲线代表音频信号的包络曲线,线性预测分析过程为:计算原始音频信号的自相关函数,然后结合计算的自相关函数使用列文森-杜宾算法求得原始音频信号的AR系数a=[a(1),……,a(Na)]和增益常数σ2,即可通过下式获取原始音频信号的功率谱:
其中,E(ejΩ)表示原始音频信号的功率谱,A(ejΩ)表示AR系数a的频谱,Na为线性预测分析的阶数。
在具体实施的时候,如图3所示,在发送端,数据经过RS编码交织之后,再通过OQAM调制然后进行串并转换,接着通过IFFT与原始音频信号的功率谱相乘,然后经过多相网络PPN;串并转换后,乘以一个中心频率为12.5KHz的载波将数据频率范围限制在10-15KHz之间,再与原始音频信号的功率谱相加,通过扬声器发出一段接近原始音频信号功率谱的声波。在接收端,麦克风接收到声波信号后,经过带通滤波器得到数据频段的信号,然后乘以一个12.5KHz的载波进行相干解调,串并转换后,经过PPN-FFT结构,再通过并串转换,OQAM解调,RS译码恢复原始数据。
图1为IFFT-PPN方法示意图,根据上面公式推导可知发送端滤波器组可以简化为一个IFFT和PPN,比直接设计原型滤波器组的复杂度低很多。图2为OQAM调制示意图,在QAM调制信号中,虚部被延迟半个符号周期被传送。用以实现全速率通信,避免了只能用相隔子载波进行传输。
图4为子载波功率控制原理图,普通多载波信号的频谱是均匀的脉冲,提取获得一段原始音频信号的功率谱后,使用低通滤波器去除原始音频信号功率谱中的高频部分,修改后的各路子载波的频谱与PPN矩阵相乘后再与去除高频部分的原始音频信号功率谱相加。
图6为整个子载波功率控制方法原理图,原始音频的功率谱与原始信号相乘。需要注意的是,原始音频信号中有一部分可能为零,在与各路子载波的频谱相乘的时候会把各路子载波的频谱变为零,因此需要设定一个最低阈值,避免音频信号的功率谱为零,从而降低信号的误码率。这样就得到了一个接近音频信号的高频功率谱,在进入扬声器之前同低频部分相加,得到一个和原始音频近似的信号,从而掩盖噪音影响。
本发明提供的方法通过对多载波声波通信的各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号,该声波信号能够掩盖FBMC声波通信***通信过程中产生的噪声。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (3)
1.一种音频播放中嵌入多载波声波通信的方法,其特征在于:将发送端滤波器组发送各路子载波的方法简化为相乘的IFFT和PPN矩阵;然后提取一段原始音频信号的功率谱包络信息,将原始音频信号的功率谱和经过IFFT后的各路子载波的频谱进行相乘从而对各路子载波的频谱功率进行修改,修改后的各路子载波的频谱与PPN矩阵相乘后再与原始音频信号的功率谱相加,得到一段接近原始音频信号功率谱的声波信号并将该声波信号向接收端滤波器组进行发射;接收端滤波器组接收到声波信号后进行解调,获得原始数据;
对音频信号进行线性预测分析从而提取原始音频信号的功率谱,线性预测分析过程为:计算原始音频信号的自相关函数,然后结合计算的自相关函数使用列文森-杜宾算法求得原始音频信号的AR系数a=[a(1),……,a(Na)]和增益常数σ2,即可通过下式获取原始音频信号的功率谱:
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&Omega;</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<msup>
<mi>&sigma;</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>|</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&Omega;</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mo>|</mo>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&Omega;</mi>
</mrow>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mi>a</mi>
</mrow>
</munderover>
<mi>a</mi>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>v</mi>
<mi>&Omega;</mi>
</mrow>
</msup>
</mrow>
其中,E(ejΩ)表示原始音频信号的功率谱,A(ejΩ)表示AR系数a的频谱,Na为线性预测分析的阶数。
2.根据权利要求1所述的音频播放中嵌入多载波声波通信的方法,其特征在于:将发送端滤波器组简化为相乘的IFFT和PPN矩阵的具体过程如下:
发送端滤波器组包含的滤波器的频率响应H(f)表示为:其中Z=ei2πf,hi表示脉冲响应系数,L表示脉冲响应系数hi的长度;根据扩展FFT方法,滤波器频率响应的Z变换H(Z)表示如下:
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>p</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>Z</mi>
<mi>M</mi>
</msup>
<mo>)</mo>
</mrow>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mi>p</mi>
</mrow>
</msup>
</mrow>
令滤波器的脉冲响应系数hi的长度为两个因子的乘积即L=kM,则Z变换表示为:
<mrow>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>Z</mi>
<mi>M</mi>
</msup>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>K</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>h</mi>
<mrow>
<mi>k</mi>
<mi>M</mi>
<mo>+</mo>
<mi>p</mi>
</mrow>
</msub>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mi>k</mi>
<mi>M</mi>
</mrow>
</msup>
</mrow>
Hp(ZM)包含了频率响应的相位偏移,设发送端滤波器组第一个滤波器频率响应为B0(f),那么发送端滤波器组第k个滤波器的频率响应由B0(f)经过k/M个单位频偏得到,如下所示:
<mrow>
<msub>
<mi>B</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<mfrac>
<mi>k</mi>
<mi>M</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>f</mi>
<mo>-</mo>
<mi>k</mi>
<mo>/</mo>
<mi>M</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
相对应的Z变换为:
<mrow>
<msub>
<mi>B</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>Z</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>i</mi>
<mi>k</mi>
<mo>/</mo>
<mi>M</mi>
</mrow>
</msup>
<msup>
<mi>Z</mi>
<mrow>
<mo>-</mo>
<mi>i</mi>
</mrow>
</msup>
</mrow>
用多相结构表示:
令W=e-j2π/M,得到如下所示的各路子载波的频谱矩阵:
为各路子载波的频谱矩阵,为IFFT,为多相结构PPN矩阵。
3.根据权利要求1所述的音频播放中嵌入多载波声波通信的方法,其特征在于:提取获得一段原始音频信号的功率谱后,使用低通滤波器去除原始音频信号功率谱中的高频部分,修改后的各路子载波的频谱与PPN矩阵相乘后再与去除高频部分的原始音频信号功率谱相加。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329624.2A CN105049128B (zh) | 2015-06-15 | 2015-06-15 | 一种音频播放中嵌入多载波声波通信的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510329624.2A CN105049128B (zh) | 2015-06-15 | 2015-06-15 | 一种音频播放中嵌入多载波声波通信的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105049128A CN105049128A (zh) | 2015-11-11 |
CN105049128B true CN105049128B (zh) | 2017-09-19 |
Family
ID=54455359
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510329624.2A Active CN105049128B (zh) | 2015-06-15 | 2015-06-15 | 一种音频播放中嵌入多载波声波通信的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105049128B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105791197B (zh) * | 2015-12-25 | 2019-04-26 | 中国科学院上海微***与信息技术研究所 | 一种偏移正交多载波***的脉冲成形滤波器 |
CN105635022B (zh) * | 2015-12-29 | 2018-09-18 | 中国科学院上海微***与信息技术研究所 | 一种偏移正交多载波基带*** |
CN106302279A (zh) * | 2016-08-04 | 2017-01-04 | 成都极比特通信技术有限公司 | 基于干扰方差统计的fbmc***均衡方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101933242A (zh) * | 2008-08-08 | 2010-12-29 | 雅马哈株式会社 | 调制装置以及解调装置 |
CN104243388B (zh) * | 2014-09-25 | 2017-10-27 | 陈景竑 | 基于ofdm的声波通信*** |
-
2015
- 2015-06-15 CN CN201510329624.2A patent/CN105049128B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105049128A (zh) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106576092B (zh) | Fbmc信号的发送方法、接收方法和发射机以及接收机 | |
US20140362934A1 (en) | Multi transform ofdm systems and methods with low peak to average power ratio signals | |
RU2016146120A (ru) | Асинхронная связь с несколькими несущими | |
JP2017521900A5 (zh) | ||
CN108990142B (zh) | 一种低信令开销的多载波非正交多址接入***的传输方法 | |
CN105357160A (zh) | 发送参考信号的方法及装置、接收参考信号的方法及装置 | |
CN105049128B (zh) | 一种音频播放中嵌入多载波声波通信的方法 | |
EP3148141B1 (en) | Method, apparatus, and device for transmitting data | |
JP2005304040A (ja) | 直交周波数分割多重接続システムにおける高速周波数ホッピングのための送受信装置 | |
US10491445B2 (en) | Data modulation for use in multi-carrier system, demodulation method, frame generation method, and node | |
CN105187350A (zh) | 一种基于滤波多音调制的时反水声通信方法 | |
CN111884974B (zh) | 一种基于频谱分解的宽带信号射频指纹特征提取方法 | |
Bandari et al. | PAPR analysis of wavelet based multitaper GFDM system | |
CN107438041B (zh) | 一种发送信号和接收信号的方法及装置 | |
CN106797361A (zh) | 循环卷积的***和方法 | |
RU2599930C1 (ru) | Способ и устройство восстановления передаваемой информации при ортогональном частотном уплотнении | |
US9893923B2 (en) | Method for transmitting and receiving QAM signal in filter bank-based multicarrier communication system, and apparatus therefor | |
US20160112235A1 (en) | Signal modulation method, signal demodulation method, signal modulation apparatus, signal demodulation apparatus and signal transmission system | |
JP6157722B2 (ja) | 受信装置 | |
CN116566780A (zh) | 基于稀疏码多址接入的多用户多载波dcsk通信方法及*** | |
RU2702258C1 (ru) | Способ ортогонального частотного уплотнения и его реализующее устройство | |
CN112202697B (zh) | 一种信号的处理方法、装置、存储介质以及电子装置 | |
WO2014107973A1 (zh) | 一种lte***的基带消峰方法及装置 | |
CN106656892B (zh) | 发送数据的方法和设备 | |
CN109768823B (zh) | 一种窄带多通道卫星通信*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |