Nonlinear multi-scale dynamic system control mechanism inference method based on Koopman operator theory

Abstract

The invention discloses a nonlinear multi-scale dynamic system control mechanism inference method based on a Koopman operator theory, and the method comprises the steps: building an order reduction mechanism model of a dynamic system from observed time series data based on the Koopman operator theory; firstly, dynamic sampling is carried out on observation time sequence data based on a self-adaptive sampling strategy designed by the method, the sampling rate is adjusted according to data changes, and information of different time scales is effectively separated. Thirdly, performing local data smoothing and numerical differentiation on the discrete sampling points by applying cubic B-spline interpolation; on the basis, a Block-Hankel matrix is innovatively constructed by using the nonlinear modeling capability of a Koopman operator theory so as to capture the multivariate dynamic characteristics of the system; and a main dynamic mode is extracted by using a kernel singular value decomposition technology, and system identification is expanded to a nonlinear power system. And finally, the SINDy algorithm is used for identifying the reduced-order model, and accurate modeling of the complex nonlinear multi-scale system is realized. According to the method, the advantages of Koopman embedding and sparse regression are fused on the basis of the multi-scale data processing capability of adaptive sampling, so that the method has relatively high robustness and universality, and has theoretical and practical significance for mechanism modeling of a nonlinear and multi-scale system.

Description

Technical Field

[0001] The invention belongs to the technical field of nonlinear, multi-scale dynamic system mechanism modeling, and specifically is a nonlinear multi-scale dynamic system control mechanism inference method based on Koopman operator theory. Background Art

[0002] In the increasingly complex industrial field, systems often exhibit various nonlinear behaviors under different operating conditions and show unique dynamics across multiple time scales. The nonlinear and multi-scale properties of industrial systems pose significant challenges to understanding and predicting their dynamics. Deep learning algorithms based on neural networks have made significant progress in predicting industrial parameters. However, the rapid development of technology has led to the pursuit of more complex and sophisticated analysis in industrial process modeling. Traditional deep learning models, as black boxes, have difficulty revealing the inherent mechanisms of dynamic systems, which limits their application in applications requiring high-precision interpretation capabilities. On the other hand, traditional physical models and mathematical methods based on first principles have also encountered difficulties in simulating complex industrial systems, especially in the absence of a comprehensive theoretical framework and precise parameters. Therefore, researchers are exploring new methods with the goal of establishing models that can describe the essence of industrial processes based on observed time series data, providing more accurate theoretical support for the analysis, design, and optimization of industrial systems.

[0003] In recent years, extensive research efforts have focused on data-driven identification of nonlinear systems, advancing the research progress of both nonparametric and parametric models. Symbolic regression based on genetic algorithms has been widely studied in various fields, but this approach may overfit the target system with a large number of false positives, resulting in high computational costs and identification errors. Sparse identification of nonlinear dynamics (SINDy) based methods have been widely used to reveal the physical principles of low-dimensional systems, but the application of this approach to industrial systems with complex nonlinear characteristics is still in the exploratory stage.

[0004] The study showed that the Koopman operator-based analytical framework shows great potential in dealing with nonlinear systems. When combined with the sparsity of SINDy, this approach can identify the most influential low-dimensional variables that affect the system dynamics, thereby building a more compact reduced-order model. However, this combination still faces many challenges in modeling multiscale systems, especially in dealing with uneven data sampling and cross-scale dynamics. Therefore, a system mechanism modeling method that can simultaneously be nonlinear and multi-scale is of great research significance. Summary of the invention

[0005] In view of the defects or deficiencies of the prior art described above, the present invention provides a control mechanism inference method for nonlinear multi-scale dynamic systems based on Koopman operator theory, which is used for data-driven modeling of dynamic systems, aiming to generate a generalizable and interpretable reduced-order model. First, the observed time series data is downsampled based on an adaptive sampling strategy, which dynamically adjusts the sampling rate according to data changes and effectively isolates information of different scales. Subsequently, cubic B-spline interpolation is applied for local data smoothing and numerical differentiation. Utilizing the nonlinear modeling capability of Koopman operator theory, system identification is extended to nonlinear dynamic systems by implementing kernel SVD on the newly constructed Block-Hankel matrix. In addition, the integrated enhanced SINDy method is used to achieve accurate modeling of complex systems.

[0006] The technical solution adopted by the present invention to achieve the above-mentioned purpose is:

[0007] A mechanism modeling method for nonlinear, multi-scale dynamic systems based on Koopman embedding theory. This method effectively separates multi-time scale information through a designed adaptive sampling strategy and uses the nonlinear modeling ability of Koopman operator theory to generate a generalizable and interpretable reduced-order model from the observed time series data. The method includes the following steps:

[0008] S1. Obtain system observation data;

[0009] S2, perform adaptive sampling according to the dynamic changes of data to effectively capture multi-scale information;

[0010] S3, the sampled data {t i ,x i}Apply cubic B-spline interpolation for local smoothing;

[0011] S4, analyze the dynamic characteristics of the system, construct the Block-Hankel matrix for the interpolated time series data, and use its high-dimensional structure to extract the spatiotemporal characteristics of the system;

[0012] S5. Decompose the constructed Block-Hankel matrix using Kernel SVD to extract the main dynamic mode and characteristic subspace of the system;

[0013] S6. Combine the extracted low-dimensional feature subspace with the enhanced SINDy method to construct a sparse dynamic model Used to reduce the order of the model by calculation Obtain an approximate representation of the Koopman operator that characterizes the incremental value of the next time point of the current attribute data, and superimpose it on the value of the previous time point as a predicted value of the future state or dynamic characteristics of the current attribute data;

[0014] Step 7: Use the obtained reduced-order mechanism model to analyze the future state changes of system data information, so as to achieve precise control.

[0015] The observation data is actually collected long-term historical observation information, or time series information of the system evolving over time simulated in a simulation environment under a variety of different initial conditions.

[0016] The adaptive sampling is to dynamically adjust the sampling rate within a range according to the change rate of the data values ​​of the continuous sampling points of the observed data, and obtain the long time series sampling point data {t i ,x i}.

[0017] The cubic B-spline interpolation includes: for a given data point {t i ,x i}(i=0,1,...,n), by recursively defining the B-spline basis function N i (t), we get a piecewise cubic B-spline polynomial Use its smooth interpolation data point t i , which is used to maintain the smoothness and stability of the curve while ensuring the interpolation accuracy.

[0018] The recursive definition of the B-spline basis function N i (t) specifically include:

[0019]

[0020] where t i are the elements of the knot vector. For cubic B-spline k=3, the above recursive process produces cubic polynomial B-spline basis functions;

[0021] The interpolation function S(t) is a linear combination of basis functions and is expressed as:

[0022]

[0023] where c i is the control point coefficient to be determined, which is solved by introducing boundary conditions;

[0024] To ensure that the interpolation curve S(t) passes through the given data point x i , the following conditions must be met:

[0025] S(t i )=x i ,for i=0,1,...,n.

[0026] Therefore, the first-order derivative of the cubic B-spline basis function is expressed as:

[0027]

[0028] The method of constructing a Block-Hankel matrix for the interpolated time series data includes: constructing a new Block-Hankel matrix different from the traditional standard embedding matrix by introducing a window length L The Block-Hankel matrix is the Hankel matrix Vertical connection formation; the Block-Hankel matrix It is not only used to capture the dynamic characteristics of each time series, but also preserves the interdependencies between different time series.

[0029] The Hankel matrix for:

[0030]

[0031] Where K = N-L+1, N is the total length of the time series.

[0032] Decomposing the constructed Block-Hankel matrix by using kernel singular value decomposition Kernel SVD includes: using Gaussian kernel function Map the Block-Hankel matrix to a high-dimensional feature space:

[0033]

[0034] in and are any two columns of the Block-Hankel matrix;

[0035] Kernel matrix K pq Perform SVD decomposition:

[0036] K=UΣV T

[0037] Then, after KSVD decomposition, the eigenvector matrix V is obtained, which consists of a set of principal singular vectors in the kernel space and is used to represent the principal dynamic mode and the eigensubspace.

[0038] The SINDy algorithm specifically includes:

[0039] Construct a function library Θ(V) that encapsulates the eigenvector matrix V using one of multiple potential dynamical function forms;

[0040]

[0041] Where m is the number of time points and p is the number of functions in the library;

[0042] Obtain a reduced-order model of the system through sparse regression:

[0043]

[0044] Reduced Order Model The calculation result of approximately represents the Koopman operator of the incremental value of the next time point of the current attribute data.

[0045] The dynamic function forms include polynomials and trigonometric functions.

[0046] The present invention has the following beneficial effects and advantages:

[0047] A new mechanism modeling method based on Koopman embedding theory is used to establish a reduced-order mechanism model of nonlinear, multiscale dynamic systems based on noisy multidimensional observational time series data. From the perspective of separation time scales, the superiority of this method is mainly reflected in the designed dynamic sampling strategy, which adjusts the sampling rate according to the rate of change (gradient) of the data. That is, in areas where the data changes drastically (high gradient areas), the sampling rate is increased to capture more details and changes. In areas where the data changes gently (low gradient areas), reducing the sampling rate and redundant data can improve the efficiency of data processing. This dynamic sampling can adaptively adjust the sampling strategy according to the data characteristics without the need to set a fixed sampling rate in advance, and is suitable for data with various different characteristics.

[0048] In addition, from the perspective of dynamic extraction of nonlinear systems, cubic B-spline interpolation is particularly effective for fitting such systems due to its smoothness and high-order continuity, avoiding the sawtooth effect that may occur in general interpolation methods such as linear interpolation. The innovative construction of the Block-Hankel matrix retains the interrelationships between different time series and captures the dynamic characteristics of multivariate systems. By applying kernel technology, the Block-Hankel matrix is ​​nonlinearly expanded, allowing kernel singular value decomposition to retain the main dynamic information of the system in high-dimensional space. Constructing the SINDy model on the extracted features can reduce the order modeling of nonlinear systems and reduce the impact of noise on the model to a certain extent, thereby improving the accuracy and robustness of the model, so that it can more accurately guide system predictive control. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 The present invention is a flow chart of the method.

[0050] Figure 2 Schematic diagram of a jacketed continuous stirred tank reactor system.

[0051] Figure 3 It is a comparison chart of the real data of the non-isothermal CSTR system with a jacket, the SINDy results and the results of the method proposed in the present invention.

[0052] Figure 4 It is a comparison chart of the real data of the multi-reaction isothermal batch reactor system, the SINDy results and the results of the method proposed in the present invention.

[0053] Figure 5 It is a comparison chart of the SINDy result of the actual grinding-classification system and the mean square error result of the method proposed in the present invention. DETAILED DESCRIPTION

[0054] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the specific implementation method of the present invention is described in detail below in conjunction with the accompanying drawings. In the following description, many specific details are set forth to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways different from those described herein, and those skilled in the art can make similar improvements without violating the connotation of the invention, so the present invention is not limited by the specific implementation disclosed below.

[0055] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art of the present invention. The terms used in the specification of the invention herein are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The present invention is further described in detail below in conjunction with the accompanying drawings and examples.

[0056] The present invention is a unified framework for system mechanism modeling based on the control requirements of industrial systems. The present invention has conducted experiments on two simulated reactor systems and a real grinding and classification process. The experimental results all show that the present invention can accurately model nonlinear, multi-scale complex systems, and has strong robustness and universality. Therefore, the present invention can be extended to general industrial systems, establish a reduced-order mechanism model to describe the system, and guide actual industrial prediction and control.

[0057] like Figure 1As shown, it is a flowchart of an example method of the present invention. A control mechanism inference method for a nonlinear multi-scale dynamic system based on the Koopman operator theory. First, the observed time series data is dynamically sampled based on the adaptive sampling strategy designed by the present invention, and the sampling rate is adjusted according to the data change to effectively separate information of different time scales. Subsequently, the discrete sampling points are locally smoothed and numerically differentiated using cubic B-spline interpolation. On this basis, the nonlinear modeling capability of the Koopman operator theory is utilized to innovatively construct a Block-Hankel matrix to capture the multivariate dynamic characteristics of the system; and the kernel singular value decomposition technique is used to extract the main dynamic modes, extending the system identification to nonlinear dynamic systems. Finally, the SINDy algorithm is used to identify the reduced-order model to achieve accurate modeling of complex nonlinear multi-scale systems. The programming language used in the program execution step of the present invention is not limited to MATLAB, Python, etc.

[0058] The specific steps of the present invention are as follows:

[0059] Step 1: Obtain system observation data.

[0060] Example: A continuous stirred tank reactor (CSTR) is a chemical system widely used in process industries. The present invention considers a non-isothermal CSTR in which a first order exothermic irreversible reaction occurs, which is controlled by a rate constant k. The schematic diagram of the reaction system is shown in Figure 2 shown.

[0061] In this process, the meaning and corresponding values ​​of each parameter are listed in Table 1 below. Based on these parameters, the material balance of substance A and the energy balance of the reactor and jacket can be formulated as follows:

[0062]

[0063]

[0064] where the reaction rate r is given by:

[0065]

[0066] The sampling details of the system are summarized in Table 2. Inspection of the sampled datasets shows that C A The concentration changes slowly, while the temperature variable T r and T j Shows rapid dynamic changes.

[0067] Table 1 Parameter values ​​and meanings of the jacketed continuous stirred tank reactor system

[0068]

[0069]

[0070] Table 2 Sampling details of the jacketed continuous stirred tank reactor system

[0071]

[0072] The method of the present invention can be used to reduce the concentration C of the reactant A in the reaction system. A , reactor temperature T r and jacket temperature T j Make predictions.

[0073] Step 2: Perform adaptive sampling to dynamically adjust the sampling rate according to the dynamic changes of the data to ensure that multi-scale information can be effectively captured;

[0074] Step 2-1, for data x(t), calculate the rate of change as follows:

[0075]

[0076] where x t and x t+1 are the data values ​​of two consecutive sampling points.

[0077] Step 2-2, set the initial sampling rate r0, and define the gradient threshold θ to trigger the sampling rate adjustment. If g>θ, increase the sampling rate r. Otherwise, reduce the sampling rate r. The sampling rate adjustment function f(r,g) can be defined as:

[0078] r new =r old ·f(r old ,g)=r old ·(1+α·(g-θ)) (6)

[0079] Where α is the adjustment coefficient that controls the sensitivity of sampling rate changes. Record the index of the selected sampling point to form a set of sampling points t s .

[0080] Step 3: Apply cubic B-spline interpolation to the sampled data for local smoothing;

[0081] For the sampled data points {t i ,x i}(i=0,1,...,n), define the B-spline basis function N i (t):

[0082]

[0083] where t iis the element of the knot vector. For cubic B-spline, using k=3, the above recursive process produces cubic polynomial B-spline basis functions. The interpolation function S(t) as a linear combination of basis functions can be expressed as:

[0084]

[0085] where c i is the control point coefficient to be determined.

[0086] In order to ensure that the interpolation curve S(t) passes through the given data points {t i ,x i}, the following conditions must be met:

[0087] S(t i )=x i ,for i=0,1,...,n.(10)

[0088] By introducing natural boundary conditions (the second-order derivative of the curve at the boundary is zero), the control point coefficient c can be solved i .

[0089] The cubic B-spline interpolation has a significant advantage, namely that it provides C 2 Continuity, which means that the interpolation curve is analytically differentiable and has continuous first-order and second-order derivatives across the segment. Its first-order derivative can be expressed as:

[0090]

[0091] Step 4: Construct a Block-Hankel matrix based on the interpolated time series data, and use its high-dimensional structure to extract the spatiotemporal characteristics of the system;

[0092] The traditional Hankel alternative view of Koopman analysis method (HAVOK) includes:

[0093] The Koopman operator theory and time series analysis technology are combined to reveal the dynamic characteristics of the system by constructing the Hankel matrix. The core idea of ​​this method is to use the linear characteristics of the Koopman operator to analyze and predict the behavior of nonlinear dynamic systems. Specifically, the Koopman operator is an infinite-dimensional linear operator acting on the state space of the dynamic system, which transforms the nonlinear evolution of the system into a linear evolution in the function space. t+1 =f(x t ) describes a nonlinear dynamic system, Koopman operator The Koopman operator operates on the observable function g and provides the evolution of the function in the state space. The operation of the Koopman operator is as follows:

[0094]

[0095] For time series data For the generated dynamic system, constructing the Hankel matrix requires choosing an embedding dimension m such that The rows of the Hankel matrix consist of consecutive subsequences in the time series, and the columns are formed by delayed copies of the subsequences. This design effectively captures the temporal dependencies present in the time series. The standard Hankel matrix is ​​defined as:

[0096]

[0097] The traditional HAVOK algorithm uses singular value decomposition (SVD) to analyze the Hankel matrix and decompose it into three parts:

[0098] H=U∑V * (14)

[0099] Where U and V are orthogonal matrices, ∑ is a diagonal matrix whose diagonal elements correspond to the singular values ​​of the Hankel matrix. SVD is used to select the largest singular value and the corresponding singular vector, which form the basis of the Koopman invariant subspace. By analyzing the columns of the matrix V, the main dynamic modes of the system can be identified. When the Koopman operator is applied, these main modes are retained, forming a Koopman invariant subspace. If r is the number of singular values ​​retained, the traditional standard embedding matrix HAVOK model can be approximated as:

[0100]

[0101] Here, U r ,∑ r and V r is a matrix consisting of the first r singular values ​​and their corresponding singular vectors.

[0102] The present invention constructs a new Block-Hankel matrix different from the traditional standard embedding matrix by introducing an appropriate window length L:

[0103]

[0104] Where K = N-L+1, N is the total length of the time series. Hankel matrix Connect vertically to form a Block-Hankel matrix This Block-Hankel matrix not only captures the dynamic characteristics of each time series, but also preserves the interdependencies between different time series.

[0105] Step 5: Decompose the constructed Block-Hankel matrix using kernel singular value decomposition (Kernel SVD) to extract the main dynamic mode and characteristic subspace of the system;

[0106] Step 5-1: Use Gaussian kernel function Map the Block-Hankel matrix to a high-dimensional feature space:

[0107]

[0108] in and are any two columns of the Block-Hankel matrix.

[0109] Step 5-2: Kernel matrix K pq Perform SVD:

[0110] K=UΣV T (18)

[0111] After KSVD decomposition, the eigenvector matrix V is obtained, which consists of a set of principal singular vectors in the kernel space, corresponding to the main dynamic modes in the data, and is also regarded as the basis of the Koopman operator approximate invariant subspace. By utilizing these principal singular vectors, the SINDy algorithm can be used to identify the functions associated with these principal delay coordinates.

[0112] Step 6: Combine the extracted low-dimensional feature subspace with the enhanced SINDy method to construct a sparse dynamic model.

[0113] Step 6-1. Construct a function library Θ(V) containing various potential dynamic forms, such as polynomials, trigonometric functions, etc.:

[0114]

[0115] Where m is the number of time points and p is the number of functions in the library.

[0116] Step 6-2: Obtain the reduced-order model of the system through sparse regression:

[0117]

[0118] The resulting model constitutes an approximate representation of the Koopman operator and can be used to predict the future state of the system or analyze its dynamic characteristics.

[0119] Thus, the present invention obtains a reduced-order mechanism model that accurately describes the system. Figure 3 The real data of the jacketed non-isothermal CSTR system, the SINDy results and the results of the method proposed in the present invention are compared. In addition to this example, the present invention also tests a simulated multi-reaction isothermal batch reactor system and a real grinding-classification system. Figure 4 It is a comparison chart of the real data of the multi-reaction isothermal batch reactor system, the SINDy results and the results of the method proposed in the present invention. Figure 5 The figure is a comparison of the SINDy result of the real grinding-classification system and the mean square error result of the method proposed in the present invention. It can be seen that the present invention contributes to the efficient and accurate modeling process and effectively captures the dynamic characteristics of the system.

[0120] Step 7: Analyze the system state changes based on the obtained reduced-order mechanism model to achieve precise control of the system.

[0121] The present invention establishes a reduced-order mechanism model that describes the dynamic characteristics of nonlinear, multi-scale complex systems from observed time series data. The model is universal and can be generalized to the field of actual engineering control.

[0122] In summary, in solving the complex and challenging issue of nonlinear multi-scale industrial system modeling, the present invention proposes a new method for inferring the control mechanism of the system. The present invention adopts dynamic sampling technology to adaptively adjust the sampling strategy according to the data characteristics, effectively isolating the dynamics of different time scales. In order to obtain continuous time series data and its corresponding derivatives from discrete sample points, the present invention uses cubic B-spline interpolation to achieve high-order smoothing, providing more accurate data support for subsequent model identification. A new Block-Hankel matrix is ​​constructed, and kernel SVD analysis is applied to effectively capture the dynamic characteristics of the system. This study uses the SINDy method to identify reduced-order dynamic models, so that complex nonlinear multi-scale systems can be accurately modeled. The reduced-order model obtained by the present invention can guide actual industrial control.

[0123] The above description is only a preferred embodiment of the present invention and does not limit the present invention in any way. Any simple modification, change and equivalent structural change made to the above embodiment based on the technical essence of the present invention still falls within the protection scope of the technical solution of the present invention.

Claims

1. A mechanism modeling method for nonlinear, multi-scale dynamic systems based on Koopman embedding theory, characterized in that: This method effectively separates multi-time scale information through a designed adaptive sampling strategy and uses the nonlinear modeling capability of Koopman operator theory to generate a generalizable and interpretable reduced-order model from the observed time series data. The method includes the following steps: S1. Obtain system observation data; S2, perform adaptive sampling according to the dynamic changes of data to effectively capture multi-scale information; S3, the sampled data {t i ,x i }Apply cubic B-spline interpolation for local smoothing; S4, analyze the dynamic characteristics of the system, construct the Block-Hankel matrix for the interpolated time series data, and use its high-dimensional structure to extract the spatiotemporal characteristics of the system; S5. Decompose the constructed Block-Hankel matrix using Kernel SVD to extract the main dynamic mode and characteristic subspace of the system; S6. Combine the extracted low-dimensional feature subspace with the enhanced SINDy method to construct a sparse dynamic The state model V is used to obtain an approximate representation of the Koopman operator representing the incremental value of the next time point of the current attribute data by calculating the reduced-order model V, and superimpose the value of the previous time point on it as the predicted value of the future state or dynamic characteristics of the current attribute data; Step 7: Use the obtained reduced-order mechanism model to analyze the future state changes of system data information, so as to achieve precise control.

2. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The observation data is actually collected long-term historical observation information, or time series information of the system evolving over time simulated in a simulation environment under a variety of different initial conditions.

3. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The adaptive sampling is to dynamically adjust the sampling rate within a range according to the change rate of the data values ​​of the continuous sampling points of the observed data, and obtain the long time series sampling point data {t i ,x i }.

4. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The cubic B-spline interpolation includes: for a given data point {t i ,x i }(i=0,1,...,n), by recursively defining the B-spline basis function N i (t), we get a piecewise cubic B-spline polynomial Use its smooth interpolation data point t i , which is used to maintain the smoothness and stability of the curve while ensuring the interpolation accuracy.

5. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 4, characterized in that: The recursive definition of the B-spline basis function N i (t) specifically include: where t i are the elements of the knot vector. For cubic B-spline k=3, the above recursive process produces cubic polynomial B-spline basis functions; The interpolation function S(t) is a linear combination of basis functions and is expressed as: where c i is the control point coefficient to be determined, which is solved by introducing boundary conditions; To ensure that the interpolation curve S(t) passes through the given data point x i , the following conditions must be met: S(t i )=x i ,for i=0,1,...,n. Therefore, the first-order derivative of the cubic B-spline basis function is expressed as:

6. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The Block-Hankel matrix is ​​constructed for the interpolated time series data, including: constructing a new Block-Hankel matrix H different from the traditional standard embedding matrix by introducing a window length L, wherein the Block-Hankel matrix is the Hankel matrix Vertical connection formation; the Block-Hankel matrix It is not only used to capture the dynamic characteristics of each time series, but also preserves the interdependencies between different time series.

7. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 6, characterized in that: The Hankel matrix for: Where K = N-L+1, N is the total length of the time series.

8. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: Decomposing the constructed Block-Hankel matrix by using kernel singular value decomposition Kernel SVD includes: using Gaussian kernel function Map the Block-Hankel matrix to a high-dimensional feature space: in and are any two columns of the Block-Hankel matrix; Kernel matrix K pq Perform SVD decomposition: K=UΣV T Then, after KSVD decomposition, the eigenvector matrix V is obtained, which consists of a set of principal singular vectors in the kernel space and is used to represent the principal dynamic mode and the eigensubspace.

9. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The SINDy algorithm specifically includes: Construct a function library Θ(V) that encapsulates the eigenvector matrix V using one of multiple potential dynamical function forms; Where m is the number of time points and p is the number of functions in the library; Obtain a reduced-order model of the system through sparse regression: Reduced Order Model The calculation result of approximately represents the Koopman operator of the incremental value of the next time point of the current attribute data.

10. The mechanism modeling method for nonlinear, multi-scale dynamic systems according to claim 1, characterized in that: The dynamic function forms include polynomials and trigonometric functions.

Images