CN104601080A - 定子磁链电压模型的偏移计算和反馈补偿方法 - Google Patents

定子磁链电压模型的偏移计算和反馈补偿方法 Download PDF

Info

Publication number
CN104601080A
CN104601080A CN201510014465.7A CN201510014465A CN104601080A CN 104601080 A CN104601080 A CN 104601080A CN 201510014465 A CN201510014465 A CN 201510014465A CN 104601080 A CN104601080 A CN 104601080A
Authority
CN
China
Prior art keywords
beta
alpha
magnetic linkage
psi
output
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
CN201510014465.7A
Other languages
English (en)
Other versions
CN104601080B (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.)
Xian Polytechnic University
Original Assignee
Xian Polytechnic 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 Xian Polytechnic University filed Critical Xian Polytechnic University
Priority to CN201510014465.7A priority Critical patent/CN104601080B/zh
Publication of CN104601080A publication Critical patent/CN104601080A/zh
Application granted granted Critical
Publication of CN104601080B publication Critical patent/CN104601080B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明公开了一种定子磁链电压模型的偏移计算和反馈补偿方法,实时检测积分器的输出偏移量,对输出偏移量进行滤波得到平滑的偏移量,把平滑的偏移量从输出中减掉,对平滑的偏移量设计调节器,将调节器输出反馈到输入端与输入偏移量相抵消,此时积分器的输出就是实际磁链,积分器的输入就是实际的定子反电势。本发明定子磁链电压模型的偏移计算和反馈补偿方法,相对于使用惯性或其他环节替代纯积分的磁链计算而言,该方法不会造成磁链幅值和相位的畸变,磁链检测的结果无偏移,提高了磁链辨识的精度,减少了磁链幅值波动量,直接从输出中减去偏移量,大大缩短了检测时间。

Description

定子磁链电压模型的偏移计算和反馈补偿方法
技术领域
本发明属于交流电机变频调速方法技术领域,具体涉及一种定子磁链电压模型的偏移计算和反馈补偿方法。
背景技术
电机磁链的在线估计是实现磁场定向控制技术、直接转矩控制技术的关键。
在基于电压模型定子磁链辨识算法中,由于纯积分环节的存在,在辨识的定子磁链中引入了直流偏置误差或初始相位积分误差,从而影响了磁链辨识的精度。为了解决纯积分环节存在的问题,一种常用的方法是用一阶低通滤波器来代替纯积分环节。一阶低通滤波器虽然会减小直流偏置误差及消除初始相位积分误差,但它又会在辨识的磁链中引入幅值误差和相位超前,尤其当电机运行在低通截止频率以下时,使电机控制精度很低。
论文《定子磁链电压模型的偏移计算和反馈补偿新方法》(陈增禄,马鑫姬弘扬,刘娟娟,汪冰.定子磁链电压模型的偏移计算和反馈补偿新方法[J].中国电机工程学报.2014(6)),提出了计算磁链输出偏移量的新方法,并通过反馈调节来抑制并消除这个偏移。该方法虽然能得到较为满意的磁链检测结果,但是检测结果存在波动;另外该方法需要较长的检测时间,不利于电机的动态控制。
发明内容
本发明的目的是提供一种定子磁链电压模型的偏移计算和反馈补偿方法,磁链检测结果无偏移,减少了磁链幅值波动量,大大缩短检测时间。
本发明所采用的技术方案是:定子磁链电压模型的偏移计算和反馈补偿方法,实时检测积分器的输出偏移量,对输出偏移量进行滤波得到平滑的偏移量,把平滑的偏移量从输出中减掉,对平滑的偏移量设计调节器,将调节器输出反馈到输入端与输入偏移量相抵消,此时积分器的输出就是实际磁链,积分器的输入就是实际的定子反电势。
具体按照以下步骤实施:
步骤1:在α-β坐标下建立包含偏移量的输入为:
e sα * ( t ) = e sα ( t ) + q α e sβ * ( t ) = e sβ ( t ) + q β - - - ( 1 )
式中,e、e表示定子反电势;e* 、e* 表示含有偏移量qα、qβ的定子反电势;下标α和β分别表示α和β轴分量;
包含偏移量的输出为:
ψ sα * ( t ) = ψ sα ( t ) + d α ψ sβ * ( t ) = ψ sβ ( t ) + d β - - - ( 2 )
式中,ψ、ψ表示磁链;ψ* 、ψ* 表示含有偏移量dα、dβ的磁链;
设上一次定子反电势采样值为esα1、esβ1,下一次定子反电势采样值为e 2、esβ2;上一次定子磁链采样值为ψ* sα1、ψ* sβ1,下一次定子磁链采样值为ψ* sα2、ψ* sβ2;则输出偏移量dα和dβ为:
d α = - e sα 1 e sβ 2 ψ sα 1 * + e sα 2 e sβ 1 ψ sα 2 * - e sβ 1 e sβ 2 ( ψ sβ 1 * - ψ sβ 2 * ) e sα 2 e sβ 1 - e sα 1 e sβ 2 d β = - e sα 2 e sβ 1 ψ sβ 1 * + e sα 1 e sβ 2 ψ sβ 2 * - e sα 1 e sα 2 ( ψ sα 1 * - ψ sα 2 * ) e sα 1 e sβ 2 - e sα 2 e sβ 1 - - - ( 3 )
步骤2:由于步骤1计算得到的偏移量是不连续变化的,直接减掉会造成磁链幅值的突变,因此加一个惯性滤波环节,表达式为
G 1 ( s ) = 1 T 1 s + 1 - - - ( 4 )
式(4)是对该惯性滤波环节的一种复数域描述,其中T1表示该惯性滤波的时间常数、其单位是秒,s表示微分运算的符号;
用式(4)对步骤1得到的输出偏移量dα和dβ进行惯性滤波,从而得到平滑的偏移量dα 和dβ
dα =dα*G1(s),dβ =dβ*G1(s)    (5)
步骤3:对步骤2得到的平滑偏移量dα 和dβ 设计调节器,表达式为
G 2 ( s ) = H ( T 2 s + 1 ) s - - - ( 6 )
式(6)是对该调节器的一种复数域描述,其中H表示该调节器的放大倍数,T2表示该调节器的时间常数、其单位是秒,s表示微分运算的符号;
则调节器输出qα 和qβ 为:
qα =dα *G2(s),qβ =dβ *G2(s)   (7)
步骤4:将步骤3得到的qα 和qβ 反馈到输入端,从输入e* 、e* 中减掉,从而抵消输入中的偏移量qα、qβ,得到纯净的定子反电势e、e
步骤5:将步骤2得到的平滑偏移量dα 和dβ ,直接从输出ψ* 、ψ* 中减掉,得到不含偏移量的磁链ψ、ψ,即ψ=ψ* -dα ,ψ=ψ* -dβ ,缩短了***调节时间。
本发明的特点还在于,
步骤1中的磁链ψ、ψ是采用一个纯积分环节对定子反电势e、e积分得到。
本发明的有益效果是:本发明定子磁链电压模型的偏移计算和反馈补偿方法,从定子反电势和磁链的相位关系出发,计算α-β轴磁链输出偏移量的公式并对其进行惯性滤波,磁链幅值平滑,减少了磁链幅值波动量,设计了偏移量反馈调节的调节器,使得闭环调节后偏移量趋于0,最终使得磁链检测的结果无偏移,相对于使用惯性或其他环节替代纯积分的磁链计算而言,该方法不会造成磁链幅值和相位的畸变,提高了磁链辨识的精度,直接从输出中减去偏移量,大大缩短了检测时间。
附图说明
图1是本发明定子磁链电压模型的偏移计算和反馈补偿方法的原理框图;
图2是采用论文《定子磁链电压模型的偏移计算和反馈补偿新方法》得到的仿真图;
图3是采用本发明定子磁链电压模型的偏移计算和反馈补偿方法得到的仿真图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明定子磁链电压模型的偏移计算和反馈补偿方法,具体按照以下步骤实施:
步骤1:在α-β坐标下建立包含偏移量的输入为:
e sα * ( t ) = e sα ( t ) + q α e sβ * ( t ) = e sβ ( t ) + q β - - - ( 1 )
式中,如图1所示,e、e表示定子反电势;e* 、e* 表示含有偏移量qα、qβ的定子反电势;下标α和β分别表示α和β轴分量;
包含偏移量的输出为:
ψ sα * ( t ) = ψ sα ( t ) + d α ψ sβ * ( t ) = ψ sβ ( t ) + d β - - - ( 2 )
式中,如图1所示,ψ、ψ表示磁链;ψ* 、ψ* 表示含有偏移量dα、dβ的磁链;其中磁链ψ、ψ是采用一个纯积分环节对定子反电势e、e积分得到。
设上一次定子反电势采样值为esα1、esβ1,下一次定子反电势采样值为e 2、esβ2;上一次定子磁链采样值为ψ* sα1、ψ* sβ1,下一次定子磁链采样值为ψ* sα2、ψ* sβ2;如图1所示,经过偏移量计算模块,计算出输出偏移量dα和dβ为:
d α = - e sα 1 e sβ 2 ψ sα 1 * + e sα 2 e sβ 1 ψ sα 2 * - e sβ 1 e sβ 2 ( ψ sβ 1 * - ψ sβ 2 * ) e sα 2 e sβ 1 - e sα 1 e sβ 2 d β = - e sα 2 e sβ 1 ψ sβ 1 * + e sα 1 e sβ 2 ψ sβ 2 * - e sα 1 e sα 2 ( ψ sα 1 * - ψ sα 2 * ) e sα 1 e sβ 2 - e sα 2 e sβ 1 - - - ( 3 )
步骤2:由于步骤1计算得到的偏移量是不连续变化的,直接减掉会造成磁链幅值的突变,因此加一个如图1所示惯性滤波环节G1(s),其表达式为:
G 1 ( s ) = 1 T 1 s + 1 - - - ( 4 )
式(4)是对该惯性滤波环节的一种复数域描述,其中T1表示该惯性滤波的时间常数、其单位是秒,s表示微分运算的符号;
用式(4)对步骤1得到的输出偏移量dα和dβ进行惯性滤波,从而得到平滑的偏移量dα 和dβ
dα =dα*G1(s),dβ =dβ*G1(s)      (5)
步骤3:对步骤2得到的平滑偏移量dα 和dβ 设计调节器,如图1所示的G2(s),其表达式为:
G 2 ( s ) = H ( T 2 s + 1 ) s - - - ( 6 )
式(6)是对该调节器的一种复数域描述,其中H表示该调节器的放大倍数,T2表示该调节器的时间常数、其单位是秒,s表示微分运算的符号;
则调节器输出qα 和qβ 为:
qα =dα *G2(s),qβ =dβ *G2(s)    (7)
步骤4:将步骤3得到的qα 和qβ 反馈到输入端,从输入e* 、e* 中减掉,从而抵消输入中的偏移量qα、qβ,得到纯净的定子反电势e、e
步骤5:将步骤2得到的平滑偏移量dα 和dβ ,直接从输出ψ* 、ψ* 中减掉,得到不含偏移量的磁链ψ、ψ,即ψ=ψ* -dα ,ψ=ψ* -dβ ,缩短了***调节时间。
本发明的原理为:采用一个纯积分环节对定子反电势es积分以得到定子磁链ψs。输入反电势中除了定子反电势es外还包含输入偏移量q;积分器输出磁链中除了定子磁链ψs以外,还包含输出偏移量d。本发明实时检测出积分器输出中的偏移量d;然后将偏移量d通过调节器反馈到纯积分环节的输入端,***动态平衡后调节器的输出就等于q,恰好与输入偏移量q相抵消。此时积分器的输出就是实际磁链ψs,积分器的输入就是实际的定子反电势es
图2是采用论文《定子磁链电压模型的偏移计算和反馈补偿新方法》得到的仿真图,横轴是时间,纵轴是幅值,图中波形分别是磁链ψ、ψ。从图2可以看出,磁链ψ、ψ在20秒左右达到稳态,波动量约为110%。图3是采用本发明定子磁链电压模型的偏移计算和反馈补偿方法得到的仿真图,横轴是时间,纵轴是幅值,图中波形分别是磁链ψ、ψ。从图3可以看出,磁链ψ、ψ在2秒左右达到稳态,检测时间比论文中减小了90%;波动量约为20%,比论文中减小了90%。
本发明的有益效果在于:
(1)对偏移量进行了滤波,磁链幅值平滑,与论文《定子磁链电压模型的偏移计算和反馈补偿新方法》中的方法相比,磁链幅值波动量减少了90%;
(2)直接从输出中减去偏移量,极大地缩短了检测时间,与论文《定子磁链电压模型的偏移计算和反馈补偿新方法》中的方法相比,检测时间缩短了90%;
(3)相对于使用惯性或其他环节替代纯积分的磁链计算而言,该方法不会造成磁链幅值和相位的畸变;
(4)该方法不需要高阶滤波环节;
(5)该方法无论在稳态还是动态过程中,都能得到满意的磁链检测结果。
实施例
假设某一时刻定子反电势为:
e=sinθ,e=cosθ
输入的偏移量为10,即
qα=10,qβ=10
根据式(1),含偏移量的输入为:
e sα * ( t ) = sin θ + 10 e sβ * ( t ) = cos θ + 10
***动态平衡以后,由式(3)计算得到输出的偏移量dα和dβ
根据步骤2,取T1=0.02,则惯性滤波环节为:
G 1 ( s ) = 1 0.02 s + 1
对输出偏移量dα和dβ进行惯性滤波,得到平滑的偏移量dα 和dβ
根据步骤3,取T2=3,H=0.2,对平滑的偏移量dα 和dβ 设计调节器:
G 2 ( s ) = 0.2 ( 3 s + 1 ) s
***动态平衡后,调节器的输出为:
qα′=10,qβ′=10
将得到的qα 和qβ 反馈到输入端,从输入e* 、e* 中减掉,得到纯净的定子反电势为:
e=sinθ,e=cosθ
将得到的平滑偏移量dα 和dβ ,直接从输出ψ* 、ψ* 中减掉,得到不含偏移量的磁链:ψ=ψ* -10,ψ=ψ* -10。

Claims (3)

1.定子磁链电压模型的偏移计算和反馈补偿方法,其特征在于,实时检测积分器的输出偏移量,对输出偏移量进行滤波得到平滑的偏移量,把平滑的偏移量从输出中减掉,对平滑的偏移量设计调节器,将调节器输出反馈到输入端与输入偏移量相抵消,此时积分器的输出就是实际磁链,积分器的输入就是实际的定子反电势。
2.如权利要求1所述的定子磁链电压模型的偏移计算和反馈补偿方法,其特征在于,具体按照以下步骤实施:
步骤1:在α-β坐标下建立包含偏移量的输入为:
e sα * ( t ) = e sα ( t ) + q α e sβ * ( t ) = e sβ ( t ) + q β - - - ( 1 )
式中,e、e表示定子反电势;e* 、e* 表示含有偏移量qα、qβ的定子反电势;下标α和β分别表示α和β轴分量;
包含偏移量的输出为:
ψ sα * ( t ) = ψ sα ( t ) + d α ψ sβ * ( t ) = ψ sβ ( t ) + d β - - - ( 2 )
式中,ψ、ψ表示磁链;ψ* 、ψ* 表示含有偏移量dα、dβ的磁链;
设上一次定子反电势采样值为esα1、esβ1,下一次定子反电势采样值为e 2、esβ2;上一次定子磁链采样值为ψ* sα1、ψ* sβ1,下一次定子磁链采样值为ψ* sα2、ψ* sβ2;则输出偏移量dα和dβ为:
d α = - e sα 1 e sβ 2 ψ sα 1 * + e sα 2 e sβ 1 ψ sα 2 * - e sβ 1 e sβ 2 ( ψ sβ 1 * - ψ sβ 2 * ) e sα 2 e sβ 1 - e sα 1 e sβ 2 d β = - e sα 2 e sβ 1 ψ sβ 1 * + e sα 1 e sβ 2 ψ sβ 2 * - e sα 1 e sα 2 ( ψ sα 1 * - ψ sα 2 * ) e sα 1 e sβ 2 - e sα 2 e sβ 1 - - - ( 3 )
步骤2:由于步骤1计算得到的偏移量是不连续变化的,直接减掉会造成磁链幅值的突变,因此加一个惯性滤波环节,表达式为
G 1 ( s ) = 1 T 1 s + 1 - - - ( 4 )
式(4)是对该惯性滤波环节的一种复数域描述,其中T1表示该惯性滤波的时间常数、其单位是秒,s表示微分运算的符号;
用式(4)对步骤1得到的输出偏移量dα和dβ进行惯性滤波,从而得到平滑的偏移量d'α和d'β
d'α=dα*G1(s),d'β=dβ*G1(s)    (5)
步骤3:对步骤2得到的平滑偏移量d'α和d'β设计调节器,表达式为
G 2 ( s ) = H ( T 2 s + 1 ) s - - - ( 6 )
式(6)是对该调节器的一种复数域描述,其中H表示该调节器的放大倍数,T2表示该调节器的时间常数、其单位是秒,s表示微分运算的符号;
则调节器输出q'α和q'β为:
q'α=d'α*G2(s),q'β=d'β*G2(s)    (7)
步骤4:将步骤3得到的q'α和q'β反馈到输入端,从输入e* 、e* 中减掉,从而抵消输入中的偏移量qα、qβ,得到纯净的定子反电势e、e
步骤5:将步骤2得到的平滑偏移量d'α和d'β,直接从输出ψ* 、ψ* 中减掉,得到不含偏移量的磁链ψ、ψ,即ψ=ψ* -d'α,ψ=ψ* -d'β,缩短了***调节时间。
3.如权利要求1所述的定子磁链电压模型的偏移计算和反馈补偿方法,其特征在于,所述步骤1中的磁链ψ、ψ是采用一个纯积分环节对定子反电势e、e积分得到。
CN201510014465.7A 2015-01-12 2015-01-12 定子磁链电压模型的偏移计算和反馈补偿方法 Expired - Fee Related CN104601080B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510014465.7A CN104601080B (zh) 2015-01-12 2015-01-12 定子磁链电压模型的偏移计算和反馈补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510014465.7A CN104601080B (zh) 2015-01-12 2015-01-12 定子磁链电压模型的偏移计算和反馈补偿方法

Publications (2)

Publication Number Publication Date
CN104601080A true CN104601080A (zh) 2015-05-06
CN104601080B CN104601080B (zh) 2017-01-11

Family

ID=53126660

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510014465.7A Expired - Fee Related CN104601080B (zh) 2015-01-12 2015-01-12 定子磁链电压模型的偏移计算和反馈补偿方法

Country Status (1)

Country Link
CN (1) CN104601080B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105162425A (zh) * 2015-08-11 2015-12-16 成都思邦力克科技有限公司 差分积分器采集终端机
CN105187018A (zh) * 2015-08-11 2015-12-23 成都思邦力克科技有限公司 差分积分器采集装置的改进结构
CN107154764A (zh) * 2017-04-11 2017-09-12 卧龙电气集团股份有限公司 一种基于积分降幅和正交误差补偿模块的电机磁链观测器

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101795109B (zh) * 2010-02-11 2012-05-30 中国船舶重工集团公司第七一二研究所 一种异步电机矢量控制器
CN102624321A (zh) * 2012-03-20 2012-08-01 浙江大学 一种双馈感应风力发电机的电磁转矩消抖方法
CN103414424B (zh) * 2013-06-19 2015-12-23 山西潞安环保能源开发股份有限公司 交流电机定子磁链估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JAEBOK LEE ET AL.: "A New Integration Method to Estimate stator flux in induction motor with dc offset error", 《TENCON SPRING CONFERENCE,2013 IEEE》 *
WEI XU ET AL.: "Reduced Parameter Sensitivity Stator Flux Linkage Observer in Deadbeat-Direct Torque and Flux Control for IPMSMs", 《IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS》 *
刘忠举等: "基于改进型积分器的矢量控制***研究", 《电力电子技术》 *
陈增禄等: "定子磁链电压模型的偏移计算和反馈补偿新方法", 《中国电机工程学报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105162425A (zh) * 2015-08-11 2015-12-16 成都思邦力克科技有限公司 差分积分器采集终端机
CN105187018A (zh) * 2015-08-11 2015-12-23 成都思邦力克科技有限公司 差分积分器采集装置的改进结构
CN105187018B (zh) * 2015-08-11 2018-01-30 成都思邦力克科技有限公司 差分积分器采集装置的改进结构
CN105162425B (zh) * 2015-08-11 2018-01-30 成都思邦力克科技有限公司 差分积分器采集终端机
CN107154764A (zh) * 2017-04-11 2017-09-12 卧龙电气集团股份有限公司 一种基于积分降幅和正交误差补偿模块的电机磁链观测器
CN107154764B (zh) * 2017-04-11 2019-07-30 卧龙电气集团股份有限公司 一种基于积分降幅和正交误差补偿模块的电机磁链观测器

Also Published As

Publication number Publication date
CN104601080B (zh) 2017-01-11

Similar Documents

Publication Publication Date Title
CN103731082B (zh) 一种基于直接转矩控制的永磁同步电机定子磁链估计方法
CN103872951B (zh) 基于滑模磁链观测器的永磁同步电机转矩控制方法
CN104953915B (zh) 一种基于趋近律的永磁同步电机滑模控制策略
CN107482977B (zh) 一种永磁同步电机转子位置和转速检测方法
CN108069021B (zh) 一种舵机及其控制***
CN105680754B (zh) 一种永磁同步电机的直交轴电流矢量复合控制器
CN104935235B (zh) 电机的电流谐波抑制方法
CN105846748B (zh) 一种基于矢量变换与信号滤波的定子磁链计算方法
CN103414424B (zh) 交流电机定子磁链估计方法
CN103516284B (zh) 一种永磁同步电机电流增量预测算法
CN110323983B (zh) 一种永磁同步电机的电流解耦方法、装置、设备及介质
CN104601080A (zh) 定子磁链电压模型的偏移计算和反馈补偿方法
CN103344815B (zh) 一种宽范围变化的电气参数测量方法及***
CN108964556A (zh) 用于驱动永磁式同步电机的无速度传感器控制装置
CN112910329A (zh) 一种基于移动均值滤波器的pmsm无位置传感器控制方法
CN110620533A (zh) 一种表贴式永磁同步电机无传感器控制方法
CN103116281A (zh) 轴向混合磁轴承无模型自适应控制***及其控制方法
Zhang et al. Robust plug-in repetitive control for speed smoothness of cascaded-PI PMSM drive
CN108649851B (zh) 一种永磁同步电机最大转矩电流比控制方法
CN105227022A (zh) 基于改进无功功率模型的异步电机转子时间常数在线辨识方法
CN106452256A (zh) 基于转子磁链观测器的异步电机参数在线矫正方法
CN103607158A (zh) 低速下基于pi调节器和锁相环原理的磁链估计器
CN101132163A (zh) 具有补偿能力的dtc磁链观测滤波方法
CN108649849B (zh) 一种简单的无传感器永磁同步电机速度估测方法
CN114024382B (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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170111

Termination date: 20180112

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