CN107134281A - Adaptive filter coefficient update method during a kind of adaptive echo is eliminated - Google Patents

Adaptive filter coefficient update method during a kind of adaptive echo is eliminated Download PDF

Info

Publication number
CN107134281A
CN107134281A CN201710306285.5A CN201710306285A CN107134281A CN 107134281 A CN107134281 A CN 107134281A CN 201710306285 A CN201710306285 A CN 201710306285A CN 107134281 A CN107134281 A CN 107134281A
Authority
CN
China
Prior art keywords
mrow
adaptive
signal
frequency domain
filter coefficient
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.)
Pending
Application number
CN201710306285.5A
Other languages
Chinese (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.)
Chongqing University of Education
Original Assignee
Chongqing University of Education
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 Chongqing University of Education filed Critical Chongqing University of Education
Priority to CN201710306285.5A priority Critical patent/CN107134281A/en
Publication of CN107134281A publication Critical patent/CN107134281A/en
Pending legal-status Critical Current

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04MTELEPHONIC COMMUNICATION
    • H04M9/00Arrangements for interconnection not involving centralised switching
    • H04M9/08Two-way loud-speaking telephone systems with means for conditioning the signal, e.g. for suppressing echoes for one or both directions of traffic
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing 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/02Speech enhancement, e.g. noise reduction or echo cancellation
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing 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/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing 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/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • G10L2021/02082Noise filtering the noise being echo, reverberation of the speech

Abstract

Adaptive filter coefficient update method in being eliminated the invention discloses a kind of adaptive echo, including:Sef-adapting filter estimated echo signal y (k) is passed through according to the reference signal x (k) of pronunciation unit;Residual signal e (k) is calculated according to the signal d (k) that the y (k) and microphone are picked up;Calculate frequency domain reference signal X (f) power spectrum xPow (f);Calculate frequency domain residual signal E (f) power spectrum ePow (f);Calculate the conjugation residual signal E of the E (f)*(f) with the related power spectrum xePow (f) of the X (f);Calculate X (f) and E*(f) coefficient correlation cohxe (f);By the median cohxe in each cohxe (f)medAs sef-adapting filter in the adaptive step factor of each frequency, adaptive filter coefficient is updated.The coefficient correlation between residual signal and reference signal is updated adaptive filter coefficient by technical scheme, and effective echo cancellor can be carried out to each speech state.

Description

Adaptive filter coefficient update method during a kind of adaptive echo is eliminated
Technical field
The present invention relates to the echo processing technology in real-time speech communicating, interactive process, one kind is related specifically to Adaptive filter coefficient update method during adaptive echo is eliminated.
Background technology
During real-time speech communicating, the sound that communication terminal loudspeaker is sent always is picked up by the microphone of the terminal Get, if being sent without handling, other side can hear oneself one's voice in speech;In field of human-computer interaction, due to interaction The sound that terminal is sent is gone back by microphone pickup again, while the voice of effector has been picked up, if in microphone pickup signals In do not eliminate the sound that interactive terminal is sent, then interactive terminal is introduced into very strong dry when recognizing effector's sound of speaking Disturb, reduce the success rate of identification, ultimately cause difficult interface.Therefore, it is necessary to enter to echo during real-time speech communicating Row is adaptive to be eliminated, and it is typically what is realized using sef-adapting filter that adaptive echo, which is eliminated, in order to preferably realize echo Eliminate, it is necessary to according to actual echo disturbed condition, be updated to the adaptive filter coefficient of sef-adapting filter.
The problem of adaptive filter coefficient updates in being eliminated for adaptive echo, prior art it is also proposed some skills Art scheme, it is typical such as a kind of Chinese patent application " sef-adapting filter and filtering method for being used to eliminate acoustic echo " (Shen Please number:201410073711.1) proposed in a kind of adaptive-filtering scheme for echo cancellor.The technical scheme is a kind of Frequency domain adaptive filtering method, by the use of remote signaling x (n) as reference signal, by adaptive-filtering to propagation path h (n) System estimation is carried out, and reference signal is filtered using the systematic parameter w (n) of estimation, and then eliminates terminal Mike's elegance Echo in the signal d (n) collected obtains residual signals e (n).On the one hand e (n) is directly exported as filter result, it is another Aspect is iterated control to adaptive filter coefficient using e (n).Control is being iterated to filter coefficient in the program When, sef-adapting filter is once updated for R times per iteration to filter coefficient limitation parameter, and wherein R is set in advance solid Definite value.
The problem of existing scheme, is:On the one hand, the renewal of adaptive filter coefficient and filter effect (residual signal) Do not associate, and adjust not according to echo intensity the adaptive step factor in real time, it is impossible to ensure well it is various not With if under sound-like state echo cancellor accuracy.
The content of the invention
In order to solve the above-mentioned problems in the prior art, the present invention proposes adaptive during a kind of adaptive echo is eliminated Filter coefficient update method is answered, to realize the efficient adaptive echo cancellor under various speech states.
To achieve these goals, present invention employs following technical scheme:
Adaptive filter coefficient update method during a kind of adaptive echo is eliminated, including:
Sef-adapting filter estimated echo signal y (k) is passed through according to the reference signal x (k) of pronunciation unit;
Residual signal e (k) is calculated according to the signal d (k) that the y (k) and microphone are picked up:
E (k)=d (k)-y (k)
Calculate frequency domain reference signal X (f) power spectrum xPow (f);
Calculate frequency domain residual signal E (f) power spectrum ePow (f);
Wherein, E (f)=FFT [0M 0,e(k)];
Calculate the conjugation residual signal E of the E (f)*(f) with the related power spectrum xePow (f) of the X (f);
XePow (f)=| | X (f) E*(f)||2
Calculate X (f) and E*(f) coefficient correlation cohxe (f);
By the median cohxe in each cohxe (f)medAs sef-adapting filter each frequency adaptive step because Son, updates adaptive filter coefficient:
Wherein, the k be the moment, k=1,2 ... M, M be the sef-adapting filter length;The f is frequency;Wk+1 (f) it is the adaptive filter coefficient after frequency f renewal;Wk(f) it is the adaptive filter coefficient before frequency f renewal.
Further, the estimated echo signal y (k) includes:
The reference signal x (k) of loudspeaker unit is converted into frequency domain reference signal X (f);
X (f)=FFT [x (k-M) ..., x (k) ..., x (k+M-1)]
Estimate frequency domain echo signal Y (f),
Fast fourier inverse transformation is carried out to the Y (f) and obtains y (k);
Y (k)=IFFT [Y (f)]
Wherein, the IFFT is fast fourier inverse transformation,Operated for frequency domain convolution;W (f) works as sef-adapting filter It is preceding in adaptive filter coefficient used in frequency f.
Further, the renewal adaptive filter coefficient includes:
The frequency domain residual signal E (f) is normalized, the frequency domain residual signal E after being normalizednor (f);
According to the Enor(f) adaptive filter coefficient is updated.
Technical scheme regard the coefficient correlation between residual signal and reference signal as the adaptive step factor To realize the renewal of adaptive filter coefficient, greatly echo cancellor can be carried out to each speech state.Talking state Four kinds are generally divided into, distal end is singly said, near-end list is said, both-end is talked, both-end is silent.State is singly said distal end, microphone pickup It is echo signal, in the starting stage of sef-adapting filter, sef-adapting filter has a convergence process, residual signal and Mike The echo signal intensity of wind pickup is suitable, then reference signal and the direct correlation of residual signal are very strong, then corresponding adaptive Answer step factor also larger, when sef-adapting filter enter stable state when, residual signal is almost 0, now echo signal with The correlation of residual signal is very weak, and only adaptive filter coefficient is finely adjusted.State is said to near-end list, residual signal is The local terminal speech of microphone pickup, reference signal is 0, and the coefficient correlation now calculated is 0, without entering to sef-adapting filter Row updates, and does not have echo now actually, is not required to carry out echo cancellor really.Both-end silent status, residual signal and ginseng Examine signal all very weak, coefficient correlation also can very little.Meanwhile, the median of each frequency coefficient correlation is used in technical solution of the present invention As the adaptive step factor, correlation central tendency is embodied, the effect of adaptive-filtering can be further improved.The present invention's Technical scheme can effectively improve the accuracy of echo cancellor.
Brief description of the drawings
Adaptive filter coefficient update method flow chart during Fig. 1 adaptive echos of the present invention are eliminated.
Embodiment
In order to better illustrate technical scheme, the embodiment to the present invention is carried out below in conjunction with the accompanying drawings It is described in detail.
Specific embodiment
The present embodiment is a kind of preferred reality of adaptive filter coefficient update method during adaptive echo of the present invention is eliminated Apply mode.
Referring to Fig. 1, adaptive filter coefficient update method flow such as Fig. 1 institutes during the adaptive echo of the present embodiment is eliminated Show, including:
S1, sef-adapting filter estimated echo signal y (k) passed through according to the reference signal x (k) of pronunciation unit;
As a kind of preferred implementation scheme of the present embodiment, this step can also include through one:
S101, the reference signal x (k) of loudspeaker unit is converted into frequency domain reference signal X (f).
X (f)=FFT [x (k-M) ..., x (k) ..., x (k+M-1)]
S102, estimation frequency domain echo signal Y (f),
Wherein, W (f) is current adaptive filter coefficient;
S103, the y (k) obtained according to the Y (f):
Y (k)=IFFT [Y (f)]
Wherein, the FFT is Fourier transform, and the IFFT is inverse fourier transform;
Originally it is that this step can also be estimated echo signal by other means in example, the present invention does not have to this Limitation.
S2, the signal d (k) picked up according to the y (k) and microphone calculate residual signal e (k):
E (k)=d (k)-y (k)
S3, the power spectrum xPow (f) for calculating frequency domain reference signal X (f):
XPow (f)=| | X (f) | |2
S4, the power spectrum ePow (f) for calculating frequency domain residual signal E (f):
EPow (f)=| | E (f) | |2
Wherein,
S5, the conjugation residual signal E for calculating the E (f)*(f) with the related power spectrum xePow (f) of the X (f);
XePow (f)=| | X (f) E*(f)||2
S6, calculating X (f) and E*(f) coefficient correlation cohxe (f);
S7, by the median cohxe in each cohxe (f)medAs sef-adapting filter each frequency adaptive step Factor mu (f), updates adaptive filter coefficient;
In the present embodiment,
In the present embodiment, this step can directly update adaptive according to the X (f), the E (f) and the cohxe (f) Answer filter coefficient:
Wherein, the k be the moment, k=1,2 ... M, M be the sef-adapting filter length;The f is frequency;Wk+1 (f) it is the adaptive filter coefficient after frequency f renewal;Wk(f) it is the adaptive filter coefficient before frequency f renewal.
As a kind of preferred implementation scheme of the present embodiment, this step can also include:
S701, the frequency domain residual signal E (f) is normalized, the frequency domain residual signal after being normalized Enor(f);
S702, according to the Enor(f) adaptive filter coefficient is updated.
Wherein, the k be the moment, k=1,2 ... M, M be the sef-adapting filter length;The f is frequency;Wk+1 (f) it is the adaptive filter coefficient after frequency f renewal;Wk(f) it is the adaptive filter coefficient before frequency f renewal.
During both-end talk situation, sef-adapting filter is to need to update, and the correlation of reference signal and residual signal It is very strong, but include local terminal speech as the residual signal of feedback signal, certainly will if being directly fed back to adaptive filter coefficient Filter divergence can be caused, frequency domain residual signal is normalized first for this preferred implementation scheme, recycle normalization Rear frequency domain residual signal updates adaptive filter coefficient, the problem of solving filter divergence well.
It should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted, although with reference to compared with The present invention is described in detail good embodiment, it will be understood by those within the art that, can be to skill of the invention Art scheme is modified or equivalent substitution, and without departing from the objective and scope of technical solution of the present invention, it all should cover at this Among the right of invention.

Claims (3)

  1. Adaptive filter coefficient update method during 1. a kind of adaptive echo is eliminated, including:
    Sef-adapting filter estimated echo signal y (k) is passed through according to the reference signal x (k) of pronunciation unit;
    Residual signal e (k) is calculated according to the signal d (k) that the y (k) and microphone are picked up:
    E (k)=d (k)-y (k)
    Calculate frequency domain reference signal X (f) power spectrum xPow (f);
    Calculate frequency domain residual signal E (f) power spectrum ePow (f);
    Calculate the conjugation residual signal E of the E (f)*(f) with the related power spectrum xePow (f) of the X (f);
    Calculate X (f) and E*(f) coefficient correlation cohxe (f);
    <mrow> <mi>c</mi> <mi>o</mi> <mi>h</mi> <mi>x</mi> <mi>e</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>x</mi> <mi>e</mi> <mi>P</mi> <mi>o</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>x</mi> <mi>P</mi> <mi>o</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>e</mi> <mi>P</mi> <mi>o</mi> <mi>w</mi> <mrow> <mo>(</mo> <mi>f</mi> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
    By the median cohxe in each cohxe (f)medAs sef-adapting filter each frequency the adaptive step factor, more New adaptive filter coefficient;
    Wherein, the k be the moment, the f be frequency, k=1,2 ... M, M be the sef-adapting filter length.
  2. 2. according to the method described in claim 1, it is characterised in that the estimated echo signal y (k) includes:
    The reference signal x (k) of loudspeaker unit is converted into frequency domain reference signal X (f);
    Estimate frequency domain echo signal Y (f),
    Fast fourier inverse transformation is carried out to the Y (f) and obtains y (k);
    Wherein,Operated for frequency domain convolution;W (f) is that sef-adapting filter is current in sef-adapting filter system used in frequency f Number.
  3. 3. method according to claim 1 or 2, it is characterised in that the renewal adaptive filter coefficient includes:
    The frequency domain residual signal E (f) is normalized, the frequency domain residual signal E after being normalizednor(f);
    According to the Enor(f) adaptive filter coefficient is updated.
CN201710306285.5A 2017-05-04 2017-05-04 Adaptive filter coefficient update method during a kind of adaptive echo is eliminated Pending CN107134281A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710306285.5A CN107134281A (en) 2017-05-04 2017-05-04 Adaptive filter coefficient update method during a kind of adaptive echo is eliminated

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710306285.5A CN107134281A (en) 2017-05-04 2017-05-04 Adaptive filter coefficient update method during a kind of adaptive echo is eliminated

Publications (1)

Publication Number Publication Date
CN107134281A true CN107134281A (en) 2017-09-05

Family

ID=59715691

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710306285.5A Pending CN107134281A (en) 2017-05-04 2017-05-04 Adaptive filter coefficient update method during a kind of adaptive echo is eliminated

Country Status (1)

Country Link
CN (1) CN107134281A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108986837A (en) * 2018-09-05 2018-12-11 科大讯飞股份有限公司 A kind of filter update method and device
CN109547655A (en) * 2018-12-30 2019-03-29 广东大仓机器人科技有限公司 A kind of method of the echo cancellation process of voice-over-net call
CN109672970A (en) * 2018-11-01 2019-04-23 恒玄科技(上海)有限公司 A kind of more microphone/ear-headphone echo cancellation filter coefficient calibration methods
CN110992923A (en) * 2019-12-19 2020-04-10 浙江大华技术股份有限公司 Echo cancellation method, electronic device, and storage device
CN111294473A (en) * 2019-01-28 2020-06-16 展讯通信(上海)有限公司 Signal processing method and device
WO2020124325A1 (en) * 2018-12-17 2020-06-25 深圳市汇顶科技股份有限公司 Echo elimination adaptive filtering method, apparatus, device and storage medium
CN111355855A (en) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 Echo processing method, device, equipment and storage medium
CN111742541A (en) * 2017-12-08 2020-10-02 华为技术有限公司 Acoustic echo cancellation method and apparatus
CN112017679A (en) * 2020-08-05 2020-12-01 海尔优家智能科技(北京)有限公司 Method, device and equipment for updating adaptive filter coefficient
WO2020252629A1 (en) * 2019-06-17 2020-12-24 深圳市汇顶科技股份有限公司 Residual acoustic echo detection method, residual acoustic echo detection device, voice processing chip, and electronic device

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104883462A (en) * 2014-02-28 2015-09-02 中国科学院声学研究所 Adaptive filter and filtering method used for eliminating acoustics echo
CN105323395A (en) * 2014-07-31 2016-02-10 重庆重邮信科通信技术有限公司 Echo eliminator and echo elimination method
CN107123430A (en) * 2017-04-12 2017-09-01 广州视源电子科技股份有限公司 Echo cancel method, device, meeting flat board and computer-readable storage medium

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104883462A (en) * 2014-02-28 2015-09-02 中国科学院声学研究所 Adaptive filter and filtering method used for eliminating acoustics echo
CN105323395A (en) * 2014-07-31 2016-02-10 重庆重邮信科通信技术有限公司 Echo eliminator and echo elimination method
CN107123430A (en) * 2017-04-12 2017-09-01 广州视源电子科技股份有限公司 Echo cancel method, device, meeting flat board and computer-readable storage medium

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111742541A (en) * 2017-12-08 2020-10-02 华为技术有限公司 Acoustic echo cancellation method and apparatus
CN111742541B (en) * 2017-12-08 2021-11-30 华为技术有限公司 Acoustic echo cancellation method, acoustic echo cancellation device and storage medium
CN108986837B (en) * 2018-09-05 2021-08-17 科大讯飞股份有限公司 Filter updating method and device
CN108986837A (en) * 2018-09-05 2018-12-11 科大讯飞股份有限公司 A kind of filter update method and device
CN109672970A (en) * 2018-11-01 2019-04-23 恒玄科技(上海)有限公司 A kind of more microphone/ear-headphone echo cancellation filter coefficient calibration methods
CN111868826A (en) * 2018-12-17 2020-10-30 深圳市汇顶科技股份有限公司 Adaptive filtering method, device, equipment and storage medium in echo cancellation
WO2020124325A1 (en) * 2018-12-17 2020-06-25 深圳市汇顶科技股份有限公司 Echo elimination adaptive filtering method, apparatus, device and storage medium
CN109547655A (en) * 2018-12-30 2019-03-29 广东大仓机器人科技有限公司 A kind of method of the echo cancellation process of voice-over-net call
CN111294473A (en) * 2019-01-28 2020-06-16 展讯通信(上海)有限公司 Signal processing method and device
CN111294473B (en) * 2019-01-28 2022-01-04 展讯通信(上海)有限公司 Signal processing method and device
WO2020252629A1 (en) * 2019-06-17 2020-12-24 深圳市汇顶科技股份有限公司 Residual acoustic echo detection method, residual acoustic echo detection device, voice processing chip, and electronic device
CN110992923A (en) * 2019-12-19 2020-04-10 浙江大华技术股份有限公司 Echo cancellation method, electronic device, and storage device
CN110992923B (en) * 2019-12-19 2022-06-07 浙江大华技术股份有限公司 Echo cancellation method, electronic device, and storage device
CN111355855A (en) * 2020-03-12 2020-06-30 紫光展锐(重庆)科技有限公司 Echo processing method, device, equipment and storage medium
CN112017679A (en) * 2020-08-05 2020-12-01 海尔优家智能科技(北京)有限公司 Method, device and equipment for updating adaptive filter coefficient
CN112017679B (en) * 2020-08-05 2024-01-26 海尔优家智能科技(北京)有限公司 Method, device and equipment for updating adaptive filter coefficients

Similar Documents

Publication Publication Date Title
CN107134281A (en) Adaptive filter coefficient update method during a kind of adaptive echo is eliminated
CN109727604B (en) Frequency domain echo cancellation method for speech recognition front end and computer storage medium
CN107026950B (en) A kind of frequency domain adaptive echo cancel method
CN107123430B (en) Echo cancel method, device, meeting plate and computer storage medium
CN104883462B (en) A kind of sef-adapting filter and filtering method for eliminating acoustic echo
US20190222691A1 (en) Data driven echo cancellation and suppression
EP3791565B1 (en) Method and apparatus utilizing residual echo estimate information to derive secondary echo reduction parameters
CN107636758A (en) Acoustic echo eliminates system and method
CN107071196B (en) A kind of adaptive echo cancellation method
CN105810202B (en) It is a kind of to drop hypoechoic method, apparatus and communication apparatus
US10880440B2 (en) Echo canceller and method therefor
TW200948030A (en) Apparatus and method for computing filter coefficients for echo suppression
CN111768796A (en) Acoustic echo cancellation and dereverberation method and device
CN110992923B (en) Echo cancellation method, electronic device, and storage device
CN107483761A (en) A kind of echo suppressing method and device
CN111031448B (en) Echo cancellation method, echo cancellation device, electronic equipment and storage medium
EP4071757A1 (en) Echo cancellation method and device
Ivry et al. Nonlinear acoustic echo cancellation with deep learning
EP2930917B1 (en) Method and apparatus for updating filter coefficients of an adaptive echo canceller
EP2939405B1 (en) Method and apparatus for audio processing
US8600037B2 (en) Audio quality and double talk preservation in echo control for voice communications
CN115457928A (en) Echo cancellation method and system based on neural network double-talk detection
CN106161820A (en) A kind of interchannel decorrelation method for stereo acoustic echo canceler
CN107005268A (en) Echo eliminator and method for echo cancellation
CN111370016B (en) Echo cancellation method and electronic equipment

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170905