CN103675011A - Soft industrial melt index measurement instrument and method of optimal support vector machine - Google Patents
Soft industrial melt index measurement instrument and method of optimal support vector machine Download PDFInfo
- Publication number
- CN103675011A CN103675011A CN201310435358.2A CN201310435358A CN103675011A CN 103675011 A CN103675011 A CN 103675011A CN 201310435358 A CN201310435358 A CN 201310435358A CN 103675011 A CN103675011 A CN 103675011A
- Authority
- CN
- China
- Prior art keywords
- support vector
- training sample
- particle
- sigma
- iter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
The invention discloses a soft industrial melt index measurement instrument and method of an optimal support vector machine. The soft measurement method comprises the following steps: carrying out fuzzy processing on outputs of a plurality of weighting least square support vector machines; and optimizing a whole fuzzy equation system by adopting a particle swarm optimization to obtain an optimal soft measurement result. According to the soft industrial melt index measurement instrument and method of the optimal support vector machine, a field intelligent instrument for measuring easily-measured variables and a control station for measuring operation variables are connected with a DCS (Data Communication System) database; a soft measurement value displayer comprises a soft industrial melt index measurement model of the optimal support vector machine; the DCS database is connected with the input end of the soft measurement model; the output end of the soft industrial melt index measurement model of the optimal support vector machine is connected with the soft melt index measurement value displayer; the soft industrial melt index measurement instrument and method of the optimal support vector machine have the characteristics of on-line parameter optimization, strong noise immunity and good popularization performance, and the automatic updating of the model can be realized.
Description
Technical field
The present invention relates to soft measuring instrument and method, relate in particular to a kind of industrial melting index soft measuring instrument and method of optimum support vector machine.
Background technology
Polypropylene is a kind of hemicrystalline thermoplastics being formed by propylene polymerization, has higher resistance to impact, and engineering properties is tough, and anti-multiple organic solvent and acid and alkali corrosion, be widely used in industry member, is one of usual modal macromolecular material.Melting index (MI) is to determine one of important quality index of the final products trade mark during polypropylene is produced, and it has determined the different purposes of product.Measuring accurately, timely of melting index, to producing and scientific research, has very important effect and directive significance.Yet the on-line analysis of melting index is measured and is still difficult at present accomplish, the in-line analyzer that lacks melting index is a subject matter of restriction polypropylene product quality.MI can only obtain by hand sampling, off-line assay, and analyzes once for general every 2-4 hour, and time lag is large, is difficult to meet the requirement of producing real-time control.
Research work major part about the online forecasting of MI all concentrates on above artificial neural network in recent years, has obtained good effect.But artificial neural network also has the shortcoming of himself, for example the interstitial content of over-fitting, hidden layer and parameter are bad determines.Secondly, the DCS data that industry spot collects also because noise, manual operation error etc. with certain uncertain error, so use the general Generalization Ability of forecasting model of the artificial neural network that determinacy is strong or not.
First nineteen sixty-five U.S. mathematician L.Zadeh has proposed the concept of fuzzy set.Fuzzy logic, in the mode of its problem closer to daily people and meaning of one's words statement, starts to replace adhering to the classical logic that all things can represent with binary item subsequently.Fuzzy logic so far successful Application industry a plurality of fields among, fields such as household electrical appliances, Industry Control.2003, Demirci proposed the concept of fuzzifying equation, and by using fuzzy membership matrix and building a new input matrix with its distortion, the gravity model appoach of then usining in local equation in Anti-fuzzy method show that analytic value is as last output.For the soft measurement of melting index in propylene polymerization production process, consider noise effect and operate miss in industrial processes, can use the fuzzy performance of fuzzy logic to reduce the impact of error on whole forecast precision.
Support vector machine, is introduced in 1998 by Vapnik, due to its good Generalization Ability, is widely used in pattern-recognition, matching and classification problem.Because standard support vector machine is responsive to isolated point and noise, so proposed again afterwards Weighted Least Squares Support Vector Machines.Weighted Least Squares Support Vector Machines can be processed the sample data with noise better than standard support vector machine, is selected as the local equation in fuzzifying equation here.
Particle cluster algorithm, Particle Swarm Optimization, is a kind of a kind of biological intelligence optimizing algorithm of seeking global optimum by imitating Bird Flight behavior being put forward by Kennedy and professor Eberhart, is called for short PSO.This algorithm, by interparticle influencing each other in colony, has reduced searching algorithm and has been absorbed in the risk of locally optimal solution, has good global search performance.Particle cluster algorithm is used to the best parameter group of search weighted least square method supporting vector machine, to reach the object of Optimized model.
Summary of the invention
In order to overcome the deficiency that the measuring accuracy of existing propylene polymerization production process is not high, low to noise sensitivity, promote poor performance, the invention provides a kind of on-line measurement, computing velocity is fast, model upgrades automatically, noise resisting ability strong, promote industrial melting index soft measuring instrument and the method for the optimum support vector machine that performance is good.
A kind of industrial melting index soft measuring instrument of optimum support vector machine, comprise propylene polymerization production process, for measuring the field intelligent instrument of easy survey variable, for measuring the control station of performance variable, the DCS database of store data and melt index flexible measured value display instrument, described field intelligent instrument, control station is connected with propylene polymerization production process, described field intelligent instrument, control station is connected with DCS database, described soft measuring instrument also comprises the industrial melting index soft-sensing model of optimum support vector machine, described DCS database is connected with the input end of the industrial melting index soft-sensing model of described optimum support vector machine, the output terminal of the industrial melting index soft-sensing model of described optimum support vector machine is connected with melt index flexible measured value display instrument, the industrial melting index soft-sensing model of described optimum support vector machine comprises:
Data preprocessing module, for by carrying out pre-service from the model training sample of DCS database input, to training sample centralization, deducts the mean value of sample, then it is carried out to standardization:
Computation of mean values:
Calculate variance:
Standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
Fuzzifying equation module, the training sample X to from data preprocessing module passes the standardization of coming, carries out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1func(μ
ik)X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is m component of corresponding Lagrange multiplier.μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
Particle cluster algorithm is optimized module, and for adopting particle cluster algorithm to be optimized the penalty factor of fuzzifying equation Weighted Least Squares Support Vector Machines local equation and error margin value, specific implementation step is as follows:
7. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
8. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
9. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
10. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
if the individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
As preferred a kind of scheme, the industrial melting index soft-sensing model of described optimum support vector machine also comprises: model modification module, for the online updating of model, regularly off-line analysis data is input in training set, and upgrade fuzzifying equation model.
An industrial melt index flexible measurement method for optimum support vector machine, described flexible measurement method specific implementation step is as follows:
1), to propylene polymerization production process object, according to industrial analysis and Operations Analyst, select performance variable and easily survey variable as the input of model, performance variable and easily survey variable and obtained by DCS database;
2), the model training sample from DCS database input is carried out to pre-service, to training sample centralization, deduct the mean value of sample, then it is carried out to standardization, making its average is 0, and variance is 1.This processing adopts following formula process to complete:
2.1) computation of mean values:
2.2) calculate variance:
2.3) standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
3), to pass the training sample come from data preprocessing module, carry out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1func(μ
ik)X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is m component of corresponding Lagrange multiplier.μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
Wherein,
for the output of fuzzy group k at training sample i.
4), adopt particle cluster algorithm to be optimized the penalty factor of Weighted Least Squares Support Vector Machines local equation in fuzzifying equation and error margin value, specific implementation step is as follows:
7. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
8. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
9. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
10. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
if the individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
As preferred a kind of scheme: described flexible measurement method is further comprising the steps of: 5), regularly off-line analysis data is input in training set, upgrade fuzzifying equation model.
Technical conceive of the present invention is: the important quality index melting index to propylene polymerization production process is carried out online soft sensor, overcome the deficiency that existing polypropylene melting index measurement instrument measuring accuracy is not high, low to noise sensitivity, promote poor performance, introduce particle cluster algorithm fuzzifying equation model is carried out to Automatic Optimal, do not need artificial experience repeatedly to adjust the parameter of Weighted Least Squares Support Vector Machines local equation in fuzzifying equation.This model has following advantage with respect to existing melting index soft-sensing model: (1) has reduced noise and the impact of manual operation error on model prediction precision; (2) strengthened the popularization performance of model, over-fitting has effectively been suppressed; (3) parameter of model is carried out to automatic optimal, improved the stability of model, reduced the possibility that model is absorbed in local optimum.
Beneficial effect of the present invention is mainly manifested in: 1, on-line measurement; 2, on-line parameter Automatic Optimal; 3, model upgrades automatically; 4, anti-noise jamming ability strong, 5, precision is high; 6, Generalization Ability is strong.
Accompanying drawing explanation
Fig. 1 is the industrial melting index soft measuring instrument of optimum support vector machine and the basic structure schematic diagram of method;
Fig. 2 is the industrial melting index soft-sensing model structural representation of optimum support vector machine.
Embodiment
Below in conjunction with accompanying drawing, the invention will be further described.The embodiment of the present invention is used for the present invention that explains, rather than limits the invention, and in the protection domain of spirit of the present invention and claim, any modification and change that the present invention is made, all fall into protection scope of the present invention.
Embodiment 1
With reference to Fig. 1, Fig. 2, a kind of industrial melting index soft measuring instrument of optimum support vector machine, comprise propylene polymerization production process 1, for measuring the field intelligent instrument 2 of easy survey variable, for measuring the control station 3 of performance variable, the DCS database 4 of store data and melt index flexible measured value display instrument 6, described field intelligent instrument 2, control station 3 is connected with propylene polymerization production process 1, described field intelligent instrument 2, control station 3 is connected with DCS database 4, described soft measuring instrument also comprises the soft-sensing model 5 of particle cluster algorithm optimization Weighted Least Squares Support Vector Machines fuzzifying equation, described DCS database 4 is connected with the input end of the industrial melting index soft-sensing model 5 of described optimum support vector machine, the output terminal of the industrial melting index soft-sensing model 5 of described optimum support vector machine is connected with melt index flexible measured value display instrument 6, the industrial melting index soft-sensing model of described optimum support vector machine comprises:
Data preprocessing module, for by carrying out pre-service from the model training sample of DCS database input, to training sample centralization, deducts the mean value of sample, then it is carried out to standardization:
Computation of mean values:
Calculate variance:
Standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
Fuzzifying equation module, the training sample X to from data preprocessing module passes the standardization of coming, carries out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1func(μ
ik)X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is corresponding Lagrange multiplier, μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
Particle cluster algorithm is optimized module, and for adopting particle cluster algorithm to be optimized the penalty factor of fuzzifying equation Weighted Least Squares Support Vector Machines local equation and error margin value, specific implementation step is as follows:
1. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
2. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
3. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
4. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
If the 5. individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
6. judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
As preferred a kind of scheme, the industrial melting index soft-sensing model of described optimum support vector machine also comprises: model modification module, for the online updating of model, regularly off-line analysis data is input in training set, and upgrade fuzzifying equation system model.
According to reaction mechanism and flow process analysis, consider the various factors in polypropylene production process, melting index being exerted an influence, get nine performance variables conventional in actual production process and easily survey variable as modeling variable, have: three strand of third rare feed flow rates, major catalyst flow rate, cocatalyst flow rate, temperature in the kettle, pressure, liquid level, hydrogen volume concentration in still.Table 1 has been listed 9 modeling variablees as soft-sensing model 5 inputs, is respectively liquid level (L) in temperature in the kettle (T), still internal pressure (p), still, the interior hydrogen volume concentration (X of still
v), 3 bursts of propylene feed flow rates (first strand of third rare feed flow rates f1, second strand of third rare feed flow rates f2, the 3rd strand of third rare feed flow rates f3), 2 bursts of catalyst charge flow rates (major catalyst flow rate f4, cocatalyst flow rate f5).Polyreaction in reactor is that reaction mass mixes rear participation reaction repeatedly, so mode input variable relates to the mean value in front some moment of process variable employing of material.The mean value of last hour for data acquisition in this example.Melting index off-line laboratory values is as the output variable of soft-sensing model 5.By hand sampling, off-line assay, obtain, within every 4 hours, analyze and gather once.
Field intelligent instrument 2 and control station 3 are connected with propylene polymerization production process 1, are connected with DCS database 4; Soft-sensing model 5 is connected with DCS database and soft measured value display instrument 6.Field intelligent instrument 2 is measured the easy survey variable that propylene polymerization is produced object, will easily survey variable and be transferred to DCS database 4; Control station 3 is controlled the performance variable that propylene polymerization is produced object, and performance variable is transferred to DCS database 4.In DCS database 4, the variable data of record is as the input of the industrial melting index soft-sensing model 5 of optimum support vector machine, and soft measured value display instrument 6 is for showing the output of the industrial melting index soft-sensing model 5 of optimum support vector machine, i.e. soft measured value.
Table 1: the required modeling variable of industrial melting index soft-sensing model of optimum support vector machine
Variable symbol | Variable implication | Variable symbol | Variable implication |
T | Temperature in the kettle | f1 | First strand of third rare feed flow rates |
p | Pressure in still | f2 | Second strand of third rare feed flow rates |
[0139]
L | Liquid level in still | f3 | The 3rd strand of third rare feed flow rates |
X v | Hydrogen volume concentration in still | f4 | Major catalyst flow rate |
f5 | Cocatalyst flow rate |
The industrial melting index soft-sensing model 5 of optimum support vector machine, comprises following 4 parts:
Data preprocessing module 7, for by carrying out pre-service from the model training sample of DCS database input, to training sample centralization, deducts the mean value of sample, then it is carried out to standardization:
Computation of mean values:
Calculate variance:
Standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1func(μ
ik)X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is corresponding Lagrange multiplier, μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
Particle cluster algorithm is optimized module 9, and for adopting particle cluster algorithm to be optimized the penalty factor of fuzzifying equation Weighted Least Squares Support Vector Machines local equation and error margin value, specific implementation step is as follows:
1. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
2. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
3. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
4. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
If the 5. individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
6. judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
Embodiment 2
With reference to Fig. 1, Fig. 2, a kind of industrial polypropylene producing melt index flexible measurement method of optimizing Weighted Least Squares Support Vector Machines fuzzifying equation model based on particle cluster algorithm, described flexible measurement method concrete methods of realizing is as follows:
1), to propylene polymerization production process object, according to industrial analysis and Operations Analyst, select performance variable and easily survey variable as the input of model, performance variable and easily survey variable and obtained by DCS database;
2), the model training sample from DCS database input is carried out to pre-service, to training sample centralization, deduct the mean value of sample, then it is carried out to standardization, making its average is 0, and variance is 1.This processing adopts following formula process to complete:
2.1) computation of mean values:
2.2) calculate variance:
2.3) standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
3), to pass the training sample come from data preprocessing module, carry out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1func(μ
ik)Xi] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is corresponding Lagrange multiplier, μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
4), adopt particle cluster algorithm to be optimized the penalty factor of Weighted Least Squares Support Vector Machines local equation in fuzzifying equation and error margin value, specific implementation step is as follows:
1. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
2. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
3. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
4. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
If the 5. individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
6. judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
As preferred a kind of scheme: described flexible measurement method is further comprising the steps of: 4), regularly off-line analysis data is input in training set, upgrade fuzzifying equation model.
The concrete implementation step of method of the present embodiment is as follows:
Step 1: to propylene polymerization production process object 1, according to industrial analysis and Operations Analyst, select performance variable and easily survey variable as the input of model.Performance variable and easily survey variable are obtained by DCS database 4.
Step 2: sample data is carried out to pre-service, completed by data preprocessing module 7.
Step 3: set up initial fuzzy equation model 8 based on model training sample data.Input data obtain as described in step 2, and output data are obtained by off-line chemical examination.
Step 4: the local weighted least square method supporting vector machine equation parameter of being optimized initial fuzzy equation model 8 by particle cluster algorithm 9.
Step 5: model modification module 10 is regularly input to off-line analysis data in training set, upgrades fuzzifying equation model, optimizes the soft-sensing model 5 of Weighted Least Squares Support Vector Machines fuzzifying equation model set up based on particle cluster algorithm.
Step 6: melt index flexible measured value display instrument 6 shows the output of the industrial melting index soft-sensing model 5 of optimum support vector machine, completes the demonstration that industrial polypropylene producing melt index flexible is measured.
Claims (2)
1. the industrial melt index flexible of an optimum support vector machine is measured soft measuring instrument, comprise for measuring the field intelligent instrument of easy survey variable, for measuring the control station of performance variable, the DCS database of store data and melt index flexible measured value display instrument, described field intelligent instrument, control station is connected with DCS database, it is characterized in that: described soft measuring instrument also comprises the industrial melting index soft-sensing model of optimum support vector machine, described DCS database is connected with the input end of the industrial melting index soft-sensing model of described optimum support vector machine, the output terminal of the industrial melting index soft-sensing model of described optimum support vector machine is connected with melt index flexible measured value display instrument, the industrial melting index soft-sensing model of described optimum support vector machine comprises:
Data preprocessing module, for by carrying out pre-service from the model training sample of DCS database input, to training sample centralization, deducts the mean value of sample, then it is carried out to standardization:
Computation of mean values:
Calculate variance:
Standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
Fuzzifying equation module, the training sample X to from data preprocessing module passes the standardization of coming, carries out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1 func(μ
ik) X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is m component of corresponding Lagrange multiplier.μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
Particle cluster algorithm is optimized module, and for adopting particle cluster algorithm to be optimized the penalty factor of fuzzifying equation Weighted Least Squares Support Vector Machines local equation and error margin value, specific implementation step is as follows:
1. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
2. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
3. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
4. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
If the 5. individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
6. judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
The industrial melting index soft-sensing model of described optimum support vector machine also comprises:
Model modification module, for the online updating of model, is regularly input to off-line analysis data in training set, upgrades fuzzifying equation model.
2. a flexible measurement method of realizing with the industrial melting index soft measuring instrument of optimum support vector machine as claimed in claim 1, is characterized in that: described flexible measurement method specific implementation step is as follows:
1), to propylene polymerization production process object, according to industrial analysis and Operations Analyst, select performance variable and easily survey variable as the input of model, performance variable and easily survey variable and obtained by DCS database;
2), the model training sample from DCS database input is carried out to pre-service, to training sample centralization, deduct the mean value of sample, then it is carried out to standardization, making its average is 0, and variance is 1.This processing adopts following formula process to complete:
2.1) computation of mean values:
2.2) calculate variance:
2.3) standardization:
Wherein, TX
ibe i training sample, N is number of training,
for the average of training sample, X is the training sample after standardization.σ
xthe standard deviation that represents training sample, σ
2 xthe variance that represents training sample.
3), to pass the training sample come from data preprocessing module, carry out obfuscation.If have c in fuzzifying equation system
*individual fuzzy group, the center of fuzzy group k, j is respectively v
k, v
j, i training sample X after standardization
idegree of membership μ for fuzzy group k
ikfor:
In formula, n is the partitioned matrix index needing in fuzzy classification process, conventionally get and do 2, || || be norm expression formula.
Use above degree of membership value or its distortion to obtain new input matrix, for fuzzy group k, its input matrix is deformed into:
Φ
ik(X
i,μ
ik)=[1 func(μ
ik) X
i] (5)
Func (μ wherein
ik) be degree of membership value μ
ikwarping function, generally get
exp (μ
ik) etc., Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.
Weighted Least Squares Support Vector Machines, as the local equation of fuzzifying equation system, is optimized matching to each fuzzy group.If i target of model training sample is output as O
i, the support vector machine of weighted is equivalent to following quadratic programming problem by conversion fitting problems:
Wherein, R (w, ξ) is the objective function of optimization problem, and minR (w, ξ) is the minimum value of the objective function of optimization problem, and N is number of training, ξ={ ξ
1..., ξ
nslack variable, ξ
ibe i component of slack variable, w is the normal vector of support vector machine lineoid, and b is corresponding side-play amount, and ω
i, i=1 ..., N and γ are respectively weight and the penalty factors of Weighted Least Squares Support Vector Machines,
i component ξ of Weighted Least Squares Support Vector Machines slack variable
ithe estimation of standard deviation, c
1for constant, get 2.5, c here
2for constant, get 3 here, the transposition of subscript T representing matrix, μ
ikrepresent i training sample X after standardization
ifor the degree of membership of fuzzy group k, Φ
ik(X
i, μ
ik) i input variable X of expression
iand the degree of membership μ of fuzzy group k
ikcorresponding new input matrix.By (6) (7) (8) formula, can derive fuzzy group k is output as at training sample i:
Wherein,
for the output of fuzzy group k at training sample i.K<> is the kernel function of Weighted Least Squares Support Vector Machines, and K<> gets linear kernel function here; α
m, m=1 ..., N is m component of corresponding Lagrange multiplier.μ
mkrepresent m training sample X
mfor the degree of membership of fuzzy group k, Φ
mk(X
m, μ
mk) m input variable X of expression
mand the degree of membership μ of fuzzy group k
mkcorresponding new input matrix.
Gravity model appoach in Anti-fuzzy method obtains the output of last fuzzifying equation system:
4), adopt particle cluster algorithm to be optimized the penalty factor of Weighted Least Squares Support Vector Machines local equation in fuzzifying equation and error margin value, specific implementation step is as follows:
1. the Optimal Parameters of determining population is the penalty factor of Weighted Least Squares Support Vector Machines local equation and error margin value, population individual amount popsize, largest loop optimizing number of times iter
max, a p particle initial position r
p, initial velocity v
p, local optimum Lbest
pand the global optimum Gbest of whole population.
2. set optimization aim function, be converted into fitness, each On Local Fuzzy equation is evaluated; By corresponding error function, calculate fitness function, and think that the large particle fitness of error is little, the fitness function of particle p is expressed as:
f
p=1/(E
p+1) (11)
In formula, E
pbe the error function of fuzzifying equation system, be expressed as:
In formula,
the prediction output of fuzzifying equation system, O
itarget output for fuzzifying equation system;
3. according to following formula, speed and the position of each particle upgraded in circulation,
v
p(iter+1)=ω×v
p(iter)+m
1a
1(Lbest
p-r
p(iter))+m
2a
2(Gbest-r
p(iter))
(13)
r
p(iter+1)=r
p(iter)+v
p(iter+1) (14)
In formula, v
prepresent the more speed of new particle p, r
prepresent the more position of new particle p, Lbest
prepresent the more individual optimal value of new particle p, Gbest is the global optimum of whole population, and iter represents cycle index, and ω is the inertia weight in particle cluster algorithm, m
1, m
2corresponding accelerator coefficient, a
1, a
2it is the random number between [0,1];
4. for particle p, if new fitness is greater than original individual optimal value, the individual optimal value of new particle more:
Lbest
p=f
p (15)
If the 5. individual optimal value Lbest of particle p
pbe greater than original population global optimum Gbest, upgrade original population global optimum Gbest:
Gbest=Lbest
p (16)
6. judge whether to meet performance requirement, if so, finish optimizing, obtain the local equation parameter of one group of fuzzifying equation of optimizing; Otherwise return to step 3., continue iteration optimizing, until reach maximum iteration time iter
max.
Described flexible measurement method is further comprising the steps of: 5), regularly off-line analysis data is input in training set, upgrade fuzzifying equation model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310435358.2A CN103675011B (en) | 2013-09-22 | 2013-09-22 | The industrial melt index soft measurement instrument of optimum support vector machine and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310435358.2A CN103675011B (en) | 2013-09-22 | 2013-09-22 | The industrial melt index soft measurement instrument of optimum support vector machine and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103675011A true CN103675011A (en) | 2014-03-26 |
CN103675011B CN103675011B (en) | 2015-09-30 |
Family
ID=50313216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310435358.2A Expired - Fee Related CN103675011B (en) | 2013-09-22 | 2013-09-22 | The industrial melt index soft measurement instrument of optimum support vector machine and method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103675011B (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680428A (en) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | Soft measuring method for denitration control system |
CN108388113A (en) * | 2018-02-07 | 2018-08-10 | 浙江大学 | Least square method supporting vector machine soft-measuring modeling method based on distribution estimation local optimum |
CN108803525A (en) * | 2018-06-28 | 2018-11-13 | 浙江大学 | A kind of quick high-precision propylene polymerization production process optimal soft survey instrument of chaos |
CN108804851A (en) * | 2018-06-28 | 2018-11-13 | 浙江大学 | A kind of high-precision propylene polymerization production process optimal soft survey instrument of chaos gunz optimizing |
CN110378461A (en) * | 2019-06-14 | 2019-10-25 | 上海交通大学 | Error self-adaptation control method and system based on space-time error separate |
CN110750756A (en) * | 2019-10-01 | 2020-02-04 | 深圳市行健自动化股份有限公司 | Method for checking and diagnosing real-time online instrument by optimal support vector machine algorithm |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050197994A1 (en) * | 2004-03-03 | 2005-09-08 | Shigeru Fujii | Intelligent robust control system for motorcycle using soft computing optimizer |
CN1916791A (en) * | 2006-09-12 | 2007-02-21 | 浙江大学 | Method of soft measuring fusion index of producing propylene through polymerization in industrialization |
CN101021723A (en) * | 2006-12-22 | 2007-08-22 | 浙江大学 | Melt index detection fault diagnozing system and method in propylene polymerization production |
CN201035377Y (en) * | 2006-12-22 | 2008-03-12 | 浙江大学 | Failure diagnosis device of melt index detecting in polymerization of propylene produce |
SG141218A1 (en) * | 2003-10-07 | 2008-04-28 | Nanyang Polytechnic | Method for prediction of single nucleotide polymorphisms |
CN102662324A (en) * | 2012-04-28 | 2012-09-12 | 江南大学 | Non-linear model predication control method of tank reactor based on on-line support vector machine |
-
2013
- 2013-09-22 CN CN201310435358.2A patent/CN103675011B/en not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SG141218A1 (en) * | 2003-10-07 | 2008-04-28 | Nanyang Polytechnic | Method for prediction of single nucleotide polymorphisms |
US20050197994A1 (en) * | 2004-03-03 | 2005-09-08 | Shigeru Fujii | Intelligent robust control system for motorcycle using soft computing optimizer |
CN1916791A (en) * | 2006-09-12 | 2007-02-21 | 浙江大学 | Method of soft measuring fusion index of producing propylene through polymerization in industrialization |
CN101021723A (en) * | 2006-12-22 | 2007-08-22 | 浙江大学 | Melt index detection fault diagnozing system and method in propylene polymerization production |
CN201035377Y (en) * | 2006-12-22 | 2008-03-12 | 浙江大学 | Failure diagnosis device of melt index detecting in polymerization of propylene produce |
CN102662324A (en) * | 2012-04-28 | 2012-09-12 | 江南大学 | Non-linear model predication control method of tank reactor based on on-line support vector machine |
Non-Patent Citations (2)
Title |
---|
JIAN SHI 等: "Melt index prediction by neural networks based on independent component analysis and multi-scale analysis", 《NEUROCOMPUTING》 * |
夏陆岳 等: "基于SNNs-RR的聚丙烯熔融指数软测量", 《化工学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680428A (en) * | 2016-12-19 | 2017-05-17 | 华北电力大学(保定) | Soft measuring method for denitration control system |
CN108388113A (en) * | 2018-02-07 | 2018-08-10 | 浙江大学 | Least square method supporting vector machine soft-measuring modeling method based on distribution estimation local optimum |
CN108803525A (en) * | 2018-06-28 | 2018-11-13 | 浙江大学 | A kind of quick high-precision propylene polymerization production process optimal soft survey instrument of chaos |
CN108804851A (en) * | 2018-06-28 | 2018-11-13 | 浙江大学 | A kind of high-precision propylene polymerization production process optimal soft survey instrument of chaos gunz optimizing |
CN110378461A (en) * | 2019-06-14 | 2019-10-25 | 上海交通大学 | Error self-adaptation control method and system based on space-time error separate |
CN110750756A (en) * | 2019-10-01 | 2020-02-04 | 深圳市行健自动化股份有限公司 | Method for checking and diagnosing real-time online instrument by optimal support vector machine algorithm |
CN110750756B (en) * | 2019-10-01 | 2023-06-20 | 深圳市行健自动化股份有限公司 | Real-time on-line instrument checksum diagnosis method through optimal support vector machine algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN103675011B (en) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101315557B (en) | Propylene polymerization production process optimal soft survey instrument and method based on genetic algorithm optimization BP neural network | |
CN103674778B (en) | The industrial melt index soft measurement instrument of RBF particle group optimizing and method | |
CN103675011B (en) | The industrial melt index soft measurement instrument of optimum support vector machine and method | |
CN103675006B (en) | The industrial melt index soft measurement instrument of least square and method | |
CN101382801B (en) | Optimum soft measuring instrument based on EGA-optimized polymerization of propylene production process and method | |
CN110348075A (en) | A kind of grinding surface roughness prediction technique based on improvement algorithm of support vector machine | |
CN103472865B (en) | The pesticide waste liquid incinerator furnace temperature optimization system of intelligence least square and method | |
CN103675005B (en) | The industrial melt index soft measurement instrument of optimum FUZZY NETWORK and method | |
CN103823430A (en) | Intelligent weighing propylene polymerization production process optimal soft measurement system and method | |
CN103675012B (en) | The industrial melt index soft measurement instrument of BP particle group optimizing and method | |
CN103675010B (en) | The industrial melt index soft measurement instrument of support vector machine and method | |
CN103839103B (en) | Propylene polymerization production process BP Optimal predictor system and method | |
CN103675009B (en) | The industrial melt index soft measurement instrument of fuzzifying equation and method | |
CN103955170A (en) | Propylene polymerization production process online forecasting system and method based on group intelligent optimization | |
CN103838206A (en) | Optimal soft measurement meter and method in optimal BP multi-mode propylene polymerization production process | |
CN103838209B (en) | Propylene polymerization production process adaptive optimal forecast system and method | |
CN103675008B (en) | The industrial melt index soft measurement instrument of Weighted Fuzzy and method | |
CN103630568B (en) | The industrial melt index soft measurement instrument of BP network and method | |
CN103675007B (en) | The industrial melt index soft measurement instrument of RBF network and method | |
CN103838957A (en) | Propylene polymerization production process radial basis optimum soft measurement instrument and method | |
CN103838958A (en) | Vague intelligent optimal soft measuring instrument and method used in propylene polymerization production process | |
CN103838205A (en) | Optimum soft measurement instrument and method in BP global optimum propylene polymerization production process | |
CN103472867A (en) | Pesticide production waste liquid incinerator temperature optimization system and method of support vector machine | |
CN103823369A (en) | Propylene polymerization production process prediction system based on BP multimode network and method | |
CN103823963A (en) | Soft measuring system and method for polypropylene melt index of robusts |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150930 Termination date: 20180922 |
|
CF01 | Termination of patent right due to non-payment of annual fee |