基于BP神经网络的PID参数自适应整定

基于BP 神经网络的PID 参数自适应整定

曾正1，蔡容容2，詹立新2

1 武汉大学电气工程学院，43007

2 2 武汉大学动力与机械学院，430072 联系方式：zengerzheng@

摘 要：针对简单单入单出（SISO ）系统中PID 控制的参数整定问题进行了仿真研究，利用BP 神经网络进行PID 参数自适应整定。

首先，得到了问题的传递函数模型，并建立了对应的离散化传递函数模型作为仿真研究的对象，并对未校正系统进行了相关的理论分析。

然后，利用BP 神经网络算法在线进行PID 参数自适应整定。同时，为了形成参照，给出了运用模拟退火算法离线整定PID 参数的仿真过程。建立了3层BP 网络结构，以δ规则为学习规则，控制器算法为有监督的Hebb 算法，教师信号为给定与被控对象的输出间的偏差信号。

最后，针对系统在多工况下的复杂情况，从静态稳定和动态稳定的角度，对系统稳定性进行了仿真分析。

仿真结果表明，所设计的基于BP 神经网络参数自适应PID 控制系统，控制品质高、鲁棒性强。

为了支持更加复杂情况下的仿真，还搭建GUI 仿真界面。 关键字：PID 控制；参数整定；BP 神经网络；稳定分析；GUI ；

1 问题背景分析 1.1 准备知识

1.1.1 控制系统的性能指标

当系统的时间响应()y t 中的瞬态分量较大而不能忽视时，称系统处于动态或过渡过程中，这时系统的特性称为动态特性。动态特性指标通常根据系统的阶跃响应曲线定义。设系统的阶跃响应曲线如图1所示，图中()lim ()x y y t →∞

∞=称为稳态值。动态性能指标主要有以

下几种[1]

。

图1 系统的阶跃响应曲线

（1）上升时间r t ：阶跃响应曲线从零第一次上升到稳态值所需的时间为上升时间。若阶跃曲线不超过稳态值（称为过阻尼系统），则定义阶跃响应曲线从稳态值的10%上升到90%所对应的时间为上升时间。

（2）最大超调p σ：设阶跃响应曲线的最大值为()p y t ，则对大超调p σ为

()()100%()

p p y t y y σ-∞=

⨯∞ (1)

p σ大，称系统阻尼小。

（3）调节时间s t ：阶跃响应曲线进入并保持在允许的误差范围内所对应的时间称为调节时间。这个误差范围通常为稳态值的∆倍，∆称为误差带，∆为5%或2%。

当系统t →∞时，系统的响应趋于稳定，此时系统跟踪阶跃信号的误差

()lim ()lim(()())ss ss t t e e t r t y t →∞

→∞

∞==- (2)

定义为稳态误差，它反映了系统跟踪阶跃信号的能力，也是最重要静态性能指标。 1.1.2 系统的数学模型

系统的数学模型在经典控制理论中有微分方程、传递函数、动态框图和信号流图，在现代控制理论中广泛采用状态空间模型。

在初始条件为零时，线性定常系统或元件输出信号的拉氏变换()Y s 与输入信号的拉氏变换()R s 之比，称为该系统或元件的传递函数，通常记为()G s 。因此，有

()

()()

Y s G s R s =

(3) 同时，应注意到在已知系统传递函数的情况下，由()()()Y s G s R s =，经过反拉氏变换，可以很容易地知道系统在给定输入信号()r t 下的响应()y t 。

在现代控制理论中，系统的数学模型是用状态空间法来表示的，对于一个简单的单输入单输出系统来说，可以描述为

x Ax Bu

y Cx Du

=+⎧⎨

=+⎩ (4) 其中，y 为输出变量，x 为n 维状态向量。

状态空间模型和传递函数模型是统一的，它们之间可以相互转换。在本问题中，我们将用到将状态空间模型转换为传递函数模型的方法，其逆转换在本文中不涉及，可参考任何一本控制理论相关的书籍。

由状态空间模型转换为传递函数模型的公式为

1()

()[]()

Y s G s C sI A B R s -=

=- (5) p

r s

(y t (y ∞(p y t

1.2 问题分析

对于SISO 系统

12

21225 6.5 6.5x x x x x u

=⎧⎨

=--+⎩ (6) 输入为u ，输出为1y x =，用PID 控制方法进行设计，并用数值仿真确定控制器参数，使系统输出跟踪期望的输出：3ref y =。

由1.1可知，系统的传递函数模型为

1

12

10() 6.5

()[][10]25 6.5 6.5() 6.525

s C s G s C sI A B s R s s s ---⎡⎤⎡⎤==-==⎢⎥⎢⎥+++⎣⎦⎣⎦ (6) 在MA TLAB 的Simulink 工具箱中建立如图2所示的仿真模型，系统跟踪给定阶跃输入信号

的结果如图3所示。

图2 系统在Simulink 下的仿真模型

从图3可以看出：系统是稳定的，但是不能够实现对给定阶跃信号的跟踪。稳态响应为0.78，

模型的稳态误差也可以通过理论计算得到

2

6.53

()lim ()lim ()()lim 0.786.525t s s y y t sG s R s s

s s s

→∞

→→∞====++ (7)

()()30.78 2.22ss ref e y y ∞=-∞=-= (8)

理论的分析结果与仿真结果完全一致。

为了实现对给定阶跃信号的跟踪，必须对系统进行校正或控制。对这样一个本身稳定的系统，可以采用预测控制、模糊控制、自适应控制等先进控制技术，也可以采用传统的PID 控制（或称为校正）技术。问题要求采用PID 控制，且整定合适的PID 参数。我们将采用BP 神经网络自适应的方法给出PID 参数的整定值，这些参数必须是使得控制系统的性能最优的，我们拭目以待。

为了便于计算机数值仿真，我们将模型进行了离散化。采用零阶保持的方法将模型离散化，采样周期为0.01s 。离散化的仿真模型为

-3-3

2

0.31810 z + 0.311210()z - 1.935 z + 0.9371

G z ⨯⨯= (9) 1.3 PID 基本原理

在模拟控制系统中，控制器最常用的控制规律是PID 控制。模拟PID 控制系统原理框图如图4所示。系统由模拟PID 控制器和被控制对象组成。

图4 模拟PID 控制系统原理框图

PID 控制器是一种线性控制器，它根据给定值rin 与实际输出值yout 构成控制偏差

()()ref e t y y t =- (10)

PID 的控制律为

()

1()(()())t

d p i

T de t u t k e t e t dt T dt

=+

+

⎰

(11) 或写成传递函数的形式

()1

()(1)()p d i U s G s k T s E s T s

=

=++ (12) 式中，p k 为比例系数，i T 为积分常数，d T 为微分时间常数。

简单说，PID 控制器各校正环节的作用为

（1）比例环节：成比例地反映控制系统的偏差信号()e t ，偏差一旦产生，控制器立即产生控制作用，减少偏差；

（2）积分环节：主要用于消除静差，提高系统的无差度。积分作用的强弱取决于积分时间常数i T ，i T 越大，积分作用越弱，反之则越强。

（3）微分环节：反映偏差信号的变化趋势（变化速率），并且能在偏差信号变化太大之前，在系统中引入一个有效的早期修正信号，从而加快系统的动作速度，减少调节时间。

2 符号说明

()y t ：系统的输出；

ref y ：系统待跟踪的信号；

()y ∞：系统稳态输出；