1.镇定单输入单输出多时滞系统的PID控制器设计方法，其特征在于所述控制方法包括以下步骤：(1)先利用工控系统对以输送到流浆箱中的低浓纸浆为输入，以制造出的成品纸厚度为输出的造纸生产过程进行辨识，给出具有如下形式的多时滞模型G(s)：令G(s)的分子和分母分别为V(s)和U(s)；式(1)中，u和v分别表示U(s)和V(s)中i is项的最高阶次，u>v，τij和αij是U(s)中s 项所对应的时滞和系数，hi是U(s)中s 项i i所对应的时滞块的个数，θij和βij是V(s)中s 项所对应的时滞和系数，fi是V(s)中s项所对应的时滞块的个数，βv为V(s)的最高阶次项的系数。然后，将辨识出的模型参数送到主机的存储单元RAM中；

(2)建立前馈回路中具有PID控制器C(s)和多时滞被控模型G(s)的单位反馈系统，其中，C(s)为具有以下形式的PID控制器：这里，kp、ki和kd分别为控制器的比例、积分和微分增益；

(3)确定系统的闭环特征函数δ(s)为

δ(s)=sU(s)+V(s)(kds2+kps+ki) (3)令ω为频率，Ur(ω)和Ui(ω)为U(jω)的实部和虚部，Vr(ω)和Vi(ω)为V(jω)的实部和虚部，将闭环特征函数乘以V(-s)且令s=jω，得到δ'(jω)=p(ω,ki,kd)+jq(ω,kp)其中,

p(ω,ki,kd)=p1(ω)+(ki-ω2kd)[Vr2(ω)+Vi2(ω)]+jq(ω,kp) (4)q(ω,kp)=ω{q1(ω)+kp[Vr2(ω)+Vi2(ω)]} (5)p(ω,ki,kd)和q(ω,kp)中的p1(ω)和q1(ω)分别为：p1(ω)=ω[Ur(ω)Vi(ω)-Ui(ω)Vr(ω)] (6)q1(ω)=Ur(ω)Vr(ω)+Ui(ω)Vi(ω) (7)(4)选取一个足够大的频率值ω*；

(5)令Re[V(jω)]和Im[V(jω)]分别为V(jω)的实部与虚部，根据下式计算ω由0变化到ω*时所对应的V(jω)的幅角变化范围*

这里，当v为偶数时，w0

(6)确定kp的最大可允许稳定范围：

2 2 *

令Q表示f1(ω)=kp与f2(ω)=-q1(ω)/[Vr(ω)+Vi(ω)]在(0,ω)上的交点数，给出满足下式的kp范围，即为kp的最大可允许稳定范围，*

其中，j(V)表示V(s)在正虚轴上的零点数， 为ω由0变化到ω 时所对应的V(jω)的幅角变化范围；

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

(8)对于其中一个遍历点 根据以下步骤确定能够保证闭环系统稳定的(kd,ki)二维稳定域：*

(a)计算方程q(ω,kp)=0在区间[0,ω)内的实零点，将这些实零点按照从小到大的顺序表示为ω0,ω1,ω2,…,ωc-1，其中，c为实零点的个数；

(b)令t=0,1,2,…,c-1，选择整数it的值：(i).如果V(-jωt)=0，那么it=0；

(ii).如果V(-s)在原点处存在零点，那么

(iii).若不满足(i)和(ii)，it=1或-1，其准确值需根据下述稳定性条件而定；

(iv).令I={i0,i1,…}，确定能够满足下述等式的所有I：γ(I)由下式给出：

其中， 这里ε是一足够小的正实数；

(c)对于遍历点 (kd,ki)二维稳定域由下式决定:[ki-A(ωt)kd+B(ωt)]it>0 (12)*

其中，ωt为式(5)在区间[0,ω)内的实零点，t=0,1,…,c-1， 以及

2 2

B(ωt)=p1(ωt)/[Vr(ωt)+Vi(ωt)]，通过求取所有ωt所对应的由式(12)所决定的不等式组的交集，即可确定具有凸多边形特性的(kd,ki)二维稳定域；

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

(10)通过工控系统将辨识出的模型参数输送到PID控制器参数的计算单元，由步骤(2)-(9)计算出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) (13)2

其中，b1=(kpR+kd+Rki)/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转换后经模拟量输出通道送到调节造纸过程中的烘缸的进气蒸汽压力的进气蒸汽阀，通过控制蒸汽压力来调节纸张的水分，使得系统在达到稳态时，纸张中的水分能够稳定在水分的设定值，以稳定成纸纤维分布的均匀度，提高纸张质量，实现在线控制。