CN110262231A - 基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 - Google Patents
基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 Download PDFInfo
- Publication number
- CN110262231A CN110262231A CN201910456527.8A CN201910456527A CN110262231A CN 110262231 A CN110262231 A CN 110262231A CN 201910456527 A CN201910456527 A CN 201910456527A CN 110262231 A CN110262231 A CN 110262231A
- Authority
- CN
- China
- Prior art keywords
- disturbance
- cstr
- input
- cstr system
- observation method
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000010354 integration Effects 0.000 claims abstract description 5
- 238000009499 grossing Methods 0.000 claims description 12
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000004088 simulation Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 3
- 238000001311 chemical methods and process Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012824 chemical production Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 239000000975 dye Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003094 perturbing effect Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000003756 stirring Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- 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)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法，并在各种条件下验证了该方法在CSTR系统中的观测效果。主要步骤为：1、建立CSTR系统动力学模型，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达；2、利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数，从而得到系统输入的估计值；3、将系统输入估计值与实际值作差来获得扰动的估计值；4、将所得到的扰动估计值前馈给CSTR系统的控制输入端，补偿掉实际扰动。本发明的主要优点是仅需要系统输出和控制输入信号，不需要对系统状态进行测量。
Description
技术领域
本发明是一种基于微分平滑的CSTR系统的扰动观测方法，属于自动控制技术领域。
背景技术
连续搅拌釜式反应器(Continuous Stirred Tank Reactor，简称为CSTR)在实际工业过程中是一种常见的化学反应器，由于其热交换能力强、产品质量稳定和低成本等优点，在化工生产的核心设备中占有重要地位，被广泛应用于染料、医药试剂、食品及合成材料等工业生产中。
注意到，CSTR系统工作范围大，设定值变化大，呈现着很强的非线性特性。对于这类非线性化工过程，除非线性调节器的参数整定得很保守，否则控制品质将会显著恶化。鉴于此，对于大部分的非线性化工过程，必须采用非线性控制方法，例如滑模控制、智能算法和鲁棒自适应控制等等。但是这些控制方法往往会受到外界扰动的影响，导致系统的动态和稳态性能差。
发明内容
本发明的目的是提出一种基于微分平滑的CSTR系统的扰动观测方法。该方法通过构造一类只需要输出信息和控制输入信息的新型扰动观测器来估算CSTR系统中的扰动，再将此得到的扰动估计量前馈给控制输入端并将其补偿掉，从而来削弱扰动对CSTR系统的稳定性影响。
本发明的技术方案为：一种基于微分平滑的CSTR系统的扰动观测方法，包括以下步骤：
建立CSTR系统动力学模型，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达；利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数，从而得到系统输入的估计值；将系统输入估计值与实际值作差来获得扰动的估计值；将所得到的扰动估计值前馈给CSTR系统的控制输入端，补偿掉实际扰动。
进一步，步骤1中，建立CSTR系统动力学模型表达式为：
其中，x1是化学反应浓度，x2是反应器温度，u是初始时的反应输入浓度，γ、Da和β均为系统系数，d(t)是外界扰动，e是自然常数。
进一步，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达的方法为：
步骤1.1.选择平滑输出变量z＝x2；
步骤1.2.由步骤1.1和CSTR系统的数学模型表达式得到：
步骤1.3.由步骤1.2的式子求导得到：
步骤1.4.再由CSTR系统的数学模型表达式得到：
步骤1.5.将步骤1.2的式子和步骤1.3的式子代入步骤1.4的式子得到：
其中，由平滑输出以及其导数构成，表示系统输入。
进一步，利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数为从而系统输入的估计值为：
进一步，将系统输入估计值与实际值作差来获得扰动的估计值为
其中，是由积分方法与比例增益k1和k2所获得。
进一步，将扰动观测值前馈至CSTR系统控制输入u端，形成闭环回路，从而将CSTR系统的实际扰动补偿掉。
本发明具有的有益效果在于：
(1)该新型扰动观测器的构造基于单输入单输出系统，仅需要系统输出和控制输入的测量，不需要系统状态的测量。
(2)本发明方法相比于现有方法大大简化和优化，易于实现，具有较强的抗扰动能力，提高了系统性能，实现了较好的控制效果。
附图说明
图1是CSTR系统设计方案的结构图；
图2是CSTR系统示意图；
图3是CSTR系统扰动观测的结构图；
图4是CSTR系统扰动观测图；
图5是CSTR系统扰动的观测误差；
图6是CSTR系统无扰动补偿时的状态图；
图7是CSTR系统有扰动补偿时的状态图。
具体实施方式
下面结合附图、CSTR系统的具体参数和仿真波形对本发明作进一步说明。
为了更清楚地说明本发明的技术方案，将对照CSTR的一组具体参数和其在MATLAB软件中的仿真结果具体说明本发明的具体实施方式。显而易见地，下面描述的仅仅是针对的一个具体的CSTR系统，对于本领域技术人员来讲，若使用该方法进行扰动观测，还应针对具体应用的CSTR系统的参数进行分析运算。进一步地，由于实际扰动是与控制输入u一起被作用于CSTR系统中的，所以将此扰动观测值直接反馈给CSTR系统控制输入u端，形成闭环控制，将实际扰动补偿掉，从而削弱扰动对CSTR系统稳定性的影响，其原理框图如图1所示。
仿真实验采用的CSTR系统示意图如图2所示。此CSTR系统主要由反应室1和搅拌机2构成。其中，反应室1主要由筒体、筒体内的内构件和换热元件组成；搅拌机2由搅拌轴及其密封装置、传动装置和搅拌器连接组成。设置该CSTR系统的参数分别是γ＝20,Da＝1,β＝0.3,u＝0,e＝2.71828。在此仿真实验中，CSTR系统的初始状态分别被选为x1(0)＝0.625,x2(0)＝0.48，并假设扰动d(t)被给为如下形式：
基于上述系统，下面通过具体实施解释本发明对CSTR系统的扰动观测方法：
步骤一、根据已给CSTR系统的具体参数，在此基础上确立采用的CSTR系统的数学模型为：
步骤二、确立应用在该具体CSTR系统上的系统输入估计量；
进一步地，(5)式整理为：
进一步地，满足：在此式中，k1和k2均为比例增益，是对积分而得到，是对积分而得到，从而估计出该系统输入；
步骤三、确定扰动的观测值；
进一步地，(6)式整理为：
其中，图3为该CSTR系统扰动观测的结构图。
步骤四、通过仿真对比，验证该新型扰动观测器的观测效果以及其对CSTR系统稳定性的影响。
进一步地，新型扰动观测器的增益K1值越小，K2值越大，其扰动观测的精度越高，估计的扰动越接近真实的扰动，但其带来的抖振也越大。本发明选取K1的取值范围为2-15，K2的取值范围为35-85，优选K1值为5，K2值为48，在提高观测扰动精度的同时，能克服抖振。假设给定以下形式的扰动：
观测到的波形如图4所示。为进一步验证观测到的扰动是否与实际扰动相近，给出扰动的观测误差图，如图5所示。
进一步地，观测出扰动，且引入扰动补偿后，通过图6和图7的比较可以发现，有扰动补偿的CSTR系统的稳定性得到了显著提高。
综上，本发明的一种基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法，并在各种条件下验证了该方法在CSTR系统中的观测效果。主要步骤为：1、建立CSTR系统动力学模型，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达；2、利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数，从而得到系统输入的估计值；3、将系统输入估计值与实际值作差来获得扰动的估计值；4、将所得到的扰动估计值前馈给CSTR系统的控制输入端，补偿掉实际扰动。本发明的主要优点是仅需要系统输出和控制输入信号，不需要对系统状态进行测量。
再次强调，本发明拟采用新型扰动观测技术对CSTR系统进行扰动观测，并对CSTR系统进行扰动补偿，从而达到削弱扰动的目的。需要注意的是，传统状态观测器一般需要状态量为条件，才能得到观测量。而本发明提出的扰动观测器只需要输出量和控制输入量即可，设计简单，容易实现。
所述实施例为本发明优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能做出的任何显而易见的改进、替换或变形均属于本发明的保护范围。
Claims (7)
1.一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，包括以下步骤：
建立CSTR系统动力学模型，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达；利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数，从而得到系统输入的估计值；将系统输入估计值与实际值作差来获得扰动的估计值；将所得到的扰动估计值前馈给CSTR系统的控制输入端，补偿掉实际扰动。
2.根据权利要求1所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，步骤1中，建立CSTR系统动力学模型表达式为：
其中，x1是化学反应浓度，x2是反应器温度，u是初始时的反应输入浓度，γ、Da和β均为系统系数，d(t)是外界扰动，e是自然常数。
3.根据权利要求2所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，选择一个平滑输出量，使得系统输入能由平滑输出量和它的各阶导数所表达的方法为：
步骤1.1.选择平滑输出变量z＝x2；
步骤1.2.由步骤1.1和CSTR系统的数学模型表达式得到：
步骤1.3.由步骤1.2的式子求导得到：
步骤1.4.再由CSTR系统的数学模型表达式得到：
步骤1.5.将步骤1.2的式子和步骤1.3的式子代入步骤1.4的式子得到：
其中，由平滑输出以及其导数构成，表示系统输入。
4.根据权利要求2所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，利用积分方法与比例增益k1和k2估计平滑输出量和它的各阶导数为满足：在此式中，k1和k2均为比例增益，是对积分而得到，是对积分而得到，从而系统输入的估计值为：
5.根据权利要求2所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，将系统输入估计值与实际值作差来获得扰动的估计值为
其中，是由积分方法与比例增益k1和k2所获得。
6.根据权利要求1所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，将扰动观测值前馈至CSTR系统控制输入u端，形成闭环回路，从而将CSTR系统的实际扰动补偿掉。
7.根据权利要求1所述的一种基于微分平滑的CSTR系统的扰动观测方法，其特征在于，外界扰动d(t)为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910456527.8A CN110262231B (zh) | 2019-05-29 | 2019-05-29 | 基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910456527.8A CN110262231B (zh) | 2019-05-29 | 2019-05-29 | 基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110262231A true CN110262231A (zh) | 2019-09-20 |
| CN110262231B CN110262231B (zh) | 2022-09-16 |
Family
ID=67915881
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910456527.8A Active CN110262231B (zh) | 2019-05-29 | 2019-05-29 | 基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110262231B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030191728A1 (en) * | 2002-03-27 | 2003-10-09 | Kulkarni Bhaskar Dattatray | Performance of artificial neural network models in the presence of instrumental noise and measurement errors |
| US20070078238A1 (en) * | 2005-09-30 | 2007-04-05 | Chevron Phillips Chemical Company, Lp | Multiple component feed methods and systems |
| CN105404144A (zh) * | 2014-08-20 | 2016-03-16 | 上海交通大学 | 连续搅拌反应釜的多模型自适应控制方法及系统 |
| CN105893654A (zh) * | 2016-03-11 | 2016-08-24 | 中南大学 | 一种一阶连续搅拌釜式反应器的鲁棒预测控制方法 |
| CN107732959A (zh) * | 2017-09-18 | 2018-02-23 | 上海电力学院 | 用于分布式光储并网系统的非线性微分平滑前馈控制方法 |
| CN109491245A (zh) * | 2018-10-30 | 2019-03-19 | 江苏大学 | 一种cstr系统的扰动补偿控制方法 |
-
2019
- 2019-05-29 CN CN201910456527.8A patent/CN110262231B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030191728A1 (en) * | 2002-03-27 | 2003-10-09 | Kulkarni Bhaskar Dattatray | Performance of artificial neural network models in the presence of instrumental noise and measurement errors |
| US20070078238A1 (en) * | 2005-09-30 | 2007-04-05 | Chevron Phillips Chemical Company, Lp | Multiple component feed methods and systems |
| CN105404144A (zh) * | 2014-08-20 | 2016-03-16 | 上海交通大学 | 连续搅拌反应釜的多模型自适应控制方法及系统 |
| CN105893654A (zh) * | 2016-03-11 | 2016-08-24 | 中南大学 | 一种一阶连续搅拌釜式反应器的鲁棒预测控制方法 |
| CN107732959A (zh) * | 2017-09-18 | 2018-02-23 | 上海电力学院 | 用于分布式光储并网系统的非线性微分平滑前馈控制方法 |
| CN109491245A (zh) * | 2018-10-30 | 2019-03-19 | 江苏大学 | 一种cstr系统的扰动补偿控制方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王维: "连续搅拌反应釜(CSTR)控制方法研究", 《中国优秀博硕士学位论文全文数据库(硕士)》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110262231B (zh) | 2022-09-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104267616B (zh) | 一种时滞系统的自抗扰控制系统的设计及整定方法 | |
| CN108303895B (zh) | 一种电液位置伺服系统多模型鲁棒自适应控制方法 | |
| CN103064282A (zh) | 非线性参数变化模型辨识方法（npv） | |
| CN110764417B (zh) | 一种基于闭环辨识模型的线性二次型最优动态前馈-反馈pid控制系统及其控制方法 | |
| CN100462877C (zh) | 工业过程中非方系统的解耦控制方法 | |
| CN111781835B (zh) | 一种镇定二阶惯性加纯滞后系统的线性自抗扰控制器设计方法 | |
| CN105807632A (zh) | 一种基于改进大林算法的加热炉温度控制器设计方法 | |
| CN105676915A (zh) | 基于自抗扰技术的生产线喷涂烘干过程的温度控制方法 | |
| CN105316613A (zh) | 一种基于变时滞偏差校正技术的锌层厚度控制方法和系统 | |
| CN109519957B (zh) | 一种超超临界锅炉闭环燃烧优化控制方法 | |
| Kurtz et al. | State and disturbance estimation for nonlinear systems affine in the unmeasured variables | |
| CN105676645A (zh) | 一种基于函数型权rbf-arx模型的双回路水箱液位预测控制方法 | |
| CN105867125B (zh) | 炼化装置耦合单元的优化控制方法及装置 | |
| CN110262231B (zh) | 基于微分平滑的连续搅拌釜式反应器系统的扰动观测方法 | |
| CN117608197B (zh) | 一种基于自适应鲁棒滑模技术的低温风洞总温控制方法 | |
| CN113159647B (zh) | 一种基于delta算子的二级化学反应器故障估计方法 | |
| CN116627028B (zh) | 交联电缆生产线控制方法 | |
| CN107703761B (zh) | 一种气动调节阀粘滞特性参数的估计方法 | |
| CN113671833A (zh) | 不确定重复控制系统协调设计方法及系统 | |
| CN109491245B (zh) | 一种cstr系统的扰动补偿控制方法 | |
| JPWO2006088072A1 (ja) | モデル構造のパラメータ決定方法、パラメータ決定装置、制御装置および温度調節器 | |
| CN117130274A (zh) | 一种用于实时混合试验的多模型补偿控制方法 | |
| CN109725526B (zh) | 一种多变量半自适应预测控制方法 | |
| CN109884884A (zh) | 一种系统控制品质的调整方法及相关装置 | |
| CN112776817B (zh) | 履带式车辆、控制方法、控制器和计算机可读存储介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240624 Address after: 518000 1104, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province Patentee after: Shenzhen Hongyue Enterprise Management Consulting Co.,Ltd. Country or region after: China Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301 Patentee before: JIANGSU University Country or region before: China |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240626 Address after: 274900 Coal Chemical Park (Dongguantun), Juye County, Heze City, Shandong Province Patentee after: Juye Fengrui Chemical Co.,Ltd. Country or region after: China Address before: 518000 1104, Building A, Zhiyun Industrial Park, No. 13, Huaxing Road, Henglang Community, Longhua District, Shenzhen, Guangdong Province Patentee before: Shenzhen Hongyue Enterprise Management Consulting Co.,Ltd. Country or region before: China |