63 results about "Linear problem" patented technology

The invention relates to the technical field of expressway monitoring, and provides an expressway abnormal state detection method and system based on a support vector machine, which processes vehicle trajectory data and ETC gantry data through Kalman filtering and adaptive weighted average to improve the data accuracy, introduces a Gaussian kernel function to construct an SVM model, and improves the detection accuracy of the expressway abnormal state. An original feature space is mapped to a high-dimensional space to deal with a nonlinear problem, samples with labels are used to train a multi-classification SVM model, and accurate detection of various abnormal states such as traffic congestion, accidents and vehicle faults is realized through a voting mechanism. According to the method, the accuracy of abnormal state detection is improved, powerful decision support is provided for expressway managers, the traffic flow is optimized, the road traffic capacity is improved, and the development of an intelligent traffic system is promoted.

The invention discloses a circulating fluidized bed temperature intelligent prediction method based on a physical information neural network, and mainly relates to the technical field of industrial process intelligent control and energy power engineering crossing. Comprising the following steps: S1, acquiring operation data of the circulating fluidized bed, and preprocessing the operation data to obtain a data set; s2, constructing a physical information residual neural network; s3, constructing a loss function of adaptive weight adjustment; s4, training the physical information residual neural network by using the data set, and performing iterative optimization in combination with the constructed physical information residual neural network and the loss function to obtain a bed temperature prediction model; s5, inputting operation data of the circulating fluidized bed into the bed temperature prediction model to obtain a bed temperature prediction value; the method can solve the problem of complexity and nonlinearity of dynamic prediction of the bed temperature of the circulating fluidized bed, can optimize the operation efficiency of the boiler, improves the combustion stability, and reduces the emission of pollutants.

According to the plateau area load prediction method and system based on the signal decomposition and the neural network, an original signal is decomposed into a plurality of intrinsic mode components with physical significance by adopting completely adaptive noise set empirical mode decomposition, mode aliasing is suppressed, and multi-scale features are extracted to effectively suppress load demand oscillation; the problem of non-stationarity of the load of the service area is solved; then, K-means clustering grouping is carried out on the components, similar modes are combined to reduce calculation redundancy, and the problems of heterogeneity and local feature redundancy in complex time series data are effectively solved; performing nonlinear dimension reduction processing on the grouped reconstructed signals by adopting kernel principal component analysis, eliminating noise interference and extracting key features, thereby effectively solving the multiple nonlinear problem; and finally, inputting the dimensionality-reduced features into a long-short-term memory network modeling time sequence dependency relationship, and realizing high-precision prediction through adaptive weighted integration of component prediction results.

The invention provides a structure reliability analysis method based on an adaptive variable fidelity model, and the method comprises the steps: generating an initial sample point set which is uniformly distributed and has representativeness through an improved random sampling method KMODMC, enabling the initial sample point set to comprise a low-fidelity sample set and a high-fidelity sample set, training a BP neural network through employing the low-fidelity sample set, and carrying out the training of the BP neural network through employing the high-fidelity sample set; a low-fidelity BP neural network model is obtained; meanwhile, based on an error training Kriging model of a high-fidelity sample set and a low-fidelity model predicted value, an error correction Kriging model is constructed, the low-fidelity BP neural network model and the error correction Kriging model are combined to form a multi-fidelity mixed agent model, and adaptive iterative optimization is performed through a double-model alternate point adding sampling strategy. And finally obtaining a high-precision multi-fidelity hybrid agent model for structural reliability evaluation. According to the method, the problems of high cost of high-fidelity simulation calculation and insufficient precision of a low-fidelity model are solved, and the adaptive capacity and prediction precision of the model in a complex nonlinear problem are effectively improved.

The invention discloses a formation resistivity determination method and device for a post-casing transient electromagnetic logging curve, equipment and a medium, and relates to the technical field of oil exploration and development, and the formation resistivity determination method comprises the following steps: obtaining a post-casing electromagnetic logging curve, and carrying out time sequence conversion and data preprocessing to obtain time sequence data; converting the data into frequency domain data by adopting wavelet transform, and performing multi-scale refinement splitting to obtain low-frequency features and high-frequency features; fitting the low-frequency features and the high-frequency features by using a symbol regression method to obtain a nonlinear high-dimensional symbol regression formula for the low-frequency data and a nonlinear high-dimensional symbol regression formula for the high-frequency data, and generating a post-casing electromagnetic logging curve nonlinear fitting formula; the mapping relation between the post-casing transient electromagnetic curve and the formation resistivity is determined through the post-casing electromagnetic logging curve nonlinear fitting formula, the formation resistivity is obtained, the uncertainty of manual selection of model parameters is avoided, the capacity of processing nonlinear problems is improved, and the accuracy and reliability of high-dimensional relation mapping and nonlinear fitting are improved.

The invention relates to the technical field of electric power, in particular to an electricity price fluctuation quantification method, device, equipment and medium based on regression analysis, and the method comprises the steps: taking electric power price fluctuation in peak load of electric power as a dependent variable, and taking new energy power generation, thermal power generation, electric power demand, generating capacity and capacity electricity price as independent variables; establishing a multiple linear regression model; verifying whether the independent variable and the dependent variable have a linear relationship or not, preliminarily carrying out multiple linear regression analysis, finding an optimal regression coefficient by minimizing a residual sum of squares, and carrying out multiple linear problem judgment; a ridge regression analysis method is adopted, a regularization item is added into a loss function to limit the size of a regression coefficient, overfitting of the model is reduced, and a ridge function is confirmed through cross validation; performing overall significance test and regression coefficient significance test on the multiple linear regression model, and analyzing the fitting degree; and performing quantitative analysis on electricity price fluctuation by using the fitted multiple linear regression model.

The invention provides a dynamic state estimation method for minimum information loss of a heat network and a gas network of a multi-energy system, and the method comprises the steps: building all subsystems of a comprehensive energy system and a multi-energy system coupling equipment model based on different physical characteristics and different time scales of an electric system, a heat system and a gas system; establishing a heat supply network state estimation model and a pseudo-dynamic state estimation model based on each subsystem and the multi-energy system coupling equipment model; substituting the parameters of the pre-estimated moment into a heat supply network state estimation model to obtain a heat supply network state estimation value; and substituting the parameters of the pre-estimated moment into the quasi-dynamic state estimation model, and solving the quasi-dynamic state estimation model by adopting an alternating direction multiplier method to obtain a gas network state estimation value. According to the method, a nonlinear problem in state estimation is converted into a linear problem, the accuracy of state estimation is greatly improved, the contradiction between the calculation complexity and the solution accuracy in the dynamic calculation of the multi-energy system network is well balanced, the numerical stability is good, and the calculation efficiency is high.

The invention provides an all-electric energy station operation optimization method and system considering multiple objectives. The method comprises the following steps: acquiring equipment parameters and operation parameters of an all-electric energy station; based on the equipment parameters and the operation parameters, a comprehensive multi-objective hierarchical optimization model of the photovoltaic-hydrogen-natural gas comprehensive energy system model is solved by adopting a genetic-simulated annealing algorithm, and an optimal operation strategy of the full-electric energy station is obtained; according to the method, when the genetic-simulated annealing algorithm is adopted to solve the comprehensive multi-objective hierarchical optimization model, global search can be realized through group search and natural evolution operation, the capability of solving the complex nonlinear problem in the solving process can be improved, multi-objective optimization can still be implemented when the complex nonlinear problem exists, and the method is suitable for large-scale popularization and application. Therefore, the optimal operation strategy obtained by solving the model genetic-simulated annealing algorithm can significantly improve the flexible adjustment capability and the environmental change adaptation capability of the full-electric energy station.

The invention provides a topological optimization structure design method based on space control points and a space point control technology of structure topological optimization. The SPCT constructs a field function by utilizing spatial arrangement of control points, so that the number of design variables of topological optimization is effectively reduced, the optimization calculation efficiency is improved, and a chessboard effect and medium-density elements in a traditional topological optimization method are avoided. The control point placement mode is visual, the initial configuration or sample data can be conveniently controlled, and the optimization efficiency is further improved. Besides, sensitivity information based on an SPCT topological optimization model is deduced, a gradient-based optimization framework is established, and a non-gradient algorithm solution scheme can be used to solve the complex nonlinear problem without the help of gradient information. The technology has obvious advantages in the aspects of reducing design variables and improving calculation efficiency, and particularly has wide application prospects in the field of structure optimization design.

The invention relates to the technical field of loss reduction optimization of a power distribution network, in particular to a power grid loss reduction method and a power grid loss reduction system considering Lorenna series power flow approximation under data driving, and the method comprises the following steps: converting a line loss optimization task of the power distribution network into an intelligent agent interaction environment; the intelligent agent interaction environment is used for converting a line loss optimization task of the power distribution network into a Markov decision process under a data driving framework; and by adopting the load flow approximation method based on the Loran series, a nonlinear problem in the traditional load flow calculation is converted into a linear problem, so that the calculation complexity is greatly reduced. According to the method, the load flow calculation speed can be remarkably increased on the premise that the calculation precision is guaranteed, and power grid loss reduction optimization can meet the real-time requirement. Especially in a high-permeability distributed new energy access complex power grid, the method can quickly respond to the dynamic change of the power grid state, and effective support is provided for loss reduction optimization in actual operation.

According to the source-load-storage coupled multi-energy micro-grid double-layer optimization method and system for the oil and gas well field, an upper-layer model, namely a load optimization layer model with the purpose of maximizing the green power consumption rate, is constructed respectively, the working system and production parameters of production equipment are optimized, and the wind and light consumption level is enhanced while the production requirement is guaranteed; a lower-layer model, namely a power generation dispatching layer model taking the lowest daily operation cost of the multi-energy micro-grid as a target, is constructed, the energy storage peak regulation capability is fully exerted, and the economic benefit of system operation is improved; meanwhile, due to the fact that the coupling relation exists between the upper layer and the lower layer of the double-layer model, the double-layer nonlinear problem is converted into single-layer mixed integer linear programming through the KKT condition and Big-M linearization, the solving time is greatly shortened, and the method is suitable for a real-time scheduling scene. And the green electricity consumption rate and the well group daily fluid production capacity are respectively improved by 9.34% and 9.61%. And the double-target balance of economy and environmental protection is realized.

The invention discloses a semi-analytic solution-based energy flow calculation method and system for an electricity-gas integrated energy system, and the method comprises the steps: firstly building an ordinary differential-algebraic model of an electricity-gas network based on a numerical boundary extrapolation and semi-discrete method; based on a semi-analytical theory of differential transformation, establishing a linearized electricity-gas network model; carrying out block calculation on a matrix and self-adaptive adjustment on a time window; and finally solving the model to obtain an energy flow calculation result. According to the method, the original nonlinear partial differential equation of the gas network is converted into the linear algebraic equation, the complex nonlinear problem in traditional gas network modeling is simplified, errors caused by local linearization and time difference are avoided, the calculation scale is reduced, and the calculation robustness is improved.

The invention discloses a damaged road network recovery sequence optimization method with maximum toughness as a target based on user balance flow distribution, and aims to solve the problems of low recovery efficiency of a post-disaster traffic road network and lack of scientific basis for decision making. The method comprises the following implementation steps: firstly, constructing a damaged road network toughness evaluation framework, and defining a toughness index as a reciprocal of a product of recovery time and total transit time in combination with a traffic network topological structure and a flow demand; secondly, establishing damaged road network toughness optimization scene parameters, abstracting a road network into a connected graph, and defining a recovery scheme set; then, a mixed integer linear programming model with the maximum toughness as the target is constructed, a nonlinear problem is converted into a linear problem through secant approximation, and constraint conditions such as flow balance and travel requirements are set; and finally, performing global optimization solution on the model by using a high-performance solver to generate an optimal recovery scheme. According to the method, the road network toughness can be scientifically quantified, the recovery sequence is optimized, blindness of experience decision making is avoided, the passing efficiency of a post-disaster traffic system is rapidly improved, efficient and scientific decision support is provided for post-disaster traffic recovery, and the method has important practical application value.

The invention relates to a low-carbon economic dispatching method and system for an integrated energy system, and the method comprises the steps: accurately describing the variable efficiency characteristics of energy conversion equipment through constructing a dynamic energy concentrator model based on a BP neural network, and processing a nonlinear problem through employing a linearization method; a typical source load prediction scene is generated in combination with a scene generation and reduction technology, and an energy flow-carbon flow double-layer optimization model is established. Wherein the energy flow layer takes the minimum operation cost as a target and considers equipment operation, multi-energy flow balance and transaction constraints; a stepped carbon transaction mechanism and demand response guided by carbon potential are introduced into a carbon flow layer, and the purpose is to minimize the cost of carbon transaction and demand response. And through iterative solution of the double-layer model, collaborative optimization of a system energy purchase and sale plan, equipment output and load change is realized until an optimal scheduling scheme is converged. According to the method, collaborative optimization of the energy flow and the carbon flow is realized, and double improvement of the scheduling of the comprehensive energy system in economy and low carbon is promoted.

The invention discloses an electricity-gas integrated energy system scheduling method based on a global linearization gas network model, and the method comprises the steps: generating a gas network pipeline operation data set, and building a global linearization gas network model based on a Koopman operator theory; then, a Koopman operator is approximated, and an integrated energy system scheduling model is constructed; and finally solving the model to obtain an optimal scheduling result. According to the method, an original nonlinear partial differential equation of the gas network can be converted into a linear algebraic equation, the complex nonlinear problem in traditional gas network modeling is simplified, errors caused by local linearization and spatial difference are avoided, the linear algebraic equation can be easily integrated into a scheduling model of a comprehensive energy system, and the method is suitable for large-scale popularization and application. And the dynamic behavior of the natural gas system can be accurately captured, and the scheduling accuracy and reliability are remarkably improved. In addition, according to the method, errors possibly generated by a traditional model in a complex gas dynamics scene are effectively reduced, and system potential safety hazards caused by inaccurate models are avoided.

This invention discloses a real-time attitude estimation system for non-cooperative industrial targets based on strong tracking filtering, relating to the fields of industrial machine vision and automation control technology. This system decouples the complete six-degree-of-freedom attitude estimation problem into a linear translational sub-state space and a nonlinear rotational sub-state space. This decomposes the nonlinear calculations that would normally be performed in six dimensions into a three-dimensional linear problem and a three-dimensional nonlinear problem. The translational sub-state only performs a simplified Kalman gain iteration with closed-form solutions, and matrix operations are simplified to low-order linear recursion. Although the rotational sub-state uses strong tracking filtering, it only needs to perform nonlinear updates in the reduced low-dimensional space. The Jacobian matrix sparsification pre-computation and incremental update unit pre-stores the constant part of the observation equation's derivative calculation offline and only incrementally updates the sparsely changing part online. This allows the system to be fully deployed on ARM architecture edge devices and stably meet millisecond-level processing requirements.

The application discloses a kind of space target ultra-short arc initial orbit intelligent determination method, belong to aerospace technology field.The application is based on artificial intelligence technology to realize space target ultra-short arc initial orbit determination, there is no strict requirement to complex nonlinear problem, and can develop more complex perturbation condition and initial orbit determination under the operation law of multiple body, coverage is wider;The application generates multiple groups of target orbits randomly and calculates the observation and state quantity of each visible arc segment, constructs the sample data set with ultra-large capacity, which is helpful for neural network to learn the nonlinear relationship between observation and state quantity using rich samples;The application builds and trains long short-term memory network model, determines the optimal network model weight parameter by designing and adjusting hyperparameter, can obtain the fast prediction solution of unknown space target ultra-short arc initial orbit according to observation;The application can meet more space target ultra-short arc initial orbit determination tasks by considering different inputs and outputs.

The present invention provides an active suspension control method, which includes the following steps: constructing a state space equation of the model according to the established dynamic model of the multi-variable non-linear active suspension to obtain an active suspension system model; linearizing the multi-variable non-linear links in the active suspension system model to obtain a corresponding quasi-linear system model of the active suspension; optimizing the quasi-linear system model of the active suspension by the Lagrange multiplier method, and solving by the bisection method to obtain the feedback gain K of the multi-variable non-linear active suspension and the feedback control law of the multi-variable non-linear active suspension. By adopting the above method, not only the non-linear problem of the inherent components is solved, but also the network non-linear problem is considered, improving the control accuracy and the stability of the system. At the same time, because the method is simple to reproduce, it is easy to implement in actual work.

The invention discloses an electrocatalyst activity attenuation prediction method based on microdynamics, and relates to the field of electrochemical catalyst.According to the electrocatalyst activity attenuation prediction method based on microdynamics, an electrocatalytic reaction and an inactivation reaction are dynamically coupled in a microdynamics model for the first time, the limitation of traditional thermodynamic stability analysis is broken through, and the activity attenuation process under the working condition can be truly reflected. By introducing a residual active site coverage variable and adopting quasi-steady-state approximation, the modeling and solving problems of a dynamic process of continuous loss of active sites are successfully solved, and a complex rigid differential equation set is simplified into an analyzable linear problem. The method not only can predict a stability attenuation curve, but also can quantitatively diagnose key steps influencing the service life through half-life period control degree analysis, realizes crossing from'phenomenon prediction 'to'mechanism interpretation' and'reverse design ', and provides a powerful theoretical tool for rationally designing a high-stability electrocatalyst.

The invention discloses a quantum linear solution method, device and equipment based on a polynomial preprocessor. The method comprises: obtaining information of a first matrix A and a first vector b in a linear system to be processed, calculating a polynomial p(A) of the first matrix A used for preprocessing the linear system, preprocessing the linear system according to the polynomial p(A), obtaining a second matrix A′ and a second vector b′, constructing a quantum circuit corresponding to an HHL algorithm, performing evolution and measurement operations of a quantum state according to the second matrix A′ and the second vector b′, and obtaining the final quantum state of the evolved quantum circuit, so as to solve the deficiencies in the prior art. The method can reduce the time complexity and the amount of calculation for solving linear problems, and at the same time reduce the occupation of hardware resources.

The invention provides an original real-time electricity price prediction method based on a quantum complex neural network and Hilbert-Huang Transform (HHT), and the real-time electricity price prediction method based on the quantum complex neural network and the Hilbert-Huang Transform (HHT). Aiming at the non-stationarity of electricity price data on the historical level, firstly, the time internal correlation of each feature channel is extracted through an HHT time sequence analysis method, and complex multi-dimensional time sequence prediction is simplified into a simple regression task, so that a tedious time sequence modeling process is avoided; for the non-linear problem of electricity price data, a quantum neural network is used for capturing the coupling relation between different factors, and the calculation speed and the model efficiency are further improved by means of the parallelism of quantum calculation. Based on the two advantages, the method can realize accurate real-time electricity price prediction.

The application provides a topology optimization structure design method based on a space control point, and a space point control technology for structure topology optimization. The SPCT constructs a field function by using the space arrangement of the control points, effectively reduces the number of design variables of the topology optimization, improves the optimization calculation efficiency, and avoids the appearance of the chessboard effect and the medium density elements in the traditional topology optimization method. The control point placement method is intuitive and convenient for controlling the initial configuration or sample data, and further improves the optimization efficiency. In addition, the application derives the sensitivity information of the SPCT topology optimization model, establishes a gradient-based optimization framework, and can use a non-gradient algorithm to solve the scheme without the aid of gradient information to solve complex nonlinear problems. The technology has obvious advantages in reducing design variables and improving calculation efficiency, and has broad application prospects in the field of structure optimization design.

A method and apparatus for simulating the operation of a cascade reversible hydropower station based on piecewise McCormick relaxation is disclosed. The method includes: constructing a simulation model of the reversible hydropower station operation based on the operating conditions of the station and aiming to maximize the total output of the station; linearizing the nonlinear part of the simulation model using piecewise McCormick relaxation and piecewise linearization to simplify it into a mixed-integer linear programming model; and iteratively relaxing the simulation model based on upstream runoff data and operating parameters to obtain the simulation results. This invention uses piecewise McCormick relaxation to transform the nonlinear problem into a linear problem, while employing an iterative relaxation algorithm to improve the solution accuracy. It does not require narrowing the feasible boundary of McCormick relaxation, has fewer iterations, and a faster convergence speed, thus improving the speed and stability of the solution while ensuring simulation accuracy.

The application provides a structure reliability analysis method based on an adaptive variable fidelity model, and the method comprises the following steps: generating an initial sample point set with uniform distribution and representativeness through an improved random sampling method KMODMC, the initial sample point set comprising a low-fidelity sample set and a high-fidelity sample set; subsequently, training a BP neural network by using the low-fidelity sample set to obtain a low-fidelity BP neural network model; at the same time, training a Kriging model based on the error of the high-fidelity sample set and the predicted value of the low-fidelity model to construct an error correction Kriging model; combining the low-fidelity BP neural network model and the error correction Kriging model to form a multi-fidelity hybrid proxy model; and performing adaptive iterative optimization through a double-model alternating point sampling strategy to finally obtain a high-precision multi-fidelity hybrid proxy model for structure reliability evaluation. The application solves the problems of high cost of high-fidelity simulation calculation and insufficient precision of a low-fidelity model, and effectively improves the adaptability and prediction precision of the model in complex nonlinear problems.

The invention discloses a structure nonlinear dynamic response incremental modal overlay analysis method and system, and belongs to the field of aerospace structure dynamics calculation. The method comprises the steps of constructing a finite element model and obtaining an initial matrix; dispersing the time-varying temperature field into a plurality of time increment intervals, updating material parameters in each interval, and calculating an instantaneous matrix; carrying out eigenvalue analysis and converting the kinetic equation into a modal coordinate system; solving a modal response increment by adopting a numerical integration method and converting the modal response increment to a physical coordinate system for superposition; and circularly executing until a full-time-history response is output. According to the method, the problem of rigidity nonlinearity caused by the temperature is solved through incremental linearization and modal order reduction, high calculation precision and high efficiency are achieved, and the method is suitable for dynamic design and analysis of aerospace structures in the high-temperature variable-temperature environment.

The invention belongs to the technical field of electric vehicle load optimization, and particularly discloses an electric vehicle load multi-target double-layer optimization method considering peak load shifting, comprising the following steps: establishing an electric vehicle load multi-target double-layer optimization model; converting a double-layer optimization problem of the multi-target double-layer optimization model into a single-layer multi-target linear problem through a KKT condition and a maximum M method; further converting a single-layer multi-target linear problem into a single-layer single-target problem; solving a single-layer single-target problem to obtain a Parteo leading edge; and determining an ideal solution at a Parteo leading edge, and taking the ideal solution as a final result of an original double-layer multi-objective optimization problem to realize electric vehicle load multi-objective double-layer optimization. According to the method, the problem that the upper-layer model in the double-layer model is single-target optimization due to the fact that comprehensive indexes such as power distribution network side peak-valley difference and load fluctuation are not considered in an existing double-layer optimization operation model is solved, and the peak-valley difference and the load fluctuation can be reduced while the income of an electric vehicle user is guaranteed.

The present invention discloses a distribution network reconstruction method and device based on operation uncertainty. An interval uncertainty set is constructed according to the uncertainty of distribution network operation. A linearized power flow constraint is established for the distribution network, and a maximum slack variable is used to replace the nonlinear power flow over-limit of the distribution network, so that the original nonlinear problem can be converted into a linear problem. Minimization of the over-limit frequency and degree of distribution network state variables is used as an objective function, and the objective function is solved in combination with the interval uncertainty set and the linearized power flow constraint to reconstruct the distribution network with the minimum risk. Therefore, the influence of the uncertainty of distribution network operation on the network reconstruction model can be considered, which is conducive to the rapid network reconstruction of the distribution network and is suitable for application in the scenario of distribution network network reconstruction.

This invention discloses a method, apparatus, and medium for solving variable quantum linear problems in subspaces. The method includes: determining the system of linear equations to be solved, constructing a variable quantum circuit and a Krylov subspace, and using the generalized minimum residual method and the variable quantum circuit to calculate an approximate solution of the system of linear equations to be solved in the subspace. This addresses the shortcomings of existing technologies and can reduce the time complexity and computational load of solving linear problems.

The invention discloses an underwater multi-target tracking method based on kernel mean value embedding belief propagation, which adopts a kernel mean value embedding method, and overcomes the problem of incomplete observability of a system with lower calculation complexity by approaching a probability density function of a target state through a small number of kernel particles. A strong nonlinear problem of a passive tracking low-dimensional measurement space is converted into a linear problem of a high-dimensional feature space, the problems of strong nonlinearity of an underwater multi-target measurement function and incomplete observability of a target state are solved, and the tracking precision of underwater passive tracking is improved. Meanwhile, by means of the kernel mean belief propagation method under Gaussian hypothesis, the problem that the calculation complexity is too high when random particles or Gaussian mixed density is used for approximating probability density distribution in a traditional BP algorithm is solved, the problem of matching of a measurement result and a target is solved with low calculation complexity, and the accuracy of the measurement result is improved. And the expandability of the system is greatly improved.

The application discloses a kind of based on user balance flow distribution with maximum toughness as target damaged road network recovery order optimization method, to solve the problem of low efficiency, decision lacks scientific basis after disaster traffic network recovery.The method comprises the following implementation steps: first, construct damaged road network toughness evaluation framework, in combination with traffic network topological structure and flow demand, define toughness index as the reciprocal of recovery time and total travel time product;Second, establish damaged road network toughness optimization scene parameter, abstract road network into connected graph and define recovery scheme set;Then, construct the mixed integer linear programming model with maximum toughness as target, convert nonlinear problem into linear problem by cutting line approximation, and set flow balance, travel demand and other constraint conditions;Finally, the model is globally optimized and solved using high-performance solver, to generate the best recovery scheme.The method of the application can scientifically quantify road network toughness and optimize recovery order, avoid the blindness of experience decision, quickly improve the traffic efficiency of post-disaster traffic system, provide efficient and scientific decision support for post-disaster traffic recovery, and have important practical application value.