CN110989349A - Multivariate system model identification method - Google Patents

Multivariate system model identification method Download PDF

Info

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
Application number
CN201911263779.5A
Other languages
Chinese (zh)
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.)
Shanghai University of Electric Power
Shanghai Electric Power University
University of Shanghai for Science and Technology
Original Assignee
Shanghai Electric Power University
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 Shanghai Electric Power University filed Critical Shanghai Electric Power University
Priority to CN201911263779.5A priority Critical patent/CN110989349A/en
Publication of CN110989349A publication Critical patent/CN110989349A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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

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

Multivariate system model identification method
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)
Figure BDA0002304102470000021
Wherein N is the total number of sampled data.
Step T3, calculating the function formula of the optimized performance index:
Figure BDA0002304102470000022
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:
Figure BDA0002304102470000031
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:
Figure BDA0002304102470000071
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:
Figure BDA0002304102470000072
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)
Figure BDA0002304102470000081
Where N is the total number of sampled data.
Step T3, calculating the function formula of the optimized performance index:
Figure BDA0002304102470000082
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.
Figure BDA0002304102470000091
Figure BDA0002304102470000101
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)
Figure FDA0002304102460000011
Wherein N is the total number of the sampling data;
step T3, calculating the function formula of the optimized performance index:
Figure FDA0002304102460000012
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:
Figure FDA0002304102460000021
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.
CN201911263779.5A 2019-12-05 2019-12-05 Multivariate system model identification method Pending CN110989349A (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (8)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN110532674B (en) Method for measuring temperature of boiler furnace of coal-fired power plant
CN108897286B (en) Fault detection method based on distributed nonlinear dynamic relation model
CN108647808B (en) Production parameter optimization prediction method, device, equipment and storage medium
CN110413601B (en) Generator data screening method based on Gauss naive Bayes and prediction error method
US20200175435A1 (en) System and method for optimizing combustion of boiler
CN105955026B (en) Fuzzy PID control method, device and system
CN110597070B (en) Method for identifying model parameters of thermal power generating unit system
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
WO2020258826A1 (en) Industrial equipment operation, maintenance and optimization method and system based on complex network model
CN108960309B (en) Dynamic process monitoring method based on RBF neural network autocorrelation elimination
CN105629766B (en) Multivariable time delay system discrimination method based on step test
CN115510904A (en) Boiler heating surface ash deposition monitoring method based on time sequence prediction
CN112765746B (en) Turbine blade top gas-thermal performance uncertainty quantification system based on polynomial chaos
US11475186B2 (en) Apparatus and method for deriving boiler combustion model
CN110989349A (en) Multivariate system model identification method
Muravyova et al. Development of the intellectual complex for parallel work of steam boilers
CN111709585A (en) Air conditioner load prediction method and device and storage medium
CN115526433A (en) Power plant reheat flue gas baffle operation prediction method based on integrated hybrid model
CN110245400B (en) Identification method of oxygen object model of boiler combustion system
CN112348696A (en) Heat supply unit peak regulation upper limit evaluation method and system based on BP neural network
CN112670976A (en) Multi-energy flow calculation and analysis of electric-gas coupling system under attack of controllable load
Łabęda-Grudziak Identification of dynamic system additive models by KDD methods
CN113705932B (en) Short-term load prediction method and device

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