CN114488780A - 一种pid控制方法以及pid控制器 - Google Patents

一种pid控制方法以及pid控制器 Download PDF

Info

Publication number
CN114488780A
CN114488780A CN202210338582.9A CN202210338582A CN114488780A CN 114488780 A CN114488780 A CN 114488780A CN 202210338582 A CN202210338582 A CN 202210338582A CN 114488780 A CN114488780 A CN 114488780A
Authority
CN
China
Prior art keywords
state
loop
output feedback
feedback value
pid
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
CN202210338582.9A
Other languages
English (en)
Other versions
CN114488780B (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.)
Zhejiang Supcon Technology Co Ltd
Original Assignee
Zhejiang Supcon Technology Co Ltd
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 Zhejiang Supcon Technology Co Ltd filed Critical Zhejiang Supcon Technology Co Ltd
Publication of CN114488780A publication Critical patent/CN114488780A/zh
Application granted granted Critical
Publication of CN114488780B publication Critical patent/CN114488780B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00—Automatic controllers
    • G05B11/01—Automatic controllers electric
    • G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42—Automatic 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控制方法及PID控制器，根据PID闭环控制回路的输入设定值和滤波输出反馈值计算控制偏差；根据PID闭环控制回路的输出反馈值和滤波输出反馈值计算噪声水平；对滤波输出反馈值的一阶导数和二阶导数进行滤波处理，得到第一输出反馈值和第二输出反馈值；根据控制偏差、第一输出反馈值、第二输出反馈值以及PID参数确定状态特征向量；根据状态特征向量和噪声水平识别PID闭环控制回路的回路状态并基于回路状态对PID参数进行控制。本申请基于状态识别的方法，根据回路控制状态给出PID参数的的调整方向，实现了PID参数的高效和自适应整定，避免人工整定消耗人力和时间，同时对回路初始参数不敏感。

Description

一种PID控制方法以及PID控制器
本申请要求于2022年3月21日提交中国专利局、申请号为202210277390.1、发明名称为“一种PID控制方法以及PID控制器”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。
技术领域
本申请涉及电气自动化控制领域，尤其涉及一种PID控制方法以及PID控制器。
背景技术
比例-积分-微分(Proportional-Integral-Derivative, PID)控制器广泛地应用于各种自动化控制的场合，PID控制器就是根据系统的误差，利用比例、积分、微分计算出控制量进行控制。
现有的PID整定方法分为用户手动整定方式和参数自动整定方式，对于用户手动整定方式，需要依靠专业技术人员的经验，对于整定响应慢的回路，往往比较耗时。对于参数自动整定方式，大多需要给定一个或若干个初始参数，且对于初始参数有一定的精度要求，随着智能控制技术的发展，如何高效地对PID参数进行整定是目前待解决的问题。
发明内容
有鉴于此，本申请提供了一种PID控制方法以及PID控制器，用以实现PID参数的高效整定，其方案如下：
一种PID控制方法，应用于PID闭环控制回路，所述方法包括：
对所述PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值，所述PID闭环控制回路的输出反馈值通过预设采样周期采样得到；
根据所述PID闭环控制回路的输入设定值和所述滤波输出反馈值计算控制偏差；
对所述滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，所述滤波输出反馈值的一阶导数通过对所述滤波输出反馈值求导得到；
对所述滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，所述滤波输出反馈值的二阶导数通过对所述第一输出反馈值求导得到；
根据所述控制偏差、所述第一输出反馈值、所述第二输出反馈值以及PID参数确定所述PID闭环控制回路的状态特征向量，其中，所述PID参数包括所述PID闭环控制回路的比例分量、积分分量以及微分分量；
根据所述PID闭环控制回路的输出反馈值和所述滤波输出反馈值计算噪声水平；
根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态；
根据所述回路状态对所述PID参数进行控制。
优选的，所述根据所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数确定所述PID闭环控制回路的状态特征向量，包括：
根据所述控制偏差，所述第一输出反馈值以及所述第二输出反馈值确定所述PID闭环控制回路所属的特征段；
根据采样时间和所述预设采样周期确定采样点总数；
根据所述特征段，所述采样点总数，所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数构建所述PID闭环控制回路的状态特征向量。
优选的，所述根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态，包括：
构建所述PID闭环控制回路的状态特征矩阵，其中，所述状态特征矩阵包括N个第一状态特征向量，所述第一状态特征向量包括所述状态特征向量，所述N为大于或等于2的正整数；
当从第1个所述第一状态特征向量对应的最大控制偏差到第N个所述第一状态特 征向量对应的最大控制偏差逐渐变小，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE001
，确定所述回路状态为振荡衰 减状态；
当所述N个第一状态特征向量对应的最大控制偏差维持在预设波动范围内，且所 述控制偏差和所述噪声水平满足
Figure 283779DEST_PATH_IMAGE002
与
Figure DEST_PATH_IMAGE003
，确定所述回路 状态为等幅振荡状态；
当从第1个所述第一状态特征对应的最大控制偏差向量到第N个所述第一状态特 征向量对应的最大控制偏差逐渐变大，且所述控制偏差和所述噪声水平满足
Figure 528816DEST_PATH_IMAGE004
，确定所述回路状态为振荡发 散状态；
当第N个所述第一状态特征向量对应的特征段属于第一预设特征段，且所述第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定所述回路状态为发散状态；
当第1个所述第一状态特征向量对应的特征段属于第二预设特征段，且所述采样点总数大于预设采样阈值，确定所述回路状态为调节慢状态；
其中，所述
Figure DEST_PATH_IMAGE005
为预设阈值，所述
Figure 370870DEST_PATH_IMAGE006
为所述噪声水平，所述
Figure DEST_PATH_IMAGE007
为第i个 第一状态特征向量对应的最大控制偏差，所述
Figure 562817DEST_PATH_IMAGE008
为第i+1个第一状态特征向 量对应的最大控制偏差。
优选的，所述根据所述回路状态对所述PID参数进行控制，包括：
根据所述PID参数确定第一调整分量和第二调整分量，所述第一调整分量是所述比例分量、所述积分分量以及所述微分分量中的最大值，所述第二调整分量是所述比例分量以及所述积分分量中的较小值；
当所述回路状态为振荡衰减状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的85%；
当所述回路状态为等幅振荡状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的70%；
当所述回路状态为振荡发散状态或发散状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的50%；
当所述回路状态为调节慢状态，将所述第二调整分量调整至所述第二调整分量的当前值的125%。
优选的，所述PID闭环控制回路的状态特征向量具体为：
Figure DEST_PATH_IMAGE009
，其中，
Figure 675392DEST_PATH_IMAGE010
为大于或等于1的正整数，
Figure DEST_PATH_IMAGE011
表示
Figure 560171DEST_PATH_IMAGE012
时刻对应的特征段，
Figure DEST_PATH_IMAGE013
为当前采样周期起始时 刻，
Figure 155100DEST_PATH_IMAGE012
表示当前时刻，
Figure 150738DEST_PATH_IMAGE014
表示
Figure 852240DEST_PATH_IMAGE013
时刻到
Figure 376763DEST_PATH_IMAGE012
时刻的采样点总数，
Figure DEST_PATH_IMAGE015
表示
Figure 927830DEST_PATH_IMAGE013
时刻到
Figure 196000DEST_PATH_IMAGE012
时刻的比例分量作用面积的大小，
Figure 516123DEST_PATH_IMAGE016
表示
Figure 742705DEST_PATH_IMAGE013
时刻到
Figure 282533DEST_PATH_IMAGE012
时刻的积分分量作用面积的大 小，
Figure DEST_PATH_IMAGE017
表示
Figure 354394DEST_PATH_IMAGE013
时刻到
Figure 794603DEST_PATH_IMAGE012
时刻的微分分量作用面积的大小，
Figure 660927DEST_PATH_IMAGE018
表示
Figure 452166DEST_PATH_IMAGE013
时刻到
Figure 62139DEST_PATH_IMAGE012
时 刻的控制偏差的面积的大小，
Figure DEST_PATH_IMAGE019
表示
Figure 655056DEST_PATH_IMAGE013
时刻到
Figure 957862DEST_PATH_IMAGE012
时刻的最大控制偏差值，
Figure 705238DEST_PATH_IMAGE020
表示
Figure 650060DEST_PATH_IMAGE013
时刻到
Figure 533703DEST_PATH_IMAGE012
时刻的最大第一输出反馈值，
Figure DEST_PATH_IMAGE021
表示
Figure 571191DEST_PATH_IMAGE013
时刻到
Figure 337022DEST_PATH_IMAGE012
时 刻的最大第二输出反馈值，
Figure 819956DEST_PATH_IMAGE022
表示转置运算。
本申请另一实施例还提供了一种PID控制器，所述PID控制器包括：
输出滤波模块，用于对所述PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值，所述PID闭环控制回路的输出反馈值通过预设采样周期采样得到；
计算模块，用于根据所述PID闭环控制回路的输入设定值和所述滤波输出反馈值计算控制偏差；
第一滤波模块，用于对所述滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，所述滤波输出反馈值的一阶导数通过对所述滤波输出反馈值求导得到；
第二滤波模块，用于对所述滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，所述滤波输出反馈值的二阶导数通过对所述第一输出反馈值求导得到；
向量确定模块，用于根据所述控制偏差、所述第一输出反馈值、所述第二输出反馈值以及PID参数确定所述PID闭环控制回路的特征向量，其中，所述PID参数包括所述PID闭环控制回路的比例分量、积分分量以及微分分量；
噪声计算模块，用于根据所述PID闭环控制回路的输出反馈值和所述滤波输出反馈值计算噪声水平；
状态确定模块，用于根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态；
控制模块，用于根据所述回路状态对所述PID参数进行控制。
优选的，所述向量确定模块包括：
特征段确定子模块，用于根据所述控制偏差，所述第一输出反馈值以及所述第二输出反馈值确定所述PID闭环控制回路所属的特征段；
状态特征向量子模块，用于根据所述特征段，采样点总数，所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数构建所述PID闭环控制回路的状态特征向量，其中，所述采样点总数根据采样时间和所述预设采样周期确定。
优选的，所述状态确定模块包括：
矩阵构建子模块，用于构建所述PID闭环控制回路的状态特征矩阵，其中，所述状态特征矩阵包括N个第一状态特征向量，所述第一状态特征向量包括所述状态特征向量，所述N为大于或等于2的正整数；
状态确定子模块，用于当从第1个所述第一状态特征向量对应的最大控制偏差到 第N个所述第一状态特征向量对应的最大控制偏差逐渐变小，且所述控制偏差和所述噪声 水平满足
Figure DEST_PATH_IMAGE023
，确定所述回路状 态为振荡衰减状态；
所述状态确定子模块，还用于当所述N个第一状态特征向量对应的最大控制偏差 维持在预设波动范围内，且所述控制偏差和所述噪声水平满足
Figure 844588DEST_PATH_IMAGE024
与
Figure DEST_PATH_IMAGE025
，确定所述回路状态 为等幅振荡状态；
所述状态确定子模块，还用于当第1个所述第一状态特征向量对应的最大控制偏 差到第N个所述第一状态特征向量对应的最大控制偏差逐渐变大，且所述控制偏差和所述 噪声水平满足
Figure 20355DEST_PATH_IMAGE004
，确定所述回路 状态为振荡发散状态；
所述状态确定子模块，还用于当第N个所述第一状态特征向量对应的特征段属于第一预设特征段，且所述第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定所述回路状态为发散状态；
所述状态确定子模块，还用于当第1个所述第一状态特征向量对应的特征段属于第二预设特征段，且所述采样点总数大于预设采样阈值，确定所述回路状态为调节慢状态；
其中，所述
Figure 7902DEST_PATH_IMAGE005
为预设阈值，所述
Figure 763369DEST_PATH_IMAGE006
为所述噪声水平，所述
Figure 621603DEST_PATH_IMAGE007
为第i个 第一状态特征向量对应的最大控制偏差，所述
Figure 204157DEST_PATH_IMAGE008
为第i+1个第一状态特征向 量对应的最大控制偏差。
优选的，所述控制模块包括：
调整分量确定子模块，用于根据所述PID参数确定第一调整分量和第二调整分量，所述第一调整分量是所述比例分量、所述积分分量以及所述微分分量中的最大值，所述第二调整分量是所述比例分量以及所述积分分量中的较小值；
控制子模块，用于当所述回路状态为振荡衰减状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的85%；
所述控制子模块，还用于当所述回路状态为等幅振荡状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的70%；
所述控制子模块，还用于当所述回路状态为振荡发散状态或发散状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的50%；
所述控制子模块，还用于当所述回路状态为调节慢状态，将所述第二调整分量调整至所述第二调整分量的当前值的125%。
优选的，所述PID闭环控制回路的状态特征向量具体为：
Figure 147842DEST_PATH_IMAGE026
，其中，
Figure 972578DEST_PATH_IMAGE010
为大于或等于1的正整数，
Figure 216478DEST_PATH_IMAGE011
表示
Figure 202888DEST_PATH_IMAGE012
时刻对应的特征段，
Figure 135335DEST_PATH_IMAGE013
为当前采样周期起始时 刻，
Figure 967024DEST_PATH_IMAGE012
表示当前时刻，
Figure 65430DEST_PATH_IMAGE014
表示
Figure 488321DEST_PATH_IMAGE013
时刻到
Figure 203337DEST_PATH_IMAGE012
时刻的采样点总数，
Figure 871341DEST_PATH_IMAGE015
表示
Figure 293095DEST_PATH_IMAGE013
时刻到
Figure 621308DEST_PATH_IMAGE012
时刻的比例分量作用面积的大小，
Figure 761302DEST_PATH_IMAGE016
表示
Figure 731532DEST_PATH_IMAGE013
时刻到
Figure 804530DEST_PATH_IMAGE012
时刻的积分分量作用面积的大 小，
Figure 336268DEST_PATH_IMAGE017
表示
Figure 229138DEST_PATH_IMAGE013
时刻到
Figure 471900DEST_PATH_IMAGE012
时刻的微分分量作用面积的大小，
Figure 133826DEST_PATH_IMAGE018
表示
Figure 803841DEST_PATH_IMAGE013
时刻到
Figure 449586DEST_PATH_IMAGE012
时 刻的控制偏差的面积的大小，
Figure 997505DEST_PATH_IMAGE019
表示
Figure 513937DEST_PATH_IMAGE013
时刻到
Figure 354854DEST_PATH_IMAGE012
时刻的最大控制偏差值，
Figure DEST_PATH_IMAGE027
表示
Figure 487895DEST_PATH_IMAGE013
时刻到
Figure 72460DEST_PATH_IMAGE012
时刻的最大第一输出反馈值，
Figure 708978DEST_PATH_IMAGE028
表示
Figure 487840DEST_PATH_IMAGE013
时刻到
Figure 311440DEST_PATH_IMAGE012
时 刻的最大第二输出反馈值，
Figure 434116DEST_PATH_IMAGE022
表示转置运算。
本申请实施例提供的一种PID控制方法及PID控制器，对PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值；根据所述PID闭环控制回路的输入设定值和所述滤波输出反馈值计算控制偏差；对所述滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，所述滤波输出反馈值的一阶导数通过对所述滤波输出反馈值求导得到；对所述滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，所述滤波输出反馈值的二阶导数通过对所述第一输出反馈值求导得到；根据所述控制偏差、所述第一输出反馈值、所述第二输出反馈值以及PID参数确定所述PID闭环控制回路的状态特征向量，其中，所述PID参数包括所述PID闭环控制回路的比例分量、积分分量以及微分分量；根据所述PID闭环控制回路的输出反馈值和所述滤波输出反馈值计算噪声水平；根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态；根据所述回路状态对所述PID参数进行控制。本申请实施例提供的PID控制方法利用状态识别的方法，根据输出反馈值和控制偏差的特性确定PID闭环回路控制状态，根据回路控制状态的来判断PID参数的合理性，并给出PID参数的的调整方向，实现了PID参数的高效和自适应整定。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
图1为本申请实施例提供的一种PID闭环控制回路示意图；
图2为本申请实施例提供的一种回路状态响应曲线示意图；
图3为本申请实施例提供的一种回路状态响应特征段示意图；
图4为本申请实施例提供的一种PID控制方法的流程示意图；
图5至图8为本申请实施例相关仿真波形示意图；
图9为本申请实施例提供的一种PID控制器的模块示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1所示的一种PID闭环控制回路，图1中
Figure DEST_PATH_IMAGE029
代表被控对象模型，
Figure 925141DEST_PATH_IMAGE030
和
Figure DEST_PATH_IMAGE031
分别为被控对象的输出反馈值和输入值，
Figure 904598DEST_PATH_IMAGE032
为干扰。
被控对象可以表示为
Figure DEST_PATH_IMAGE033
，控制偏差
Figure 779275DEST_PATH_IMAGE034
定义为
Figure DEST_PATH_IMAGE035
； 其中，
Figure 236801DEST_PATH_IMAGE012
为时间，
Figure 316753DEST_PATH_IMAGE036
为设定值。
定义输出反馈值
Figure 670374DEST_PATH_IMAGE030
的一阶导数为：
Figure 734145DEST_PATH_IMAGE038
定义输出反馈 值
Figure 464203DEST_PATH_IMAGE030
的二阶导数为：
Figure 177424DEST_PATH_IMAGE040
其中，
Figure DEST_PATH_IMAGE041
为采样周期。
对于线性被控对象，当回路干扰
Figure 498684DEST_PATH_IMAGE032
较小且PID参数保持不变时，则认为其单位闭 环阶跃响应的特征保持不变，响应特征包括超调量、上升时间、峰值时间、调节时间等，这一 组PID参数对应一个状态，定义PID回路的状态为
Figure 784172DEST_PATH_IMAGE042
，其中
Figure DEST_PATH_IMAGE043
表示PID参数，状态简 写为
Figure 317921DEST_PATH_IMAGE044
。
PID作用包括比例Proportional、积分Integral和微分Derivative 三部分，被控对象的输入值可以表示为：
Figure 372465DEST_PATH_IMAGE046
比例
Figure DEST_PATH_IMAGE047
、积分
Figure 366091DEST_PATH_IMAGE048
和微分
Figure 100002_DEST_PATH_IMAGE049
三部分的作用可以表示为：
比例作用
Figure 935613DEST_PATH_IMAGE050
；
积分作用
Figure 741895DEST_PATH_IMAGE051
；
微分作用
Figure DEST_PATH_IMAGE052
；
此时，被控对象可以进一步表示为
Figure 447683DEST_PATH_IMAGE053
。
其中，
Figure DEST_PATH_IMAGE054
为比例调节系数，
Figure 845166DEST_PATH_IMAGE041
为采样周期，
Figure 606711DEST_PATH_IMAGE055
为积分时间，
Figure DEST_PATH_IMAGE056
为微分时间。
根据测量到的输入值
Figure 13421DEST_PATH_IMAGE031
、
Figure 42557DEST_PATH_IMAGE012
时刻的控制偏差值
Figure 814204DEST_PATH_IMAGE034
、
Figure 561580DEST_PATH_IMAGE057
时刻的控制偏差值
Figure DEST_PATH_IMAGE058
以及当前的PID参数，可以计算得到
Figure 506402DEST_PATH_IMAGE012
时刻的比例作用
Figure 891510DEST_PATH_IMAGE047
和微分作用
Figure 630796DEST_PATH_IMAGE049
， 则积分作用
Figure 865468DEST_PATH_IMAGE059
，当需要对三部分参数进行调整时，可以对各 部分中的参数进行调整，例如，当需要增大比例作用时，可以增大比例调节系数；当需要增 大积分作用时，可以增大比例调节系数或减少积分时间；当需要增大微分作用时，可以增大 比例调节系数或增大微分时间。
图2示出了一种常见的闭环曲线，输出
Figure 551664DEST_PATH_IMAGE030
曲线呈现出振荡衰减的趋势，AC段为 上升段，控制偏差
Figure DEST_PATH_IMAGE060
且逐渐减小，其中AB段为加速段，BC段为减速段。CD段为减速上 升段，CD段的控制偏差
Figure 86551DEST_PATH_IMAGE061
且绝对值逐渐增大，代表着超调。DE段和HI段为下降段，控 制偏差
Figure 200000DEST_PATH_IMAGE061
且
Figure DEST_PATH_IMAGE062
，但是DE段为加速段，HI段为减速段。EF段为减速下降段，且 控制偏差
Figure 954592DEST_PATH_IMAGE060
。FH段为上升段。
如图3所示，将闭环相应曲线划分为8种特征段类型和分界点P，其控制偏差
Figure 975637DEST_PATH_IMAGE034
、 输出反馈值一阶导
Figure 833872DEST_PATH_IMAGE063
和输出反馈值二阶导
Figure DEST_PATH_IMAGE064
的特征对应参见表1，其中分界点P表 示控制偏差
Figure 180540DEST_PATH_IMAGE034
方向变化的分界点，其中图3中的横轴为时间，代表曲线的变化方向，虚线 代表设定值
Figure 858646DEST_PATH_IMAGE036
，是为了体现输出反馈值
Figure 152224DEST_PATH_IMAGE030
对应的特征段与设定值的相对位置。
特征段
Figure DEST_PATH_IMAGE066
Figure DEST_PATH_IMAGE068
Figure DEST_PATH_IMAGE070
 ID
RA1 <0 ≥0 ≥0 1
RA2 >0 ≥0 ≥0 2
RD1 <0 ≥0 ≤0 3
RD2 >0 ≥0 ≤0 4
DA1 <0 ≤0 ≤0 5
DA2 >0 ≤0 ≤0 6
DD1 <0 ≤0 ≥0 7
DD2 >0 ≤0 ≥0 8
P =0 0
表1
为了实现PID参数的高效整定，本申请实施例提供了一种PID控制方法，应用于PID闭环控制回路，所述方法参阅图4，该方法可以包括：
步骤S100：通过预设采样周期采样得到PID闭环控制回路的输出反馈值，对PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值。
可以理解的是，当PID参数或设定值
Figure 491064DEST_PATH_IMAGE036
发生改变时，回路状态
Figure 477474DEST_PATH_IMAGE044
随之发生变化， 此时按照预设采样周期开始识别和采集PID闭环控制回路的输出反馈值。
在实际计算中，由于干扰的影响，回路的输出反馈值
Figure 908455DEST_PATH_IMAGE030
存在毛刺，参见图2中的 DI段以及图5输出反馈值仿真波形所示，如果直接根据回路的输出反馈值确定特征段，会导 致特征段很短，影响后续分析，因此，可以对输出
Figure 5724DEST_PATH_IMAGE030
进行宽度为(2L+1)个采样周期的移 动平均滤波得到
Figure 307393DEST_PATH_IMAGE071
，本申请实施例所提供的滤波方式仅为示例，不因此限定本申请的 保护范围。
具体的，滤波输出反馈值为
Figure DEST_PATH_IMAGE072
。
步骤S200：根据PID闭环控制回路的输入设定值和滤波输出反馈值计算控制偏差。
具体的，控制偏差为
Figure 497328DEST_PATH_IMAGE073
，可以理解的是，由滤波后的输出反馈 值计算得到控制偏差，消除了毛刺，计算得到控制偏差更为稳定，参见图6示出的由滤波后 的输出反馈值计算得到的控制偏差仿真波形图。
步骤S300：对滤波输出反馈值求导得到滤波输出反馈值的一阶导数，对滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值。
具体的，滤波输出反馈值的一阶导数为
Figure DEST_PATH_IMAGE074
，同样对
Figure 477922DEST_PATH_IMAGE063
进行 宽度为(2L+1)个采样周期的移动平均滤波得到第一输出反馈值
Figure 113303DEST_PATH_IMAGE075
，第一输出反馈值为
Figure DEST_PATH_IMAGE076
，参见图7示出的第一输出反馈值的仿真波形图。
步骤S400：对第一输出反馈值求导得到滤波输出反馈值的二阶导数，对滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值。
具体的，滤波输出反馈值的二阶导数为
Figure 66215DEST_PATH_IMAGE077
，同样对
Figure 895893DEST_PATH_IMAGE064
进 行宽度为(2L+1)个采样周期的移动平均滤波得到第二输出反馈值
Figure 770309DEST_PATH_IMAGE064
，第二输出反馈值 为
Figure DEST_PATH_IMAGE078
，参见图8示出的第二输出反馈值的仿真波形图。
步骤S500：根据控制偏差、第一输出反馈值、第二输出反馈值以及PID参数确定PID闭环控制回路的状态特征向量，其中，PID参数包括PID闭环控制回路的比例分量、积分分量以及微分分量。
可选的，步骤S500可以包括步骤S501和步骤S502，具体如下。
S501：根据控制偏差，第一输出反馈值以及第二输出反馈值确定PID闭环控制回路所属的特征段。
可以理解的是，特征段可以根据表1所示的规则进行确定。
S502：根据特征段，采样点总数，控制偏差，第一输出反馈值，第二输出反馈值以及PID参数构建PID闭环控制回路的状态特征向量。
作为一种选择，本申请实施例中状态特征向量可以为：
Figure 537276DEST_PATH_IMAGE079
，其中，
Figure 344695DEST_PATH_IMAGE010
为大于或等于1的正整数，
Figure 843810DEST_PATH_IMAGE011
表示
Figure 471100DEST_PATH_IMAGE012
时刻对应的特征段，
Figure 182704DEST_PATH_IMAGE013
为当前采样周期起始时 刻，
Figure 334376DEST_PATH_IMAGE012
表示当前时刻，
Figure 4392DEST_PATH_IMAGE014
表示
Figure 118978DEST_PATH_IMAGE013
时刻到
Figure 899852DEST_PATH_IMAGE012
时刻的采样点总数，
Figure 681863DEST_PATH_IMAGE015
表示
Figure 522781DEST_PATH_IMAGE013
时刻到
Figure 124663DEST_PATH_IMAGE012
时刻的比例分量作用面积的大小，
Figure 210693DEST_PATH_IMAGE016
表示
Figure 316052DEST_PATH_IMAGE013
时刻到
Figure 327871DEST_PATH_IMAGE012
时刻的积分分量作用面积的大 小，
Figure 417049DEST_PATH_IMAGE017
表示
Figure 805305DEST_PATH_IMAGE013
时刻到
Figure 30750DEST_PATH_IMAGE012
时刻的微分分量作用面积的大小，
Figure 479049DEST_PATH_IMAGE018
表示
Figure 789945DEST_PATH_IMAGE013
时刻到
Figure 217777DEST_PATH_IMAGE012
时 刻的控制偏差的面积的大小，
Figure 297729DEST_PATH_IMAGE019
表示
Figure 916929DEST_PATH_IMAGE013
时刻到
Figure 449542DEST_PATH_IMAGE012
时刻的最大控制偏差值，
Figure 914021DEST_PATH_IMAGE020
表示
Figure 379637DEST_PATH_IMAGE013
时刻到
Figure 435318DEST_PATH_IMAGE012
时刻的最大第一输出反馈值，
Figure 455227DEST_PATH_IMAGE021
表示
Figure 490441DEST_PATH_IMAGE013
时刻到
Figure 279406DEST_PATH_IMAGE012
时 刻的最大第二输出反馈值，
Figure 974829DEST_PATH_IMAGE022
表示转置运算。
可以理解的是，状态特征向量中记录有回路的特征段信息，不同特征段对应不同 的特征向量，各个特征向量可以根据分界点P进行区分，参见表1，状态特征向量
Figure DEST_PATH_IMAGE080
可以为前 一个分界点时刻
Figure 278771DEST_PATH_IMAGE013
到当前时刻
Figure 616212DEST_PATH_IMAGE012
的状态特征向量。
步骤S600：根据PID闭环控制回路的输出反馈值和滤波输出反馈值计算噪声水平。
具体的，可以根据最近L个采样周期来估计噪声水平，定义噪声水平为原输出
Figure 790841DEST_PATH_IMAGE030
与滤波输出
Figure 424210DEST_PATH_IMAGE071
的差值的标准差，即
Figure 418711DEST_PATH_IMAGE081
。
可以理解的是，本申请实施例提供的方法中，不对步骤S600的执行顺序进行限定，步骤S600可以在计算获得滤波输出反馈值之后的任一执行顺序进行。
步骤S700：根据状态特征向量和噪声水平确定PID闭环控制回路的回路状态。
可选的，步骤S700可以包括步骤S701和步骤S702，具体如下。
步骤S701：构建PID闭环控制回路的状态特征矩阵，其中，状态特征矩阵包括N个第一状态特征向量，第一状态特征向量包括状态特征向量，N为大于或等于2的正整数。
具体的，状态特征矩阵包括第一状态特征向量，第一状态特征矩阵包括步骤S502中构建的状态特征向量，通过各个状态特征向量构建PID闭环控制回路的状态特征矩阵，便于通过状态特征矩阵确定回路的状态。
可选的，状态特征矩阵还可以包括第二状态特征向量，第二状态特征矩阵包括PID闭环控制回路在同一状态下的所有数据统计。
作为一种选择，本申请实例第二状态特征矩阵可以为：
Figure 294263DEST_PATH_IMAGE083
其中，
Figure DEST_PATH_IMAGE084
为状态发生变化的时刻，
Figure 588978DEST_PATH_IMAGE012
表示当前时刻，
Figure 626204DEST_PATH_IMAGE043
表示当前状态的PID参数，
Figure 373581DEST_PATH_IMAGE085
表示当前状态下的采样点总数，
Figure DEST_PATH_IMAGE086
表示
Figure 85447DEST_PATH_IMAGE084
时刻到
Figure 969089DEST_PATH_IMAGE012
时刻的比例分量作用面积的大 小，
Figure 177217DEST_PATH_IMAGE087
表示
Figure 880730DEST_PATH_IMAGE084
时刻到
Figure 98085DEST_PATH_IMAGE012
时刻的积分分量作用面积的大小，
Figure DEST_PATH_IMAGE088
表示
Figure 367392DEST_PATH_IMAGE084
时刻到
Figure 12000DEST_PATH_IMAGE012
时刻 的微分分量作用面积的大小，
Figure 235434DEST_PATH_IMAGE089
表示
Figure 990900DEST_PATH_IMAGE084
时刻到
Figure 849135DEST_PATH_IMAGE012
时刻的控制偏差的面积的大小，
Figure DEST_PATH_IMAGE090
表示
Figure 195802DEST_PATH_IMAGE084
时刻到
Figure 139488DEST_PATH_IMAGE012
时刻的最大控制偏差值，
Figure 698645DEST_PATH_IMAGE091
表示
Figure 411386DEST_PATH_IMAGE084
时刻到
Figure 899261DEST_PATH_IMAGE012
时刻的最 大第一输出反馈值，
Figure DEST_PATH_IMAGE092
表示
Figure 595822DEST_PATH_IMAGE084
时刻到
Figure 693091DEST_PATH_IMAGE012
时刻的最大第二输出反馈值，
Figure 525918DEST_PATH_IMAGE022
表示转置 运算。
可选的，状态特征矩阵可以为
Figure 683230DEST_PATH_IMAGE093
，状态
Figure 601507DEST_PATH_IMAGE044
随着PID参数的 改变而发生变化。
可以理解的是，当状态特征矩阵中状态特征向量的数量不足以确定PID闭环控制回路的回路状态时，需要继续记录回路数据直到识别出回路状态。
步骤S702：通过状态特征矩阵的变化可以确定PID闭环控制回路的回路状态，包括：
当从第1个第一状态特征向量对应的最大控制偏差到第N个第一状态特征向量对 应的最大控制偏差逐渐变小，且控制偏差和噪声水平满足
Figure 15651DEST_PATH_IMAGE001
，确定回路状态为振荡衰减状 态；
当N个第一状态特征向量对应的最大控制偏差维持在预设波动范围内，且控制偏 差和噪声水平满足
Figure 702984DEST_PATH_IMAGE002
与
Figure 31197DEST_PATH_IMAGE003
，确定回路状态为 等幅振荡状态；
当从第1个第一状态特征向量对应的最大控制偏差到第N个第一状态特征向量对 应的最大控制偏差逐渐变大，且控制偏差和噪声水平满足
Figure DEST_PATH_IMAGE094
，确定回路状态为振荡发散状 态；
当第N个第一状态特征向量对应的特征段属于第一预设特征段，且第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定回路状态为发散状态；
具体的，第一预设特征段包括RA1特征段或DA2特征段，预设偏差阈值可以根据实际情况设置；
当第1个第一状态特征向量对应的特征段属于第二预设特征段，且采样点总数大于预设采样阈值，确定回路状态为调节慢状态；
具体的，第二预设特征段包括RD2特征段或DD1特征段，预设采样阈值可以根据实际情况设置；
其中，所述
Figure 436770DEST_PATH_IMAGE005
为预设阈值，所述
Figure 875842DEST_PATH_IMAGE006
为所述噪声水平，所述
Figure 417682DEST_PATH_IMAGE007
为第i个 第一状态特征向量对应的最大控制偏差，所述
Figure 447955DEST_PATH_IMAGE008
为第i+1个第一状态特征向 量对应的最大控制偏差。
步骤S800：根据回路状态对PID参数进行控制。
根据前述判断的回路状态对PID参数进行整定，当回路状态为振荡衰减状态、等幅振荡状态、振荡发散状态、发散状态或调节慢状态时，对PID参数进行整定，直到回路状态不再出现控制不佳的状态，即正常状态，结束整定。
具体的，步骤S800包括步骤S801和步骤802，具体如下。
步骤S801：根据PID参数确定第一调整分量和第二调整分量，第一调整分量是比例分量、积分分量以及微分分量中的最大值，第二调整分量是比例分量以及积分分量中的较小值；
具体的，可以根据
Figure 576710DEST_PATH_IMAGE086
，
Figure 553893DEST_PATH_IMAGE087
，
Figure 481398DEST_PATH_IMAGE088
计算比例分量、积分分量以及微分 分量的大小，由此判断各分量之间的关系。
步骤S802：根据回路状态对PID参数进行控制，具体的，包括：
当回路状态为振荡衰减状态，将第一调整分量的值调整至第一调整分量的当前值的85%；
当回路状态为等幅振荡状态，将第一调整分量的值调整至第一调整分量的当前值的70%；
当回路状态为振荡发散状态或发散状态，将第一调整分量的值调整至第一调整分量的当前值的50%；
当回路状态为调节慢状态，将第二调整分量调整至第二调整分量的当前值的125%。
可以理解的是，上述对PID参数的调整数值是本申请实例给到的一种方案，技术人员可以具体更具实际情况进一步做出调整，此处不对分量调整比例做出具体限制。
通过对不同回路状态的识别，针对不同回路状态，当回路处于需整定状态时，如振荡衰减状态、等幅振荡状态、振荡发散状态、发散状态或调节慢状态时，结合当前PID参数的具体情况和当前回路状态给出PID参数的调整方向。
综上所述，本申请实施例提供了一种PID控制方法，基于状态识别的方法，根据输出反馈值和控制偏差的特性确定PID闭环回路控制状态，根据回路控制状态的来判断PID参数的合理性，并给出PID参数的的调整方向，实现了PID参数的高效和自适应整定，避免人工整定消耗人力和时间，同时对回路初始参数不敏感。
与上述方法相对应，本发明实施例还提供了一种PID控制器，请参阅图9，示出了PID控制器的结构示意图，可以包括：
输出滤波模块100，用于对PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值，PID闭环控制回路的输出反馈值通过预设采样周期采样得到；
计算模块200，用于根据PID闭环控制回路的输入设定值和滤波输出反馈值计算控制偏差；
第一滤波模块300，用于对滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，滤波输出反馈值的一阶导数通过对滤波输出反馈值求导得到；
第二滤波模块400，用于对滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，滤波输出反馈值的二阶导数通过对第一输出反馈值求导得到；
向量确定模块500，用于根据控制偏差、第一输出反馈值、第二输出反馈值以及PID参数确定PID闭环控制回路的特征向量，其中，PID参数包括PID闭环控制回路的比例分量、积分分量以及微分分量；
噪声计算模块600，用于根据PID闭环控制回路的输出反馈值和滤波输出反馈值计算噪声水平；
状态确定模块700，用于根据状态特征向量和噪声水平确定PID闭环控制回路的回路状态；
控制模块800，用于根据回路状态对PID参数进行控制。
优选的，向量确定模块500可以进一步包括：
特征段确定子模块501，用于根据控制偏差，第一输出反馈值以及第二输出反馈值确定PID闭环控制回路所属的特征段；
状态特征向量子模块502，用于根据特征段，采样点总数，控制偏差，第一输出反馈值，第二输出反馈值以及PID参数构建PID闭环控制回路的状态特征向量，其中，采样点总数根据采样时间和预设采样周期确定。
优选的，状态确定模块700可以进一步包括：
矩阵构建子模块701，用于构建PID闭环控制回路的状态特征矩阵，其中，状态特征矩阵包括N个第一状态特征向量，第一状态特征向量包括状态特征向量，N为大于或等于2的正整数；
状态确定子模块702，用于当从第1个第一状态特征向量对应的最大控制偏差到第 N个第一状态特征向量对应的最大控制偏差逐渐变小，且控制偏差和噪声水平满足
Figure 885835DEST_PATH_IMAGE001
，确定回路状态为振荡衰减状 态；
状态确定子模块702，还用于当第一状态特征向量对应的最大控制偏差维持在预 设波动范围内，且控制偏差和噪声水平满足
Figure 734842DEST_PATH_IMAGE002
与
Figure 781295DEST_PATH_IMAGE003
，确定回路状态为 等幅振荡状态；
状态确定子模块702，还用于当从第1个第一状态特征向量对应的最大控制偏差到 第N个第一状态特征向量对应的最大控制偏差逐渐变大，且控制偏差和噪声水平满足
Figure 297727DEST_PATH_IMAGE004
，确定回路状态为振荡发散状 态；
状态确定子模块702，还用于当第N个第一状态特征向量对应的特征段属于第一预设特征段，且第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定回路状态为发散状态；
具体的，第一预设特征段包括RA1特征段或DA2特征段，预设偏差阈值可以根据实际情况设置；
状态确定子模块702，还用于当第1个第一状态特征向量对应的特征段属于第二预设特征段，且采样点总数大于预设采样阈值，确定回路状态为调节慢状态；
具体的，第二预设特征段包括RD2特征段或DD1特征段，预设采样阈值可以根据实际情况设置；
其中，所述
Figure 404224DEST_PATH_IMAGE005
为预设阈值，所述
Figure 507571DEST_PATH_IMAGE006
为所述噪声水平，所述
Figure 92136DEST_PATH_IMAGE007
为第i个 第一状态特征向量对应的最大控制偏差，所述
Figure 197495DEST_PATH_IMAGE008
为第i+1个第一状态特征向 量对应的最大控制偏差。
优选的，控制模块800进一步可以包括：
调整分量确定子模块801，用于根据PID参数确定第一调整分量和第二调整分量，第一调整分量是比例分量、积分分量以及微分分量中的最大值，第二调整分量是比例分量以及积分分量中的较小值；
控制子模块802，用于当回路状态为振荡衰减状态，将第一调整分量的值调整至第一调整分量的当前值的85%；
控制子模块803，还用于当回路状态为等幅振荡状态，将第一调整分量的值调整至第一调整分量的当前值的70%；
控制子模块803，还用于当回路状态为振荡发散状态或发散状态，将第一调整分量的值调整至第一调整分量的当前值的50%；
控制子模块803，还用于当回路状态为调节慢状态，将第二调整分量调整至第二调整分量的当前值的125%。
优选的，PID状态特征向量可以为：
Figure 209314DEST_PATH_IMAGE095
，其中，
Figure 32913DEST_PATH_IMAGE010
为大于或等于1的正整数，
Figure 155590DEST_PATH_IMAGE011
表示
Figure 912193DEST_PATH_IMAGE012
时刻对应的特征段，
Figure 861957DEST_PATH_IMAGE013
为当前采样周期起始时 刻，
Figure 172853DEST_PATH_IMAGE012
表示当前时刻，
Figure 99221DEST_PATH_IMAGE014
表示
Figure 444751DEST_PATH_IMAGE013
时刻到
Figure 798372DEST_PATH_IMAGE012
时刻的采样点总数，
Figure 330985DEST_PATH_IMAGE015
表示
Figure 61043DEST_PATH_IMAGE013
时刻到
Figure 261080DEST_PATH_IMAGE012
时刻的比例分量作用面积的大小，
Figure 552647DEST_PATH_IMAGE016
表示
Figure 103714DEST_PATH_IMAGE013
时刻到
Figure 637463DEST_PATH_IMAGE012
时刻的积分分量作用面积的大 小，
Figure 692007DEST_PATH_IMAGE017
表示
Figure 387431DEST_PATH_IMAGE013
时刻到
Figure 894635DEST_PATH_IMAGE012
时刻的微分分量作用面积的大小，
Figure 966496DEST_PATH_IMAGE018
表示
Figure 141126DEST_PATH_IMAGE013
时刻到
Figure 40074DEST_PATH_IMAGE012
时 刻的控制偏差的面积的大小，
Figure 300154DEST_PATH_IMAGE019
表示
Figure 175706DEST_PATH_IMAGE013
时刻到
Figure 673684DEST_PATH_IMAGE012
时刻的最大控制偏差值，
Figure 976489DEST_PATH_IMAGE020
表示
Figure 458286DEST_PATH_IMAGE013
时刻到
Figure 606370DEST_PATH_IMAGE012
时刻的最大第一输出反馈值，
Figure 257057DEST_PATH_IMAGE021
表示
Figure 996343DEST_PATH_IMAGE013
时刻到
Figure 965436DEST_PATH_IMAGE012
时 刻的最大第二输出反馈值，
Figure 182791DEST_PATH_IMAGE022
表示转置运算。
需要说明的是，本申请实施例提供的一种PID控制器中各模块所执行的步骤以及相关技术特征与申请实施例所提供方法相对应，装置部分的描述可以参见前述方法部分的实施例，此处不赘述。
综上所述，本申请实施例提供了一种PID控制器，能够基于状态识别的方法，根据输出反馈值和控制偏差的特性确定PID闭环回路控制状态，根据回路控制状态的来判断PID参数的合理性，并给出PID参数的的调整方向，实现了PID参数的高效和自适应整定，避免人工整定消耗人力和时间，同时对回路初始参数不敏感。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。
本领域技术人员可以理解，图所示的流程图仅是本申请的实施方式可以在其中得以实现的一个示例，本申请实施方式的适用范围不受到该流程图任何方面的限制。
在本申请所提供的几个实施例中，应该理解到，所揭露的方法、装置和设备，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
所述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机，服务器，或者网络设备等）执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器（ROM，Read-Only Memory）、随机存取存储器（RAM，Random Access Memory）、磁碟或者光盘等各种可以存储程序代码的介质。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种PID控制方法，其特征在于，所述方法应用于PID闭环控制回路，所述方法包括：
对所述PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值，所述PID闭环控制回路的输出反馈值通过预设采样周期采样得到；
根据所述PID闭环控制回路的输入设定值和所述滤波输出反馈值计算控制偏差；
对所述滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，所述滤波输出反馈值的一阶导数通过对所述滤波输出反馈值求导得到；
对所述滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，所述滤波输出反馈值的二阶导数通过对所述第一输出反馈值求导得到；
根据所述控制偏差、所述第一输出反馈值、所述第二输出反馈值以及PID参数确定所述PID闭环控制回路的状态特征向量，其中，所述PID参数包括所述PID闭环控制回路的比例分量、积分分量以及微分分量；
根据所述PID闭环控制回路的输出反馈值和所述滤波输出反馈值计算噪声水平；
根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态；
根据所述回路状态对所述PID参数进行控制。
2.根据权利要求1所述的方法，其特征在于，所述根据所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数确定所述PID闭环控制回路的状态特征向量，包括：
根据所述控制偏差，所述第一输出反馈值以及所述第二输出反馈值确定所述PID闭环控制回路所属的特征段；
根据采样时间和所述预设采样周期确定采样点总数；
根据所述特征段，所述采样点总数，所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数构建所述PID闭环控制回路的状态特征向量。
3.根据权利要求2所述的方法，其特征在于，所述根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态，包括：
构建所述PID闭环控制回路的状态特征矩阵，其中，所述状态特征矩阵包括N个第一状态特征向量，所述第一状态特征向量包括所述状态特征向量，所述N为大于或等于2的正整数；
当从第1个所述第一状态特征向量对应的最大控制偏差到第N个所述第一状态特征向量对应的最大控制偏差逐渐变小，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE002
，确定所述回路状态为振荡衰减状态；
当所述N个第一状态特征向量对应的最大控制偏差维持在预设波动范围内，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE004
与
Figure DEST_PATH_IMAGE006
，确定所述回路状态为等幅振荡状态；
当从第1个所述第一状态特征对应的最大控制偏差向量到第N个所述第一状态特征向量对应的最大控制偏差逐渐变大，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE008
，确定所述回路状态为振荡发散状态；
当第N个所述第一状态特征向量对应的特征段属于第一预设特征段，且所述第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定所述回路状态为发散状态；
当第1个所述第一状态特征向量对应的特征段属于第二预设特征段，且所述采样点总数大于预设采样阈值，确定所述回路状态为调节慢状态；
其中，所述
Figure DEST_PATH_IMAGE010
为预设阈值，所述
Figure DEST_PATH_IMAGE012
为所述噪声水平，所述
Figure DEST_PATH_IMAGE014
为第i个第一状态特征向量对应的最大控制偏差，所述
Figure DEST_PATH_IMAGE016
为第i+1个第一状态特征向量对应的最大控制偏差。
4.根据权利要求3所述的方法，其特征在于，所述根据所述回路状态对所述PID参数进行控制，包括：
根据所述PID参数确定第一调整分量和第二调整分量，所述第一调整分量是所述比例分量、所述积分分量以及所述微分分量中的最大值，所述第二调整分量是所述比例分量以及所述积分分量中的较小值；
当所述回路状态为振荡衰减状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的85%；
当所述回路状态为等幅振荡状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的70%；
当所述回路状态为振荡发散状态或发散状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的50%；
当所述回路状态为调节慢状态，将所述第二调整分量调整至所述第二调整分量的当前值的125%。
5.根据权利要求2所述的方法，其特征在于，所述PID闭环控制回路的状态特征向量具体为：
Figure DEST_PATH_IMAGE018
，其中，
Figure DEST_PATH_IMAGE020
为大于或等于1的正整数，
Figure DEST_PATH_IMAGE022
表示
Figure DEST_PATH_IMAGE024
时刻对应的特征段，
Figure DEST_PATH_IMAGE026
为当前采样周期起始时刻，
Figure 14431DEST_PATH_IMAGE024
表示当前时刻，
Figure DEST_PATH_IMAGE028
表示
Figure 790626DEST_PATH_IMAGE026
时刻到
Figure 101522DEST_PATH_IMAGE024
时刻的采样点总数，
Figure DEST_PATH_IMAGE030
表示
Figure 90206DEST_PATH_IMAGE026
时刻到
Figure 438667DEST_PATH_IMAGE024
时刻的比例分量作用面积的大小，
Figure DEST_PATH_IMAGE032
表示
Figure 854605DEST_PATH_IMAGE026
时刻到
Figure 652796DEST_PATH_IMAGE024
时刻的积分分量作用面积的大小，
Figure DEST_PATH_IMAGE034
表示
Figure 445172DEST_PATH_IMAGE026
时刻到
Figure 176367DEST_PATH_IMAGE024
时刻的微分分量作用面积的大小，
Figure DEST_PATH_IMAGE036
表示
Figure 294365DEST_PATH_IMAGE026
时刻到
Figure 579853DEST_PATH_IMAGE024
时刻的控制偏差的面积的大小，
Figure DEST_PATH_IMAGE038
表示
Figure 178849DEST_PATH_IMAGE026
时刻到
Figure 967813DEST_PATH_IMAGE024
时刻的最大控制偏差值，
Figure DEST_PATH_IMAGE040
表示
Figure 256712DEST_PATH_IMAGE026
时刻到
Figure 763917DEST_PATH_IMAGE024
时刻的最大第一输出反馈值，
Figure DEST_PATH_IMAGE042
表示
Figure 632516DEST_PATH_IMAGE026
时刻到
Figure 869462DEST_PATH_IMAGE024
时刻的最大第二输出反馈值，表示转置运算。
6.一种PID控制器，其特征在于，所述PID控制器包括：
输出滤波模块，用于对所述PID闭环控制回路的输出反馈值进行滤波处理，以得到滤波输出反馈值，所述PID闭环控制回路的输出反馈值通过预设采样周期采样得到；
计算模块，用于根据所述PID闭环控制回路的输入设定值和所述滤波输出反馈值计算控制偏差；
第一滤波模块，用于对所述滤波输出反馈值的一阶导数进行滤波处理，以得到第一输出反馈值，所述滤波输出反馈值的一阶导数通过对所述滤波输出反馈值求导得到；
第二滤波模块，用于对所述滤波输出反馈值的二阶导数进行滤波处理，以得到第二输出反馈值，所述滤波输出反馈值的二阶导数通过对所述第一输出反馈值求导得到；
向量确定模块，用于根据所述控制偏差、所述第一输出反馈值、所述第二输出反馈值以及PID参数确定所述PID闭环控制回路的特征向量，其中，所述PID参数包括所述PID闭环控制回路的比例分量、积分分量以及微分分量；
噪声计算模块，用于根据所述PID闭环控制回路的输出反馈值和所述滤波输出反馈值计算噪声水平；
状态确定模块，用于根据所述状态特征向量和所述噪声水平确定所述PID闭环控制回路的回路状态；
控制模块，用于根据所述回路状态对所述PID参数进行控制。
7.根据权利要求6所述的PID控制器，其特征在于，所述向量确定模块包括：
特征段确定子模块，用于根据所述控制偏差，所述第一输出反馈值以及所述第二输出反馈值确定所述PID闭环控制回路所属的特征段；
状态特征向量子模块，用于根据所述特征段，采样点总数，所述控制偏差，所述第一输出反馈值，所述第二输出反馈值以及所述PID参数构建所述PID闭环控制回路的状态特征向量，其中，所述采样点总数根据采样时间和所述预设采样周期确定。
8.根据权利要求7所述的PID控制器，其特征在于，所述状态确定模块包括：
矩阵构建子模块，用于构建所述PID闭环控制回路的状态特征矩阵，其中，所述状态特征矩阵包括N个第一状态特征向量，所述第一状态特征向量包括所述状态特征向量，所述N为大于或等于2的正整数；
状态确定子模块，用于当从第1个所述第一状态特征向量对应的最大控制偏差到第N个所述第一状态特征向量对应的最大控制偏差逐渐变小，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE044
，确定所述回路状态为振荡衰减状态；
所述状态确定子模块，还用于当所述N个第一状态特征向量对应的最大控制偏差维持在预设波动范围内，且所述控制偏差和所述噪声水平满足
Figure DEST_PATH_IMAGE046
与
Figure DEST_PATH_IMAGE048
，确定所述回路状态为等幅振荡状态；
所述状态确定子模块，还用于当第1个所述第一状态特征向量对应的最大控制偏差到第N个所述第一状态特征向量对应的最大控制偏差逐渐变大，且所述控制偏差和所述噪声水平满足
Figure 382790DEST_PATH_IMAGE008
，确定所述回路状态为振荡发散状态；
所述状态确定子模块，还用于当第N个所述第一状态特征向量对应的特征段属于第一预设特征段，且所述第一状态特征向量对应的最大控制偏差大于预设偏差阈值，确定所述回路状态为发散状态；
所述状态确定子模块，还用于当第1个所述第一状态特征向量对应的特征段属于第二预设特征段，且所述采样点总数大于预设采样阈值，确定所述回路状态为调节慢状态；
其中，所述
Figure 439608DEST_PATH_IMAGE010
为预设阈值，所述
Figure 49581DEST_PATH_IMAGE012
为所述噪声水平，所述
Figure 813137DEST_PATH_IMAGE014
为第i个第一状态特征向量对应的最大控制偏差，所述
Figure 115943DEST_PATH_IMAGE016
为第i+1个第一状态特征向量对应的最大控制偏差。
9.根据权利要求8所述的PID控制器，其特征在于，所述控制模块包括：
调整分量确定子模块，用于根据所述PID参数确定第一调整分量和第二调整分量，所述第一调整分量是所述比例分量、所述积分分量以及所述微分分量中的最大值，所述第二调整分量是所述比例分量以及所述积分分量中的较小值；
控制子模块，用于当所述回路状态为振荡衰减状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的85%；
所述控制子模块，还用于当所述回路状态为等幅振荡状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的70%；
所述控制子模块，还用于当所述回路状态为振荡发散状态或发散状态，将所述第一调整分量的值调整至所述第一调整分量的当前值的50%；
所述控制子模块，还用于当所述回路状态为调节慢状态，将所述第二调整分量调整至所述第二调整分量的当前值的125%。
10.根据权利要求7所述的PID控制器，其特征在于，所述PID闭环控制回路的状态特征向量具体为：
Figure DEST_PATH_IMAGE049
，其中，
Figure 191215DEST_PATH_IMAGE020
为大于或等于1的正整数，
Figure 136037DEST_PATH_IMAGE022
表示
Figure 816417DEST_PATH_IMAGE024
时刻对应的特征段，
Figure 293054DEST_PATH_IMAGE026
为当前采样周期起始时刻，
Figure 996568DEST_PATH_IMAGE024
表示当前时刻，
Figure 948343DEST_PATH_IMAGE028
表示
Figure 748809DEST_PATH_IMAGE026
时刻到
Figure 127838DEST_PATH_IMAGE024
时刻的采样点总数，
Figure 380964DEST_PATH_IMAGE030
表示
Figure 870852DEST_PATH_IMAGE026
时刻到
Figure 260245DEST_PATH_IMAGE024
时刻的比例分量作用面积的大小，
Figure 606912DEST_PATH_IMAGE032
表示
Figure 81756DEST_PATH_IMAGE026
时刻到
Figure 375334DEST_PATH_IMAGE024
时刻的积分分量作用面积的大小，
Figure 887742DEST_PATH_IMAGE034
表示
Figure 608574DEST_PATH_IMAGE026
时刻到
Figure 773976DEST_PATH_IMAGE024
时刻的微分分量作用面积的大小，
Figure 871245DEST_PATH_IMAGE036
表示
Figure 500809DEST_PATH_IMAGE026
时刻到
Figure 126963DEST_PATH_IMAGE024
时刻的控制偏差的面积的大小，
Figure 45240DEST_PATH_IMAGE038
表示
Figure 742938DEST_PATH_IMAGE026
时刻到
Figure 899113DEST_PATH_IMAGE024
时刻的最大控制偏差值，
Figure 555222DEST_PATH_IMAGE040
表示
Figure 229304DEST_PATH_IMAGE026
时刻到
Figure 668376DEST_PATH_IMAGE024
时刻的最大第一输出反馈值，
Figure 210216DEST_PATH_IMAGE042
表示
Figure 506068DEST_PATH_IMAGE026
时刻到
Figure 867779DEST_PATH_IMAGE024
时刻的最大第二输出反馈值，
Figure DEST_PATH_IMAGE051
表示转置运算。
CN202210338582.9A 2022-03-21 2022-04-01 一种pid控制方法以及pid控制器 Active CN114488780B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2022102773901 2022-03-21
CN202210277390 2022-03-21

Publications (2)

Publication Number Publication Date
CN114488780A true CN114488780A (zh) 2022-05-13
CN114488780B CN114488780B (zh) 2022-07-26

Family

ID=81488924

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210338582.9A Active CN114488780B (zh) 2022-03-21 2022-04-01 一种pid控制方法以及pid控制器

Country Status (1)

Country Link
CN (1) CN114488780B (zh)

Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5990104A (ja) * 1982-11-15 1984-05-24 Mitsubishi Electric Corp オ−トチユ−ナ
JPH07210205A (ja) * 1994-01-18 1995-08-11 Mitsubishi Heavy Ind Ltd Ｐｉｄ調節器のパラメータ調整装置
US5488560A (en) * 1993-06-30 1996-01-30 West Instruments Limited Apparatus for and method of controlling a process
TW201217924A (en) * 2010-10-22 2012-05-01 Chan Li Machinery Co Ltd for enabling effective control of PID control loop for non-continuous (batch type) process control during start-up process
CN106774243A (zh) * 2016-12-02 2017-05-31 浙江中控软件技术有限公司 Pid性能评估方法
CN106950826A (zh) * 2017-01-22 2017-07-14 航天长征化学工程股份有限公司 一种基于蚁群算法的气化装置pid参数整定方法
CN108873702A (zh) * 2018-07-17 2018-11-23 太原科技大学 一种电液位置伺服控制系统的线性自抗扰控制方法及装置
CN108958021A (zh) * 2018-07-07 2018-12-07 江苏迅达电磁线有限公司 一种工业机器人的直线电机pid控制器的设计方法
CN109001973A (zh) * 2018-07-05 2018-12-14 广东工业大学 一种积分分离式pid控制方法
CN109375499A (zh) * 2018-08-29 2019-02-22 江苏坚诺机电测试技术有限公司 Pid闭环控制系统及其控制方法
CN109856959A (zh) * 2019-02-28 2019-06-07 滨州学院 一种基于pid的直线电机控制方法及系统
US10516403B1 (en) * 2019-02-27 2019-12-24 Ciena Corporation High-order phase tracking loop with segmented proportional and integral controls
CN110764417A (zh) * 2019-11-13 2020-02-07 东南大学 一种基于闭环辨识模型的线性二次型最优动态前馈-反馈pid控制系统及其控制方法
CN111812967A (zh) * 2020-05-27 2020-10-23 浙江中控技术股份有限公司 基于稳定裕度和动态响应指标的pid控制参数整定方法
CN113110029A (zh) * 2021-04-16 2021-07-13 北京黑蚁兄弟科技有限公司 一种基于混合滤波的pid控制方法和装置、一种pid控制设备
CN113138552A (zh) * 2021-03-31 2021-07-20 国网浙江省电力有限公司电力科学研究院 基于阶跃响应数据和临界比例度法的pid参数整定方法
CN113189862A (zh) * 2021-04-29 2021-07-30 中国大唐集团科学技术研究院有限公司中南电力试验研究院 一种火电厂pid控制回路品质评价方法
CN113359413A (zh) * 2021-05-26 2021-09-07 广东电网有限责任公司 一种实际微分器的参数调节系统及方法
CN113433819A (zh) * 2021-06-09 2021-09-24 浙江中控技术股份有限公司 一种筛选pid控制回路中数据的方法和系统辨识方法

Patent Citations (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5990104A (ja) * 1982-11-15 1984-05-24 Mitsubishi Electric Corp オ−トチユ−ナ
US5488560A (en) * 1993-06-30 1996-01-30 West Instruments Limited Apparatus for and method of controlling a process
JPH07210205A (ja) * 1994-01-18 1995-08-11 Mitsubishi Heavy Ind Ltd Ｐｉｄ調節器のパラメータ調整装置
TW201217924A (en) * 2010-10-22 2012-05-01 Chan Li Machinery Co Ltd for enabling effective control of PID control loop for non-continuous (batch type) process control during start-up process
CN106774243A (zh) * 2016-12-02 2017-05-31 浙江中控软件技术有限公司 Pid性能评估方法
CN106950826A (zh) * 2017-01-22 2017-07-14 航天长征化学工程股份有限公司 一种基于蚁群算法的气化装置pid参数整定方法
CN109001973A (zh) * 2018-07-05 2018-12-14 广东工业大学 一种积分分离式pid控制方法
CN108958021A (zh) * 2018-07-07 2018-12-07 江苏迅达电磁线有限公司 一种工业机器人的直线电机pid控制器的设计方法
CN108873702A (zh) * 2018-07-17 2018-11-23 太原科技大学 一种电液位置伺服控制系统的线性自抗扰控制方法及装置
CN109375499A (zh) * 2018-08-29 2019-02-22 江苏坚诺机电测试技术有限公司 Pid闭环控制系统及其控制方法
US10516403B1 (en) * 2019-02-27 2019-12-24 Ciena Corporation High-order phase tracking loop with segmented proportional and integral controls
CN109856959A (zh) * 2019-02-28 2019-06-07 滨州学院 一种基于pid的直线电机控制方法及系统
CN110764417A (zh) * 2019-11-13 2020-02-07 东南大学 一种基于闭环辨识模型的线性二次型最优动态前馈-反馈pid控制系统及其控制方法
CN111812967A (zh) * 2020-05-27 2020-10-23 浙江中控技术股份有限公司 基于稳定裕度和动态响应指标的pid控制参数整定方法
CN113138552A (zh) * 2021-03-31 2021-07-20 国网浙江省电力有限公司电力科学研究院 基于阶跃响应数据和临界比例度法的pid参数整定方法
CN113110029A (zh) * 2021-04-16 2021-07-13 北京黑蚁兄弟科技有限公司 一种基于混合滤波的pid控制方法和装置、一种pid控制设备
CN113189862A (zh) * 2021-04-29 2021-07-30 中国大唐集团科学技术研究院有限公司中南电力试验研究院 一种火电厂pid控制回路品质评价方法
CN113359413A (zh) * 2021-05-26 2021-09-07 广东电网有限责任公司 一种实际微分器的参数调节系统及方法
CN113433819A (zh) * 2021-06-09 2021-09-24 浙江中控技术股份有限公司 一种筛选pid控制回路中数据的方法和系统辨识方法

Also Published As

Publication number Publication date
CN114488780B (zh) 2022-07-26

Similar Documents

Publication Publication Date Title
CN111766777B (zh) 一种pid控制器及pid控制方法
CN108489015B (zh) 基于极点配置和帕德近似的空调系统温度控制方法
CN104360823A (zh) 动态调整固态硬盘带宽的方法及其固态硬盘
CN114488780B (zh) 一种pid控制方法以及pid控制器
EP1850474B1 (en) Method and system for amplifier thermal management
EP3575893A1 (en) Control parameter computation method and control parameter computation device
CN114523966A (zh) 车速控制方法、装置、电子设备以及可读存储介质
CN1854627A (zh) 一种变风量空调系统变静压与总风量双重控制方法
CN112003588A (zh) 基于多态可变步长归一化均方的自适应信号滤波方法
CN114074661A (zh) 一种车速控制方法及装置
CN113110029B (zh) 一种基于混合滤波的pid控制方法和装置、一种pid控制设备
CN114243690A (zh) 一种电网有功安全校正方法、装置、电子设备及存储介质
CN109695893B (zh) 一种锅炉系统中的氧气浓度控制方法、装置、设备及系统
Freiling et al. Monotonicity and convexity properties of matrix Riccati equations
CN109634220B (zh) 一种六自由度机器人运动控制方法及系统
CN107145154B (zh) 一种用于控制飞行器姿态角的控制系统及其控制方法
Kim et al. Controller design for time domain specifications
CN111058952A (zh) 基于自适应模糊技术的电子节气门控制方法、系统及车辆
EP0422809A2 (en) Adaptive apparatus
CN109026711A (zh) 螺杆压缩机无级能量调节自动控制方法、装置及介质
JPH01213701A (ja) 自己増殖型制御装置および同型制御方法，ならびに同制御装置で使用される増殖型コントローラ，その動作方法，その制御方法およびスーパーバイザ
CN110409987A (zh) 一种感光式百叶窗自适应控制方法及系统
CN115016250B (zh) 一种改进型抗饱和结构的伺服速度控制方法
CN112810395B (zh) 基于模糊指令和用户喜好的车载空调智能控制系统及方法
JPH11506553A (ja) ダイナミックプロセスのモデル化装置

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