CN112835286A - Pid参数自动整定方法及系统 - Google Patents

Pid参数自动整定方法及系统 Download PDF

Info

Publication number
CN112835286A
CN112835286A CN202011355027.4A CN202011355027A CN112835286A CN 112835286 A CN112835286 A CN 112835286A CN 202011355027 A CN202011355027 A CN 202011355027A CN 112835286 A CN112835286 A CN 112835286A
Authority
CN
China
Prior art keywords
pid parameter
controlled system
pid
parameter combination
value
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
CN202011355027.4A
Other languages
English (en)
Other versions
CN112835286B (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.)
Jiangxi University of Science and Technology
Original Assignee
Jiangxi 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 Jiangxi University of Science and Technology filed Critical Jiangxi University of Science and Technology
Priority to CN202011355027.4A priority Critical patent/CN112835286B/zh
Publication of CN112835286A publication Critical patent/CN112835286A/zh
Application granted granted Critical
Publication of CN112835286B publication Critical patent/CN112835286B/zh
Active 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)
  • Feedback Control In General (AREA)

Abstract

本发明提供一种PID参数自动整定方法及系统,属于自动控制领域。所述方法包括:S1)设置超参数取值和PID参数组合初值;S2)令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,计算N次往返阶跃的平均误差积分;其中N为正整数,且N>1;S3)将所述PID参数组合作为优化变量,将获取所述平均误差积分的最小值作为优化目标,构成一个非约束最优化问题,判断所述PID参数组合是否满足最优化收敛条件,若不满足,根据所述求解算法更新所述PID参数组合,重复步骤S2),直到所述PID参数组合满足所述求解算法的最优化收敛条件;S4)输出满足最优化收敛条件的PID参数组合。保证具有强烈非线性、快速动态响应或对噪声敏感的受控系统的控制及时性和准确性。

Description

PID参数自动整定方法及系统
技术领域
本发明涉及自动控制领域,具体地涉及一种PID参数自动整定方法及一种PID参数自动整定系统。
背景技术
PID参数整定在进行受控系统的精准控制时具有很重要的作用。经过相对漫长的技术发展,目前进行PID参数整定的成熟方法包括临界比例度法、反应曲线法、衰减曲线法和经验法等。由于以上常规PID参数整定方法来源于普遍经验之总结,所以针对某些特殊受控系统不一定完全适用。例如,具有强烈非线性的受控系统由于不稳定性而导致临界比例度很窄,进而在参数微调时,临界比例度法往往达不到理想的整定效果;反应曲线法有很强的适用广泛性,但对具有强烈非线性的受控系统进行PID参数整定时,很难准确套用反应曲线模型,影响参数整定效果;衰减曲线法在整定噪声敏感的受控系统时,由于干扰频繁、反应曲线不规则或不断有小摆动,无法得到准确的衰减比例度和衰减周期,导致整定效果不理想;经验法适用于动态响应适中的受控系统,若受控系统的动态响应较快,工程师无法做出及时和准确的参数整定,容易导致控制系统失去稳定性。
所以,以上常规PID参数整定方法并不适用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,例如,磁悬浮系统、无人飞行器系统和平衡车等。为了解决这种问题,需要创造一种适用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统的PID参数自动整定方法。
发明内容
本发明实施方式的目的是提供一种PID参数自动整定方法,以至少解决上述的现有PID参数整定方法不适用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统的问题。
为了实现上述目的,本发明第一方面提供一种PID参数自动整定方法,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,所述方法包括:S1)根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值;S2)令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据,并计算N次往返阶跃的平均误差积分;其中N为正整数,且N>1;S3)将所述PID参数组合作为优化变量,将获取所述平均误差积分的最小值作为优化目标,构成一个非约束最优化问题,按照非约束最优化问题的已有求解算法,判断所述PID参数组合是否满足最优化收敛条件,若不满足,根据所述求解算法更新所述PID参数组合,重复步骤S2),直到所述PID参数组合满足所述求解算法的最优化收敛条件;S4)输出满足最优化收敛条件的PID参数组合,完成PID参数整定。
可选的,步骤S1)中,所述超参数包括:PID控制器的采样周期、积分增益、参考设定值X、参考设定值阶跃幅度ΔX、误差积分时间和阶跃时间间隔;其中,所述受控系统的所述参考设定值为用户设定的参考状态值。
可选的,步骤S1)中,所述PID参数组合包括:比例增益和微分增益。
可选的,步骤S2)中,往返阶跃过程指所述被控系统的所述参考设定值在原参考设定值X和所述参考设定值阶跃幅度ΔX的基础上,以阶跃时间间隔在X-ΔX和X+ΔX之间阶跃变化;其中,所述阶跃时间间隔不小于所述误差积分时间。
可选的,步骤S2)中,所述误差积分为平方误差积分、时间乘平方误差积分、绝对误差积分准则或时间乘绝对误差积分中的任意一种。
可选的,步骤S3)中,所述平均误差积分的计算方法包括:令受控系统按照所述阶跃幅度和所述阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据;其中,所述往返阶跃数据之构成为所述受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系;将所述N组往返阶跃数据分为N条上升阶跃响应曲线和N条下降阶跃响应曲线,其中所述上升阶跃响应曲线的所述参考设定值为X+ΔX,所述下降阶跃响应曲线的所述参考设定值为X-ΔX;分别计算所述上升阶跃响应曲线的平均响应曲线和所述下降阶跃响应曲线的平均响应曲线;根据所述上升阶跃平均响应曲线和所述下降阶跃平均响应曲线,分别计算上升阶跃的误差积分和下降阶跃的误差积分,并求两者的平均值作为所述N次往返阶跃的平均误差积分。
可选的,步骤S3)中,步骤S3)中,所述非约束最优化问题的已有求解算法至少包括:最速下降法、梯度下降法、牛顿法、拟牛顿法和人工智能法。
可选的,步骤S3)中,所述根据所述求解算法更新所述PID参数组合,包括:获取所述非约束最优化问题的收敛情况;根据所述求解算法和所述收敛情况计算所述PID参数组合的修正值;根据所述修正值更新所述PID参数组合。
本发明第二方面提供一种PID参数自动整定系统,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,所述系统包括:传感器单元,用于获取所述受控系统的实时状态;PID控制器单元,用于根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值;PID参数整定单元,用于根据所述超参数和所述PID参数组合使受控系统进行若干次往返阶跃,获得每一次阶跃数据并计算平均误差积分;所述PID参数整定单元还被用于判断所述PID参数组合是否满足最优化收敛条件和根据所述求解算法更新所述PID参数组合;执行单元,用于输出满足最优化收敛条件的PID参数组合。
另一方面,本发明提供一种计算机可读储存介质,该计算机可读存储介质上储存有指令,其在计算机上运行时使得计算机执行上诉中任一项所述的PID参数自动整定方法。
通过上述技术方案,根据受控系统的实际运行效果总结受控系统的受控数学模型,根据数学模型设置受控系统的超参数和PID整定参数组合,通过在理想受控系统的参考设定值上下波动,反复在理想参考设定值附件进行整定控制并获得每次整定过程的误差积分,将多组误差积分的平均值作为该组PID整定参数组合的误差积分,并将误差积分最小值作为一个全局最优化问题,判断是否满足最优化收敛条件,当满足最优化收敛条件时判断设定的该组PID整定参数组合为最优化参数组合。保证具有强烈非线性、快速动态响应或对噪声敏感的受控系统的控制及时性和准确性。
本发明实施方式的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施方式的进一步理解,并且构成说明书的一部分,与下面的具体实施方式一起用于解释本发明实施方式,但并不构成对本发明实施方式的限制。在附图中:
图1是本发明一种实施方式提供的PID参数自动整定方法的方法流程图;
图2是本发明一种实施方式提供的平均误差积分获取方法的方法流程图;
图3是本发明一种实施方式提供的往返阶跃示意图;
图4是本发明一种实施方式提供的梯度下降法求解最优化问题的方法示意图;
图5是本发明一种实施方式提供的PID参数自动整定方法的控制流程图;
图6是本发明一种实施方式提供的PID参数自动整定系统的系统结构图。
附图标记说明
10-传感器单元;20-PID控制单元;30-PID参数整定单元;40-执行单元。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明,并不用于限制本发明。
图6是本发明一种实施方式提供的PID参数自动整定系统。如图6所述,本发明实施方式提供一种PID参数自动整定系统,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,所述系统包括:传感器单元10,用于获取所述受控系统的实时状态;PID控制器单元20,用于根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值;PID参数整定单元30,用于根据所述超参数和所述PID参数组合使受控系统进行若干次往返阶跃,获得每一次阶跃数据并计算平均误差积分;所述PID参数整定单元还被用于判断所述PID参数组合是否满足最优化收敛条件和根据所述求解算法更新所述PID参数组合。执行单元40,用于输出满足最优化收敛条件的PID参数组合。
图1是本发明一种实施方式提供的PID参数自动整定方法的方法流程图。如图1所示,本发明实施方式提供一种PID参数自动整定的方法,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,所述方法包括:
步骤S10:根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值。
具体的,超参数包括:PID控制器的采样周期、积分增益、参考设定值X、参考设定值阶跃幅度ΔX、误差积分时间和阶跃时间间隔;其中,受控系统的参考设定值为用户设定的参考状态值;PID参数组合初值包括比例增益和微分增益。其中,比例增益、积分增益和微分增益三个参数(kp,ki,kd)服务于增量式数字PID控制器,数学表达模型为:
Δuk=uk-uk-1=kp(ek-ek-1)+kiek+kd(ek-2ek-1+ek-2)
其中,ek=xk-X为受控系统的第k个采样时刻的误差值,xk为受控系统的第k个采样时刻的状态值,uk为控制器的第k个采样时刻的输出值。
本方法的核心目的就是求解最优PID参数组合(kp,kd)。在PID参数整定时,根据控制经验将积分增益ki设定为常数,积分增益的数值大小与控制器优化后的鲁棒性直接相关,即ki越大,控制器经优化后的鲁棒性越好。所以积分增益的设定值由受控系统控制器的鲁棒性需求决定,例如,磁悬浮系统和无人飞行器具有非常不同的系统特性:磁悬浮系统的控制目标是维持稳定悬浮、具备抗干扰能力,需要保证控制器的鲁棒性,则可以在PID参数整定过程中采用较大的积分增益;而无人飞行器的稳定性要求相对较低,可以在PID参数整定过程中采用较小的积分增益。
除了上述的积分增益,参考设定值阶跃幅度ΔX与优化后控制器的鲁棒性也直接相关。参考设定值阶跃幅度ΔX越大,表示扰动幅度越大,则控制器经过优化后的鲁棒性也就越好;但是,参考设定值阶跃幅度ΔX不宜过大,避免扰动幅过大与强烈非线性二者叠加,造成整定结果与实际参考设定值X的情况偏差过大。考虑应用本发明提出的PID参数自动整定方法的受控系统具有的非线性特性,优选的,将设定值阶跃幅度ΔX设置为参考设定值最大变化范围的1~10%。
所以,根据受控系统优化后控制器的鲁棒性需求,综合积分增益和参考设定值阶跃幅度两者的影响因素,完成积分增益ki和参考设定值阶跃幅度ΔX设定。之后,根据设定的ki和ΔX进行误差积分时间设定,设定原则为确保受控系统有足够调整时间而达到新的稳定平衡,优选设定为ΔX/ki的1~10倍。
另外,根据受控系统特性和控制经验进行PID控制器的采样周期、参考设定值X、和阶跃时间间隔设定,该部分数据优选通过相关人员录入PID参数自动整定系统,保证设置的参数合理。
完成超参数设置后,根据上述增量式数字PID控制器的数学表达式和设定的积分增益ki,尝试并设置一组可行的PID参数组合,作为PID参数组合初值。
在本发明实施例中,根据受控系统特性和控制经验设置受控系统的超参数和PID参数组合,在保证满足强烈非线性、快速动态响应或对噪声敏感的受控系统稳定控制的前提下,兼顾了优化后的控制鲁棒性。
步骤S20:令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据,并计算N次往返阶跃的平均误差积分;其中N为正整数,且N>1,具体的,如图2,包括以下步骤:
步骤S201:令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据。
具体的,如图3,根据预设的超参数和PID参数组合,令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,在多次往返阶跃过程中,参考设定值X根据参考设定值阶跃幅度ΔX交替更新,相邻两个更新时间点的间隔等于阶跃时间间隔,其中,参考设定值XSP更新规则满足以下关系:
Figure BDA0002802317730000081
其中,XSP,上升为上升阶跃的参考设定值;XSP,下降为下降阶跃的参考设定值。一次往返阶跃过程包括一次上升阶跃和一次下降阶跃,即在一次往返阶跃过程中,参考设定值需要根据上述关系更新二次。根据阶跃时间间隔使参考设定值在XSP,上升和XSP,下降之间反复阶跃,从而令受控系统进行多次往返阶跃,优选的,往返阶跃次数N由控制人员依据经验设定,根据受控系统对噪声的敏感程度进行适应性得增减。例如,预设往返阶跃次数为N次,则在整个往返阶跃过程中,共进行N次上升阶跃和N次下降阶跃。每完成一次上升阶跃或一次下降阶跃便获得一组阶跃数据,其中,阶跃数据为受控系统从阶跃开始至误差积分时间内受控系统状态与时间的对应关系,其构成为受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系。所以最终共获得N组上升阶跃数据和N组下降阶跃数据。
在本发明实施例中,通过在参考设定值初值附近进行小幅往返阶跃可以有效避免强烈非线性对受控系统运行的干扰,误差积分可以适应受控系统的快速动态响应,多次往返提高了PID参数整定在噪声敏感条件下的精度和效率。
步骤S202:计算N次往返阶跃的平均误差积分;其中N为正整数,且N>1。
具体的,如步骤S201所述,经过N次往返阶跃后,共获得N组上升阶跃数据和N组下降阶跃数据,即N个上升阶跃过程中受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系和N个下降阶跃过程中受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系,将每一个函数关系转换为一个对应的二维坐标曲线,则共获得N条上升阶跃响应曲线和N条下降阶跃响应曲线。分别计算N条上升阶跃响应曲线的平均响应曲线和N条下降阶跃响应曲线的平均响应曲线,计算规则为:
Figure BDA0002802317730000091
Figure BDA0002802317730000092
为受控系统的在第i条上升(或下降)阶跃响应曲线第k个采样时刻的状态值。
通过计算获得上升阶跃平均响应曲线和下降阶跃平均响应曲线后,通过上升阶跃平均响应曲线和下降阶跃平均响应曲线分别计算上升阶跃的误差积分和下降阶跃的误差积分。计算规则为:
Figure BDA0002802317730000093
其中,
Figure BDA0002802317730000094
为时间乘平方误差积分,
Figure BDA0002802317730000095
为每条阶跃响应曲线的数据总数,Δt为采样时间间隔。获得对应的上升阶跃的误差积分和下降阶跃的误差积分,最后求取二者的平均误差积分
Figure BDA0002802317730000096
计算公式为:
Figure BDA0002802317730000097
所获得的误差积分平均值
Figure BDA0002802317730000098
就作为本次设定的比例增益kp和微分增益kd的PID参数组合(kp,kd)的平均误差积分,记为
Figure BDA0002802317730000099
在另一种可能的实施方式中,关于N数值的确定,优选的,每完成一次往返阶跃,根据上述关系式更新平均误差积分,将获得的新的平均误差积分与前一次平均误差积分进行比较,计算二者的差值,当差值小于预设值时,记为一次趋于稳定关系,若连续几个新的平均误差积分均满足趋于稳定关系,且次数达到预设值时,判定为整定趋于稳定,自动停止往返阶跃,输出最终的平均误差积分。例如,在完成5轮整定后,误差积分平均值更新到了前5轮的误差积分平均值
Figure BDA0002802317730000101
进行第6轮整定后,获得前六轮新的平均误差积分
Figure BDA0002802317730000102
判断
Figure BDA0002802317730000103
Figure BDA0002802317730000104
之间的差值不大于预设值,记为一次趋于稳定关系,即
Figure BDA0002802317730000105
继续进行整定,获得
Figure BDA0002802317730000106
和,
Figure BDA0002802317730000107
判断连续第7轮和第8轮均为趋于稳定关系,在连续次数达到预设值时,例如预设连续3次均为趋于稳定关系时,系统在自动在预设轮数后结束新一轮整定,获得误差积分平均值。系统进行适应性自动设定整定轮数,提高系统的智能性。
步骤S30:将所述PID参数组合作为优化变量,将获取所述平均误差积分的最小值作为优化目标,构成一个非约束最优化问题,按照非约束最优化问题的已有求解算法,判断所述PID参数组合是否满足最优化收敛条件,若不满足,根据所述求解算法更新所述PID参数组合,重复步骤S20,直到所述PID参数组合满足所述求解算法的最优化收敛条件。
常用的求解算法很多,例如,最速下降法、梯度下降法、牛顿法、拟牛顿法和人工智能法。在一种可能的实施方式中,本实施例以梯度下降法来进行举例说明。
具体的,为了有效评估控制系统的鲁棒性,优选的,如图4,将平均误差积分的最小值作为一个非约束最优化问题,按照非约束最优化问题的已有求解算法,判断本次PID参数组合是否满足最优化收敛条件。
在一种可能的实施方式中,使用梯度下降法求解非约束最优化问题,通过梯度下降法来一步步的迭代求解,得到最小化的平均误差积分和最优PID参数组合。梯度下降法的迭代公式为,
Figure BDA0002802317730000108
其中,Kj=(kp,j,kd,j),是第j次迭代采用的PID参数组合,αj是迭代步长,
Figure BDA0002802317730000111
是平均误差积分在Kj处的梯度。其中,
Figure BDA0002802317730000112
的表达式为,
Figure BDA0002802317730000113
Figure BDA0002802317730000114
Figure BDA0002802317730000115
换句话说,为了求
Figure BDA0002802317730000116
我们还需要做实验,求另外四个相邻PID参数组合下的平均误差积分。而Δkp和Δkd取决于收敛条件。例如,若|kp,j+1-kp,j|<Δkp,收敛和|kd,j+1-kd,j|<Δkd,收敛为收敛条件,则Δkp≈Δkp,收敛和Δkd≈Δkd,收敛
在一种可能的实施方式中,如图5,PID参数整定系统运用于磁悬浮系统,首先传感器单元10获取受控系统的控制特性,根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值,然后令受控系统按照阶跃幅度和阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据;其中,往返阶跃数据之构成为所述受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系;将N组往返阶跃数据分为N条上升阶跃响应曲线和N条下降阶跃响应曲线,其中上升阶跃响应曲线的所述参考设定值为X+ΔX,下降阶跃响应曲线的所述参考设定值为X-ΔX;分别计算上升阶跃响应曲线的平均响应曲线和下降阶跃响应曲线的平均响应曲线;根据上升阶跃平均响应曲线和下降阶跃平均响应曲线,分别计算上升阶跃的误差积分和下降阶跃的误差积分,并求两者的平均值作为N次往返阶跃的平均误差积分。通过梯度下降法获取非约束最优化问题的收敛情况;根据梯度下降法算法和收敛情况计算所述PID参数组合的修正值;根据修正值更新所述PID参数组合。重复进行上述步骤,直到非约束最优化问题满足收敛条件,输出本次(kp,kd)作为最优参数组合。
本发明实施方式还提供一种计算机可读储存介质,该计算机可读存储介质上储存有指令,其在计算机上运行时使得计算机执行上诉中任一项所述的PID参数自动整定方法。
本领域技术人员可以理解实现上述实施方式的方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序存储在一个存储介质中,包括若干指令用以使得单片机、芯片或处理器(processor)执行本发明各个实施方式所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上结合附图详细描述了本发明的可选实施方式,但是,本发明实施方式并不限于上述实施方式中的具体细节,在本发明实施方式的技术构思范围内,可以对本发明实施方式的技术方案进行多种简单变型,这些简单变型均属于本发明实施方式的保护范围。另外需要说明的是,在上述具体实施方式中所描述的各个具体技术特征,在不矛盾的情况下,可以通过任何合适的方式进行组合。为了避免不必要的重复,本发明实施方式对各种可能的组合方式不再另行说明。
此外,本发明的各种不同的实施方式之间也可以进行任意组合,只要其不违背本发明实施方式的思想,其同样应当视为本发明实施方式所公开的内容。

Claims (10)

1.一种PID参数自动整定方法,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,其特征在于,所述方法包括:
S1)根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值;
S2)令受控系统以阶跃幅度和阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据,并计算N次往返阶跃的平均误差积分;其中N为正整数,且N>1;
S3)将所述PID参数组合作为优化变量,将获取所述平均误差积分的最小值作为优化目标,构成一个非约束最优化问题,按照非约束最优化问题的已有求解算法,判断所述PID参数组合是否满足最优化收敛条件,若不满足,根据所述求解算法更新所述PID参数组合,重复步骤S2),直到所述PID参数组合满足所述求解算法的最优化收敛条件;
S4)输出满足最优化收敛条件的PID参数组合,完成PID参数整定。
2.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S1)中,所述超参数包括:PID控制器的采样周期、积分增益、参考设定值X、参考设定值阶跃幅度ΔX、误差积分时间和阶跃时间间隔;其中,所述受控系统的所述参考设定值为用户设定的参考状态值。
3.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S1)中,所述PID参数组合包括:比例增益和微分增益。
4.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S2)中,往返阶跃过程指所述被控系统的所述参考设定值在原参考设定值X和所述参考设定值阶跃幅度ΔX的基础上,以阶跃时间间隔在X-ΔX和X+ΔX之间阶跃变化;其中,所述阶跃时间间隔不小于所述误差积分时间。
5.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S2)中,所述误差积分为平方误差积分、时间乘平方误差积分、绝对误差积分准则或时间乘绝对误差积分中的任意一种。
6.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S3)中,所述平均误差积分的计算方法包括:
令受控系统按照所述阶跃幅度和所述阶跃时间间隔进行N次往返阶跃,获得受控系统每次的阶跃数据;其中,所述往返阶跃数据之构成为所述受控系统从阶跃开始至误差积分时间内的状态与时间的函数关系;
将所述N组往返阶跃数据分为N条上升阶跃响应曲线和N条下降阶跃响应曲线,其中所述上升阶跃响应曲线的所述参考设定值为X+ΔX,所述下降阶跃响应曲线的所述参考设定值为X-ΔX;
分别计算所述上升阶跃响应曲线的平均响应曲线和所述下降阶跃响应曲线的平均响应曲线;
根据所述上升阶跃平均响应曲线和所述下降阶跃平均响应曲线,分别计算上升阶跃的误差积分和下降阶跃的误差积分,并求两者的平均值作为所述N次往返阶跃的平均误差积分。
7.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S3)中,所述非约束最优化问题的已有求解算法至少包括:最速下降法、梯度下降法、牛顿法、拟牛顿法和人工智能法。
8.根据权利要求1所述的PID参数自动整定方法,其特征在于,步骤S3)中,所述根据所述求解算法更新所述PID参数组合,包括:
获取所述非约束最优化问题的收敛情况;
根据所述求解算法和所述收敛情况计算所述PID参数组合的修正值;
根据所述修正值更新所述PID参数组合。
9.一种PID参数自动整定系统,应用于具有强烈非线性、快速动态响应或对噪声敏感的受控系统,其特征在于,所述系统包括:
传感器单元,用于获取所述受控系统的实时状态;
PID控制器单元,用于根据受控系统特性和控制经验,设置超参数取值和PID参数组合初值;
PID参数整定单元,用于根据所述超参数和所述PID参数组合使受控系统进行若干次往返阶跃,获得每一次阶跃数据并计算平均误差积分;所述PID参数整定单元还被用于判断所述PID参数组合是否满足最优化收敛条件和根据所述求解算法更新所述PID参数组合;
执行单元,用于输出满足最优化收敛条件的PID参数组合。
10.一种计算机可读储存介质,该计算机可读存储介质上储存有指令,其在计算机上运行时使得计算机执行权利要求1至8中任一项权利要求所述的PID参数自动整定方法。
CN202011355027.4A 2020-11-27 2020-11-27 Pid参数自动整定方法及系统 Active CN112835286B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011355027.4A CN112835286B (zh) 2020-11-27 2020-11-27 Pid参数自动整定方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011355027.4A CN112835286B (zh) 2020-11-27 2020-11-27 Pid参数自动整定方法及系统

Publications (2)

Publication Number Publication Date
CN112835286A true CN112835286A (zh) 2021-05-25
CN112835286B CN112835286B (zh) 2023-02-17

Family

ID=75923389

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011355027.4A Active CN112835286B (zh) 2020-11-27 2020-11-27 Pid参数自动整定方法及系统

Country Status (1)

Country Link
CN (1) CN112835286B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115903453A (zh) * 2022-09-30 2023-04-04 荣耀终端有限公司 一种马达的控制方法、控制系统以及电子设备

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013128214A1 (en) * 2012-02-28 2013-09-06 Aristole University Of Thessaloniki-Research Committee A method for auto-tuning of pid controllers and apparatus therefor
US20140081427A1 (en) * 2011-05-27 2014-03-20 Politechnika Swietokrzyska Tuning methods of the pi and pid controllers parameters
CN107728481A (zh) * 2017-11-14 2018-02-23 江西理工大学 一种基于模型预测控制的闭环建模方法及装置
CN108226679A (zh) * 2018-01-10 2018-06-29 云南电网有限责任公司电力科学研究院 一种接力器反应时间常数可调的测试方法及装置
CN108241288A (zh) * 2016-12-23 2018-07-03 重庆邮电大学 一种基于改进tlbo算法的pid控制器参数整定方法
CN109709795A (zh) * 2018-12-24 2019-05-03 东华大学 一种基于天牛须搜索算法的pid控制器参数整定方法
CN110045597A (zh) * 2019-05-30 2019-07-23 青岛科技大学 用于刀具检测机械臂精准工作的改进模糊pid控制方法
CN110377970A (zh) * 2019-06-27 2019-10-25 中国南方电网有限责任公司 弱阻尼低频振荡模式的水轮机调速器参数优化方法及系统
WO2020029217A1 (zh) * 2018-08-10 2020-02-13 曾喆昭 一种智慧pid控制方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140081427A1 (en) * 2011-05-27 2014-03-20 Politechnika Swietokrzyska Tuning methods of the pi and pid controllers parameters
WO2013128214A1 (en) * 2012-02-28 2013-09-06 Aristole University Of Thessaloniki-Research Committee A method for auto-tuning of pid controllers and apparatus therefor
CN108241288A (zh) * 2016-12-23 2018-07-03 重庆邮电大学 一种基于改进tlbo算法的pid控制器参数整定方法
CN107728481A (zh) * 2017-11-14 2018-02-23 江西理工大学 一种基于模型预测控制的闭环建模方法及装置
CN108226679A (zh) * 2018-01-10 2018-06-29 云南电网有限责任公司电力科学研究院 一种接力器反应时间常数可调的测试方法及装置
WO2020029217A1 (zh) * 2018-08-10 2020-02-13 曾喆昭 一种智慧pid控制方法
CN109709795A (zh) * 2018-12-24 2019-05-03 东华大学 一种基于天牛须搜索算法的pid控制器参数整定方法
CN110045597A (zh) * 2019-05-30 2019-07-23 青岛科技大学 用于刀具检测机械臂精准工作的改进模糊pid控制方法
CN110377970A (zh) * 2019-06-27 2019-10-25 中国南方电网有限责任公司 弱阻尼低频振荡模式的水轮机调速器参数优化方法及系统

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
EUGENIO TURCO NETO等: "A sequential optimization procedure to tune MIMO PID controllers", 《2017 INDIAN CONTROL CONFERENCE (ICC)》 *
NAIDU,K等: "Multiobjective optimization using weighted sum Artificial Bee Colony algorithm for Load Frequency Control", 《INTERNATIONAL JOURNAL OF ELECTRICAL POWER & ENERGY SYSTEMS》 *
于兰英等: "基于改进鸡群算法的比例-积分-微分参数整定", 《科学技术与工程》 *
冯新强等: "基于逐维黄金分割法的PID控制器参数优化", 《国网技术学院学报》 *
可寅: "基于粒子群算法的热工过程PID参数优化", 《仪器仪表用户》 *
杨奕飞等: "基于IEM的动力定位PID控制器参数整定", 《电子技术应用》 *
王晓天等: "基于遗传算法和神经网络的PID参数自整定", 《吉林大学学报(理学版)》 *
王欠欠等: "基于遗传算法整定的PID控制", 《中国仪器仪表》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115903453A (zh) * 2022-09-30 2023-04-04 荣耀终端有限公司 一种马达的控制方法、控制系统以及电子设备
CN115903453B (zh) * 2022-09-30 2024-04-09 荣耀终端有限公司 一种马达的控制方法、控制系统以及电子设备

Also Published As

Publication number Publication date
CN112835286B (zh) 2023-02-17

Similar Documents

Publication Publication Date Title
CN108489015B (zh) 基于极点配置和帕德近似的空调系统温度控制方法
CN110103709B (zh) 一种车辆踏板开度的控制方法、系统、设备和介质
CN112835286B (zh) Pid参数自动整定方法及系统
CN106354013B (zh) 攻角的线性自抗扰控制方法
CN110942483B (zh) 函数快速收敛模型构建方法、装置和终端
WO2021213177A1 (zh) 一种船舶分段建造曲面变形动态补偿方法
CN109062040B (zh) 基于系统嵌套优化的预测pid方法
CN113276833B (zh) 车辆的横向运动控制方法、控制终端及存储介质
CN113183957A (zh) 车辆控制方法、装置、设备及自动驾驶车辆
CN111880413B (zh) 一种针对船舶航向保持的自适应动态面算法
Kamenko et al. A novel fuzzy logic scheme for PID controller auto-tuning
CN112757299A (zh) 质心轨迹的获取方法、装置、机器人及存储介质
KR20160119511A (ko) 퍼지 pid 제어기 생성 장치 및 그 생성 방법
CN110262221B (zh) 一种热工过程中对象的pid控制器参数控制方法
Wang et al. Intelligent chaos synchronization of fractional order systems via mean-based slide mode controller
CN109828586B (zh) 基于非线性不确定性的无人船舶鲁棒h∞航向控制方法
US11086277B2 (en) System and method for determining the parameters of a controller
CN111221250B (zh) 一种具有参数不确定性和多个外部扰动的非线性系统及其设计方法
US20230324885A1 (en) Control assistance device, control system, and control assistance method
CN114371618A (zh) 基于神经网络的扩张状态观测器补偿方法及自抗扰控制器
CN111898726B (zh) 一种电动汽车控制系统参数优化方法、设备及存储介质
CN115128948A (zh) Rst史密斯预测器
CN114614490A (zh) 无功电压控制方法和装置、介质以及计算装置
Sadalla et al. Analysis of simple anti-windup compensation in approximate pole-placement control of a second order oscillatory system with time-delay
CN111585292B (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