BP神经网络整定PID

小羊学习记录

已于 2023-02-21 17:24:10 修改

阅读量1w

点赞数 18

文章标签：神经网络

于 2022-07-01 19:19:32 首次发布

本文链接：https://blog.csdn.net/m0_66341377/article/details/125554784

版权

BP神经网络整定PID：

传统增量式数字PID控制算法为：

$u(k) = u(k - 1) + {k_p}[e(k) - e(k - 1)] + {k_i}e(k) + {k_d}[e(k) - 2e(k - 1) + e(k - 2)]$

现引入三层BP神经网络，网络输入层为： ${O_j}^{(1)} = x(j)$ j=1,2,...,M

隐含层的输入输出为： $\begin{array}{l} ne{t_i}^{(2)}(k) = \sum\limits_{j = 0}^M {w_{ij}^{(2}} O_j^{(1)}\\ O_i^{(2)}(k) = f(net_i^{(2)}(k)) \end{array}$ i=1,2,...,Q

其中， $w_{ij}^{(2)}$ 为隐含层加权系数。隐含层是网络的内部信息处理层，负责信息变换。隐含层神经元的活化函数取正负对称的Sigmoid函数：

$f(x) = \tanh (x) = \frac{{{e^x} - {e^{ - x}}}}{{{e^x} + {e^{ - x}}}}$

网络输出层的输入和输出为：

$\begin{array}{l} net_l^{(3)}(k) = \sum\limits_{i = 0}^Q {w_{li}^{(3)}O_i^{(2)}(k)} \\ O_l^{(3)}(k) = g(net_l^{(3)}(k))\\ O_1^{(3)}(k) = {k_p}\\ O_2^{(3)}(k) = {k_i} \end{array}$

其中， $O_l^{(3)}(k)$ 表示神经网络输出层的三个输出节点，输出层向外界输出信息处理结果。这里的三个输出节点分别对应PID控制器的三个可调参数。而输出层神经元的活化函数为非负的Sigmoid函数：

$g(x) = \frac{{{e^x}}}{{{e^x} + {e^{ - x}}}}$

性能指标函数为： $E(k) = \frac{1}{2}{(rin(k) - yout(k))^2}$

整个神经网络按E(k)对加权系数的负梯度方向搜索调整，修正网络的权系数，并附加一个使搜索调整快速收敛于全局极小的惯性项：

$\Delta w_{li}^{(3)}(k) = - \eta \frac{{\partial E(k)}}{{\partial \Delta w_{li}^{(3)}}} + \alpha \Delta w_{li}^{(3)}(k - 1)$ ，式中 $\eta$ 为学习速率， $\alpha$ 为惯性系数。

$\frac{{\partial E(k)}}{{\partial \Delta w_{li}^{(3)}}} = \frac{{\partial E(k)}}{{\partial y(k)}} \bullet \frac{{\partial y(k)}}{{\partial u(k)}} \bullet \frac{{\partial u(k)}}{{\partial O_l^{(3)}(k)}} \bullet \frac{{\partial O_l^{(3)}(k)}}{{\partial net_l^{(3)}(k)}} \bullet \frac{{\partial net_l^{(3)}(k)}}{{\partial w_{li}^{(3)}(k)}}$

其中y(k)为被控对象的传递函数。

且 $\frac{{\partial net_l^{(3)}(k)}}{{\partial w_{li}^{(3)}(k)}} = O_i^{(2)}(k)$ ，由于 $\frac{{\partial y(k)}}{{\partial u(k)}}$ 未知，这里近似用符号函数 ${\mathop{\rm sgn}} (\frac{{\partial y(k)}}{{\partial u(k)}})$ 取代，由此产生的计算不精确的影响通过调整学习速率来补偿。

$\begin{array}{l} \frac{{\partial u(k)}}{{\partial O_1^{(3)}(k)}} = e(k) - e(k - 1)\\ \frac{{\partial u(k)}}{{\partial O_2^{(3)}(k)}} = e(k)\\ \frac{{\partial u(k)}}{{\partial O_3^{(3)}(k)}} = e(k) - 2e(k - 1) + e(k - 2) \end{array}$

推导出网络输出层权的学习算法为：

$\Delta w_{li}^{(3)}(k) = \alpha \Delta w_{li}^{(3)}(k - 1) + \eta \delta _l^{(3)}O_i^{(2)}(k)$

$\delta _{l}^{(3)}=e(k)sgn(\frac{\partial y(k)}{\partial u(k)})\frac{\partial u(k)}{\partial O_{l}^{(3)}(k)}g'(net_{l}^{(3)}(k))$

l=1,2,3...

式中, g'(x)=g(x)(1-g(x))

同理可得隐含层加权系数的学习算法：

$\Delta w_{ij}^{(2)}(k) = \alpha \Delta w_{ij}^{(2)}(k - 1) + \eta \delta _i^{(2)}O_j^{(1)}(k)$

$\delta _{i}^{(2)}=f'(net_{i}^{(2)}(k))\sum_{l=1}^{3}\delta _{l}^{(3)}w_{li}^{(3)}(k)$ i=1,2,3...Q

式中, $f'(x)=(1-f(x)_{}^{2})/2$ 。

S函数：

function [sys,x0,str,ts,simStateCompliance] = sfun_BP(t,x,u,flag,T,j,xite,alfa)
switch flag,
  case 0,
    [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes(T,j);
%初始化函数
  case 3,
    sys=mdlOutputs(t,x,u,xite,alfa);
%输出函数
  case {1,2,4,9},
    sys=[];
  otherwise
    DAStudio.error('Simulink:blocks:unhandledFlag', num2str(flag));
end
function [sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes(T,j)
%调用初始化函数，T步长，j隐含层神经元数
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 0;
sizes.NumOutputs     = 4;
%定义输出变量，包括控制变量u,三个PID参数：Kp,Ki,Kd
sizes.NumInputs      = 8;
%定义输入变量，包括7个参数[e(k);e(k-1);e(k-2);y(k);y(k-1);r(k);u(k-1)]和偏置量u(8)=1
sizes.DirFeedthrough = 1;
sizes.NumSampleTimes = 1; 
sys = simsizes(sizes);
x0  = [];
str = [];
ts  = [T 0];
global wi_2 wi_1 wo_2 wo_1
wi_2 =rand(j,4).*2-1;
%隐含层（k-2)权值系数矩阵，维数j*4，范围【-1，1】
wo_2 = rand(3,j);
%输出层（k-2）权值系数矩阵，维数3*j，范围【0，1】
wi_1 = wi_2;
%隐含层（k-1)权值系数矩阵，维数j*4
wo_1 = wo_2;
%输出层（k-1）权值系数矩阵，维数3*j
simStateCompliance = 'UnknownSimState';
function sys=mdlOutputs(t,x,u,xite,alfa)
%调用输出函数
M=[50;2;25];
%PID权值
global wi_2 wi_1 wo_2 wo_1
xi = [u(6),u(4),u(1),u(8)];
%神经网络的输入xi=[u(6),u(4),u(1),u(8)]=[r(k),y(k),e(k),1] 维数1*4
xx = [u(1)-u(2);u(1);u(1)+u(3)-2*u(2)];
%xx=[u(1)-u(2);u(1);u(1)+u(3)-2*u(2)]=[e(k)-e(k-1);e(k);e(k)+e(k-2)-2*e(k-1)]3*1
I = xi*wi_1';
%计算隐含层的输入，I=神经网络的输入*隐含层权值系数矩阵的转置wi_1'
Oh = (exp(I)-exp(-I))./(exp(I)+exp(-I));
%激活函数，计算隐含层的输出,为1*j的矩阵
O = wo_1*Oh';
%计算输出层的输入，维数3*1
K =exp(O)./(exp(O)+exp(-O));
%激活函数,计算输出层的输出K=[Kp,Ki,Kd]，维数为3*1
K(1)=M(1)*K(1);K(2)=M(2)*K(2);K(3)=M(3)*K(3);
uu = u(7)+K'*xx;
%根据增量式PID控制算法计算控制变量u(k) 维数1*1
if uu>15
    uu=15;
end
if uu<-15
    uu=-15;
end
%限制输出u
dyu = sign((u(4)-u(5))/(uu-u(7)+0.0001));
%计算输出层加权系数修正公式的sgn
%sign((y(k)-y(k-1))/(u(k)-u(k-1)+0.0001）)近似代表偏导 维数1*1
dO = (2./(exp(O)+exp(-O))).^2;
%激活函数，维数3*1
delta3 = u(1)*dyu*xx.*dO;
wo = wo_1+xite*delta3*Oh+alfa*(wo_1-wo_2);
%输出层加权系数矩阵的修正
dI = 2./(exp(I)+exp(-I)).^2;
%激活函数，维数1*j
wi = wi_1+xite*(dI.*(delta3'*wo))'*xi+alfa*(wi_1-wi_2);
%隐含层加权系数修正
wi_2=wi_1;
wi_1=wi;
wo_2=wo_1;
wo_1=wo;
sys = [uu;K(:)];
%输出层输出sys=[uu;Kp;Ki;Kd]

simulink：

得到的阶跃响应为：

其中，

KP:

KI:

KD:

参考原文链接：https://blog.csdn.net/weixin_42650162/article/details/90678503

参考文献

[1]李捷菲. 基于BP神经网络的PID控制系统研究与设计[D].吉林大学,2019.

[2]K. Jiangming and L. Jinhao, "Self-Tuning PID Controller Based on Improved BP Neural Network," 2009 Second International Conference on Intelligent Computation Technology and Automation, 2009, pp. 95-98, doi: 10.1109/ICICTA.2009.32.

[3]Gan Jialiang, Li Zhimin and Tan Huaijiang, "Research on self-tuning PID control strategy based on BP neural network," Proceedings of 2011 International Conference on Electronics and Optoelectronics, 2011, pp. V2-16-V2-21, doi: 10.1109/ICEOE.2011.6013163.