1.一种基于数据驱动的时滞系统PID控制器镇定方法，其特征在于所述控制方法包括以下步骤：(1) 给被控对象一个阶跃输入信号u，对控制系统输出进行滤波，测量其输出端开始出现响应信号y的时间，该时间为被控对象的时滞，用θ表示；

(2) 建立包括被控时滞对象G(s)、PID控制器C(s)和转换开关的继电反馈系统，C(s)的形式为：其中，kp、ki和kd分别为控制器的比例、积分和微分增益；将继电反馈系统中的开关拨至继电器输入端；

(3) 由继电反馈下的输入和输出数据，利用继电反馈辨识，获得被控对象的多点频率响应数据，并根据频率响应数据绘制被控对象的Bode图和Nyquist图；下面分别考虑被控对象为稳定系统和不稳定系统两种情况：(a)稳定的被控时滞对象

把被控过程在继电反馈下的响应分为稳态部分和暂态部分：u(t)=Δu(t)+us(t)，y(t)=Δy(t)+ys(t)(2)其中，u(t)和y(t)分别为输入和输出响应，us(t)和ys(t)分别为输入和输出响应的稳态部分，Δu(t)和Δy(t)分别为输入和输出响应的暂态部分；设系统在t=Tf时达到稳态的周期振荡，达到振荡后的振荡周期为Tc，则被控对象的频率响应数据为：其中，T是采样周期，Δy(kT)和ys(kT)分别是Δy(t)和ys(t)在t=KT时刻的采样值，Δu(kT)和us(kT)分别是Δu(t)和us(t)在t=KT时刻的采样值，N和Nc是分别满足等式(N-1)T=Tf和Nc=(Tc-T)/T的正整数，k=0,1,2,…,N-1，ωl=2πl/(NT)；被控对象G(s)的频域响应在频域控制理论中可被描述为G(jω)=Gr(ω)+jGi(ω) (4)其中，Gr(ω)和Gi(ω)分别为被控对象的频域响应数据的实部和虚部，从而可绘出被控对象的Bode图和Nyquist图；

(b)不稳定的被控时滞对象

如果被控对象G(s)是不稳定的，将首先通过手动调节给出一个能使G(s)稳定的PID控制器G0(s)，接着根据稳定的被控对象的频率响应数据的获取方法，给出由控制器G0(s)和被控对象G(s)组成的闭环系统的频率响应数据T(jωl)，最后由下式给出被控时滞对象G(s)的频率响应数据并由频率响应数据绘出被控对象的Bode图和Nyquist图；

(4) 基于步骤(3)中所获得的Bode图和Nyquist图，计算确定PID控制器稳定域所必需的特性参数。将G(s)写成如下的形式其中，其中，N(s)和D(s)是关于s的多项式，θ表示时滞；令n和m分别为式(6)中D(s)和N(s)的最高阶次，r(N)、l(N)和j(N)分别为G(s)在右半平面、左半平面和虚轴上的零点，r(D)、l(D)和j(D)分别为G(s)在右半平面、左半平面和虚轴上的极点；由Bode图和Nyquist图给出确定PID稳定域所需的被控对象的特性参数值n-m，r(N)，r(D)，j(N)和j(D)：(a) 根据下式确定n-m的值

其中，Pdb(ω)=20log10|G(jω)|，dPdb(ω)/d(log10ω)表示函数Pdb(ω)关于log10ω的一阶导数；

(b) 计算j(N)和j(D)的值

若Bode图中的幅值曲线有突变，则必存在虚轴零点或极点，若是上升突变，则为虚轴极点，若为下降突变，则为虚轴零点；令突变点的个数为V，所对应的ω值为ωv，其中v=1,2,…,V。接着，绘制 的幅值曲线图，其中，uv=1,2,…，直至幅值曲线图上的突变情况消失。假定突变情况消失时，uv=Uv；对于上升突变的情况，则虚轴极点的个数利用上述计算上升突变情况下的虚轴极点个数的方法，可得到下降突变时的虚轴零点个数j(N)；

(c) 计算r(D)和r(N)的值

(i) 稳定的被控时滞对象

*

由于被控对象是稳定的，因此r(D)=0；令L 为一足够大的正整数，由下式计算r(N)：*

其中， 为取整函数， 为ω由0变化到2 Lπ/θ的相位差；

(ii) 不稳定的被控时滞对象

令能使闭环系统稳定的控制器C0(s)的传递函数的分子和分母多项式分别为NC(s)和DC(s)，nc为多项式DC(s)的最高阶次，l(NC)和r(NC)分别为NC(s)的左半平面零点个数和右半平面零点个数，η为区间[0,π/4]内的实数，包含控制器C0(s)和被控对象G(s)的闭环传递函数T(s)为*

根据T(jω)的Bode图和Nyquist图可计算出T(jω)在ω由(-2Lπ+η)/θ变化*到(2Lπ+η)/θ的相位差 ，从而利用下式获取r(N)的值由

G(jω)的Bode图和Nyquist图可计算出 ，进而根据式(9)获取r(D)的值；

(5)令z=θω，由频域响应数据Gr(ω)和Gi(ω)获得Gr(z)，Gi(z)和 ，其中，Gr(z)和Gi(z)分别为G(jz/θ)的实部和虚部；

(6) 将继电反馈控制系统的开关拨至控制器输入端；

(7) 选取足够大的正整数l*和在区间[0,π/4]内的η值，根据下述条件确定能使闭环系统稳定的kp的最大可允许稳定范围：*

其中，W为函数f1=kp 和 在(0,2lπ+η)间的交点个数；

(8) 令kp的最大可允许稳定范围为[kpmin,kpmax]，将kp值在该范围内进行等间隔的遍历，即每个遍历点为 ，其中F为遍历点之间的间隔，w=0,1,…,F；

(9) 对于其中一个遍历点 ，根据以下步骤确定能够保证闭环系统稳定的(kd,ki)二维稳定域：(a) 根据系统的闭环特征函数得到

令s=jz/θ，可得到 的实部 和 分别为

*

确定 在区间 [0,2lπ+η)的实根，按从小到大的顺序表示为z0,z1,z2,…*,zd-1，其中，z0=0，d为 在区间 [0,2lπ+η)的实数个数；

(b) 令I={i0,i1,…,id-1}，根据下述稳定性的充分必要条件确定I：*

γ(I)=4l+(n-m)+j(N)+2r(N) (16)其中，

，ε是一足够小的正实数，t=0,1,2,…d时所对应的it需满足下述条件之一：

(i)若被控对象在虚轴上有零点zt/θ，则it=0；

(ii)若被 控 对象 在原 点处 有 零点，则 ，其 中，表示函数 关于z的一阶导数；

(iii) 对其它的t=0,1,2,…d-1，it=1或-1；

(c) 对于遍历点 ，(kd,ki)二维稳定域由下式决定:其中，zt为 在区间 [0,2l*π+η)的实根，t=0,1,2,…d-1。通过求取所有zt所对应的由(17)所决定的不等式组的交集，即可确定具有凸多边形特性的(kd,ki)二维稳定域；

(10) 对于步骤(8)中所给出的kp的每个遍历点，都重复步骤(9)，确定能使闭环系统稳定的所有PID控制器集合；

(11) 在所获得的PID控制器稳定集合中选取控制参数，执行被控对象的PID镇定控制：首先对控制系统输出采样滤波，经模拟量输入通道传输信号，并将信号接入检测变送装置，再经A/D转换后得到数字量输入信号与此时的系统设定值比较后得到不同时刻的跟踪误差，基于跟踪误差，按照离散域PID控制算式计算控制信号增量Δu(n)的值，与前一时刻的控制信号u(n-1)通过加法器进行加法运算就得到当前时刻的输出控制信号u(n)，其中，n为当前时刻的采样步数。Δu(n)计算公式如下：Δu(n)=b1e(n)+b2e(n-1)+b3e(n-2) (18)其中，b1=(kpR+kd+R2ki)/R，b2=-(kpR+2kd)/R，b3=kd/R ，R为系统采样周期，Δu(n)为当前采样步数为n时控制器输出信号增量，e(n)为当前采样步数为n时的跟踪误差，e(n-1)为采样步数为n-1时的跟踪误差，e(n-2)为采样步数为n-2时的跟踪误差；输出控制信号u(n)由D/A转换后输出至执行器，由执行器作用到被控对象，使被控对象运行在稳定状态，从而实现被控对象的PID镇定控制。