CN112151047B - 一种应用于语音数字信号的实时自动增益控制方法 - Google Patents
一种应用于语音数字信号的实时自动增益控制方法 Download PDFInfo
- Publication number
- CN112151047B CN112151047B CN202011030786.3A CN202011030786A CN112151047B CN 112151047 B CN112151047 B CN 112151047B CN 202011030786 A CN202011030786 A CN 202011030786A CN 112151047 B CN112151047 B CN 112151047B
- Authority
- CN
- China
- Prior art keywords
- gain
- voice
- input
- time
- smoothing
- 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
Images
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
- G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
- G10L19/08—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters
- G10L19/083—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters the excitation function being an excitation gain
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
- G10L21/00—Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Quality & Reliability (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种应用于语音数字信号的实时自动增益控制方法,该方法包括:使用平滑滤波器求语音数字信号的近似上包络;将包络作为NLMS自适应滤波器的输入,NLMS的目标参考值根据输入信号为语音或非语音分别设置,将NLMS的自适应权重作为快速增益输出;若输入信号长时间非语音且时间超过预设的阀值,系统将切换为静默模式,快速增益设为1;若快速增益超过预设的最大增益值,则将快速增益设置为最大增益值;通过平滑滤波器将快速增益处理成平滑增益;使用平滑增益乘以输出语音信号得到输出语音信号。该方法具有快速、有效、实时的特点,可应用在实时语音通话,网络视频会议等场景,能有效解决语音忽大忽小的问题,改善用户的使用体验。
Description
技术领域
本发明涉及语音信号处理领域,具体是一种应用于语音数字信号的实时自动增益控制方法。
背景技术
随着计算机和互联网技术的发展,实时语音通话和视频通话的应用越来越广泛,音视频技术也越来越得到重视。语音自动增益控制(Automatic Gain Control简称AGC)是音视频技术中音频处理的重要环节,语音自动增益控制主要目标是解决语音通话时声音忽大忽小的情况,语音自动增益控制处理效果将直接影响着用户的体验。近些年来,出现了许多实用的语音实时自动增益控制方法,他们大多都是基于语音的峰值比较和门限比较的方法,这类方法的增益响应速度往往不够快,对音量极小的语音信号往往增益不够明显。
发明内容
本发明的目的是针对现有技术的不足,而提供一种应用于语音数字信号的实时自动增益控制方法。这种方法能自适应调整语音的增益、增益小音量的语音、减小语音音量的动态范围,使得语音的音量变得相对统一,从而提高用户的听感体验。
实现本发明目的的技术方案是:
一种应用于语音数字信号的实时自动增益控制方法,与现有技术不同处在于,所述方法是采用归一化最小均方(Normalized Least Mean Square简称NLMS)自适应滤波器为核心控制单元,通过增益控制信息、包络提取和增益平滑滤波器控制整个增益过程,包括如下步骤:
1)在语音流的处理过程中,接收到经过降噪、消回声和去混响处理过后的数字语音信号x(n)和语音活动检测VAD(Voice Activity Detection简称VAD)的结果v(n),每次处理的语音流以采样点为单位或以帧为单位,若语音流处理单位为样点,则系统的输入语音为x(n)=x(t)、VAD信息为v(n)=v(t),n代表当前时刻,t代表当前采样点;若语音流处理单位为帧,则输入语音和VAD信息分别为:
x(n)=max[X(l)],
v(n)=v(l),
其中,l代表当前帧,X(l)表示当前帧的采样点向量,max表示取向量的最大值;
2)对步骤1)接受到的语音信号x(n)做绝对值处理,并求得语音信号x(n)近似包络xenv(n);
3)将步骤2)中近似包络xenv(n)作为NLMS自适应滤波器的输入,根据步骤1)得到的VAD结果v(n)选择NLMS的目标参考值,将NLMS的自适应权重作为系统的快速增益gf(n),包括如下过程:
3-1)将步骤2)中得到的近似包络作为NLMS自适应滤波器的输入,根据步骤1)中接收到的VAD结果v(n)选取NLMS的参考值Ref,假设语音采样点的值为-1到1之间,当语音活动时,参考值Ref取0.001;当非语音时,参考值Ref取0.5;
3-2)所述NLMS自适应滤波器为一阶自适应滤波器,滤波器权重gf(n)即为系统的快速增益,快速增益的迭代公式为:
gf(n)=gf(n-1){1+μxenv(n)[Ref-yf(n-1)]},
yf(n)=gf(n)xenv(n),
步长根据自适应滤波器的输入的大小自适应调节,表达式为:
μ=μ0/(r+xenv(n)),
其中,初始步长μ0根据输入单位的采样点数设置,当以采样点输入时,μ0取0.1,r为避免分子过小而设置的,取0.0001;
4)根据步骤1)中的VAD信息和预设的最大增益阀值,调整步骤3)中系统的快速增益gf(n),过程为:
4-1)系统设置一个静默的限制时间Tl,根据步骤1)中接收的VAD结果v(n),若连续非语音时间超过Tl,则系统将切换为静默模式,Tl设置为1秒;
4-2)若系统处于静默模式,则将步骤3)中的快速增益gf(n)将直接设为1;
4-3)若系统不处于静默模式,且步骤3)得到的快速增益gf(n)大于最大增益阀值gmax,则将快速增益gf(n)设置为gmax;
5)将步骤4)中调整后的快速增益gf(n)输入自适应平滑滤波器得到系统的平滑增益g(n),自适应平滑滤波器的递归表达式为:
g(n)=αsg(n-1)+(1-αs)gf(n),
自适应平滑滤波器的平滑时间ts设置的越小,增益的变化就越快,底噪变化也会变快,使得听感变得突兀,ts根据底噪的大小,设置在0.01秒到0.5秒之间;
6)将步骤5)中得到的系统平滑增益g(n)和步骤1)中的输入语音信号x(n)相乘,得到输出语音信号y(n),表达式如下:
y(n)=x(n)*g(n),
若输入单位不是采样点而是帧,则整帧乘以g(n)并输出;
7)重复上述步骤1)至步骤6),将新的语音连续输入,最终实现实时的自动增益控制,在第一次运行时,需要初始值,将步骤2)中的近似包络xenv初值设为0、步骤3)的快速增益gf(n)初值设为1、步骤5)的平滑增益gf(n)初值设为1。
步骤2)中所述的求得语音信号x(n)近似包络xenv(n)的过程为:先对步骤1)接收到的语音信号x(n)先求绝对值xabs(n),然后采用Attack/Release滤波器求得近似包络xenv(n),公式如下:
其中,n-1代表上一时刻,上升沿即Attack和下降沿即Release分别采用不同的滤波因子αEA和αER,滤波因子由平滑时间获得:
上升平滑时间tEA和下降平滑时间tER分别取0.001秒和1秒,fs为输入语音单位的采样率,若输入为采样点,则取值为语音采样率Fs;若输入为帧,则取值为:fs=Fs/L。
本技术方案能有效地解决语音通话中声音忽大忽小的问题,处理过后的语音听感提升明显,相较于其他现有的技术,本技术方案对极小音量的语音信号的增益效果更加明显,增益的自适应速度更快,此外,增益能根据预设的平滑时间平滑,使得系统当输入信号底噪较高时也能保持良好的听感。
这种方法能自适应调整语音的增益、增益小音量的语音、减小语音音量的动态范围,使得语音的音量变得相对统一,从而提高用户的听感体验。
附图说明
图1为实施例的原理流程框图;
图2为实施例中输入语音及输入语音的近似包络示意图;
图3为实施例中快速增益和输入VAD信息示意图;
图4为实施例中平滑增益示意图;
图5为实施例中处理过后的输出语音示意图。
具体实施方式
下面结合附图和实施例对本发明的内容做进一步阐述,但不是对本发明的限定。
实施例:
参照图1,一种应用于语音数字信号的实时自动增益控制方法,所述方法是采用NLMS自适应滤波器为核心控制单元,通过增益控制信息、包络提取和增益平滑滤波器控制整个增益过程,包括如下步骤:
1)在语音流的处理过程中,接收到经过降噪、消回声和去混响处理过后的数字语音信号x(n)和语音活动检测(Voice Activity Detection简称VAD)的结果v(n),每次处理的语音流以采样点为单位或以帧为单位,若语音流处理单位为样点,则系统的输入语音为x(n)=x(t)、VAD信息为v(n)=v(t),n代表当前时刻,t代表当前采样点;若语音流处理单位为帧,则输入语音和VAD信息分别为:
x(n)=max[X(l)],
v(n)=v(l),
其中,l代表当前帧,X(l)表示当前帧的采样点向量,max表示取向量的最大值;
2)对步骤1)接受到的语音信号x(n)做绝对值处理,并求得语音信号x(n)近似包络xenv(n);
3)将步骤2)中近似包络xenv(n)作为NLMS自适应滤波器的输入,根据步骤1)得到的VAD结果v(n)选择NLMS的目标参考值,将NLMS的自适应权重作为系统的快速增益gf(n),包括如下过程:
3-1)将步骤2)中得到的近似包络作为NLMS自适应滤波器的输入,根据步骤1)中接收到的VAD结果v(n)选取NLMS的参考值Ref,假设语音采样点的值为-1到1之间,当语音活动时,参考值Ref取0.001;当非语音时,参考值Ref取0.5;
3-2)所述NLMS自适应滤波器为一阶自适应滤波器,滤波器权重gf(n)即为系统的快速增益,快速增益的迭代公式为:
gf(n)=gf(n-1){1+μxenv(n)[Ref-yf(n-1)]},
yf(n)=gf(n)xenv(n),
步长根据自适应滤波器的输入的大小自适应调节,表达式为:
μ=μ0/(r+xenv(n)),
其中,初始步长μ0根据输入单位的采样点数设置,当以采样点输入时,μ0取0.1,r为避免分子过小而设置的,取0.0001;
4)根据步骤1)中的VAD信息和预设的最大增益阀值,调整步骤3)中系统的快速增益gf(n),过程为:
4-1)系统设置一个静默的限制时间Tl,根据步骤1)中接收的VAD结果v(n),若连续非语音时间超过Tl,则系统将切换为静默模式,Tl设置为1秒;
4-2)若系统处于静默模式,则将步骤3)中的快速增益gf(n)将直接设为1;
4-3)若系统不处于静默模式,且步骤3)得到的快速增益gf(n)大于最大增益阀值gmax,则将快速增益gf(n)设置为gmax;
5)将步骤4)中调整后的快速增益gf(n)输入自适应平滑滤波器得到系统的平滑增益g(n),自适应平滑滤波器的递归表达式为:
g(n)=αsg(n-1)+(1-αs)gf(n),
自适应平滑滤波器的平滑时间ts设置的越小,增益的变化就越快,底噪变化也会变快,使得听感变得突兀,ts根据底噪的大小,设置在0.01秒到0.5秒之间;
6)将步骤5)中得到的系统平滑增益g(n)和步骤1)中的输入语音信号x(n)相乘,得到输出语音信号y(n),表达式如下:
y(n)=x(n)*g(n),
若输入单位不是采样点而是帧,则整帧乘以g(n)并输出;
7)重复上述步骤1)至步骤6),将新的语音连续输入,最终实现实时的自动增益控制,在第一次运行时,需要初始值,将步骤2)中的近似包络xenv初值设为0、步骤3)的快速增益gf(n)初值设为1、步骤5)的平滑增益gf(n)初值设为1。
步骤2)中所述的求得语音信号x(n)近似包络xenv(n)的过程为:先对步骤1)接收到的语音信号x(n)先求绝对值xabs(n),然后采用Attack/Release滤波器求得近似包络xenv(n),公式如下:
其中,n-1代表上一时刻,上升沿即Attack和下降沿即Release分别采用不同的滤波因子αEA和αER,滤波因子由平滑时间获得:
上升平滑时间tEA和下降平滑时间tER分别取0.001秒和1秒,fs为输入语音单位的采样率,若输入为采样点,则取值为语音采样率Fs;若输入为帧,则取值为:fs=Fs/L。
如图2-图5所示,由图2和图5的实验变化结果可以明显看出,本例方法能够对忽大忽小的语音自适应的调整增益,使输入语音能够增益为一段音量相对统一的语音并输出,另外,由图3可以看出,滤波过后的增益曲线相对平滑,这使得本例方法即使是在有少量底噪的情况下,最终输出的语音也能保证良好的听感。
Claims (2)
1.一种应用于语音数字信号的实时自动增益控制方法,其特征在于,所述方法是采用NLMS自适应滤波器为核心控制单元,通过增益控制信息、包络提取和增益平滑滤波器控制整个增益过程,包括如下步骤:
1)在语音流的处理过程中,接收到经过降噪、消回声和去混响处理过后的数字语音信号x(n)和语音活动检测VAD的结果v(n),每次处理的语音流以采样点为单位或以帧为单位,若语音流处理单位为样点,则系统的输入语音为x(n)=x(t)、VAD信息为v(n)=v(t),n代表当前时刻,t代表当前采样点;若语音流处理单位为帧,则输入语音和VAD信息分别为:
x(n)=max[X(l)],
v(n)=v(l),
其中,l代表当前帧,X(l)表示当前帧的采样点向量,max表示取向量的最大值;
2)对步骤1)接受到的语音信号x(n)做绝对值处理,并求得语音信号x(n)近似包络xenv(n);
3)将步骤2)中近似包络xenv(n)作为NLMS自适应滤波器的输入,根据步骤1)得到的VAD结果v(n)选择NLMS的目标参考值,将NLMS的自适应权重作为系统的快速增益gf(n),包括如下过程:
3-1)将步骤2)中得到的近似包络作为NLMS自适应滤波器的输入,根据步骤1)中接收到的VAD结果v(n)选取NLMS的参考值Ref,假设语音采样点的值为-1到1之间,当语音活动时,参考值Ref取0.001;当非语音时,参考值Ref取0.5;
3-2)所述NLMS自适应滤波器为一阶自适应滤波器,滤波器权重gf(n)即为系统的快速增益,快速增益的迭代公式为:
gf(n)=gf(n-1){1+μxenv(n)[Ref-yf(n-1)]},
yf(n)=gf(n)xenv(n),
步长根据自适应滤波器的输入的大小自适应调节,表达式为:
μ=μ0/(r+xenv(n)),
其中,初始步长μ0根据输入单位的采样点数设置,当以采样点输入时,μ0取0.1,r为避免分子过小而设置的,取0.0001;
4)根据步骤1)中的VAD信息和预设的最大增益阀值,调整步骤3)中系统的快速增益gf(n),过程为:
4-1)系统设置一个静默的限制时间Tl,根据步骤1)中接收的VAD结果v(n),若连续非语音时间超过Tl,则系统将切换为静默模式,Tl设置为1秒;
4-2)若系统处于静默模式,则将步骤3)中的快速增益gf(n)将直接设为1;
4-3)若系统不处于静默模式,且步骤3)得到的快速增益gf(n)大于最大增益阀值gmax,则将快速增益gf(n)设置为gmax;
5)将步骤4)中调整后的快速增益gf(n)输入自适应平滑滤波器得到系统的平滑增益g(n),自适应平滑滤波器的递归表达式为:
g(n)=αsg(n-1)+(1-αs)gf(n),
其中,fs为输入语音单位的采样率,自适应平滑滤波器的平滑时间ts设置的越小,增益的变化就越快,底噪变化也会变快,使得听感变得突兀,ts根据底噪的大小,设置在0.01秒到0.5秒之间;
6)将步骤5)中得到的系统平滑增益g(n)和步骤1)中的输入语音信号x(n)相乘,得到输出语音信号y(n),表达式如下:
y(n)=x(n)*g(n),
若语音流处理单位不是采样点而是帧,则整帧的采样点乘以g(n)并输出;
7)重复上述步骤1)至步骤6),将新的语音连续输入,最终实现实时的自动增益控制,在第一次运行时,需要初始值,将步骤2)中的近似包络xenv初值设为0、步骤3)的快速增益gf(n)初值设为1、步骤5)的平滑增益gf(n)初值设为1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011030786.3A CN112151047B (zh) | 2020-09-27 | 2020-09-27 | 一种应用于语音数字信号的实时自动增益控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011030786.3A CN112151047B (zh) | 2020-09-27 | 2020-09-27 | 一种应用于语音数字信号的实时自动增益控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112151047A CN112151047A (zh) | 2020-12-29 |
CN112151047B true CN112151047B (zh) | 2022-08-05 |
Family
ID=73894262
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011030786.3A Active CN112151047B (zh) | 2020-09-27 | 2020-09-27 | 一种应用于语音数字信号的实时自动增益控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112151047B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1809105A (zh) * | 2006-01-13 | 2006-07-26 | 北京中星微电子有限公司 | 适用于小型移动通信设备的双麦克语音增强方法及系统 |
CN101009099A (zh) * | 2007-01-26 | 2007-08-01 | 北京中星微电子有限公司 | 数字自动增益控制方法及装置 |
WO2013124712A1 (en) * | 2012-02-24 | 2013-08-29 | Nokia Corporation | Noise adaptive post filtering |
CN105322904A (zh) * | 2015-09-08 | 2016-02-10 | 湖南国科微电子股份有限公司 | 一种数字音频信号的电平控制方法 |
CN106782504A (zh) * | 2016-12-29 | 2017-05-31 | 百度在线网络技术(北京)有限公司 | 语音识别方法和装置 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1238489B1 (en) * | 1999-12-13 | 2008-03-05 | Broadcom Corporation | Voice gateway with downstream voice synchronization |
US9173025B2 (en) * | 2012-02-08 | 2015-10-27 | Dolby Laboratories Licensing Corporation | Combined suppression of noise, echo, and out-of-location signals |
EP2806424A1 (en) * | 2013-05-20 | 2014-11-26 | ST-Ericsson SA | Improved noise reduction |
GB201501791D0 (en) * | 2015-02-03 | 2015-03-18 | Microsoft Technology Licensing Llc | Non-linear echo path detection |
-
2020
- 2020-09-27 CN CN202011030786.3A patent/CN112151047B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1809105A (zh) * | 2006-01-13 | 2006-07-26 | 北京中星微电子有限公司 | 适用于小型移动通信设备的双麦克语音增强方法及系统 |
CN101009099A (zh) * | 2007-01-26 | 2007-08-01 | 北京中星微电子有限公司 | 数字自动增益控制方法及装置 |
WO2013124712A1 (en) * | 2012-02-24 | 2013-08-29 | Nokia Corporation | Noise adaptive post filtering |
CN105322904A (zh) * | 2015-09-08 | 2016-02-10 | 湖南国科微电子股份有限公司 | 一种数字音频信号的电平控制方法 |
CN106782504A (zh) * | 2016-12-29 | 2017-05-31 | 百度在线网络技术(北京)有限公司 | 语音识别方法和装置 |
Non-Patent Citations (4)
Title |
---|
Hands-free voice communication with TV;Istvan I. Papp等;《IEEE Transactions on Consumer Electronics》;20110722;第57卷(第2期);全文 * |
Noise reduction and echo cancellation front-end for speech codecs;Filiz Basbug等;《IEEE TRANSACTIONS ON SPEECH AND AUDIO PROCESSING》;20030131;第11卷(第1期);全文 * |
带VoIP的ZigBee电话的语音处理技术研究;黎水芬;《中国优秀硕士学位论文全文数据库》;20070315;全文 * |
高清视频通信云服务若干关键技术研究;李铎等;《电信工程技术与标准化》;20170815(第08期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112151047A (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8170879B2 (en) | Periodic signal enhancement system | |
JP3484801B2 (ja) | 音声信号の雑音低減方法及び装置 | |
KR101120970B1 (ko) | 이동 오디오 디바이스에 대한 자동 음량 및 다이나믹 레인지 조절 | |
US7610196B2 (en) | Periodic signal enhancement system | |
US6023674A (en) | Non-parametric voice activity detection | |
JP5265056B2 (ja) | 雑音抑圧装置 | |
EP2898510B1 (en) | Method, system and computer program for adaptive control of gain applied to an audio signal | |
US20080008311A1 (en) | Methods and apparatus for improving voice quality in an environment with noise | |
US20060089958A1 (en) | Periodic signal enhancement system | |
JPH09503590A (ja) | 会話の品質向上のための背景雑音の低減 | |
JPH07306695A (ja) | 音声信号の雑音低減方法及び雑音区間検出方法 | |
CN112242147B (zh) | 一种语音增益控制方法及计算机存储介质 | |
WO2008104446A2 (en) | Method for reducing noise in an input signal of a hearing device as well as a hearing device | |
CN108133712B (zh) | 一种处理音频数据的方法和装置 | |
KR20160113224A (ko) | 오디오 신호를 압축하는 오디오 압축 시스템 | |
CN111199751B (zh) | 一种麦克风的屏蔽方法、装置和电子设备 | |
JP4843691B2 (ja) | 信号特性変化装置 | |
US20200279575A1 (en) | Automatic gain control for speech recognition engine in far field voice user interface | |
US8639294B2 (en) | System and method for performing automatic gain control in mobile phone environments | |
WO2017196382A1 (en) | Enhanced de-esser for in-car communication systems | |
CN112151047B (zh) | 一种应用于语音数字信号的实时自动增益控制方法 | |
US9614486B1 (en) | Adaptive gain control | |
CN111243631B (zh) | 一种自动增益控制方法及电子设备 | |
CA2524162C (en) | Periodic signal enhancement system | |
JP2001236090A (ja) | 音声入力装置 |
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 |