CN109782587B - MISO different-factor compact-format model-free control method - Google Patents

MISO different-factor compact-format model-free control method Download PDF

Info

Publication number
CN109782587B
CN109782587B CN201910103040.1A CN201910103040A CN109782587B CN 109782587 B CN109782587 B CN 109782587B CN 201910103040 A CN201910103040 A CN 201910103040A CN 109782587 B CN109782587 B CN 109782587B
Authority
CN
China
Prior art keywords
miso
control
control method
factor
different
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.)
Active
Application number
CN201910103040.1A
Other languages
Chinese (zh)
Other versions
CN109782587A (en
Inventor
卢建刚
陈晨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201910103040.1A priority Critical patent/CN109782587B/en
Publication of CN109782587A publication Critical patent/CN109782587A/en
Application granted granted Critical
Publication of CN109782587B publication Critical patent/CN109782587B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

The invention discloses a MISO different-factor compact format model-free control method, which aims at the limitation of the existing MISO compact format model-free control method adopting the same-factor structure, namely: aiming at the limitation that different control inputs in control input vectors can only adopt penalty factors with the same value and step factors with the same value at the moment k, the MISO compact-format model-free control method adopting the different-factor structure is provided, and the penalty factors with different values and/or the step factors with different values can be adopted aiming at different control inputs in the control input vectors at the moment k, so that the control problem that the characteristics of each control channel are different in complex objects such as a strong nonlinear MISO system and the like can be solved. Compared with the existing control method, the method has the advantages of higher control precision, better stability and wider applicability.

Description

MISO different-factor compact-format model-free control method
Technical Field
The invention belongs to the field of automatic control, and particularly relates to a MISO (multiple input single output) different-factor compact-format model-free control method.
Background
Controlled objects in industries such as oil refining, petrifaction, chemical industry, pharmacy, food, papermaking, water treatment, thermal power, metallurgy, cement, rubber, machinery, electricity and the like comprise reactors, rectifying towers, machines, equipment, devices, production lines, workshops and factories, wherein a plurality of controlled objects are MISO (Multiple Input and Single Output) systems. The control of high precision, high stability and high applicability of the MISO system is realized, and the method has important significance for energy conservation, consumption reduction, quality improvement and efficiency improvement of the industry. However, the control of MISO systems, especially of strongly non-linear MISO systems, has long been a significant challenge in the field of automation control.
The existing control methods of the MISO system comprise a MISO compact format model-free control method. The MISO compact format model-free control method is a novel data-driven control method, does not depend on any mathematical model information of a controlled object, only depends on input and output data measured by the MISO controlled object in real time to analyze and design a controller, and has the advantages of simplicity, low calculation burden, strong robustness and good application prospect. The theoretical basis of the MISO compact-format model-free control method is proposed by Houzhong and Jinshangtai in the 'model-free adaptive control-theory and application' (scientific publishing house, 2013, page 95) of the union of Houzhong and Jinshangtai, and the control algorithm is as follows:
Figure BDA0001966075900000011
where u (k) is a control input vector at time k, and u (k) is [ u (k) ]1(k),…,um(k)]TM is the total number of control inputs (m is a positive integer greater than 1); e (k) is the error at time k; phi (k) is an estimated value of a pseudo Jacobian matrix of the MISO system at the moment k, | phi (k) | is a 2 norm of the matrix phi (k); λ is a penalty factor; ρ is the step factor.
The existing MISO compact format model-free control method adopts the same factor structure, that is to say: at time k, for different control inputs u in a control input vector u (k)1(k),…,um(k) Only the same value of penalty factor λ and the same value of step factor ρ can be used. When the existing MISO same-factor compact-format model-free control method is applied to complex objects such as a strong nonlinear MISO system, ideal control effect is often difficult to realize due to different control channel characteristics, and popularization and application of the MISO compact-format model-free control method are restricted.
Therefore, in order to break the application bottleneck of the existing MISO same-factor compact-format model-free control method, the invention provides the MISO different-factor compact-format model-free control method.
Disclosure of Invention
In order to solve the problems in the background art, an object of the present invention is to provide a MISO different-factor compact-format model-free control method, comprising:
when the controlled object is a Multiple Input and Single Output (MISO) system, the MISO different-factor compact format model-free control method calculates the ith control Input u at the time ki(k) The mathematical formula of (a) is as follows:
Figure BDA0001966075900000021
wherein k is a positive integer; i represents the ith of the total number of the MISO system control inputs, i is a positive integer, i is more than or equal to 1 and less than or equal to m, m is the total number of the MISO system control inputs, and m is a positive integer more than 1; u. ofi(k) The ith control input at time k; e (k) is the error at time k; phi (k) is the estimated value of the pseudo Jacobian matrix of the MISO system at the moment of k, phij,i(k) Is the ith row and the ith column element of the matrix phi (k), and | phi (k) | is the 2 norm of the matrix phi (k); lambda [ alpha ]iA penalty factor for the ith control input; rhoiStep size factor for the ith control input;
aiming at the MISO system, the MISO different-factor compact format model-free control method traverses the value of i through the positive integer interval [1, m]All the values in the time point k are calculated to obtain the control input vector u (k) ═ u1(k),…,um(k)]T
The MISO different-factor compact-format model-free control method has different-factor characteristics; the different factor characteristic refers to that for any two unequal positive integers i and x in a positive integer interval [1, m ], during the control of the MISO system by adopting the control method, at least one of the following two inequalities is true:
λi≠λx;ρi≠ρx
while adopting the above technical scheme, the present invention can also adopt or combine the following further technical schemes:
the k time error e (k) is obtained by calculation through an error calculation function; the independent variables of the error calculation function include an output expected value and an output actual value.
The error calculation function employs e (k) y*(k) -y (k), wherein y*(k) Outputting expected values for k time, and outputting actual values for k time; or using e (k) ═ y*(k +1) -y (k), wherein y*(k +1) outputting an expected value at the moment k + 1; or using e (k) ═ y (k) — y*(k) (ii) a Or using e (k) ═ y (k) — y*(k+1)。
The controlled object comprises a reactor, a rectifying tower, a machine, equipment, a device, a production line, a workshop and a factory.
The hardware platform for operating the control method comprises any one or any combination of an industrial control computer, a single chip microcomputer controller, a microprocessor controller, a field programmable gate array controller, a digital signal processing controller, an embedded system controller, a programmable logic controller, a distributed control system, a field bus control system, an industrial Internet of things control system and an industrial Internet control system.
The MISO different-factor compact-format model-free control method provided by the invention can adopt penalty factors with different values or step factors with different values aiming at different control inputs in the control input vector, and can solve the control problem of different control channel characteristics in complex objects such as a strong nonlinear MISO system and the like. Therefore, compared with the existing MISO same-factor compact-format model-free control method, the MISO different-factor compact-format model-free control method provided by the invention has higher control precision, better stability and wider applicability.
Drawings
FIG. 1 is a functional block diagram of the present invention;
FIG. 2 is a control effect diagram of a two-input single-output MISO system adopting the MISO different-factor compact-format model-free control method of the present invention;
FIG. 3 is a control input curve for a two-input single-output MISO system employing the MISO different-factor compact-format model-free control method of the present invention;
FIG. 4 is a diagram of the control effect of a two-input single-output MISO system using the existing MISO same-factor compact-format model-free control method;
fig. 5 is a control input curve of a two-input single-output MISO system when the existing MISO same-factor compact-format model-free control method is adopted.
Detailed Description
The invention is further described with reference to the following figures and specific examples.
Fig. 1 shows a schematic block diagram of the present invention. For a MISO system with m control inputs (m is a positive integer greater than 1), a MISO different-factor compact-format model-free control method is adopted for control; needleFor the ith control input ui(k) (i 1, …, m), determining a MISO differential factor compact format model-free control method for calculating ui(k) All parameters in the mathematical formula (2) contain a penalty factor lambdaiStep factor piThe value of (d); recording the current time as k time; will output the expected value y*(k) The difference between the output actual value y (k) and the output actual value y (k) is used as the k time error e (k); based on the k time error vector e (k) ═ e1(k),…,en(k)]TAnd a penalty factor lambdaiStep factor piThe ith control input u at the k moment is calculated by adopting an MISO different-factor compact-format model-free control methodi(k) (ii) a For the MISO system, the MISO different-factor compact format model-free control method traverses the value of i through the positive integer interval [1, m]All the values in the time point k are calculated to obtain the control input vector u (k) ═ u1(k),…,um(k)]T(ii) a And (4) after the control input vector u (k) acts on the controlled object, obtaining an output actual value of the controlled object at the later moment, and then carrying out MISO (MISO) different-factor compact-format model-free control at the later moment.
The following are specific examples of the present invention.
The two-input single-output MISO system adopted by the controlled object has the complex characteristic of strong nonlinearity, and belongs to a typical difficult-to-control MISO system:
Figure BDA0001966075900000051
desired value y of system output*(k) The following were used:
y*(k)=(-1)round((k-1)/100)
in a specific embodiment, m is 2.
For the above specific examples, two sets of experiments were performed for comparative validation. To more clearly compare the control performance of the two sets of tests, the Root Mean Square Error (RMSE) was used as a control performance evaluation index:
Figure BDA0001966075900000052
wherein e (k) ═ y*(k)-y(k),y*(k) The expected value is output for time k, and y (k) the actual value is output for time k. The smaller the value of RMSE (e), the more the actual output value y (k) and the desired output value y*(k) Is smaller overall and the control performance is better.
The hardware platform for operating the control method of the invention adopts an industrial control computer.
For the first set of experiments: by adopting the MISO different-factor compact-format model-free control method, the penalty factor lambda of the 1 st control input is set1Step factor p 3.0011.05; setting penalty factor lambda of 2 nd control input2Step factor p 4.0021.35; then, the two-input single-output MISO system is controlled, wherein FIG. 2 is a control effect diagram of the output, and FIG. 3 is a control input curve; from the control performance evaluation index, the rmse (e) output in fig. 2 is 0.3325.
For the second set of experiments: directly adopting the existing MISO same-factor compact-format model-free control method, setting a penalty factor lambda to be 3.00 and setting a step factor rho to be 1.00, and controlling the two-input single-output MISO system, wherein FIG. 4 is an output control effect diagram, and FIG. 5 is a control input curve; from the control performance evaluation index, the rmse (e) output in fig. 4 is 0.3396.
The comparison results of the two groups of test control performance evaluation indexes are listed in table 1, the results of the first group of tests adopting the control method of the invention are superior to the results of the second group of tests adopting the existing MISO same-factor compact format model-free control method, wherein the improvement effect can be found to be remarkable by comparing the control performance indexes of the two groups of tests, and the MISO different-factor compact format model-free control method provided by the invention has higher control precision, better stability and wider applicability.
TABLE 1 control Performance comparison
Figure BDA0001966075900000061
Further, the following three points should be particularly pointed out:
(1) controlled objects in industries such as oil refining, petrifaction, chemical industry, pharmacy, food, papermaking, water treatment, thermal power, metallurgy, cement, rubber, machinery, electricity and the like comprise reactors, rectifying towers, machines, equipment, devices, production lines, workshops and factories, wherein a plurality of controlled objects are MISO systems, have complex characteristics of strong nonlinearity and are typical difficult-to-control objects; for example, a continuous stirred tank reactor CSTR commonly used in oil refining, petrochemical, chemical, pharmaceutical industries, etc. is a common two-input single-output MISO system, in which two control inputs are a feed flow and a cooling water flow, respectively, and an output is a reaction temperature; when the chemical reaction has strong exothermic effect, the MISO system of the continuous stirred reactor CSTR has complex characteristics of strong nonlinearity, and is a typical uncontrollable object. In the above specific embodiment, the two-input single-output MISO system adopted by the controlled object also has the complex characteristic of strong nonlinearity, and belongs to a MISO system which is particularly difficult to control; the invention can realize high-precision, high-stability and high-applicability control on the controlled object, and the control method can also realize high-precision, high-stability and high-applicability control on a reactor, a rectifying tower, a machine, equipment, a device, a production line, a workshop, a factory and other complex MISO systems.
(2) In the above embodiment, the hardware platform for operating the control method of the present invention is an industrial control computer; in practical application, any one or any combination of a single-chip microcomputer controller, a microprocessor controller, a field programmable gate array controller, a digital signal processing controller, an embedded system controller, a programmable logic controller, a distributed control system, a field bus control system, an industrial internet of things control system and an industrial internet control system can be selected as a hardware platform for operating the control method according to specific conditions.
(3) In the above-described embodiment, the desired value y will be output*(k) The difference from the output actual value y (k) is used as the k time error e (k), i.e. e (k) y*(k) -y (k), only said errorOne method in the difference calculation function; the expected value y can also be output at the moment k +1*The difference between (k +1) and the time k output y (k) is taken as the error e (k), i.e., e (k) y*(k +1) -y (k); the error calculation function may also employ other calculation methods in which the arguments include an output desired value and an output actual value, for example,
Figure BDA0001966075900000071
for the controlled object of the above embodiment, good control effects can be achieved by using the different error calculation functions.
The above-described embodiments are intended to illustrate the present invention, but not to limit the present invention, and any modifications, equivalents, improvements, etc. made within the spirit of the present invention and the scope of the claims fall within the scope of the present invention.

Claims (5)

  1. The MISO different-factor compact format model-free control method is characterized by comprising the following steps:
    when the controlled object is a Multiple Input and Single Output (MISO) system, the MISO different-factor compact format model-free control method calculates the ith control Input u at the time ki(k) The mathematical formula of (a) is as follows:
    Figure FDA0002725128280000011
    wherein k is a positive integer; i represents the ith of the total number of the MISO system control inputs, i is a positive integer, i is more than or equal to 1 and less than or equal to m, m is the total number of the MISO system control inputs, and m is a positive integer more than 1; u. ofi(k) The ith control input at time k; e (k) is the error at time k; phi (k) is the estimated value of the pseudo Jacobian matrix of the MISO system at the moment of k, phij,i(k) Is the ith row and the ith column element of the matrix phi (k), and | phi (k) | is the 2 norm of the matrix phi (k); lambda [ alpha ]iA penalty factor for the ith control input; rhoiStep size factor for the ith control input;
    aiming at the MISO system, the MISO different-factor compact format model-free control method traverses the value of i through the positive integer interval [1, m]All the values in the time point k are calculated to obtain the control input vector u (k) ═ u1(k),…,um(k)]T
    The MISO different-factor compact-format model-free control method has different-factor characteristics; the different factor characteristic refers to that for any two unequal positive integers i and x in a positive integer interval [1, m ], during the control of the MISO system by adopting the control method, at least one of the following two inequalities is true:
    λi≠λx;ρi≠ρx
  2. 2. the MISO different-factor compact-format model-free control method according to claim 1, characterized in that: the k time error e (k) is obtained by calculation through an error calculation function; the independent variables of the error calculation function include an output expected value and an output actual value.
  3. 3. The MISO different-factor compact-format model-free control method according to claim 2, characterized in that: the error calculation function employs e (k) y*(k) -y (k), wherein y*(k) Outputting expected values for k time, and outputting actual values for k time; or using e (k) ═ y*(k +1) -y (k), wherein y*(k +1) outputting an expected value at the moment k + 1; or using e (k) ═ y (k) — y*(k) (ii) a Or using e (k) ═ y (k) — y*(k+1)。
  4. 4. The MISO different-factor compact-format model-free control method according to claim 1, characterized in that: the controlled object comprises a reactor, a rectifying tower, a machine, equipment, a device, a production line, a workshop and a factory.
  5. 5. The MISO different-factor compact-format model-free control method according to claim 1, characterized in that: the hardware platform for operating the control method comprises any one or any combination of an industrial control computer, a single chip microcomputer controller, a microprocessor controller, a field programmable gate array controller, a digital signal processing controller, an embedded system controller, a programmable logic controller, a distributed control system, a field bus control system, an industrial Internet of things control system and an industrial Internet control system.
CN201910103040.1A 2019-02-01 2019-02-01 MISO different-factor compact-format model-free control method Active CN109782587B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910103040.1A CN109782587B (en) 2019-02-01 2019-02-01 MISO different-factor compact-format model-free control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910103040.1A CN109782587B (en) 2019-02-01 2019-02-01 MISO different-factor compact-format model-free control method

Publications (2)

Publication Number Publication Date
CN109782587A CN109782587A (en) 2019-05-21
CN109782587B true CN109782587B (en) 2020-12-22

Family

ID=66504087

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910103040.1A Active CN109782587B (en) 2019-02-01 2019-02-01 MISO different-factor compact-format model-free control method

Country Status (1)

Country Link
CN (1) CN109782587B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7142626B2 (en) * 2003-05-30 2006-11-28 George Shu-Xing Cheng Apparatus and method of controlling multi-input-single-output systems
CN105759603A (en) * 2016-03-23 2016-07-13 东北大学 Voltage transformation circuit control system and method based on automatic optimizing model-free controller
CN107633692A (en) * 2017-09-29 2018-01-26 河南理工大学 A kind of city expressway Entrance ramp MFA control method
CN108287471A (en) * 2017-12-04 2018-07-17 浙江大学 Inclined methods of self-tuning of the format Non-Model Controller based on systematic error of MIMO

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7142626B2 (en) * 2003-05-30 2006-11-28 George Shu-Xing Cheng Apparatus and method of controlling multi-input-single-output systems
CN105759603A (en) * 2016-03-23 2016-07-13 东北大学 Voltage transformation circuit control system and method based on automatic optimizing model-free controller
CN107633692A (en) * 2017-09-29 2018-01-26 河南理工大学 A kind of city expressway Entrance ramp MFA control method
CN108287471A (en) * 2017-12-04 2018-07-17 浙江大学 Inclined methods of self-tuning of the format Non-Model Controller based on systematic error of MIMO

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Controller design and experiment for autothermal reforming of methanol in miniature reactor;Jiangang Lu 等;《ISATransactions》;20140104;全文 *
无模型控制器参数在线自整定研究;李雪园;《中国优秀硕士学位论文全文数据库信息科技辑》;20180815;第47-69页 *
无模型自适应控制器参数的实验研究;张晨晖 等;《自动化应用》;20151231;全文 *

Also Published As

Publication number Publication date
CN109782587A (en) 2019-05-21

Similar Documents

Publication Publication Date Title
CN111522230B (en) MIMO different factor tight format model-free control method
CN111522233B (en) Parameter self-tuning MIMO different factor full-format model-free control method
CN111522232B (en) MIMO different factor full-format model-free control method
CN111522229B (en) Parameter self-tuning MIMO different factor offset format model-free control method
CN109782586B (en) Parameter self-tuning MISO different-factor compact-format model-free control method
CN114240243B (en) Rectifying tower product quality prediction method and device based on dynamic system identification
CN113489014A (en) Rapid and flexible full-pure embedded type power system optimal power flow evaluation method
CN111522235B (en) MIMO different factor tight format model-free control method with self-setting parameters
CN109782587B (en) MISO different-factor compact-format model-free control method
CN109581865B (en) MISO different-factor partial-format model-free control method
CN109814390B (en) MISO different-factor full-format model-free control method
CN109814388B (en) Parameter self-tuning MISO different-factor offset format model-free control method
CN110298767B (en) Thermal power plant time sequence variable abnormity monitoring method and system
CN107038489B (en) Multi-target unit combination optimization method based on improved NBI method
CN109581992B (en) Parameter self-tuning MISO different-factor full-format model-free control method
CN111522231B (en) MIMO different factor offset format model-free control method
CN113553232B (en) Technology for carrying out unsupervised anomaly detection on operation and maintenance data through online matrix image
CN111582751B (en) Time-weighted electricity purchasing risk early warning method
CN117606591A (en) Reservoir dispatching data deviation processing method
CN111752152A (en) LSTM neural network-based parameter self-tuning method for SISO (SISO) compact format model-free controller

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