发明名称 幂函数与对数嵌套的PID变目标算法
摘要 一种幂函数与对数嵌套的PID变目标算法,主要针对的是工控设备中利用夹套控温的问题而提出,在实际生产中通过控制被控设备夹套的蒸汽压力与设备内温度的连锁,从而将被控设备的温度控制在指定的区间范围内。借助离散化的PID算法将可变的夹套蒸汽压力SpPt作为目标值进行设定,从而实现控温的功能。本算法具有正反双向调节的功能,可在升温阶段与降温阶段分别计算夹套压力的设定值,保证了大系统、滞后性强的系统可以快速响应降温过程。升温与降温曲线的参数可以做到根据不同控制设备的参数进行人工设定,结合设备特性,达到最好的调节效果。算法实现利用PLC的有限资源,实现了高精度的控制算法,做到系统开销与系统性能的平衡。算法稳定性好,控制系统运行稳定,控制曲线良好。
申请公布号 CN103513572B 申请公布日期 2016.01.20
申请号 CN201310508850.8 申请日期 2013.10.25
申请人 中国船舶重工集团公司第七〇三研究所 发明人 冯贵墨;栾洋;王震;刘庆阁;刘亮亮;王明;金磊;关纯刚;高斌;唐晓斌;佟欣
分类号 G05B11/42(2006.01)I;G05B19/05(2006.01)I 主分类号 G05B11/42(2006.01)I
代理机构 大庆禹奥专利事务所 23208 代理人 朱士文;杨晓梅
主权项 一种幂函数与对数嵌套的PID变目标算法,其特征在于:系统的关键控制参数如下:<img file="dest_path_image001.GIF" wi="50" he="22" />:控制算法中被控设备夹套蒸汽压力上限,由设备铭牌可知其夹套压力的上限值,一般在其安全允许范围内较大即可;以下仿真时<img file="dest_path_image002.GIF" wi="117" he="22" />;<img file="dest_path_image003.GIF" wi="45" he="22" />:控制算法中被控设备夹套蒸汽压力下限,保证设备到达温控点后,夹套蒸汽压力的下限值,热损失与夹套提供的热量达到平衡;以下仿真时<img file="dest_path_image004.GIF" wi="129" he="22" />;<img file="dest_path_image005.GIF" wi="54" he="22" />:被控设备在升温过程中提供给离散化PID算法的设定值,因为在实际工作中夹套蒸汽压力的调节响应速度快、随动性好,可实现实时调节;<img file="dest_path_image006.GIF" wi="38" he="20" />:被控设备在工作时实时的夹套蒸汽压力值;<img file="dest_path_image007.GIF" wi="38" he="22" />:被控设备工作时的目标保温温度,一般低于沸点;以下仿真时<img file="dest_path_image008.GIF" wi="88" he="25" />;<img file="dest_path_image009.GIF" wi="48" he="22" />:被控设备进入幂函数与对数嵌套的PID变目标算法的启动温度,在此温度之前,被控设备夹套蒸汽压力设定为<img file="836056dest_path_image001.GIF" wi="50" he="22" />,保证系统的快速升温;以下仿真时<img file="566246dest_path_image009.GIF" wi="48" he="22" />=45℃;<img file="dest_path_image010.GIF" wi="40" he="20" />:被控设备当前工作的温度,由当前温度自动计算出当前夹套蒸汽压力<img file="408300dest_path_image006.GIF" wi="38" he="20" />的设定值;公式一:<img file="dest_path_image011.GIF" wi="458" he="49" />借用了<img file="dest_path_image012.GIF" wi="52" he="21" />函数的性质,在<img file="616558dest_path_image007.GIF" wi="38" he="22" />点,即为整个控温函数的终点,希望压力保持恒定,希望控温函数的导数为0;所以幂函数的求导可知,在整个函数的最低点,其函数的导数为零;所以建立坐标系,横坐标为温度,纵坐标为压力设定值,可知现在应该求得一个幂函数,其函数的最低点为(<img file="493247dest_path_image007.GIF" wi="38" he="22" />,<img file="dest_path_image013.GIF" wi="46" he="22" />);根据函数图像平移的原理,可以求得由<img file="128759dest_path_image012.GIF" wi="52" he="21" />函数的图像平移后的方程为:<img file="dest_path_image014.GIF" wi="281" he="34" />,此函数必经过(<img file="458109dest_path_image009.GIF" wi="48" he="22" />,<img file="657010dest_path_image001.GIF" wi="50" he="22" />),带入公式可得:<img file="dest_path_image015.GIF" wi="289" he="34" />,求得<img file="dest_path_image016.GIF" wi="154" he="46" />所以温控公式变形为:<img file="dest_path_image017.GIF" wi="344" he="64" />这里值得注意的是,在项目实践中,大部分使用PLC来实现功能,以西门子系统的PLC为例,系统中只支持以e为底的幂函数;使用幂函数的基本性质:<img file="dest_path_image018.GIF" wi="80" he="30" />,所以原公式变形为:<img file="dest_path_image019.GIF" wi="307" he="33" />,其中<img file="dest_path_image020.GIF" wi="38" he="20" />,且<img file="dest_path_image021.GIF" wi="36" he="20" />时为直线方程;公式二:<img file="dest_path_image022.GIF" wi="550" he="63" />目标夹套压力设定值应具有以下的通用形式:<img file="dest_path_image023.GIF" wi="367" he="31" />由实际工况可知,k应该具有以下的性质:1:k应是以<img file="886740dest_path_image010.GIF" wi="40" he="20" />为变量的函数,形式为<img file="dest_path_image024.GIF" wi="88" he="22" />;2:<img file="dest_path_image025.GIF" wi="62" he="22" />在区间<img file="dest_path_image026.GIF" wi="142" he="24" />内为单调递减;3:<img file="dest_path_image027.GIF" wi="102" he="24" />,且<img file="dest_path_image028.GIF" wi="92" he="22" />,<img file="dest_path_image029.GIF" wi="85" he="22" />;4:<img file="817787dest_path_image025.GIF" wi="62" he="22" />在<img file="dest_path_image030.GIF" wi="88" he="22" />点的斜率尽量小;于是构造满足以上条件的函数<img file="405674dest_path_image025.GIF" wi="62" he="22" />,选用对数函数<img file="142686dest_path_image025.GIF" wi="62" he="22" />=<img file="dest_path_image031.GIF" wi="57" he="25" />为基础进行构造;由对数函数的性质可知,<img file="728388dest_path_image031.GIF" wi="57" he="25" />当<img file="dest_path_image032.GIF" wi="36" he="20" />时,<img file="971281dest_path_image026.GIF" wi="142" he="24" />区间内是单调递增的,令<img file="212907dest_path_image031.GIF" wi="57" he="25" />通过<img file="dest_path_image033.GIF" wi="70" he="22" />和<img file="dest_path_image034.GIF" wi="58" he="22" />两点;可知此时的函数是由<img file="550347dest_path_image031.GIF" wi="57" he="25" />的图像平移<img file="dest_path_image035.GIF" wi="78" he="22" />个单位得到的,所以<img file="6868dest_path_image025.GIF" wi="62" he="22" />=<img file="dest_path_image036.GIF" wi="101" he="25" />=<img file="dest_path_image037.GIF" wi="94" he="25" />;将<img file="217400dest_path_image034.GIF" wi="58" he="22" />带入<img file="946322dest_path_image025.GIF" wi="62" he="22" />可得<img file="25136dest_path_image025.GIF" wi="62" he="22" />=<img file="585431dest_path_image037.GIF" wi="94" he="25" />=<img file="dest_path_image038.GIF" wi="94" he="25" />=1;可得<img file="dest_path_image039.GIF" wi="134" he="22" />,所以<img file="167197dest_path_image025.GIF" wi="62" he="22" />=<img file="dest_path_image040.GIF" wi="84" he="28" />;但是此时注意<img file="180152dest_path_image025.GIF" wi="62" he="22" />并不满足所需函数的设定条件,继续对<img file="62658dest_path_image025.GIF" wi="62" he="22" />进行变换,利用函数性质将<img file="228191dest_path_image025.GIF" wi="62" he="22" />的图像数值反转,并向上平移1个单位,<img file="170739dest_path_image025.GIF" wi="62" he="22" />变化为:<img file="dest_path_image041.GIF" wi="176" he="28" /><img file="dest_path_image042.GIF" wi="316" he="45" />且有:<img file="484040dest_path_image026.GIF" wi="142" he="24" />,所以:<img file="dest_path_image043.GIF" wi="89" he="22" />;<img file="966974dest_path_image025.GIF" wi="62" he="22" />在区间<img file="252593dest_path_image026.GIF" wi="142" he="24" />内为单调递减,并且当<img file="366042dest_path_image030.GIF" wi="88" he="22" />时,<img file="dest_path_image044.GIF" wi="66" he="22" />取得最小值;当<img file="353590dest_path_image007.GIF" wi="38" he="22" />足够大,<img file="577898dest_path_image009.GIF" wi="48" he="22" />足够小时,<img file="dest_path_image045.GIF" wi="65" he="22" />可以尽量趋近于0;且有:<img file="449515dest_path_image028.GIF" wi="92" he="22" />,<img file="733865dest_path_image029.GIF" wi="85" he="22" />;进行图像仿真;此时构造了一条满足设定的函数,且并不是唯一满足;<img file="dest_path_image046.GIF" wi="128" he="28" />,代入夹套压力通用公式可得:<img file="dest_path_image047.GIF" wi="518" he="38" />但是这里要注意的是,西门子系列的PLC不支持<img file="dest_path_image048.GIF" wi="42" he="25" />的计算,根据<img file="21758dest_path_image048.GIF" wi="42" he="25" />的运算性质可得:<img file="dest_path_image049.GIF" wi="86" he="42" />,带入公式可得:<img file="112074dest_path_image022.GIF" wi="550" he="63" />对<img file="293657dest_path_image046.GIF" wi="128" he="28" />分析可知,对于任意给定的系统,<img file="296379dest_path_image005.GIF" wi="54" he="22" />的图像应是唯一的,满足系统比较好的变化曲线,在调温终点<img file="196202dest_path_image007.GIF" wi="38" he="22" />系统的变化率很小;通过仿真图像图即可看出,该函数曲线也存在着不足,系统的灵活性小,变化曲线唯一;并且大系统的滞后性强,在升温阶段可以快速响应;但是一旦系统工作温度低于设定值,需要继续加热时,系统的响应缓慢,不易控;公式三:<img file="dest_path_image050.GIF" wi="533" he="50" />由以上的算法推导和仿真可知,公式一和公式二部分满足了实际的使用要求,所以将公式一与公式二相结合,证明如下:同上,目标夹套压力设定值应具有以下的通用形式:<img file="824629dest_path_image023.GIF" wi="367" he="31" /><img file="dest_path_image051.GIF" wi="14" he="20" />也应符合以上的构造原则,设<img file="dest_path_image052.GIF" wi="86" he="22" />=<img file="dest_path_image053.GIF" wi="104" he="28" />,<img file="dest_path_image054.GIF" wi="61" he="22" />=<img file="dest_path_image055.GIF" wi="21" he="22" />,<img file="811783dest_path_image032.GIF" wi="36" he="20" />,<img file="dest_path_image056.GIF" wi="14" he="16" />称为弯曲系数;由上面的推导可知:1:<img file="dest_path_image057.GIF" wi="28" he="20" /><img file="dest_path_image058.GIF" wi="82" he="22" />,<img file="234675dest_path_image051.GIF" wi="14" he="20" />是一个关于<img file="434843dest_path_image010.GIF" wi="40" he="20" />的函数;2:<img file="804644dest_path_image052.GIF" wi="86" he="22" />,当<img file="570606dest_path_image026.GIF" wi="142" he="24" />时,<img file="dest_path_image059.GIF" wi="14" he="16" />是关于<img file="429978dest_path_image010.GIF" wi="40" he="20" />单调递减的,可得<img file="dest_path_image060.GIF" wi="54" he="24" />;<img file="851863dest_path_image057.GIF" wi="28" he="20" /><img file="dest_path_image061.GIF" wi="37" he="22" />=<img file="25355dest_path_image055.GIF" wi="21" he="22" />,可得<img file="98354dest_path_image060.GIF" wi="54" he="24" />时,<img file="66310dest_path_image051.GIF" wi="14" he="20" />是关于<img file="972561dest_path_image059.GIF" wi="14" he="16" />单调递减的;由函数单调性的传递关系可知,<img file="684165dest_path_image057.GIF" wi="28" he="20" /><img file="549353dest_path_image058.GIF" wi="82" he="22" />,当<img file="750527dest_path_image026.GIF" wi="142" he="24" />时,<img file="333955dest_path_image051.GIF" wi="14" he="20" />是关于<img file="396720dest_path_image010.GIF" wi="40" he="20" />单调递减的;3:<img file="381994dest_path_image057.GIF" wi="28" he="20" /><img file="691753dest_path_image058.GIF" wi="82" he="22" />,<img file="dest_path_image062.GIF" wi="44" he="24" /><img file="dest_path_image063.GIF" wi="105" he="22" />,<img file="dest_path_image064.GIF" wi="45" he="25" /><img file="dest_path_image065.GIF" wi="112" he="22" />,<img file="dest_path_image066.GIF" wi="54" he="24" />;4:<img file="dest_path_image067.GIF" wi="28" he="20" /><img file="dest_path_image068.GIF" wi="52" he="22" />=<img file="dest_path_image069.GIF" wi="101" he="22" />,将<img file="513209dest_path_image061.GIF" wi="37" he="22" />与<img file="628933dest_path_image025.GIF" wi="62" he="22" />带入得,<img file="468713dest_path_image067.GIF" wi="28" he="20" /><img file="dest_path_image070.GIF" wi="285" he="45" />在曲线调节的终点,<img file="dest_path_image071.GIF" wi="74" he="22" /><img file="dest_path_image072.GIF" wi="342" he="45" />=0通过以上的推导和构造,可知<img file="821809dest_path_image057.GIF" wi="28" he="20" /><img file="114250dest_path_image058.GIF" wi="82" he="22" />是一条符合构造条件的曲线,并在调节终点导数为零,保证了系统的平滑过渡;<img file="dest_path_image073.GIF" wi="522" he="38" />同样为了方便公式的工程化实现,对公式进行变形可得:<img file="dest_path_image074.GIF" wi="554" he="52" />经仿真,得到曲线图;通过分析升温阶段初期曲线变化快,后段变化缓慢,温控终点变化率为零,符合系统的滞后性特点;同时当系统处于降温阶段时,可以将弯曲系数<img file="784397dest_path_image056.GIF" wi="14" he="16" />取值为<img file="dest_path_image075.GIF" wi="60" he="20" />,此时调节目标迅速变大,增加系统的供气量,迅速将系统加热,转换进入升温阶段;将此算法命名为幂函数与对数嵌套的PID变目标算法。
地址 150000 黑龙江省哈尔滨市道里区洪湖路35号