CN101975584B - 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法 - Google Patents

一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法 Download PDF

Info

Publication number
CN101975584B
CN101975584B CN2010102715778A CN201010271577A CN101975584B CN 101975584 B CN101975584 B CN 101975584B CN 2010102715778 A CN2010102715778 A CN 2010102715778A CN 201010271577 A CN201010271577 A CN 201010271577A CN 101975584 B CN101975584 B CN 101975584B
Authority
CN
China
Prior art keywords
sampling
signal
add
sampled value
unit
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
CN2010102715778A
Other languages
English (en)
Other versions
CN101975584A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN2010102715778A priority Critical patent/CN101975584B/zh
Publication of CN101975584A publication Critical patent/CN101975584A/zh
Application granted granted Critical
Publication of CN101975584B publication Critical patent/CN101975584B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Gyroscopes (AREA)

Abstract

本发明公开了一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,本测试方法采用数字式逐点平均扫描的方法恢复其原始波形,将调制串扰信号经过前置滤波电路放大,A/D转换器转换成数字干扰信号后,进入中心处理器中进行处理并传送到信息融合计算,计算结果通过D/A转换器转换成模拟信号后输出,即为恢复出的干扰原始波形。此种测量方法能够有效地测量干涉式光纤陀螺内部的周期性串扰信号。

Description

一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法
技术领域
本发明涉及一种对误差信号的测试方法,更特别地说,是指一种适用于干涉式光纤陀螺的信号处理装置(检测电路)中调制串扰误差的开环测试方法。
背景技术
干涉式光纤陀螺是一种测量角速度的仪器,其硬件包括光源1、耦合器2、Y波导3、光纤环4、光电探测器5和信号处理装置6(也称角速度信息的检测电路)组成(请参见图1所示)。所述的信号处理装置6包括用于对光电探测器5输出的光功率信号f1进行滤波放大的前放滤波电路61,然后滤波放大的光功率信号顺次经A/D转换器62、中心处理器63、D/A转换器64和放大调理电路65后输出调制信号作用到Y波导3上(请参见图2所示);所述的晶振电路66提供中心处理器63所需的时钟信息;所述中心处理器63可以由DSP芯片实现、也可以由FPGA芯片实现、也可以由DSP芯片+FPGA芯片实现。干涉式光纤陀螺对角速度的测量是通过在光纤环4中传播的两束相向的光在光纤陀螺自身的转动中,引起的非互易相位差的大小来表征的。陀螺是一种敏感相对于惯性空间角运动的装置。它作为一种重要的惯性敏感器,用于测量运载体的姿态角和角速度,是构成惯性系统的核心器件。可以应用在飞行器导航、舰船导航和陆用导航中。
干涉式光纤陀螺环形干涉仪中,光波在Y波导3和耦合器2之间的两路光路的群传输时间之差倒数的二分之一称为光纤陀螺的本征频率(eigen frequency)。光纤陀螺最小互易性结构的光功率响应是一个隆起的余弦函数,为了获得较高的灵敏度,故给该信号施加一个偏置,使之工作在一个响应斜率不为零的点附近。而Y波导3中的寄生非线性或振幅调制可能会削弱偏置的质量。在Y波导3非线性的情况下,一种简单的解决方法是使光纤陀螺环工作在本征频率(或其奇次谐波)上,因此,光纤陀螺的信号处理装置6通常都是基于其本征频率来设计其控制时序的。
在数字闭环光纤陀螺中,调制、解调技术是提高光纤陀螺精度、改善陀螺抗干扰能力的重要手段,闭环反馈控制技术是提高陀螺动态范围、改善标度因数线性度的重要方法,但这两项技术也导致信号检测电路设计复杂,引入多项检测误差。其中最主要的一项检测误差为调制串扰误差,是由于调制信号通过光路和电路串扰到闭环前向通道中产生,是数字闭环光纤陀螺特有的一种误差。
发明内容
为了测试干涉式光纤陀螺的信号处理装置(检测电路)中的调制串扰信号f2并恢复其原始波形,本发明提出一种适用于干涉式光纤陀螺的检测电路调制串扰误差的开环测试方法。该开环测试方法采用数字式逐点平均扫描方式对调制串扰信号f2进行扫描,然后对获取的扫描信号进行累加平均的融合处理,最后复原得到调制串扰信号f2的原始波形。本发明方法利用了调制串扰信号f2的干扰源为调制方波同频的周期信号,因此能够有效地测量出光电转换后的探测器信号受到同频周期信号干扰的情况。
本发明的一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,所述检测电路(信号处理装置)包括有前放滤波电路、A/D转换器、中心处理器、D/A转换器、放大调理电路和晶振电路;前放滤波电路对光电探测器输出的光功率信号f1进行滤波放大处理,然后滤波放大后的光功率信号经A/D转换器后形成调制串扰信号f2输出给中心处理器;为了实现对调制串扰信号f2依照图3所示的周期性波形进行数字式逐点平均扫描,并将扫描获得的信息进行处理的步骤为:
第一步:采样时序控制单元控制时钟信息t来获得调制串扰信号f2进行数字式逐点平均扫描的采样时间点t0
第二步:在采样时间点t0上,对调制串扰信号f2在一个干扰信号的周期T内均匀采样M次,并将采样值xij储存在相应的采样寄存单元中;
所述采样值xij中i表示循环次数,i=1,2,...N,N表示设定的最大循环次数;
所述采样值xij中j表示采样次数,j=1,2,...M,M表示设定的最大采样次数;
第三步:多路选择单元依据顺时方式对采样寄存单元A、采样寄存单元B至采样寄存单元M中的采样值进行通道选通;
第四步:移位累加存取单元对接收到的采样寄存单元A中的采样值进行移位累加处理,得到第一个累加和;
移位累加存取单元对接收到的采样寄存单元B中的采样值进行移位累加处理,得到第二个累加和;
移位累加存取单元对接收到的采样寄存单元M中的采样值进行移位累加处理,得到第M个累加和;
第五步:累加和平均单元对接收的第一个累加和进行平均处理,得到第一个采样平均值;
累加和平均单元对接收的第二个累加和进行平均处理,得到第二个采样平均值;
累加和平均单元对接收的第M个累加和进行平均处理,得到第M个采样平均值;
第六步:将第一个采样平均值、第二个采样平均值、……、第M个采样平均值进行时间-电压的拟和,得到恢复后的串扰噪声。
本发明是一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,其优点在于:(一)能够有效的测量光电转换后的探测器信号受到同频周期信号干扰的情况;(二)使用检测电路(信号处理装置6)本身对干扰信号进行采样、累加,并恢复出原始干扰信号波形,不使用外部仪器(如锁相放大器),简单易行;(三)信号处理过程完全在中心处理器63内部进行,既可以恢复周期干扰信号的原始波形,也可以恢复周期干扰信号解调后的波形。
附图说明
图1是背景技术中干涉式光纤陀螺的结构框图。
图2是背景技术中干涉式光纤陀螺的信号处理装置的结构框图。
图3是本发明数字式逐点平均扫描原理示意图。
图4是本发明数字式逐点平均扫描法恢复调制串扰信号原始波形结构示意图。
图5A是实测含有周期串扰噪声信号的时间-电压关系曲线图。
图5B是采用本发明方法提取出的串扰噪声的时间-电压关系曲线图。
具体实施方式
下面将结合附图对本发明作进一步的详细说明。
为了测试干涉式光纤陀螺的信号处理装置(检测电路)中的调制串扰信号f2并恢复其原始波形,本发明提出的一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,此种测量方法可以有效地测量干涉式光纤陀螺内部的周期性串扰信号。因数字闭环光纤陀螺中信号处理装置所检测到的系统误差干扰源是与调制方波同频的周期信号(如图5A所示),故本发明测试方法采用数字式逐点平均扫描的方法来扫描调制串扰信号f2,然后对扫描获得的采样点信息进行累加平均处理,从而恢复其原始波形。
在本发明中,数字式逐点平均扫描方法的扫描是按一个周期内的调制方波的波形(如图3所示)进行的。
在本发明中,图5A中表示了实测含有调制串扰噪声的信号,信号幅值为1.5mV,信号为高斯白噪声分布,内部叠加了幅值为微伏量级的串扰噪声,频率和调制方波信号同频,调制串扰是数字闭环光纤陀螺的主要误差源,将引起陀螺输出产生死区、周期性噪声干扰、小角速度漂移、小角速度非线性等误差,严重降低陀螺精度。
为了实现对调制串扰信号f2依照图3所示的周期性波形进行“数字式逐点平均扫描”,本发明在中心处理器63选取的载体(芯片)上运用VerilogHDL编程设计进行模式配置得以实现。参见图4所示,在本发明中的中心处理器63包括有采样时序控制单元、移位累加存取单元、累加和平均单元、结果寄存单元、多路选择单元和采样寄存器;所述的采样寄存器由多个采样寄存单元构成。
在本发明中,采样时序控制单元依据晶振电路66提供的时钟信息t,对调制串扰信号f2采用数字式逐点平均扫描方式进行采样时间点t0的控制;
在本发明中,多路选择单元用于选择当前所需处理的采样寄存单元内的采样值xij
在本发明中,移位累加存取单元对接收到的采样值xij进行移位累加处理,得到累加和;
在本发明中,累加和平均单元对接收的累加和进行平均处理,得到采样平均值A(tj);
在本发明中,结果寄存单元用于对采样平均值进行储存。
在本发明中,若被测干扰信号的周期为T,在每个周期的起始处触发采样过程,每个周期内均匀采样M次,采样时间间隔为Δt,对于第j通道采样信号的运算过程为
Figure BSA00000255607600041
且j=1,2,…,M,式中,tj是最后一次对第j通道采样的时刻。按照通道采样信号的运算关系分别计算出各个j所对应的采样平均值A(tj),将采样平均值按时间顺序排列进行拟和(时间-电压)就能够恢复被测信号波形。
为了实现对干涉式光纤陀螺的信号处理装置中调制串扰误差干扰源的采集及复原,本发明采用下列的处理步骤:
第一步:采样时序控制单元控制时钟信息t来获得调制串扰信号f2进行数字式逐点平均扫描的采样时间点t0
第二步:在采样时间点t0上,对调制串扰信号f2在一个干扰信号的周期T内均匀采样M次,并将采样值xij(i=1,2,...N,j=1,2,...M)储存在相应的采样寄存单元中;
在本发明中,所述采样值xij中i表示循环次数,N表示设定的最大循环次数,j表示采样次数,M表示设定的最大采样次数。
参见图4所示,图中采样寄存单元A中储存j=1,且i=1,2,...N的采样值xi1;采样寄存单元2中储存j=2,且i=1,2,...N的采样值xi2;采样寄存单元M中储存j=M,且i=1,2,...N的采样值xiM;同理可得,任意一采样寄存单元中储存的采样值记为xij
第三步:多路选择单元依据顺时方式对采样寄存单元A、采样寄存单元B至采样寄存单元M中的采样值进行通道选通;
第四步:移位累加存取单元对接收到的采样寄存单元A中的采样值进行移位累加处理,得到第一个累加和;
移位累加存取单元对接收到的采样寄存单元B中的采样值进行移位累加处理,得到第二个累加和;
移位累加存取单元对接收到的采样寄存单元M中的采样值进行移位累加处理,得到第M个累加和;
第五步:累加和平均单元对接收的第一个累加和进行平均处理,得到第一个采样平均值;
累加和平均单元对接收的第二个累加和进行平均处理,得到第二个采样平均值;
累加和平均单元对接收的第M个累加和进行平均处理,得到第M个采样平均值;
第六步:将第一个采样平均值、第二个采样平均值、……、第M个采样平均值进行时间-电压的拟和,得到恢复后的串扰噪声(如图5B所示)。
在本发明中,图5B表示了使用数字式逐点平均扫描方式恢复得到的调制串扰噪声信号,该调制串扰信号f2和调制方波信号同频率,为167KHz,幅值为5μV,通过本发明的方法复原得到的串扰信号能够为后续的补偿提供数据支持,从而提高光纤陀螺的测量精度。
本发明测试方法采用数字式逐点平均扫描的方法恢复周期性的检测系统误差干扰信号的原始波形或解调后的波形。首先断开信号处理装置6与Y波导3的引线,信号处理装置6不向Y波导3施加调制信号,使测试波形中包含除了非理想方波调制信号以外的全部调制串扰信号。所述光功率信号f1经过前置滤波电路61放大,A/D转换器62转换成数字干扰信号f2后,进入中心处理器63中,中心处理器63内部的采样时序控制单元控制时钟,对调制串扰信号在每个周期内均匀采样后,储存在采样寄存器中,中心处理器63内部的多路选择单元对储存在采样寄存器中的采样值xij(i=1,2,...N,j=1,2,...M)进行选择并传送到信息融合计算,计算结果通过D/A转换器64转换成模拟信号后输出,即为恢复出的干扰原始波形。

Claims (2)

1.一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,所述检测电路包括有前放滤波电路(61)、A/D转换器(62)、中心处理器(63)、D/A转换器(64)、放大调理电路(65)和晶振电路(66);前放滤波电路(61)对光电探测器(5)输出的光功率信号f1进行滤波放大处理,然后滤波放大后的光功率信号经A/D转换器(62)后形成调制串扰信号f2输出给中心处理器(63);首先断开信号处理装置(6)与Y波导(3)的引线,信号处理装置(6)不向Y波导(3)施加调制信号,使测试波形中包含除了非理想方波调制信号以外的全部调制串扰信号;其特征在于包括有下列处理步骤:
所述中心处理器(63)包括有采样时序控制单元、移位累加存取单元、累加和平均单元、结果寄存单元、多路选择单元和采样寄存器;所述的采样寄存器由多个采样寄存单元构成;采样时序控制单元依据晶振电路(66)提供的时钟信息t,对调制串扰信号f2采用数字式逐点平均扫描方式进行采样时间点t0的控制;多路选择单元用于选择当前所需处理的采样寄存单元内的采样值xij;移位累加存取单元对接收到的采样值xij进行移位累加处理,得到累加和;累加和平均单元对接收的累加和进行平均处理,得到采样平均值A(tj);结果寄存单元用于对采样平均值进行储存;
第一步:采样时序控制单元控制时钟信息t来获得调制串扰信号f2进行数字式逐点平均扫描的采样时间点t0
第二步:在采样时间点t0上,对调制串扰信号f2在一个干扰信号的周期T内均匀采样M次,并将采样值xij储存在相应的采样寄存单元中;
所述采样值xij中i表示循环次数,i=1,2,...N,N表示设定的最大循环次数;
所述采样值xij中j表示采样次数,j=1,2,...M,M表示设定的最大采样次数;
第三步:多路选择单元依据顺时方式对第一采样寄存单元、第二采样寄存单元、......、第M采样寄存单元中的采样值进行通道选通;
第四步:移位累加存取单元对接收到的第一采样寄存单元中的采样值进行移位累加处理,得到第一个累加和;
移位累加存取单元对接收到的第二采样寄存单元中的采样值进行移位累加处理,得到第二个累加和;
移位累加存取单元对接收到的第M采样寄存单元中的采样值进行移位累加处理,得到第M个累加和; 
第五步:累加和平均单元对接收的第一个累加和进行平均处理,得到第一个采样平均值;
累加和平均单元对接收的第二个累加和进行平均处理,得到第二个采样平均值;
累加和平均单元对接收的第M个累加和进行平均处理,得到第M个采样平均值;
第六步:将第一个采样平均值、第二个采样平均值、......、第M个采样平均值进行时间-电压的拟合,得到恢复后的串扰噪声。
2.根据权利要求1所述的一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法,其特征在于:数字式逐点平均扫描是内嵌在中心处理器(63)选取的芯片上,并运用Verilog HDL编程语言实现。 
CN2010102715778A 2010-09-03 2010-09-03 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法 Expired - Fee Related CN101975584B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102715778A CN101975584B (zh) 2010-09-03 2010-09-03 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102715778A CN101975584B (zh) 2010-09-03 2010-09-03 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法

Publications (2)

Publication Number Publication Date
CN101975584A CN101975584A (zh) 2011-02-16
CN101975584B true CN101975584B (zh) 2012-05-30

Family

ID=43575485

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102715778A Expired - Fee Related CN101975584B (zh) 2010-09-03 2010-09-03 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法

Country Status (1)

Country Link
CN (1) CN101975584B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2482450C1 (ru) * 2011-12-16 2013-05-20 Федеральное государственное унитарное предприятие "Центр эксплуатации объектов наземной космической инфраструктуры" (ФГУП "ЦЭНКИ") Устройство тестирования электронного блока волоконно-оптического гироскопа
CN103401743B (zh) * 2013-07-02 2017-01-25 深圳市九洲电器有限公司 一种i2c总线消除干扰信号的方法和装置
CN109061230B (zh) * 2018-08-30 2022-06-14 衡阳市衡山科学城科技创新研究院有限公司 一种多轴检测电路板串扰特性的评估方法及系统
CN109387223B (zh) * 2018-09-27 2022-04-19 北京计算机技术及应用研究所 一种光纤陀螺信号检测电路的死区性能检测系统
CN109211274B (zh) * 2018-09-30 2020-07-14 北京控制工程研究所 一种光纤陀螺检测电路自主验证方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1118684C (zh) * 2001-06-13 2003-08-20 北京航空航天大学 开环光纤陀螺信号检测方法及装置
CN1664538A (zh) * 2004-03-01 2005-09-07 北京航空航天大学 光纤陀螺用集成光学调制器在线测试方法及其测试装置
CN101270991A (zh) * 2008-04-25 2008-09-24 北京航空航天大学 采用方波调制测量干涉式光纤陀螺本征频率和半波电压的系统
CN100498223C (zh) * 2006-09-30 2009-06-10 北京航空航天大学 具有信号差分放大的光纤陀螺前置放大和调制输出电路
CN100541127C (zh) * 2006-11-24 2009-09-16 北京航空航天大学 采用不对称方波调制方法测量干涉式光纤陀螺本征频率的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6212811A (ja) * 1985-07-10 1987-01-21 Japan Aviation Electronics Ind Ltd 光干渉角速度計
JPH09280871A (ja) * 1996-04-08 1997-10-31 Hitachi Cable Ltd 光ファイバジャイロ

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1118684C (zh) * 2001-06-13 2003-08-20 北京航空航天大学 开环光纤陀螺信号检测方法及装置
CN1664538A (zh) * 2004-03-01 2005-09-07 北京航空航天大学 光纤陀螺用集成光学调制器在线测试方法及其测试装置
CN100498223C (zh) * 2006-09-30 2009-06-10 北京航空航天大学 具有信号差分放大的光纤陀螺前置放大和调制输出电路
CN100541127C (zh) * 2006-11-24 2009-09-16 北京航空航天大学 采用不对称方波调制方法测量干涉式光纤陀螺本征频率的方法
CN101270991A (zh) * 2008-04-25 2008-09-24 北京航空航天大学 采用方波调制测量干涉式光纤陀螺本征频率和半波电压的系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JP昭62-12811A 1987.01.21
JP特开平9-280871A 1997.10.31
于海成等.利用DSP实现光纤陀螺数字信号处理中的多采样率转换.《测控技术》.2003,第22卷(第3期),全文. *
谭呈明等.数字闭环光纤陀螺复位电压控制技术研究.《中国惯性技术学报》.2005,第13卷(第1期),全文. *

Also Published As

Publication number Publication date
CN101975584A (zh) 2011-02-16

Similar Documents

Publication Publication Date Title
CN101270991B (zh) 采用方波调制测量干涉式光纤陀螺本征频率和半波电压的系统
CN107389097B (zh) 光纤陀螺Sagnac光纤环本征频率跟踪测量方法
CN111174773B (zh) 一种多相调制解调的光纤陀螺多闭环算法
CN101975584B (zh) 一种适用于干涉式光纤陀螺的检测电路系统误差的开环测试方法
CN104682187B (zh) 一种基于闭环反馈的拉曼激光系统相位噪声自动补偿装置及方法
CN100541127C (zh) 采用不对称方波调制方法测量干涉式光纤陀螺本征频率的方法
CN100458367C (zh) 一种对光纤陀螺2π电压进行自动跟踪的四态调制及解调方法
CN101709971B (zh) 一种抑制光纤陀螺振动误差的信号解调方法
CN112697125B (zh) 一种光纤陀螺本征频率的在线跟踪方法及系统
CN102650526A (zh) 一种基于相位比较调频连续波光纤陀螺的开环检测电路
CN111220142B (zh) 一种新型超高精度光纤陀螺仪的信号处理方法及装置
CN115077567B (zh) 一种基于波导复位误差的标度因数补偿系统及方法
CN202947744U (zh) 一种三轴一体化高精度光纤陀螺的光电控制系统
CN106871931B (zh) 一种闭环光纤陀螺温度补偿方法
CN104713575A (zh) 一种闭环光纤陀螺频率特性的测试方法
CN102183249B (zh) 一种光纤陀螺仪的萨格奈克相移跟踪方法
CN113310483B (zh) 一种数字闭环光纤陀螺本征频率实时跟踪装置及方法
CN113720321A (zh) 一种光纤陀螺温度补偿系统
CN114993281A (zh) 基于低相干光的全数字化闭环角速度检测系统及方法
Sun et al. Study on the Drift of Modulated Phase in Interference Fiber Optic Gyroscope.
CN109029412A (zh) 一种光纤陀螺闭环反馈回路非线性度测试方法及系统
CN101101213B (zh) 一种补偿数字闭环光纤陀螺标度因数非线性度的方法
CN103591946B (zh) 一种可消除尖峰信号的干涉式光纤陀螺
CN105674976A (zh) 光纤陀螺调制解调方法、提高标度因数稳定性方法及装置
CN204944505U (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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20110216

Assignee: General aviation (Beijing) Photoelectric Technology Co Ltd

Assignor: Beihang University

Contract record no.: 2013990000773

Denomination of invention: Open loop measuring method applicable to detection circuit system error of interference optical fiber gyroscope

Granted publication date: 20120530

License type: Exclusive License

Record date: 20131120

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
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: 20120530

Termination date: 20200903