CN103338036A - 基于相位群处理的原子钟频率信号链接控制方法 - Google Patents

基于相位群处理的原子钟频率信号链接控制方法 Download PDF

Info

Publication number
CN103338036A
CN103338036A CN2013102945186A CN201310294518A CN103338036A CN 103338036 A CN103338036 A CN 103338036A CN 2013102945186 A CN2013102945186 A CN 2013102945186A CN 201310294518 A CN201310294518 A CN 201310294518A CN 103338036 A CN103338036 A CN 103338036A
Authority
CN
China
Prior art keywords
signal
frequency
output
atomic clock
phase
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
CN2013102945186A
Other languages
English (en)
Other versions
CN103338036B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN201310294518.6A priority Critical patent/CN103338036B/zh
Publication of CN103338036A publication Critical patent/CN103338036A/zh
Application granted granted Critical
Publication of CN103338036B publication Critical patent/CN103338036B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

本发明公开了一种基于相位群处理的原子钟频率信号链接控制方法,该控制方法分两路闭环控制完成:第一路为内环闭环控制,用于铯原子钟中9192631770Hz激励信号的产生并使晶体振荡器输出稳定和准确的信号;第二路为外环闭环控制,利用10MHz信号的分频输出与14.591479MHz高稳晶振的信号进行群相位同步和连续的处理实现最终的高精度10MHz振荡器信号输出。该控制方法不但具有很好的准确度指标而且频标输出信号的短期稳定度和相位噪声指标也完全保留了10MHz高稳定度晶体振荡器的高指标。

Description

基于相位群处理的原子钟频率信号链接控制方法
技术领域
本发明属于频率信号链接技术领域,尤其涉及一种基于相位群处理的原子钟频率信号链接控制方法。
背景技术
频率链接主要研究的是差异很大的频率信号之间实现相位或频率的准确和稳定的传递,所体现出的是利用相位群同步现象和相位群连续的特点,通过相位群处理实现频率信号之间的精密链接。目前这种频率链接已经在原子钟性能改进中得到了实现。能通过广义相位处理的频率链接方法来改造钟的系统结构并且改善其性能指标。目前,无论是主动型原子钟还是被动型原子钟,都是在频率处理的基础上,采用频率归一化进行处理的。
现有技术复杂的频率变换导致了噪声的增加,无法既保留原子钟本身的高稳定度指标,又充分发挥晶振的优良指标。
发明内容
本发明的目的在于提供一种基于相位群处理的原子钟频率信号链接控制方法,旨在解决现有技术因频率变换复杂而导致的噪声增加和某些指标下降的问题。
本发明是这样实现的,一种基于相位群处理的原子钟频率信号链接控制方法,分两路闭环控制完成:
第一路为内环闭环控制,用于铯原子钟中激励信号的产生并使晶体振荡器输出准确和稳定的信号;
第二路为外环闭环控制,利用信号的分频输出与高稳晶振的信号进行群相位同步和连续的处理实现最终的振荡器信号输出。
进一步,第一路为内环闭环控制,用于铯原子钟中9192631770Hz激励信号的产生并使晶体振荡器输出稳定的信号;
第二路为外环闭环控制,用于实现最终的高精度10MHz振荡器信号输出。
进一步,第一路内环闭环控制的实现方法如下:
14.591479MHz压控晶体振荡器(14.591479MHzVCOCXO)输出信号经过630倍倍频和调相后得到铯原子的特定的能级跃迁频率的激励信号,即9192631770Hz的激励信号;
在9192631770Hz的激励信号激励下,铯束管内的铯原子产生相应的能级跃迁,该跃迁信号的频率与激励信号的频率完全相等时有最大输出信号。否则,在探测器收到的跃迁信号偏离最大值时,会有一个误差电压输出给14.591479MHz压控晶体振荡器,调节其振荡频率至出现跃迁信号的最大值,使14.591479MHz晶体振荡器输出一个稳定的频率。
进一步,第二路外环闭环控制的实现方法如下:
10MHz压控晶体振荡器(10MHzVCOCXO)经过1/10000分频的信号与14.591479MHz的晶体振荡器输出进行相位群同步,实现了异频检相。这两个信号间的最小公倍数周期Tminc为1秒,在这一秒内测量1000次的相位信息,将这个信息被送入CPU;
CPU计算出对10MHz压控晶体振荡器的压控电压,达到对其调整和锁定的目的,具体实施方法为:
在1秒内有1000次对信号间的相位差进行测量,每次测量和1秒后的那一次测量间的间隔为一个最小公倍数周期,即1秒。这1000个成组的数据中有一些数据可能超出时间测量模块的测量范围,按照每一个Tminc根据相位差采集的结果用公式算出10MHz压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T,其中ΔT是每组对应数据间隔T的相位差的变化值,T是Tminc或者其倍数值。在所有的测量结果中剔除那些超出可测范围的、测量过程出现局部满周期交替的情况、以及噪声引起的起伏数值太大的情况等。将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差。从这个频差值和10MHzVCOCXO的压控灵敏度等,计算出对该振荡器的压控电压,达到对其调整和锁定的目的。
CPU的输出进行D/A转换,成为电压信号以调节10MHz压控晶体振荡器的输出。
进一步,内环采用了对14.591479MHz高稳定度晶体振荡器的锁定环路,对其输出信号进行630倍频输出9192.631770MHz的信号;它用作为铯束管的激励信号,产生的检测信号控制14.591479MHz高稳定度晶体振荡器的输出,将其锁定在铯原子能级跃迁对应信号的高稳定度和准确度上;
外环则是用被铯原子跃迁信号准确度和稳定度锁定的14.591479MHz振荡器信号来锁定10MHz振荡器信号,10MHz信号经1/10000分频后的1ms信号与14.591479MHz进行比对,测量其群相位差的变化,每秒测量1000次,在10MHz信号1/10000乘以1/1000分频之后,得到1s信号,以此1s信号作为时间标记控制CPU对其内的1000个相位差测量数据进行处理;
按照每一个Tminc内相位差的出现规律用公式算出压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T;其中ΔT是每组对应数据间隔T的相位差的变化值,T是Tminc或者其倍数值;其中在测量结果中剔除那些超出时间测量模块可测范围的、局部满周期交替的情况以及噪声引起的起伏数值太大的情况;将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差;从频差值和10MHzVCOCXO的压控灵敏度,计算出对该振荡器的压控电压,达到对其调整和锁定的结果;其中,Tminc为最小公倍数周期,若fx=Afmacx,fref=Bfmacx,则,Tminc=1/fmacx.将Tminc定义为一个群周期,在Tminc内包含多个比对的相位差。
本发明通过第一路内环闭环控制,铯原子能级跃迁的高稳定度传递到14.591479MHz压控晶体振荡器输出信号;通过第二路外环闭环控制,使得输出信号的短期稳定度和相位噪声指标也完全保留了10MHz高稳定度晶体振荡器的高指标。
附图说明
图1是本发明提供的基于相位群处理的原子钟频率信号链接控制方法原理实现框图;
图2是本发明提供的基于相位群处理的原子钟频率信号链接控制方法的鉴相结果示意图。
具体实施方式
本发明是这样实现的,一种基于相位群处理的原子钟频率信号链接控制方法,分两路闭环控制完成,如图1所示:
第一路为内环闭环控制,用于铯原子钟中9192631770Hz激励信号的产生并使14.591479MHz晶体振荡器输出稳定的信号;
第二路为外环闭环控制,用于实现最终的高精度10MHz振荡器信号输出。
进一步,第一路内环闭环控制的实现方法如下:
14.591479MHz压控晶体振荡器(14.591479MHzVCOCXO)输出信号经过630倍倍频和调相后得到铯原子的特定的能级跃迁频率的激励信号,即9192631770Hz的激励信号;
在9192631770Hz的激励信号的激励下,铯束管内的铯原子产生相应的能级跃迁,该跃迁信号的频率与激励信号的频率完全相等时有最大输出信号。否则,在探测器收到的跃迁信号偏离最大值时,会有一个误差电压输出给14.591479MHz压控晶体振荡器,调节其振荡频率至出现跃迁信号的最大值,使14.591479MHz晶体振荡器输出一个稳定的频率。
进一步,第二路外环闭环控制的实现方法如下:10MHz压控晶体振荡器(10MHzVCOCXO)经过1/10000分频后的信号与14.591479MHz的晶体振荡器输出进行相位群同步,实现了异频检相。这两个信号间的最小公倍数周期Tminc为1秒,在这一秒内测量1000次相位信息,将这个信息送入CPU;
CPU计算出对10MHz压控晶体振荡器的压控电压,达到对其调整和锁定的目的,具体实施方法为:在1秒内有1000次对信号间的相位差进行测量,每次测量和1秒后的那一次测量间的间隔为一个最小公倍数周期,即1秒。这1000个成组的数据中有一些数据可能超出时间测量模块的测量范围,按照每一个Tminc根据相位差采集的结果用公式算出10MHz压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T,其中ΔT是每组对应数据间隔T的相位差的变化值,T是Tminc或者其倍数值。在所有的测量结果中剔除那些超出可测范围的、测量过程出现局部满周期交替的情况、以及噪声引起的起伏数值太大的情况等。将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差。从这个频差值和10MHzVCOCXO的压控灵敏度等,计算出对该振荡器的压控电压,达到对其调整和锁定的结果。
CPU的输出进行D/A转换,成为电压信号以调节10MHz压控晶体振荡器的输出。
14.591479MHz压控晶体振荡器(14.591479MHzVCOCXO)输出信号经过630倍倍频和调相后得到铯原子的特定的能级跃迁频率的激励信号,即9192631770Hz的激励信号;在9192631770Hz的激励信号电磁波激励下,铯束管内的铯原子产生相应的能级跃迁,该跃迁信号的频率与激励信号的频率完全相等时有最大输出信号。否则,在探测器收到的跃迁信号偏离最大值时,会有一个误差电压输出给14.591479MHz压控晶体振荡器,调节其振荡频率至出现跃迁信号的最大值,使14.591479MHz晶体振荡器输出一个稳定的频率。10MHz压控晶体振荡器(10MHzVCOCXO)依次经过1/10000分频后的1ms信号与14.591479MHz的晶体振荡器输出进行相位群同步,也就是实现了异频检相,并以相关的1s信号做做测量时间,这两个信号间的最小公倍数周期Tminc为1秒,这个信息被送入CPU;在1秒内有1000次对信号间的相位差进行测量,每次测量和1秒后的那一次测量间的间隔为一个最小公倍数周期,即1秒。这1000个成组的数据中有一些数据可能超出时间测量模块的测量范围,按照每一个Tminc根据相位差采集的结果用公式算出10MHz压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T,其中ΔT是每组对应数据间隔T的相位差的变化值,T是Tminc或者其倍数值。在所有的测量结果中剔除那些超出可测范围的、测量过程出现局部满周期交替的情况、以及噪声引起的起伏数值太大的情况等。将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差。从这个频差值和10MHzVCOCXO的压控灵敏度等,计算出对该振荡器的压控电压,即CPU的输出,达到对其调整和锁定的结果。CPU的输出进行D/A转换,成为电压信号以调节10MHz压控晶体振荡器的输出。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (7)

1.一种基于相位群处理的原子钟频率信号链接控制方法,其特征在于,所述基于相位群处理的原子钟频率信号链接控制方法,分两路闭环控制完成:
第一路为内环闭环控制,用于铯原子钟中激励信号的产生并使晶体振荡器输出准确和稳定的信号;
第二路为外环闭环控制,利用信号的分频输出与高稳晶振的信号进行群相位同步和连续的处理实现最终的振荡器信号输出。
2.如权利要求1所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,第一路为内环闭环控制,用于铯原子钟中9192631770Hz激励信号的产生并使14.591479MHz晶体振荡器输出准确和稳定的信号;
第二路为外环闭环控制,利用10MHz信号的分频输出与14.591479MHz高稳晶振的信号进行群相位同步和连续的处理实现最终的10MHz振荡器信号输出。
3.如权利要求2所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,铯原子钟物理部分的激励信号的产生,只是直接倍频;所以引起的附加噪声很小;所述第一路内环闭环控制的实现方法如下:
14.591479MHz压控晶体振荡器输出信号经过630倍倍频和调相后得到铯原子的特定的能级跃迁频率的激励信号,即9192631770Hz的激励信号;
在9192631770Hz的激励信号激励下,铯束管内的铯原子产生相应的能级跃迁,该跃迁信号的频率与激励信号的频率完全相等时有最大输出信号;否则,在探测器收到的跃迁信号偏离最大值时,会有一个误差电压输出给14.591479MHz压控晶体振荡器,调节其振荡频率至出现跃迁信号的最大值,使14.591479MHz晶体振荡器输出一个稳定的频率。
4.如权利要求2所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,14.591479MHz晶体振荡器输出信号与10MHz压控晶体振荡器信号之间的锁相控制采用的是不同频率信号之间的不经过复杂频率变换的直接锁相控制;所述第二路外环闭环控制的实现方法如下:
压控晶体振荡器经过1/10000分频后作为开门信号与晶体振荡器输出关门信号进行相位群同步,实现了异频检相;这两个信号间的最小公倍数周期Tminc为1s,并以相关的1s信号作为测量时间标记,最后测量信息被送入CPU。
5.如权利要求4所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,CPU计算出对压控晶体振荡器的压控电压,达到对其调整和锁定的结果,具体实施方法为:
在1秒内有1000次对信号间的相位差进行测量,每次测量和1秒后的那一次测量间的间隔为1秒;这1000个成组的数据中有一些数据可能超出时间测量模块的测量范围,按照每一个Tminc根据相位差采集的结果用公式算出压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T,其中ΔT是每组对应数据间隔τ的相位差的变化值,T是Tminc或者其倍数值;在测量结果中剔除那些超出可测范围的、测量过程出现局部满周期交替的情况、以及噪声引起的起伏数值太大的情况;将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差;从这个频差值和10MHzVCOCXO的压控灵敏度,计算出对该振荡器的压控电压,达到对其调整和锁定的结果。
6.如权利要求5所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,对CPU的输出进行D/A转换,成为电压信号以调节10MHz压控晶体振荡器的输出。
7.如权利要求2所述基于相位群处理的原子钟频率信号链接控制方法,其特征在于,内环采用了对14.591479MHz高稳定度晶体振荡器的锁定环路,对其输出信号进行630倍频输出9192.631770MHz的信号,它用作为铯束管的激励信号,铯束管产生的检测信号控制14.591479MHz高稳定度晶体振荡器的输出,将其锁定在铯原子能级跃迁对应信号的高稳定度和准确度上;
外环则是用被铯原子跃迁信号准确度和稳定度锁定的14.591479MHz振荡器信号来锁定10MHz振荡器信号,10MHz信号经1/10000分频后为1ms信号与14.591479MHz进行比对,测量其群相位差的变化,每秒测量1000次;1ms信号再经过1/1000分频之后,得到1s信号,以此1s信号作为时间标记控制CPU对其内的1000个相位差测量数据进行处理;
按照每一个Tminc内相位差的出现规律用公式算出压控晶体振荡器相对于其标称值的相对平均频率变化偏差为Δf/f=ΔT/T;其中ΔT是每组对应数据间隔T的相位差的变化值,T是Tminc或者其倍数值;其中在测量结果中剔除那些超出时间测量模块可测范围的、局部满周期交替的情况以及噪声引起的起伏数值太大的情况;将其他有效的测量数据先按照式Δf/f=ΔT/T分别计算出各个相对频差,然后按照需要的平均时间值求得对应的频差;从这个频差值和10MHzVCOCXO的压控灵敏度,计算出对该振荡器的压控电压,达到对其调整和锁定的结果;其中,Tminc为最小公倍数周期,若fx=Afmacx,fref=Bfmacx,则,Tminc=1/fmacx.将Tminc定义为一个群周期,在Tminc内包含多个比对的相位差。
CN201310294518.6A 2013-07-12 2013-07-12 基于相位群处理的原子钟频率信号链接控制方法 Expired - Fee Related CN103338036B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310294518.6A CN103338036B (zh) 2013-07-12 2013-07-12 基于相位群处理的原子钟频率信号链接控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310294518.6A CN103338036B (zh) 2013-07-12 2013-07-12 基于相位群处理的原子钟频率信号链接控制方法

Publications (2)

Publication Number Publication Date
CN103338036A true CN103338036A (zh) 2013-10-02
CN103338036B CN103338036B (zh) 2016-04-13

Family

ID=49246150

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310294518.6A Expired - Fee Related CN103338036B (zh) 2013-07-12 2013-07-12 基于相位群处理的原子钟频率信号链接控制方法

Country Status (1)

Country Link
CN (1) CN103338036B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107493102A (zh) * 2017-09-18 2017-12-19 周渭 一种对被动型原子钟数字化锁相处理的新方法
CN109743022A (zh) * 2019-01-09 2019-05-10 西安电子科技大学 一种改善晶体振荡器老化漂移率的方法
CN110988463A (zh) * 2019-11-07 2020-04-10 西安电子科技大学 数字化相位比对精准获取信号频率及频率稳定度的方法
CN112702060A (zh) * 2020-12-11 2021-04-23 兰州空间技术物理研究所 一种用于铯原子钟的微波幅度控制方法
CN113395069A (zh) * 2021-06-15 2021-09-14 郑州轻大产业技术研究院有限公司 一种基于模糊区脉冲检测的高精度异频数字锁相环系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0351780A2 (en) * 1988-07-19 1990-01-24 Fujitsu Limited Cesium oscillator
CN201937572U (zh) * 2011-02-22 2011-08-17 合肥威师智能电子电器厂 一种cpt原子钟伺服电路
CN202998067U (zh) * 2012-12-17 2013-06-12 江汉大学 基于双微波探测系统频标

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0351780A2 (en) * 1988-07-19 1990-01-24 Fujitsu Limited Cesium oscillator
CN201937572U (zh) * 2011-02-22 2011-08-17 合肥威师智能电子电器厂 一种cpt原子钟伺服电路
CN202998067U (zh) * 2012-12-17 2013-06-12 江汉大学 基于双微波探测系统频标

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周渭等: "《周期性信号之间的相位量子特征及应用》", 《西安电子科技大学学报(自然科学版)》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107493102A (zh) * 2017-09-18 2017-12-19 周渭 一种对被动型原子钟数字化锁相处理的新方法
CN107493102B (zh) * 2017-09-18 2020-06-30 周渭 一种对被动型原子钟数字化锁相处理的方法
CN109743022A (zh) * 2019-01-09 2019-05-10 西安电子科技大学 一种改善晶体振荡器老化漂移率的方法
CN110988463A (zh) * 2019-11-07 2020-04-10 西安电子科技大学 数字化相位比对精准获取信号频率及频率稳定度的方法
CN112702060A (zh) * 2020-12-11 2021-04-23 兰州空间技术物理研究所 一种用于铯原子钟的微波幅度控制方法
CN113395069A (zh) * 2021-06-15 2021-09-14 郑州轻大产业技术研究院有限公司 一种基于模糊区脉冲检测的高精度异频数字锁相环系统
CN113395069B (zh) * 2021-06-15 2023-03-03 郑州轻大产业技术研究院有限公司 一种基于模糊区脉冲检测的高精度异频数字锁相环系统

Also Published As

Publication number Publication date
CN103338036B (zh) 2016-04-13

Similar Documents

Publication Publication Date Title
CN102045062B (zh) 一种基于Cordic算法的数字锁相环
CN103338036A (zh) 基于相位群处理的原子钟频率信号链接控制方法
CN104199278B (zh) 基于多导航系统的抗遮挡的高精度同步时钟系统及其同步方法
CN105091835B (zh) 一种转台动态角位置误差测量精度提高装置及校正方法
CN105049040B (zh) 一种利用gnss校正相干布居囚禁原子钟输出频率的方法
CN102540216A (zh) 一种自适应跟踪环路及实现方法
CN203377841U (zh) 基于卫星授时的晶振驯服设备
CN106885998B (zh) 一种提高铯光泵磁共振信号测频精度和速度的方法和电路
US11965919B2 (en) Phase frequency detector-based high-precision feedback frequency measurement apparatus and method
CN104090160A (zh) 一种高精度频率测量装置
CN103983957A (zh) 一种多普勒偏移测量方法及其装置
Roncagliolo et al. Optimized Carrier Tracking Loop Design for Real‐Time High‐Dynamics GNSS Receivers
CN105182069A (zh) 一种异频架构下的高分辨群量子化相位处理方法
CN110995250A (zh) 一种高分辨率任意频率信号的直接数字锁相电路及方法
CN104076373A (zh) 一种基于多信息融合辅助的载波跟踪方法与系统
CN105021904A (zh) 一种基于dds移相技术的快速相位噪声测量系统及测量方法
CN114578679A (zh) 一种基于授时压控技术应用于隧道的时间同步方法
CN202075347U (zh) 一种用于原子频标的环路振荡周期测量设备
CN106291102A (zh) 一种频标比对装置及方法
CN109856956A (zh) 基于比对的时钟授时装置
CN212850459U (zh) 一种基于fpga的高精度数字锁相环系统
CN205377841U (zh) 时间标准频率源
CN101207436A (zh) 一种反馈延时锁相装置及方法以及相位误差检测单元
CN109120260A (zh) 一种基于asic-tdc的时钟模块高精度鉴相系统及方法
CN207133596U (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

Granted publication date: 20160413

CF01 Termination of patent right due to non-payment of annual fee