A SUCCESSIVE LEAST SQUARES METHOD FOR STRUCTURED TOTAL LEAST SQUARES

A new method for Total Least Squares (TLS) problems is presented. It differs from previous approaches and is based on the solution of successive Least Squares problems.The method is quite suitable for Structured TLS (STLS) problems. We study mostly the case of Toeplitz matrices in this paper. The numerical tests illustrate that the method converges to the solution fast for Toeplitz STLS problems. Since the method is designed for general TLS problems, other structured problems can be treated similarly.     

A NOTE ON TWO-DIMENSIONAL SYMPLECTIC SCHEME-SHOOTING METHOD FOR SOLVING THE SCHRDINGER EQUATION

The symplectic scheme-shooting method (SSSM) for solving the two-dimensional time-independent Schrodinger equation is presented. The generalized time-independent Schrodinger equation based on the active and cyclic coordinates is evaluated. The procedure of the separation of the active and cyclic coordinates of A2B type molecule (C2(?) symmetry) is given.     

A GENERALIZED SCALAR AUXILIARY VARIABLE METHOD FOR THE TIME-DEPENDENT GINZBURG-LANDAU EQUATIONS

This paper develops a generalized scalar auxiliary variable(SAV) method for the time-dependent Ginzburg-Landau equations.The backward Euler method is used for discretizing the temporal derivative of the time-dependent Ginzburg-Landau equations.In this method,the system is decoupled and linearized to avoid solving the non-linear equation at each step.The theoretical analysis proves that the generalized SAV method can preserve the maximum bound principle and energy stability,and this is confirmed ...     

AN ADAPTIVE MULTI-SCALE CONJUGATE GRADIENT METHOD FOR DISTRIBUTED PARAMETER ESTIMATION OF 2-D WAVE EQUATION

An adaptive multi-scale conjugate gradient method for distributed parameter estimations (or inverse problems) of wave equation is presented. The identification of the coefficients of wave equations in two dimensions is considered. First, the conjugate gradient method for optimization is adopted to solve the inverse problems. Second, the idea of multi-scale inversion and the necessary conditions that the optimal solution should be the fixed point of multi-scale inversion method is considered. An adaptive multi-scale inversion method for the inoerse problem is developed in conjunction with the conjugate gradient method. Finally, some numerical results are shown to indicate the robustness and effectiveness of our method.     

GL method for solving equations in math-physics and engineering

In this paper, we propose a GL method for solving the ordinary and the partial differential equation in mathematical physics and chemics and engineering. These equations govern the acustic, heat, electromagnetic, elastic, plastic, flow, and quantum etc. macro and micro wave field in time domain and frequency domain. The space domain of the differential equation is infinite domain which includes a finite inhomogeneous domain. The inhomogeneous domain is divided into finite sub domains. We present the solution of the differential equation as an explicit recursive sum of the integrals in the inhomogeneous sub domains. Actualy, we propose an explicit representation of the inhomogeneous parameter nonlinear inversion. The analytical solution of the equation in the infinite homogeneous domain is called as an initial global field. The global field is updated by local scattering field successively subdomaln by subdomain. Once all subdomains are scattered and the updating process is finished in all the sub domains, the solution of the equation is obtained. We call our method as Global and Local field method, in short , GL method. It is different from FEM method, the GL method directly assemble inverse matrix and gets solution. There is no big matrix equation needs to solve in the GL method. There is no needed artificial boundary and no absorption boundary condition for infinite domain in the GL method. We proved several theorems on relationships between the field solution and Green＇s function that is the theoretical base of our GL method. The numerical discretization of the GL method is presented. We proved that the numerical solution of the GL method convergence to the exact solution when the size of the sub domain is going to zero. The error estimation of the GL method for solving wave equation is presented. The simulations show that the GL method is accurate, fast, and stable for solving elliptic, parabolic, and hyperbolic equations. The GL method has advantages and wide applications in the 3D electromagnetic （EM）     

A conservative local discontinuous Galerkin method for the solution of nonlinear Schrdinger equation in two dimensions

In this study, we present a conservative local discontinuous Galerkin(LDG) method for numerically solving the two-dimensional nonlinear Schrdinger(NLS) equation. The NLS equation is rewritten as a firstorder system and then we construct the LDG formulation with appropriate numerical flux. The mass and energy conserving laws for the semi-discrete formulation can be proved based on different choices of numerical fluxes such as the central, alternative and upwind-based flux. We will propose two kinds of time discretization methods for the semi-discrete formulation. One is based on Crank-Nicolson method and can be proved to preserve the discrete mass and energy conservation. The other one is Krylov implicit integration factor(IIF) method which demands much less computational effort. Various numerical experiments are presented to demonstrate the conservation law of mass and energy, the optimal rates of convergence, and the blow-up phenomenon.     

A MIXED FINITE ELEMENT METHOD FOR THE TRANSIENT BEHAVIOR OF A SEMICONDUCTOR DEVICE

The transient behavior of a semiconductor device is described by a system of three quasilinear partial differential equations. One is elliptic in form for the electric potential and the other two are parabolic in form for the conservation of electron and hole concentrations. The electric potential equation is discretized by a mixed finite element method. The electron and hole density equations are treated by a Galerkin method that applies a variant of the method of characteristics to the transport terms. Optimal order convergence analysis in L2 is given for the proposed method.     

关于带W权Drazin逆的表示定理及迭代方法

To study singular linear system, Cline and Greville[8] proposed the concept of W-weighted Drazin inverse for the rectangular matrices,where the properties were also discussed. The computation for the W-weighted Drazin inverse is of much interest, which is mainly divided into two kinds of methods: direct method[2,4,6] and iterative method[3,5,7,9,12,13]. In this paper, we study the iterative method and successive matrix squaring(SMS) method for the W-weighted Drazin inverse and generalize the main results in [12,13].     

The Variable Coefficient ABE-I Method and Its Stability

The ABE-I (Alternating Block Explicit-Implicit) method for diffusion problem is extended to solve the variable coefficient problem and the unconditional stability of the ABE-I method is proved by the energy method.     

A new periodic solution to Jacobi elliptic functions of MKdV equation and BBM equation

Based on the homogeneous balance method,the Jacobi elliptic expansion method and the auxiliary equation method,the first elliptic function equation is used to get a new kind of solutions of nonlinear evolution equations.New exact solutions to the Jacobi elliptic function of MKdV equations and Benjamin-Bona-Mahoney (BBM) equations are obtained with the aid of computer algebraic system Maple.The method is also valid for other (1+1)-dimensional and higher dimensional systems.     

求解离散不适定问题的正则化GMERR方法

迭代极小残差方法是求解大型线性方程组的常用方法, 通常用残差范数控制迭代过程.但对于不适定问题, 即使残差范数下降, 误差范数未必下降. 对大型离散不适定问题,组合广义最小误差(GMERR)方法和截断奇异值分解(TSVD)正则化方法, 并利用广义交叉校验准则(GCV)确定正则化参数,提出了求解大型不适定问题的正则化GMERR方法.数值结果表明, 正则化GMERR方法优于正则化GMRES方法.     

基于非结构自适应网格的复合有限体积法

利用文献[1]中将Lax-Wendroff格式和Lax-Friedrichs格式整体复合作用构成二维无结构网格上的复合型有限体积法,同时利用Delaunay方法,根据流场流动特性变化的梯度值为指示器对网格进行加密和粗化,实现自适应,并将此方法应用到二维浅水波方程的求解上,进行了二维部分溃坝,倾斜水跃的数值实验.结果表明,该方法是一个计算稳定、能适应复杂的求解域、能很好地捕捉激波、且计算速度快的算法.     

Interactive bicriterion solution method and its application to critical path method problems

An interactive solution method is developed for bicriterion mathematical programming (BCMP) problems. The new method, called the dichotomous bicriterion mathematical programming (DBCMP) method, combines Tchebycheff theory and the existing paired comparison method (PCM). The DBCMP method is then compared with the PCM method based on critical path method problems with two conflicting objectives: minimizing the total crashing cost and minimizing the total project completion time. The extension of the DBCMP method to BCMP problems with multiple decision makers is also discussed.     

Dual reciprocity hybrid radial boundary node method for the analysis of Kirchhoff plates

A meshless method of dual reciprocity hybrid radial boundary node method (DHRBNM) for the analysis of arbitrary Kirchhoff plates is presented, which combines the advantageous properties of meshless method, radial point interpolation method (RPIM) and BEM. The solution in present method comprises two parts, i.e., the complementary solution and the particular solution. The complementary solution is solved by hybrid radial boundary node method (HRBNM), in which a three-field interpolation scheme is employed, and the boundary variables are approximated by RPIM, which is applied instead of moving least square (MLS) and obtains the Kronecker’s delta property where the traditional HBNM does not satisfy. The internal variables are interpolated by two groups of symmetric fundamental solutions. Based on those, a hybrid displacement variational principle for Kirchhoff plates is developed, and a meshless method of HRBNM for solving biharmonic problems is obtained, by which the complementary solution can be solved.     

一类凸优化的混合下降算法

邻近点算法(PPA)是一类求解凸优化问题的经典算法, 但往往需要精确求解隐式子问题,于是近似邻近点算法(APPA)在满足一定的近似规则下非精确求解PPA的子问题, 降低了求解难度. 本文利用近似规则的历史信息和随机数扩张预测校正步产生了两个方向, 通过随机数组合两个方向获得了一类凸优化的混合下降算法.在近似规则满足的情况下, 给出了混合下降算法的收敛性证明. 一系列的数值试验表明了混合下降算法的有效性和效率性.     

Stabilization method for continuous approximation of transient convection problem

The aim of this work is to study a new finite element (FE) formulation for the approximation of nonsteady convection equation. Our approximation scheme is based on the Streamline Upwind Petrov Galerkin (SUPG) method for space variable, x, and a modified of the Euler implicit method for time variable, t. The most interest for this scheme lies in its application to resolve by continuous (FE) method the complex of viscoelastic fluid flow obeying an Oldroyd‐B differential model; this constituted our aim motivation and allows us to treat the constitutive law equation, which expresses the relation between the stress tensor and the velocity gradient and includes tensorial transport term. To make the analysis of the method more clear, we first study, in this article this modified method for the advection equation. We point out the stability of this new method and the error estimate of the approximation solution is discussed. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

重根牛顿变换的Julia集

作者分析了重根牛顿变换的Julia集理论,并利用迭代法构造了标准牛顿变换、松弛牛顿变换和重根牛顿变换的Julia集．采用实验数学方法,作者得出如下结论:(1)函数f(z)=zα(zβ-1) 的三种牛顿变换Julia集的中心为原点目具有β倍的旋转对称性; (2)三种牛顿变换Julia集的重根吸引域对α具有敏感的依赖性;(3)由于的零点是松弛牛顿变换的中性或斥性不动点,故松弛牛顿变换的Julia集中不存在单根吸引域;(4)由于∞点不是重根牛顿变换的不动点,故重根牛顿变换的Julia集中多为重根和单根吸引域;(5)重根牛顿法受计算误差影响最小,松弛牛顿法次之, 标准牛顿法最大．     

SOLVING TRUST REGION PROBLEM IN LARGE SCALEOPTIMIZATION

1.IntroductionLetHbeagivennxnsymmetricpositivesemidefinitematrixandcER".Inthispaperweconsiderthefollowingquadraticprogrammingwithasimplequadraticalconstraint1r--XTH. CTX=min2(1)s.tIlxlls5a,wheretheparameteraisprescribed.Thisproblemoccursfrequentlyintrustrationmethodforunconstrainedoptimization[11.Anumberofapproachesforsolving(1)havebeenproposedintheliterature[2,3,6,12--16].OnetechniqueistoapproximateaLagrangemultiplierAbyNewton'smethod.Theapproximationofthisparametermaybequitedelicate,ho…     

奇异杂交边界点方法求解泊松方程

As a boundary-type meshless method,the singular hybrid boundary node method(SHBNM)is based on the modified variational principle and the moving least square(MLS)approximation,so it has the advantages of both boundary element method(BEM)and meshless method.In this paper,the dual reciprocity method(DRM)is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution.The general solution is achieved by means of SHBNM,and the particular solution is approximated by using the radial basis function(RBF).Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain.The postprocess is very simple.Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.