EP0774750B1 - Détermination des fréquences du spectre de raies pour utilisation dans un radiotéléphone - Google Patents

Détermination des fréquences du spectre de raies pour utilisation dans un radiotéléphone Download PDF

Info

Publication number
EP0774750B1
EP0774750B1 EP19960308081 EP96308081A EP0774750B1 EP 0774750 B1 EP0774750 B1 EP 0774750B1 EP 19960308081 EP19960308081 EP 19960308081 EP 96308081 A EP96308081 A EP 96308081A EP 0774750 B1 EP0774750 B1 EP 0774750B1
Authority
EP
European Patent Office
Prior art keywords
polynomials
lpc
zero
zeros
polynomial
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.)
Expired - Lifetime
Application number
EP19960308081
Other languages
German (de)
English (en)
Other versions
EP0774750A2 (fr
EP0774750A3 (fr
Inventor
Vesa T. Ruoppila
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.)
Nokia Oyj
Original Assignee
Nokia Oyj
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 Nokia Oyj filed Critical Nokia Oyj
Publication of EP0774750A2 publication Critical patent/EP0774750A2/fr
Publication of EP0774750A3 publication Critical patent/EP0774750A3/fr
Application granted granted Critical
Publication of EP0774750B1 publication Critical patent/EP0774750B1/fr
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; 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
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/24—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being the cepstrum

Definitions

  • This invention relates generally to speech encoding methods and apparatus and, in particular, to linear predictive coding (LPC) speech and audio coding techniques that employ a line spectrum frequency representation of an LPC filter.
  • LPC linear predictive coding
  • Linear predictive coding is a known technique for analyzing a speech signal and for characterizing the signal in terms of coefficients which are encoded, broadcast, received and decoded to recover an approximation of the original signal.
  • the parameters of a LPC filter are coded and sent as a part of the information stream.
  • the use of line spectrum frequencies is an alternative to the use of, for example, polynomial coefficients or reflection coefficients for representing the LPC filter.
  • the line spectrum frequencies have useful properties for quantization and interpolation which make them a more attractive representation than polynomial or reflection coefficients.
  • a method for determining line spectrum frequencies of a LPC filter is disclosed.
  • the predictor polynominal of the LPC filter is decomposed into symmetric and antisymmetric auxiliary polynomials, the zeros of which determine the line spectrum frequencies of the filter.
  • the line spectrum frequency representation is determined by solving the zeros of the two auxiliary polynomials. Due to the symmetry of the auxiliary polynomials, their zeros are preferably solved from two cosine series. In speech codecs this is usually done by a bisection algorithm, and by employing the definition of Chebyshev polynomials to evaluate the cosine series.
  • a potential aim of this invention is to reduce the complexity of the line spectrum frequency computation. This may be obtained by rewriting the cosine series as polynomials using explicit forms of Chebyshev polynomials. This enables an evaluation of the series by nested multiplications. Moreover, the already-computed zeros are successively eliminated from the polynomial by polynomial deflation. This procedure and the properties of the auxiliary polynomials can enable the initial values to be chosen in the zero finding algorithm such that a zero is found by only a few iterations using the zero of the other polynomial. Thus, the invention may educe considerably the arithmetic operations required to compute the line spectrum frequencies. The method of this invention thus has the potential to be implemented with relatively low complexity, and furthermore to be accomplished using fixed-point arithmetic.
  • FIG. 1 A simplified block diagram of a speech encoder 10 employing the spectrum frequency representation of the LPC filter is presented in Fig. 1.
  • the speech encoder 10 may form a portion of a radiotelephone, such as a digital cellular user terminal or a personal communicator device.
  • An input audio signal such as a speech signal obtained from a speech transducer or microphone 5, is converted into a digital form by an analog-to-digital (A/D) converter 12.
  • A/D analog-to-digital
  • the digital output of the A/D converter is preprocessed by separating the signal into frames of convenient length, typically of the order of tens of milliseconds. It should be noted that the A/D conversion is not necessary if the signal is already in digital form. After preprocessing, the signal is applied to an LPC analysis block 14.
  • the LPC analysis produces coefficients for an LPC filter, also referred to herein as an LPC-analysis filter 16.
  • the output of the LPC analysis block 14 is transformed into a line spectrum frequency (LSF) representation in block 18.
  • LSF line spectrum frequency
  • the LSF coefficients may be quantized in block 20 and then interpolated in block 22 in order to construct a LPC analysis filter for each speech subframe.
  • the analysis filter is constructed separately for each subframe.
  • the residual is encoded in the excitation encoding block 26, and the resulting encoded residual, i.e., an encoded excitation signal, is multiplexed (block 28) with the quantized LSF coefficients into a bit stream transmitted to a speech decoder (not shown) via a communication channel 30.
  • the communication channel 30 is a radio channel linking the mobile terminal to a base station (not shown) by a transmitter 32 and an antenna 34.
  • the "side information" input to the multiplexer 28 determines, for example, the operational mode of the speech coder, particularly for variable rate codecs such as QCELP. For a speech coder operated in a fixed rate mode, this input may not be required.
  • a n+1 (z) A n (z) + k n+1 z -n-1 A n (z -1 ), wherein k 1 , k 2 ,...,k n+1 are reflection coefficients.
  • the recurrence formula (2) is called the Levinson-Durbin solution to the Yule-Walker equations. It expresses the relationship between the (n+1)th and the nth degree predictor polynomials. For the purposes of this description it is assumed that all roots of the predictor polynomial A n (z) are inside the unit circle, i.e., that the predictor polynomial is minimum phase.
  • P n+1 (z) and Q n+1 (z) can be factored as follows: where ⁇ 1 , ⁇ 2 ,..., ⁇ n are the phase angles of the zeros of the polynomials such that (10) 0 ⁇ ⁇ 1 ⁇ ⁇ 2 ⁇ ⁇ n ⁇ ⁇ ⁇ 1 , ⁇ 2 ,..., ⁇ n are the line spectrum frequencies of A n (z).
  • the degree of the polynomial Q(z) is 2m Q , where
  • the symmetric polynomial Q(z) can be rewritten similarly to the form where q 1 , q 2 ,..., q mQ are the coefficients of the polynomial Q(z).
  • Kang and Fransen (“Application of line spectrum pairs to low bit rate speech encoders", Proceedings of IEEE International Conference on Acoustics, Speech, and Signal Processing, Tampa, FL, pp. 7.3.1-7.3.4, March 1984) have proposed autocorrelation based and ratio-filter based methods for finding line spectrum frequencies. However, these two methods also require an evaluation of a large number of trigonometric functions.
  • Saoudi et al. ("A new efficient algorithm to compute the LSP parameters for speech coding", Signal Processing, vol. 28, pp. 201-212, 1992) reformulated the problem. That is, they introduced an algorithm which solves line spectrum frequencies from eigenvalues of tridiagonal symmetric matrices without computing the predictor polynomial. The eigenvalues are computed by the bisection method. Saoudi et al. also compared the complexity of several algorithms introduced for computing the line spectrum frequencies. Unfortunately, the results of the comparison are not comprehensive.
  • the algorithm (22) is known to be used in speech codecs, it has several shortcomings. For example, if some zeros of the polynomials are sufficiently near each other, the algorithm has a tendency to miss zeros, since the sign change is not detected. This shortcoming can be circumvented by making the grid denser. In other words, the intervals are made sufficiently small so that two or more roots do not occur in the same interval. The definitive choice of the grid interval is always a compromise between a reliability and computational burden. That is, as the grid interval is made more dense, the overall processing burden increases as well.
  • the method successively eliminates already-found zeros from the polynomials.
  • the procedure is known as polynomial deflation, or as synthetic division in numerical analysis, see, e.g., Kincaid and Chaney.
  • polynomial deflation allows for the elimination of the sign-change detection procedure from the algorithm.
  • Another advantage is that efficient algorithms (which have better convergence properties than the bisection method) can be exploited together with the property (21) to locate the zeros. For example, when zeros are computed by Newton's method, the algorithm is guaranteed to find all line spectrum frequencies if sufficient numerical accuracy is used.
  • step (d) of algorithm (29) i.e., in the block Solve the last zero in Fig. 2.
  • the initial values can be chosen differently.
  • Figs. 4A-4D illustrate the progress of the technique (algorithm (29)) employing polynomial deflation.
  • the zeros of the polynomials R(x) and S(x) corresponding to A 10 (z) are presented in Fig. 4A by crosses (x) and circles (o), respectively.
  • (x - x 1 ) is divided from the polynomial R(x).
  • the polynomial R(x) after the first deflation is shown in Fig. 4C.
  • the search is continued from S(x) using x 1 as an initial value.
  • Fig. 4D illustrates the polynomial R(x) after two deflations.
  • the is invention teaches in one aspect a method for determining line spectrum frequencies of a linear predictive coder (LPC) filter that is expressed as symmetric and antisymmetric polynomials, the zeros of which determine the line spectrum frequencies of the LPC filter.
  • the method includes the steps of (a) expressing the polynomials using explicit forms of Chebyshev polynomials; (b) iteratively solving a zero of a first of the polynomials using a zero of the other one of the polynomials; and (c) successively eliminating zeros from the polynomials by polynomial division so as to determine the line spectrum frequencies.
  • the first polynomial can be either the symmetric or the antisymmetric polynomial, in which case the other polynomial is then the antisymmetric or the symmetric polynomial, respectively.
  • any of the zeros of the first polynomial can be solved using any of the zeros of the other polynomial.
  • the immittance spectrum frequency representation of the LPC filter is based on a similar polynomial decomposition of a predictor polynomial as the line spectrum frequency representation.
  • the symmetric polynomial P n (z) and the antisymmetric polynomial Q n (z) also have similar properties as polynomials (3) and (4).
  • the roots of P n (z) and Q n (z) are on the unit circle, and they are simple and separate from each other.
  • phase angles ⁇ 1 , ⁇ 2 , ..., ⁇ n-1 and the parameter k give an unique parametrization for the LPC filter.
  • the properties of the immittance spectrum frequencies and their relation to line spectrum frequencies have been discussed in more detail by Bistritz and Peller.
  • polynomials (35) and (36) are not identical to the polynomials defined in (3) and (4).
  • the degrees of the symmetric polynomials P(z) and Q(z) defined in (35) and (36) are 2m p and 2m Q , respectively, where and
  • Equations (38) and (39) can be written to the form (25) and (26) by proceeding as described above.
  • the immittance spectrum frequencies ⁇ 1 , ⁇ 2 , ..., ⁇ n-1 can be solved in accordance with this invention by the techniques summarized in the algorithm (29).
  • the audio encoder can be used in a PC or workstation connected to a network.
  • the communication channel 30 may then be a wired network (e.g., Internet).

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Computational Linguistics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Noise Elimination (AREA)

Claims (11)

  1. Procédé pour déterminer des fréquences du spectre de raies d'un filtre de codeur à prédiction linéaire (CPL) exprimées sous la forme de polynômes symétrique et antisymétrique, dont les zéros déterminent les fréquences du spectre de raies du filtre CPL, comprenant les étapes consistant à:
    exprimer des polynômes utilisant des formes explicites de polynômes de Tchebyshev;
    résoudre de manière itérative un zéro d'un premier des polynômes en utilisant un zéro de l'autre des polynômes; et
    éliminer successivement des zéros des polynômes par division de polynôme afin de déterminer les fréquences du spectre de raies.
  2. Procédé selon la revendication 1, comprenant en outre une étape consistant à transmettre un signal à codage CPL à une voie de communication.
  3. Procédé selon la revendication 2, dans lequel l'étape de transmission transmet le signal à codage CPL à une voie de communication radioélectrique.
  4. Procédé selon l'une quelconque des revendications précédentes, dans lequel le zéro de l'autre des polynômes est utilisé comme valeur initiale lors de la solution itérative du zéro du premier des polynômes.
  5. Station mobile permettant des communications sans fil sur une voie de communication, ladite station mobile comprenant un transducteur de paroles pour délivrer un signal de paroles et comprenant en outre:
    un codeur à prédiction linéaire (CPL) dont une entrée est couplée au signal de parole et une sortie est couplée à la voie de communication; ledit CPL comportant
    un filtre CPL ayant une première entrée couplée au signal de parole et une sortie;
    un bloc d'analyse de CPL ayant une entrée couplée au signal de parole et une sortie pour générer des coefficients CPL pour ledit filtre CPL; et
    un bloc de transformation ayant une entrée couplée à ladite sortie dudit bloc d'analyse CPL pour transformer lesdits coefficients CPL en une représentation des fréquences de spectre de raies (FSR), ledit bloc de transformation ayant une sortie couplée à une seconde entrée dudit filtre CPL; dans laquelle
    ledit filtre CPL comporte des polynômes auxiliaires symétrique et antisymétrique dont les zéros déterminent les fréquences de spectre de raies du filtre CPL, et dans laquelle
    ledit bloc de transformation comporte un premier moyen pour exprimer les polynômes auxiliaires en utilisant des formes explicites de polynômes de Tchebyshev;
    un second moyen pour résoudre de manière itérative un zéro d'un premier des polynômes en utilisant un zéro de l'autre des polynômes; et
    ledit bloc de transformation comporte en outre un troisième moyen pour éliminer successivement des zéros des polynômes par division de polynômes de façon à déterminer les fréquences de spectre d'immittance.
  6. Station mobile selon la revendication 5, comprenant en outre un bloc de quantification et un bloc d'interpolation couplé en série entre ladite sortie dudit bloc de transformation et ladite seconde entrée dudit filtre d'analyse CPL.
  7. Station mobile selon la revendication 5 ou 6, dans laquelle ledit troisième moyen réagit à un degré deux de l'un, donné, des polynômes, pour résoudre analytiquement le dernier zéro.
  8. Procédé pour déterminer les fréquences de spectre d'immittance d'un filtre de codeur à prédiction linéaire (CPL) exprimées sous la forme de polynômes symétrique et antisymétrique dont les zéros déterminent les fréquences de spectre d'immittance, comprenant les étapes consistant à:
    exprimer les polynômes à l'aide de formes explicites de polynômes de Tchebyshev;
    résoudre de manière itérative un zéro d'un premier des polynômes à l'aide d'un zéro de l'autre des polynômes; et
    éliminer successivement des zéros des polynômes par division de polynôme afin de déterminer les fréquences de spectre d'immittance.
  9. Procédé selon la revendication 8, comprenant en outre une étape consistant à transmettre un signal à codage CPL à une voie de communication.
  10. Procédé selon la revendication 9, dans lequel l'étape de transmission transmet le signal à codage CPL à une voie de communication radioélectrique.
  11. Procédé selon l'une quelconque des revendications 8 à 10, dans lequel le zéro de l'autre des polynômes est utilisé comme valeur initiale lors de la solution itérative du zéro du premier des polynômes.
EP19960308081 1995-11-15 1996-11-07 Détermination des fréquences du spectre de raies pour utilisation dans un radiotéléphone Expired - Lifetime EP0774750B1 (fr)

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US678795P 1995-11-15 1995-11-15
US6787 1995-11-15
US73582796A 1996-10-23 1996-10-23
US735827 1996-10-23

Publications (3)

Publication Number Publication Date
EP0774750A2 EP0774750A2 (fr) 1997-05-21
EP0774750A3 EP0774750A3 (fr) 1998-08-05
EP0774750B1 true EP0774750B1 (fr) 2003-02-05

Family

ID=26676063

Family Applications (1)

Application Number Title Priority Date Filing Date
EP19960308081 Expired - Lifetime EP0774750B1 (fr) 1995-11-15 1996-11-07 Détermination des fréquences du spectre de raies pour utilisation dans un radiotéléphone

Country Status (2)

Country Link
EP (1) EP0774750B1 (fr)
DE (1) DE69626088T2 (fr)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2783651A1 (fr) * 1998-09-22 2000-03-24 Koninkl Philips Electronics Nv Dispositif et procede de filtrage d'un signal de parole, recepteur et systeme de communications telephonique
EP2916319A1 (fr) * 2014-03-07 2015-09-09 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Concept pour le codage d'informations
EP3349212A1 (fr) * 2017-01-13 2018-07-18 Nokia Technologies Oy Procede de determination de frequences spectrales lineaires
CN110070894B (zh) * 2019-03-26 2021-08-03 天津大学 一种改进的多个病理单元音识别方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9002282D0 (en) * 1990-02-01 1990-03-28 Psion Plc Encoding speech
SE467806B (sv) * 1991-01-14 1992-09-14 Ericsson Telefon Ab L M Metod att kvantisera linjespektralfrekvenser (lsf) vid beraekning av parametrar foer ett analysfilter ingaaende i en talkodare

Also Published As

Publication number Publication date
DE69626088T2 (de) 2003-10-09
DE69626088D1 (de) 2003-03-13
EP0774750A2 (fr) 1997-05-21
EP0774750A3 (fr) 1998-08-05

Similar Documents

Publication Publication Date Title
US6122608A (en) Method for switched-predictive quantization
EP0666557B1 (fr) Interpolation de formes d'onde par décomposition en bruit et en signaux périodiques
KR100361883B1 (ko) 오디오신호압축방법,오디오신호압축장치,음성신호압축방법,음성신호압축장치,음성인식방법및음성인식장치
US5339384A (en) Code-excited linear predictive coding with low delay for speech or audio signals
US6681204B2 (en) Apparatus and method for encoding a signal as well as apparatus and method for decoding a signal
EP0501421B1 (fr) Système de codage de parole
KR101733326B1 (ko) 개선된 확률 분포 추정을 이용한 선형 예측 기반 오디오 코딩
KR100408911B1 (ko) 선스펙트럼제곱근을발생및인코딩하는방법및장치
KR20070051872A (ko) 음성 부호화 장치, 음성 복호화 장치 및 이들의 방법
KR100603167B1 (ko) 시간 동기식 파형 보간법을 이용한 피치 프로토타입파형으로부터의 음성 합성
EP0865029B1 (fr) Interpolation de formes d'onde par décomposition en bruit et en signaux périodiques
US6889185B1 (en) Quantization of linear prediction coefficients using perceptual weighting
EP0658876B1 (fr) Codeur pour des paramètres de parole
EP0819303B1 (fr) Quantification des parametres spectraux pour un codage efficace de la parole, utilisant une matrice de prediction scindee
US5233659A (en) Method of quantizing line spectral frequencies when calculating filter parameters in a speech coder
EP0774750B1 (fr) Détermination des fréquences du spectre de raies pour utilisation dans un radiotéléphone
EP0724252B1 (fr) Codeur de parole du type CELP comprenant un prédicteur à long terme amélioré
EP0899720A2 (fr) Quantisation des coefficients de prédiction linéaire
US5704001A (en) Sensitivity weighted vector quantization of line spectral pair frequencies
US4908863A (en) Multi-pulse coding system
EP0729132A2 (fr) Codeur de signaux sur canal large
EP0203940A1 (fr) Vocoder relp pour processeurs de signaux numeriques
JPH09212198A (ja) 移動電話装置における線スペクトル周波数決定方法及び移動電話装置
JP3194930B2 (ja) 音声符号化装置
Nagarajan et al. Efficient implementation of linear predictive coding algorithms

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): DE FR GB SE

PUAL Search report despatched

Free format text: ORIGINAL CODE: 0009013

AK Designated contracting states

Kind code of ref document: A3

Designated state(s): DE FR GB SE

17P Request for examination filed

Effective date: 19990205

GRAG Despatch of communication of intention to grant

Free format text: ORIGINAL CODE: EPIDOS AGRA

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: NOKIA CORPORATION

RIC1 Information provided on ipc code assigned before grant

Free format text: 7G 10L 19/06 A

17Q First examination report despatched

Effective date: 20020320

GRAG Despatch of communication of intention to grant

Free format text: ORIGINAL CODE: EPIDOS AGRA

GRAG Despatch of communication of intention to grant

Free format text: ORIGINAL CODE: EPIDOS AGRA

GRAH Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOS IGRA

GRAH Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOS IGRA

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Designated state(s): DE FR GB SE

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

REF Corresponds to:

Ref document number: 69626088

Country of ref document: DE

Date of ref document: 20030313

Kind code of ref document: P

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

ET Fr: translation filed
PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20031106

REG Reference to a national code

Ref country code: GB

Ref legal event code: 732E

REG Reference to a national code

Ref country code: FR

Ref legal event code: TP

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: SE

Payment date: 20081107

Year of fee payment: 13

EUG Se: european patent has lapsed
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20091108

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20121130

Year of fee payment: 17

Ref country code: FR

Payment date: 20121113

Year of fee payment: 17

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20121025

Year of fee payment: 17

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20131107

REG Reference to a national code

Ref country code: FR

Ref legal event code: ST

Effective date: 20140731

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20140603

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 69626088

Country of ref document: DE

Effective date: 20140603

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20131202

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20131107