CN101557641B - 一种适用于蜂窝中继系统的子载波和功率分配方法 - Google Patents

一种适用于蜂窝中继系统的子载波和功率分配方法 Download PDF

Info

Publication number
CN101557641B
CN101557641B CN2009100983079A CN200910098307A CN101557641B CN 101557641 B CN101557641 B CN 101557641B CN 2009100983079 A CN2009100983079 A CN 2009100983079A CN 200910098307 A CN200910098307 A CN 200910098307A CN 101557641 B CN101557641 B CN 101557641B
Authority
CN
China
Prior art keywords
subcarrier
station
relay
access
link
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
CN2009100983079A
Other languages
English (en)
Other versions
CN101557641A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN2009100983079A priority Critical patent/CN101557641B/zh
Publication of CN101557641A publication Critical patent/CN101557641A/zh
Application granted granted Critical
Publication of CN101557641B publication Critical patent/CN101557641B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radio Relay Systems (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种适用于蜂窝中继系统的子载波和功率分配方法。基站通过专用控制信道收集系统中所有链路的信道状态信息,根据蜂窝中继系统的帧结构特点以及子载波和功率资源配置特点,分别完成下行接入域和下行中继域的子载波和功率联合分配,并通过专用控制信道将分配结果反馈给中继站,然后基站和中继站使用分配的子载波和发射功率完成数据传输。采用本发明所公开的方法,在下行接入域子载波分配过程中同时考虑了链路在子载波上的信道状态和接入站可供分配的发射功率,有效地提高了资源的利用效率;下行中继域在下行接入域分配基础上最大化系统有效吞吐量,系统总有效吞吐量有较大提高。

Description

一种适用于蜂窝中继系统的子载波和功率分配方法
技术领域
本发明涉及无线数字通信技术领域,具体涉及一种适用于蜂窝中继系统的子载波和功率分配方法。
背景技术
未来无线通信系统的一个显著特点是高速的数据传输能力。为了满足高速数据传输需求,未来无线通信系统分配的频率更高,高频电磁波的传播衰减大,小区覆盖范围相对现有网络变小,需要更多的基站满足系统吞吐量和覆盖的需求,由此就带来了高额的网络基础设备资金投入。为了解决这个问题,不仅需要采用更先进的物理层传输技术,更需要采用新的网络体系结构来支持。在传统蜂窝系统中引入智能中继技术,不仅能够提升用户的数据传输速率,还能有效扩大无线通信的覆盖范围,从而弥补现有系统不足,降低网络基础设施投资和运营成本,是一种极具吸引力和竞争力的未来无线通信网络架构。
蜂窝中继系统的网络元素,包括基站、中继站和移动用户,网络元素的增加,带来了更加复杂的网络拓扑结构。中继站将基站与用户之间的无线链路拆分为距离较短的两跳链路,有效地减少了总的路径损耗,但是由于中继链路的资源开销,如果不能对子载波和功率进行有效利用,则可能无法达到提升系统性能的初衷,针对这一问题已有大量文献就蜂窝中继系统中的子载波和功率分配方法进行了研究。现有方法,一般采用两步完成无线资源分配:子载波分配和功率调整。其中,对于子载波分配,采用接入站(包括基站和中继站)与用户之间链路在具体子载波上的信道增益作为子载波分配的决策依据,然后将发射功率在子载波分配结果上进行分配,并在问题建模中以所有接入站总发射功率大小受限作为约束条件。而实际系统中,基站与中继站在空间上彼此间隔一定距离分布在小区内,各自具有独立的电力能源供应,因此基站和中继站具有各自独立的发射功率资源,且无法实现共享。另一方面,小区内所有接入站都可以占用所有分配给该小区的子载波资源。在这样的资源结构特点条件下,仅以信道增益作为子载波分配决策依据,无法充分利用系统资源,考虑这样一种情况,对于某一信道条件足够好的接入站与用户之间的链路,在极端情况下则将会获得所有子载波,但由于该接入站发射功率有限,此时获得系统总吞吐量,与将部分子载波分配给其它信道条件较差链路的分配结果相比,后者更好地利用了系统的功率资源,也更可能获得较大的系统总吞吐量增益。所以,蜂窝中继系统的子载波和功率分配,需要同时考虑子载波的信道增益和可供分配的发射功率量。
发明内容
本发明的目的是为了解决上述方法中无线资源没有得到充分利用的问题,提供一种适用于蜂窝中继系统的子载波和功率分配方法。
适用于蜂窝中继系统的子载波和功率分配方法包括下行接入域分配阶段和下行中继域分配阶段:
下行接入域分配阶段包括如下步骤:
1)参数初始化
未分配子载波集合Nblank={1,2,…,N},接入站m,包括基站和中继,m=0表示基站,m≠0表示中继站,子载波分配集合
Figure G2009100983079D00021
d m , k , n = 0 , ∀ m , ∀ k , ∀ n ;
2)子载波分配包括如下步骤:
a)接入站m将子载波集合 N m available = N m ∪ N blank 内的子载波n分配给接入站m与用户k之间的链路(m,k)=argmaxgm,k,n,直至Nm available内子载波被全部分配,并将接入站总功率Pm在子载波分配结果基础上进行分配,即:
p m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线,gm,k,n为接入站m与用户k之间链路在子载波n上的信道增益,σ2为热噪声功率谱密度;
b)选择 ( m * , k * , n * ) = arg max ( m , k , n ) p m , k , n g m , k , n ;
c)将子载波n*分配给接入站m*与用户k*之间的链路,置 d m * , k * , n * = 1 , 同时更新未分配子载波集合Nblank=Nblank-n*,接入站m*子载波分配集合 N m * = N m * + n * ;
d)重复过程步骤a)到步骤c),直至
Figure G2009100983079D00028
时转入步骤3)。
3)功率分配
将接入站m的发射功率Pm分配到子载波分配集合Nm的所有子载波上,接入站m与用户之k间链路在子载波n上发射功率为:
P m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线;
4)中继站m,在下行中继域的速率需求:
r m = Σ k ∈ U m Σ n = 1 M d m , k , n B N log ( 1 + p m , k , n g m , k , n B N Γ σ 2 )
其中,Γ为信噪比差,Um为附着在中继站m上的用户集合,转入下行中继域分配阶段;
下行中继域分配阶段包括如下步骤:
5)参数初始化
未分配子载波集合Nblank={1,2,…,N},中继站m的子载波分配集合
Figure G2009100983079D00031
当前分配过程已经获得的速率Rm=0,基站的发射功率在每个子载波上平均分配;
6)子载波分配包括如下步骤:
e)中继站m使用未分配子载波集合Nblank内子载波n时,获得的传输速率:
R m , n = B N log ( 1 + p 0 , m , n g 0 , m , n B N Γ σ 2 )
f)将子载波n分配给基站与中继站 m * = arg max m ∈ { 1,2 , . . . , M } ( min ( R m , n , r m - R m ) ) 之间的链路;
g)更新Nblank=Nblank-n, N m * = N m * + n , Rm=Rm+Rm,n
h)重复步骤e)到步骤g),直至未分配子载波集合
Figure G2009100983079D00035
本发明给出了一种适用于蜂窝中继系统的子载波和功率分配方法,分别完成下行接入域和下行中继域的子载波分配。在子载波分配过程中,同时考虑了接入站与用户之间链路的信道增益和接入站可供分配的发射功率,提高了资源的利用效率,并在子载波分配基础上完成了发射功率的最优分配,实现了系统的吞吐量增益。
附图说明
图1是蜂窝中继系统网络下行传输部分结构图;
图2是本发明的性能曲线。
具体实施方式
下面根据附图以及一个实例对本发明所公开方法做一个详细的描述。
适用于蜂窝中继系统的子载波和功率分配方法包括下行接入域分配阶段和下行中继域分配阶段:
下行接入域分配阶段包括如下步骤:
1)参数初始化
未分配子载波集合Nblank={1,2,…,N},接入站m,包括基站和中继,m=0表示基站,m≠0表示中继站,子载波分配集合
Figure G2009100983079D00036
d m , k , n = 0 , ∀ m , ∀ k , ∀ n ;
2)子载波分配包括如下步骤:
a)接入站m将子载波集合 N m available = N m ∪ N blank 内的子载波n分配给接入站m与用户k之间的链路(m,k)=argmaxgm,k,n,直至Nm available内子载波被全部分配,并将接入站总功率Pm在子载波分配结果基础上进行分配,即:
p m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线,gm,k,n为接入站m与用户k之间链路在子载波n上的信道增益,σ2为热噪声功率谱密度;
b)选择 ( m * , k * , n * ) = arg max ( m , k , n ) p m , k , n g m , k , n ;
c)将子载波n*分配给接入站m*与用户k*之间的链路,置 d m * , k * , n * = 1 , 同时更新未分配子载波集合Nblank=Nblank-n*,接入站m*子载波分配集合 N m * = N m * + n * ;
d)重复过程步骤a)到步骤c),直至
Figure G2009100983079D00045
时转入步骤3)。
3)功率分配
将接入站m的发射功率Pm分配到子载波分配集合Nm的所有子载波上,接入站m与用户之k间链路在子载波n上发射功率为:
P m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线;
4)中继站m,在下行中继域的速率需求:
r m = Σ k ∈ U m Σ n = 1 M d m , k , n B N log ( 1 + p m , k , n g m , k , n B N Γ σ 2 )
其中,Γ为信噪比差,Um为附着在中继站m上的用户集合,转入下行中继域分配阶段;
下行中继域分配阶段包括如下步骤:
5)参数初始化
未分配子载波集合Nblank={1,2,…,N},中继站m的子载波分配集合
Figure G2009100983079D00048
当前分配过程已经获得的速率Rm=0,基站的发射功率在每个子载波上平均分配;
6)子载波分配包括如下步骤:
e)中继站m使用未分配子载波集合Nblank内子载波n时,获得的传输速率:
R m , n = B N log ( 1 + p 0 , m , n g 0 , m , n B N Γ σ 2 )
f)将子载波n分配给基站与中继站 m * = arg max m ∈ { 1,2 , . . . , M } ( min ( R m , n , r m - R m ) ) 之间的链路;
g)更新Nblank=Nblank-n, N m * = N m * + n , Rm=Rm+Rm,n
h)重复步骤e)到步骤g),直至未分配子载波集合
Figure G2009100983079D000412
实施例
蜂窝中继系统网络结构图如图1所示,基站部署在小区中央位置,M个中继站分布在小区内,并且摆放在距离基站2/3小区半径的位置上,K个用户均匀分布在小区内,根据附着接入站的不同,分别为通过直连传输方式接入基站的直连用户和通过中继传输方式接入中继站的中继用户。本实例中,蜂窝中继系统采用IEEE 802.16j协议草案中带内非透明中继模式的帧结构,小区内所有用户在下行接入域内接入系统,基站在下行中继域将中继用户的下行数据转发到中继用户附着的中继站;系统所使用的总带宽为B,分成N个子载波。基站的总发射功率为P0,中继站的总发射功率为Pm,m=1,…,M。在本实例中,M=6,K∈{1,5,10,15,20,25,30,35,40,45},N=96,P0=43dBm,Pm=33dBm,接收机热噪声功率谱密度σ2=-174dBm,信噪比差Γ≈5.48dB,在小区内随机产生K用户,并由路径损耗模型PLBS-MS=128.1+37.6logd(dB)和PLRS-MS=128.1+28.8logd(dB)计算接入站m与用户k之间链路的路径增益gm,k,n。根据本发明提供的方法,具体步骤如下:
下行接入域分配阶段包括如下步骤:
1)参数初始化
未分配子载波集合Nblank={1,2,…,N},接入站m,包括基站和中继,m=0表示基站,m≠0表示中继站,子载波分配集合
Figure G2009100983079D00051
d m , k , n = 0 , ∀ m , ∀ k , ∀ n ;
2)子载波分配包括如下步骤:
a)接入站m将子载波集合 N m available = N m ∪ N blank 内的子载波n分配给接入站m与用户k之间的链路(m,k)=argmaxgm,k,n,直至Nm available内子载波被全部分配,并将接入站总功率Pm在子载波分配结果基础上进行分配,即:
p m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线,gm,k,n为接入站m与用户k之间链路在子载波n上的信道增益,σ2为热噪声功率谱密度;
b)选择 ( m * , k * , n * ) = arg max ( m , k , n ) p m , k , n g m , k , n ;
c)将子载波n*分配给接入站m*与用户k*之间的链路,置 d m * , k * , n * = 1 , 同时更新未分配子载波集合Nblank=Nblank-n*,接入站m*子载波分配集合 N m * = N m * + n * ;
d)重复过程步骤a)到步骤c),直至
Figure G2009100983079D00058
时转入步骤3)。
3)功率分配
将接入站m的发射功率Pm分配到子载波分配集合Nm的所有子载波上,接入站m与用户之k间链路在子载波n上发射功率为:
P m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线;
4)中继站m,在下行中继域的速率需求:
r m = Σ k ∈ U m Σ n = 1 M d m , k , n B N log ( 1 + p m , k , n g m , k , n B N Γ σ 2 )
其中,Γ为信噪比差,Um为附着在中继站m上的用户集合,转入下行中继域分配阶段;
下行中继域分配阶段包括如下步骤:
5)参数初始化
未分配子载波集合Nblank={1,2,…,N},中继站m的子载波分配集合当前分配过程已经获得的速率Rm=0,基站的发射功率在每个子载波上平均分配;
6)子载波分配包括如下步骤:
e)中继站m使用未分配子载波集合Nblank内子载波n时,获得的传输速率:
R m , n = B N log ( 1 + p 0 , m , n g 0 , m , n B N Γ σ 2 )
f)将子载波n分配给基站与中继站 m * = arg max m ∈ { 1,2 , . . . , M } ( min ( R m , n , r m - R m ) ) 之间的链路;
g)更新Nblank=Nblank-n, N m * = N m * + n , Rm=Rm+Rm,n
h)重复步骤e)到g),直至未分配子载波集合
Figure G2009100983079D00067
图2中给出了本发明所提供的蜂窝中继系统子载波和功率分配方法在不同用户数时,系统总吞吐量性能变化情况,可以看出随着用户数目的增加,系统总吞吐量也随之增加,并且相对于对比方法,本发明所提供方法均取得了较明显的吞吐量增益,其中作为对比算法的是蜂窝中继系统基于子载波固定分配的方法。

Claims (1)

1.一种适用于蜂窝中继系统的子载波和功率分配方法,其特征在于包括下行接入域分配阶段和下行中继域分配阶段:
下行接入域分配阶段包括如下步骤:
1)参数初始化
未分配子载波集合Nblank={1,2,…,N},其中N为子载波的个数,接入站m,包括基站和中继站,m=0表示基站,m≠0表示中继站,子载波分配集合
Figure FSB00000234219400011
Figure FSB00000234219400012
2)子载波分配包括如下步骤:
a)接入站m将子载波集合
Figure FSB00000234219400013
内的子载波n分配给接入站m与用户k之间的链路(m,k)=argmaxgm,k,n,直至
Figure FSB00000234219400014
内子载波被全部分配,并将接入站总发射功率Pm在子载波分配结果基础上进行分配,即:
p m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中pm,k,n为接入站m与用户k之间链路在子载波n上的发射功率,λm为注水线,gm,k,n为接入站m与用户k之间链路在子载波n上的信道增益,σ2为热噪声功率谱密度;
b)选择 ( m * , k * , n * ) = arg max ( m , k , n ) p m , k , n g m , k , n ;
c)将子载波n*分配给接入站m*与用户k*之间的链路,置同时更新未分配子载波集合Nblank=Nblank-n*,接入站m*子载波分配集合
d)重复过程步骤a)到步骤c),直至
Figure FSB00000234219400019
时转入步骤3);
3)功率分配
将接入站m的总发射功率Pm分配到子载波分配集合Nm的所有子载波上,接入站m与用户k之间链路在子载波n上发射功率为:
P m , k , n = max ( λ m - σ 2 g m , k , n , 0 )
其中λm为注水线;
4)中继站m,在下行中继域的速率需求:
r m = Σ k ∈ U m Σ n = 1 N d m , k , n B N log ( 1 + P m , k , n g m , k , n B N Γ σ 2 )
其中,B为系统所使用的总带宽,Г为信噪比差,Um为附着在中继站m上的用户集合,转入下行中继域分配阶段;
下行中继域分配阶段包括如下步骤:
5)参数初始化
未分配子载波集合Nblank={1,2,…,N},中继站m的子载波分配集合当前分配过程已经获得的速率Rm=0,基站的发射功率在每个子载波上平均分配;
6)子载波分配包括如下步骤:
e)中继站m使用未分配子载波集合Nblank内子载波n时,获得的传输速率:
R m , n = B N log ( 1 + p 0 , m , n g 0 , m , n B N Γ σ 2 )
其中p0,m,n为基站与中继站m之间链路在子载波n上的发射功率,g0,m,n为基站与中继站m之间链路在子载波n上的信道增益;
f)将子载波n分配给基站与中继站
Figure FSB00000234219400023
之间的链路,M为中继站的个数;
g)更新Nblank=Nblank-n,
Figure FSB00000234219400024
Rm=Rm+Rm,n
h)重复步骤e)到步骤g),直至未分配子载波集合
Figure FSB00000234219400025
CN2009100983079A 2009-05-07 2009-05-07 一种适用于蜂窝中继系统的子载波和功率分配方法 Expired - Fee Related CN101557641B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009100983079A CN101557641B (zh) 2009-05-07 2009-05-07 一种适用于蜂窝中继系统的子载波和功率分配方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009100983079A CN101557641B (zh) 2009-05-07 2009-05-07 一种适用于蜂窝中继系统的子载波和功率分配方法

Publications (2)

Publication Number Publication Date
CN101557641A CN101557641A (zh) 2009-10-14
CN101557641B true CN101557641B (zh) 2010-12-08

Family

ID=41175521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009100983079A Expired - Fee Related CN101557641B (zh) 2009-05-07 2009-05-07 一种适用于蜂窝中继系统的子载波和功率分配方法

Country Status (1)

Country Link
CN (1) CN101557641B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101697637B (zh) * 2009-10-22 2012-11-07 上海交通大学 Ofdma系统中资源分配方法
US8812006B2 (en) * 2009-11-02 2014-08-19 Telefonaktiebolaget L M Ericsson (Publ) Methods and arrangements in a wireless communication system
CN101730109B (zh) * 2009-11-09 2012-06-06 中国人民解放军理工大学 基于博弈论的正交频分多址接入中继系统资源分配方法
CN102075998A (zh) * 2009-11-24 2011-05-25 中兴通讯股份有限公司 载波聚合系统中功率控制命令通知方法及系统
CN102123507A (zh) * 2010-01-08 2011-07-13 中兴通讯股份有限公司 一种载波聚合系统中载波的使用方法和装置
EP2638751B1 (en) 2010-11-09 2015-01-14 Telefonaktiebolaget L M Ericsson (PUBL) Method and apparatuse for composite reporting of wireless relay power capability
CN102186213B (zh) * 2011-05-05 2013-10-23 浙江大学 基于异构中继无线网络系统架构的上行无线资源管理方法
CN103298154B (zh) * 2013-05-31 2016-04-13 天津大学 一种基于温差电池的无线传感器网络动态能量分配方法
CN107659967A (zh) * 2017-08-25 2018-02-02 浙江工业大学 一种集能型无线中继网络吞吐量最大化的可再生能量优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060209671A1 (en) * 2005-03-18 2006-09-21 Samsung Electronics Co., Ltd. System and method for subcarrier allocation in a wireless multihop relay network
US20080117854A1 (en) * 2006-10-16 2008-05-22 Nokia Corporation Bandwidth allocation for relay networks
CN101399799A (zh) * 2007-09-25 2009-04-01 中兴通讯股份有限公司 一种正交频分复用系统的数据中继方法及中继站

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060209671A1 (en) * 2005-03-18 2006-09-21 Samsung Electronics Co., Ltd. System and method for subcarrier allocation in a wireless multihop relay network
US20080117854A1 (en) * 2006-10-16 2008-05-22 Nokia Corporation Bandwidth allocation for relay networks
CN101399799A (zh) * 2007-09-25 2009-04-01 中兴通讯股份有限公司 一种正交频分复用系统的数据中继方法及中继站

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
余官定等.OFDM中继系统的中继子载波对选择和功率分配.电子与信息学报.2007,29(2),265-268. *

Also Published As

Publication number Publication date
CN101557641A (zh) 2009-10-14

Similar Documents

Publication Publication Date Title
CN101557641B (zh) 一种适用于蜂窝中继系统的子载波和功率分配方法
Siddique et al. Wireless backhauling of 5G small cells: Challenges and solution approaches
Pabst et al. Relay-based deployment concepts for wireless and mobile broadband radio
CN101248691B (zh) 用于控制从线性或树型拓扑结构网络的节点接入tdma无线信道的方法
CN103442369B (zh) 一种广域覆盖场景下的灵活组网方法
CN104284407B (zh) 嵌入d2d的蜂窝网络中基于全双工中继的功率控制方法
CN101971521B (zh) 分布式转发器的布置方法及其系统
CN101420702A (zh) 中继系统部署方法和中继系统
CN101505482A (zh) 半分布式资源分配方法及系统
CN107222907A (zh) 基于上行非正交多址中继辅助的蜂窝网络传输方法
CN102820954A (zh) 一种降低异构网络小区间干扰的方法
CN101951645B (zh) 一种蜂窝中继网络中的下行自适应传输方法
CN103369568B (zh) Lte-a中继系统中基于博弈论的无线资源优化方法
Ancans et al. Analysis of characteristics and requirements for 5G mobile communication systems
CN103701737B (zh) 一种云无线接入网络通信的干扰消除方法
CN107241799A (zh) 异构网络中能量有效的资源分配和干扰对齐联合方法
CN106358204A (zh) 一种蜂窝异构网络下的全双工小区间干扰协调方法
Celik et al. On the feasibility of blind dynamic TDD in ultra-dense wireless networks
CN106102173A (zh) 基于多播波束成形的无线回传及基站分簇联合优化方法
CN101389128B (zh) 一种用户终端通过无线中继接入网络的方法
Psilopanagiotis et al. On relay-based subcarrier allocation and power management in 5G multicellular networks
CN103379648B (zh) 中继协调调度方法和装置
CN106130609B (zh) 一种波束通信的协作增强传输方法
CN101635973A (zh) 中继增强型蜂窝系统及其路由选择与比例公平性调度方法
Ayyadurai et al. Multihop cellular network optimization using genetic algorithms

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: 20101208

Termination date: 20120507