A numerical approach for solving weakly singular partial integro‐differential equations via two‐dimensional‐orthonormal Bernstein polynomials with the convergence analysis

In this paper, we develop an efficient matrix method based on two‐dimensional orthonormal Bernstein polynomials (2D‐OBPs) to provide approximate solution of linear and nonlinear weakly singular partial integro‐differential equations (PIDEs). First, we approximate all functions involved in the considerable problem via 2D‐OBPs. Then, by using the operational matrices of integration, differentiation, and product, the solution of Volterra singular PIDEs is transformed to the solution of a linear or nonlinear system of algebraic equations which can be solved via some suitable numerical methods. With a small number of bases, we can find a reasonable approximate solution. Moreover, we establish some useful theorems for discussing convergence analysis and obtaining an error estimate associated with the proposed method. Finally, we solve some illustrative examples by employing the presented method to show the validity, efficiency, high accuracy, and applicability of the proposed technique.     

A finite difference scheme for nonlinear ultra-parabolic equations

In this paper, our aim is to study a numerical method for an ultraparabolic equation with nonlinear source function. Mathematically, the bibliography on initial–boundary value problems for ultraparabolic equations is not extensive although the problems have many applications related to option pricing, multi-parameter Brownian motion, population dynamics and so forth. In this work, we present the approximate solution by virtue of finite difference scheme and Fourier series. For the nonlinear case, we use an iterative scheme by linear approximation to get the approximate solution and obtain error estimates. A numerical example is given to justify the theoretical analysis.     

An approximate solution for linear boundary-value problems with slowly varying coefficients

In this paper, an approximate closed-form solution for linear boundary-value problems with slowly varying coefficient matrices is obtained. The derivation of the approximate solution is based on the freezing technique, which is commonly used in analyzing the stability of slowly varying initial-value problems as well as solving them. The error between the approximate and the exact solutions is given, and an upper bound on the norm of the error is obtained. This upper bound is proportional to the rate of change of the coefficient matrix of the boundary-value problem. The proposed approximate solution is obtained for a two-point boundary-value problem and is compared to its solution obtained numerically. Good agreement is observed between the approximate and the numerical solutions, when the rate of change of the coefficient matrix is small.     

线性常微分方程的全局误差估计和优化求解方法

线性常微分方程初值问题求解在许多应用中起着重要作用.目前,已存在很多的数值方法和求解器用于计算离散网格点上的近似解,但很少有对全局误差(global error)进行估计和优化的方法.本文首先通过将离散数值解插值成为可微函数用来定义方程的残差;再给出残差与近似解的关系定理并推导出全局误差的上界;然后以最小化残差的二范数为目标将方程求解问题转化为优化求解问题;最后通过分析导出矩阵的结构,提出利用共轭梯度法对其进行求解.之后将该方法应用于滤波电路和汽车悬架系统等实际问题.实验分析表明,本文估计方法对线性常微分方程的初值问题的全局误差具有比较好的估计效果,优化求解方法能够在不增加网格点的情形下求解出线性常微分方程在插值解空间中的全局最优解.     

A new method for solving the nonlinear second-order boundary value differential equations

In this paper we use measure theory to solve a wide range of second-order boundary value ordinary differential equations. First, we transform the problem to a first order system of ordinary differential equations (ODE’s) and then define an optimization problem related to it. The new problem is modified into one consisting of the minimization of a linear functional over a set of Radon measures; the optimal measure is then approximated by a finite combination of atomic measures and the problem converted approximatly to a finite-dimensional linear programming problem. The solution to this problem is used to construct the approximate solution of the original problem. Finally we get the error functionalE (we define in this paper) for the approximate solution of the ODE’s problems.     

A cascadic multigrid algorithm for semilinear elliptic problems

Summary. We propose a cascadic multigrid algorithm for a semilinear elliptic problem. The nonlinear equations arising from linear finite element discretizations are solved by Newton's method. Given an approximate solution on the coarsest grid on each finer grid we perform exactly one Newton step taking the approximate solution from the previous grid as initial guess. The Newton systems are solved iteratively by an appropriate smoothing method. We prove that the algorithm yields an approximate solution within the discretization error on the finest grid provided that the start approximation is sufficiently accurate and that the initial grid size is sufficiently small. Moreover, we show that the method has multigrid complexity. Received February 12, 1998 / Revised version received July 22, 1999 / Published online June 8, 2000     

Superconvergence and asymptotic expansions for linear finite element approximations on criss-cross mesh

In this paper, we discuss the error estimation of the linear finite element solution on criss-cross mesh. Using space orthogonal decomposition techniques, we obtain an asymptotic expansion and superconvergence results of the finite element solution. We first prove that the asymptotic expansion has different forms on the two kinds of nodes and then derive a high accuracy combination formula of the approximate derivatives.     

Superconvergence and asymptotic expansions for linear finite element approximations on crisscross mesh

In this paper, we discuss the error estimation of the linear finite element solution on criss-cross mesh. Using space orthogonal decomposition techniques, we obtain an asymptotic expansion and superconvergence results of the finite element solution. We first prove that the asymptotic expansion has different forms on the two kinds of nodes and then derive a high accuracy combination formula of the approximate derivatives.     

Error estimates for large-scale ill-posed problems

The computation of an approximate solution of linear discrete ill-posed problems with contaminated data is delicate due to the possibility of severe error propagation. Tikhonov regularization seeks to reduce the sensitivity of the computed solution to errors in the data by replacing the given ill-posed problem by a nearby problem, whose solution is less sensitive to perturbation. This regularization method requires that a suitable value of the regularization parameter be chosen. Recently, Brezinski et al. (Numer Algorithms 49, 2008) described new approaches to estimate the error in approximate solutions of linear systems of equations and applied these estimates to determine a suitable value of the regularization parameter in Tikhonov regularization when the approximate solution is computed with the aid of the singular value decomposition. This paper discusses applications of these and related error estimates to the solution of large-scale ill-posed problems when approximate solutions are computed by Tikhonov regularization based on partial Lanczos bidiagonalization of the matrix. The connection between partial Lanczos bidiagonalization and Gauss quadrature is utilized to determine inexpensive bounds for a family of error estimates. In memory of Gene H. Golub. This work was supported by MIUR under the PRIN grant no. 2006017542-003 and by the University of Cagliari.     

Convergence of a Projection-Difference Method for the Approximate Solution of a Parabolic Equation with a Weighted Integral Condition on the Solution

In a Hilbert space, we consider an abstract linear parabolic equation defined on an interval with a nonlocal weighted integral condition imposed on the solution. This problem is solved approximately by a projection-difference method with the use of the implicit Euler method in the time variable. The approximation to the problem in the spatial variables is developed with the finite element method in mind. An estimate of the approximate solution is obtained, the convergence of the approximate solutions to the exact solution is proved, and the error estimates, as well as the orders of the rate of convergence, are established.     

Solving a variational parabolic equation with the periodic condition by a projection-difference method with the Crank–Nicolson scheme in time

A solution to a smoothly solvable linear variational parabolic equation with the periodic condition is sought in a separable Hilbert space by an approximate projection-difference method using an arbitrary finite-dimensional subspace in space variables and the Crank–Nicolson scheme in time. Solvability, uniqueness, and effective error estimates for approximate solutions are proven. We establish the convergence of approximate solutions to a solution as well as the convergence rate sharp in space variables and time.     

An iterative method with error estimators

Iterative methods for the solution of linear systems of equations produce a sequence of approximate solutions. In many applications it is desirable to be able to compute estimates of the norm of the error in the approximate solutions generated and terminate the iterations when the estimates are sufficiently small. This paper presents a new iterative method based on the Lanczos process for the solution of linear systems of equations with a symmetric matrix. The method is designed to allow the computation of estimates of the Euclidean norm of the error in the computed approximate solutions. These estimates are determined by evaluating certain Gauss, anti-Gauss, or Gauss–Radau quadrature rules.     

Scheme of the finite element method with multiplicative separation of the singularity for a spectral boundary value problem for a degenerate differential operator

The paper deals with the numerical solution of a generalized spectral boundary value problem for an elliptic operator with degenerating coefficients. We suggest an approximate method based on the multiplicative separation of the singularity, whereby the eigenfunctions are approximated by piecewise linear functions multiplied by a weight specially chosen depending on the order of degeneration of the coefficients. For this method, we obtain error estimates justifying its optimality.     

Solving deconvolution type problems by wavelet decomposition methods

The paper considers an inverse problem associated with equations of the form Kf = g, where K is a convolution-type operator. The aim is to find a solution f for given function g. We construct approximate solutions by applying a wavelet basis that is well adapted to this problem. For this basis we calculate the elementary solutions that are the approximate preimages of the wavelets. The solution for the inverse problem is then constructed as an appropriate finite linear combination of the elementary solutions. Under certain assumptions we estimate the approximation error and discuss the advantages of the proposed scheme.     

一阶线性隐式微分算子方程的近似解

本文采用投影法讨论了Hilbert空间中一阶线性隐式微分算子方程的Cauchy问题的近似解的收敛性及误差估计.     

A recurrence method for a special class of continuous time linear programming problems

This article studies a numerical solution method for a special class of continuous time linear programming problems denoted by (SP). We will present an efficient method for finding numerical solutions of (SP). The presented method is a discrete approximation algorithm, however, the main work of computing a numerical solution in our method is only to solve finite linear programming problems by using recurrence relations. By our constructive manner, we provide a computational procedure which would yield an error bound introduced by the numerical approximation. We also demonstrate that the searched approximate solutions weakly converge to an optimal solution. Some numerical examples are given to illustrate the provided procedure.     

Roundoff error analysis of algorithms based on Krylov subspace methods

We study the roundoff error propagation in an algorithm which computes the orthonormal basis of a Krylov subspace with Householder orthonormal matrices. Moreover, we analyze special implementations of the classical GMRES algorithm, and of the Full Orthogonalization Method. These techniques approximate the solution of a large sparse linear system of equations on a sequence of Krylov subspaces of small dimension. The roundoff error analyses show upper bounds for the error affecting the computed approximated solutions.This work was carried out with the financial contribution of the Human Capital and Mobility Programme of the European Union grant ERB4050PL921378.     

Numerical Solutions of Matrix Differential Models Using Higher-Order Matrix Splines

This paper deals with the construction of approximate solution of first-order matrix linear differential equations using higher-order matrix splines. An estimation of the approximation error, an algorithm for its implementation and some illustrative examples are included.     

Improved componentwise verified error bounds for least squares problems and underdetermined linear systems

Recently Miyajima presented algorithms to compute componentwise verified error bounds for the solution of full-rank least squares problems and underdetermined linear systems. In this paper we derive simpler and improved componentwise error bounds which are based on equalities for the error of a given approximate solution. Equalities are not improvable, and the expressions are formulated in a way that direct evaluation yields componentwise and rigorous estimates of good quality. The computed bounds are correct in a mathematical sense covering all sources of errors, in particular rounding errors. Numerical results show a gain in accuracy compared to previous results.     

Error analysis of approximate solutions to parametric vector quasiequilibrium problems

In this paper, we give some results on error estimates of approximate solutions to parametric vector quasiequilibrium problems in metric linear spaces. Under some special cases, the error estimates are equivalent to H?lder stability or Lipschitz stability of the set-valued solution map at a given point. An application to variational inequalities is also presented.