用随机奇异值分解算法求解矩阵恢复问题

本文研究了大型低秩矩阵恢复问题.利用随机奇异值分解（RSVD）算法,对稀疏矩阵做奇异值分解.该算法与Lanczos方法相比,在误差精度一致的同时运算时间大大降低,且该算法对相对低秩矩阵也有效.     

用随机奇异值分解算法求解矩阵恢复问题（英文）

本文研究了大型低秩矩阵恢复问题.利用随机奇异值分解(RSVD)算法,对稀疏矩阵做奇异值分解.该算法与Lanczos方法相比,在误差精度一致的同时运算时间大大降低,且该算法对相对低秩矩阵也有效.     

非线性发展方程的小模板简化Pade格式

在有理逼近的紧致格式的理论基础上，采用特别的统一的Pade逼近形式，构造了针对高阶非线性发展方程的、简单小模板的差商格式．不仅保持了格式的四阶精度，而且还可以采用追赶法求解得到的3对角矩阵，或者采用三阶Runge-Kuna法直接求解积分．计算效果通过多种算例表明是十分令人满意的．相对于其他差分格式，此方法具有模板较小而精度保持四阶的优点．     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

一类min-max-min问题的区间算法

讨论了一类由一阶连续可微函数构成的无约束min-max-min问题．通过构造目标函数的区间扩张、无解区域删除原则，建立了求解min-max-min问题的区间算法，证明了算法的收敛性，给出了数值算例．理论证明和数值结果表明方法是可靠和有效的．     

部分变量带非负约束的严格凸二次规划问题的新算法

本文将正交校正共轭梯度法推广来解只有部分变量带非负约束而其它变量无约束的严格凸二次规划，所建立的新算法的优点是：在迭代过程中，不用求逆矩阵，这样能保持矩阵的稀疏性，数值结果表明：算法对大规模稀疏二次规划问题是可行和有效的．     

交替方向的扩散方程精细积分并行算法

给出了交替方向的二维扩散方程的精细积分算法，将一个时间步积分分为两个方向，使大规模矩阵的计算转化为一些小矩阵的计算，减小了每一步求解的计算量．对于方形区域的齐次方程，计算结果与全城精细积分完全相同，而计算量和存储量都要小得多．算例表明了算法具有较高的并行计算加速比和计算效率．     

三维热传导方程的高精度有限差分方法

提出了数值求解三维热传导方程的一个四阶精度的有限差分格式,首先对三个空间方向上的二阶导数项,采用四次样条函数来近似,从而得到半离散的常微分方程.然后利用常微分方程的解析解表达式,时间矩阵利用Padé近似,得到时间和空间均为四阶精度的差分格式.最后利用方法计算了两个数值算例,并与文献中结果进行了对比,从而验证了高精度格式的性能.     

两个带变时间步的两重网格算法的稳定性分析

对用于求解非线性发展方程的两个带变时间步的两重网格算法，对空间变量用有限元离散，对时间变量分别用一阶精度Euler显式和二阶精度半隐式差分格式离散，然后构造两重网格算法，通过深入的稳定性分析，得出本文的算法优于标准全离散有限元算法。     

Winograd矩阵乘法算法用于任意阶矩阵时的一种新处理方法

摘要t矩阵乘法StraSsen算法及其变形winograd算法用分而治之的方法把矩阵乘法时间复杂性由传统的D(n。)改进到0(佗kg。n.但是对于奇数阶矩阵,在划分子矩阵时,要作特殊处理才能继续使用此算法.本文提出了一种非等阶“十”字架划分方法,可以最少化填零,最大化性能,使得奇数阶矩阵乘法的时间复杂性更加接近偶数阶矩阵乘法的效果.计算实例显示该方法是有效的.     

四阶杆振动方程的两类高精度辛格式

In this paper, we present two classes of symplectic schemes with high order accuracy for solving four-order rod vibration equation utt uxxxx=0 via the third type generating function method. First, the equation of four order rod vibration is written into the canonical Hamilton system; second, overcoming successfully the essential difficult on the calculus of high order variations derivative, we get the semi-discretization with arbitrary order of accuracy in time direction for the PDEs by the third type generating function method. Furthermore the discretization of the related modified equation of original equation is obtained. Finally, arbitrary order accuracy symplectic schemes are obtained. Numerical results are also presented to show the effectiveness of the scheme, high order accuracy and properties of excellent long-time numerical behavior.     

Numerical computation of branch points in nonlinear equations

Summary The numerical computation of branch points in systems of nonlinear equations is considered. A direct method is presented which requires the solution of one equation only. The branch points are indicated by suitable testfunctions. Numerical results of three examples are given.     

Connection coefficients on an interval and wavelet solutions of Burgers equation

A definition of connection coefficients is introduced and techniques of computation are presented. We use semi-implicit time difference scheme to solve Burgers equation by applying the evaluations of connection coefficients in calculating the integrals of the variational form. Comparisons of accuracy and robustness of numerical solutions are mentioned in the examples.     

Finite difference spectral approximation for the time‐space fractional telegraph equation and its parameter estimation

The object of this paper is to present the numerical solution of the time‐space fractional telegraph equation. The proposed method is based on the finite difference scheme in temporal direction and Fourier spectral method in spatial direction. The fast Fourier transform (FFT) technique is applied to practical computation. The stability and convergence analysis are strictly proven, which shows that this method is stable and convergent with (2?α) order accuracy in time and spectral accuracy in space. Moreover, the Levenberg‐Marquardt (L‐M) iterative method is employed for the parameter estimation. Finally, some numerical examples are given to confirm the theoretical analysis.     

Crank–Nicolson Finite Element Scheme and Modified Reduced-Order Scheme for Fractional Sobolev Equation

Abstract

In this paper, a Crank–Nicolson finite difference/finite element method is considered to obtain the numerical solution for a time fractional Sobolev equation. Firstly, the classical finite element method is presented. Stability and error estimation for the fully discrete scheme are rigorously established. However, the amount of calculation and computing time are too large due to many degrees of freedom of classical finite element scheme and nonlocality of fractional differential operator. And then the modified reduced-order finite element scheme with low dimensions and sufficiently high accuracy, which is based on proper orthogonal decomposition technique, is provided. Stability and convergence for the reduced-order scheme are also studied. At last, numerical examples show that the results of numerical computation are consistent with previous theoretical conclusions.     

微分本构粘弹性轴向运动弦线横向振动分析的差分法

给出了微分本构粘弹性轴向运动弦线横向振动数值仿真的一种差分法．文中建立了具有微分本构的粘弹性运动弦线的横向振动模型;通过对系统的控制方程和本构方程在不同的分数节点离散,得到一种新的差分方法．利用这一方法,弦线振动方程的数值计算过程可以交替地显式进行,且有较小的截断误差和好的数值稳定性．与通用的方法比较,新的方法计算简单、方便．文中利用方程的不变量检验了数值结果的可靠性,并利用这一方法给出了一类弦线模型的参数振动分析．     

A high‐order ADI finite difference scheme for a 3D reaction‐diffusion equation with neumann boundary condition

In this article, we extend the fourth‐order compact boundary scheme in Liao et al. (Numer Methods Partial Differential Equations 18 (2002), 340–354) to a 3D problem and then combine it with the fourth‐order compact alternating direction implicit (ADI) method in Gu et al. (J Comput Appl Math 155 (2003), 1–17) to solve the 3D reaction‐diffusion equation with Neumann boundary condition. First, the reaction‐diffusion equation is solved with a compact fourth‐order finite difference method based on the Padé approximation, which is then combined with the ADI method and a fourth‐order compact scheme to approximate the Neumann boundary condition, to obtain fourth order accuracy in space. The accuracy in the temporal dimension is improved to fourth order by applying the Richardson extrapolation technique, although the unconditional stability of the numerical method is proved, and several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed new algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

四阶杆振动方程的一族高稳定的十字架格式

用辛几何的观点得到了四阶杆振动方程的一族十字架辛格式，对于四阶杆振动方程的稳定条件不一定随时间方向的精度的提高而放宽，而随空间方向精度的提高稳定范围缩小．数值例子表明单辛算法具有良好的数值稳定性．     

Rationalized Haar functions method for solving Fredholm and Volterra integral equations

This paper presents a computational technique for Fredholm integral equation of the second kind and Volterra integral equation of the second kind. The method is based upon Haar functions approximation. Properties of Rationalized Haar functions are first presented, the operational matrix of integration together with product operational matrix and Newton–Cotes nodes are utilized to reduce the computation of integral equations into some algebraic equations. The method is computationally attractive and applications are demonstrated through illustrative examples.     

Method for solving the Boltzmann kinetic equation for polyatomic gases

A method is proposed for computing the collision operator of a generalized Boltzmann kinetic equation with allowance for energy transfer from translational to vibrational or rotational degrees of freedom. The collision operator is computed using a projection method on a uniform velocity grid. The operator satisfies the mass, momentum, and energy conservation laws and vanishes for an equilibrium velocity distribution function. Approximate models are suggested that provide savings on the computation of rotational-translational relaxation. Numerical examples are presented.