阻尼牛顿法迭代中病态问题趋于稳定的渐变(瞬变)过程分析

<正>1引言设映像F:DR~n→R~n,考虑非线性方程组F(x)=0,x∈DR~n,其中F(x)=(f_1(x),f_2(x),…,f_n(x))T,分量f_i(x):R~n→R(i=1,2,…,n)是连续可微实值函数.目前,非线性方程组求解的数值方法有牛顿法、同伦型法、单纯形法与胞腔排除法等[1]~[3]牛顿法是一种非常实用的计算方法,迭代公式如下x=x+p,(2)其中x为前次迭代近似,x为紧接着x后的迭代近似,p=-[F'(x)]~(-1)F(x)为牛顿修正,F'(x)为x处的雅可比矩阵.     

函数因子法——非线性方程组求解中处理奇异问题的一种新方法

§1.引言 非线性方程组 F(x)=0,F:D?R~n→R~n (1.1)嵌入参数t,构成同伦H:[0,T]×D?R~(n 1)→R~n,使得 H(0,x~0)=0,H(T,x)=F(x),(1.2)这里T可以是有限的或 ∞,当T为 ∞时以极限过程代替求值.若 H(t.x)=0(1.3)存在连续解x(t):[0,T]→D,则非线性方程组(1.1)的解x~*=x(T).若(1.3)的解     

高阶稀疏局部非线性方程组的一种拟牛顿方法

§1.引言 当用有限元法或有限差分法分析非线性偏微分方程问题时,必然会导致求解非线性方程组的问题,即求 F(x)=0 (1.1)的解.其中,x=(x_1,x_2,…,Xx_n)~T∈D,D?R~n;F:D→R~n是一个非线性映射.因此,有效地求解非线性方程组(1.1),是分析相应的非线性问题的关键. 不管这些非线性问题是来自流体力学、固体力学,还是其他的物理范畴,它们所对应     

对称矩阵分解因子的直接换元修正法

§1.引言 考虑非线性方程组 F(x)=0, (1)其中F:Ω?R~n→R~n使F′(x)对称.本文给出求解(1)的一种分解修正法,这种方法始于Jacobian F′(x)的初始对称三角分解,然后利用换元技巧直接修正上三角分解因子,进而前代与回代求迭代点.本文分析了分解修正法的运算量,证明了这个算法不用重新启动仍具有局部超线性收敛性和大范围收敛性.此外,这个算法自然保持分解因子的稀疏传递性和修正矩阵的对称传递性,特别当Jacobian正定时,还具有正定传递性.由此本文完成了[1]和[2]无法完成的工作.本算法特别适于大规模带状方程组和最优化问题,数值例子也表明了这一点.     

一种保持对称性、稀疏性的拟Newton法

考虑非线性方程组 (1)F(x)=0其中,F(x)=(f_1(x),…,f_n(x))~T,x∈R~n.当F(x)为梯度算子时,F(x)的Jacobian是对称的.这类问题在实际计算工作中大量存在,比如近些年来研究很多的非线性泛函极小化问题就是如此,因此,人们自然地想到要把对于解非线性方程组很有效的Broyden方法发展到对称情形.1970年Powell提出PSB修正:     

换元Newton型方法的计算效能和最佳换元周期

考虑非线性方程组: F(x)=0, (1.1)其中F:R~n→R~n是二次连续可微函数.一般地说,解方程组(1.1)的拟Newton法较Newton法更为有效.我们可以将拟Newton法解释为逐次在R~n的子空间上构造F′(x)的近似(割线近似)得到的算法.按照这种思想,如果将子空间依次循环取成F′(x)的例     

计算简单分歧点的Weber辅助方程离散Newton法

1.引言 本文考虑单参数有限维非线性方程组G:??R~n×R~1→R~n, G(x,λ)=0 x∈R~n,λ∈R (1.1)的数值求解.方程组(1.1)的解在集合     

Newton-like方法的收敛性及其应用

1 引言我们考虑非线性方程组 F(x)=0, (1．1) 其中F:Rn→Rn是给定的非线性向量函数,并具有如下性质: (1)存在x*使得F(x*)=0; (2)F(x)在x*的邻域内是连续可微的; (3)F′(x*)是非奇异的． Newton法是求(1．1)的数值解的经典算法:     

弱半光滑函数总体极小的广义填充函数法

设F:R~n→R为目标函数,并设F存在极小点。我们的目的是求出x∈R~n使得对所有的x∈R~n有 F(X)≤ F(x). (1.1)即求解F的总体极小. 关于求总体极小问题,到目前为止尚无理论上较为成熟、实际计算中又较为有效的方法.葛人溥在[1]中提出一种求解(1.1)的填充函数法.其基本想法是利用填充函数逐次求     

求解大规模非线性方程组的分层多元谱梯度算法

1 引 言 考虑非线性方程组问题: F(x)=0, x∈Rn (1) 其中,F:Rn→Rn为连续可微的非线性映射.我们讨论大规模情形,并假设F(x)的Jacobian矩阵无法获取,或存储量太大无法承受.     

A novel numerical method and its convergence for nonlinear delay Volterra integro-differential equations

In this paper, we suggest a convergent numerical method for solving nonlinear delay Volterra integro-differential equations. First, we convert the problem into a continuous-time optimization problem and then use a shifted pseudospectral method to discrete the problem. Having solved the last problem, we can achieve the pointwise and continuous approximate solutions for the main delay Volterra integro-differential equations. Here, we analyze the convergence of the method and solve some numerical examples to show the efficiency of the method.     

Spectral gradient projection method for solving nonlinear monotone equations

An algorithm for solving nonlinear monotone equations is proposed, which combines a modified spectral gradient method and projection method. This method is shown to be globally convergent to a solution of the system if the nonlinear equations to be solved is monotone and Lipschitz continuous. An attractive property of the proposed method is that it can be applied to solving nonsmooth equations. We also give some preliminary numerical results to show the efficiency of the proposed method.     

A NUMERICAL EMBEDDING METHOD FOR SOLVING THE NONLINEAR COMPLEMENTARITY PROBLEM(I) ——THEORY

AbstractIn this paper, we extend the numerical embedding method for solving the smooth equations to the nonlinear complementarity problem. By using the nonsmooth theory, we prove the existence and the continuation of the following path for the corresponding homotopy equations. Therefore the basic theory of the numerical embedding method for solving the nonlinear complementarity problem is established. In part II of this paper, we will further study the implementation of the method and give some numerical exapmles.     

Numerical method for the simultaneous determination of the hydraulic properties of unsaturated porous media

This paper presents a numerical algorithm for solving the inverse coefficient problem for nonlinear parabolic equations. This problem arises in simultaneous determination of the hydraulic properties of unsaturated porous media from a simple outflow experiment. The novel feature of the method is that it is not based on output least squares. In this method, the unknown functions are represented as polygons (continuous and piecewise linear functions) every new linear pieces that are determined in each time step by using information based only on previous time intervals. The results of some numerical experiments are displayed.     

A globally convergent BFGS method for nonlinear monotone equations without any merit functions

Since 1965, there has been significant progress in the theoretical study on quasi-Newton methods for solving nonlinear equations, especially in the local convergence analysis. However, the study on global convergence of quasi-Newton methods is relatively fewer, especially for the BFGS method. To ensure global convergence, some merit function such as the squared norm merit function is typically used. In this paper, we propose an algorithm for solving nonlinear monotone equations, which combines the BFGS method and the hyperplane projection method. We also prove that the proposed BFGS method converges globally if the equation is monotone and Lipschitz continuous without differentiability requirement on the equation, which makes it possible to solve some nonsmooth equations. An attractive property of the proposed method is that its global convergence is independent of any merit function.We also report some numerical results to show efficiency of the proposed method.

     

在再生核空间中方程AuBu+Cu＝f的求解方法

$ 1 引言 本文研究下面一类非线性算子方程求解问题 AμBμ Cμ=f, (1.1)其中f,μ∈W(Ω),μ(O)=1,||f ||=1,A,B,C∈(W(Ω)→W(Ω)),(W(Ω)→W(Ω))是W(Ω)到W(Ω)的连续线性算子空间,W(Ω)是定义在Ω域上的(Ω是实数域R的有界域)再生核空间。 本文是在再生核空间上,通过将一维非线性算子方程(1.1)转化为二维线性算子方     

Multigrid method for coupled semilinear elliptic equation

This paper introduces a kind of multigrid finite element method for the coupled semilinear elliptic equations. Instead of the common way of directly solving the coupled semilinear elliptic problems on some fine spaces, the presented method transforms the solution of the coupled semilinear elliptic problem into a series of solutions of the corresponding decoupled linear boundary value problems on the sequence of multilevel finite element spaces and some coupled semilinear elliptic problems on a very low dimensional space. The decoupled linearized boundary value problems can be solved by some multigrid iterations efficiently. The optimal error estimate and optimal computational work are proved theoretically and demonstrated numerically. Moreover, the requirement of bounded second‐order derivatives of the nonlinear term in the existing multigrid method is reduced to a Lipschitz continuous condition in the proposed method.     

A parametrization method for solving nonlinear two-point boundary value problems

A sharper version of the local Hadamard theorem on the solvability of nonlinear equations is proved. Additional parameters are introduced, and a two-parameter family of algorithms for solving nonlinear two-point boundary value problems is proposed. Conditions for the convergence of these algorithms are given in terms of the initial data. Using the right-hand side of the system of differential equations and the boundary conditions, equations are constructed from which initial approximations to the unknown parameters can be found. A criterion is established for the existence of an isolated solution to a nonlinear two-point boundary value problem. This solution is shown to be a continuous function of the data specifying the problem.     

A new analytical technique to solve Fredholm’s integral equations

This paper shows that the homotopy analysis method, the well-known method to solve ODEs and PDEs, can be applied as well as to solve linear and nonlinear integral equations with high accuracy. Comparison of the present method with Adomian decomposition method (ADM), which is well-known in solving integral equations, reveals that the ADM is only special case of the present method. Also, some linear and nonlinear examples are presented to show high efficiency and illustrate the steps of the problem resolution.     

Integrability and beyond

An algebraic approach to solving nonlinear functional equations in the Riemann theta functions is stated. By the inverse scattering method and some general methods of the theory of partial differential equations, the solution of the initial boundary value problem for the nonlinear Schrödinger equation is presented. Bibliography:17 titles.