CN102468845A - 一种消除数控振荡器的频率偏移的方法 - Google Patents

一种消除数控振荡器的频率偏移的方法 Download PDF

Info

Publication number
CN102468845A
CN102468845A CN201010532092XA CN201010532092A CN102468845A CN 102468845 A CN102468845 A CN 102468845A CN 201010532092X A CN201010532092X A CN 201010532092XA CN 201010532092 A CN201010532092 A CN 201010532092A CN 102468845 A CN102468845 A CN 102468845A
Authority
CN
China
Prior art keywords
controlled oscillator
frequency
digital controlled
common factor
nco
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
Application number
CN201010532092XA
Other languages
English (en)
Other versions
CN102468845B (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.)
TD Tech Ltd
Original Assignee
TD Tech Ltd
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 TD Tech Ltd filed Critical TD Tech Ltd
Priority to CN201010532092.XA priority Critical patent/CN102468845B/zh
Publication of CN102468845A publication Critical patent/CN102468845A/zh
Application granted granted Critical
Publication of CN102468845B publication Critical patent/CN102468845B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Stabilization Of Oscillater, Synchronisation, Frequency Synthesizers (AREA)

Abstract

本发明提供了一种消除数控振荡器的频率偏移的方法,该方法包括:根据数控振荡器中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子;从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率;根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器。通过使用上述的方法,可有效地消除NCO的频率偏移,实现完全无频偏的NCO,提高NCO的杂散性能,降低NCO的设计复杂度,节约***的资源,降低生产成本。

Description

一种消除数控振荡器的频率偏移的方法
技术领域
本发明涉及通信技术领域,尤其是指一种消除数控振荡器的频率偏移的方法。
背景技术
目前,随着通信技术的飞速发展,数字技术已经越来越多地应用到各种通信产品之中。而随着数字频率合成技术的发展,数控振荡器(NCO)由于具有频率分辨率高、切换时间短、相位噪声低等特点,已经在通信***中得到了非常广泛的应用。
在现有技术中,NCO是一种高精度的正、余弦数字信号发生器,可以将第一频率的输入信号转换成第二频率的输出信号,从而产生正交的正弦信号和余弦信号。一般来说,目前NCO的实现方法有三种:查表法、实时计算法以及上述两种方法的结合。以下将以查表法为例对NCO的结构和原理进行简单地介绍。
图1为现有技术中的数控振荡器的结构示意图。如图1所示,现有技术中的NCO一般包括:相位累加器101和波形存储器(ROM)102。NCO的工作原理一般为:首先根据各个正、余弦波相位计算好相位的正余弦函数值,并以相位为地址将上述正余弦函数值样本存储在ROM中;然后,相位累加器101接收外界输入的频率控制字,并在每一个时钟脉冲输入时,把所接收的频率控制字进行累加,得到每个时钟周期所对应的数字相位;接着,再以上述数字相位为地址,在ROM中所存储的正余弦函数值进行查找,得到与上述数字相位相对应的正余弦函数值,并最终输出与所查找的正余弦函数值相应的、正交的正弦信号和余弦信号。
在一般通信***中,实际使用中的NCO一般都存在一定的频率偏移(可简称为频偏,下同),其原因在于:由于对于实际使用中的NCO的频率精度的要求很高,因此,在对NCO进行设计时,NCO的最小频率分辨率一般都被设计得很小(例如,一般都小于0.1Hz)。当NCO的最小频率分辨率很小时,必然要求NCO中的ROM中最好是能存储有所有的与各种相位对应的正余弦函数值。然而,由于***资源有限,因此NCO中的ROM的存储容量不能太大,相位累加器的长度也不能很大,所以当NCO的最小频率分辨率很小时,ROM中不可能存储与所有可能的相位相对应的正余弦函数值,而只能存储一定数量的正余弦函数值。因此,此时NCO的输出信号的频率只能是与数字相位对应的相似的频率,而并不是准确的所需输出的频率,从而使得NCO的输出频率与所需输出的频率之间存在频率偏移。
由上可知,在现有技术中,实际使用中的NCO如果要达到***的要求,都难以有效地消除频率偏移,以实现完全无频偏的NCO,从而降低了***的性能,而且还需占用大量的***资源;同时,现有技术中的NCO由于最小分辨率很小,因此设计复杂度也很大,从而也提高了相应的生产成本。
发明内容
本发明提供了一种消除数控振荡器的频率偏移的方法,从而消除NCO的频率偏移,实现了完全无频偏的NCO。
为达到上述目的,本发明中的技术方案是这样实现的:
一种消除数控振荡器的频率偏移的方法,该方法包括:
根据数控振荡器中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子;
从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率;
根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器。
所述中频频率和本振频率为所述数控振荡器中所使用的最小中频频率和最小本振频率;所述***的输出频率为***的最小输出频率。
所述计算所述三种频率的公因子包括:计算所述三种频率的最大公因子。
所述从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率包括:
从所计算的公因子中选择任意一个公因子的值作为所述数控振荡器的最小频率分辨率。
所述从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率包括:
从所计算的公因子中选择最大公因子的值作为所述数控振荡器的最小频率分辨率。
所述根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器包括:
根据所述最小频率分辨率和所述数控振荡器的采样频率,计算一个完整的采样周期中的所有采样点的总数目;
在波形存储器中设置一个存储有预设数目的采样点数据的波形查找表,并根据所存储的采样点数据的数目设置所述相位累加器的长度。
所述预设数目的采样点数据为一个完整的采样周期中的所有的采样点数据。
所述预设数目的采样点数据为一个完整的采样周期中的角度在[0,π/4]区间内的采样点数据。
综上可知,本发明中提供了一种消除数控振荡器的频率偏移的方法。在所述消除数控振荡器的频率偏移的方法中,由于将根据数控振荡器中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子;然后将所计算的公因子的值作为数控振荡器的最小频率分辨率;并根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器,从而可有效地消除NCO的频率偏移,实现了完全无频偏的NCO,提高了NCO的杂散性能,降低了NCO的设计复杂度,节约了***的资源,降低了生产成本。
附图说明
图1为现有技术中的数控振荡器的结构示意图。
图2为本发明中的消除数控振荡器的频率偏移的方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点表达得更加清楚明白,下面结合附图及具体实施例对本发明再作进一步详细的说明。
本发明提供了一种消除数控振荡器的频率偏移的方法,在该方法中,将根据数控振荡器中所使用的中频频率、本振频率以及***中输出频率,计算所述三种频率的公因子;然后将所计算的公因子的值作为数控振荡器的最小频率分辨率;并根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器,从而有效地消除了NCO的频率偏移,实现了完全无频偏的NCO;此外,还可以大大降低NCO的设计复杂度,节约***的资源,从而降低生产成本。
图2为本发明中的消除数控振荡器的频率偏移的方法的流程示意图。如图2所示,在本发明中的消除数控振荡器的频率偏移的方法中,主要包括如下所述的步骤:
步骤201,根据NCO中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子。
在本步骤中,将首先获取NCO中所使用的中频频率、本振频率以及***的输出频率,然后再计算上述这三种频率的公因子。
例如,在TD-SCDMA***中,空口输出的频率并不是任意的。根据协议的要求,空口输出的频率必须是200KHz的整数倍,由此可知,***的最小输出频率为200KHz。例如,在一般的TD-SCDMA***中的射频拉远单元(RRU,RemoteRadio Unit)***中,从基带到射频的频率转换过程一般包括:先通过数字上变频(DUC)中的NCO将各个载波信号从基带的频谱上以零频为中心搬开,再通过一次频率变换将上述分离后的载波信号的中心频率从零频变换为中频频率fIF;由于通信***中一般采用带通采样,所以中频频率fIF的值一般为
Figure BSA00000332998300051
(其中,fs为采样频率);最后,再根据***的本振频率fL0将上述已变换为中频频率的载波信号变换成射频输出端的射频信号,该射频信号的输出频率frf。因此,最后的***输出频率frf可以用如下所述的公式表示:
frf=fnco+fIF+fL0                        (1)
其中,由于根据协议要求,TD-SCDMA***中***输出频率frf必须是200KHz的整数倍,因此必然有:frf=0.2a(MHz),其中,a为正整数;另外,在TD-SCDMA***中,中频频率fIF和***的本振频率fL0一般均为预先设置的固定值,例如,由于TD-SCDMA***中的码片速率为1.28MHz,因此fIF和fL0一般都被设置为1.28MHz的整数倍,因此必然有:fIF=1.28c(MHz)以及fL0=1.28d(MHz),其中,b和c均为正整数,也就是说,最小中频频率和最小本振频率均为1.28MHz。
所以,在本发明的具体实施例中,可根据上述的中频频率、本振频率和***的输出频率,计算这三种频率的公因子(即计算1.28c、1.28d和0.2a的公因子,单位为MHz),计算结果的单位可以是MHz、KHz或其它便于计算上述公因子的单位。例如,当***的输出频率、中频频率、本振频率分别为***的最小输出频率0.2MHz、最小中频频率1.28MHz和最小本振频率1.28MHz时(即a=c=d=1时),如果以1KHz为最小的单位,则通过计算得到的这三种频率的公因子将分别为40KHz、20KHz、10KHz、8KHz、5KHz、4KHz、2KHz和1KHz,其中,40KHz为最大公因子。
此外,在本发明的具体实施例中,所述计算所述三种频率的公因子也可以是只计算这三种频率的最大公因子,此时,上述实施例中的计算结果将为40KHz。
步骤202,从所计算的公因子中选择一个公因子的值作为所述NCO的最小频率分辨率。
在本步骤,可以从步骤201中所计算的上述三种频率的公因子中任意选择一个公因子的值作为数控振荡器的最小频率分辨率。其中,在本发明的较佳实施例中,所选择的公因子可以是最大公因子。另外,如果在步骤201中仅计算了上述三种频率的最大公因子,则在本步骤中可以直接将所计算的最大公因子的值作为数控振荡器的最小频率分辨率。
例如,当上述所计算的最大公因子的值为40KHz时,所述数控振荡器的最小频率分辨率即为40KHz,此时,根据上述公式(1)可得:
frf=0.04MHz*Arf=0.04MHz*(Bnco+CIF+DL0)=fnco+fIF+fL0     (2)
其中,Arf、Bnco、CIF和DL0均为正整数。
步骤203,根据所述最小频率分辨率和所述NCO的采样频率,设置所述NCO中的波形存储器和相位累加器。
由于数控振荡器的最小频率分辨率已经被确定,因此在本步骤中,可以根据上述NCO的最小频率分辨率以及NCO的采样频率,设置该NCO中的波形存储器和相位累加器。
具体来说,在本发明的具体实施例中,可先根据所述最小频率分辨率和所述数控振荡器的采样频率,计算一个完整的采样周期中的所有采样点的总数目;然后再在ROM中设置一个存储有预设数目的采样点数据的波形查找表,并根据所存储的采样点数据的数目设置相位累加器的长度。其中,在本发明的具体实施例中,上述预设数目的采样点数据可以是一个完整的采样周期中的所有的采样点数据。
例如,当NCO的最小频率分辨率为40KHz,采样频率为30.72MHz时,根据上述最小频率分辨率和采样频率可知,在一个完整的采样周期中,所有的采样点的总数目为:30.72/0.04=768(个)。因此,只需在ROM中设置一个存储有768个40KHz的采样点数据的波形查找表,并设置一个长度为n=Roundup(log2768)=10比特(bit)的相位累加器(其中,“Roundup”表示向上取整的运算),即可通过查表的方法使得该NCO可以输出具有40KHz的任意整数倍的频率的输出信号。具体的查表方法可以使用现有技术中常用的查表方法,在此不再赘述。
由上可知,由于上述NCO的最小频率分辨率远远大于现有技术中的NCO的最小频率分辨率(一般小于0.1Hz),一个完整采样周期中的采样点的总数目较少,因此可以在上述NCO的波形存储器中存储所有采样点的数据,从而使得该NCO的输出信号的频率必然是所需要的40KHz的整数倍的频率,而不会出现其它数值的频率,从而有效地消除了NCO的频率偏移,实现了完全无频偏的NCO;同时,也大大简化了该NCO的设计过程,大大降低NCO的设计复杂度;另外,由于一个完整采样周期中的采样点的总数目较小,因此波形存储器的存储空间也比较小,从而使得相位累加器的长度也可以大为减小,例如,当采样点的总数目为768时,相位累加器的长度可以仅为10比特(而现有技术中的相位累加器的长度一般为32比特以上),从而可以大大节约***的资源,降低生产成本。
进一步地,在本发明的具体实施例中,还可以利用正弦函数和余弦函数的特性,对波形存储器的存储空间的大小进行进一步的优化。例如,由于NCO可以同时输出正弦信号和余弦信号,而根据正弦函数和余弦函数的联合图象可知,在[0,π/4]的区间内,正、余弦函数的幅值都已经出现过一次,至于其它的角度在[π/4,2π]区间内的各个采样点的正弦值和余弦值,则可以通过相应的三角变换公式变换到[0,π/4]区间内的正弦值和余弦值,因此,上述波形存储器中的波形查找表中所存储的预设数目的采样点数据可以仅仅是角度在[0,π/4]区间内的采样点的正弦值和余弦值,也就是说,上述预设数目的采样点数据为一个完整的采样周期中的角度在[0,π/4]区间内的采样点数据;此时,所存储的采样点的数目为一个完整的采样周期中的所有采样点的总数目的八分之一。例如,当所有采样点的总数目为768时,所存储的采样点的数目为768/8=96(个)。如果使用一个16比特的无符号数来量化采样点的正弦值或余弦值,则该波形查找表所需要的存储空间仅为(96×32)bit。因此,相对于现有技术来说,本发明中的波形存储器所占用的***资源将非常小。
此外,由于本发明中的NCO可以输出完全无频偏的频率,因此,该NCO的最终输出误差仅仅取决于该NCO输出的量化位数,从而可使得该NCO具有很好的杂散性能。例如,当使用16比特的无符号数来进行量化时,所导致的量化误差为:
10 log ( 1 12 ( 2 - 16 ) 2 ) = - 107.12 dBFs
其中,dbFs表示以满刻度为零电平的分贝数,可简称为满刻度分贝值。例如,低于满刻度20db的信号便是-20dbFs。
相对于现有技术来说,上述量化误差远远小于现有技术中的输出误差,因此,本发明中的NCO还将具有很好的杂散性能。
综上可知,在本发明的实施例中提出了上述的消除数控振荡器的频率偏移的方法。在上述消除数控振荡器的频率偏移的方法中,由于将根据数控振荡器中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子;然后将所计算的公因子的值作为数控振荡器的最小频率分辨率;并根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器,从而有效地消除了NCO的频率偏移,实现了完全无频偏的NCO,并大大提高了NCO的杂散性能,因而大大提高了***的性能;同时,还可以大大降低NCO的设计复杂度,节约***的资源,从而降低生产成本。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明保护的范围之内。

Claims (8)

1.一种消除数控振荡器的频率偏移的方法,其特征在于,该方法包括:
根据数控振荡器中所使用的中频频率、本振频率以及***的输出频率,计算所述三种频率的公因子;
从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率;
根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器。
2.根据权利要求1所述的方法,其特征在于:
所述中频频率和本振频率为所述数控振荡器中所使用的最小中频频率和最小本振频率;所述***的输出频率为***的最小输出频率。
3.根据权利要求1所述的方法,其特征在于,所述计算所述三种频率的公因子包括:
计算所述三种频率的最大公因子。
4.根据权利要求1所述的方法,其特征在于,所述从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率包括:
从所计算的公因子中选择任意一个公因子的值作为所述数控振荡器的最小频率分辨率。
5.根据权利要求1所述的方法,其特征在于,所述从所计算的公因子中选择一个公因子的值作为所述数控振荡器的最小频率分辨率包括:
从所计算的公因子中选择最大公因子的值作为所述数控振荡器的最小频率分辨率。
6.根据权利要求1所述的方法,其特征在于,所述根据所述最小频率分辨率和所述数控振荡器的采样频率,设置所述数控振荡器中的波形存储器和相位累加器包括:
根据所述最小频率分辨率和所述数控振荡器的采样频率,计算一个完整的采样周期中的所有采样点的总数目;
在波形存储器中设置一个存储有预设数目的采样点数据的波形查找表,并根据所存储的采样点数据的数目设置所述相位累加器的长度。
7.根据权利要求6所述的方法,其特征在于:
所述预设数目的采样点数据为一个完整的采样周期中的所有的采样点数据。
8.根据权利要求6所述的方法,其特征在于:
所述预设数目的采样点数据为一个完整的采样周期中的角度在[0,π/4]区间内的采样点数据。
CN201010532092.XA 2010-10-29 2010-10-29 一种消除数控振荡器的频率偏移的方法 Expired - Fee Related CN102468845B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010532092.XA CN102468845B (zh) 2010-10-29 2010-10-29 一种消除数控振荡器的频率偏移的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010532092.XA CN102468845B (zh) 2010-10-29 2010-10-29 一种消除数控振荡器的频率偏移的方法

Publications (2)

Publication Number Publication Date
CN102468845A true CN102468845A (zh) 2012-05-23
CN102468845B CN102468845B (zh) 2013-10-16

Family

ID=46072091

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010532092.XA Expired - Fee Related CN102468845B (zh) 2010-10-29 2010-10-29 一种消除数控振荡器的频率偏移的方法

Country Status (1)

Country Link
CN (1) CN102468845B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104639160A (zh) * 2013-11-13 2015-05-20 上海联影医疗科技有限公司 数字控制振荡器及磁共振成像***
CN106095706A (zh) * 2015-04-08 2016-11-09 美高森美半导体无限责任公司 具有主设备和从设备的pll***

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001041381A2 (en) * 1999-12-02 2001-06-07 Silicon Wave, Inc. Gaussian filtered fsk modulator that includes a gaussian filter, an integrator and a look-up table
WO2008036613A1 (en) * 2006-09-18 2008-03-27 National Semiconductor Corporation Numerically controlled oscillator (nco) output clock phase smoothing
CN101276002A (zh) * 2007-03-26 2008-10-01 中国海洋石油总公司 高温单片相位可编程直接数字频率合成源
CN101807089A (zh) * 2010-04-02 2010-08-18 广西大学 输出信号偏移量任意可调的波形信号发生器

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001041381A2 (en) * 1999-12-02 2001-06-07 Silicon Wave, Inc. Gaussian filtered fsk modulator that includes a gaussian filter, an integrator and a look-up table
WO2008036613A1 (en) * 2006-09-18 2008-03-27 National Semiconductor Corporation Numerically controlled oscillator (nco) output clock phase smoothing
CN101276002A (zh) * 2007-03-26 2008-10-01 中国海洋石油总公司 高温单片相位可编程直接数字频率合成源
CN101807089A (zh) * 2010-04-02 2010-08-18 广西大学 输出信号偏移量任意可调的波形信号发生器

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈绍宽, 李明: "NCO的研究及其FPGA实现", 《仪表技术》, no. 5, 31 May 2006 (2006-05-31), pages 30 - 32 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104639160A (zh) * 2013-11-13 2015-05-20 上海联影医疗科技有限公司 数字控制振荡器及磁共振成像***
CN106095706A (zh) * 2015-04-08 2016-11-09 美高森美半导体无限责任公司 具有主设备和从设备的pll***
CN106095706B (zh) * 2015-04-08 2019-01-15 美高森美半导体无限责任公司 具有主设备和从设备的pll***

Also Published As

Publication number Publication date
CN102468845B (zh) 2013-10-16

Similar Documents

Publication Publication Date Title
CN101776935B (zh) 一种基于dds的数字调制信号发生器
CN102325058B (zh) 一种变频***群时延测试方法
CN101510756A (zh) 基于mimo实时测试平台的数字信号下变频处理***
CN103630743B (zh) 一种外差式频谱分析仪频率修正的方法
CN101572558B (zh) 一种中频收发芯片
CN101458329A (zh) 时域并行采样率变换方法
CN101908896B (zh) 一种多频段射频接收机
CN101378263A (zh) 基于数字中频的多载波数字接收机及多载波数字接收方法
CN102882517A (zh) 一种生成低失真度的低频正弦信号的装置及方法
CN104201991A (zh) 实现捷变频米波雷达的数字下变频***
CN107196673B (zh) 一种基于任意可变数字中频的信号分析装置及方法
US6690748B2 (en) Receiver with improved digital intermediate to base band demodulator
CN102751998A (zh) 一种基于软件无线电接收机的数字中频模块
CN106130571A (zh) 一种基于频域和空域压缩感知的信号采样接收方法及装置
CN102468845B (zh) 一种消除数控振荡器的频率偏移的方法
CN102368690A (zh) 微纳卫星测控数字中频与基带处理方法及装置
CN101814940B (zh) 数字中频光纤直放站的多信道数字选频数字信号处理方法
CN103916199A (zh) 一种天线信号的时延和相位调整装置与方法
CN201509196U (zh) 一种中频收发芯片
CN103856278B (zh) 调频数模同播的实现方法和装置
CN102594361A (zh) 音频异步采样率转换处理方法
CN209281395U (zh) 一种宽频率范围可重构信号处理平台
CN106972833A (zh) 一种可任意倍数重采样的数字上变频器
CN110198194B (zh) 无线电调制信号发生装置
CN101848038A (zh) 基于软件无线电滤波的无线传输装置及其方法

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

Termination date: 20211029