CN101156318A - 预测器 - Google Patents
预测器 Download PDFInfo
- Publication number
- CN101156318A CN101156318A CNA200680011587XA CN200680011587A CN101156318A CN 101156318 A CN101156318 A CN 101156318A CN A200680011587X A CNA200680011587X A CN A200680011587XA CN 200680011587 A CN200680011587 A CN 200680011587A CN 101156318 A CN101156318 A CN 101156318A
- Authority
- CN
- China
- Prior art keywords
- value
- fallout predictor
- matrix
- calculating
- tri
- 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
Links
- 239000011159 matrix material Substances 0.000 claims description 24
- 230000017105 transposition Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000005236 sound signal Effects 0.000 description 21
- 230000003044 adaptive effect Effects 0.000 description 11
- 230000006835 compression Effects 0.000 description 6
- 238000007906 compression Methods 0.000 description 6
- 238000005303 weighing Methods 0.000 description 4
- 238000000034 method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000002203 pretreatment Methods 0.000 description 2
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
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
- 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/0017—Lossless audio signal coding; Perfect reconstruction of coded audio signal by transmission of coding error
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)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Error Detection And Correction (AREA)
Abstract
公开了一种预测器,用于计算针对多个样本值x(n)的预测值e(n),其中n是时间索引,设定P(0)=δI,δ是小的正的常数,I是M×M单位矩阵,M是所述预测器的阶,而且设定W(0)=0;以及针对每个时间索引n=1,2,…,进行上面的计算:其中,K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分并使用与所述上(或下)三角部分相同的值填充矩阵的其余部分的运算;针对每个n,确定m是否小于或等于预定值;如果m小于或等于所述预定值,则将P(n)设定为预定的矩阵。
Description
本发明涉及预测器。
无损音频编码器是这样一种音频编码器,其从原始音频信号来产生编码的音频信号以使得相应的音频解码器能够从该编码的音频信号来产生原始音频信号的准确拷贝。
在MPEG-4标准化形成过程中,产生了针对音频无损编码(ALS)的标准。无损音频编码器典型地包括两部分:线性预测器,其通过减小包含于原始音频信号中的音频样本的相关性来从原始音频信号产生残留信号;和熵编码器,其对残留信号进行编码以形成编码的音频信号。在产生残留信号的过程中,预测器将相关性减小得越多,可实现的对原始音频信号的压缩就越多,即所编码的音频信号相对于原始音频信号的压缩比越高。
如果原始音频信号是立体声的信号,即包括针对第一通道和第二通道的音频样本,那么就既有通道内的相关性,即相同通道的音频样本之间的相关性,也有通道间的相关性,即不同通道的音频样本之间的相关性。
通常用于无损音频编码的线性预测器是基于RLS(递归最小平方)算法的预测器。
经典的RLS算法可以总结如下:
通过设定P(0)=δI来初始化该算法,
其中δ是小的正的常数,I是M×M单位矩阵,其中M是预测器的阶。
此外,被定义为W(n)=[w0(n),w1(n),…wM-1(n)]T的M×1加权向量W(n)通过W(0)=0来初始化。
对于每个时刻,n=1,2,...,进行下面的计算:
V(n)=P(n-1)*X(n)
其中X(n)是M×1矩阵(即M维向量)形式的输入信号,其被定义为X(n)=[x(n-1),...,x(n-M)]T
(P(n)是M×M矩阵,因此V(n)是M×1矩阵)
K(n)=m*V(n)
e(n)=x(n)-WT(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*VT(n)]}
K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分并使用与上(或下)三角部分相同的值填充矩阵的其余部分的运算。
使用定点数学来实现上述的经典RLS算法存在两个问题。
第一,由于固定点的有限动态范围,变量m往往容易被舍入为零。如果m是零,则K(n)将为零,P(n)将根据λ-1(略大于1)缓慢增加,并且除非以减小XT V(n)(XT V(n)的高值导致m为零)的方式改变该输入X(n),否则P(n)将最终溢出。
第二,V(n)的动态范围非常大(有时大于232),同时需要高精度(至少32比特)以保持高预测增益。然而,因为对于大多数32比特定点实现,用于上述方程的变量的动态范围太大。所以,当使用类似用于该算法中的其它变量的定点实现来编码V(n)时,存在精确度损失。
本发明的一个目的是为了解决当使用固定点实现的RLS算法时出现的发散问题和精确度问题。
该目的通过具有根据独立权利要求的特征的预测器来实现。
提供了一种预测器,用于计算针对多个样本值x(n)的预测值e(n),其中n是时间索引,其中
设定P(0)=δI,其中δ是小的正的常数,I是M×M单位矩阵,其中M是预测器的阶,并且设定W(0)=0;
针对每个时间索引n=1,2,...,进行下面的计算:
V(n)=p(n-1)*X(n)where X(n)=[x(n-1),...,x(n-M)]T
K(n)=m*V(n)
e(n)=x(n)-WT(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*VT(n)]}
其中K(n)是M×1矩阵(即M维向量),λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分,并使用与上(或下)三角部分相同的值填充矩阵的其余部分的运算;
并且其中进一步针对每个n,确定m是否小于或等于预定值,如果m小于或等于预定值,则将P(n)设定为预定的矩阵。
进一步提供一种预测器,用于计算针对多个样本值x(n)的预测值e(n),其中n是时间索引,其中
设定P(0)=δI,其中δ是小的正的常数,I是M×M单位矩阵,其中M是预测器的阶,并且设定W(0)=0;
而且针对每个时间索引n=1,2,…,进行下面的计算:
V(n)=P(n-1)*X(n)whereX(n)=[x(n-1),...,x(n-M)]T
K(n)=m*V(n)
e(n)=x(n)-WT(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*VT(n)]}
其中K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分并使用与上(或下)三角部分相同的值填充矩阵的其余部分的操作,而且其中更进一步,变量V(n)被编码为标量乘以变量V′(n)的积,该标量以使V′(n)位于预定区间内的方式被预先确定。
作为解释,当值m已经在预测算法的一个步骤中变得非常小时,P(n)被重新初始化。这样,***因为P(n)将不溢出而保持稳定。
此外,使用一标量,即缩放因子(下面用vscale表示)对V(N)进行缩放,以使得V(N)=vscale*V′(N)。这样,与V(N)相比,缩放的变量V′(N)的范围被缩减。因此,当定点实现被用于编码V′(N)时,不会有精确度损失。
例如,根据MPEG-4标准规范,可以使用小的常数0.0001初始化P(0)。在另一实施例中,设置P(0)=δ-1 I,其中δ是小的正的常数。
本发明的优选实施例由从属权利要求来体现。
在一个实施例中,预定值为0。预定值也可以是小的正的常数。预定的向量是例如P(0)=δI。预定向量也可以是P(0)=δ-1 I。在一个实施例中,将定点实现用于计算。特别地,在一个实施例中,使用定点实现来编码V′(n)。
下面参考附图来解释本发明的说明性的实施例。
图1示出根据本发明实施例的编码器。
图2示出根据本发明实施例的解码器。
图1示出根据本发明实施例的编码器100。
编码器100接收原始音频信号101作为输入。
该原始音频信号包括多个帧。将每个帧分成块。每个块包括多个样本。该音频信号可以包括用于多个音频通道的音频信息。在这种情况下,通常帧包括用于每个通道的一个块,即帧中的每个块对应一个通道。
原始音频信号101是数字音频信号,并且例如是通过以某采样率(例如48kHz、96kHz和192kHz)来采样模拟音频信号而产生的,其中每个样本具有某分辨率(例如8比特、16比特、10比特和14比特)。
提供缓冲器102来存储一个帧,即包含在一个帧中的音频信息。
通过自适应预测器103来处理原始音频信号101(即处理原始音频信号101的所有样本),所述自适应预测器103根据原始音频信号101的过去样本的过去样本值来计算原始音频信号101的当前(即当前处理的)样本的当前样本值的预测(估计)104。为此,自适应的预测器103使用自适应算法。这个过程将在下面详细描述。
针对当前样本值的预测104通过减法单元106从当前样本值减去,以产生当前残留105。
当前残留105然后通过熵编码器107被熵编码。熵编码器107能够例如执行Rice编码或BGMC(块式Gilbert-Moore代码)编码。
所编码的当前残留、描述由熵编码器107执行的当前残留105的编码的编码索引、由自适应预测器用来产生预测104的预测器系数、以及任选地其它信息通过多路复用器108被多路复用,以使得当处理原始信号101的所有样本时,形成比特流109,所述比特流109持有无损编码的原始信号101和把它解码的信息。
对于编码和压缩原始音频信号,编码器100可以提供几个具有不同复杂度的压缩等级。然而,对于高压缩等级,编码效率方面的差异通常是相当小的,因此放弃最高的压缩以便减小计算工作量可能是合适的。
通常,比特流109以某种方式例如经由计算机网络传送给在下面所描述的解码器。
图2示出根据本发明实施例的解码器200。
解码器200接收对应于比特流109的比特流201作为输入。
作为解释,解码器200执行编码器的逆功能。
正如所说明的,比特流201持有编码的残留、代码索引和预测器系数。该信息通过多路分离器202从比特流201中被多路分离。
使用各自的编码索引,当前(即:当前处理的)编码的残留通过熵编码器203被解码以形成当前残留206。
因为假设已经处理了对应于当前残留206的样本之前的样本的样本值,所以类似于自适应预测器103的自适应预测器204能够产生当前样本值的预测205,即要根据当前残留206来无损重构的样本值,该预测205通过加法器单元207与当前残留206相加。
加法单元207的输出是无损重构的当前样本,该无损重构的当前样本与由编码器100来处理以形成当前编码的残留的样本是相同的。
解码器200的计算量依赖于自适应预测器204的阶,自适应预测器204的阶通过编码器100来选择。除了自适应预测器204的阶以外,解码器200的复杂度与编码器100的复杂度相同。
编码器100在一个实施例中也提供CRC(循环冗余校验)校验和,该校验和在比特流109中被提供给解码器200,以便解码器200能够校验被解码的数据。在编码器100侧,CRC校验和能够被用于确保压缩的文件是可无损解码的。
下面说明根据本发明的一个实施例的自适应预测器103的功能性。
通过设定P(0)=δI来初始化预测器,
其中δ是小的正的常数,I是M×M单位矩阵,其中M是预测器阶。
此外,通过设定W(0)=0来对定义为:
W(n)=[w0(n),w1(n),…wM-1(n)]T的M×1加权向量W(n)进行初始化,该M×1加权向量W(n)示例性地作为初始滤波器加权向量。
针对每个时刻,即针对要被预测器处理的每个样本值x(n),其中n=1,2,...,是相应的时间索引,进行下面的计算:
V(n)=P(n-1)*X(n)
其中X(n)是M×1矩阵形式的输入信号,其被定义为
X(n)=[x(n-1),...,x(n-M)]T
(P(n)是M×M矩阵,因此V(n)是M×1矩阵)
向量X(n)是当前样本值x(n)之前的样本值的向量。作为解释,向量X(n)持有用于预测当前值的过去值。
K(n)=m*V(n)
e(n)=x(n)-W T(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*VT(n)]}
K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号(表示转置运算),Tri表示计算P(n)的上(或下)三角部分并使用与上(或下)三角部分相同的值填充矩阵的其余部分的运算。
为了防止出现由于m被舍入为零而导致发散的问题,在每个步骤中检验m是否为零。如果是这种情况,则例如根据P(n)=δI来重新初始化P(n)。
为了解决由变量、特别是V(n)被以定点格式编码的所导致的精确度损失的问题,引入缩放因子vscale。
严格地选择缩放因子vscale以用于V(n)。利用与移位参数有关vscale,比例因子vscale使其它变量能够以32比特的形式来简单表示。这样,该算法能够主要用32比特定点运算而不是仿效浮点数学运算来操作。
V(n)被编码为vscale和变量V′(n)的积。选择vscale以使得V′(n)可以按照定点格式编码而不损失(或很少损失)精确度,例如与浮点实现相比。
Claims (6)
1.一种预测器,用于计算针对多个样本值x(n)的预测值e(n),其中n是时间索引,其中
设定P(0)=δI,其中δ是小的正的常数,I是M×M单位矩阵,其中M是所述预测器的阶,而且设定W(0)=0;
以及针对每个时间索引n=1,2,...,进行下面的计算:
V(n)=P(n-1)*X(n)where X(n)=[x(n-1),...,x(n-M)]T
K(n)=m*V(n)
e(n)=x(n)-W T(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*V T(n)]}
其中,K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分并使用与所述上(或下)三角部分相同的值填充矩阵的其余部分的运算;
以及其中进一步
针对每个n,确定m是否小于或等于预定值;
如果m小于或等于所述预定值,则将P(n)设定为预定的矩阵。
2.根据权利要求1所述的预测器,其中所述预定值是小的正的常数。
3.根据权利要求1或2所述的预测器,其中所述预定的向量是δI。
4.根据权利要求1至3中的任何一项所述的预测器,其中采用定点实现进行所述计算。
5.一种预测器,其用于计算针对多个样本值x(n)的预测值e(n),其中n是时间索引,其中
设定P(0)=δI,其中δ是小的正的常数,I是M×M单位矩阵,其中M是所述预测器的阶,而且设定W(0)=0;
以及针对每个时间索引n=1,2,...,进行下面的计算:
V(n)=P(n-1)*X(n)whereX(n)=[x(n-1),...,x(n-M)]T
K(n)=m*V(n)
e(n)=x(n)-W T(n-1)X(n)
W(n)=W(n-1)+K(n)e(n)
P(n)=Tri{λ-1[P(n-1)-K(n)*V T(n)]}
K(n)是M×1矩阵,λ是略小于1的正值,T是转置符号,Tri表示计算P(n)的上(或下)三角部分,并使用与所述上(或下)三角部分相同的值填充矩阵的其余部分的运算;
以及其中进一步
所述变量V(n)被编码为一标量乘以变量V′(n)的乘积
所述标量以使得V′(n)位于预定区间内的方式来预先确定。
6.根据权利要求5所述的预测器,其中使用定点实现来编码所述变量V′(n)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US66066905P | 2005-03-11 | 2005-03-11 | |
US60/660,669 | 2005-03-11 | ||
PCT/SG2006/000049 WO2006096137A2 (en) | 2005-03-11 | 2006-03-09 | Predictor |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101156318A true CN101156318A (zh) | 2008-04-02 |
CN101156318B CN101156318B (zh) | 2012-05-09 |
Family
ID=36953767
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200680011587XA Expired - Fee Related CN101156318B (zh) | 2005-03-11 | 2006-03-09 | 编码器 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20100023575A1 (zh) |
EP (1) | EP1859531A4 (zh) |
CN (1) | CN101156318B (zh) |
SG (1) | SG160390A1 (zh) |
TW (1) | TW200703940A (zh) |
WO (1) | WO2006096137A2 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021246A (zh) * | 2014-05-28 | 2014-09-03 | 复旦大学 | 一种应用于低功耗容错电路的自适应长度预测器 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101167124B (zh) * | 2005-04-28 | 2011-09-21 | 松下电器产业株式会社 | 语音编码装置和语音编码方法 |
RU2007139784A (ru) * | 2005-04-28 | 2009-05-10 | Мацусита Электрик Индастриал Ко., Лтд. (Jp) | Устройство кодирования звука и способ кодирования звука |
JP6148811B2 (ja) * | 2013-01-29 | 2017-06-14 | フラウンホーファーゲゼルシャフト ツール フォルデルング デル アンゲヴァンテン フォルシユング エー.フアー. | 周波数領域におけるlpc系符号化のための低周波数エンファシス |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5568378A (en) * | 1994-10-24 | 1996-10-22 | Fisher-Rosemount Systems, Inc. | Variable horizon predictor for controlling dead time dominant processes, multivariable interactive processes, and processes with time variant dynamics |
US5664053A (en) * | 1995-04-03 | 1997-09-02 | Universite De Sherbrooke | Predictive split-matrix quantization of spectral parameters for efficient coding of speech |
US5923711A (en) * | 1996-04-02 | 1999-07-13 | Zenith Electronics Corporation | Slice predictor for a signal receiver |
US6463410B1 (en) * | 1998-10-13 | 2002-10-08 | Victor Company Of Japan, Ltd. | Audio signal processing apparatus |
JP3387089B2 (ja) * | 2000-10-20 | 2003-03-17 | 日本ビクター株式会社 | 音声符号化装置 |
-
2006
- 2006-03-09 SG SG201001709-3A patent/SG160390A1/en unknown
- 2006-03-09 EP EP06717173A patent/EP1859531A4/en not_active Withdrawn
- 2006-03-09 CN CN200680011587XA patent/CN101156318B/zh not_active Expired - Fee Related
- 2006-03-09 WO PCT/SG2006/000049 patent/WO2006096137A2/en active Application Filing
- 2006-03-09 US US11/908,300 patent/US20100023575A1/en not_active Abandoned
- 2006-03-10 TW TW095108102A patent/TW200703940A/zh unknown
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104021246A (zh) * | 2014-05-28 | 2014-09-03 | 复旦大学 | 一种应用于低功耗容错电路的自适应长度预测器 |
Also Published As
Publication number | Publication date |
---|---|
SG160390A1 (en) | 2010-04-29 |
WO2006096137A2 (en) | 2006-09-14 |
TW200703940A (en) | 2007-01-16 |
EP1859531A2 (en) | 2007-11-28 |
US20100023575A1 (en) | 2010-01-28 |
EP1859531A4 (en) | 2008-04-09 |
CN101156318B (zh) | 2012-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
USRE49363E1 (en) | Variable bit rate LPC filter quantizing and inverse quantizing device and method | |
RU2541864C2 (ru) | Аудио или видео кодер, аудио или видео и относящиеся к ним способы для обработки многоканальных аудио или видеосигналов с использованием переменного направления предсказания | |
KR100469002B1 (ko) | 오디오 코딩 방법 및 장치 | |
KR100936498B1 (ko) | 스테레오 호환성의 멀티채널 오디오 코딩 | |
EP2543038B1 (en) | Decoding of multi-channel audio signals using complex prediction | |
RU2005137320A (ru) | Способ и устройство для квантования усиления в широкополосном речевом кодировании с переменной битовой скоростью передачи | |
JP3254687B2 (ja) | 音声符号化方式 | |
US8665945B2 (en) | Encoding method, decoding method, encoding device, decoding device, program, and recording medium | |
CN101484937A (zh) | 使用缓冲器调节对已预测编码的数据进行解码 | |
JP2009513993A (ja) | 変換表現に転換するか、または変換表現を逆変換するための装置および方法 | |
EP1847022B1 (en) | Encoder, decoder, method for encoding/decoding, computer readable media and computer program elements | |
EP1837997B1 (en) | Long-term prediction encoding method, long-term prediction decoding method, devices thereof, program thereof, and recording medium | |
JP3557255B2 (ja) | Lspパラメータ復号化装置及び復号化方法 | |
JPH08179795A (ja) | 音声のピッチラグ符号化方法および装置 | |
CA2147394C (en) | Quantization of input vectors with and without rearrangement of vector elements of a candidate vector | |
CN101156318A (zh) | 预测器 | |
US20050075888A1 (en) | Fast codebook selection method in audio encoding | |
JP6780108B2 (ja) | 符号化装置、復号装置、平滑化装置、逆平滑化装置、それらの方法、およびプログラム | |
WO2023198383A1 (en) | Method for quantizing line spectral frequencies | |
JP3557414B2 (ja) | Lspパラメータ符号化装置及び符号化方法 | |
KR20230023560A (ko) | 부호화 방법 및 복호화 방법, 상기 방법을 수행하는 부호화기 및 복호화기 | |
AU2021436955A1 (en) | Encoding device, decoding device, encoding method, and decoding method | |
JP2001298367A (ja) | オーディオ信号符号化方法、オーディオ信号復号化方法およびオーディオ信号符号化/復号化装置と前記方法を実施するプログラムを記録した記録媒体 | |
JPH04107600A (ja) | コード励振線形予測符号化器及び復号化器 | |
JPH05341800A (ja) | 音声符号化装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120509 Termination date: 20130309 |