发明名称 一种反应釜温度混合控制方法
摘要 本发明涉及一种反应釜温度混合控制方法。现有方法控制效果不好,不适应节能减排及环保的需求。本发明方法首先采用典型的响应曲线法设计反应釜温度过程的比例积分微分控制器,然后利用反应釜温度实时过程数据建立非最小化模型预测函数控制器所需的过程模型,再设计非最小化模型预测函数比例积分微分控制器,根据目标函数得到当前的控制参数值。本发明提出的控制方法可以有效减少反应釜温度工艺参数与实际反应釜温度工艺参数之间的误差,进一步弥补了传统控制器的不足,同时保证控制装置操作在最佳状态,使生产过程的反应釜温度工艺参数达到严格控制。
申请公布号 CN102156496A 申请公布日期 2011.08.17
申请号 CN201110086036.2 申请日期 2011.04.07
申请人 杭州电子科技大学 发明人 张日东;薛安克;王建中;葛铭;孔亚广
分类号 G05D23/30(2006.01)I 主分类号 G05D23/30(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 1.一种反应釜温度混合控制方法,其特征在于该方法的具体步骤是:步骤(1).采用响应曲线法设计反应釜温度过程的比例积分微分控制器,具体方法是:Ⅰ.将反应釜温度过程的比例积分微分控制器停留在手动操作状态,操作拨盘使其输出有阶跃变化,由记录仪表记录反应釜温度的输出值<img file="2011100860362100001DEST_PATH_IMAGE002.GIF" wi="45" he="25" />,将其转换成无量纲形式<img file="2011100860362100001DEST_PATH_IMAGE004.GIF" wi="44" he="26" />,具体是:<img file="2011100860362100001DEST_PATH_IMAGE006.GIF" wi="150" he="26" />其中,<img file="2011100860362100001DEST_PATH_IMAGE008.GIF" wi="46" he="25" />是反应釜温度输出值<img file="9421DEST_PATH_IMAGE002.GIF" wi="45" he="25" />的稳态值;Ⅱ.选取满足<img file="2011100860362100001DEST_PATH_IMAGE010.GIF" wi="189" he="26" />的两个计算点<img file="2011100860362100001DEST_PATH_IMAGE012.GIF" wi="17" he="25" />和<img file="2011100860362100001DEST_PATH_IMAGE014.GIF" wi="18" he="25" />,依据下式计算比例积分微分控制器所需要的参数<img file="2011100860362100001DEST_PATH_IMAGE016.GIF" wi="18" he="18" />、<img file="2011100860362100001DEST_PATH_IMAGE018.GIF" wi="41" he="21" />:<img file="2011100860362100001DEST_PATH_IMAGE020.GIF" wi="92" he="73" />其中,<img file="2011100860362100001DEST_PATH_IMAGE022.GIF" wi="14" he="18" />为反应釜温度过程的比例积分微分控制器输出的阶跃变化幅度,<img file="2011100860362100001DEST_PATH_IMAGE024.GIF" wi="60" he="24" />分别是计算比例积分微分控制器所需要的增益、时间常数和滞后参数;Ⅲ.计算反应釜温度过程的比例积分微分控制器的参数,具体是:<img file="2011100860362100001DEST_PATH_IMAGE026.GIF" wi="98" he="73" />其中<img file="2011100860362100001DEST_PATH_IMAGE028.GIF" wi="22" he="25" />为比例积分微分控制器的比例参数,<img file="2011100860362100001DEST_PATH_IMAGE030.GIF" wi="16" he="25" />为比例积分微分控制器的积分参数,<img file="2011100860362100001DEST_PATH_IMAGE032.GIF" wi="18" he="25" />分别为比例积分微分控制器的微分参数;步骤(2).利用反应釜温度实时过程数据建立非最小化模型预测函数控制器所需的过程模型,具体方法是:首先建立反应釜温度实时运行数据库,将过程的比例积分微分控制器停留在自动操作状态,操作拨盘使其输入有阶跃变化,通过数据采集装置采集<img file="2011100860362100001DEST_PATH_IMAGE034.GIF" wi="20" he="20" />组实时过程运行数据,将采集的实时过程运行数据作为数据驱动的样本集合,表示为<img file="2011100860362100001DEST_PATH_IMAGE036.GIF" wi="84" he="36" />,<img file="2011100860362100001DEST_PATH_IMAGE038.GIF" wi="85" he="22" />,<img file="2011100860362100001DEST_PATH_IMAGE040.GIF" wi="20" he="25" />表示第<img file="2011100860362100001DEST_PATH_IMAGE042.GIF" wi="10" he="18" />组工艺参数的输入数据,<img file="2011100860362100001DEST_PATH_IMAGE044.GIF" wi="40" he="33" />表示第<img file="942391DEST_PATH_IMAGE042.GIF" wi="10" he="18" />组工艺参数的输出值;然后以该反应釜温度实时过程运行数据集合为基础建立基于最小二乘法的离散差分方程形式的局部受控自回归滑动平均模型:<img file="2011100860362100001DEST_PATH_IMAGE046.GIF" wi="318" he="26" /><img file="2011100860362100001DEST_PATH_IMAGE048.GIF" wi="356" he="25" />其中,<img file="2011100860362100001DEST_PATH_IMAGE050.GIF" wi="34" he="22" />表示当前时刻过程模型的工艺参数的输出值,<img file="2011100860362100001DEST_PATH_IMAGE052.GIF" wi="17" he="17" />表示过程模型的工艺参数的过去时刻的输入和输出数据的集合,<img file="2011100860362100001DEST_PATH_IMAGE054.GIF" wi="34" he="22" />表示当前过程模型工艺参数对应的控制变量,<img file="2011100860362100001DEST_PATH_IMAGE056.GIF" wi="14" he="20" />为当前的递推步数,<img file="2011100860362100001DEST_PATH_IMAGE058.GIF" wi="18" he="17" />表示通过辨识得到的模型参数的集合,<img file="2011100860362100001DEST_PATH_IMAGE060.GIF" wi="17" he="18" />表示矩阵的转置,<img file="2011100860362100001DEST_PATH_IMAGE062.GIF" wi="62" he="20" />分别为对应实际过程的输出变量阶次、输入变量阶次、实际过程的时滞;采用的辨识手段为:<img file="2011100860362100001DEST_PATH_IMAGE064.GIF" wi="209" he="26" /><img file="2011100860362100001DEST_PATH_IMAGE066.GIF" wi="253" he="26" /><img file="2011100860362100001DEST_PATH_IMAGE068.GIF" wi="204" he="26" />其中,<img file="2011100860362100001DEST_PATH_IMAGE070.GIF" wi="20" he="21" />和<img file="2011100860362100001DEST_PATH_IMAGE072.GIF" wi="16" he="17" />为辨识中的两个矩阵,<img file="2011100860362100001DEST_PATH_IMAGE074.GIF" wi="132" he="25" />为遗忘因子,<img file="2011100860362100001DEST_PATH_IMAGE076.GIF" wi="13" he="20" />为单位矩阵;步骤(3).设计非最小化模型预测函数比例积分微分控制器,具体方法是:a.将步骤(2)建立的过程模型转化为差分模型形式:<img file="2011100860362100001DEST_PATH_IMAGE078.GIF" wi="554" he="22" />其中,<img file="2011100860362100001DEST_PATH_IMAGE080.GIF" wi="18" he="20" />是差分算子,<img file="2011100860362100001DEST_PATH_IMAGE082.GIF" wi="184" he="25" />为通过转换模型得到的相关系数,<img file="2011100860362100001DEST_PATH_IMAGE084.GIF" wi="14" he="16" />为对应实际过程的输入和输出变量的统一阶次;b.选取<img file="2011100860362100001DEST_PATH_IMAGE086.GIF" wi="529" he="26" />进一步将步骤a的模型转化为状态空间输入输出模型:<img file="2011100860362100001DEST_PATH_IMAGE088.GIF" wi="286" he="49" />其中具体的转化参数矩阵为<img file="2011100860362100001DEST_PATH_IMAGE090.GIF" wi="376" he="220" /><img file="2011100860362100001DEST_PATH_IMAGE092.GIF" wi="186" he="28" /><img file="2011100860362100001DEST_PATH_IMAGE094.GIF" wi="182" he="26" />c.依据步骤b建立用于非最小化模型预测函数比例积分微分控制器设计的非最小化参数模型:<img file="2011100860362100001DEST_PATH_IMAGE096.GIF" wi="232" he="22" /><img file="723573DEST_PATH_IMAGE050.GIF" wi="34" he="22" />=<img file="2011100860362100001DEST_PATH_IMAGE098.GIF" wi="42" he="22" />其中,<img file="2011100860362100001DEST_PATH_IMAGE100.GIF" wi="109" he="49" />,<img file="2011100860362100001DEST_PATH_IMAGE102.GIF" wi="286" he="52" /><img file="2011100860362100001DEST_PATH_IMAGE104.GIF" wi="17" he="25" />,<img file="2011100860362100001DEST_PATH_IMAGE106.GIF" wi="20" he="25" />,<img file="2011100860362100001DEST_PATH_IMAGE108.GIF" wi="18" he="26" />为相应的参数矩阵;d.依据步骤c计算出的模型参数整定非最小化模型预测函数比例积分微分控制器的参数,具体方法是:①建立多步最优预测输出<img file="2011100860362100001DEST_PATH_IMAGE110.GIF" wi="20" he="26" /><img file="2011100860362100001DEST_PATH_IMAGE112.GIF" wi="232" he="28" />其中,<img file="2011100860362100001DEST_PATH_IMAGE114.GIF" wi="329" he="86" />参数<img file="2011100860362100001DEST_PATH_IMAGE116.GIF" wi="482" he="90" />,<img file="2011100860362100001DEST_PATH_IMAGE118.GIF" wi="44" he="28" />为实际反应釜温度的测量值;②建立非最小化模型预测函数比例积分微分控制器的参考轨迹<img file="2011100860362100001DEST_PATH_IMAGE120.GIF" wi="16" he="20" />和目标函数<img file="2011100860362100001DEST_PATH_IMAGE122.GIF" wi="17" he="20" />;<img file="2011100860362100001DEST_PATH_IMAGE124.GIF" wi="228" he="46" />其中,<img file="DEST_PATH_IMAGE126.GIF" wi="128" he="22" />为各个时刻的输出参考轨迹;③依据②的目标函数得到当前的控制参数值<img file="DEST_PATH_IMAGE128.GIF" wi="22" he="25" /><img file="DEST_PATH_IMAGE130.GIF" wi="346" he="52" />。
地址 310018 浙江省杭州市下沙高教园区2号大街