CN109507871B - 用于两轮平衡车车身平衡控制的pid参数整定方法及产品 - Google Patents

用于两轮平衡车车身平衡控制的pid参数整定方法及产品 Download PDF

Info

Publication number
CN109507871B
CN109507871B CN201811513257.1A CN201811513257A CN109507871B CN 109507871 B CN109507871 B CN 109507871B CN 201811513257 A CN201811513257 A CN 201811513257A CN 109507871 B CN109507871 B CN 109507871B
Authority
CN
China
Prior art keywords
pid parameter
real
inclination angle
pid
current
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
CN201811513257.1A
Other languages
English (en)
Other versions
CN109507871A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of 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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201811513257.1A priority Critical patent/CN109507871B/zh
Publication of CN109507871A publication Critical patent/CN109507871A/zh
Application granted granted Critical
Publication of CN109507871B publication Critical patent/CN109507871B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Motorcycle And Bicycle Frame (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定方法及产品,该PID参数整定方法,能够根据两轮平衡车的在输入PID参数后的输出值,对PID参数样本进行判定,而且利用拟定输入、输出之间的关系,利用理论PID参数与实际输出之间的关系,对最终整定的PID参数进行确定,在实施方程中,不需要人为操作,能够方便、快捷地对两轮平衡车车身平衡控制的PID参数进行整定。本发明提供的一种用于两轮平衡车车身平衡控制的PID参数整定系统、两轮平衡车、计算机设备、计算机可读存储介质也具有上述有益效果,在此不再赘述。

Description

用于两轮平衡车车身平衡控制的PID参数整定方法及产品
技术领域
本发明涉及自动控制技术领域,特别涉及一种用于两轮平衡车车身平衡控制的PID参数整定方法及产品。
背景技术
现如今,随着现代科技的高速发展,人们的出行活动变得越来越频繁,对出行工具的需求也越来越高。在众多的出行工具中,电动平衡车,这种代步工具的出现,凭借着体形小便于携带、绿色环保的独特优势,成功的吸引了众人的眼球,迎合了大量80、90后对现代出行方式的需求,形成较好的市场反应。
在对平衡车控制系统的设计过程中,比较关键的就是对两轮平衡车车身平衡的控制。一个好的车身平衡系统不仅能够提高用户体验,还能提高用户在使用过程中的安全系数,有效减少意外的发生。
目前,市面上比较广泛的是根据车身的姿态传感器反馈的车身倾角信息作为反馈信号,采用PID闭环控制的方法来设计控制平衡车的车身平衡控制系统。PID控制系统包含了三个控制环节,分别为比例环节、积分环节、微分环节。
在设计PID控制系统过程中,最关键的就是对PID的比例、积分、微分三个控制环节的权重比参数进行调节,简称PID参数整定。目前比较广泛PID参数整过程定是基于在调试过程中根据控制系统的输出波形变化情况来进行调整。这就要求调试者需要有一定的调试经验。除了上述方法,还有一种PID参数整定方法叫做Ziegler-Nichols整定法,它的主要方法是仅在比例作用下将系统调节到临界震荡状态,得到临界增益和震荡周期,根据这两个量进行相应的计算便得到相对合适的PID参数。上述Ziegler-Nichols整定法被广泛运用到PID调节中,但也存在临界状态难以把握的缺点。同时,该方法具有一般性,不能对不同控制系统在实际应用中的差异性作出弥补。
因此,如何提供一种用于两轮平衡车车身平衡控制的PID整定方案,能够方便、快捷地对两轮平衡车车身平衡控制的PID参数进行整定,是本领域技术人员亟待解决的技术问题。
发明内容
有鉴于此,本发明的目的在于提供一种用于两轮平衡车车身平衡控制的PID参数整定方法及产品,能够方便、快捷地对两轮平衡车车身平衡控制的PID参数进行整定。其具体方案如下:
第一方面,本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定方法,包括:
A1:获取N组初始PID参数样本,以择一作为当前PID参数;
A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
A3:判断实时倾角值是否小于预设倾角阈值;
A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
A6:利用当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
A8:判断实时倾角值是否小于预设倾角阈值;
A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。
优选地,
所述获取所述初始PID参数样本输入两轮平衡车后车身的实时倾角值,包括:
将PID参数输入两轮平衡车控制系统;
测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
优选地,
所述拟合实时倾角值关于PID参数的当前关系式,包括:
根据所述实时倾角值与所述预设倾角阈值,计算两者的均方误差;
拟合所述均方误差关于PID参数的当前关系式。
优选地,
所述获取N组初始PID参数样本,包括:
获取PID参数的参数区间;
随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
第二方面,本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定系统,包括:
参数样本获取模块,用于实现A1:获取N组初始PID参数样本,以择一作为当前PID参数;
实时倾角获取模块,用于实现A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块,用于实现A3:判断实时倾角值是否小于预设倾角阈值;
参数样本循环模块,用于实现A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
关系拟合模块,用于实现A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
理论参数求取模块,用于实现A6:利用所述当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
理论倾角获取模块,用于实现A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块,用于实现A8:判断实时倾角值是否小于预设倾角阈值;
拟合循环模块,用于实现A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到得到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,该当前理论最小PID参数即为最终整定PID参数。
优选地,
所述实时倾角获取模块,包括:
PID参数输入单元,用于将PID参数输入两轮平衡车控制系统;
实时倾角测取单元,用于测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
优选地,
所述参数样本获取模块,包括:
参数区间获取单元,用于获取PID参数的参数区间;
参数样本随机模块,用于随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
第三方面,本发明提供一种两轮平衡车,应用于上述第一方面所述的用于两轮平衡车车身平衡控制的PID参数整定方法,包括:
参数样本获取模块,用于获取当前PID参数样本;
参数样本控制模块,用于利用所述当前PID参数样本,控制所述两轮平衡车;
实时倾角获取单元,用于获取当前PID参数样本输入两轮平衡车后车身的实时倾角值。
第四方面,本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定计算机设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如上述第一方面所述的用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
第五方面,本发明提供一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如上述第一方面所述用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定方法,包括:A1:获取N组初始PID参数样本,以择一作为当前PID参数;A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;A3:判断实时倾角值是否小于预设倾角阈值;A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;A6:利用当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;A8:判断实时倾角值是否小于预设倾角阈值;A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。本发明提供的用于两轮平衡车车身平衡控制的PID参数整定方法,能够根据两轮平衡车的在输入PID参数后的输出值,对PID参数样本进行判定,而且利用拟定输入、输出之间的关系,利用理论PID参数与实际输出之间的关系,对最终整定的PID参数进行确定,在实施方程中,不需要人为操作,能够方便、快捷地对两轮平衡车车身平衡控制的PID参数进行整定。
本发明提供的一种用于两轮平衡车车身平衡控制的PID参数整定系统、两轮平衡车、计算机设备、计算机可读存储介质也具有上述有益效果,在此不再赘述。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明一种具体实施方式所提供的一种用于两轮平衡车车身平衡控制的PID参数整定方法的流程图;
图2为本发明一种具体实施方式所提供的实时倾角值获取方法的流程图;
图3为本发明一种具体实施方式所提供的一种关系式拟合方式的流程图;
图4为本发明一种具体实施方式所提供的一种初始PID参数样本的获取流程图;
图5为本发明一种具体实施方式所提供的一种用于两轮平衡车车身平衡控制的PID参数整定方法的流程图;
图6为本发明一种具体实施方式所提供的一种两轮平衡车的硬件组成示意图;
图7为本发明一种具体实施方式所提供的一种两轮平衡车的PID调节示意图;
图8为本发明又一种具体实施方式中所提供一种用于两轮平衡车车身平衡控制的PID参数整定系统的组成示意图;
图9为本发明一种具体实施方式所提供的两轮平衡车的组成结构示意图;
图10为本发明又一种具体实施方式所提供的计算机设备的结构示意图.
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请参考图1,图1为本发明一种具体实施方式所提供的一种用于两轮平衡车车身平衡控制的PID参数整定方法的流程图。
在本发明一种具体实施方式中,本发明实施例提供一种用于两轮平衡车车身平衡控制的PID参数整定方法,包括:
步骤A1:获取N组初始PID参数样本,以择一作为当前PID参数;
在本发明实施例中,首先需要获取初始PID参数样本,在该初始PID参数样本的基础上,进行一系列的处理,例如在实际中,可以根据经验,设定PID参数,也就是对Kp、Ki、Kd三个参数的范围进行设定,可以在三个参数范围内进行均匀的取值,例如Kp可以取值Kp1、Kp2……Kpn,一共n个取值,;而Ki有m个取值、Kd有q个取值,n、m、q为大于1的整数对这些取值进行任意可能的组合。并且选择一个组合作为当前PID参数,以进行下一步处理。当然,也可以采取其他的方式进行初始PID参数样本的获取。
步骤A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
在获取当前PID参数样本后,需要知道当前PID参数样本的输入两轮平衡车后的效果,因此,可以利用当前PID参数样本去输入到两轮平衡车,并且利用当前PID参数样本去控制两轮平衡车,对于车站的平衡,可以选定一个基准,例如可以选定竖直方向作为基准,来测定车身踏板法线方向的角度,作为实时倾角值,当然,也可以选择车身踏板与水平方向的夹角值作为实时倾角值。值得说明的是,该实时倾角值,可以为一个数值,也可以为一组数值,也就是说,在将当前PID参数样本输入到两轮平衡车后,可以测取车身在一组预设时刻的每个实时倾角值。
步骤A3:判断实时倾角值是否小于预设倾角阈值;
步骤A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
对于实时倾角值,一般有一个预判,希望该实时倾角值不要太大,以免影响使用者的安全,因为可以设定一个预设倾角阈值,该预设倾角阈值对于两轮平衡车的使用者来说,应该是安全的,如果当前PID参数小于预设倾角阈值,那么可以判定该当前PID参数可以作为最终整定PID参数。如果当前PID参数大于或等于预设倾角阈值,那么可以判定该当前PID参数不可以作为最终整定PID参数,则应该继续在PID参数样本中,选定当前PID参数样本进行循环处理,如果所有的PID参数样本都不符合要求,那么可以获取N组所述初始PID参数样本对应的实时倾角值。当然,如果存在符合要求的PID参数样本,那么可以确定该PID参数样本为最终整定PID参数。
步骤A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
如果所有的PID参数样本都不符合要求,那么可以获取N组所述初始PID参数样本对应的实时倾角值,对于这些输入和输出可以进行拟合,从而得到输出与输入之间的关系式。
步骤A6:利用当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
对于该关系式,可以通过数学方法求取使得输出值最小的当前理论最小PID参数。
步骤A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
在获取当理论最小PID参数后,还要输入到两轮平衡中进行严重,也就是说,需要利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值。
步骤A8:判断实时倾角值是否小于预设倾角阈值;
步骤A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。
并且验证该实时倾角值是否符合要求,如果还不符合要求,那么可以将此事的当前理论最小PID参数及对应的实时倾角值作为一组输入输出,替换掉原本的N组初始PID参数样本所对应的输入中输出误差最大的一组,从而拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。
请参考图2,图2为本发明一种具体实施方式所提供的实时倾角值获取方法的流程图。
进一步地,为了获取所述初始PID参数样本输入两轮平衡车后车身的实时倾角值,具体地可以:
A21:将PID参数输入两轮平衡车控制系统;
A22:测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
请参考图3,图3为本发明一种具体实施方式所提供的一种关系式拟合方式的流程图。
也就是说,本发明是将PID参数样本实际输入到真实的需要控制的两轮平衡车中去,而不是采用的模拟两轮车系统,从而可以的得到更加符合实际的输入输出值的对应关系。进一步地,为了拟合实时倾角值关于PID参数的当前关系式,具体地,可以:
A51:根据所述实时倾角值与所述预设倾角阈值,计算两者的均方误差;
A52:拟合所述均方误差关于PID参数的当前关系式。
请参考图4,图4为本发明一种具体实施方式所提供的一种初始PID参数样本的获取流程图。
在上述具体实施方式的基础上,本具体实施方式中,为了获取N组初始PID参数样本,可以采用随机取值的方式,增加获取初始PID参数的快捷性。具体地可以采用以下步骤:
A11:获取PID参数的参数区间;
A12:随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
本发明提供一种用于两轮平衡车车身平衡控制的PID参数整定方法,,能够根据两轮平衡车的在输入PID参数后的输出值,对PID参数样本进行判定,而且利用拟定输入、输出之间的关系,利用理论PID参数与实际输出之间的关系,对最终整定的PID参数进行确定,在实施方程中,不需要人为操作,能够方便、快捷地对两轮平衡车车身平衡控制的PID参数进行整定。
请参考图5、图6、图7,图5为本发明一种具体实施方式所提供的一种用于两轮平衡车车身平衡控制的PID参数整定方法的流程图;图6为本发明一种具体实施方式所提供的一种两轮平衡车的硬件组成示意图;图7为本发明一种具体实施方式所提供的一种两轮平衡车的PID调节示意图。
在本发明的又一种具体实施方式中,本发明实施例可以利用MATLAB分析系统对PID参数的Kp、Ki、Kd进行自动整定。
本发明包括:两个系统,一个是两轮平衡车的硬件控制系统,包括了车身结构、MCU控制器,车轮与驱动电机、倾角检测传感器、电机转速传感器以及通讯接口;软件系统包括了基于MCU控制器的车身平衡控制系统与基于MATLAB的参数分析系统。车身平衡控制系统与MATLAB参数分析系统之间的数据传输用USB转TTL模块进行通讯。通讯方式可以采用UART方式。
车身平衡控制的硬件控制系统信息如下:
使用STM32F4103C8T6作为主控芯片,通过IIC接口读出MPU6050的平衡角速度和平衡角加速度,获取车身的倾角信息,最终通过PID输出PWM电机控制信号,由驱动完成对电机的控制。
PID参数自整定步骤如下:
车身平衡控制系统运行步骤如下:
1、完成车身平衡系统的初始化并等待接收MATLAB分析系统命令,转步骤2;
2、若接收到参数设定命令,则根据参数设置Kp、Ki、Kd,设置完成之后向MATALB分析系统发送设置成功应答;
3、若收到启动测试命令,则启动电机,进行车身平衡测试,并采用定时中断方式,每隔100MS通过MPU6050获取车身倾角数据,将数据传输至MATLAB参数分析系统,测试过程持续10S,测试完成转步骤4;
4、测试完成之后停止电机输出,等待接收MATLAB命令;
MATLAB的参数分析系统运行步骤如下:
1.完成MATLAb参数分析系统初始化并与车身平衡控制系统建立通讯连接,转步骤2;
2.给定PID参数初始区间,设Kp的参数区间为、Ki的参数区间为、Kd的参数区间为,设定样本集的个数为N,在参数区间内随机生成N个参数作为初始输入样本集X,
设定期望倾角;设定期望均方误差指标;对输入样本集的每组参数进行实际测试,转步骤3;
3.将待测试参数样本Xi发送至车身控制系统,待接收到车身平衡控制系统的应答信号之后发送启动测试命令,并实时接收车身平衡控制系统反馈倾角信息Gj,测试完成之后获得车身倾角集G,也就是在将PID参数样本输入到两轮平衡车进行PID控制后,一个预设时刻组(100个时刻)苏哦对应的车身倾角值,也就是实时倾角值:
将集G与期望倾角(也就是预设倾角阈值)比较计算测试期望均方误差作为样本输出值Yi
a.若不大于指标,则输出最终整定参数,参数自整定过程完成;
b.否则转至步骤4
4.此时获得一个完整的输入输出样本;加入到测试样本集去,判断样本集个数,
a.若个数没达到设定值N,则说明初始参数输入样本还未全部完成,则将下一个待测试参数发送到车身控制系统,转步骤3;
b.若样本集个数达到N,则跳转至步骤5
c.若样本个数达到N+1,则将新的输入输出样本代替原先样本集中样本输出值Y最大的样本,转至步骤5;
5.开始拟合模型:
建立模型:
利用MATLAB软件的lsqcurvefit()函数对N个样本集进行参数估计,求出之后获得数学模型F;
根据样本集X各个参数的最大最小值重新设定Kp,Ki,Kd的约束条件:
求三元函数模型F在约束条件下使得最小值点Xi,过程如下:
求F函数一次偏导方程组的解,代入公式:,
求出的值,
若大于0则令,
若小于0则在约束条件边界取使成立,令
将Xi作为待测试参数,转至步骤3;
依次循环进行直至步骤3.a,则整个参数自整定过程完成。
发明实施例提供一种基于MATLAB软件对控制系统的PID参数进行自动调节的方法,整个PID参数整定过程中无需人为干预即可完成。在参数自整定过程中,对实际的车身平衡控制系统不断进行实验测试,直到平衡指标达到要求为止,从而得到与当前实际的车身平衡控制系统契合度最佳的PID参数,即Kp、Ki、Kd。
请参考图8,图8为本发明又一种具体实施方式中所提供一种用于两轮平衡车车身平衡控制的PID参数整定系统的组成示意图。
在本发明又一种具体实施方式中,本发明实施例提供一种用于两轮平衡车车身平衡控制的PID参数整定系统800,包括:
参数样本获取模块810,用于实现A1:获取N组初始PID参数样本,以择一作为当前PID参数;
实时倾角获取模块820,用于实现A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块830,用于实现A3:判断实时倾角值是否小于预设倾角阈值;
参数样本循环模块840,用于实现A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
关系拟合模块850,用于实现A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
理论参数求取模块860,用于实现A6:利用所述当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
理论倾角获取模块870,用于实现A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块880,用于实现A8:判断实时倾角值是否小于预设倾角阈值;
拟合循环模块890,用于实现A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到得到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,该当前理论最小PID参数即为最终整定PID参数。
优选地,
所述实时倾角获取模块,包括:
PID参数输入单元,用于将PID参数输入两轮平衡车控制系统;
实时倾角测取单元,用于测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
优选地,
所述参数样本获取模块,包括:
参数区间获取单元,用于获取PID参数的参数区间;
参数样本随机模块,用于随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
请参考图9,图9为本发明一种具体实施方式所提供的两轮平衡车的组成结构示意图。
在本发明一种具体实施方式中,本发明实施例提供一种两轮平衡车900,应用于上述任一种具体实施方式中所述的用于两轮平衡车车身平衡控制的PID参数整定方法,包括:
参数样本获取模块910,用于获取当前PID参数样本;
参数样本控制模块920,用于利用所述当前PID参数样本,控制所述两轮平衡车;
实时倾角获取单元930,用于获取当前PID参数样本输入两轮平衡车后车身的实时倾角值。
请参考图10,图10为本发明又一种具体实施方式所提供的计算机设备的结构示意图。
在本发明的又一种具体实施方式中,本发明实施例提供一种计算机设备,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现上述任一种具体实施方式所述的用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
下面参考图10,其示出了适于用来实现本申请实施例的计算机设备的结构示意图。图10示出的计算机设备仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图10所示,计算机系统1000包括处理器(CPU)1001,其可以根据存储在只读存储器(ROM)1002中的程序或者从存储部分1008加载到随机访问存储器(RAM)1003中的程序而执行各种适当的动作和处理。在RAM 1003中,还存储有系统1000操作所需的各种程序和数据。
CPU 1001、ROM 1002以及RAM 1003通过总线1004彼此相连。输入/输出(I/O)接口1003也连接至总线1004。
以下部件连接至I/O接口1005:包括键盘、鼠标等的输入部分1006;包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分1007;包括硬盘等的存储部分1008;以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分1009。通信部分1009经由诸如因特网的网络执行通信处理。驱动器1010也根据需要连接至I/O接口1007。可拆卸介质1011,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器1010上,以便于从其上读出的计算机程序根据需要被安装入存储部分1008。
特别地,根据本公开的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本公开的实施例包括一种计算机程序产品,其包括承载在计算机可读介质上的计算机程序,该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分1009从网络上被下载和安装,和/或从可拆卸介质1011被安装。在该计算机程序被处理器(CPU)1001执行时,执行本申请的方法中限定的上述功能。需要说明的是,本申请所述的计算机可读介质可以是计算机可读信号介质或者计算机可读介质或者是上述两者的任意组合。计算机可读介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读介质的更具体的例子可以包括但不限于:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本申请中,计算机可读介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本申请中,计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括但不限于:无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本申请的操作的计算机程序代码,所述程序设计语言包括面向目标的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如”C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。
在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
附图中的流程图和框图,图示了按照本申请各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
作为本发明的又一具体实施方式,本发明实施例提供一种计算机可读存储介质所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现上述任意具体实施方式中的用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
该计算机可读介质可以是上述实施例中描述的计算机或终端设备中所包含的;也可以是单独存在,而未装配入该计算机设备中。上述计算机可读介质承载有一个或者多个程序,当上述一个或者多个程序被该计算机设备执行时,使得该计算机设备:A1:获取N组初始PID参数样本,以择一作为当前PID参数;A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;A3:判断实时倾角值是否小于预设倾角阈值;A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;A6:利用当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;A8:判断实时倾角值是否小于预设倾角阈值;A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实施倾角值与预设倾角阈值相差最大的初始PID参数样本做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上对本发明所提供的一种用于两轮平衡车车身平衡控制的PID参数整定方法及产品进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种用于两轮平衡车车身平衡控制的PID参数整定方法,其特征在于,包括:
A1:获取N组初始PID参数样本,以择一作为当前PID参数;
A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
A3:判断实时倾角值是否小于预设倾角阈值;
A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
A6:利用当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
A8:判断实时倾角值是否小于预设倾角阈值;
A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实时倾角值与预设倾角阈值相差最大的初始PID参数样本及其对应的实时倾角值做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,确定该当前理论最小PID参数为最终整定PID参数。
2.根据权利要求1所述的用于两轮平衡车车身平衡控制的PID参数整定方法,其特征在于,
所述获取所述初始PID参数样本输入两轮平衡车后车身的实时倾角值,包括:
将PID参数输入两轮平衡车控制系统;
测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
3.根据权利要求2所述的用于两轮平衡车车身平衡控制的PID参数整定方法,其特征在于,
所述拟合实时倾角值关于PID参数的当前关系式,包括:
根据所述实时倾角值与所述预设倾角阈值,计算两者的均方误差;
拟合所述均方误差关于PID参数的当前关系式。
4.根据权利要求1至3任一项所述的用于两轮平衡车车身平衡控制的PID参数整定方法,其特征在于,
所述获取N组初始PID参数样本,包括:
获取PID参数的参数区间;
随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
5.一种用于两轮平衡车车身平衡控制的PID参数整定系统,其特征在于,包括:
参数样本获取模块,用于实现A1:获取N组初始PID参数样本,以择一作为当前PID参数;
实时倾角获取模块,用于实现A2:利用当前PID参数样本,获取当前PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块,用于实现A3:判断实时倾角值是否小于预设倾角阈值;
参数样本循环模块,用于实现A4:如果否,则从N组初始PID参数样本确定当前PID参数,并进入步骤A2,直到获取N组所述初始PID参数样本对应的实时倾角值;
关系拟合模块,用于实现A5:根据N组所述初始PID参数样本及其对应的实时倾角值,拟合实时倾角值关于PID参数的当前关系式;
理论参数求取模块,用于实现A6:利用所述当前关系式,求取使得实时倾角值最接近所述预设倾角阈值的当前理论最小PID参数;
理论倾角获取模块,用于实现A7:利用当前理论最小PID参数,获取当前理论最小PID参数样本输入两轮平衡车后车身的实时倾角值;
倾角阈值判断模块,用于实现A8:判断实时倾角值是否小于预设倾角阈值;
拟合循环模块,用于实现A9:如果否,则将当前理论最小PID参数及其对应的实时倾角值与N组所述初始PID参数样本及其对应的实时倾角值中实时倾角值与预设倾角阈值相差最大的初始PID参数样本及其对应的实时倾角值做替换,得到替换后的N组PID参数样本及其对应的实时倾角值,并拟合实时倾角值关于PID参数的当前关系式,进入步骤A6,直到得到当前理论最小PID参数对应的实时倾角值小于预设倾角阈值,该当前理论最小PID参数即为最终整定PID参数。
6.根据权利要求5所述的用于两轮平衡车车身平衡控制的PID参数整定系统,其特征在于,
所述实时倾角获取模块,包括:
PID参数输入单元,用于将PID参数输入两轮平衡车控制系统;
实时倾角测取单元,用于测取PID参数输入两轮平衡车控制系统后,M个预设时刻的实时倾角值;其中,所述M为正整数。
7.根据权利要求5或6所述的用于两轮平衡车车身平衡控制的PID参数整定系统,其特征在于,
所述参数样本获取模块,包括:
参数区间获取单元,用于获取PID参数的参数区间;
参数样本随机模块,用于随机在所述参数区间内生成N组PID参数作为初始PID参数样本。
8.一种两轮平衡车,应用于如权利要求1至4任一项所述的用于两轮平衡车车身平衡控制的PID参数整定方法,其特征在于,包括:
参数样本获取模块,用于获取当前PID参数样本;
参数样本控制模块,用于利用所述当前PID参数样本,控制所述两轮平衡车;
实时倾角获取单元,用于获取当前PID参数样本输入两轮平衡车后车身的实时倾角值。
9.一种用于两轮平衡车车身平衡控制的PID参数整定计算机设备,其特征在于,包括:
存储器,用于存储计算机程序;
处理器,用于执行所述计算机程序时实现如权利要求1至4任一项所述用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至4任一项所述用于两轮平衡车车身平衡控制的PID参数整定方法的步骤。
CN201811513257.1A 2018-12-11 2018-12-11 用于两轮平衡车车身平衡控制的pid参数整定方法及产品 Expired - Fee Related CN109507871B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811513257.1A CN109507871B (zh) 2018-12-11 2018-12-11 用于两轮平衡车车身平衡控制的pid参数整定方法及产品

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811513257.1A CN109507871B (zh) 2018-12-11 2018-12-11 用于两轮平衡车车身平衡控制的pid参数整定方法及产品

Publications (2)

Publication Number Publication Date
CN109507871A CN109507871A (zh) 2019-03-22
CN109507871B true CN109507871B (zh) 2022-03-25

Family

ID=65752196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811513257.1A Expired - Fee Related CN109507871B (zh) 2018-12-11 2018-12-11 用于两轮平衡车车身平衡控制的pid参数整定方法及产品

Country Status (1)

Country Link
CN (1) CN109507871B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110793720B (zh) * 2019-11-05 2021-05-04 腾讯科技(深圳)有限公司 用于标定机械系统平衡位置的方法、装置及存储介质
CN111103790A (zh) * 2019-12-05 2020-05-05 珠海格力电器股份有限公司 Pid控制器的参数整定方法、装置、存储介质、终端及系统
CN110941280A (zh) * 2019-12-16 2020-03-31 华南理工大学广州学院 一种激光追踪平衡车控制方法
CN110941281A (zh) * 2019-12-16 2020-03-31 华南理工大学广州学院 一种激光追踪平衡车控制系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833314A (zh) * 2010-03-30 2010-09-15 深圳达实智能股份有限公司 污水处理控制系统及污水处理控制方法
CN103529850A (zh) * 2013-10-28 2014-01-22 华南理工大学广州学院 两轮自平衡车控制方法
CN104298113A (zh) * 2014-10-22 2015-01-21 五邑大学 两轮机器人自适应模糊平衡控制器
CN105807689A (zh) * 2016-04-29 2016-07-27 上海海洋大学 一种两轮自平衡车控制系统
CN106444363A (zh) * 2016-12-14 2017-02-22 浙江中控技术股份有限公司 一种pid参数整定方法及整定系统
CN107340298A (zh) * 2017-06-29 2017-11-10 广西师范大学 基于摄像头路面检测的平衡车系统测控方法
CN108931986A (zh) * 2018-06-19 2018-12-04 福建海源自动化机械股份有限公司 一种两轮汽车自平衡控制方法、装置及存储介质

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9536149B1 (en) * 2016-02-04 2017-01-03 Proxy Technologies, Inc. Electronic assessments, and methods of use and manufacture thereof

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833314A (zh) * 2010-03-30 2010-09-15 深圳达实智能股份有限公司 污水处理控制系统及污水处理控制方法
CN103529850A (zh) * 2013-10-28 2014-01-22 华南理工大学广州学院 两轮自平衡车控制方法
CN104298113A (zh) * 2014-10-22 2015-01-21 五邑大学 两轮机器人自适应模糊平衡控制器
CN105807689A (zh) * 2016-04-29 2016-07-27 上海海洋大学 一种两轮自平衡车控制系统
CN106444363A (zh) * 2016-12-14 2017-02-22 浙江中控技术股份有限公司 一种pid参数整定方法及整定系统
CN107340298A (zh) * 2017-06-29 2017-11-10 广西师范大学 基于摄像头路面检测的平衡车系统测控方法
CN108931986A (zh) * 2018-06-19 2018-12-04 福建海源自动化机械股份有限公司 一种两轮汽车自平衡控制方法、装置及存储介质

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
Genetic fuzzy self-tuning PID controllers for antilock braking systems;Abdel Badie Sharkawy;《Web of Science》;20101231;第23卷;第1041–1052页 *
PID控制原理及参数整定方法;金奇,邓志杰;《重庆工学院学报(自然科学)》;20080531;第22卷(第5期);第91-94页 *
PID控制参数整定方法;刘镇,姜学智,李东海;《电子技术应用》;19971231(第5期);第4-6页 *
Robust task-space control of an autonomous underwater vehicle-manipulator system by PID-like fuzzy control scheme with disturbance estimator;P.S. Londhe,等;《Ocean Engineering》;20171231;第139卷;第1-13页 *
双轮自平衡车的双闭环式PID控制系统设计与实现;林枫,蔡延光;《工业控制计算机》;20171231;第30卷(第6期);第73-75页 *
基于PLC的调速控制系统算法及机械调控系统PID参数仿真整定的分析与研究;严继超;《机械化工》;20161231(第13期);第45页 *
基于STM32的智能平衡车控制系统设计;高正中,龚群英,宋森森;《现代电子技术》;20160715;第39卷(第14期);第46-48页 *
浅析PID参数整定方法;戴雅馨;《纯碱工业》;20091231(第6期);第15-17页 *

Also Published As

Publication number Publication date
CN109507871A (zh) 2019-03-22

Similar Documents

Publication Publication Date Title
CN109507871B (zh) 用于两轮平衡车车身平衡控制的pid参数整定方法及产品
US9940223B2 (en) Human-machine interface test system
CN110658542B (zh) 自动驾驶汽车定位识别方法、装置、设备及存储介质
EP3913727A1 (en) Battery residual value determination system
Mollov et al. μ-synthesis and hardware-in-the-loop simulation of miniature helicopter control system
CN107220170A (zh) 一种测试方法、装置及电子设备
CN112140894B (zh) 续驶里程处理方法、装置及车辆
CN104699594A (zh) 一种测试方法和测试设备
CN114936107A (zh) 车辆整车数据校验方法、装置、电子设备及存储介质
CN101713636B (zh) 三维坐标系补正及绘制系统与方法
CN110717315A (zh) 系统数据批量修改方法、装置、存储介质及电子设备
CN113001546B (zh) 一种提高工业机器人运动速度安全性的方法及系统
GB2583589A (en) Interactive guidance system for selecting thermodynamics methods in process simulations
CN112365883A (zh) 座舱系统语音识别测试方法、装置、设备及存储介质
CN109813934A (zh) 加速轴向的校准方法、装置及计算机可读存储介质
CN112051074B (zh) 车辆标定参数获取方法、系统及计算机可读存储介质
Chapoutot et al. Static analysis of Simulink programs
US20230286621A1 (en) Ship state estimation device, ship state estimation system, ship state estimation method, and ship state estimation program
CN116008379B (zh) 基于模型拟合和机器学习的自动滴定系统、方法和装置
CN102043759A (zh) 一种数学模型数值计算程序的验证方法
CN112036565B (zh) 一种数字工作记忆广度的自适应训练方法及装置
CN116279535A (zh) 车辆控制方法和装置、电子设备、计算机存储介质
Zhmud et al. Numerical Optimization of Locked Automatic Control System in Software VisSim: New Structures and Methods
CN116879579B (zh) 车辆加速度计算方法、装置、计算机设备及存储介质
US20210026999A1 (en) Method and device for validating a simulation of a technical system

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220325