一种延时系统中的PID控制器参数整定方法

技术领域

本发明涉及扰动估计及抑制领域，具体涉及一种延时系统中的PID控制器参数整定方法，进一步提升系统上升时间、调节时间等动态响应性能指标以及系统的抗干扰能力。

背景技术

在工程控制应用中，系统经常受到外部扰动和自身内部不确定性影响，这些扰动严重影响系统稳定性能和控制效果，甚至可能造成闭环系统不稳定。以惯性稳定系统为例，文献[3](Tian J,Yang W,Peng Z,et al.Application of MEMS accelerometers andgyroscopes in fast steering mirror control systems[J].Sensors,2016,16(4):440.)中采用光纤陀螺仪、加速度计、高分辨率位置探测器组成了多回路反馈控制系统，系统总扰动抑制能力是各回路作用的叠加，但这种方法需要在惯性稳定平台上安装额外的惯性传感器，不利于实现惯性稳定平台小惯性和快速性的要求，同时也增加了实验空间和经济成本。文献[4](Glück M,Pott J U,Sawodny O.Piezo-actuated vibrationdisturbance mirror for investigating accelerometer-based tip-tiltreconstruction in large telescopes[J].IFAC-PapersOnLine,2016,49(21):361-366.)中对于由基座传感器测量到的外部振动，采用基于测量的直接前馈方法进行抑制，但需要准确识别从基座到倾斜镜的扰动传递特性。文献[5](Tang T,Niu S,Chen X,etal.Disturbance observer-based control of tip-tilt mirror for mitigatingtelescope vibrations[J].IEEE Transactions on Instrumentation and Measurement,2018,68(8):2785-2791.)将DOB引入到惯性稳定系统来增强系统抗干扰能力。然而，被控对象模型的特性往往不能被准确识别，补偿器的设计受到限制。

为了进一步提高系统的扰动抑制能力，一种可以用来观测外部扰动和系统内部不确定性的扩张状态观测器被广泛使用。传统的线性扩张状态观测器设计方法将外部扰动和系统内部不确定归为总扰动，并将其扩张为一个状态变量进行观测和补偿。这种设计方法下的系统将被简化为双积分串联标准型，降低了控制的难度。但是，在传统扩张状态观测器的设计中与已知模型信息没有关系，这会使观测器的观测精度的下降，导致系统抗干扰能力下降。本发明对传统的扩张状态观测器设计进行改进，使观测器的观测精度上升。

随着经典控制理论的不断发展，线性二次型最优控制在现代控制理论中有着广泛的应用。线性二次型最优控制方法得到的最优控制律具有许多优良的特性，包括闭环稳定性以及系统过程为单输入单输出时，线性二次型最优控制下的系统相位裕度至少为60°，系统幅值裕度为无穷大。此外，线性二次型最优控制方法中通过选择加权矩阵Q和R，可以控制状态调节要求和控制能量消耗之间的权衡。这种优良的性质促使控制设计师将其用于PID控制器参数的整定。本发明中采用线性二次型最优控制与主导极点配置技术相结合的方法对带有改进后的扩张状态观测器的延时系统进行PID控制器参数进行整定。

发明内容

本发明目的在于：由于工程控制应用中系统延时因素对系统稳定性的影响不可忽略，针对延时系统中存在的建模误差以及不确定扰动问题，本发明提供了一种延时系统中的PID控制器参数整定方法。

本发明采用的技术方案为：一种延时系统中的PID控制器参数整定方法，具体实施步骤如下：

针对二阶模型系统G(s)：

其中，a

步骤1：建立系统微分方程：

其中，y,u为系统输出和输入，w为外部扰动，

将微分方程转化为扩张状态空间方程形式，系统扩张状态空间方程如下：

其中，

步骤2：根据线性系统理论中状态观测器的设计，连续扩展状态观测器ESO如下所示：

其中，z＝[z

步骤3：扩张状态观测器可在一定频率范围内准确估计不确定扰动f′，并对扩展状态z

u(t-L)＝u

将公式(1.49)带入公式(1.46)中，可以得到：

/>

将

将公式(1.51)改写为传递函数为：

通过与公式(1.45)比较，公式(1.51)中的系统传递函数与公式(1.45)中的模型一致。这是本发明中对线性扩张状态观测器的改进。传统的线性扩张状态观测器设计方法与已知模型信息无关，系统将被简化为双积分串联标准型。

步骤4：采用线性二次型最优控制(LQR)以及主导极点技术相结合的方法来整定图1控制框架中PID控制器参数，实现对扰动补偿后的系统控制。图1中u

u

其中，

在外部扰动w(t)以及系统参考信号r(t)作用下，被控对象位置输出信号为y(t)，假设参考信号r(t)＝0，则跟踪误差e(t)＝r(t)-y(t)＝-y(t)，在此条件下，被控系统公式(1.52)可以表示为：

根据状态空间的形式，被控系统公式(1.52)中状态变量的导数可以写成：

其中，

A,B,X,L分别是状态转移矩阵、控制矩阵、状态矩阵和时滞项。从公式(1.55)中可以看出：当t＝L时，控制信号有效。因此，将公式(1.55)分成两部分：如下：

公式(1.58)中u

u

通过公式(1.59)中的变化，即u

为了使公式(1.56)系统具备LQR的性能，需要将如下的二次型代价函数最小化：

其中，Q为半正定状态权矩阵，R为正定控制权矩阵。标准LQR方法给出了最优控制向量u

u

其中，P为对称正定Riccati系数矩阵，可通过求解如下连续代数Riccati方程得到：

A

从公式(1.59)中得到：

u

公式(1.63)给出了u

其中A

A

由于公式(1.66)中得到的系统矩阵不包含任何时间延迟，因此直接应用极点配置的方法来得到期望的闭环时间性能。将公式(1.61)带入到公式(1.58)，可以得到：

然后通过建立闭环系统的特征方程△(s)＝|sI-A

当A

其中

当A

其中非主导极点p

步骤5：在线性二次型最优控制中，标准做法是通过改变加权矩阵Q和保持加权矩阵R不变来进一步设计控制器参数。

假设，

通过公式(1.56)、(1.66)以及(1.70)，可以得到相应的闭环系统特性方程如下：

△(s)＝|sI-A

其中，η＝r

通过对比公式(1.69)与公式(1.73)右边相同状态变量前的系数，可以得到：

矩阵P和矩阵Q的剩余元素可通过求解Riccati方程式(1.62)得到，如下所示：

步骤6：为了得到最优控制信号u

其中，p

通过使用部分分式法，f

步骤7：为了得到最优控制信号u

其中，γ＝ηp

的计算如下：

通过使用部分分式法，f

其中，

p

步骤8：联立公式(1.79)、(1.76)以及(1.64)，计算0≤t

对比公式(1.53)与公式(1.81)，可以得到PID控制器参数：

步骤9：联立公式(1.79)以及(1.65)，计算t≥L的PID参数。

对比公式(1.53)与公式(1.83)，可以得到PID控制器参数：

其中，K

步骤10：设计扩张状态观测器增益矩阵。假设误差状态变量为e(t)＝x(t)-z(t)，由公式(1.48)减去公式(1.47)，则观测器误差矩阵方程为：

从上式可以看出，观测器误差矩阵方程中，(A

|sI-(A

根据相关文献[2](Herbst G.A simulative study on active disturbancerejection control(ADRC)as a control tool for practitioners[J].Electronics,2013,2(3):246-279.)结论，扩张状态观测器经过参数化，可将对应的特征方程的极点放到同一位置上(-w

|sI-(A

将公式(1.86)与公式(1.87)右边相同变量前的系数进行比较，可以得到：

根据以上技术方案，可以实现以下有益效果：

1.与未带有改进扩张状态观测器的方法进行相比，本发明设计方法下系统的抗干扰能力得到进一步提升，即本发明中的改进扩张状态观测器方法有效。

2.与其他时域优化方法(性能指标函数为平方误差积分(ISE)、时间平方误差积分(ITSE)、绝对误差积分(IAE)、时间绝对误差积分(ITAE))相比，本发明设计方法下系统上升时间、调节时间等动态响应性能指标以及系统抗干扰能力得到显著改善。

附图说明

图1为本发明中基于线性二次型最优控制以及改造扩张状态观测器的控制框图；

图2为阶跃扰动作用下系统在不同方法下的阶跃响应对比图；

图3为阶跃扰动作用下的本发明方法下的带有改进扩张状态观测器的系统以及未带有改进扩张状态观测器的系统的阶跃响应对比图。

具体实施方式

以下结合附图和非最小相位系统为例，对本发明的具体实施步骤做详细说明：

本发明一种延时系统中的PID控制器参数整定方法，包括如下步骤：

步骤1：针对非最小相位系统：

将上述非最小相位被控对象转化为扩张状态空间方程形式，系统扩张状态空间方程如下：

其中，

步骤2：根据系统扩张状态空间方程设计改造后的线性扩张状态观测器：

其中，z＝[z

根据公式(1.88)，改造后的扩张状态观测器的增益矩阵L为：

步骤3：在改造后的扩张状态观测器基础上，设计PID控制器，实现对扰动补偿后的系统控制。如图1所示，采取LQR方法以及主导极点配置技术实现对PID控制器参数的整定。通常，加权矩阵R＝1。取系统期望阻尼比ζ

通过使用公式(1.72)，(1.73)，LQR方法整定的加权矩阵P、Q参数如下：

Q＝diag(5.7158 1.4889 9.829) (1.89)

公式(1.89)，(1.90)的特征值为：

eig(P)＝[1.405 3.6579 29.3282]

从公式(1.91)可以看出：加权矩阵P、Q的特征值大于零，满足LQR方法的使用条件。

步骤4：0≤tL的PID控制器的参数可以通过公式(1.84)计算得到。考虑到0≤tL的恒定的PID控制器参数，这样的优势为便于实际实施操作，需要较少的控制力。

t>L的PID控制器的参数为：

[K

步骤5：将本发明中的方法与其他时域优化方法(性能指标函数为平方误差积分(ISE)、时间平方误差积分(ITSE)、绝对误差积分(IAE)、时间绝对误差积分(ITAE))进行比较。仿真中使用MATLAB优化工具箱的fmincon()函数来寻找基于性能指标函数的给定时域性能指标下的优化PID控制器参数集，所有优化方法都是从相同的PID参数初始值0.3开始的，即

[K

将本发明中的方法作用于带有改进扩张状态观测器的系统以及未带有改进扩张状态观测器的系统进行比较，仿真中扩张状态观测器的带宽取为40HZ。

图2为在40s时，给定阶跃幅值为1的扰动作用下的系统在不同方法下的阶跃响应对比图。从图2中可以看出，与其他时域优化方法(性能指标函数为平方误差积分(ISE)、时间平方误差积分(ITSE)、绝对误差积分(IAE)、时间绝对误差积分(ITAE))相比，本发明设计方法下系统上升时间、调节时间等动态响应性能指标以及系统抗干扰能力得到显著改善。

图3为在40s时，给定阶跃幅值为1的扰动作用下的本发明方法下的带有改进扩张状态观测器的系统以及未带有改进扩张状态观测器的系统的阶跃响应对比图。从图3中可以看出，本发明设计方法下系统抗干扰能力得到进一步提升。

以上结合附图和实例对本发明的具体实施方式、过程和效果做详细说明，但所述内容仅为本方法的一个实施例，不能只限定该方法的实施范围。