CN101346675B - 自动调整用于空间分布过程的大规模多元模型预测控制器 - Google Patents
自动调整用于空间分布过程的大规模多元模型预测控制器 Download PDFInfo
- Publication number
- CN101346675B CN101346675B CN2006800488517A CN200680048851A CN101346675B CN 101346675 B CN101346675 B CN 101346675B CN 2006800488517 A CN2006800488517 A CN 2006800488517A CN 200680048851 A CN200680048851 A CN 200680048851A CN 101346675 B CN101346675 B CN 101346675B
- Authority
- CN
- China
- Prior art keywords
- model
- parameter
- predictive controller
- model predictive
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- D—TEXTILES; PAPER
- D21—PAPER-MAKING; PRODUCTION OF CELLULOSE
- D21G—CALENDERS; ACCESSORIES FOR PAPER-MAKING MACHINES
- D21G9/00—Other accessories for paper-making machines
- D21G9/0009—Paper-making control systems
- D21G9/0045—Paper-making control systems controlling the calendering or finishing
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
- Paper (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
用于多阵列造纸机横向(CD)过程的大规模多元模型预测控制器的自动调整方法可以显著改善相对于传统控制器的控制器的性能。造纸机CD过程是大规模空间分布动态系统。由于这些系统的(几乎全部)空间变量性质,由具有矩形轮换矩阵块的转移矩阵近似闭环转移函数,转移矩阵的输入和输出奇异向量是维数与致动器的数量或测量的数量相等的傅立叶分量。该近似使用于这些系统的模型预测控制器能够通过数字搜索最优化权重来进行调整,以便在用于性能和鲁棒性的二维频率域中使闭环转移函数具有适当形状。一种新颖的定标方法用于在空间频率域中对多元系统的输入和输出进行定标。
Description
技术领域
本发明通常涉及用于监控和控制连续造纸(continuoussheetmaking)系统的技术,并特别是涉及自动调整用于例如造纸机横向过程的空间分布过程的大规模多元模型预测控制器。
背景技术
在利用现代高速机器造纸的技术中,必须连续监控和控制纸性质以保证纸张质量并在制造过程出现扰动时使不合格的成品量最小化。最经常要测量的纸张变量包括制造过程中不同阶段时纸张的基本重量、水分含量和纸厚度测径器(即厚度)。通常通过例如在过程开始时调节给料供给速度、在过程中期附近调节施加给纸张的蒸汽量、或在过程结束时改变压光辊之间的夹压力,来控制这些过程变量。
造纸机采用跨越连续移动卷筒展开的大型致动器阵列来控制如由从致动器的下游的一个(或数个)扫描传感器测量的纸性质的横向(CD)分布。通常,设计者正使用配对规则来选择用于控制一个纸张性质的一个CD致动器阵列,而多阵列CD过程的相互作用在传统CD控制中通常被忽略。
不幸的是大多数良好设计的单阵列CD系统是病态的(ill-conditioned)。即使在稳态时,一些它们的奇异值也小得难以察觉。大维数和病态使这些过程向控制提出挑战。近来已经提出对于多阵列CD过程而言过程的病态可能是由于多阵列测量和致动器之间的相互作用。这意味着多阵列CD系统比单阵列CD系统可能更难以控制。
已经对CD过程中的模型预测控制(MPC)的应用考虑过一段时间。虽然大多数刊发论文仅考虑一个致动器阵列和一个受控性质,从而不注意控制多个纸张性质的并列多CD致动器阵列的问题,但是多阵列CD控制系统变得更加普遍。工业模型预测控制实施可以采用从补充工业模型识别工具(complementary industrial modelidentification tool)获得的CD过程的多阵列模型。多阵列控制的优点在已经报道过的改善的性能方面是明显的。在线最优化的主要缺点是所需的巨大的计算负担,如约束二次方程式编程(QP)问题可能需要产生如600个那么多的致动器选点,从频繁如每15秒上至6000次测量上经受上至1800个约束。最优化问题是高结构化的,而开发这种结构的最优化算法已经被开发出来。可能的补充技术是使用模型缩小技术来减少最优化问题的大小。
在R.Shridhar和D.J.Cooper的“A tuning strategy forunconstrained multivariable model predictive control”,Industrial & Engineering Chemistry & Research,vol.37,no.10,pp 4003-4016,1998,和D.Dougherty和D.J.Cooper的“Tuningguidelines of a dynamic matrix controller for integrating(non-self-regulating)processes”,Industrial & EngineeringChemistry & Research,vol.42,no.8,pp 1739-1752,2003中,作者提出用于多元动态矩阵控制器的一些调整指导方针。在K.Y.Rani和H.Unbehauen的“Study of predictive controller tuningmethods”,Automatica,vol.33,no.12,pp 2243-2248,1997中,作者提出基于调整规则和闭环防真的用于预测控制器的调整程序。
在J.H.Lee和Z.Yu的“Tuning of model predictivecontrollers for robust performance”,Computers & ChemicalEngineering,vol.18,no.1,pp.15-37,1994中,提出基于MPC控制器的闭环性质的频率域分析的调整规则。在A.A1-Ghazzawi等的“On-line tuning strategy for model predictive controller”,Journal of Process Control,vol.11,no.3,pp 265-284,2001中,基于闭环预测的输出和MPC调整参数之间的线性近似提出用于线性模型预测控制算法的在线调整策略。J.Trierweiler和L.A.Farina的“RPN tuning strategy for model predictive control”,Journal of Process Control,vol.13,no.7,pp.591-598,2003,基于用于多输入多输出(MIMO)MPC的鲁棒性能数量(robustperformance number)提出一种调整策略。不过,上述调整策略不可直接用于调整大规模二维工业CD MPC,尤其是在空间域中。
用于实施造纸机CD MPC控制系统的程序在图1中示出。选择预测水平和最优化权重的第三步通常是ad hoc,并通常经由闭环系统的仿真来估计。该技术的状态是迅速成为对于大规模MPC问题无法解决的问题的试误法。一阶静态最优化器已经被提出,用于在准确预测稳态性能的同时大大减少计算时间。对于CD过程而言,已经看出可以在空间域和时间域中分别调整CD控制器。不过,即使利用静态最优化器,通过试误法调整多元CD预测控制器也是非常困难的,即使对于有经验的工程师也是如此,这是因为存在如此多的(通常大于10)调整参数。另一个实际问题是未指出对于不可避免的模型不确定性来说控制器是否是鲁棒稳定的。
发明内容
本发明针对用于大规模多元MPC的新颖自动调整方法。该方法特别适于应用于造纸机或机外超级压光机中的空间分布的多CD阵列过程。
单阵列CD过程可以由矩形轮换矩阵近似,并接着转换到空间和时间频率域中。本发明部分基于认识到多CD阵列过程可以由矩形轮换矩阵块近似。在适当转换和置换之后,可以将大维数的对象模型(plant model)分离成跨越空间频率的小MIMO系统的族。将在二维频率域中进行调整。过程模型的定标保证调整的参数不随测量和促动的物理单位变化而改变。自动调整的多阵列CD过程在二维频率域中提供所需的闭环行为。
在一个实施例中,本发明针对用于空间分布过程的多元模型预测控制器(MPC)的自动调整方法,包括下列步骤:
(a)识别用于空间分布过程的过程模型;
(b)对过程模型的输入和输出进行定标;
(c)将过程模型转换到二维频率域中;
(d)将MPC转换到二维频率域中;
(e)确定空间调整参数关于模型不确定性的最佳解;和
(f)确定时间调整参数关于模型不确定性的最佳解。
步骤(e)和(f)的每一个都可以采用(i)非结构性模型不确定性,例如乘法(multiplicative)输入不确定性、乘法输出不确定性、逆加法模型不确定性、逆乘法输入不确定性、逆乘法输出不确定性,或(ii)结构性模型不确定性,来满足鲁棒稳定条件。
在另一个实施例中,本发明针对一种过程控制系统,其具有多元模型预测控制器(MPC),用于提供对具有至少一个被操纵致动器阵列和至少一个受控测量阵列的多阵列横向(CD)过程的控制,一种用于提供过程控制的方法,包括下列步骤:
(a)通过下列步骤自动调整MPC:
(i)识别用于CD过程的过程模型;
(ii)对过程模型的输入和输出进行定标;
(iii)将过程模型转换到二维频率域中;
(iv)将MPC转换到二维频率域中;
(v)计算用于MPC的空间和时间调整参数的可能范围;
(vi)确定空间调整参数关于模型不确定性的最佳解;和
(vii)确定时间调整参数关于模型不确定性的最佳解;
(b)将最佳调整参数输入到MPC;和
(c)利用MPC对多阵列CD过程进行控制。
在另一个实施例中,本发明针对在多元阵列横向过程中形成材料的系统,其中该系统包括:
至少两组致动器阵列,每个在横向(CD)上都毗连材料分布,其中每组致动器阵列都是可控的以改变材料的性质;
用于测量和获得关于材料在横向上的性质的性质数据的装置;和
用于为多阵列横向过程提供CD控制的多元模型预测控制器(MPC),其中MPC包括用于自动调整MPC的装置,其中MPC响应于指示性质数据的信号,给至少两组致动器阵列提供信号以改变材料的性质。
附图说明
图1是用于实施工业造纸机CD控制的程序;
图2是造纸厂中使用的机外(off-machine)超级压光机过程。
图3是具有方程14中的无约束MPC的闭环系统框图;
图4是空间频率域中sub子对象模型(subplant model)(z)的对角线性质;
图6a是CD MPC的调整程序;
图6b是CD MPC系统的自动调整程序;
图7(a)和7(b)描绘出空间频率域中的灵敏度函数(a)和控制灵敏度函数(b)的图形(仅示出一个测量和一个致动器阵列);
图8(a)和8(b)描绘出时间频率域中的灵敏度函数(a)和控制灵敏度函数(b)的图形(仅示出一个测量和一个致动器阵列);
图9是空间频率域中的定标静态模型的奇异值的图形;
图10(a)-(f)是空间频率域中稳态时的((a)、(b)和(c))和时间频率域中零空间频率时的((d)、(e)和(f))灵敏度函数的图形;
图11(a)-(c)是通过使用不同控制策略测量分布的2-σ趋势图的曲线。
具体实施方式
本发明的自动调整的方法将通过把该技术实施在控制造纸机的高度压光过程的大规模多元MPC中来说明。高度压光是造纸业中的公共过程。高度压光过程是机上连续过程或机外分批过程。利用它来改善纸张的表面性质,例如对于高质量印刷纸张关键的平滑和光泽(发亮)。纸张可以是无涂层的、单面涂层的或双面涂层的。
应该明白本发明的技术足以灵活到可应用于受到多元MPC控制的任何大规模工业多元横向(CD)过程。例如,除了机外和机上高度压光过程之外,其它造纸机过程也可以通过多元MPC进行控制。由潮湿原料连续造纸的适当造纸机过程在下列多个专利中被进一步介绍,例如,美国专利MacHattie等人的6,805,899,Heaven等人的6,466,839,Hu等人的6,149,770,Hagart-Alexander等人的6,092,003,Heaven等人的6,080,278,Hu等人的6,059,931,Hu等人的6,853,543和He的5,892,679,它们都转让给了Honeywell International,Inc.并通过参考并入在此。本发明技术还可以应用于塑料纸制造(plasticsheetmaking)、橡胶纸制造(rubber sheetmaking)和板金操作(sheetmetal operation)。
超级压光机的多CD阵列过程模型
图2示出造纸厂中采用的机外超级压光机过程。将超级压光机CD过程用作典型的多CD阵列过程模型来说明本发明。也就是,象其它造纸机CD过程那样,超级压光过程是二维(空间和时间)过程。过程开始于在开卷机10处由卷轴开卷纸张。接着在垂直布置的一系列辊12、14、16、18、20、22、24、26、28和30之间装入纸张。纸张经过的两个辊之间的点被称为辊隙。将辊布置成利用堆垛中间的两个连续的软辊18、20来交替改变硬和软。包含两个连续软辊的辊隙被称为换向辊隙。就是在这点处由硬辊进行的工作从纸张的一侧转移到另一侧。纸张经过该堆垛底部出来,通过测量纸张性质的扫描器40,并接着在重绕架42处的卷轴上卷起。用于造纸制造的在线扫描传感器在美国专利Dahlquist的4,879,471、Dablquist等人的5,094,535和Dahlquist的5,166,748中公开,它们全部转让给了Honeywell International,Inc.并通过参考并入在此。可选地,为了CD测量,传感器阵列可以沿着CD定位在毗连移动的纸张的任何适当位置处。
利用热水加热辊16和22。该过程还配备有两个蒸汽盒致动器阵列48和50以及两个感应加热致动器阵列44和46。对于两个蒸汽盒致动器阵列而言,纸正面(top side)的那个指定为u1而纸反面(wire side)的那个指定为u4。对于两个感应加热致动器阵列而言,纸正面的那个指定为u2而纸反面的那个指定为u3。每个蒸汽盒致动器阵列都有n1=n4=31个致动器而每个感应加热致动器阵列都有n2=n3=64个致动器。被控制的纸张性质是由y1指示的纸正面光泽(平滑)、由y2指示的纸厚度测径器(厚度)和由y3指示的纸反面光泽。在扫描器的传感器的信号处理之后,测量阵列的数量为m=192。当该示例性超级压光机具有维数192的3个测量和4个致动器阵列时,本发明的调整方法可应用于具有其它仪器配置的超级压光机,并且实际上,该方法通常可应用于其它多元CD过程。
多CD阵列系统过程模型由下式给出:
Y(z)=G(z)U(z)+D(z),(1)
Gij(z)=Pijhij(z), (3)
其中
Y(z)=[y1(z),y2(z),y3(z)]T∈C576×1,y1(z)∈C192×1,
U(z)=[u1(z),u2(z),u3(z),u4(z),]T∈C190×1,
D(z)=[d1(z),d2(z),d3(z)]T∈C576×1,di(z)∈C192×1,
这里Y(z)、U(z)和D(z)分别是扰动分布、致动器分布和测量分布的z-变换;Pij是描述(i,j)th过程的空间响应的(i,j)th相互作用矩阵;hij(z)是(i,j)th过程的时间响应的z-变换,整数Tdij表示(i,j)th过程死时间,而aij与(i,j)th过程时间常数和采样时间有关。空间相互作用矩阵Pij和时间响应hij(z)通过标准工业软件工具由输入-输出数据识别,例如在D.M.Gorinevsky和C.Gheorghe的“Identification tool forcross-directional processes”,IEEE Transactions on ControlSystems Technology,vol.11,no.5,pp 629-640,Sept,2003.中所描述的。阵列中的每个致动器都被认为具有相同的空间响应形状。另外,空间相互作用矩阵Pij被近似为矩形轮换矩阵(RCMs)。在J.Fan等人的“Two dimensional frequency analysis for unconstrainedmodel predictive control of cross directional processes”,Automa tica,vol.40,no.11,pp.1891-1903,2004.中进一步介绍了RCMs的概念以及它们的性质。实际上,这种近似通常是较小的。
通常,为了正确调整CD MPC控制系统,有必要要求调整参数不随测量和促动的物理单位改变而改变。还要求调整参数不随测量和促动的分辨率改变而改变。测量和促动的分辨率改变意味着仅仅输入和输出的维数被增大/减小而空间响应形状不改变(也就是,Pij的列向量的维数和Pij的列向量的数量被增大/减小但是它们的形状不改变)。为了完成这种任务,有必要对(1)中的过程模型G(z)定标(scaled)。在正确的预处理之后,可以通过在例如Lu的美国专利5,574,638中介绍的标准技术来实现定标,该专利转让给了Honeywell International,Inc.,其通过参考被并入这里。下面给出用于多CD过程模型G(z)的新颖定标函数。
定标
可以将多横向过程模型G(z)转换成一系列小矩阵模型,其根据如下“时间和空间频率”部分所示的空间频率进行分类。稳态时(z=1),这些模型是:
这里vk是第k阶空间频率(k=0,1,...p),vp是(除了那些零相互作用subplant之外的)所有subplant的最小“截止”频率。可以将一个subplant vc的“截止”频率限定为
这里(4)中的全部“可控制的”频率模型g(vk,1)的平均模型用于表现非定标模型,
这里
p是来自全部subplant vp的上述最小“截止”频率。
令定标函数为对角线矩阵dy∈R3×3和du∈R4×4,于是“可控制的”频率范围中的定标平均模型gs∈R3×4由下式定义:
dy和du用于使g的条件数最小化并可以通过美国专利5,574,638中介绍的迭代方法来解决。
用于(1)中的输入U(z)和输出Y(z)的定标函数为:
Dy=diag(dy(I)Im,dy(2)Im,dy(3)Im),(8)
这里算子“diag”表示产生对角线矩阵,而Ix意味着x乘x单位矩阵。因此,可以通过下式计算定标的输入、输出和多CD过程模型
Us(z)=DuU(z),(10)
Ys(z)=DyY(z),(11)
为了简化符号,从现在开始,将模型G(z)认为是不受输入和输出的物理单位影响的定标模型。
也已经假设真实的过程响应属于一组可行的响应模型,其可以描述为
Gp(z)∈∏:={G(z)+Δ(z):σ(Δ(eiω))<α||G(z)||∞}。 (13)
这里α∈(0,1)是将扰动转移矩阵Gp(z)限制到(1)中的标称模型G(z)附近的正标量数量,Δ(z)是非结构附加的模型不确定性,σ表示最大奇异值,z=eiω和符号ω表示时间频率。
工业CD模型预测控制器
本发明方法可以自动调整空间分布过程使用的任何适当的大规模多元模型预测控制器。在该示例中,该方法在用于CD过程的近来开发的工业MPC控制器中被采用,所述CD过程在转让给了HoneywellASCa Inc的Backstrom和He的美国专利6,807,510中被介绍,该专利通过参考并入在此。对于这种工业MPC,二次方程式编程问题如下,
受到(1)和致动器阵列上的下列约束:
ΩΔU(k)≤b-φU(k-1) (15)
这里 是预测输出分布,Ysp是用于受控纸张性质的选点,Hp是预测水平,并将控制水平选择为1,这是由于相对简单的动态和计算考虑,致动器阵列ΔU(k)(=U(k)-U(k-1))和Ui中的变化是用于致动器分布的目标位置,而权重
Q1=diag(q11Im,q12Im,q13Im), (16)
Qj=diag(qj1In1,...,qj4In4),j=2,3,(17)
这里diag(x1,...,xn)意味着(块)对角线矩阵具有(块)对角线元素x1,...,xn,Ix表示x乘x单位矩阵, 被称为“弯曲力矩矩阵”;qii(i=1至3)、qjk(j=2,3,4)和k=1至4是正标量数量,其需要被调整,(15)中的Ω、φ和b分别是约束矩阵和向量,其由致动器上的物理限制的考虑、或有时由过程的热考虑分别得出。
闭环转移函数
严格地讲,(14)-(15)中的模型预测控制器是时间变化的非线性控制器。非常难以直接分析闭环系统的性质。因此,在该示例中,策略是首先忽略约束(15),然后计算闭环转移矩阵并基于这些矩阵分析系统的性能和鲁棒性,并最终为了实施而再次引入约束(15)。当再次引入约束时这种技术不保证性能,所以它将提供对于通常情况下可接受的性能,此时没有约束起作用。它以这种方式提供电流试误加闭环仿真方法上的优点。目标是设计权重Q1至Q4,使用于非约束情况的闭环转移函数具有二维频率域中的所需性质。
图3示出具有非约束MPC的闭环系统的框图。来自(1)-(3)和(4)的常数预滤器矩阵Kr∈R576×576和复转移矩阵K(z)∈C190×576的导出是标准的。
图3中的线性闭环系统对于(13)中的全部对象(plant)Gp(z)都是鲁棒稳定的,如果它是标称稳定的并且
R(z)=K(z)·[I-G(z)K(z)]-1, (21)
这里Δ(z)是(13)中的非结构不确定性,而结果(20)是由小增益定理得出的,R(z)∈C190×576(也被称为控制敏感度函数)是使致动器选点U(z)与输出扰动D(z)相连接的标称闭环转移矩阵。
可以由图3的使被测量的分布Y(z)与输出扰动D(z)相联系的标称转移矩阵将CD控制性能标准定义为,
S(z)=[I-G(z)K(z)]-1, (22)
这里S(z)∈C576×576也被称为灵敏度函数。
工业模型预测控制器(14)中可用的设计自由度是权重矩阵Q(1)至Q(4)和预测水平Hp。这些参数应被自动调整使得:
(i)由(14)和(1)给出的闭环控制系统对于(13)中的所有Gp(z)∈∏都是稳定的;
(ii)稳态时,(14)中的误差分布E(k)的2-范数被保持得尽量小;和
(iii)二维频率域中的所需闭环性质。
时间和空间频率
这部分概述关于矩形轮换矩阵(RCMs)和矩形轮换块矩阵提出的部分工作。
已知这些单CD阵列系统Gij(z)是空间频率限带的。因此,(2)中的Gij(z)近似为空间频率限带RCMs。令 表示空间频率域中的模型,这里Fm∈C192×192和 是复傅立叶矩阵而H意味着共轭转置。注意 具有如图4中所示的对角线非零元素,其表示由于其RCM特性导致的空间频率增益。在图4中,gij(v0,z)是零空间频率v0时的增益,gij(vk,z)和是共轭的并表示具有不同相位的相同空间频率增益,对于偶数nj有p=q-1=nj/2-1,而对于奇数nj有p=q-1=(nj-1)/2。
二维频率域中的多阵列CD过程模型可以通过下式得到
这里Py∈R576×576和Pu∈R190×190是置换(酉)矩阵,而Fy和Fu是Fy=diag(Fm,Fm,Fm),
对于多阵列系统,如果(2)中的标称对象G(z)是RCM-块矩阵并且其子对象(subplant)Gij(z)是空间限带的,则图(3)中的非约束MPC控制器K(z)、(22)中的灵敏度函数S(z)和(21)中的转移矩阵R(z)也具有空间限带的RCM-块结构。这使闭环转移函数的分析进入到二维频率域。
自动调整
由(14)式,调整参数是预测水平Hp和权重Q1至Q4。通常(14)中预测水平Hp的设计是基于过程的动态并被选择成最大时间延迟Tdij与子对象的最大时间常数的3倍之和。
可以通过注意不失去一般性来减少调整参数选择中自由度的数量。(14)中的Q1可以被固定并通过Q2、Q3和Q4修改闭环性能。通常在CD MPC中,权重矩阵Q2被认为是影响闭环系统的动态而Q3影响更慢的行为和稳态性能,Q4影响高空间频率行为。
从J.Fan等人的“Two-dimensional frequency analysis forunconstrained model predictive control of cross-directionalprocesses”,Automatica,vol.40,no.11,pp.1891-1903,2004,中的定理3,已知二维频率域中(22)中的灵敏度函数S(z)在稳态(z=1)时不依赖于Q2,而在空间频率v=0时不依赖于Q4。这意味着可以分别调整闭环控制系统,首先在稳态时用于空间频率域中的空间调整,接着在空间频率v=0时用于时间频率域中的时间调整。在图6a中概述调整策略。调整目标是令控制器尽可能主动(aggressive),同时仍遵守(20)中的鲁棒稳定条件,也就是,令(17)、(18)中的q2j、q3j、q4j(j=1)尽可能小。
在这种设计中,已经假设用户将基于产品的要求和致动器的物理考虑来提供相对重要的被控制性质和相对主动的被操纵致动器阵列。对于默认而言,每个被控制的性质都被认为是相同重要的并且每个致动器阵列都同样主动地工作。令qi1,j和qik,j(k=2,3,4)分别表示第i次测量的相对重要参数和第j次致动器阵列的相对主动参数。对于这些参数的默认值都为1。
可以在如图6b中所示的CD MPC系统的自动调整方法中采用如图6a中所示的调整策略,这里用于实施模型识别工具以识别过程模型的优选技术是执行标准撞击测试和模型识别步骤,其例如在美国专利Gorinevsky和Heaven的6,086,237中被介绍,该专利转让给了Honeywell International,Inc,通过参考并入在此。使用过程模型,可以在可能的预计算范围内自动搜索易于调节的调整参数。被确定的调整参数被接着加载给控制CD过程的MPC,例如超级压光机。本发明方法保证控制器是最佳的和加鲁棒的,也就是,控制系统能够在多种操作条件下始终如一地操作。
空间调整
为了简化自动调整过程,用于空间调整参数Q3和Q4的全局空间比例参数ρ被使用。也就是,
由上述讨论,除了相对重要参数之外的测量权重被固定,
Q1=Qi1=diag(qi11Im,qi12Im,qi13Im). (26)
稳态时(z=1)CD MPC的灵敏度函数S(z)和控制灵敏度函数R(z)为:
这里Gss是定标稳态模型,Gsum是与预测水平和过程动态相关的矩阵。
空间调整目标是找到最小的ρ以满足(20)中的鲁棒稳定条件。首先得出使高限制ρhigh被保证以满足鲁棒稳定条件的范围ρ∈(ρlow,ρhigh)。保守性能(conservative performance)和自动调整程序中的该结果接着使用关于ρ的分半搜索方法,以便降低保守性(conservativeness)。
对于ρ而言最主动的调整数量是ρ=0。如果ρ→0,则S(1)→0而σ(R(1))→+∞,这是由于Gss的最小奇异值几乎为0的事实。因此,ρ的下限为ρlow=0。
最鲁棒的调整为ρ→+∞。如果ρ→+∞,则σ(R)→0。不过,为了满足鲁棒稳定条件(20),需要ρ足够大。由(28),如果令
ρ>(1+α)σ((Qi3+Qi4)-1)σ(Gsum)σ(Qi1)σ(Gss),,(29)
则,
这里α表示(13)中模型不确定性的程度。这意味着(20)中的鲁棒稳定条件在稳态时被满足。注意Qi3和Qi4对角线元素qi3j和qi4j(j=1至4)是正数。
由(29),对于Qi1、Qi3和Qi4的所有可能值,ρ的高限(稳态时其满足(20)中的鲁棒稳定条件)为:
如前面部分中所述的那样,作为RCM-块矩阵,可以通过类似于g(vk,z)的r(vk,1)来计算(28)中的控制灵敏度函数R(1)的奇异值。对于每个空间频率,为了满足鲁棒条件,从0到ρhigh搜索最主动的ρ。
搜索方法是分半方法。该方法将保证ρ是收敛的并进一步满足鲁棒稳定条件。
图7中示出的简单示例用于说明减小的Q3改善了空间性能并减小了鲁棒稳定裕度。为了清楚起见该示例仅示出一个致动器阵列和一个测量阵列的情况。注意图7a中的空间灵敏度函数的最大峰值是1,这不同于其时间对应部分(见图8a)。原因在于控制器在空间域中是非因果的,而控制器在时间域中是因果的。
时间调整
在上述空间调整程序之后,得知空间调整参数Q3和Q4。在这部分中,我们将介绍怎样设计时间调整参数Q2,其被定义为
这里qi2j(j=1至4)是用于4个致动器阵列的用户定义的相对时间主动调整权重(默认值为1),而γ是全局时间比例参数。时间调整目标是找到适当的γ使面对关于最大灵敏度的约束时鲁棒稳定条件被满足。首先得出γ∈(γlow,γhigh)的范围,接着在该范围中搜索γ,用于满足鲁棒稳定条件和所需的闭环性能。
空间频率v=0时,灵敏度函数S(z)和控制灵敏度函数R(z)为:
s(v0,z)=[I3+g(v0,z)k(v0,z)]-1, (34)
r(v0,z)=[I4+k(v0,z)g(v0,z)]-1k(v0,z),(35)
这里g(v0,z)∈C3×4和k(v0,z)∈C4×3分别是空间频率v0=0时的对象模型G(z)和控制器K(z),k(v0,z)具有下列结构,
k(v0,z)=[γΘ(z)+Φ(z)]-1ψ(z),(36)
这里Θ(z)、Φ(z)和ψ(z)是权重q12j、q13j、qi1j、预测水平Hp和对象模型的时间参数的函数。
γ的最小(或最大的动态主动)值为γlow=0(也就是Q2=0)。γ的最大(或最大动态惰性(sluggish))值为+∞。但是,我们希望找到不同于+∞的γ的合理最大值。γ的该合理最大值可以使灵敏度函数的最高峰值为ηlow(例如,ηlow=1.01)并满足鲁棒稳定条件。也就是,我们希望令γ=γhigh,以使
||s(v0,eiω,γ)||max≤σ(s(v0,eiω,γ))≤ηlow,(37)
这里||s||max表示矩阵s的最大元素。
由(37)和(38)可直接获得γhigh。在获得γ的范围之后,设立目标峰值用于图8中的灵敏度函数和控制灵敏度函数(见ηsp和ksp)。然后,从γhigh到0按顺序搜索全局时间比例以达到灵敏度函数的目标峰值ηsp而不超过控制灵敏度函数的峰值ksp和违反(38)中的鲁棒稳定条件。如果通过任何γ都不能达到ηsp,则该方法努力找到γ以使控制灵敏度函数的峰值等于ksp而不违反(38)中的鲁棒稳定条件。在搜索全局时间比例γ的过程中,首先使用更小范围(γhigh,0.1γhigh)来替代整个范围(γhigh,0),接着使用分半方法来搜索γ,并看峰值是否足够接近它们的目标(也就是说目标的±1%)。如果找到解,则停止搜索。否则,移动到下一个更小范围并再次搜索直到找到解。这种搜索程序的原因是找到相对更大的γ,也就是更鲁棒的控制器是优选的,如果更鲁棒的控制器能够满足性能要求,即最大灵敏度函数峰值和鲁棒稳定条件的话。
图8示出一个简单示例,其说明空间频率v=0时,时间频率域中的灵敏度函数和控制灵敏度函数。该示例仅示出对于一个受控的纸张性质和一个致动器阵列系统的情况。
示例
来自造纸厂的超级压光机中的CD过程被用作难以通过传统CD控制器进行控制的示例。图9示出空间频率域中稳态时过程模型的奇异值。很明显过程模型的病态是由于低空间频率范围中多致动器阵列和测量阵列的方向性和高空间频率范围中的低增益两方面的缘故。在调整之后,对于空间频率域中稳态时每个测量的灵敏度函数在图10a、b和c中示出,而时间频率域中零空间频率时的灵敏度函数在图10d、e和f中说明。
在工业中,测量分布和目标分布(2-σ)之间的误差分布的2-范数的两倍被认为是控制器的空间性能指数。图11示出通过运行半实物(hardware-in-loop)仿真器中的传统控制器(TC)和MPC的时间性能。对于传统控制器,纸正面和纸反面蒸汽盒致动器阵列用于分别控制纸正面和纸反面的光泽,而两个感应加热致动器阵列用于控制纸厚度测径器分布。很明显MPC控制器的时间性能比传统控制器更好。在使用传统控制器之后,纸正面光泽的2-σ变坏。在表1中总结了两种控制器的空间性能。很明显MPC控制器的空间性能比传统控制器的空间性能更好。
表1
控制之前2-σ | (TC)控制之后2-σ | (MPC)控制之后2-σ | |
纸正面光泽 | 2.8711 | 4.0326 | 1.5450 |
纸厚度测径器 | 0.0882 | 0.0758 | 0.0408 |
纸反面光泽 | 3.5333 | 2.7613 | 2.3109 |
上面已经介绍了本发明的原理、优选实施例和操作模式。不过,本发明不应该被理解为限于所讨论的特殊实施例。而是,应该将上述实施例认为是示例性的而非限制性的,并且应该意识到在不背离由所附权利要求限定的本发明的范围的情况下,本领域技术人员可以对这些实施例进行更改。
Claims (7)
1.一种用于空间分布过程的多元模型预测控制器的自动调整方法,包括下列步骤:
(a)识别用于作为闭环控制系统的一部分的空间分布过程的过程模型;
(b)对过程模型的输入和输出进行定标;
(c)将过程模型转换到二维频率域中;
(d)将多元模型预测控制器转换到二维频率域中;
(e)确定空间调整参数关于模型不确定性的最佳解,包括由于模型不确定性搜索空间频率域用于稳定闭环系统,以及满足空间频率域中的鲁棒稳定条件;
(f)确定时间调整参数关于模型不确定性的最佳解,包括由于模型不确定性搜索时间频率域用于稳定闭环系统,以及满足时间频率域中的鲁棒稳定条件;以及
其中,步骤(e)和(f)在迭代过程中被重复直到适当的空间调整参数和时间调整参数被确定,以使作为部分闭环控制系统的空间分布过程显示出鲁棒稳定性。
2.根据权利要求1所述的方法,其中步骤(e)和(f)每一个都包括计算用于多元模型预测控制器的空间和时间调整参数的可能范围的步骤。
3.根据权利要求1所述的方法,其中步骤(e)包括计算至少一个权重矩阵的调整参数值和二次方程式函数中的预测水平,而步骤(f)包括计算至少一个权重矩阵的最佳调整参数值和二次方程式函数中的预测水平。
4.根据权利要求1所述的方法,其中步骤(e)包括在空间频率域中搜索全局空间调整参数,而步骤(f)包括在时间频率域中搜索全局时间调整参数。
5.根据权利要求1所述的方法,其中步骤(e)包括通过空间频率域中的空间灵敏度函数搜索最佳空间调整权重。
6.根据权利要求1所述的方法,其中步骤(f)包括通过时间频率域中的时间灵敏度函数和控制灵敏度函数搜索最佳空间调整权重。
7.一种在空间分布的多元阵列横向过程中形成材料的系统,其中该系统包括:
至少两组致动器阵列,每组致动器阵列在横向上都毗连所述材料分布,其中每组致动器阵列都是可控的以改变材料的性质;
用于测量和获得关于材料在横向上的性质的性质数据的装置;和
用于为多阵列横向过程提供横向控制的多元模型预测控制器MPC,其中多元模型预测控制器包括用于自动调整多元模型预测控制器的装置,其中多元模型预测控制器响应于指示性质数据的信号,给至少两组致动器阵列提供信号以改变材料的性质,以及其中所述用于自动调整多元模型预测控制器的装置包括:
(a)识别用于横向过程的过程模型的装置;
(b)对过程模型的输入和输出进行定标的装置;
(c)将过程模型转换到二维频率域中的装置;
(d)将多元模型预测控制器转换到二维频率域中的装置;
(e)计算用于多元模型预测控制器的空间和时间调整参数的可能范围的装置;
(f)确定空间调整参数关于模型不确定性的最佳解的装置,通过由于模型不确定性搜索空间频率域用于稳定闭环系统,以及满足空间频率域中的鲁棒稳定条件,而确定空间调整参数关于模型不确定性的最佳解;和
(g)确定时间调整参数关于模型不确定性的最佳解的装置,通过由于模型不确定性搜索时间频率域用于稳定闭环系统,以及满足时间频率域中的鲁棒稳定条件,而确定时间调整参数关于模型不确定性的最佳解,以及
其中,在迭代过程中装置(f)和装置(g)重复运行,直到适当的空间调整参数和时间调整参数被确定,以使作为闭环控制系统的一部分的空间分布过程显示出鲁棒稳定性。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/260,809 US7650195B2 (en) | 2005-10-27 | 2005-10-27 | Automated tuning of large-scale multivariable model predictive controllers for spatially-distributed processes |
US11/260,809 | 2005-10-27 | ||
PCT/US2006/041641 WO2007050692A2 (en) | 2005-10-27 | 2006-10-25 | Automated tuning of large-scale multivariable model predictive controllers for spatially-distributed processes |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101346675A CN101346675A (zh) | 2009-01-14 |
CN101346675B true CN101346675B (zh) | 2011-08-03 |
Family
ID=37781810
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2006800488517A Expired - Fee Related CN101346675B (zh) | 2005-10-27 | 2006-10-25 | 自动调整用于空间分布过程的大规模多元模型预测控制器 |
Country Status (7)
Country | Link |
---|---|
US (1) | US7650195B2 (zh) |
EP (1) | EP1941328B1 (zh) |
CN (1) | CN101346675B (zh) |
AT (1) | ATE471540T1 (zh) |
CA (1) | CA2627777C (zh) |
DE (1) | DE602006014984D1 (zh) |
WO (1) | WO2007050692A2 (zh) |
Families Citing this family (54)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070078873A1 (en) * | 2005-09-30 | 2007-04-05 | Avinash Gopal B | Computer assisted domain specific entity mapping method and system |
US8036760B2 (en) * | 2005-10-04 | 2011-10-11 | Fisher-Rosemount Systems, Inc. | Method and apparatus for intelligent control and monitoring in a process control system |
US7444191B2 (en) | 2005-10-04 | 2008-10-28 | Fisher-Rosemount Systems, Inc. | Process model identification in a process control system |
US7738975B2 (en) * | 2005-10-04 | 2010-06-15 | Fisher-Rosemount Systems, Inc. | Analytical server integrated in a process control network |
US20070225835A1 (en) * | 2006-03-23 | 2007-09-27 | Yucai Zhu | Computer method and apparatus for adaptive model predictive control |
US7949417B2 (en) | 2006-09-22 | 2011-05-24 | Exxonmobil Research And Engineering Company | Model predictive controller solution analysis process |
US7987005B2 (en) * | 2006-12-19 | 2011-07-26 | Chevron U.S.A. Inc. | System, method and program for dynamic control and optimization of a process having manipulated and controlled variables |
CN100465825C (zh) * | 2007-06-13 | 2009-03-04 | 袁璞 | 变结构非线性模型预估控制器 |
KR101022661B1 (ko) * | 2008-02-28 | 2011-03-22 | 한국전자통신연구원 | 금속-절연체 전이(mit) 소자를 구비한 대전류 제어회로,그 대전류 제어회로를 포함하는 시스템 |
US8155932B2 (en) | 2009-01-08 | 2012-04-10 | Jonas Berggren | Method and apparatus for creating a generalized response model for a sheet forming machine |
US8209048B2 (en) | 2009-01-12 | 2012-06-26 | Abb Automation Gmbh | Method and apparatus for creating a comprehensive response model for a sheet forming machine |
US9760073B2 (en) * | 2010-05-21 | 2017-09-12 | Honeywell International Inc. | Technique and tool for efficient testing of controllers in development |
US8862249B2 (en) * | 2010-05-27 | 2014-10-14 | Honeywell Asca Inc. | Apparatus and method for modeling and control of cross-direction fiber orientation processes |
US9020029B2 (en) * | 2011-01-20 | 2015-04-28 | Alcatel Lucent | Arbitrary precision multiple description coding |
US9436179B1 (en) | 2013-03-13 | 2016-09-06 | Johnson Controls Technology Company | Systems and methods for energy cost optimization in a building system |
US10418833B2 (en) | 2015-10-08 | 2019-09-17 | Con Edison Battery Storage, Llc | Electrical energy storage system with cascaded frequency response optimization |
US9852481B1 (en) | 2013-03-13 | 2017-12-26 | Johnson Controls Technology Company | Systems and methods for cascaded model predictive control |
US9235657B1 (en) | 2013-03-13 | 2016-01-12 | Johnson Controls Technology Company | System identification and model development |
DK177915B1 (en) * | 2013-05-28 | 2015-01-05 | Core As | Process control method |
US10386820B2 (en) | 2014-05-01 | 2019-08-20 | Johnson Controls Technology Company | Incorporating a demand charge in central plant optimization |
EP3304223A4 (en) * | 2015-06-03 | 2019-01-23 | Honeywell Limited/Honeywell Limitée | METHOD AND DEVICE FOR ROBUST VOTING MODEL-BASED PROCESS CONTROL DEVICES WITH UNMISTAKED MIMO PROCESSES |
US10190789B2 (en) | 2015-09-30 | 2019-01-29 | Johnson Controls Technology Company | Central plant with coordinated HVAC equipment staging across multiple subplants |
US10418832B2 (en) | 2015-10-08 | 2019-09-17 | Con Edison Battery Storage, Llc | Electrical energy storage system with constant state-of charge frequency response optimization |
US10190793B2 (en) | 2015-10-08 | 2019-01-29 | Johnson Controls Technology Company | Building management system with electrical energy storage optimization based on statistical estimates of IBDR event probabilities |
US10222083B2 (en) | 2015-10-08 | 2019-03-05 | Johnson Controls Technology Company | Building control systems with optimization of equipment life cycle economic value while participating in IBDR and PBDR programs |
US10742055B2 (en) | 2015-10-08 | 2020-08-11 | Con Edison Battery Storage, Llc | Renewable energy system with simultaneous ramp rate control and frequency regulation |
US10222427B2 (en) | 2015-10-08 | 2019-03-05 | Con Edison Battery Storage, Llc | Electrical energy storage system with battery power setpoint optimization based on battery degradation costs and expected frequency response revenue |
US10250039B2 (en) | 2015-10-08 | 2019-04-02 | Con Edison Battery Storage, Llc | Energy storage controller with battery life model |
US11210617B2 (en) | 2015-10-08 | 2021-12-28 | Johnson Controls Technology Company | Building management system with electrical energy storage optimization based on benefits and costs of participating in PDBR and IBDR programs |
US10389136B2 (en) | 2015-10-08 | 2019-08-20 | Con Edison Battery Storage, Llc | Photovoltaic energy system with value function optimization |
US10564610B2 (en) | 2015-10-08 | 2020-02-18 | Con Edison Battery Storage, Llc | Photovoltaic energy system with preemptive ramp rate control |
US10554170B2 (en) | 2015-10-08 | 2020-02-04 | Con Edison Battery Storage, Llc | Photovoltaic energy system with solar intensity prediction |
US10283968B2 (en) | 2015-10-08 | 2019-05-07 | Con Edison Battery Storage, Llc | Power control system with power setpoint adjustment based on POI power limits |
US10700541B2 (en) | 2015-10-08 | 2020-06-30 | Con Edison Battery Storage, Llc | Power control system with battery power setpoint optimization using one-step-ahead prediction |
US10197632B2 (en) | 2015-10-08 | 2019-02-05 | Taurus Des, Llc | Electrical energy storage system with battery power setpoint optimization using predicted values of a frequency regulation signal |
US10459428B2 (en) * | 2016-03-08 | 2019-10-29 | Honeywell Limited | Optimal closed-loop input design for identification of flat-sheet process models |
US10594153B2 (en) | 2016-07-29 | 2020-03-17 | Con Edison Battery Storage, Llc | Frequency response optimization control system |
US10778012B2 (en) | 2016-07-29 | 2020-09-15 | Con Edison Battery Storage, Llc | Battery optimization control system with data fusion systems and methods |
US11449046B2 (en) * | 2016-09-16 | 2022-09-20 | Honeywell Limited | Model-plant mismatch detection with support vector machine for cross-directional process behavior monitoring |
US10358771B2 (en) | 2016-09-23 | 2019-07-23 | Honeywell Limited | Method of designing model predictive control for cross directional flat sheet manufacturing processes to guarantee spatial robustness and to prevent actuator picketing |
US10309059B2 (en) | 2016-09-23 | 2019-06-04 | Honeywell International Inc. | Method of designing model predictive control for cross directional flat sheet manufacturing processes to guarantee temporal robust stability and performance |
US10174456B2 (en) | 2016-09-23 | 2019-01-08 | Honeywell Limited | Technique to improve paper machine cross-directional model predictive control performance by creating a measurement profile reference trajectory |
WO2018161147A1 (en) | 2017-03-06 | 2018-09-13 | Honeywell Limited | Method and apparatus for designing model-based control having spatial robustness for multiple-array cross- direction (cd) web manufacturing or processing systems or other systems |
WO2018161148A1 (en) | 2017-03-06 | 2018-09-13 | Honeywell Limited | Method and apparatus for designing model-based control having temporally robust stability and performance for multiple-array cross-direction (cd) web manufacturing or processing systems or other systems |
US10678197B2 (en) * | 2017-03-06 | 2020-06-09 | Honeywell Limited | Method and apparatus for designing model-based control having temporally robust stability and performance for multiple-array cross-direction (CD) web manufacturing or processing systems or other systems |
CN107423811B (zh) * | 2017-07-04 | 2018-12-14 | 河海大学 | 基于bp人工神经网络和情景模拟组合的径流变化归因识别方法 |
US10838441B2 (en) | 2017-11-28 | 2020-11-17 | Johnson Controls Technology Company | Multistage HVAC system with modulating device demand control |
US10838440B2 (en) | 2017-11-28 | 2020-11-17 | Johnson Controls Technology Company | Multistage HVAC system with discrete device selection prioritization |
US11163271B2 (en) | 2018-08-28 | 2021-11-02 | Johnson Controls Technology Company | Cloud based building energy optimization system with a dynamically trained load prediction model |
US11159022B2 (en) | 2018-08-28 | 2021-10-26 | Johnson Controls Tyco IP Holdings LLP | Building energy optimization system with a dynamically trained load prediction model |
CN109519957B (zh) * | 2018-09-29 | 2019-11-12 | 东南大学 | 一种超超临界锅炉闭环燃烧优化控制方法 |
CN111835309B (zh) * | 2020-07-27 | 2024-03-26 | 中国工程物理研究院总体工程研究所 | 一种持续载荷模拟器运动奇异性滤波控制方法 |
CN112034713B (zh) * | 2020-09-07 | 2021-10-19 | 山东大学 | 非理想网络环境下移动目标的最优状态估计方法及系统 |
CN114351496B (zh) * | 2021-12-17 | 2023-07-18 | 浙江华章科技有限公司 | 一种网压部真空系统压力自动整定方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5400258A (en) * | 1993-09-03 | 1995-03-21 | Measurex Corporation | Automatic cross-directional control zone alignment for sheetmaking systems |
CN1231436A (zh) * | 1999-04-14 | 1999-10-13 | 袁璞 | 通用多变量模型预估协调控制方法 |
US6064809A (en) * | 1998-06-05 | 2000-05-16 | The Board Of Trustees Of The University Of Illinois | Fast model predictive ellipsoid control process |
CN1645276A (zh) * | 2005-02-03 | 2005-07-27 | 上海交通大学 | 化工双输入输出过程的分布式pi和pid控制器的定量整定方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5121332A (en) * | 1989-03-31 | 1992-06-09 | Measurex Corporation | Control system for sheetmaking |
EP0524317A4 (en) | 1991-02-08 | 1995-02-15 | Tokyo Shibaura Electric Co | Model forecasting controller |
US5351184A (en) | 1993-01-26 | 1994-09-27 | Honeywell Inc. | Method of multivariable predictive control utilizing range control |
US5574638A (en) | 1995-04-03 | 1996-11-12 | Lu; Zhuxin J. | Method of optimal scaling of variables in a multivariable predictive controller utilizing range control |
US6086237A (en) | 1995-12-13 | 2000-07-11 | Measurex Devron Inc. | Automated identification of web shrinkage and alignment parameters in sheet making machinery using a modeled actuator response profile |
US5893055A (en) | 1997-05-30 | 1999-04-06 | Abb Industrial Systems, Inc. | Two-dimensional web property variation modeling and control |
US6253113B1 (en) | 1998-08-20 | 2001-06-26 | Honeywell International Inc | Controllers that determine optimal tuning parameters for use in process control systems and methods of operating the same |
US6826521B1 (en) | 2000-04-06 | 2004-11-30 | Abb Automation Inc. | System and methodology and adaptive, linear model predictive control based on rigorous, nonlinear process model |
WO2002005042A2 (en) | 2000-07-12 | 2002-01-17 | Aspen Technology, Inc. | Automated closed loop step testing of process units |
US6807510B1 (en) | 2003-05-05 | 2004-10-19 | Honeywell Acsa Inc. | Model predictive controller for coordinated cross direction and machine direction control |
-
2005
- 2005-10-27 US US11/260,809 patent/US7650195B2/en active Active
-
2006
- 2006-10-25 WO PCT/US2006/041641 patent/WO2007050692A2/en active Application Filing
- 2006-10-25 DE DE602006014984T patent/DE602006014984D1/de active Active
- 2006-10-25 CN CN2006800488517A patent/CN101346675B/zh not_active Expired - Fee Related
- 2006-10-25 CA CA2627777A patent/CA2627777C/en active Active
- 2006-10-25 AT AT06817371T patent/ATE471540T1/de not_active IP Right Cessation
- 2006-10-25 EP EP06817371A patent/EP1941328B1/en not_active Not-in-force
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5400258A (en) * | 1993-09-03 | 1995-03-21 | Measurex Corporation | Automatic cross-directional control zone alignment for sheetmaking systems |
US6064809A (en) * | 1998-06-05 | 2000-05-16 | The Board Of Trustees Of The University Of Illinois | Fast model predictive ellipsoid control process |
CN1231436A (zh) * | 1999-04-14 | 1999-10-13 | 袁璞 | 通用多变量模型预估协调控制方法 |
CN1645276A (zh) * | 2005-02-03 | 2005-07-27 | 上海交通大学 | 化工双输入输出过程的分布式pi和pid控制器的定量整定方法 |
Non-Patent Citations (4)
Title |
---|
Junqiang Fan,Gregory E.Stewart ,Guy A.Dumont.Two-dimensional frequency analysisforunconstrainedmodelpredictive control of cross-directionalprocesses.AUTOMATIC40 11.2004,40(11),1891-1903. |
Junqiang Fan,Gregory E.Stewart ,Guy A.Dumont.Two-dimensional frequency response analysis and insights forweight selectioin in cross-directional model predictive control.IEEE CONFERENCE ON DECISION AND CONTROL1 1.2003,1(1),3717-3723. |
Junqiang Fan,Gregory E.Stewart,Guy A.Dumont.Two-dimensional frequency analysisforunconstrainedmodelpredictive control of cross-directionalprocesses.AUTOMATIC40 11.2004,40(11),1891-1903. * |
Junqiang Fan,Gregory E.Stewart,Guy A.Dumont.Two-dimensional frequency response analysis and insights forweight selectioin in cross-directional model predictive control.IEEE CONFERENCE ON DECISION AND CONTROL1 1.2003,1(1),3717-3723. * |
Also Published As
Publication number | Publication date |
---|---|
US20070100476A1 (en) | 2007-05-03 |
DE602006014984D1 (de) | 2010-07-29 |
US7650195B2 (en) | 2010-01-19 |
EP1941328B1 (en) | 2010-06-16 |
WO2007050692A3 (en) | 2007-06-14 |
CA2627777C (en) | 2016-01-12 |
CA2627777A1 (en) | 2007-03-05 |
CN101346675A (zh) | 2009-01-14 |
ATE471540T1 (de) | 2010-07-15 |
EP1941328A2 (en) | 2008-07-09 |
WO2007050692A2 (en) | 2007-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101346675B (zh) | 自动调整用于空间分布过程的大规模多元模型预测控制器 | |
US8600525B1 (en) | Efficient quadratic programming (QP) solver for process control and optimization | |
VanAntwerp et al. | Cross-directional control of sheet and film processes | |
CN100520645C (zh) | 用于协调横向和机器方向控制的模型预测控制器 | |
US7496413B2 (en) | Apparatus and method for coordinating controllers to control a paper machine or other machine | |
CA2647716C (en) | Fast performance prediction of multivariable model predictive controller for paper machine cross-directional processes | |
Jin et al. | Multiple model LPV approach to nonlinear process identification with EM algorithm | |
US10095206B2 (en) | Apparatus and method for modeling and control of cross-direction fiber orientation processes | |
sheng Wang et al. | D-FNN based soft-sensor modeling and migration reconfiguration of polymerizing process | |
US10678197B2 (en) | Method and apparatus for designing model-based control having temporally robust stability and performance for multiple-array cross-direction (CD) web manufacturing or processing systems or other systems | |
CN110637261B (zh) | 用于设计基于模型的控制的方法和装置 | |
CN110622079B (zh) | 用于设计具有用于多阵列交叉方向(cd)幅材制造或处理系统或其他系统的空间稳健性的基于模型的控制的方法和装置 | |
Rigopoulos et al. | Reduced order cross-directional controller design for sheet forming processes | |
Farahmand et al. | Tensor-based blind alignment of MIMO CD processes | |
Kumar et al. | Robust Cross Directional Controller Design for Paper Machine Spatial Distributed System | |
CN108884635B (zh) | 用于识别平坦片材工艺模型的最优闭环输入设计 | |
Jiang et al. | Off-policy two-dimensional reinforcement learning for optimal tracking control of batch processes with network-induced dropout and disturbances | |
Tolonen | Simulating the performance of paper machine control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20110803 Termination date: 20201025 |