CN101872182A - Batch process monitoring method based on recursive non-linear partial least square - Google Patents
Batch process monitoring method based on recursive non-linear partial least square Download PDFInfo
- Publication number
- CN101872182A CN101872182A CN 201010180651 CN201010180651A CN101872182A CN 101872182 A CN101872182 A CN 101872182A CN 201010180651 CN201010180651 CN 201010180651 CN 201010180651 A CN201010180651 A CN 201010180651A CN 101872182 A CN101872182 A CN 101872182A
- Authority
- CN
- China
- Prior art keywords
- data
- model
- matrix
- batch
- square
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a batch process monitoring method based on recursive non-linear partial least squares. The traditional method adopts the fixed statistical model and control limits, and the model can not be updated promptly to track process changes, thereby resulting in failure and errors in reporting. The method comprises the following steps: establishing a batch process monitoring model based on the non-linear partial least squares according to a process database, calculating control limits, and applying the established model to industrial process data collected on-line to monitor a batch process. When a new batch of data is obtained, the model and the control limits are updated in a recursive method so as to adapt to changes of the batch process. The method provided by the invention overcomes the defect of the traditional nonlinear multivariable statistical process monitoring method that the time varying problem can not be dealt with, and the model parameters and control limits can be updated with the new batch of data in order to adapt to changes of the batch process, thereby improving the monitoring performance.
Description
Technical field
The invention belongs to areas of information technology, relate to a kind of batch process monitoring method based on recursive non-linear partial least square.
Background technology
Along with the popularization of quick manufacturing technology, be applicable to that the batch process of production short run high value added product more and more comes into one's own.In batch process, a lot of quality index can not on-line measurement, normally after a batch of end, judges last product quality quality according to product sampling analysis value.In order better to control product quality, need set up the process monitoring model to batch process, control operation variable according to on-line measurement is monitored product quality, thereby determine whether production status is normal, can in time take its corresponding measures when unusual service condition occurring, reduce and produce down time, the security of safeguards system and reliability to greatest extent.Along with robotization, computer network and development of database, factory can directly obtain a large amount of real-time running datas from production run.But to from observation data, realize assessment, exceed slip-stick artist or operator's limit of power the process operation situation.Do not rely on the analytic model of object based on the multivariate statistics process monitoring technique of data-driven, and make full use of factory's rich data resource, effectively problems such as big, the data dimension height of deal with data amount, data collinearity are used widely industrial.Batch process has stronger non-linear usually, and as time passes, Properties of Objects and working point all may change.Traditional non-linear multivariate statistics course monitoring method all is to adopt fixing statistical model and control limit, thereby when process feature or operating conditions change, can not in time carry out the variation of model modification tracing process, thereby produce the situation of failing to report and reporting by mistake.Therefore, press for a kind of can the non-linear batch process monitoring method that can handle time variation again of processing procedure.
Summary of the invention
Purpose of the present invention is exactly the weak point at the existing process monitoring technique, and a kind of batch process monitoring method based on recursive non-linear partial least square is provided.This method has remedied the deficiency that traditional non-linear multivariate statistics course monitoring method can not be handled time variation, can utilize new data to upgrade model parameter and control limit, thereby the variation of procedure of adaptation characteristic has improved monitoring performance.
The inventive method adopts means such as data acquisition, process identification, data-driven, at first based on the process monitoring model of process database foundation based on non-linear partial least square, calculation control is limit simultaneously, then the model of setting up is applied to the real-time industrial process data of online acquisition, process is monitored.After obtaining new data, adopt recurrence method that model and control limit are upgraded, thus the variation of the procedure of adaptation.
The concrete steps of the inventive method are:
Step (1) is set up non-linear partial least square (NLPLS) process monitoring model based on process database, simultaneously the calculation control limit.Concrete grammar is:
A. by data collector gatherer process service data,, be used to set up model with the process operation data of gathering sample set as data-driven.In the batch process, each batch process repeatedly repeatability is produced, its data acquisition than the continuous process data acquisition many one dimension " batch " element, have sequential property.Therefore the data acquisition of batch production process is with three-dimensional data formation formula (batch * time * variable) expression, and X (all process datas that the expression of I * J * K) collects, wherein, I represents a batch number, and J represents the sample number, and K represents the process variable number; Quality variable generally obtains by off-line analysis when process finishes, and each batch obtains M quality variable assay value, and I batch quality variable formation matrix Y (I * M).
During modeling, at first with three-dimensional data piece X along time-axis direction cutting batch and variable data piece, each data block is horizontal to the right successively, regards the data of each batch as a data sample, forms a new two-dimensional matrix X (I * JK).Then with matrix X as input matrix, Y is used to set up model as output matrix.Wherein, the data of each batch are to being expressed as { x (i) } and { y (i) }, and x (i) expression is imported data for i batch, and y (i) represents i batch of output data.To import data constitutes input matrix X, output data is constituted output matrix Y;
B. set up the non-linear partial least square model based on inputoutput data, method is:
Matrix X and Y are carried out normalized, and making it average is 0, and variance is 1; Then input matrix is listed as expansion, the expansion item is 1 column vector 1 entirely for the latent node output matrix G of radial basis function (RBF) neural network and element, the output g of the corresponding input vector effect of each row latent node down of G wherein, the bias term coefficient that conceals node is 1; Following augmentation input matrix and output matrix are carried out partial least square (PLS) recurrence:
The NLPLS process monitoring model representation that obtains is for [1 X G], Y}:
In the formula, X
EExpression augmentation input matrix, A and H are respectively the weights matrix of coefficients of corresponding original input vector and the latent node output vector of corresponding RBF network, and b is the output offset vector, and T represents transposition.
Unknown parameter in the NLPLS process monitoring model is latent node center vector c, respective width vector σ, weights coefficient matrices A and H, model bias vector b, and these parameters are determined as follows:
1. with the k-means clustering algorithm input data are carried out cluster, obtain latent node center c; This algorithm can be determined optimum cluster centre number, can make cluster centre reasonably be distributed in the data space simultaneously;
2. adopt p neighbour rule to calculate latent node width:
Wherein N is the number of latent node center, c
iBe p nearest latent node center of j latent node center of distance.
3. adopt PLS to return and determine weights coefficient matrices A, H and bias vector b:
Calculate latent node output matrix G according to latent node center that obtains and width, then input matrix is expanded, obtain augmentation input matrix [1 X G].To data to { [1 X G], Y} carry out PLS and return, and obtain PLS model parameter matrix { T, W, P, B, Q}.In order in the model modification of back, to keep all information, extract the order that the characteristic variable number equals augmentation input matrix [1 X G], and the proper vector number a that model kept that finally is used to predict adopts cross validation method to determine, the parameter matrix that obtains is designated as { Ta, Wa, Pa, Ba, Qa} calculates PLS regression coefficient matrix β by them, thereby obtain A, H and b.
C. the control of compute statistics SPE is limit.SPE is a scalar in the value of each batch, and it has portrayed this batch measured value x (i) and y (i) departure degree to principal component model.After obtaining model, be to i batch of its SPE value with the PLS algorithm:
The Model Calculation value of representing i batch process variable,
The Model Calculation value of representing i batch quality variable, SPE control limit Q
aBe calculated as follows:
Wherein:
λ
j=‖t
i‖
2
t
iBe i the proper vector that keeps, C
αBe normal distribution at insolation level be critical value under the α, C
αWith h
0The pivot number that jack per line, k are in the principal component model to be kept, n is whole pivot numbers.
Step (2) is applied to the real-time industrial process data of online acquisition with the NLPLS process monitoring model of setting up, and calculates the SPE value of the data of newly arriving, and limits Q with control
αCompare: if SPE 〉=Q
α, declarative procedure has occurred unusually; If SPE<Q
α, declarative procedure is normal, utilizes the method for step (3) that NLPLS process monitoring model is upgraded, and upgrades control limit Q simultaneously
α
Step (3) utilizes new data in conjunction with original NLPLS process monitoring model, adopts the recursive non-linear partial least square algorithm that model is upgraded, and upgrades control limit Q simultaneously
α, concrete grammar is:
If pass through in the NLPLS process monitoring model that obtains after k-1 batch, the latent node center matrix of RBF network is
The corresponding center vector of each row, the respective width vector is
The width of the corresponding latent node of each element, { W (k-1), P (k-1), B (k-1), Q (k-1) } is PLS model parameter matrix, after k batch end, obtains new input/output variable x (k) and y (k);
A. adopt with step (1) in identical method new data is carried out the data pre-service.Calculate the output vector of the latent node of former NLPLS model, be designated as g (k) for new samples x (k).
B. judge whether to increase new latent node:
If all elements of g (k) all less than setting value, then adds new latent node.New latent node center is taken as x (k), and corresponding width σ adopts the arest neighbors rule to calculate:
σ=z
c-ησ
c
Wherein, z
cBe the distance of x (k) to nearest latent node center, η is overlapping parameter, and span is [0,1], σ
cBe width, thereby obtain new latent node center matrix and width vector from the nearest latent node of x (k):
As follows to parameter matrix P (k-1) and vectorial g (k) expansion simultaneously:
In the formula, 0 is that whole elements all are 0 row vector.
If all elements of g (k) all more than or equal to setting value, does not then need to increase latent node, C
g, σ
g, P, g remain unchanged.
C. x (k) is expanded, obtain augmentation input vector: x
E(k)
T=[1x (k)
TG (k)
T].
D. with new data x
E(k) and y (k) combine with old PLS model parameter matrix, carry out PLS then and return, form is as follows:
According to 3. method of step in the step (1), calculate PLS regression parameter A (k), H (k) and b (k).Preserve new model parameter
For prediction with use during model modification next time.
Step (4) is calculated new control limit Q based on new NLPLS process monitoring model parameter
α, be used for new data, return step (2).
The batch process monitoring method that the present invention proposes based on recursive non-linear partial least square, do not rely on the analytic model of object, make full use of factory's rich data resource, effectively the deal with data amount is big, data dimension height, problems such as data collinearity, this method has remedied the problem of failing to report and reporting by mistake that traditional non-linear multivariate statistics course monitoring method adopts fixing statistical model and control limit to be produced, can utilize new data to upgrade model parameter and control limit, thereby the variation of procedure of adaptation characteristic, can processing procedure is non-linear can handle time variation again, improve monitoring performance.
Embodiment
A kind of batch process monitoring method based on recursive non-linear partial least square, concrete steps are:
Step (1) is set up non-linear partial least square (NLPLS) process monitoring model based on process database, simultaneously the calculation control limit.Concrete grammar is:
A. by data collector gatherer process service data,, be used to set up model with the process operation data of gathering sample set as data-driven.In the batch process, each batch process repeatedly repeatability is produced, its data acquisition than the continuous process data acquisition many one dimension " batch " element, have sequential property.Therefore the data acquisition of batch production process is with three-dimensional data formation formula (batch * time * variable) expression, and X (all process datas that the expression of I * J * K) collects, wherein, I represents a batch number, and J represents the sample number, and K represents the process variable number; Quality variable generally obtains by off-line analysis when process finishes, and each batch obtains M quality variable assay value, and I batch quality variable formation matrix Y (I * M).
During modeling, at first with three-dimensional data piece X along time-axis direction cutting batch and variable data piece, each data block is horizontal to the right successively, regards the data of each batch as a data sample, forms a new two-dimensional matrix X (I * JK).Then with matrix X as input matrix, Y is used to set up model as output matrix.Wherein, the data of each batch are to being expressed as { x (i) } and { y (i) }, and x (i) expression is imported data for i batch, and y (i) represents i batch of output data.To import data constitutes input matrix X, output data is constituted output matrix Y;
B. set up the non-linear partial least square model based on inputoutput data, method is:
Matrix X and Y are carried out normalized, and making it average is 0, and variance is 1; Then input matrix is listed as expansion, the expansion item is 1 column vector 1 entirely for the latent node output matrix G of radial basis function (RBF) neural network and element, the output g of the corresponding input vector effect of each row latent node down of G wherein, the bias term coefficient that conceals node is 1; Following augmentation input matrix and output matrix are carried out partial least square (PLS) recurrence:
The NLPLS process monitoring model representation that obtains is for [1 X G], Y}:
In the formula, X
EExpression augmentation input matrix, A and H are respectively the weights matrix of coefficients of corresponding original input vector and the latent node output vector of corresponding RBF network, and b is the output offset vector, and T represents transposition.
Unknown parameter in the NLPLS process monitoring model is latent node center vector c, respective width vector σ, weights coefficient matrices A and H, model bias vector b, and these parameters are determined as follows:
1. with the k-means clustering algorithm input data are carried out cluster, obtain latent node center c;
2. adopt p neighbour rule to calculate latent node width:
Wherein N is the number of latent node center, c
iBe p nearest latent node center of j latent node center of distance.
3. adopt PLS to return and determine weights coefficient matrices A, H and bias vector b:
Calculate latent node output matrix G according to latent node center that obtains and width, then input matrix is expanded, obtain augmentation input matrix [1 X G].To data to { [1 X G], Y} carry out PLS and return, and obtain PLS model parameter matrix { T, W, P, B, Q}.Extract the order that the characteristic variable number equals augmentation input matrix [1 X G], and the proper vector number a that model kept that finally is used to predict adopts cross validation method to determine, the parameter matrix that obtains is designated as { Ta, Wa, Pa, Ba, Qa}, calculate PLS regression coefficient matrix β by them, thereby obtain A, H and b.
C. the control of compute statistics SPE is limit.SPE is a scalar in the value of each batch, and it has portrayed this batch measured value x (i) and y (i) departure degree to principal component model.After obtaining model, be to i batch of its SPE value with the PLS algorithm:
The Model Calculation value of representing i batch process variable,
The Model Calculation value of representing i batch quality variable, SPE control limit Q
aBe calculated as follows:
Wherein:
λ
j=‖t
i‖
2
t
iBe i the proper vector that keeps, C
αBe normal distribution at insolation level be critical value under the α, C
αWith h
0The pivot number that jack per line, k are in the principal component model to be kept, n is whole pivot numbers.
Step (2) is applied to the real-time industrial process data of online acquisition with the NLPLS process monitoring model of setting up, and calculates the SPE value of the data of newly arriving, and limits Q with control
αCompare: if SPE 〉=Q
α, declarative procedure has occurred unusually; If SPE<Q
α, declarative procedure is normal, utilizes the method for step (3) that NLPLS process monitoring model is upgraded, and upgrades control limit Q simultaneously
α
Step (3) utilizes new data in conjunction with original NLPLS process monitoring model, adopts the recursive non-linear partial least square algorithm that model is upgraded, and upgrades control limit Q simultaneously
α, concrete grammar is:
If pass through in the NLPLS process monitoring model that obtains after k-1 batch, the latent node center matrix of RBF network is
The corresponding center vector of each row, the respective width vector is
The width of the corresponding latent node of each element, { W (k-1), P (k-1), B (k-1), Q (k-1) } is PLS model parameter matrix, after k batch end, obtains new input/output variable x (k) and y (k);
A. adopt with step (1) in identical method new data is carried out the data pre-service.Calculate the output vector of the latent node of former NLPLS model, be designated as g (k) for new samples x (k).
B. judge whether to increase new latent node:
If all elements of g (k) all less than setting value, then adds new latent node.New latent node center is taken as x (k), and corresponding width σ adopts the arest neighbors rule to calculate:
σ=z
c-ησ
c
Wherein, z
cBe the distance of x (k) to nearest latent node center, η is overlapping parameter, and span is [0,1], σ
cBe width, thereby obtain new latent node center matrix and width vector from the nearest latent node of x (k):
As follows to parameter matrix P (k-1) and vectorial g (k) expansion simultaneously:
In the formula, 0 is that whole elements all are 0 row vector.
If all elements of g (k) is all more than or equal to setting value, C
g, σ
g, P, g remain unchanged.
C. x (k) is expanded, obtain augmentation input vector: x
E(k)
T=[1x (k)
TG (k)
T].
D. with new data x
E(k) and y (k) combine with old PLS model parameter matrix, carry out PLS then and return, form is as follows:
According to 3. method of step in the step (1), calculate PLS regression parameter A (k), H (k) and b (k).Preserve new model parameter
For prediction with use during model modification next time.
Step (4) is calculated new control limit Q based on new NLPLS process monitoring model parameter
α, be used for new data, return step (2).
Claims (1)
1. batch process monitoring method based on recursive non-linear partial least square is characterized in that the concrete steps of this method are:
Step (1) is set up non-linear partial least square process monitoring model based on process database, simultaneously the calculation control limit; Concrete grammar is:
A. by data collector gatherer process service data, with the process operation data of gathering sample set, be used to set up model as data-driven, represent X (all process datas that the expression of I * J * K) collects with three-dimensional data formation formula; Wherein I represents that a batch number, J represent that sample number, K represent the process variable number; Each batch obtains M quality variable assay value, and I batch quality variable formation matrix Y (I * M);
During modeling, at first with three-dimensional data piece X along time-axis direction cutting batch and variable data piece, each data block is horizontal to the right successively, regards the data of each batch as a data sample, forms a new two-dimensional matrix X (I * JK); Then with matrix X as input matrix, Y is used to set up model as output matrix; Wherein, the data of each batch are to being expressed as { x (i) } and { y (i) }, and x (i) expression is imported data for i batch, and y (i) represents i batch of output data; To import data constitutes input matrix X, output data is constituted output matrix Y;
B. set up the non-linear partial least square model based on inputoutput data, method is:
Matrix X and Y are carried out normalized, and making it average is 0, and variance is 1; Then input matrix is listed as expansion, the expansion item is 1 column vector 1 entirely for the latent node output matrix G of radial basis function neural network and element, the output g of the corresponding input vector effect of each row latent node down of G wherein, the bias term coefficient that conceals node is 1; Following augmentation input matrix and output matrix are carried out the partial least square recurrence:
The non-linear partial least square process monitoring model representation that obtains is for [1 X G], Y}:
In the formula, X
EExpression augmentation input matrix, A and H are respectively the weights matrix of coefficients of corresponding original input vector and the latent node output vector of corresponding radial basis function neural network, and b is the output offset vector, and T represents transposition;
Unknown parameter in the non-linear partial least square process monitoring model is latent node center vector c, respective width vector σ, weights coefficient matrices A and H, model bias vector b, and these parameters are determined as follows:
1. with the k-means clustering algorithm input data are carried out cluster, obtain latent node center c;
2. adopt p neighbour rule to calculate latent node width:
Wherein N is the number of latent node center, c
iBe p nearest latent node center of j latent node center of distance;
3. adopt partial least square to return and determine weights coefficient matrices A, H and bias vector b:
Calculate latent node output matrix G according to latent node center that obtains and width, then input matrix is expanded, obtain augmentation input matrix [1 X G]; To data to { [1 X G], Y} carry out partial least square and return, and obtain partial least square model parameter matrix { T, W, P, B, Q}; Extract the order that the characteristic variable number equals augmentation input matrix [1 X G], and the proper vector number a that model kept that finally is used to predict adopts cross validation method to determine, the parameter matrix that obtains is designated as { Ta, Wa, Pa, Ba, Qa}, calculate partial least square regression coefficient matrix β by them, thereby obtain A, H and b;
C. the control of compute statistics SPE is limit; SPE is a scalar in the value of each batch, and it has portrayed this batch measured value x (i) and y (i) departure degree to principal component model; After obtaining model, be to i batch of its SPE value with the partial least square algorithm:
The Model Calculation value of representing i batch process variable,
The Model Calculation value of representing i batch quality variable, SPE control limit Q
aBe calculated as follows:
Wherein:
λ
j=‖t
i‖
2
t
iBe i the proper vector that keeps, C
αBe normal distribution at insolation level be critical value under the α, C
αWith h
0The pivot number that jack per line, k are in the principal component model to be kept, n is whole pivot numbers;
Step (2) is applied to the real-time industrial process data of online acquisition with the non-linear partial least square process monitoring model of setting up, and calculates the SPE value of the data of newly arriving, and limits Q with control
αCompare: if SPE 〉=Q
α, declarative procedure has occurred unusually; If SPE<Q
α, declarative procedure is normal, utilizes the method for step (3) that non-linear partial least square process monitoring model is upgraded, and upgrades control limit Q simultaneously
α
Step (3) utilizes new data in conjunction with original non-linear partial least square process monitoring model, adopts the recursive non-linear partial least square algorithm that model is upgraded, and upgrades control limit Q simultaneously
α, concrete grammar is:
If pass through in the non-linear partial least square process monitoring model that obtains after k-1 batch, the latent node center matrix of radial basis function neural network is
The corresponding center vector of each row, the respective width vector is
The width of the corresponding latent node of each element, { W (k-1), P (k-1), B (k-1), Q (k-1) } is partial least square model parameter matrix, after k batch end, obtains new input/output variable x (k) and y (k);
A. adopt with step (1) in identical method new data is carried out the data pre-service; Calculate the output vector of the latent node of former non-linear partial least square process monitoring model, be designated as g (k) for new samples x (k);
B. judge whether to increase new latent node:
If all elements of g (k) all less than setting value, then adds new latent node; New latent node center is taken as x (k), and corresponding width σ adopts the arest neighbors rule to calculate:
σ=z
c-ησ
c
Wherein, z
cBe the distance of x (k) to nearest latent node center, η is overlapping parameter, and span is [0,1], σ
cBe width, thereby obtain new latent node center matrix and width vector from the nearest latent node of x (k):
As follows to parameter matrix P (k-1) and vectorial g (k) expansion simultaneously:
In the formula, 0 is that whole elements all are 0 row vector;
If all elements of g (k) is all more than or equal to setting value, C
g, σ
g, P, g remain unchanged;
C. x (k) is expanded, obtain augmentation input vector: x
E(k)
T=[1x (k)
TG (k)
T];
D. with new data x
E(k) and y (k) divide least square model parameter matrix to combine with former subordinates, carry out partial least square then and return, form is as follows:
According to 3. method of step in the step (1), calculate partial least square regression parameter A (k), H (k) and b (k); Preserve new model parameter A (k), H (k), b (k), P (k), B (k), Q (k),
For prediction with use during model modification next time;
Step (4) is calculated new control limit Q based on new non-linear partial least square process monitoring model parameter
α, be used for new data, return step (2).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010180651 CN101872182A (en) | 2010-05-21 | 2010-05-21 | Batch process monitoring method based on recursive non-linear partial least square |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010180651 CN101872182A (en) | 2010-05-21 | 2010-05-21 | Batch process monitoring method based on recursive non-linear partial least square |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101872182A true CN101872182A (en) | 2010-10-27 |
Family
ID=42997100
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010180651 Pending CN101872182A (en) | 2010-05-21 | 2010-05-21 | Batch process monitoring method based on recursive non-linear partial least square |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101872182A (en) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102662321A (en) * | 2012-03-23 | 2012-09-12 | 清华大学 | Online updating method of principal component analysis monitoring model |
CN103076741A (en) * | 2013-01-18 | 2013-05-01 | 杭州电子科技大学 | Quadric form control method for non-minimum realization of expansion state space in chemical process |
CN103076744A (en) * | 2013-01-18 | 2013-05-01 | 杭州电子科技大学 | Linear quadric form control method for non-minimum realization of state space in chemical process |
CN103092078A (en) * | 2013-01-07 | 2013-05-08 | 北京中医药大学 | Whole-process optimizing method of multi-stage batch production process |
CN105955214A (en) * | 2016-05-26 | 2016-09-21 | 华东理工大学 | Batch process fault detection method based on sample timing sequence and neighborhood similarity information |
CN108037668A (en) * | 2017-12-28 | 2018-05-15 | 杭州电子科技大学 | A kind of new Chemical Batch Process modeling and monitoring method |
CN109074051A (en) * | 2016-03-28 | 2018-12-21 | 三菱电机株式会社 | Quality management device, method for quality control and quality management program |
CN109101758A (en) * | 2018-09-03 | 2018-12-28 | 江南大学 | Batch process process conditions design method based on T-PLS model |
CN109541940A (en) * | 2018-11-13 | 2019-03-29 | 海南师范大学 | Mix fault tolerant control method based on 2D model multistage batch process constrained predictive |
CN109932908A (en) * | 2019-03-20 | 2019-06-25 | 杭州电子科技大学 | A kind of multidirectional pivot analysis process monitoring method based on alarm reliability fusion |
CN112834752A (en) * | 2021-01-05 | 2021-05-25 | 上海市疾病预防控制中心 | Blood glucose standardized measurement system, method, terminal and medium suitable for large-batch crowd |
CN113096087A (en) * | 2021-04-06 | 2021-07-09 | 清华大学深圳国际研究生院 | Soft measurement method for purity of phosphorus recovery product |
CN116702095B (en) * | 2023-06-01 | 2024-04-05 | 大连海事大学 | Modularized marine ship motion attitude real-time forecasting method |
-
2010
- 2010-05-21 CN CN 201010180651 patent/CN101872182A/en active Pending
Non-Patent Citations (2)
Title |
---|
《中国博士学位论文全文数据库》 2005 李春富 基于数据的软测量建模方法及其应用的研究 摘要、第17-18,40-43,58-62,72-77页 1 , 2 * |
《机电工程》 201003 苏金明等 递推PLS方法及其在过程监控中的应用 第93-95,99页 1 第27卷, 第3期 2 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102662321A (en) * | 2012-03-23 | 2012-09-12 | 清华大学 | Online updating method of principal component analysis monitoring model |
CN102662321B (en) * | 2012-03-23 | 2014-03-12 | 清华大学 | Online updating method of principal component analysis monitoring model |
CN103092078A (en) * | 2013-01-07 | 2013-05-08 | 北京中医药大学 | Whole-process optimizing method of multi-stage batch production process |
CN103092078B (en) * | 2013-01-07 | 2015-05-20 | 北京中医药大学 | Whole-process optimizing method of multi-stage batch production process |
CN103076741A (en) * | 2013-01-18 | 2013-05-01 | 杭州电子科技大学 | Quadric form control method for non-minimum realization of expansion state space in chemical process |
CN103076744A (en) * | 2013-01-18 | 2013-05-01 | 杭州电子科技大学 | Linear quadric form control method for non-minimum realization of state space in chemical process |
CN103076744B (en) * | 2013-01-18 | 2015-11-18 | 杭州电子科技大学 | Chemical process non-minimum realizes state space linear quadric form control method |
CN103076741B (en) * | 2013-01-18 | 2015-11-25 | 杭州电子科技大学 | Chemical process non-minimum realizes extended mode space quadric form control method |
CN109074051A (en) * | 2016-03-28 | 2018-12-21 | 三菱电机株式会社 | Quality management device, method for quality control and quality management program |
CN105955214B (en) * | 2016-05-26 | 2018-10-02 | 华东理工大学 | Batch process fault detection method based on sample time-series and neighbour's affinity information |
CN105955214A (en) * | 2016-05-26 | 2016-09-21 | 华东理工大学 | Batch process fault detection method based on sample timing sequence and neighborhood similarity information |
CN108037668A (en) * | 2017-12-28 | 2018-05-15 | 杭州电子科技大学 | A kind of new Chemical Batch Process modeling and monitoring method |
CN109101758A (en) * | 2018-09-03 | 2018-12-28 | 江南大学 | Batch process process conditions design method based on T-PLS model |
CN109541940A (en) * | 2018-11-13 | 2019-03-29 | 海南师范大学 | Mix fault tolerant control method based on 2D model multistage batch process constrained predictive |
CN109541940B (en) * | 2018-11-13 | 2022-03-29 | 海南师范大学 | Multi-stage intermittent process limited prediction hybrid fault-tolerant control method based on 2D model |
CN109932908A (en) * | 2019-03-20 | 2019-06-25 | 杭州电子科技大学 | A kind of multidirectional pivot analysis process monitoring method based on alarm reliability fusion |
CN109932908B (en) * | 2019-03-20 | 2022-03-01 | 杭州电子科技大学 | Multi-directional principal component analysis process monitoring method based on alarm reliability fusion |
CN112834752A (en) * | 2021-01-05 | 2021-05-25 | 上海市疾病预防控制中心 | Blood glucose standardized measurement system, method, terminal and medium suitable for large-batch crowd |
CN112834752B (en) * | 2021-01-05 | 2024-04-30 | 上海市疾病预防控制中心 | Blood glucose standardized measurement system, method, terminal and medium suitable for large-scale crowd |
CN113096087A (en) * | 2021-04-06 | 2021-07-09 | 清华大学深圳国际研究生院 | Soft measurement method for purity of phosphorus recovery product |
WO2022213958A1 (en) * | 2021-04-06 | 2022-10-13 | 清华大学深圳国际研究生院 | Soft measurement method for purity of phosphorus recovery product |
CN116702095B (en) * | 2023-06-01 | 2024-04-05 | 大连海事大学 | Modularized marine ship motion attitude real-time forecasting method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101872182A (en) | Batch process monitoring method based on recursive non-linear partial least square | |
CN112101480B (en) | Multivariate clustering and fused time sequence combined prediction method | |
CN101872444B (en) | Batch-to-batch optimization method of batch process by combining medium-term correction strategy | |
CN101446831B (en) | Decentralized process monitoring method | |
CN108038300A (en) | Optical fiber state evaluating method based on improved membership function combination neutral net | |
CN105004542A (en) | Online monitoring and fault diagnosing method for mixing and flavouring process of cigarette filament production based on principal component analysis | |
CN101308385A (en) | Non-linear procedure failure testing method based on two-dimensional dynamic kernel principal component analysis | |
CN104699050A (en) | Leaf-shred preparation segment on-line monitoring and fault diagnosing method for cigarette filament treatment driven by data | |
CN104865827B (en) | Oil pumping unit oil extraction optimization method based on multi-working-condition model | |
CN102357934B (en) | Quality monitor soft sensing method based on rubber mixing process | |
CN104808648A (en) | Online and real-time batch process monitoring method based on k nearest neighbor | |
CN101872163B (en) | Process monitoring method based on recursive non-linear partial least squares | |
CN117008479B (en) | Carbon emission optimization control method and system based on biomass gasification furnace | |
CN101187803A (en) | Ammonia converter production optimization method based on data excavation technology | |
CN102601881A (en) | Method for monitoring on-line quality and updating prediction model of rubber hardness | |
CN101446828A (en) | Nonlinear process quality prediction method | |
CN105259318A (en) | Foul smell OU value prediction method and system based on meteorological parameters | |
CN112417627A (en) | Power distribution network operation reliability analysis method based on four-dimensional index system | |
CN116755411A (en) | Industrial equipment data acquisition method and device | |
CN114117919B (en) | Instant learning soft measurement modeling method based on sample collaborative representation | |
CN117556366B (en) | Data abnormality detection system and method based on data screening | |
CN103389360A (en) | Probabilistic principal component regression model-based method for soft sensing of butane content of debutanizer | |
CN114253242A (en) | VPN-based Internet of things cloud equipment data acquisition system | |
CN116842841A (en) | Immersed tunnel service state prediction method based on artificial intelligence | |
CN108345726A (en) | Input variable automatic selecting method based on mutual information and soft instrument error |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20101027 |