Font Size: a A A

Some Investigations On Solving Nonlinear Equation

Posted on:2010-11-28Degree:DoctorType:Dissertation
Country:ChinaCandidate:P P TangFull Text:PDF
GTID:1100360302979888Subject:Computational Mathematics
Abstract/Summary:PDF Full Text Request
This dissertation mainly studies some problems in solving nonlinear equations. The study aims at following aspects which are very important in studyingof nonlinear equation solving. One is construction of maximal order iterationmethod based on some mount of information, one is semi-local behavior of iteration method with some kind of theory frame, the other is global behavior forpolynomials.Information that used in every step is very important when we constructan iteration method. Generally, information is determined by values of function and derivations of function at previous approximate points. For all kinds ofinformation, we choose the most natural one and give a definition of standardinformation N(xnsn,xn-1sn-1,…,xn-lsn-l;f)={f(k)(xj):k=0,…,sj-1,j=n-l,…,n}. Since it is veryhard to achieve information, the standard information used in this paper isN(xns,…,xn-l+1s,xn-ls';f).We know that qd method is a kind of root finding method basedof Pade approximation. For some mount of standard information, we obtain ageneralization of Pade approximation. It is a rational approximation which isvery important in construction of iteration method based on standard information.Given standard information, the most efficient iteration method is maximalorder method. We construct two kinds of iteration methods Hp,s,f and Mp,s,f.They are maximal order iteration methods based on standard information. Further more, they are a natural generalization of Halley's family and M(?)ller method,respectively, with great convergence property. where g(·)=1/f(·) and 0<s'≤s,l≥0.whereThen if xi(i=-l,…,-1,0) are close enough to the root z* of f, xn generated bythe two iteration methods will converge to the root nearest to x0 of f as n→∞.Their order of convergence is the unique real positive root of the polynomial(?).The two iteration methods are both iteration methods with maximal orderbased on standard information N(xns,,…,xn-l+1s,xn-ls';f). By using partial Faadi formula, Bell polynomial and cycle index of symmetric group in combinatoricswe obtain an explicit expression of both iteration method Hp,s,f and Mp,s,f.In this dissertation we also have studied semi-local behavior, which is veryimportant in the research of convergence property of iteration method, of iteration method Hp,s,f. We study the semi-local behavior of iteration method Hp,s,fin operator class Sγ(k)(Ω) of domain B(x0, r). A function f is in the operator class Sγ(k)(Ω)of domain B(x0, r) if k is a positive integer, 0<r<(?) and f satisfieswhereΩ:[0,(?)]→[0,(?)] has continuous monotone increasing derivative of order kin [0, r]The operator class Sγ(k)(Ω) of domain B(x0, r) was first obtained by XinghuaWang. Wang first gave weak condition which is a much weaker condition thanSmale's analytic condition. Furthermore, he constructed a very general conditionoperator class Sγ(k)(Ω) of domain B(x0,r). The universal constant also existsunder this kind of condition, that iswhere r0 satisfiesThe iteration method Hp,s,f will convergent and has an error estimate when thefunction f is in the operator class Sγ(k)(Ω) of domain B(x0,r) and initial condi-tionsβ=‖f'(x0)-1f(x0)‖≤b,xi∈B(x0,r),i=-1,…,-l are satisfied. Theconstant b obtained in our dissertation is the same as that of Halley's family andEuler's family, whose universal constant is obtained by Xinghua Wang. Therefore, the universal constant can be applied to much wider iteration methods.Many iterations, such as Halley's family and Euler's family, are rationalwhen f is a polynomial. In the viewpoint of discrete dynamical system, the mostimportant problem is that whether there exists any generally convergent onepoint stationary algorithm. McMullen given a negative answer to this questionin 1987. He proved that there does not exist any generally convergent onepoint stationary algorithm when the degree of polynomial is greater than 3. Furthermore, he also gave a generally convergent rational iteration method forpolynomials of degree 3. Let p(z) = z3+az + b,is generally convergent for polynomials of degree 3.Therefore, we derive a question that does there exist any rational operatorwhich is generally convergent for polynomials of degree less than or equal to 3and can also be used to as root finding algorithm for general functions. We givea positive answer and construct an algorithm. Let p(z) be a polynomial of degreed with the coefficient of zd-1 vanishes for simplicity. LetWe have proved that Tp is generally convergent for quadratic and cubic polynomials and it is an order 2 algorithm for polynomials of degree d≥4. Forany non-polynomials, d can be chosen arbitrarily. Furthermore, Tp is an order 3algorithm for every function no matter it is a polynomial or not, when we fixedd by 3.The method Tp has the following properties:1. The fixed points of Tp(z) in C are either superattracting or repelling. Thatis, the simple roots of p are superattracting and the extraneous fixed pointsin C are all repelling.2. The roots of p with multiplicities ni≥3 are all repelling fixed points ofTp(z) with multipliers 1+(?). The double roots are not fixed points ofTp(z).3. For d≥4 and the coefficient of zd-2 does not vanish, then∞is an attractingfixed point of Tp(z) with multiplier 1 -(?). Hence the Julia set of Tp(z) isbounded.It is an important implement of McMullen's result and meaningful in dynamics.
Keywords/Search Tags:Halley's family, Müller method, semi-local convergence, generally convergent, qd method
PDF Full Text Request
Related items