CN112016243B - A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface - Google Patents

A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface Download PDF

Info

Publication number
CN112016243B
CN112016243B CN202010748756.XA CN202010748756A CN112016243B CN 112016243 B CN112016243 B CN 112016243B CN 202010748756 A CN202010748756 A CN 202010748756A CN 112016243 B CN112016243 B CN 112016243B
Authority
CN
China
Prior art keywords
parameter
beta
traffic flow
model
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010748756.XA
Other languages
Chinese (zh)
Other versions
CN112016243A (en
Inventor
赵德
王炜
李东亚
周伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202010748756.XA priority Critical patent/CN112016243B/en
Publication of CN112016243A publication Critical patent/CN112016243A/en
Application granted granted Critical
Publication of CN112016243B publication Critical patent/CN112016243B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/06Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Software Systems (AREA)
  • Evolutionary Computation (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Development Economics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Computer Hardware Design (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Educational Administration (AREA)
  • Chemical & Material Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Operations Research (AREA)
  • Data Mining & Analysis (AREA)
  • Computing Systems (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于响应曲面的交通流量预测模型参数快速标定方法,包括以下步骤:S1、确定交通流量预测模型参数的解空间;S2、计算初始化参数组合下的模型误差;S3、生成参数与模型误差间的响应曲面;S4、收敛性检验;S5、确定下一个参数组合点;S6、计算在该参数组合点下的模型误差。本发明方法可以实现交通流预测模型参数的快速标定,不需要对所有参数组合下的模型预测结果进行计算,大大节省模型试算占用的时间,并且本发明方法兼顾搜索效率与公平性,可防止标定结果限于局部最优解。同时,本发明方法应用场景多、操作简单、易于编程。

Figure 202010748756

The invention discloses a method for quickly calibrating parameters of a traffic flow prediction model based on a response surface, comprising the following steps: S1, determining the solution space of the parameters of the traffic flow prediction model; S2, calculating the model error under the combination of initialization parameters; S3, generating parameters The response surface between the model error and the model error; S4, the convergence test; S5, determine the next parameter combination point; S6, calculate the model error under the parameter combination point. The method of the invention can realize the rapid calibration of the parameters of the traffic flow prediction model, and does not need to calculate the model prediction results under all parameter combinations, which greatly saves the time occupied by the model trial calculation. The calibration results are limited to local optimal solutions. Meanwhile, the method of the present invention has many application scenarios, simple operation and easy programming.

Figure 202010748756

Description

Traffic flow prediction model parameter calibration method based on response surface
Technical Field
The invention relates to the fields of traffic transportation engineering and intelligent traffic, in particular to a traffic flow prediction model parameter calibration method based on a response curved surface.
Background
The traffic flow prediction is the basis of technologies such as traffic state deduction, navigation vehicle path selection, intelligent traffic control and the like of an intelligent traffic system, and an accurate traffic flow prediction model can better serve the operation of the intelligent traffic system. In recent years, with the advance of a new technological revolution such as traffic big data, artificial intelligence, deep learning and the like, a traffic flow prediction model is greatly developed, and a more complex and more accurate traffic flow prediction model, such as a traffic flow prediction model based on deep learning, is derived.
Although the model can greatly improve the prediction accuracy of the traffic flow, the model has the problem of difficult calibration of model parameters (namely, finding optimal parameters), and the model parameters and hyper-parameters are too many, the solution space is too large, and the solution is difficult to solve by direct methods such as enumeration and the like.
Disclosure of Invention
The purpose of the invention is as follows: the invention provides a traffic flow prediction model parameter calibration method based on a response curved surface, which can realize the rapid calibration of traffic flow prediction model parameters and obtain the optimal calibration result.
The technical scheme is as follows: the invention relates to a traffic flow prediction model parameter calibration method based on a response surface, which comprises the following steps:
s1, determining a solution space of the traffic flow prediction model parameters; traffic flow prediction model f (beta)12,…,βn) With n parameters beta12,…,βnThe value range of each parameter is the solution space of each parameter, and is respectively expressed as Z1,Z2,…,Zn
S2, calculating a model error under the initialization parameter combination;
s21, solution space Z of the parameters determined in step S11,Z2,…,ZnIn the method, m groups of values of n parameters are randomly extracted by adopting a Latin hypercube sampling method to form a parameter library U which is expressed as { beta [ ]1 12 1,…,βn 1},{β1 22 2,…,βn 2},…,{β1 m2 m,…,βn m}; wherein m is>0.5n2+1.5n+1;
S22, calculating a traffic flow prediction model f (beta) under the condition of m groups of parameter values12,…,βn) M sets of predicted results X1,X2,…,Xm
S23, calculating the error between the m groups of prediction results and the corresponding real results, wherein the error calculation formula is as follows:
εi=|Yi-Xi|,i∈{1,2,...,m}
wherein, XiFor the i-th group of predictors, YiAnd the real result is corresponding to the ith group of prediction results.
S3, generating a response surface between the parameters and the model errors; for the m sets of prediction errors of step S23, fitting the prediction errors epsilon and n parameters (beta) by using a multiple quadratic regression equation g12,…,βn) The functional relationship between ═ g (. beta.), (beta.)12,…,βn) The functional relation is the current response curved surface between the parameter and the model error; solving by using maximum likelihood estimation or a least square method during fitting;
s4, testing convergence; if any one of the following convergence conditions is met or two convergence conditions are met simultaneously, the convergence test is passed, and a final calibration result of the model parameters is output: in the solution space Z of the parameters1,Z2,…,ZnSolving the current response curved surface g (beta) by adopting a Newton method in the range12,…,βn) The corresponding parameter combination of the lowest point; otherwise, the convergence does not pass, and the process advances to step S5;
convergence condition 1:
Figure BDA0002609315740000021
wherein e is the maximum allowable error;
convergence condition 2: m > M, wherein M is the maximum number of tests;
s5, determining the next parameter combination point;
s51, solving space Z in parameter1,Z2,…,ZnIn the range, the current response curved surface g (beta) is solved by adopting a Newton method12,…,βn) Corresponding to the lowest point of { beta } in the parameter set1 k2 k,…,βn k};
S52, solving space Z in parameter1,Z2,…,ZnWithin the range, a set of parameter combinations { beta } is randomly chosen1 l2 l,…,βn lIf the parameter combination { beta }1 l2 l,…,βn lRe-extracting if the element belongs to U until the extracted parameter combination
Figure BDA0002609315740000022
Figure BDA0002609315740000023
S53, randomly selecting one parameter from the two parameter combinations generated in the steps S51 and S52 according to 1/2 probability as the next parameter combination point;
s54, adding the next parameter combination point determined in the step S53 into a parameter base U, and updating m to m + 1;
s6, calculating the model error under the parameter combination point; the next parameter combination point determined in step S53 is calculated in step S22, and the traffic flow prediction model f (β)12,…,βn) The prediction error at this point is calculated by the method of step S23, and the process returns to step S3.
Has the advantages that:
1) the response surface of the invention can approximate the mapping relation between the multivariate parameter and the prediction error, and does not need to calculate the model prediction results under all parameter combinations, thereby greatly saving the time occupied by model trial calculation, being particularly effective when the model parameter combinations are more, and realizing the rapid calibration of the traffic flow prediction model parameters.
2) Compared with a conventional heuristic algorithm, the method disclosed by the invention has the advantages of both search efficiency and fairness. When the search direction of the solution is determined, on one hand, efficient search is considered, and the optimal solution under the current response curved surface is searched; on the other hand, other space ranges which are not collected are continuously explored through convergence test, the obtained calibration result is prevented from being limited to a local optimal solution, and the optimal calibration result is obtained in a plurality of groups of parameter combinations which approach to the real result.
3) The method can cope with more diversified scenes. The method of the invention can be applied to continuous parameters, discrete parameters, integer parameters and combinations thereof; further, the traffic flow prediction model may be various, such as machine learning, deep learning, time series, and the like.
4) The method is simple to operate and easy to program.
Drawings
FIG. 1 is a flow chart of the present invention;
Detailed Description
As shown in fig. 1, the traffic flow prediction model parameter calibration method based on the response surface provided by the invention includes the following steps: s1, determining a solution space of the traffic flow prediction model parameters; s2, calculating a model error under the initialization parameter combination; s3, generating a response surface between the parameters and the model errors; s4, testing convergence; s5, determining the next parameter combination point; s6, calculating the model error at the parameter combination point.
S1: determining a solution space for traffic flow prediction model parameters
Traffic flow prediction model f (beta)12,…,βn) With n parameters beta12,…,βnThe value range of each parameter is the solution space of each parameter, and is respectively expressed as Z1,Z2,…,Zn
S2: calculating model error under initialization parameter combination
The solution space Z of the parameters determined at step S11,Z2,…,ZnIn the method, m groups of values of n parameters are randomly extracted by adopting a Latin hypercube sampling method to form a parameter library U which is expressed as { beta [ ]1 12 1,…,βn 1},{β1 22 2,…,βn 2},…,{β1 m2 m,…,βn m}; wherein m is>0.5n2+1.5n+1;
Calculating a traffic flow prediction model f (beta) under the condition of m groups of parameter values12,…,βn) M sets of predicted results X1,X2,…,Xm
And calculating the error between the m groups of predicted results and the corresponding real results, wherein the error calculation formula is as follows:
εi=|Yi-Xi|,i∈{1,2,...,m}
wherein, XiIs the ith group of prediction results; y isiThe real result corresponding to the prediction result of the ith group can be obtained by field observation and collection under corresponding conditions.
S3: generating a response surface between parameters and model errors
For the m sets of prediction errors of step S2, fitting the prediction errors epsilon and n parameters (beta) by using a multiple quadratic regression equation g12,…,βn) The functional relationship between ═ g (. beta.), (beta.)12,…,βn) The functional relation is the current response curved surface between the parameter and the model error; solving by using maximum likelihood estimation or a least square method during fitting;
s4: convergence test
If any one of the following convergence conditions is met or two convergence conditions are met simultaneously, the convergence test is passed, and a final calibration result of the model parameters is output: in the solution space Z of the parameters1,Z2,…,ZnSolving the current response curved surface g (beta) by adopting a Newton method in the range12,…,βn) The corresponding parameter combination of the lowest point; otherwise, the convergence does not pass, and the process advances to step S5;
convergence condition 1:
Figure BDA0002609315740000041
wherein e is the maximum allowable error;
convergence condition 2: m > M, wherein M is the maximum number of tests;
s5: determining a next parameter combination point
In the solution space Z of the parameters1,Z2,…,ZnIn the range, the current response curved surface g (beta) is solved by adopting a Newton method12,…,βn) Corresponding to the lowest point of { beta } in the parameter set1 k2 k,…,βn k};
In the solution space Z of the parameters1,Z2,…,ZnWithin the range, a set of parameter combinations { beta } is randomly chosen1 l2 l,…,βn lIf the parameter combination { beta }1 l2 l,…,βn lRe-extracting if the element belongs to U until the extracted parameter combination
Figure BDA0002609315740000042
Figure BDA0002609315740000043
Randomly selecting one parameter from the two sets of parameter combinations generated in the previous step of S5 as a next parameter combination point according to 1/2 probability;
adding the parameter combination point to a parameter library U, and updating m to m + 1;
s6: calculating the model error at the parameter combination point
The next parameter combination point determined in step S5 is calculated in step S22, and the traffic flow prediction model f (β)12,…,βn) The prediction error at this point is calculated by the method of step S2, and the process returns to step S3.

Claims (6)

1. A traffic flow prediction model parameter calibration method based on a response surface is characterized in that: the method comprises the following steps:
s1, determining a solution space of the traffic flow prediction model parameters; traffic flow prediction model f (beta)12,…,βn) With n parameters beta12,…,βnThe value range of each parameter is the solution space of each parameter, and is respectively expressed as Z1,Z2,…,Zn
S2, calculating a model error under the initialization parameter combination;
s21, parameter solution space Z determined in step S11,Z2,…,ZnIn the method, m groups of values of n parameters are randomly extracted by adopting a Latin hypercube sampling method to form a parameter library U which is expressed as { beta [ ]1 12 1,…,βn 1},{β1 22 2,…,βn 2},…,{β1 m2 m,…,βn m};
S22, calculating a traffic flow prediction model f (beta) under the condition of m groups of parameter values12,…,βn) M sets of predicted results X1,X2,…,Xm
S23, calculating the error between the m groups of prediction results and the corresponding real results;
s3, generating a response surface between the parameters and the model errors: for the m sets of prediction errors of step S2, fitting the prediction errors epsilon and n parameters (beta) by using a multiple quadratic regression equation g12,…,βn) The functional relationship between ═ g (. beta.), (beta.)12,…,βn) The functional relation is the current response curved surface between the parameter and the model error; solving by using maximum likelihood estimation or a least square method during fitting;
s4, testing convergence; setting a convergence condition, if the convergence condition is met, checking to pass, and outputting a final calibration result of the model parameters; otherwise, the convergence does not pass, and the process advances to step S5;
s5, adding the next parameter combination point;
s51, solving space Z in parameter1,Z2,…,ZnIn the range, the current response curved surface g (beta) is solved by adopting a Newton method12,…,βn) Corresponding to the lowest point of { beta } in the parameter set1 k2 k,…,βn k};
S52, solving space Z in parameter1,Z2,…,ZnWithin the range, a set of parameter combinations { beta } is randomly chosen1 l2 l,…,βn lIf the parameter combination { beta }1 l2 l,…,βn lRe-extracting if the element belongs to U until the extracted parameter combination
Figure FDA0003390866330000011
Figure FDA0003390866330000012
S53, randomly selecting one parameter from the two parameter combinations generated in the steps S51 and S52 according to 1/2 probability as the next parameter combination point;
s54, adding the next parameter combination point determined in the step S53 into a parameter base U, and updating m to m + 1;
s6, calculating model errors under the newly added p parameter combination points; step S22 is adopted for p parameter combination points newly added in step S5 to solve the traffic flow prediction model f (beta)12,…,βn) The prediction error at this point is calculated by the method of step S23, and the process returns to step S3.
2. The traffic flow prediction model parameter calibration method based on the response surface according to claim 1, characterized in that: the convergence condition in step S4 includes two conditions:
convergence condition 1:
Figure FDA0003390866330000021
wherein e is the maximum allowable error;
convergence condition 2: m > M, wherein M is the maximum number of tests,
if any convergence condition is met or two convergence conditions are met simultaneously, the convergence check is passed, and a final calibration result of the model parameters is output.
3. The traffic flow prediction model parameter calibration method based on the response surface according to claim 1 or 2, characterized in that: the value of m in the step S21 has the following requirements:
m>0.5n2+1.5n+1。
4. the traffic flow prediction model parameter calibration method based on the response surface as claimed in claim 3, wherein: the error calculation formula in step S23 is as follows:
εi=|Yi-Xi|,i∈{1,2,...,m};
wherein, XiFor the i-th group of predictors, YiAnd the real result is corresponding to the ith group of prediction results.
5. The traffic flow prediction model parameter calibration method based on the response surface as claimed in claim 4, wherein: in step S3, ∈ ═ g (β)12,…,βn) Maximum likelihood estimation or least square method is adopted in fitting.
6. The traffic flow prediction model parameter calibration method based on the response surface as claimed in claim 5, wherein: the final calibration result of the model parameters in step S4 is: in the solution space Z of the parameters1,Z2,…,ZnSolving current response surface g (beta) in range12,…,βn) Corresponding to the lowest point of (a).
CN202010748756.XA 2020-07-30 2020-07-30 A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface Active CN112016243B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010748756.XA CN112016243B (en) 2020-07-30 2020-07-30 A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010748756.XA CN112016243B (en) 2020-07-30 2020-07-30 A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface

Publications (2)

Publication Number Publication Date
CN112016243A CN112016243A (en) 2020-12-01
CN112016243B true CN112016243B (en) 2022-02-15

Family

ID=73498894

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010748756.XA Active CN112016243B (en) 2020-07-30 2020-07-30 A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface

Country Status (1)

Country Link
CN (1) CN112016243B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111192453A (en) * 2019-12-30 2020-05-22 深圳市麦谷科技有限公司 Short-term traffic flow prediction method and system based on Bayesian optimization

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007061618A2 (en) * 2005-11-22 2007-05-31 Exxonmobil Upstream Research Company Simulation system and method
CN100561509C (en) * 2006-07-11 2009-11-18 南京大学 An Improved Hybrid Genetic Algorithm Design Method for Optimizing Water Quality Model Parameters
CN108108583A (en) * 2016-11-24 2018-06-01 南京理工大学 A kind of adaptive SVM approximate models parameter optimization method
CN106844924B (en) * 2017-01-12 2021-02-23 桂林电子科技大学 Method for optimizing PCB microstrip line structure based on response surface method and genetic algorithm

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111192453A (en) * 2019-12-30 2020-05-22 深圳市麦谷科技有限公司 Short-term traffic flow prediction method and system based on Bayesian optimization

Also Published As

Publication number Publication date
CN112016243A (en) 2020-12-01

Similar Documents

Publication Publication Date Title
CN111199016A (en) DTW-based improved K-means daily load curve clustering method
CN112115419B (en) System state estimation method and system state estimation device
CN110705029B (en) A flow field prediction method for oscillating flapping-wing energy harvesting system based on transfer learning
CN116448419A (en) Zero-sample bearing fault diagnosis method based on high-dimensional parameter multi-objective efficient optimization of deep model
CN108334898A (en) A kind of multi-modal industrial process modal identification and Fault Classification
CN110059337A (en) Gyro system life-span prediction method based on Bayes multisource data fusion
CN112016253A (en) High-fidelity chaotic polynomial correction method suitable for CFD uncertainty quantification
Wang et al. Fault feature selection based on modified binary PSO with mutation and its application in chemical process fault diagnosis
Bontempi et al. Recursive lazy learning for modeling and control
CN115391745B (en) Rainfall forecast correction method and system based on probability matching average method
CN104657442B (en) Multiple target community discovery method based on Local Search
CN114564787B (en) Bayesian optimization method, device and storage medium for target related airfoil design
Zhang et al. A right-hand side function surrogate model-based method for the black-box dynamic optimization problem
CN112016243B (en) A Parameter Calibration Method for Traffic Flow Prediction Model Based on Response Surface
US12223012B2 (en) Machine learning variable selection and root cause discovery by cumulative prediction
Yaro et al. Enhancing fingerprint localization accuracy with inverse weight-normalized context similarity coefficient-based fingerprint similarity metric
CN106203737A (en) The multistep forecasting method of a kind of wind farm wind velocity and system
JP2011008562A (en) Device and method for generating plant control information, and computer program therefor
CN109711040B (en) Intelligent industrial design reinforcement learning algorithm based on search direction learning
CN110909492A (en) Sewage treatment process soft measurement method based on extreme gradient lifting algorithm
Dan et al. Application of machine learning in forecasting energy usage of building design
CN105824994A (en) Zoning-based combined approximate model building method
CN110147614B (en) Engineering safety evaluation method based on grading difference Stacking multi-model ensemble learning
Lv et al. A dynamic programming bayesian network structure learning algorithm based on mutual information
Gao et al. Parameter estimation of Gaussian mixture model and its application in multimode process monitoring

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant