CN107643090A - 自整角机/旋转变压器的数字角转换方法 - Google Patents

自整角机/旋转变压器的数字角转换方法 Download PDF

Info

Publication number
CN107643090A
CN107643090A CN201710786295.3A CN201710786295A CN107643090A CN 107643090 A CN107643090 A CN 107643090A CN 201710786295 A CN201710786295 A CN 201710786295A CN 107643090 A CN107643090 A CN 107643090A
Authority
CN
China
Prior art keywords
way
sin
signal
rotary transformer
voltage signal
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
CN201710786295.3A
Other languages
English (en)
Other versions
CN107643090B (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201710786295.3A priority Critical patent/CN107643090B/zh
Publication of CN107643090A publication Critical patent/CN107643090A/zh
Application granted granted Critical
Publication of CN107643090B publication Critical patent/CN107643090B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Transmission And Conversion Of Sensor Element Output (AREA)
  • Measurement Of Length, Angles, Or The Like Using Electric Or Magnetic Means (AREA)

Abstract

本发明公开了一种自整角机/旋转变压器的数字角转换方法,该方法首先采集自整角机或旋转变压器两路电压信号和一路激磁信号;然后利用两路电压信号与轴角的三角函数关系,通过激磁信号简化公式,得到只有轴角为未知变量的三角函数关系;最后通过求反三角函数,得到最终的数字化轴角。本发明具有转换过程简单、易于实现、精度高、实时性强、适用场合广泛、成本低、稳定可靠等优点。

Description

自整角机/旋转变压器的数字角转换方法
技术领域
本发明涉及数字角转换方法,特别是一种自整角机/旋转变压器的数字角转换方法。
背景技术
自整角机通常应用在自动控制系统中需要指示位置、角度的数值或者远距离调节执行机构的速度等场合,需要一根或多根轴随着另外的与其无机械连接的轴同步转动。旋转变压器是一种精密角度、位置、速度检测装置,适用于所有使用旋转编码器的场合,特别是高温、严寒、潮湿、高速、高震动等旋转编码器无法正常工作的场合。
自整角机/旋转变压器作为角度测量器件,在各种工业领域得到了广泛的应用,为适应工业现代化的需求,工业控制现场的网络化、智能化逐步发展,需要将自整角机/旋转变压器输出的相电压瞬时信号转化为数字量轴角。
美国ADI公司的自整角机/旋转变压器数字转换器AD2S44,集成了输入调节器、高速sin/cos乘法器、误差放大器、相敏调节器、积分器、压控振荡器、可逆计数器和输出锁存器,结构复杂,成本高昂。国内201所推出的ADI产品的原位替代品,如10、11、14ZSZ/XSZ系列,虽简化了ADI产品SDC/RDC 1702、1700、1704系列的部分结构,但精度也随之下降,成本仍然高昂。
申请号为201310186444.4的发明专利公开了一种自整角机/旋转变压器的数字转换器转换方法,该方法很大程度上解决了现有技术的成本问题,但该发明的转换方法多次进行近似,累加误差导致其精度较低,且使用单片机进行一系列的乘除和加减运算,实时性难以保证。
发明内容
本发明的目的在于提供一种自整角机/旋转变压器的数字角转换方法。
实现本发明目的的技术方案为:一种自整角机/旋转变压器的数字角转换方法,包括以下步骤:
步骤1,采集自整角机或旋转变压器两路电压信号和一路激磁信号;
步骤2,利用两路电压信号与轴角的三角函数关系,通过激磁信号简化公式,得到只有轴角为未知变量的三角函数关系;
步骤3,通过求反三角函数,得到最终的数字化轴角。
与现有技术相比,本发明的显著优点为:
本发明的自整角机/旋转变压器数字角转换方法具有转换过程简单、易于实现、精度高、实时性强、适用场合广泛、成本低、稳定可靠等优点。
附图说明
图1为本发明的自整角机/旋转变压器数字角转换方法流程图。
图2为本发明实施例的原理框架图。
具体实施方式
结合图1,一种自整角机/旋转变压器的数字角转换方法,包括以下步骤:
步骤1,采集自整角机或旋转变压器两路电压信号和一路激磁信号;
步骤2,利用两路电压信号与轴角的三角函数关系,通过激磁信号简化公式,得到只有轴角为未知变量的三角函数关系;
步骤3,通过求反三角函数,得到最终的数字化轴角。
进一步的,对于自整角机,数字角转换方法具体为:
自整角机的线电压瞬时值为
e1=E0sin(wt)sinθ
e2=E0sin(wt)sin(θ+120°)
e3=E0sin(wt)sin(θ+240°)
E0为激磁信号幅值,w为激磁信号角频率,θ为轴角;
模数转换器选择采集两路线电压信号及一路激磁信号,任取两路线电压信号,线电压信号取ei,ej,i=1,2,3,j=1,2,3,且i≠j;
设激磁信号Es=E0sin(ωt),采集后送至FPGA,分别进行两路加、减运算:
Ub=ej-ei=-3E0sin(ωt)cosθ
根据使用场合切换模式1和模式2,采用模式1时只有一路线电压信号及一路激磁信号参与运算,采用模式2时有两路线电压信号及一路激磁信号参与运算;模式1即反余弦法,可少采集处理一路电压信号,运算少,对应高实时性要求场合,模式2即反正切法,因两路电压信号都随激磁信号变化,处理过程趋一致,可减小转换对象通道的信号误差,运用量化极限逼近判决法,可进行更高精度的转换,对应高精度场合,另在一些无差别场合,可根据实际使用效果,进行切换。
两种模式的运算过程为:
模式1:计算-Ub/(3Es)=cosθ,由数学软件生成arccos的值,将值例化到FPGA的单口ROM中,通过查表1,再结合Ua、Es的符号,确定整周期轴角θ;
表1
其中,θ0为反三角函数得出值;
模式2:计算得到tanθ后,通过量化误差计算,再通过极限逼近判决,由查找表2求结合Ub、Es的符号,确定整周期轴角θ;
表2
对于旋转变压器,数字角转换方法具体为:
旋转变压器的两绕组电压的瞬时值为:
ea=E0sin(wt)sinθ
eb=E0sin(wt)cosθ
模数转换器采集两路瞬时电压信号及激磁信号,Ua=ea,Ub=eb,设激磁信号Es=E0sin(ωt),根据使用场合切换模式1和模式2,模式1为只有绕组电压信号及一路激磁信号参与运算,模式2有两路绕组电压信号及一路激磁信号参与运算;
模式1:计算Ub/Er=cosθ,得到cosθ,由数学软件生成arccos的值,将值例化到FPGA的单口ROM中,通过查表1,再结合Ua,Er的符号,确定整周期轴角θ。
模式2:计算Ua/Ub=tanθ,得到tanθ后,通过量化误差计算,再通过极限逼近判决,由查找表2求arctan(Ua/Ub),再结合Ub、Er的符号,确定整周期轴角θ。
模式2中通过量化误差,计算最小量化误差的对应端极值,将极限值归入极限判决,通过判决将极限值在最小量化误差允许内归入有限值。
进一步的,自整角机或旋转变压器的两路电压信号和一路激磁信号通过衰减网络衰减后经隔离放大器进行隔离。
下面结合具体实施例对本发明进行详细说明。
实施例1
一种自整角机/旋转变压器的数字角转换方法,自整角机/旋转变压器的两路线电压瞬时信号及激磁电压信号经过衰减网络,按一定衰减比进行衰减,根据输入信号的幅值范围,衰减比可调节,以适应隔离放大器输入范围。
自整角机/旋转变压器的激磁信号为
e=E0sin(ωt)
(1)对于自整角机信号
e1=E0sin(ωt)sinθ
e2=E0sin(ωt)sin(θ+120°)
e3=E0sin(ωt)sin(θ+240°)
e2,e3,e经过衰减及隔离放大器隔离后为:
E2=K1e2=K1E0sin(ωt)sin(θ+120°)
E3=K1e3=K1E0sin(ωt)sin(θ+240°)
Es=K1e=K1E0sin(ωt)
K1为自整角机信号经衰减网络后的衰减系数。
(2)对于旋转变压器信号
ea=E0sin(ωt)sinθ
eb=E0sin(ωt)cosθ
ea,eb,e经过衰减及隔离放大器隔离后为:
Ea=K2ea=K2E0sin(ωt)sinθ
Eb=K2eb=K2E0sin(ωt)cosθ
Er=K2e=K2E0sin(ωt)
K2为旋转变压器信号经衰减网络后的衰减系数。
上述公式中ω为激磁电压信号角频率,θ为自整角机/旋转变压器的轴角,E0为信号的最大幅值。
E2,E3,Es和Ea,Eb,Er经ADC采样数字化后,由FPGA进行如下处理:模式1
对于自整角机,对E2,E3进行运算处理:
-Ub/(3Es)=cosθ得:θ=arccos[-Ub/(3Es)]
对于旋转变压器
Ub/Er=cosθ得:θ=arccos(Ub/Er)
由FPGA内的单口ROM中查表得出θ,-Ub/(3Es)∈[-1,1],Ub/Er∈[-1,1],则在一个最小正周期中,θ∈[0,π],结合Ua,Er、Ua,Es符号即可得出整周期的θ值。
模式2
对于自整角机,对E2,E3进行运算处理:
得:
对于旋转变压器
Ua/Ub=tanθ得:
θ=arctan(Ua/Ub)
Ua/Ub∈(-∞,∞),单周期中,θ∈[-90°,90°],由于在θ→±90°时,tanθ→±∞,使用的查找表深度是有限的,故使用量化精度在两端极限进行逼近判决,结合Ub,Es、Ub,Er符号即可得出整周期的θ值。
使用量化误差极限逼近判决法,该方法的主要内容为:θ∈[-90°,90°],对θ的量化字长为nbit,小数部分字长为fbit,则量化结果为有符号定点数,表示为Fix_n_f,整数部分最大数为90,故应满足2n-f-1-1>90,-2n-f-1<-90,Fix_n_f的表示范围为[-2n-f-1.(1-2-f),(2n-f-1-1).(1-2-f)],量化误差为Δθ=2-f
计算tan(90°-Δθ)=M,tan(-90°+Δθ)=N,当tanθ>M时,可知θ的取值范围为90°-Δθ<θ<90°,量化误差Δθ在小数部分量化字长f一定大时,可将此时的θ取值判决为90°,同理当tanθ<N时,将此时的θ取值判决为-90°,其他数值直接查表即可。
实施例2
如图2所示,一种自整角机/旋转变压器的数字角转换方法,自整角机/旋转变压器的两路线电压瞬时信号及激磁电压信号经过衰减网络,按一定衰减比进行衰减,根据输入信号的幅值范围,衰减比可调节,以适应隔离放大器输入范围。
按一定衰减比衰减后的信号,经隔离放大器以固定增益隔离放大,经ADC采样量化后,送入FPGA中,通过输入信号类型的选择,以适应输入信号由自整角机和旋转变压器输入的不同所带来的差异。
上述的衰减网络选择电阻分压衰减,选用精度系数高、温度性能好的电阻,通过选择将不同的电阻值接入衰减网络,实现不同衰减比,隔离放大器可选用低成本高精度的ISO12X系列,本实施例选用具有高共模抑制比,高隔离电压,低非线性误差的ISO121,ADC选用多相电机控制系统专用模数转换芯片AD760X系列,本实施例选用高精度,高采样信噪比的AD7608,FPGA选用Spartan-6系列的XC6SLX150。
本实施例中自整角机为SGS-1,激磁电压115V,激磁频率400Hz,次级电压90V,为解释方便,旋转变压器为45XZ4-10,激磁电压115V,频率400Hz,变比1.0。
则为适应隔离放大器的输入范围,115V的衰减为取衰减比k1=1:16.5,90V的衰减为取k2=1:22.5,为方便论述,以下计算仍以k1,k2表示,隔离放大器的固定增益为1。
对于自整角机,经衰减网络和隔离放大器后,ADC的输入信号分别为:
对于旋转变压器,经衰减网络和隔离放大器后,ADC的输入信号分别为:
模式1
对于自整角机,对E2,E3进行运算处理:
得:
对于旋转变压器
Ub/Er=cosθ得:
θ=arccos(Ub/Er)
由FPGA内的单口ROM中查表得出θ,-23k1Ub/(54k2Es)∈[-1,1],Ub/Er∈[-1,1],则在一个最小正周期中,θ∈[0,π],结合Ua,Es、Ua,Er符号即可得出整周期的θ值。
模式2
对于自整角机,对E2,E3进行运算处理:
得:
对于旋转变压器
Ua/Ub=tanθ得:
θ=arctan(Ua/Ub)
Ua/Ub∈(-∞,∞),单周期中,θ∈[-90°,90°],由于在θ→±90°时,tanθ→±∞,使用的查找表深度是有限的,故使用量化精度在两端极限进行逼近判决,结合Ub,Es、Ub,Er符号即可得出整周期的θ值。
使用量化误差极限逼近判决法,该方法具体为:θ∈[-90°,90°],对θ的量化字长为18bit,小数部分字长为10bit,则量化结果为有符号定点数,表示为Fix_18_10,整数部分最大数为90,故满足218-10-1-1>90,-218-10-1<-90,Fix_18_10的表示范围为[-128.9990234375,127.9990234375],量化误差为Δθ=0.0009765625。
tan(90°-Δθ)≈58671,tan(-90°+Δθ)=-58671,为对应查找表地址,可通过移位和加法将正切值正值化,正值化后范围[0,131071],值小于6865时,θ=-90°,值大于124207时,θ=90°,为缩小查找表空间,中间值部分减去6865后按地址查表得θ。

Claims (3)

1.一种自整角机/旋转变压器的数字角转换方法,其特征在于,包括以下步骤:
步骤1,采集自整角机或旋转变压器两路电压信号和一路激磁信号;
步骤2,利用两路电压信号与轴角的三角函数关系,通过激磁信号简化公式,得到只有轴角为未知变量的三角函数关系;
步骤3,通过求反三角函数,得到最终的数字化轴角。
2.根据权利要求1所述的自整角机/旋转变压器的数字角转换方法,其特征在于,对于自整角机,数字角转换方法具体为:
自整角机的线电压瞬时值为
e1=E0sin(wt)sinθ
e2=E0sin(wt)sin(θ+120°)
e3=E0sin(wt)sin(θ+240°)
其中,E0为激磁信号幅值,w为激磁信号角频率,θ为轴角;
模数转换器选择采集两路线电压信号及一路激磁信号,任取两路线电压信号,线电压信号取ei,ej,i=1,2,3,j=1,2,3,且i≠j;
设激磁信号Es=E0sin(ωt),采集后送至FPGA,分别进行两路加、减运算:
<mrow> <msub> <mi>U</mi> <mi>a</mi> </msub> <mo>=</mo> <msub> <mi>e</mi> <mi>j</mi> </msub> <mo>+</mo> <msub> <mi>e</mi> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <msqrt> <mn>3</mn> </msqrt> <msub> <mi>E</mi> <mn>0</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mi>&amp;omega;</mi> <mi>t</mi> <mo>)</mo> </mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;theta;</mi> </mrow>
Ub=ej-ei=-3E0sin(ωt)cosθ
根据使用场合切换模式1和模式2,采用模式1时只有一路线电压信号及一路激磁信号参与运算,采用模式2时有两路线电压信号及一路激磁信号参与运算,两种模式的运算过程为:
模式1:计算-Ub/(3Es)=cosθ,由数学软件生成arccos的值,将值例化到FPGA的单口ROM中,通过查表1,再结合Ua、Es的符号,确定整周期轴角θ;
表1
其中,θ0为反三角函数得出值;
模式2:计算得到tanθ后,通过量化误差计算,再通过极限逼近判决,由查找表2求结合Ub、Es的符号,确定整周期轴角θ;
表2
对于旋转变压器,数字角转换方法具体为:
旋转变压器的两绕组电压的瞬时值为:
ea=E0sin(wt)sinθ
eb=E0sin(wt)cosθ
模数转换器采集两路瞬时电压信号及激磁信号,Ua=ea,Ub=eb,设激磁信号Es=E0sin(ωt),根据使用场合切换模式1和模式2,模式1为只有绕组电压信号及一路激磁信号参与运算,模式2有两路绕组电压信号及一路激磁信号参与运算;
模式1:计算Ub/Er=cosθ,得到cosθ,由数学软件生成arccos的值,将值例化到FPGA的单口ROM中,通过查表1,再结合Ua,Er的符号,确定整周期轴角θ。
模式2:计算Ua/Ub=tanθ,得到tanθ后,通过量化误差计算,再通过极限逼近判决,由查找表2求arctan(Ua/Ub),再结合Ub、Er的符号,确定整周期轴角θ。
3.根据权利要求2所述的自整角机/旋转变压器的数字角转换方法,其特征在于,自整角机或旋转变压器的两路电压信号和一路激磁信号通过衰减网络衰减后经隔离放大器进行隔离,然后进行数字角转换。
CN201710786295.3A 2017-09-04 2017-09-04 自整角机/旋转变压器的数字角转换方法 Expired - Fee Related CN107643090B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710786295.3A CN107643090B (zh) 2017-09-04 2017-09-04 自整角机/旋转变压器的数字角转换方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710786295.3A CN107643090B (zh) 2017-09-04 2017-09-04 自整角机/旋转变压器的数字角转换方法

Publications (2)

Publication Number Publication Date
CN107643090A true CN107643090A (zh) 2018-01-30
CN107643090B CN107643090B (zh) 2020-02-14

Family

ID=61110899

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710786295.3A Expired - Fee Related CN107643090B (zh) 2017-09-04 2017-09-04 自整角机/旋转变压器的数字角转换方法

Country Status (1)

Country Link
CN (1) CN107643090B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108507459A (zh) * 2018-03-12 2018-09-07 连云港杰瑞电子有限公司 一种自整角机信号角度的测量方法
CN109870618A (zh) * 2018-12-17 2019-06-11 天津新星科能源技术有限公司 一种多极旋转变压器检测设备及控制方法
CN111752310A (zh) * 2020-05-20 2020-10-09 哈船光电(武汉)有限公司 一种高精度轴角信号转换模块
CN113708700A (zh) * 2021-08-26 2021-11-26 重庆长安新能源汽车科技有限公司 纯电动汽车的电机旋变初始角自学习方法及标定方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202274885U (zh) * 2011-10-26 2012-06-13 闫勇 旋转变压器-数字转换模块
CN103023497A (zh) * 2012-12-10 2013-04-03 北京航天测控技术有限公司 自整角机模块的数字信号与模拟信号转换精度测试方法
CN103036572A (zh) * 2012-12-06 2013-04-10 上海无线电设备研究所 一种用于双通道旋转变压器的解码器及其解码方法
CN103336459A (zh) * 2013-05-08 2013-10-02 安徽师范大学 一种自整角机/旋转变压器的数字转换器转换方案
CN204515404U (zh) * 2015-01-21 2015-07-29 安徽师范大学 数字转换器

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN202274885U (zh) * 2011-10-26 2012-06-13 闫勇 旋转变压器-数字转换模块
CN103036572A (zh) * 2012-12-06 2013-04-10 上海无线电设备研究所 一种用于双通道旋转变压器的解码器及其解码方法
CN103023497A (zh) * 2012-12-10 2013-04-03 北京航天测控技术有限公司 自整角机模块的数字信号与模拟信号转换精度测试方法
CN103336459A (zh) * 2013-05-08 2013-10-02 安徽师范大学 一种自整角机/旋转变压器的数字转换器转换方案
CN204515404U (zh) * 2015-01-21 2015-07-29 安徽师范大学 数字转换器

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵华: "一种自整角机信号数字化转换方法研究与实现", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108507459A (zh) * 2018-03-12 2018-09-07 连云港杰瑞电子有限公司 一种自整角机信号角度的测量方法
CN108507459B (zh) * 2018-03-12 2020-06-16 连云港杰瑞电子有限公司 一种自整角机信号角度的测量方法
CN109870618A (zh) * 2018-12-17 2019-06-11 天津新星科能源技术有限公司 一种多极旋转变压器检测设备及控制方法
CN111752310A (zh) * 2020-05-20 2020-10-09 哈船光电(武汉)有限公司 一种高精度轴角信号转换模块
CN113708700A (zh) * 2021-08-26 2021-11-26 重庆长安新能源汽车科技有限公司 纯电动汽车的电机旋变初始角自学习方法及标定方法
CN113708700B (zh) * 2021-08-26 2023-08-01 深蓝汽车科技有限公司 纯电动汽车的电机旋变初始角自学习方法及标定方法

Also Published As

Publication number Publication date
CN107643090B (zh) 2020-02-14

Similar Documents

Publication Publication Date Title
CN107643090A (zh) 自整角机/旋转变压器的数字角转换方法
CN101521480B (zh) 一种旋转变压器信号解算方法及解算器
CN104914304B (zh) 一种波动负荷下的电能精确计量方法及计量装置
CN102751991A (zh) 基于δς调制原理的旋转变压器数字转换装置及方法
CN103557830B (zh) 一种兼容多种传感器的测角板卡
CN101281041A (zh) 单电源cmos集成电路自整角机/旋转变压器-数字转换技术
CN110865238B (zh) 一种基于准谐波模型采样算法的交流电阻测量方法及装置
CN101281042A (zh) 高精度cmos集成电路自整角机/旋转变压器-数字转换技术
CN104122439A (zh) 提高相位校正精度的电能表
CN102879017B (zh) 双速轴角数字转换器粗精组合系统
CN104330096B (zh) 一种测量信号的校正补偿和自动标定的方法及装置
CN100462686C (zh) 一种自整角机/旋转变压器-模拟直流电压转换方法
US20020033694A1 (en) Method to determine fundamental and harmonic oscillations of a measured electrical quantity
CN103575976B (zh) 一种纯90度移相式无功功率测量方法
CN105227185A (zh) 一种高精度模拟旋转变压器输出电路
CN210982602U (zh) 一种具有谐波测量的高准确度三相标准电能表
CN104135284B (zh) 一种鉴相方法、装置,以及锁相方法和锁相环
CN103336459B (zh) 一种自整角机/旋转变压器的数字转换器转换方案
CN202274885U (zh) 旋转变压器-数字转换模块
CN115825546A (zh) 一种直流暂态电流宽频数字标准器实现方法
CN109581044B (zh) 应用于旋转变压器的电压检测装置及工作方法
CN103575975A (zh) 一种纯90度移相式无功功率测量电路
CN103344826A (zh) 非线性负载电能计量装置
CN108037731A (zh) 一种相位积分运算变换的频差干涉信号高分辨细分系统
CN210246636U (zh) 电机检测系统和应用该检测系统的电动汽车

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for 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: 20200214

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