中国科学院数学与系统科学研究院期刊网

14 May 2025, Volume 47 Issue 2
    

  • Select all
    |
    Reviews and Perspectives
  • Li Chenyi, Wen Zaiwen
    Mathematica Numerica Sinica. 2025, 47(2): 191-213. https://doi.org/10.12286/jssx.j2024-1273
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper provides a brief exploration of the basic principles and applications of mathematical formalization, with a focus on the formal language Lean and its application in mathematical optimization. We first review the development background of mathematical formalization, explain the construction principles of the Lean language and its correctness guarantee mechanisms, and introduce the role of the theorem library Mathlib4 in Lean. By comparing natural language with formalized expressions, we illustrate the advantages of using formalization to verify mathematics, emphasizing its important role in the accurate verification of mathematical theories. In the field of mathematical optimization, this paper discusses the current progress of formalizing mathematical optimization theory, using formalized examples of classic theorems such as the quadratic upper bound lemma, and further highlights the characteristics and advantages of formalized mathematics. Additionally, we explore the formalization goals in operations research and investigate the technologies of automated formalization and automated theorem proving, analyzing the potential and challenges of automation tools in the mathematical formalization process. Finally, we summarize the current state of research in mathematical formalization, give suggestions for further advancing the field, and discuss the significant role of formalization in the development of applied mathematical theory.
  • Articles
  • Zheng Hua, Zhang Yongxiong, Lu Xiaoping
    Mathematica Numerica Sinica. 2025, 47(2): 214-233. https://doi.org/10.12286/jssx.j2024-1198
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Vertical nonlinear complementarity problems have wide applications. The design of numerical methods for solving vertical nonlinear complementarity problems has been a hot topic among researchers in recent years. In this paper, a modulus-based synchronous multisplitting iteration method is established by matrix multisplitting and equivalent modulus equation for vertical nonlinear complementarity problem. Some convergence conditions of the proposed method are presented under H-matrix assumption. The convergence domain of the relaxation parameters of the accelerated overrelaxation iteration is obtained. By OpenACC framework, numerical tests are given to show the high parallel computational efficiency of the proposed method.
  • Zeng Ling, Chen Yumei, Xie Xiaoping
    Mathematica Numerica Sinica. 2025, 47(2): 234-254. https://doi.org/10.12286/jssx.j2024-1202
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This thesis proposes a class of space-time mixed finite element method for time fractional diffusion equations involving Riemann-Liouville derivatives of order α ∈ (0, 1). The spatial discretization employs m(m ≥ 0) order Raviart-Thomas (RT) finite elements, while the temporal discretization uses piecewise r(r ≥ 0) degree discontinuous Galerkin (DG) finite elements. To address the solution singularity near t = 0, graded meshes are used in the time direction. The well-posedness of the fully discrete scheme is analyzed. Error estimates are derived in two cases: r = 0, i.e. the temporal discretization uses the piecewise constant DG scheme, and r = 1, i.e. the temporal discretization uses the piecewise linear DG scheme. Numerical experiments are provided.
  • Chen Xin, Qin Yuefeng, Zhou Xuelin, Li Jiaofen
    Mathematica Numerica Sinica. 2025, 47(2): 255-284. https://doi.org/10.12286/jssx.j2024-1213
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Multidimensional scaling (MDS) is a technique used in multidimensional data analysis that depicts the similarities or relationships between observed objects as distances between points in a lower-dimensional space. By representing high-dimensional data within a lowdimensional framework, MDS preserves the relative distances between data points. This study focuses on developing an efficient numerical algorithm for a specific type of individual differences scaling model, known as O-INDSCAL, within symmetric multidimensional scaling, which accounts for individual differences among observed objects. Initially, leveraging the concept of the alternating least squares algorithm, the multivariable constrained matrix optimization model associated with the O-INDSCAL model is transformed into a fixed-point iteration problem. By thoroughly examining the acceleration principles and implementation processes of various polynomial extrapolation and Anderson acceleration methods in vector sequence acceleration, we have designed a matrix sequence acceleration algorithm tailored to this problem model. Numerical experiments indicate that the proposed acceleration algorithms significantly enhance the convergence speed of sequences generated by fixed-point iterations. Furthermore, when compared with existing continuous-time projection gradient flow algorithms and the first-order and second-order Riemannian algorithms available in the Manopt toolbox for manifold optimization, the proposed algorithms demonstrate notable improvements in iterative efficiency.
  • Hu Xinghua, Wang Chi
    Mathematica Numerica Sinica. 2025, 47(2): 285-303. https://doi.org/10.12286/jssx.j2024-1214
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, the Hermite polynomials are used as the hidden layer of the neural network, the initial weights of the Hermite neural network are optimized using genetic algorithm, and the inverse of the error function between the actual output and the desired output of the Hermite neural network optimized by the genetic algorithm is also chosen as the fitness function of the genetic algorithm, to construct a new genetic algorithm-optimized Hermite neural network to solve the Caputo fractal-fractional order Bagley-Torvik differential equation numerical method. The general form of the numerical solution of the Caputo fractal-fractional order Bagley-Torvik differential equation is given in conjunction with the Taylor’s formula at multiple points, and the absolute error and convergence of the algorithm are theoretically investigated. Comparison with existing numerical methods is made and the results show the effectiveness and feasibility of the method in this paper.
  • Yu Haifang, Gao Jianfang
    Mathematica Numerica Sinica. 2025, 47(2): 304-314. https://doi.org/10.12286/jssx.j2024-1215
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper mainly discusses the oscillation of numerical solutions for a class of neutral delay differential equations with multiple delays. The original equation is discretized by using the linear θ-method to obtain the corresponding difference equation. By discussing the properties of the solutions of the difference equation, the oscillation properties of the numerical solutions for the original equation are transformed into that of a non-neutral difference equation. According to the relationship between the oscillation of the difference equation and the characteristic roots of the characteristic equation, the oscillation of the numerical solutions is discussed for $0\leq \theta \leq \frac{1}{2}$ and $\frac{1}{2}<\theta \leq1$, respectively. Meanwhile, the properties of non-oscillatory numerical solutions are also investigated. Finally, numerical examples are given to illustrate the conclusions.
  • Gao Xiaonan, Long Xianjun
    Mathematica Numerica Sinica. 2025, 47(2): 315-325. https://doi.org/10.12286/jssx.j2024-1218
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, we propose an adaptive step-size rule three-operator splitting algorithm with an inertia term to solve the nonsmooth DC programming problem. Under suitable assumptions, we prove that the sequence generated by the algorithm converges to a critical point of the problem. Finally, we apply the algorithm for solving the sparse recovery problem, and numerical experiments show the effectiveness and superiority of the new algorithm.
  • Zhu Wenchang
    Mathematica Numerica Sinica. 2025, 47(2): 326-346. https://doi.org/10.12286/jssx.j2024-1221
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Single-phase multi-component flow problems are widely existed in oil and gas reservoirs, groundwater pollution, etc., and numerical simulation is of great significance. In this paper, for the traditional IMPEC (Implicit Pressure-Explicit Concentration) method, which does not follow the problem of mass conservation for all components, a physically preserving IMPEC method with mass conservation for each component is used for the time discretization, and the upwind block-centered finite-difference method is used for the spatial discretization of the pressure equation, Darcy’s equation, and the component equations. In this paper, the physics-preserving mass conservation properties of all the constructed components as well as the molar concentrations are rigorously proved under reasonable condition, respectively. Finally, numerical simulations are carried out to demonstrate the validity of the proposed algorithm by means of a single-phase multi-component flow problem.
  • Ji Min, Li Hong
    Mathematica Numerica Sinica. 2025, 47(2): 347-362. https://doi.org/10.12286/jssx.j2024-1222
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    A compact difference scheme for the two-dimensional Sobolev equation with fourth-order accuracy in space is derived, and the convergence of the compact difference scheme is proved. The difference scheme is rewritten in vector form, a reduced-order high-order compact difference scheme is constructed using the Proper Orthogonal Decomposition (POD) method, and error estimate for the approximate solution is also provided. Numerical example is presented to calculate the numerical error, spatial convergence order and temporal convergence order of both the compact difference scheme and the reduced-order compact difference scheme, verifying that the experimental results are consistent with the theoretical analysis. Furthermore, by comparing the CPU computation time before and after dimension reduction, the superiority of applying the POD method for dimension reduction in compact schemes is demonstrated.
  • Zeng Yu, Xu Weiru, Hu Wenyu
    Mathematica Numerica Sinica. 2025, 47(2): 363-384. https://doi.org/10.12286/jssx.j2024-1231
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, we consider a special kind of matrix combined by a skew-symmetric tridiagonal matrix and a periodic skew-symmetric tridiagonal matrix. It is referred to as a generalized periodic skew-symmetric tridiagonal matrix. An inverse eigenvalue problem for constructing thus a matrix is studied, that is, the matrix is reconstructed from three prescribed balanced sets and a positive number. Firstly, the matrix can be transformed into a generalized periodic symmetric tridiagonal matrix by unitary similarity, and then the relationships between the eigenvalues of this matrix and those of its leading and trailing principal submatrices are analyzed. Two aspects are respectively discussed about whether the spectra of these two principal submatrices have a disjoint set and the parity of the order of the above trailing principal submatrix. Then, the necessary and sufficient conditions for the existence of solutions to the inverse eigenvalue problems in various cases are provided, and the largest number of solutions and reconstruction algorithms are determined. Finally, the effectiveness of the algorithm is validated by two numerical examples.