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

29 January 2026, Volume 48 Issue 1
    

  • Select all
    |
    Articles
  • Xu Xiang, Zhao Yue
    Mathematica Numerica Sinica. 2026, 48(1): 1-29. https://doi.org/10.12286/jssx.j2025-1349
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper aims to investigate some recent progress on inverse source problems of timeharmonic wave equations and establish the stability in general cases. For scattering models of deterministic and stochastic wave equations, we show the methodology to obtain stability and summarize existing theoretical and numerical results.
  • Zhou Yanping, Chen Yanping, Hu Hanzhang, Qin Fangfang
    Mathematica Numerica Sinica. 2026, 48(1): 30-46. https://doi.org/10.12286/jssx.j2025-1284
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The compressible miscible displacement in a porous medium is widely used in many computational simulation fields of science and engineering. The mathematical model of the problem is the initial boundary value problem coupled by two parabolic partial differential equations. We apply the mixed finite element method to discrete pressure equation and the characteristic expanded mixed finite element method to discrete concentration equation. Next, we present the proofs of error estimates of the mixed finite element methodcharacteristic expanded mixed finite element method. Finally, the theoretical results are verified through numerical examples.
  • Luo Yiqing, Zhang Weihong
    Mathematica Numerica Sinica. 2026, 48(1): 47-61. https://doi.org/10.12286/jssx.j2025-1285
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    For the numerical solution of complex symmetric indefinite linear systems, we apply the minimal residual technique to improve the modified real and imaginary parts (MCRI) iteration method, and propose an iteration method referred to as minimum residual MCRI (MRMCRI) iteration method. Theoretically, the unconditional convergence of the method is proved by spectral theory, and the quasi-optimal iteration parameter is given, which is not affected by the scale or characteristics of the problem. Numerical experiments further validate the efficiency and robustness of the MRMCRI method, especially in solving problems dominated by the imaginary part of the coefficient matrix.
  • He Bo, Xu Jiawei, Li Shihai, Peng Zheng
    Mathematica Numerica Sinica. 2026, 48(1): 62-83. https://doi.org/10.12286/jssx.j2025-1287
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper studies a class of structured composite optimization problems, where the objective function is the sum of convex functions, subject to nonconvex equality constraints, with both the objective and constraint functions being continuously differentiable. Such problems have significant applications in electronic design automation, particularly in chiplet placement. Based on the Proximal-Perturbed Lagrangian (P-Lagrangian) method [Oper. Res. Lett. 51, 357–363, 2023], we propose the P-Lagrangian based Alternating Direction Method of Multipliers (PLADMM) for solving nonconvex constrained structured composite optimization problems. Under standard assumptions, we establish the convergence theory for PLADMM, proving that it converges to a KKT point. Numerical experiments demonstrate that PLADMM can effectively solve the MCNC benchmark chiplet placement problem.
  • Liu Yanru, Jia Junqing, Jiang Xiaoyun
    Mathematica Numerica Sinica. 2026, 48(1): 84-101. https://doi.org/10.12286/jssx.j2025-1288
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The fractional nonlinear Schr$\mathrm{\ddot{o}} $dinger equation is used to describe non-local phenomena in quantum physics and to explore the quantum behavior of remote interactions or time-dependent processes with multiple scales. In this paper, by using the second-order Strang time-splitting Fourier spectral method, the long-time dynamic improved uniform error bound of the spatial fractional nonlinear Schr$\mathrm{\ddot{o}} $dinger equation with small potential energy term is established. Firstly, the equations are semi-discretization by the second-order Strang time splitting method, and then the fully discretization scheme is derived by the spatial Fourier spectral method. Regularity compensation oscillation (RCO) technique is employed to prove the improved uniform error bounds at $O(\varepsilon^2 \tau^2)$ in temporal semi-discretization and $O(h^m+\varepsilon^2\tau^2)$ in full-discretization up to the long-time $T_\varepsilon=T/{\varepsilon^2} (T>0$ fixed), respectively. Finally, some numerical examples are given for convergence test and application analysis, which confirms the error bounds established and verifies the effectiveness of the numerical method presented in this paper.
  • Ma Yumin, Cai Xingju, Zhang Haiping, Wang Maoran
    Mathematica Numerica Sinica. 2026, 48(1): 102-122. https://doi.org/10.12286/jssx.j2025-1292
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper studies a class of nonconvex and nonsmooth two-block optimization problems, where the objective function consists of two nonconvex and nonsmooth separable functions and a smooth coupling function. To address such problems, we propose an improved algorithm based on the inexact inetial proximal gradient method and Nesterov’s acceleration method, namely the proximal alternating linearized minimization algorithm with two distinct extrapolation parameters. Building upon the traditional proximal alternating linearized minimization framework, the algorithm introduces distinct extrapolation parameter sequences to achieve dual extrapolation for one of the variables. Specifically, during each iteration, the algorithm constructs more tractable subproblems by linearizing the coupling function and adding proximal terms based on two different extrapolated points. Theoretically, under certain assumptions, we prove that any limit point of the bounded sequence generated by the algorithm is a critical point of the objective function. Furthermore, when the objective function satisfies the Kurdyka-Lojasiewicz property, we establish the global convergence of the algorithm. Notably, the proposed algorithm allows extrapolation parameters to take negative values, providing new possibilities for enhancing performance. To validate the effectiveness of the proposed algorithm, we apply it to solve the sparse principal component analysis problem. Extensive numerical experiments demonstrate that the proposed algorithm consistently outperforms existing methods in terms of both convergence speed and computational efficiency. Notably, when one of the extrapolation parameters is negative, the algorithm achieves further performance improvements, demonstrating its enhanced flexibility in parameter selection and underscoring its potential for broader applications.
  • Dai Shuling, Zhang Jianhua
    Mathematica Numerica Sinica. 2026, 48(1): 123-140. https://doi.org/10.12286/jssx.j2025-1306
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Randomized double and triple Kaczmarz algorithms are effective stochastic iterative methods for solving the extended normal equations $ A^{\mathsf{T}} A x = A^{\mathsf{T}} b - c $. However, their computational efficiency still has room for improvement. Based on the surrogate hyperplane projection technique, in this paper we propose the residual-based surrogate hyperplane double and triple Kaczmarz (RSHDK/RSHTK) algorithms for solving the extended normal equations. The new algorithms are applicable to arbitrary coefficient matrices $A\in \mathbb{R}^{m\times n}$ and significantly outperform the randomized double and triple Kaczmarz algorithms in terms of both iteration counts and computation time. For both consistent and inconsistent systems, we establish the convergence theories of the new algorithms and demonstrate that their convergence factors are smaller than those of the corresponding standard algorithms. Finally, numerical experiments further validate the effectiveness of the new algorithms.
  • Deng Anqi, Tang Lingyan
    Mathematica Numerica Sinica. 2026, 48(1): 141-159. https://doi.org/10.12286/jssx.j2025-1312
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Due to nonlinear schemes, the numerical performance of high-precision schemes is significantly influenced by sensitivity parameters and scale factor values, and they are unable to preserve the steady-state solution of the hyperbolic balance law equation. To address these issues, this paper first introduces a scaling function and constructs a class of scale and translation invariant fifth-order weighted compact nonlinear schemes. Then, through a pre-balancing scheme and source term splitting method, the proposed scheme is applied to the numerical solution of the shallow water equation with bottom topography source terms. Theoretical analysis shows that the new scheme's computational results are scale and translation invariant, and it can accurately preserve the moving water equilibrium of the shallow water equation. Numerical examples verify that the scheme achieves high-order accuracy, possesses good stability and balance, and can precisely capture small disturbances near the steady-state solution.
  • Wang Shihan, Yang Yang, Wang Wenqiang
    Mathematica Numerica Sinica. 2026, 48(1): 160-180. https://doi.org/10.12286/jssx.j2025-1313
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This work first establishes the well-posedness of solutions for Caputo tempered variableorder fractional stochastic differential equations (CTVO-FSDEs), including the existence, uniqueness and continuous dependence on initial conditions. An Euler-Maruyama scheme is further developed, with a rigorous proof of its strong convergence. Notably, when the fractional order reduces to a constant, our results align with existing theoretical findings in the literature. Finally, numerical simulations are presented to validate the theoretical analysis, demonstrating excellent agreement between computational results and analytical predictions.
  • Su Zhaogang, Tang Yuyang, Chen Shengjie, Chen Liang, Deng jiayi
    Mathematica Numerica Sinica. 2026, 48(1): 181-210. https://doi.org/10.12286/jssx.j2025-1330
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Unit Commitment (UC) is a core challenge in optimizing the operation of power systems. With the rapid expansion of power system scale, Mixed Integer Programming (MIP) methods are facing severe computational challenges. We propose an innovative two-stage column generation algorithm to efficiently solve the large-scale unit commitment. This method is based on Dantzig-Wolfe decomposition and reconstructs the original problem by introducing scheduling strategy variables, achieving effective decoupling between units. For the reconstructed model, we develop a two-stage computational approach: in the first stage, a parallel column generation method is used to solve the Linear Programming (LP) relaxation problem; in the second stage, based on the generated columns, the restricted master problem is solved to obtain high-quality integer feasible solutions. Numerical experiments conducted on large-scale unit commitment test cases involving 1000 to 1500 units show that the proposed method achieves over 2.5×average speedup in solving time compared to the commercial solver CPLEX. The relative optimality gap is around 0.01%, and the solving time is very stable with prominent scalability. This algorithm provides an effective computational tool for practical scheduling optimization of large-scale power systems, with both theoretical and practical significance.