一类线性隐式A稳定的单步法

1.引言 对于Stiff方程组初值问题的数值解法,Dahlquist在[1]中引进了 A稳定的概念,并且证明了显式的线性多步法(包括显式的Runge-Kutta方法)不可能是A稳定的.现在已经有许许多多隐式A稳定或Stiff稳定的方法,但绝大多数在数值解的过程中必须解由于隐式方法所产生的非线性方程组,而非线性方程组的求解过程往往又要采用Newton-Raphson迭代方法,因此需要计算方程y’=f(x,y)的右函数f(x,y)的Jacobi矩阵以及与此有关的逆矩阵.本文的主要思想是:既然在数值解过程中要计算f(x,y)的Jacobi矩阵,那么不妨在数值公式中明显的出现f(x,y)的一阶偏导数.我们将A稳定公式     

矩阵乘法的推广及运用

利用矩阵的Kronecker积定义了一种矩阵乘积"*积",并且对这种乘积的性质进行了研究,发现它对于任意两个矩阵都有意义而且具有通常矩阵乘积的所有性质,并且在一些特殊情况下它比通常的矩阵乘积更和谐对称,而且当在"合适维数"下它就是通常的矩阵乘积,所以可以把这种"*积"看作是对通常矩阵乘积的推广.     

块半正定矩阵Hadamard积的行列式不等式

基于分块矩阵的Schur补和Albert定理,证明了一些含有块Hadamard积的行列式不等式,并且用不同于文献的方法证明了半正定Hermitian矩阵块Hadamard积的行列式不等式的一个猜想,此结果推广了半正定Hermitian矩阵在块Hadamard积下的Oppenheim不等式.     

非奇异M-矩阵的逆矩阵和M-矩阵的Hadamard积的最小特征值下界估计

给出非奇异M-矩阵的逆矩阵和M-矩阵的Hadamard积的最小特征值下界新的估计式,这些估计式都只依赖于矩阵的元素.数值例子表明,新估计式在一定条件下改进了Fiedler和Markham的猜想,也改进了其它已有的结果.     

矩阵Hadamard积和Fan积的特征值界的一些新估计式（英文）

本文研究了非奇异M-矩阵A与B的Fan积的最小特征值下界和非负矩阵A与B的Hadamard积的谱半径上界的估计问题.利用Brauer定理,得到了一些只依赖于矩阵的元素且易于计算的新估计式,改进了文献[41现有的一些结果.     

r-置换矩阵的一些性质

本文针对正交表和置换群中的置换矩阵问题,提出了r-置换矩阵的概念,研究了其性质,并且给出这类矩阵逆的求法以及利用Hadamard积得出确定方阵为r-置换矩阵的充要条件,对于我们研究和推广置换矩阵有极其重要的意义。     

M-矩阵和逆M-矩阵的Hadamard积(英文)

A=[a_ij]∈M_n和B=[b₍ij(]∈M_n的Hadamard积可表示为AoB=[a_ijb_ij]∈M_n.如果A,B∈M_n是M-矩阵,那么AoB^-1也是M-矩阵.证明了（a）一个非奇异的M-matrix是一对M-矩阵和逆M-矩阵的Hadamard积,同时也证明了（b）一个P-矩阵是两个P-矩阵的Hadamard积.     

广义正定矩阵的Hadamard积和Kronecker积的一些性质

本文讨论了各类型广义正定矩阵的 Hadamard 积和 Kronecker 积的一些重要性质,得到了判断 n 阶实矩阵是广义正定矩阵的一些充要条件,它们是[1]-[4]中相应定理的推广,最后,我们修正了[4]中的一个错误。     

非负矩阵Hadamard积的最大特征值的新上界

关于非负矩阵A和B的Hadamard积的最大特征值的上界问题,主要利用Gerschgorin定理和Brauer定理给出了新的估计式,并把新结果与现有结果进行了比较.数值算例表明新结果在只依赖矩阵元素的条件下改进了现有的一些估计式.     

两类线性隐式方法

众所周知,在Stiff常微分方程组初值问题的数值解法中,向后微分公式(即Gear方法)是目前最通用的方法之一(见[1]).但是,Gear方法是一类隐式方法,在数值解的过程中,一般说来,每向前积分一步,需要解一个非线性方程组,它的求解是采用Newton-Raphson迭代方法,因此需要给出适当精度的预估值和计算右函数f(t,y)的Jacobi阵以     

一类求解单调变分不等式的隐式方法

1．引言变分不等式是一个非常有趣。非常困难的数学问题［"］．它具有广泛的应用（例如，数学规划中的许多基本问题都可以归结为一个变分不等式问题），因而得到深入的研究并有了不少算法［1，2，5－8，17－21］．对线性单调变分不等式，我们最近提出了一系列投影收缩算法Ig－13］．本文考虑求解单调变分不等式其中0CW是一闭凸集，F是从正p到自身的一个单调算子，一即有我们用比（·）表示到0上的投影．求解单调变分不等式的一个简单方法是基本投影法［1，6］，它的迭代式为然而，如果F不是仿射函数，只有当F一致强单调且LIPSChitZ连续…     

Periodic solution for strongly nonlinear vibration systems by He's variational iteration method

In this paper, we use variational iteration method for strongly nonlinear oscillators. This method is a combination of the traditional variational iteration and variational method. Some examples are given to illustrate the effectiveness and convenience of the method. The obtained results are valid for the whole solution domain with high accuracy. The method can be easily extended to other nonlinear oscillations and hence widely applicable in engineering and science. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical methods for nonlinear partial differential equations of fractional order

In this article, we implement relatively new analytical techniques, the variational iteration method and the Adomian decomposition method, for solving nonlinear partial differential equations of fractional order. The fractional derivatives are described in the Caputo sense. The two methods in applied mathematics can be used as alternative methods for obtaining analytic and approximate solutions for different types of fractional differential equations. In these schemes, the solution takes the form of a convergent series with easily computable components. Numerical results show that the two approaches are easy to implement and accurate when applied to partial differential equations of fractional order.     

Regularizing properties of a truncated newton-cg algorithm for nonlinear inverse problems

This paper develops truncated Newton methods as an appropriate tool for nonlinear inverse problems which are ill-posed in the sense of Hadamard. In each Newton step an approximate solution for the linearized problem is computed with the conjugate gradient method as an inner iteration. The conjugate gradient iteration is terminated when the residual has been reduced to a prescribed percentage. Under certain assumptions on the nonlinear operator it is shown that the algorithm converges and is stable if the discrepancy principle is used to terminate the outer iteration. These assumptions are fulfilled, e.g., for the inverse problem of identifying the diffusion coefficient in a parabolic differential equation from distributed data.     

The alternative Legendre Tau method for solving nonlinear multi-order fractional differential equations

In this paper, the alternative Legendre polynomials (ALPs) are used to approximate the solution of a class of nonlinear multi-order fractional differential equations (FDEs). First, the operational matrix of fractional integration of an arbitrary order and the product operational matrix are derived for ALPs. These matrices together with the spectral Tau method are then utilized to reduce the solution of the mentioned equations into the one of solving a system of nonlinear algebraic equations with unknown ALP coefficients of the exact solution. The fractional derivatives are considered in the Caputo sense and the fractional integration is described in the Riemann-Liouville sense. Numerical examples illustrate that the present method is very effective for linear and nonlinear multi-order FDEs and high accuracy solutions can be obtained only using a small number of ALPs.     

On diagonally-preconditioning the 2-step BFGS method with accumulated steps for linearly constrained nonlinear programming

We present an algorithm for very large-scale linearly constrained nonlinear programming (LCNP) based on a Limited-Storage Quasi-newton method. In large-scale programming solving the reduced Newton equation at each iteration can be expensive and may not be justified when far from a local solution; besides, the amount of storage required by the reduced Hessian matrix, and even the computing time for its Quasi-Newton approximation, may be prohibitive. An alternative based on the reduced Truncated-Newton methodology, that has proved to be satisfactory for large-scale problems, is not recommended for very large-scale problems since it requires an additional gradient evaluation and the solving of two systems of linear equations per each minor iteration. We recommend a 2-step BFGS approximation of the inverse of the reduced Hessian matrix that does not require to store any matrix since the product matrix-vector is the vector to be approximated; it uses the reduced gradient and information from two previous iterations and the so-termed restart iteration. A diagonal direct BFGS preconditioning is used.     

Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations

Perturbation methods depend on a small parameter which is difficult to be found for real-life nonlinear problems. To overcome this shortcoming, two new but powerful analytical methods are introduced to solve nonlinear heat transfer problems in this article; one is He's variational iteration method (VIM) and the other is the homotopy-perturbation method (HPM). The VIM is to construct correction functionals using general Lagrange multipliers identified optimally via the variational theory, and the initial approximations can be freely chosen with unknown constants. The HPM deforms a difficult problem into a simple problem which can be easily solved. Nonlinear convective–radiative cooling equation, nonlinear heat equation (porous media equation) and nonlinear heat equation with cubic nonlinearity are used as examples to illustrate the simple solution procedures. Comparison of the applied methods with exact solutions reveals that both methods are tremendously effective.     

三对角逆M矩阵的判定

1、引言 三对角逆M矩阵是指同时为三对角矩阵和逆M矩阵的一类特殊矩阵．文用图论方法探讨三对角逆M矩阵结构,给出了三对角矩阵为逆M矩阵的充分必要条件．此条件提供了判定三对角矩阵是逆M矩阵的方法,但较复杂．文讨论了这类矩阵在Hadamard积下的封闭性．由于三对角逆M矩阵在理论和应用上都有一定价值,所以,寻求一种简单而实用的判定方法是必要的．本文通过对这类矩阵结构特点的研究找到了这样一种方法．同时,由此证明了这类矩阵在Hadamard积下的封闭性．     

Construction of Newton-like iteration methods for solving nonlinear equations

In this paper, we present a simple, and yet powerful and easily applicable scheme in constructing the Newton-like iteration formulae for the computation of the solutions of nonlinear equations. The new scheme is based on the homotopy analysis method applied to equations in general form equivalent to the nonlinear equations. It provides a tool to develop new Newton-like iteration methods or to improve the existing iteration methods which contains the well-known Newton iteration formula in logic; those all improve the Newton method. The orders of convergence and corresponding error equations of the obtained iteration formulae are derived analytically or with the help of Maple. Some numerical tests are given to support the theory developed in this paper.     

Some iterative methods for the largest positive definite solution to a class of nonlinear matrix equation

In this paper, we propose some inversion-free iteration methods for finding the largest positive definite solution of a class of nonlinear matrix equation. Then, we consider the properties of the solution for this nonlinear matrix equation. Also, we establish Newton’s iteration method for finding the largest positive definite solution and prove its quadratic convergence. Furthermore, we derive the semi-local convergence of the Newton’s iteration method. Finally, some numerical examples are presented to illustrate the effectiveness of the theoretical results and the behavior of the considered methods.