CN101621306B - 多输入多输出系统预编码矩阵的映射方法和装置 - Google Patents

多输入多输出系统预编码矩阵的映射方法和装置 Download PDF

Info

Publication number
CN101621306B
CN101621306B CN 200810129137 CN200810129137A CN101621306B CN 101621306 B CN101621306 B CN 101621306B CN 200810129137 CN200810129137 CN 200810129137 CN 200810129137 A CN200810129137 A CN 200810129137A CN 101621306 B CN101621306 B CN 101621306B
Authority
CN
China
Prior art keywords
matrix
codebook
coding
receiver
coding matrix
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 - Fee Related
Application number
CN 200810129137
Other languages
English (en)
Other versions
CN101621306A (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.)
ZTE Corp
Original Assignee
ZTE Corp
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 ZTE Corp filed Critical ZTE Corp
Priority to CN 200810129137 priority Critical patent/CN101621306B/zh
Publication of CN101621306A publication Critical patent/CN101621306A/zh
Application granted granted Critical
Publication of CN101621306B publication Critical patent/CN101621306B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开了一种多输入多输出系统预编码矩阵的映射方法和装置,该方法包括:根据接收机的信道矩阵进行奇异值分解得到预编码矩阵;确定预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵;确定各个相关矩阵绝对值的迹;将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。所述装置包括设置在接收机的奇异值分解模块、相关矩阵确定模块、迹确定模块以及映射反馈模块,用于执行上述方法。利用本发明,可以降低多输入多输出系统接收机的运算处理的复杂度,提高处理速度,增强实时性。

Description

多输入多输出系统预编码矩阵的映射方法和装置
技术领域
本发明涉及无线通信技术领域,尤其涉及一种多输入多输出(MIMO,Multiple Input Multiple Output)系统预编码矩阵映射方法和装置。 
背景技术
多输入多输出系统的发射机和接收机分别包含多个通道(通道包括基带处理、中频处理、以及射频天线等单元),其技术核心在于将调制后的用户数据分解为多个并行数据流,发射机同时发射多个并行数据流,接收机接收各个并行数据流后,基于发射机与接收机之间各通道之间的空间特性(SS,Spatial Signature),对接收数据进行检测和解调,从而得到对应的原数据流。若发射机或接收机仅包含单个通道,则对应的系统称为单输入多输出(SIMO,Single Input Multiple Output)系统或多输入单输出(MISO,MultipleInput Single Output)系统。SIMO系统和MISO系统可认为是MIMO系统的特殊形式。利用MIMO技术,不仅能够显著提高信道容量,而且可以明显增强系统可靠性,降低误码率。 
预编码(Precoding)是MIMO的重要实现技术之一,其目的在于尽可能保证各个并行数据流之间的正交性,从而降低接收机检测和解调的复杂度。在实现上,接收机首先对信道矩阵进行奇异值分解(SVD,Singular ValueDecomposition),得到预编码矩阵,然后利用码本实现预编码矩阵的映射,获得对应的预编码矩阵索引(PMI,Precoding Matrix Index),同时将PMI索引通过无线链路反馈至发射机,发射机基于PMI索引得到对应的预编码矩阵,并对发射数据进行预编码处理。预编码矩阵码本用于对预编码矩阵进行量化,在接收机和发射机均有存储,因此接收机只需反馈对应的PMI索 引给发射机即可,从而明显降低通信系统空口传输的数据反馈量。 
显然,在预编码矩阵码本一定的情况下,预编码矩阵的映射方法将直接影响通信系统的性能。中国专利申请公开号CN101136718、公开日2008年3月5日、发明名称为《无线通信系统中多输入多输出的空间复用的预编码方法》的专利申请的公开文本中公开了一种多输入多输出的预编码矩阵的映射方法。在该现有技术中,采用有限反馈的预编码的空间复用的方法,即预先定义一组发射机和接收机都已知的预编码矩阵的集合(也就是预编码码本);接收机通过信道估计算法获得当前的信道状态信息,按照均方误差矩阵的迹最小化的原则从预编码的码本中选择一个最优的预编码矩阵,也就是说映射到一个最优的预编码矩阵,得到该与编码矩阵的PMI,然后通过反馈信道在信道的相干时间内将此映射得到的预编码矩阵的PMI反馈到发射机,发射机根据此预编码矩阵的PMI来获得预编码矩阵,然后根据此预编码矩阵对发射信号向量进行预编码操作。 
但是,上述现有技术的预编码矩阵映射方法仅适用于理论分析预编码矩阵的码本映射,在接收机获取选择最优预编码矩阵的PMI的过程中,需要按照均方误差矩阵的迹最小化的原则进行运算,并在该过程中涉及矩阵求逆等运算,这些运算处理的复杂性较高,工程实现难度较高,导致处理速度低,实时性较差。 
发明内容
有鉴于此,本发明所要解决的技术问题在于提供一种多输入多输出系统预编码矩阵的映射方法,降低处理的复杂度,增强实时性。 
本发明所要解决的另一技术问题在于提供一种多输入多输出系统预编码矩阵的映射装置,以降低预编码映射处理的复杂度,增强实时性。 
为了实现上述发明目的,本发明的主要技术方案为: 
一种多输入多输出系统预编码矩阵的映射方法,包括: 
A、根据接收机的信道矩阵进行奇异值分解得到预编码矩阵; 
B、确定步骤A所述预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵; 
C、确定各个相关矩阵绝对值的迹; 
D、将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。 
优选的,步骤A具体包括: 
A1、根据接收机接收的数据进行信道估计,得到信道矩阵; 
A2、根据式H=U·D·VH对所述信道矩阵进行奇异值分解,得到预编码矩阵;其中:H为所述信道矩阵;U为译码矩阵,为MR×Λ阶酉矩阵;D为Λ阶奇异值对角矩阵;V为所要求解的预编码矩阵,是MT×Λ阶酉矩阵;MR为接收的通道端口数,MT为发射的通道端口数;Λ为所述信道矩阵的秩或传输数据流数。 
优选的,步骤B中,所述预编码矩阵与所述各个码本矩阵之间的相关矩阵通过式R(k)=VH·W(k)或通过式R(k)=WH(k)·V得到;其中:V为步骤A所得到的预编码矩阵;R(k)为所述预编码矩阵V与第k个码本矩阵的相关矩阵;W(k)为第k个码本矩阵,k=1,...,K,K为码本数。 
优选的,步骤C中,所述相关矩阵绝对值的迹为该相关矩阵主对角线元素的绝对值之和。 
优选的,所述预知的预编码矩阵的各个码本矩阵预先存储在接收机和发射机中。 
优选的,该方法进一步包括:所述发射机接收到所述码本矩阵的索引后,从所述预知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。 
一种多输入多输出系统预编码矩阵的映射装置,包括设置在接收机的奇异值分解模块、相关矩阵确定模块、迹确定模块以及映射反馈模块,其中: 
奇异值分解模块用于对接收机的信道矩阵进行奇异值分解,得到预编码 矩阵; 
相关矩阵确定模块用于确定奇异值分解模块得到的预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵; 
迹确定模块用于确定各个相关矩阵绝对值的迹; 
映射反馈模块用于将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。 
优选的,所述相关矩阵绝对值的迹为该相关矩阵主对角线元素的绝对值之和。 
优选的,所述预知的预编码矩阵的各个码本矩阵预先存储在接收机和发射机中。 
优选的,该装置进一步包括预编码模块,设置在发射机,用于在发射机接收到所述码本矩阵的索引后,从所述预知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。 
相对于现有技术,由于本发明通过计算预编码矩阵与各个码本矩阵之间的相关矩阵,再计算各个相关矩阵绝对值的迹,最后基于迹值最大对应的预编码矩阵码本即为所求,并据此反馈PMI;在上述计算处理的过程中,所述相关性运算主要是计算内积,所述迹的运算主要是计算主对角线元素之和,这些运算处理在工程实现上比较简单,不需要涉及到矩阵求逆等复杂运算,因此本发明可以降低处理过程的复杂度,既易于工程实现,又加快了处理速度,增强了预编码的实时性。 
附图说明
图1是本发明所述MIMO系统预编码矩阵映射方法的一种实施例的流程图; 
图2是本发明所述MIMO系统预编码矩阵映射装置的一种实施例的组成结构示意图。 
具体实施方式
下面通过具体实施例和附图对本发明做进一步详细说明。 
本发明的核心思想是:对接收机的信道矩阵进行奇异值分解,得到预编码矩阵;然后确定预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵;接着确定所述各个相关矩阵绝对值的迹;最后,将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。 
图1是本发明所述MIMO系统预编码矩阵映射方法的一种实施例流程图。如图1所示,该实施例包括以下步骤: 
步骤10,接收机对信道矩阵进行奇异值分解(SVD),得到预编码矩阵。 
本步骤中,首先接收机根据接收的数据进行信道估计,得到信道矩阵。 
接收机的接收数据表示为: 
y=H·x+n    (1) 
式(1)中x为发射机的发射数据;n为加性高斯白噪声;H为信道矩阵,表示为: 
H = h 1,1 h 1,2 · · · h 1 , M T h 2,1 h 2,2 · · · h 2 , M T · · · · · · h m R , m T · · · h M R , 1 h M R , 2 · · · h M R , M T - - - ( 2 )
式(2)中 
Figure S2008101291371D00052
为第mR接收通道端口、第mT发射通道端口对应的信道估计;mR=1,2,......,MR;mT=1,2,......,MR和MT分别表示接收与发射的通道端口数。所述 可以由接收机根据目前常用的信道估计方法得到。 
其次,对信道矩阵H进行奇异值分解(SVD),即根据式(3)进行处理得到预编码矩阵: 
H=U·D·VH    (3) 
式(3)中,U为MR×Λ阶酉矩阵,称为译码矩阵,Λ为所述信道矩阵H的秩或传输数据流数;D为Λ阶奇异值对角矩阵;VH表示对V进行共轭转置运算,上标H表示共轭转置;V为MT×Λ阶酉矩阵,即所要求解的预编码矩阵,表示为V={vi,j|i=1,...,MT;j=1,...,Λ},即 
V = v 1,1 v 1,2 · · · v 1 , Λ v 2,1 v 2,2 · · · v 2 , Λ · · · · · · · · · · · · v M T , 1 v M T , 2 · · · v M T , Λ - - - ( 4 )
步骤20,计算所述预编码矩阵V与预编码矩阵的各个码本矩阵之间的相关矩阵。 
此处,预编码矩阵的码本矩阵为预先定义的一组发射机和接收机都已知的预编码矩阵的集合,预先存储在发射机和接收机中,预编码矩阵的各个码本矩阵表示为W(k)={wi,j(k)|i=1,...,MT;j=1,..,Λ},即: 
W ( k ) = w 1,1 ( k ) w 1,2 ( k ) · · · w 1 , Λ ( k ) w 2,1 ( k ) w 2,2 ( k ) · · · w 2 , Λ ( k ) · · · · · · · · · · · · w M T , 1 ( k ) w M T , 2 ( k ) · · · w M T , Λ ( k ) - - - ( 5 )
式中k=1,...,K为码本索引,K为码本数。 
所述预编码矩阵V与第k个码本矩阵之间的相关矩阵根据式(6)或式(7)得到: 
R(k)=VH·W(k)    (6) 
R(k)=WH(k)·V    (7) 
步骤30,计算所述各个相关矩阵绝对值的迹。 
预编码矩阵V与第k个码本矩阵的相关矩阵为R(k)={ri,j(k)|i,j=1,...,Λ},则R(k)绝对值的迹即为主对角线元素的绝对值之和(即模值之和),也就是R(k)对应的绝对值矩阵 R · ( k ) = { | r i , j ( k ) | | i , j = 1 , . . . , Λ } 的迹。 
所述各个相关矩阵绝对值的迹β(k)表示为 
β ( k ) = Σ l = 1 Λ | r l , l ( k ) | = trace { R · ( k ) } - - - ( 8 )
步骤40,基于最大迹值准则,将所述各个相关矩阵绝对值的迹的最大值对应的码本矩阵作为所求的预编码矩阵码本矩阵,该码本矩阵对应的索引即为反馈至发射机的预编码矩阵索引(PMI)。 
基于迹值最大准则,所求的预编码矩阵索引以及码本矩阵WPrecoding表示为: 
PMI = I PMI = arg k { max [ β ( k ) ] } - - - ( 9 )
WPrecoding=W(IPMI)    (10) 
通过上述映射方法,可以从已知的预编码矩阵码本中映射选择得到一个最优的矩阵码本,将该码本矩阵的PMI通过反馈信道在信道的相干时间内反馈到发射机。发射机根据接收到的所述PMI从已知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。 
图2为本发明所述MIMO系统预编码矩阵映射装置的一种实施例的组成结构示意图。参见图2,该装置包括设置在接收机200的奇异值分解模块201、相关矩阵确定模块202、迹确定模块203以及映射反馈模块204,其中: 
奇异值分解模块201用于对接收机的信道矩阵进行奇异值分解,得到预编码矩阵. 
相关矩阵确定模块202用于确定奇异值分解模块得到的预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵;所述预知的预编码矩阵的各个码本矩阵预先存储在接收机200和发射机100中。 
迹确定模块203用于确定所述各个相关矩阵绝对值的迹。 
映射反馈模块204用于将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。 
所述相关矩阵绝对值的迹为该相关矩阵主对角线元素的绝对值之和。并且,上述各个模块的具体计算处理过程请参见上述的MIMO系统预编码矩 阵映射方法。 
另外,本发明的所述装置还可以进一步包括预编码模块205,设置在发射机100,用于在发射机100接收到所述码本矩阵的索引后,从所述预知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。 
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。 

Claims (10)

1.一种多输入多输出系统预编码矩阵的映射方法,其特征在于,包括:
A、根据接收机的信道矩阵进行奇异值分解得到预编码矩阵;
B、确定步骤A所述预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵;
C、确定各个相关矩阵绝对值的迹;
D、将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。
2.根据权利要求1所述的方法,其特征在于,步骤A具体包括:
A1、根据接收机接收的数据进行信道估计,得到信道矩阵;
A2、根据式H=U·D·VH对所述信道矩阵进行奇异值分解,得到预编码矩阵;其中:H为所述信道矩阵;U为译码矩阵,为MR×Λ阶酉矩阵;D为Λ阶奇异值对角矩阵;V为所要求解的预编码矩阵,是MT×Λ阶酉矩阵;MR为接收的通道端口数,MT为发射的通道端口数;Λ为所述信道矩阵的秩或传输数据流数。
3.根据权利要求1所述的方法,其特征在于,步骤B中,所述预编码矩阵与所述各个码本矩阵之间的相关矩阵通过式R(k)=VH·W(k)或通过式R(k)=WH(k)·V得到;其中:V为步骤A所得到的预编码矩阵;R(k)为所述预编码矩阵V与第k个码本矩阵的相关矩阵;W(k)为第k个码本矩阵,k=1,...,K,K为码本数。
4.根据权利要求1所述的方法,其特征在于,步骤C中,所述相关矩阵绝对值的迹为该相关矩阵主对角线元素的绝对值之和。
5.根据权利要求1所述的方法,其特征在于,所述预知的预编码矩阵的各个码本矩阵预先存储在接收机和发射机中。
6.根据权利要求5所述的方法,其特征在于,该方法进一步包括:所述发射机接收到所述码本矩阵的索引后,从所述预知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。
7.一种多输入多输出系统预编码矩阵的映射装置,其特征在于,包括设置在接收机的奇异值分解模块、相关矩阵确定模块、迹确定模块以及映射反馈模块,其中:
奇异值分解模块用于对接收机的信道矩阵进行奇异值分解,得到预编码矩阵;
相关矩阵确定模块用于确定奇异值分解模块得到的预编码矩阵与预知的预编码矩阵的各个码本矩阵之间的相关矩阵;
迹确定模块用于确定各个相关矩阵绝对值的迹;
映射反馈模块用于将所述各个相关矩阵绝对值的迹中最大值对应的码本矩阵作为所映射的预编码矩阵,将该码本矩阵的索引反馈给发射机。
8.根据权利要求7所述的映射装置,其特征在于,所述相关矩阵绝对值的迹为该相关矩阵主对角线元素的绝对值之和。
9.根据权利要求7所述的映射装置,其特征在于,所述预知的预编码矩阵的各个码本矩阵预先存储在接收机和发射机中。
10.根据权利要求9所述的映射装置,其特征在于,该装置进一步包括预编码模块,设置在发射机,用于在发射机接收到所述码本矩阵的索引后,从所述预知的预编码矩阵的码本矩阵中得到对应的码本矩阵,然后根据该码本矩阵对经过空时编码的发射信号向量进行预编码操作。
CN 200810129137 2008-06-30 2008-06-30 多输入多输出系统预编码矩阵的映射方法和装置 Expired - Fee Related CN101621306B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200810129137 CN101621306B (zh) 2008-06-30 2008-06-30 多输入多输出系统预编码矩阵的映射方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200810129137 CN101621306B (zh) 2008-06-30 2008-06-30 多输入多输出系统预编码矩阵的映射方法和装置

Publications (2)

Publication Number Publication Date
CN101621306A CN101621306A (zh) 2010-01-06
CN101621306B true CN101621306B (zh) 2013-01-09

Family

ID=41514390

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200810129137 Expired - Fee Related CN101621306B (zh) 2008-06-30 2008-06-30 多输入多输出系统预编码矩阵的映射方法和装置

Country Status (1)

Country Link
CN (1) CN101621306B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102771073A (zh) * 2010-04-29 2012-11-07 富士通株式会社 生成码书的方法和装置及使用该码书的通信系统、移动台、基站和通信方法
US10044423B2 (en) * 2016-10-20 2018-08-07 Samsung Electronics Co., Ltd System and method for precoder selection in multiple-input multiple-output (MIMO) systems with discrete fourier transform (DFT)-based codebook
CN106503696B (zh) * 2016-12-05 2019-08-13 电子科技大学 一种针对视觉映射目标值的增强编码方法
CN114065122A (zh) * 2020-07-31 2022-02-18 深圳市中兴微电子技术有限公司 数据处理方法、设备和存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006049417A1 (en) * 2004-11-01 2006-05-11 Lg Electronics Inc. A method of transmitting a precoding matrix in a multi-input multi-output (mimo) system
CN1956430A (zh) * 2005-10-27 2007-05-02 财团法人工业技术研究院 Mimo-ofdm系统及其中之预编码与反馈的方法
CN101047417A (zh) * 2007-04-20 2007-10-03 哈尔滨工程大学 一种多用户mimo系统下行链路天线选择预处理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006049417A1 (en) * 2004-11-01 2006-05-11 Lg Electronics Inc. A method of transmitting a precoding matrix in a multi-input multi-output (mimo) system
CN1956430A (zh) * 2005-10-27 2007-05-02 财团法人工业技术研究院 Mimo-ofdm系统及其中之预编码与反馈的方法
CN101047417A (zh) * 2007-04-20 2007-10-03 哈尔滨工程大学 一种多用户mimo系统下行链路天线选择预处理方法

Also Published As

Publication number Publication date
CN101621306A (zh) 2010-01-06

Similar Documents

Publication Publication Date Title
CN102725967B (zh) 用于信息反馈以及预编码的方法和装置
US8644411B2 (en) Generalized reference signaling scheme for multi-user multiple input, multiple output (MU-MIMO) using arbitrarily precoded reference signals
US7912141B2 (en) Pre-coding method for MIMO system and apparatus using the method
CN101848071B (zh) 分层空时系统中有限反馈预编码非线性译码方法
CN101232317A (zh) 用于德耳塔量化器的处理通信信号的方法和系统
CN101232356A (zh) Mimo系统中的预编码方法、系统和装置
CN102088340A (zh) 一种多天线系统发送、接收信息的方法及装置
CN101207464B (zh) 广义格拉斯曼码本的反馈方法
CN101483503B (zh) 多输入多输出系统数据发送、接收方法及装置
CN101621306B (zh) 多输入多输出系统预编码矩阵的映射方法和装置
CN101094022A (zh) 发射机、通信系统以及通信方法
CN101262311B (zh) 多用户mimo系统下行链路发射端信息数据的预处理方法
CN101867399A (zh) 一种多输入多输出系统中的多用户调度方法及装置
CN106301496A (zh) 基于天线选择和预编码的空间调制系统
CN101667860B (zh) 在多天线数字无线通信系统中检测信号的方法
CN101964696A (zh) 一种应用于空分复用的预编码矩阵选择方法及系统
CN101645757B (zh) 一种预编码矩阵选择方法和装置
CN103973409A (zh) 预编码装置和包含该预编码装置的发射机
CN103368698B (zh) 一种多输入多输出mimo预编码方法、设备及系统
Wu et al. Deep joint source-channel coding for adaptive image transmission over MIMO channels
CN101572581A (zh) 一种信干噪比确定方法及装置
CN110365377B (zh) 多天线空分多址和scma非正交多址结合的下行传输方法
CN102611534B (zh) 一种码本变换更新的方法和系统
KR101049113B1 (ko) 이동 통신 시스템의 전송 다이버시티 방법
CN114050853B (zh) 基于联合非正交码本与预编码设计的多用户mimo传输方法

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130109

Termination date: 20170630