CN110989349A - Multivariate system model identification method - Google Patents
Multivariate system model identification method Download PDFInfo
- Publication number
- CN110989349A CN110989349A CN201911263779.5A CN201911263779A CN110989349A CN 110989349 A CN110989349 A CN 110989349A CN 201911263779 A CN201911263779 A CN 201911263779A CN 110989349 A CN110989349 A CN 110989349A
- Authority
- CN
- China
- Prior art keywords
- transfer function
- multivariate
- function matrix
- input
- multivariate system
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012546 transfer Methods 0.000 claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims abstract description 37
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 22
- 230000005284 excitation Effects 0.000 claims abstract description 18
- 239000002245 particle Substances 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 32
- 238000005070 sampling Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 241000544061 Cuculus canorus Species 0.000 claims description 3
- 239000004576 sand Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 2
- 238000012795 verification Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000010200 validation analysis Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
Images
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
The invention provides a multivariate system model identification method, aiming at the problem that a multivariate system is difficult to identify in engineering application, for the multivariate system with M inputs, the identification data under M groups of irrelevant excitation conditions are collected together, a transfer function matrix G(s) of the multivariate system is set, a particle swarm optimization algorithm is adopted to carry out model identification on a selected multivariate system structure model, the parameters of the transfer function matrix G(s) are solved to obtain a transfer function matrix G(s), and finally, the best fitted equivalent model of the tested system is obtained. The verification proves that the identification method is simple and easy to implement, high in identification precision and high in engineering application value.
Description
Technical Field
The invention belongs to the field of control science and engineering, and relates to a multivariate system model identification method.
Background
With the development of science and technology and social economy, modern industrial equipment is increasingly large and complex, and the production process consisting of a plurality of links is generally coupled and associated among the links. This coupling and association is manifested in that a certain input variable of the system will affect a plurality of output variables simultaneously, or a certain output variable will be affected by a plurality of input variables. This coupling and correlation has become a key difficult factor affecting multivariable system modeling and multivariable system control.
Currently, in the actual implementation of multivariate system identification engineering, a univariate identification method is also commonly used, which uses only one set of data for identification. Because the multiple input quantities of the multivariable system are not all artificially controllable, the result of using only the univariate identification method for the multivariable system identification is inevitably more failure than success, and the identification precision is difficult to improve.
In addition, the number of the manually controllable input quantities is often a small number of the total number of the input quantities of the multivariate system, which means that most of the identification data which can be recorded is multivariate excitation response data, multivariate excitation signals are different from univariate excitation signals, and whether the observed multivariate excitation signals are effective for system model identification needs to be judged before identification calculation.
Disclosure of Invention
In order to solve the above problems, the invention provides an effective identification method in engineering application aiming at the problem that multivariable system is difficult to identify, which adopts the following technical scheme:
the invention provides a multivariate system model identification method, which is characterized by comprising the following steps:
step S1, defining M dimension input vector of multivariable system as U (S) and Q dimension output vector Y (S), setting transfer function matrix G (S) of multivariable system, obtaining input and output relation Y (S) G (S) U (S), and determining transfer function matrix G (S) by carrying out model identification on multivariable system;
and step S2, solving parameters of a transfer function matrix G (S) by using an intelligent optimization algorithm on the M groups of data groups of input vectors and output vectors of the multivariable system under irrelevant excitation to obtain the transfer function matrix G (S), wherein the intelligent optimization algorithm comprises a particle swarm optimization algorithm, a differential evolution algorithm and a cuckoo algorithm.
In step S2, the process of solving the parameters of the transfer function matrix g (S) by using the intelligent optimization algorithm includes the following steps:
and step T1, taking the simulation program corresponding to the intelligent optimization algorithm as an identification program, and initializing.
Step T2, obtaining response data under the same excitation input through simulation calculation aiming at the transfer function matrix G(s)Wherein N is the total number of sampled data.
Step T3, calculating the function formula of the optimized performance index:
step T4, adopting an identification procedure, adjusting the model parameters to:
{ai,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,nij}
{bi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,mij}
wherein n isijFor the denominator order of the transfer function, mijIs the molecular order of the transfer function.
And step T5, repeating the steps from T2 to T4 until the preset counting times are reached.
And step T6, outputting the adjusted model parameters to obtain a transfer function matrix G(s).
Further, in the multivariate system model identification method provided by the present invention, the step S1 comprises the following steps:
step S1.1, let the input vector of the multivariate system be M-dimensional, and define the input vector of the multivariate system as U (S) ═ U1(s) U2(s) … UM(s)]T,Ui(s) (i ═ 1,2, …, M) is the ith input to the multivariate system, the output of the multivariate system is given by dimension Q, and the output vector of the multivariate system is defined as Y(s) ═ Y1(s) Y2(s) … YQ(s)]T,Yj(s) (j ═ 1,2, …, Q) is the jth output of the multivariable system;
step S1.2, setting a transfer function matrix of the multivariable system as G (S), and obtaining the input and output relation of the multivariable system as follows:
Y(s)=G(s)U(s)
wherein G(s) is a transfer function matrix of the multivariate system, and is developed as follows:
the transfer function matrix g(s) is determined by model identification of the multivariate system.
Further, in the multivariate system model identification method provided by the present invention, the step S2 comprises the following steps:
s2.1, selecting M groups of data groups of input vectors and output vectors of the multivariate system under irrelevant excitation, and assuming that the sampling time is TsAnd the total number of the sampling data is N, and the acquired data group is as follows:
{ui,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
{yi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
and S2.2, taking the data groups of the input vectors and the output vectors of the M groups of collected multivariable systems under irrelevant excitation as data to be identified, and solving parameters of a transfer function matrix G (S) by adopting an intelligent optimization algorithm to obtain a transfer function matrix G (S).
Action and Effect of the invention
According to the practical method for identifying the multivariate system model effectively applied to engineering, which is provided by the invention, aiming at the problem that the multivariate system is difficult to identify, for the multivariate system with M inputs, the identification data under M groups of irrelevant excitation conditions are gathered together, and the model identification is carried out on the selected multivariate system structure model by adopting an intelligent optimization algorithm, so that the best fitted equivalent model of the tested system is obtained. The verification proves that the identification method is simple and easy to implement, high in identification precision and high in engineering application value.
Drawings
FIG. 1 is a schematic diagram of a transfer function matrix of a multivariate system model of an embodiment of the invention;
FIG. 2 is a model of a low temperature reheater system of an embodiment of the present invention;
FIG. 3a) is a block diagram of the field operational data of 10:00-11:23 on 31 days 3/2015 of the distributed control system database of an embodiment of the invention, i.e., operational data 1;
FIG. 3b) is the field operational data of 15:30-16:53 on 31 days 3/2015 of the distributed control system database 2015 of an embodiment of the invention, i.e., operational data 2;
FIG. 3c) is the field operational data for distributed control system database 2015, 3, 31, 22:10-23:33, i.e., operational data 3, of an embodiment of the present invention;
FIG. 3d) is the field operation data of the centralized and distributed control system database 2016, 07:00-08:23, 1 month 4, of the present invention, i.e., operation data 4;
FIG. 4 shows the result of the multivariate system identification method of the present invention for the operation data 1-3 in the embodiment of the present invention;
FIG. 5 shows the result of identifying the running data 1 by a univariate identification method according to the embodiment of the present invention;
FIG. 6 is a model validation of the multivariate system identification method of the invention for operational data 4 in an embodiment of the invention;
FIG. 7 is a model verification for the operation data 4 identified by the univariate identification method according to the embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made with reference to the accompanying drawings and examples
< example >
FIG. 2 is a model of a low temperature reheater system of an embodiment of the present invention.
The 660MW once-through boiler of a certain power plant is taken as an embodiment object, and the structure diagram is shown in FIG. 2. Divide the reheat air temperature model into an inlet temperature T13 transfer function channels of feed water flow D and baffle opening degree F are subjected to system identification, and the output is the outlet temperature T of the low-temperature reheater2。
Screening out four sections according to field operation data of a plurality of days stored in a distributed control system databaseAs shown in the figure, 3a) is the field operation data of 10:00-11:23 on 31 days 3 and 31 months in 2015, 3b) is the field operation data of 15:30-16:53 on 31 days 3 and 31 months in 2015, 3c) is the field operation data of 22:10-23:33 on 31 days 3 and 31 months in 2015, and 3d) is the field operation data of 07:00-08:23 on 1 day 1 and 1 month in 2016. The sampling period is 5s and the sampling duration is 5000 s. Wherein: t is2Is the outlet temperature, T, of the low-temperature reheater1Is the inlet temperature of the low temperature reheater; d is water supply flow; f is the opening degree of the flue damper.
Fig. 3a) is the field operational data of 10:00-11:23 on 31 days 3/2015 of the distributed control system database 2015, i.e., operational data 1, according to an embodiment of the invention.
Fig. 3b) is the field operational data of 15:30-16:53 on 31 days 3/2015 of the distributed control system database 2015, i.e. operational data 2, according to an embodiment of the invention.
Fig. 3c) is the field operational data of 22:10-23:33 on 31/3/2015 of the distributed control system database 2015, i.e., operational data 3, according to an embodiment of the invention.
Fig. 3d) is the field operation data of 07:00-08:23 on 1/4/2016 of the distributed control system database of the embodiment of the present invention, i.e., operation data 4.
The data shown in fig. 3a), 3b), 3c) are used for model recognition and the data shown in fig. 3d) are used for model verification. According to the multivariate system model identification method provided by the embodiment, the following steps are carried out:
step S1, defining M-dimensional input vector and Q-dimensional output vector y (S) of the multivariate system, setting the transfer function matrix g (S) of the multivariate system, obtaining the input/output relational expression y (S) ═ g (S) u (S), and determining the transfer function matrix g (S) by performing model identification on the multivariate system, which is the specific steps from step S1.1 to S1.2.
Step S1.1, let the input vector of the multivariate system be M-dimensional, and define the input vector of the multivariate system as U (S) ═ U1(s) U2(s) … UM(s)]T,Ui(s) (i ═ 1,2, …, M) is the ith input to the multivariate system, the output of the multivariate system is given by dimension Q, and the output vector of the multivariate system is defined as Y(s) ═ Y1(s) Y2(s) … YQ(s)]T,Yj(s) (j ═ 1,2, …, Q) is the jth output of the multivariable system;
step S1.2, the defined input and output relationship of the multivariable system is as follows:
Y(s)=G(s)U(s)
wherein G(s) is expanded to:
FIG. 1 is a schematic diagram of a transfer function matrix of a multivariate system model according to an embodiment of the invention.
As shown in fig. 1, the multivariate system is modeled to determine the transfer function matrix g(s).
The most widely applied transfer function model in the process control of the power plant is generally a multi-order inertia link, and the transfer function of the model is as follows:
wherein: k is an amplification factor; t is a time constant; n is the order of the system.
And step S2, selecting M groups of data groups of input vectors and output vectors of the multivariable system under irrelevant excitation, solving parameters of a transfer function matrix G (S) by using an intelligent optimization algorithm to obtain a transfer function matrix G (S), and the specific steps are as described in steps S2.1-S2.2.
S2.1, selecting M groups of data groups of input vectors and output vectors of the multivariate system under irrelevant excitation, and assuming that the sampling time is TsAnd the total number of the sampling data is N, and the acquired data group is as follows:
{ui,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
{yi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
step S2.2, using the data sets of the input vectors and the output vectors of the M groups of multivariate systems under irrelevant excitation as data to be identified, and solving parameters of a transfer function matrix g (S) by using an intelligent optimization algorithm to obtain a transfer function matrix g (S), wherein the intelligent optimization algorithm includes a particle swarm optimization algorithm, a differential evolution algorithm, and a cuckoo algorithm, and the embodiment uses the particle swarm optimization algorithm, and the specific steps are as described in steps T1 to T6.
And step T1, taking a simulation program corresponding to the particle swarm optimization algorithm as an identification program, and initializing.
Step T2, obtaining response data under the same excitation input through simulation calculation aiming at the transfer function matrix G(s)Where N is the total number of sampled data.
Step T3, calculating the function formula of the optimized performance index:
step T4, adopting a particle swarm optimization algorithm to adjust model parameters:
{ai,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,nij}
{bi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,mij}
wherein n isijFor the denominator order of the transfer function, mijIs the molecular order of the transfer function.
And step T5, repeating the steps from T2 to T4 until the preset counting times are reached.
And step T6, outputting the adjusted model parameters to obtain a transfer function matrix G(s).
FIG. 4 shows the result of the multivariate system identification method of the present invention for the operation data 1-3 in the embodiment of the present invention.
As shown in fig. 4, a) operation data 1, b) operation data 2, and c) operation data 3 in the figure represent the results of identification of the outlet temperature of the low-temperature reheater in the data shown in fig. 3a), 3b), and 3c), respectively, and the dotted line in the figure represents the output of the model obtained by the identification, and the solid line represents the actual output of the low-temperature reheater. The identified model parameters are shown in table 1.
Parameter(s) | G1(s) | G2(s) | G3(s) |
K | 2.295 | -0.042 | 0.439 |
T | 409.478 | 163.685 | 401.581 |
n | 4 | 1 | 3 |
TABLE 1
FIG. 5 shows the result of the identification of the operation data 1 by the univariate identification method according to the embodiment of the present invention.
In contrast, as shown in fig. 5, the line represents the output of the identified model, and the solid line represents the actual output of the unit. The identified model parameters are shown in table 2.
TABLE 2
As shown in fig. 4 and 5, the response curve of the output value of the model identified by the univariate identification method and the method of the present invention is very close to the response curve of the true value of the actual unit, but the parameters of the model identified by the univariate identification method and the method of the present invention are very different.
Analyzing the model parameter data shown in table 1, it can be known that the outlet temperature of the reheater increases with the increase of the inlet steam temperature, decreases with the increase of the feedwater flow, and increases with the increase of the baffle opening, which completely conforms to the dynamic characteristics between the outlet temperature of the reheater and the three inputs obtained by the mechanism analysis; the model parameter data shown in table 2 is just contrary to the reheat air temperature dynamics obtained by the mechanism analysis, and therefore, the model parameters identified in table 2 are qualitatively incorrect.
FIG. 6 is a model validation of the multivariate system identification method of the present invention for operational data 4 in an embodiment of the present invention.
As shown in fig. 6, the output response curve of the identification model identified by the method of the present invention can be well fitted to the actual output response curve, and the mean square error is 0.9878.
FIG. 7 is a model verification for the operation data 4 identified by the univariate identification method according to the embodiment of the present invention.
As shown in fig. 7, the output curve of the identification model obtained by identification using the univariate identification method is greatly different from the actual output curve of the unit, and the mean square error is 3.2002. Therefore, the model identified by the method has higher reliability.
Examples effects and effects
According to the multivariate system model identification method provided by the invention, identification data under M groups of uncorrelated excitation conditions are collected together, and then an intelligent optimization algorithm is adopted to carry out model identification on the selected multivariate system structure model. In the embodiment, three groups of operation data are identified, the obtained model parameter data conform to the dynamic characteristics between the outlet temperature of the reheater and three inputs obtained through mechanism analysis, and in addition, the verification of one group of operation data proves that the method is simple and easy to implement, high in identification precision and high in engineering application value.
The above embodiments only describe the implementation and operation results of the practical method for multi-variable system identification provided by the present invention, but the present invention is not limited to the above embodiments, and the method provided by the present invention is effective for other multi-variable systems.
Claims (4)
1. A multivariate system model identification method is used for carrying out model identification on a multivariate system, and is characterized by comprising the following steps:
step S1, defining M-dimensional input vectors and Q-dimensional output vectors of the multivariate system as u (S) and y (S), setting a transfer function matrix g (S) of the multivariate system to obtain an input-output relational expression y (S) ═ g (S) u (S), and determining the transfer function matrix g (S) by performing the model identification on the multivariate system;
step S2, solving parameters of the transfer function matrix g (S) by using an intelligent optimization algorithm on M groups of data sets of the input vector and the output vector of the multivariate system under uncorrelated excitation to obtain the transfer function matrix g (S).
In step S2, the process of solving the parameters of the transfer function matrix g (S) by using the intelligent optimization algorithm includes the following steps:
step T1, taking the simulation program corresponding to the intelligent optimization algorithm as an identification program, and initializing;
step T2, obtaining response data under the same excitation input through simulation calculation aiming at the transfer function matrix G(s)Wherein N is the total number of the sampling data;
step T3, calculating the function formula of the optimized performance index:
step T4, adopting the identification procedure to adjust model parameters:
{ai,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,nij}
{bi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,mij}
wherein n isijFor the denominator order of the transfer function, mijIs the molecular order of the transfer function;
step T5, repeating the steps from T2 to T4 until reaching the preset counting times;
and T6, outputting the adjusted model parameters to obtain the transfer function matrix G(s).
2. The multivariate system model identification method as defined in claim 1, wherein the step S1 comprises the steps of:
step S1.1, assuming that the input vector of the multivariate system is M-dimensional, defining the input vector of the multivariate system as U (S) ═ U1(s) U2(s) … UM(s)]T,Ui(s) (i ═ 1,2, …, M) is the ith input to the multivariable system, assuming that the output quantities of the multivariable system are in dimension Q, the output vector defining the multivariable system is Y(s) ═ Y1(s) Y2(s) … YQ(s)]T,Yj(s) (j ═ 1,2, …, Q) is the jth output of the multivariable system;
s1.2, setting a transfer function matrix of the multivariable system as G (S), and obtaining the input and output relation of the multivariable system as follows:
Y(s)=G(s)U(s)
wherein G(s) is expanded to:
the model identification of the multivariate system is the determination of the transfer function matrix g(s).
3. The multivariate system model identification method as defined in claim 1, wherein the step S2 comprises the steps of:
step S2.1, assuming a sampling time T for the data sets of the input and output vectors of the multivariate system model under uncorrelated excitationsAnd the total number of the sampling data is N, and the acquired data group is as follows:
{ui,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
{yi,j,k,i=1,2,…,M;j=1,2,…,Q;k=1,2,…,N}
s2.2, using the acquired data groups of the input vectors and the output vectors of the M groups of multivariable systems under irrelevant excitation as data to be identified, and solving parameters of the transfer function matrix G (S) by adopting an intelligent optimization algorithm to obtain the transfer function matrix G (S).
4. The multivariate system model identification method as defined in claim 1, wherein:
the intelligent optimization algorithm comprises a particle swarm optimization algorithm, a differential evolution algorithm and a cuckoo algorithm.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911263779.5A CN110989349A (en) | 2019-12-05 | 2019-12-05 | Multivariate system model identification method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911263779.5A CN110989349A (en) | 2019-12-05 | 2019-12-05 | Multivariate system model identification method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110989349A true CN110989349A (en) | 2020-04-10 |
Family
ID=70092294
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911263779.5A Pending CN110989349A (en) | 2019-12-05 | 2019-12-05 | Multivariate system model identification method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110989349A (en) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6097410A (en) * | 1983-10-31 | 1985-05-31 | Idemitsu Petrochem Co Ltd | Identifying method of multivariable system process |
US6381504B1 (en) * | 1996-05-06 | 2002-04-30 | Pavilion Technologies, Inc. | Method for optimizing a plant with multiple inputs |
CN101924533A (en) * | 2010-07-19 | 2010-12-22 | 浙江工业大学 | Multivariable time-lag parameter estimation method based on FIR (Finite Impulse Response) model identification |
CN103760768A (en) * | 2013-12-30 | 2014-04-30 | 广东电网公司电力科学研究院 | Generator set multivariable system identification method based on normal operating data |
CN106066608A (en) * | 2016-05-30 | 2016-11-02 | 上海电力学院 | The modelling by mechanism of multivariate continuous process model merges discrimination method with Experimental modeling |
CN106773702A (en) * | 2017-01-03 | 2017-05-31 | 上海电力学院 | The multiple excitation discrimination method of Multivariable Linear continuous system |
CN107085633A (en) * | 2017-04-12 | 2017-08-22 | 华侨大学 | The device and method of multiple spot vibratory response frequency domain prediction based on SVMs |
CN108954491A (en) * | 2018-07-30 | 2018-12-07 | 东南大学 | A kind of control method of photo-thermal medium temperature offset-type electric boiler heating system |
-
2019
- 2019-12-05 CN CN201911263779.5A patent/CN110989349A/en active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6097410A (en) * | 1983-10-31 | 1985-05-31 | Idemitsu Petrochem Co Ltd | Identifying method of multivariable system process |
US6381504B1 (en) * | 1996-05-06 | 2002-04-30 | Pavilion Technologies, Inc. | Method for optimizing a plant with multiple inputs |
CN101924533A (en) * | 2010-07-19 | 2010-12-22 | 浙江工业大学 | Multivariable time-lag parameter estimation method based on FIR (Finite Impulse Response) model identification |
CN103760768A (en) * | 2013-12-30 | 2014-04-30 | 广东电网公司电力科学研究院 | Generator set multivariable system identification method based on normal operating data |
CN106066608A (en) * | 2016-05-30 | 2016-11-02 | 上海电力学院 | The modelling by mechanism of multivariate continuous process model merges discrimination method with Experimental modeling |
CN106773702A (en) * | 2017-01-03 | 2017-05-31 | 上海电力学院 | The multiple excitation discrimination method of Multivariable Linear continuous system |
CN107085633A (en) * | 2017-04-12 | 2017-08-22 | 华侨大学 | The device and method of multiple spot vibratory response frequency domain prediction based on SVMs |
CN108954491A (en) * | 2018-07-30 | 2018-12-07 | 东南大学 | A kind of control method of photo-thermal medium temperature offset-type electric boiler heating system |
Non-Patent Citations (3)
Title |
---|
汪凤泉等: "多变量转子系统逆重复序列正交多输入在线辨识方法", 《振动、测试与诊断》 * |
袁世通: "1000MW超超临界机组建模理论与方法的研究", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅱ辑》 * |
鹿振宇等: "耦合多新息随机梯度辨识方法性能分析", 《控制与决策》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109992921B (en) | On-line soft measurement method and system for thermal efficiency of boiler of coal-fired power plant | |
US20200175435A1 (en) | System and method for optimizing combustion of boiler | |
CN108647808B (en) | Production parameter optimization prediction method, device, equipment and storage medium | |
CN107194625B (en) | Wind power plant wind curtailment electric quantity evaluation method based on neural network | |
CN110413601B (en) | Generator data screening method based on Gauss naive Bayes and prediction error method | |
CN110597070B (en) | Method for identifying model parameters of thermal power generating unit system | |
CN105955026B (en) | Fuzzy PID control method, device and system | |
CN108960309B (en) | Dynamic process monitoring method based on RBF neural network autocorrelation elimination | |
CN109165418B (en) | Room temperature measuring method based on household calorimeter data | |
CN110941186B (en) | Steam temperature control optimization method based on neural network and universal gravitation search algorithm | |
CN110442911B (en) | High-dimensional complex system uncertainty analysis method based on statistical machine learning | |
CN110276460A (en) | Industrial equipment O&M and optimization method and system based on complex network model | |
CN116306221A (en) | Air pipeline system control method based on LSTM prediction model | |
US20200320237A1 (en) | Apparatus and method for deriving boiler combustion model | |
CN111931436A (en) | Burner nozzle air quantity prediction method based on numerical simulation and neural network | |
CN115510904A (en) | Boiler heating surface ash deposition monitoring method based on time sequence prediction | |
CN111709585A (en) | Air conditioner load prediction method and device and storage medium | |
CN110989349A (en) | Multivariate system model identification method | |
Muravyova et al. | Development of the intellectual complex for parallel work of steam boilers | |
CN112765746A (en) | Turbine blade top gas-thermal performance uncertainty quantification system based on polynomial chaos | |
CN101529347A (en) | Method and device for the identification of a delay-susceptible control path, control device, and computer program product | |
CN115081301B (en) | Carbon emission dynamic evolution method based on mixed PSO-MKPLS | |
Li et al. | The application of immune genetic algorithm in PID parameter optimization for level control system | |
CN115526433A (en) | Power plant reheat flue gas baffle operation prediction method based on integrated hybrid model | |
CN105117530B (en) | Method for identifying parameters of steam turbine and speed regulating system thereof by combination of thickness and thickness regulation |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20200410 |
|
WD01 | Invention patent application deemed withdrawn after publication |