On optimal iterative methods

Abstract

Let \(\left( X,\rho \right)\) be a complete matrix space, the nonlinear mapping \(\varphi:I\subset X\rightarrow X\) and the equation \(x=\varphi \left(x\right) \) with solution \(x^{\ast}\). We consider another application, \(F:X^{k}\rightarrow X\) for which we assume the diagonal coincides with \(\varphi\):  \(F(x,…,x)=\varphi(x)\). In order to solve the mentioned equation we consider the iterative method \[x_{n+1}=F\left(x_{n},x_{n-1},\cdots,x_{n-k+1}\right),\] \(n=k-1,k,…\) Let \(i_{0},i_{1},….,i_{k-1}\) be a permutation of the numbers \(0,1,…,k-1\) and therefore \(i_{0}-n-k-1,~i_{1}+n-k+1,…,i_{k-1}+n-k+1\) a  permutation of the numbers \(n-k+1,\ n-k+2,…,n\). Among the class of methods given by \[x_{n+1}=F\left( x_{i_{0}+n-k+1},x_{i_{1}+n-k+1},…,x_{i_{k-1}+n-k+1}\right) \] we determine the method for which the difference \(\mathcal{\rho}\left( x^{\ast},x_{n+1}\right)\) is the smallest.

Authors

Ion Păvăloiu
(Tiberiu Popoviciu Institute of Numerical Analysis)

Ioan Şerb
(Tiberiu Popoviciu Institute of Numerical Analysis)

Title

Original title (in French)

Sur des méthodes itératives optimales

English translation of the title

On optimal iterative methods

Keywords

multistep successive approximations; optimal iterative methods

PDF

Cite this paper as:

I. Păvăloiu, I. Şerb, Sur des méthodes itératives optimales, Seminar on functional analysis and numerical methods, Preprint no. 1 (1983), pp. 175-182 (in French).

About this paper

Journal
Seminar on Functional Analysis and Numerical Methods
Preprint
Publisher Name
“Babes-Bolyai” University
Faculty of Mathematics and Physics
Research Seminars
DOI

Not available yet.

References

[1] I. Pavaloiu, Introducere in teoria aproximarii solutiilor ecuatiilor, Ed. Dacia,  Cluj-Napoca, 1976.

[2] I. Pavaloiu, Rezolvarea ecuatiilor prin interpolare, Ed. Dacia, Cluj-Napoca, 1981.

[3] Weinischke, J. H., Uber eine Klasse von Iterationsverfahren, Num. Math., 6 (1964), 395–404.

Paper (preprint) in HTML form

/*====================================================================== Core CSS for LaTeXML documents converted to (X)HTML */ /* Generic Page layout */ .ltx_page_header, .ltx_page_footer { font-size:0.8em; } .ltx_page_header *[rel~="prev"], .ltx_page_footer *[rel~="prev"] { float:left; } .ltx_page_header *[rel~="up"], .ltx_page_footer *[rel~="up"] { display:block; text-align:center; } .ltx_page_header *[rel~="next"], .ltx_page_footer *[rel~="next"] { float:right; } /* What was I trying for here; need more selective rule! .ltx_page_header .ltx_ref, .ltx_page_footer .ltx_ref { margin:0 1em; } */ .ltx_page_header li { padding:0.1em 0.2em 0.1em 1em;} /* Main content */ .ltx_page_content { clear:both; } .ltx_page_header { border-bottom:1px solid; margin-bottom:5px; } .ltx_page_footer { clear:both; border-top:1px solid; margin-top:5px; } .ltx_page_header:after, .ltx_page_footer:after, .ltx_page_content:after { content:"."; display:block; height:0; clear:both; visibility:hidden; } .ltx_page_footer:before { content:"."; display:block; height:0; clear:both; visibility:hidden; } .ltx_page_logo { font-size:80%; margin-top: 5px; clear:both; float:right; } .ltx_page_logo a { font-variant: small-caps; } .ltx_page_logo img { vertical-align:-3px; } /* if shown */ .ltx_page_navbar li { white-space:nowrap; display:block; overflow:hidden; } /* If ref got turned into span, it's "this section"*/ .ltx_page_navbar li span.ltx_ref { white-space:normal; overflow:visible; } /* Ought to be easily removable/overridable? */ .ltx_pagination.ltx_role_newpage { height:2em; } /*====================================================================== Document Structure; Titles & Frontmatter */ /* undo bold here to remove the browser's native h# styling, at let all other styles override it (with more specific rules)*/ .ltx_title { font-size:100%; font-weight:normal; } /* Hack to simulate run-in! put class="ltx_runin" on a title or tag for it to run-into the following text. */ .ltx_runin { display:inline; } .ltx_runin:after { content:" "; } .ltx_runin + .ltx_para, .ltx_runin + .ltx_para p, .ltx_runin + p { display:inline; } .ltx_outdent { margin-left: -2em; } /* .ltx_chapter_title, etc should be in ltx-article.css etc. */ .ltx_page_main { margin:0px; padding:1em 3em 1em 2em; } .ltx_tocentry { list-style-type:none; } /* support for common author block layouts.*/ /* add class ltx_authors_1line to get authors in single line with pop-up affiliation, etc. */ .ltx_authors_1line .ltx_creator, .ltx_authors_1line .ltx_author_before, .ltx_authors_1line .ltx_author_after { display:inline;} .ltx_authors_1line .ltx_author_notes { display:inline-block; } .ltx_authors_1line .ltx_author_notes:before { content:"*"; color:blue;} .ltx_authors_1line .ltx_author_notes span { display:none; } .ltx_authors_1line .ltx_author_notes:hover span { display:block; position:absolute; z-index:10; background:white; text-align:left; border: 1px solid black; border-radius: 0 5px 5px 5px; box-shadow: 5px 5px 10px gray; } /* add class=ltx_authors_multiline to get authors & affliations on separate lines*/ .ltx_authors_multiline .ltx_creator, .ltx_authors_multiline .ltx_author_before, .ltx_authors_multiline .ltx_author_after, .ltx_authors_multiline .ltx_author_notes, .ltx_authors_multiline .ltx_author_notes .ltx_contact { display:block; } /*====================================================================== Para level */ .ltx_float { margin: 1ex 3em 1ex 3em; } td.ltx_subfigure, td.ltx_subtable, td.ltx_subfloat { width:50%; } /* theorems, figure, tables, floats captions.. */ /*====================================================================== Blocks, Lists, Floats */ .ltx_p, .ltx_quote, .ltx_block, .ltx_para { display: block; } /* alignment within blocks */ .ltx_align_left { text-align:left; } .ltx_align_right { text-align:right; } .ltx_align_center { text-align:center; } .ltx_align_justify { text-align:justify; } .ltx_align_top { vertical-align:top; } .ltx_align_bottom { vertical-align:bottom; } .ltx_align_middle { vertical-align:middle; } .ltx_align_baseline { vertical-align:baseline; } .ltx_align_floatleft { float:left; } .ltx_align_floatright { float:right; } .ltx_td.ltx_align_left, .ltx_th.ltx_align_left, .ltx_td.ltx_align_right, .ltx_th.ltx_align_right, .ltx_td.ltx_align_center, .ltx_th.ltx_align_center { white-space:nowrap; } .ltx_td.ltx_align_left.ltx_wrap, .ltx_th.ltx_align_left.ltx_wrap, .ltx_td.ltx_align_right.ltx_wrap, .ltx_th.ltx_align_right.ltx_wrap, .ltx_td.ltx_align_center.ltx_wrap, .ltx_th.ltx_align_center.ltx_wrap, .ltx_td.ltx_align_justify, .ltx_th.ltx_align_justify { white-space:normal; } .ltx_tabular .ltx_tabular { width:100%; } .ltx_inline-block { display:inline-block; } /* equations in non-aligned mode (not normally used) */ .ltx_eqn_div { display:block; width:95%; text-align:center; } /* equations in aligned mode (aligning tags, etc as well as equations) */ .ltx_eqn_table { display:table; width:100%; border-collapse:collapse; } .ltx_eqn_row { display:table-row; } .ltx_eqn_cell { display:table-cell; width:auto; } /* Padding between column pairs in ams align */ table.ltx_eqn_align tr.ltx_equation td.ltx_align_left + td.ltx_align_right, table.ltx_eqn_align tr.ltx_equation td.ltx_align_left + td.ltx_align_center, table.ltx_eqn_align tr.ltx_equation td.ltx_align_center + td.ltx_align_right, table.ltx_eqn_align tr.ltx_equation td.ltx_align_center + td.ltx_align_center { padding-left:3em; } table.ltx_eqn_eqnarray tr.ltx_eqn_lefteqn + tr td.ltx_align_right { min-width:2em; } .ltx_eqn_eqno { max-width:0em; overflow:visible; white-space: nowrap; } .ltx_eqn_eqno.ltx_align_right .ltx_tag { float:right; } .ltx_eqn_center_padleft, .ltx_eqn_center_padright { width:50%; min-width:2em;} .ltx_eqn_left_padleft, .ltx_eqn_right_padright { min-width:2em; } .ltx_eqn_left_padright, .ltx_eqn_right_padleft { width:100%; } /* Various lists */ .ltx_itemize, .ltx_enumerate, .ltx_description { display:block; } .ltx_itemize .ltx_item, .ltx_enumerate .ltx_item { display: list-item; } /* Position the tag to look like a normal item bullet. */ li.ltx_item > .ltx_tag { display:inline-block; margin-left:-1.5em; min-width:1.5em; text-align:right; } .ltx_item .ltx_tag + .ltx_para, .ltx_item .ltx_tag + .ltx_para .ltx_p { display:inline; } /* NOTE: Need to try harder to get runin appearance? */ dl.ltx_description dt { margin-right:0.5em; float:left; font-weight:bold; font-size:95%; } dl.ltx_description dd { margin-left:5em; } dl.ltx_description dl.ltx_description dd { margin-left:3em; } /* Theorems */ .ltx_theorem {margin:1em 0em 1em 0em; } .ltx_title_theorem { font-size:100%; } /* Bibliographies */ .ltx_bibliography dt { margin-right:0.5em; float:left; } .ltx_bibliography dd { margin-left:3em; } /*.ltx_biblist { list-style-type:none; }*/ .ltx_bibitem { list-style-type:none; } .ltx_bibitem .ltx_tag { font-weight:bold; margin-left:-2em; width:3em; } /*.bibitem-tag + div { display:inline; }*/ .ltx_bib_title { font-style:italic; } .ltx_bib_article .bib-title { font-style:normal !important; } .ltx_bib_journal { font-style:italic; } .ltx_bib_volume { font-weight:bold; } /* Indices */ .ltx_indexlist li { list-style-type:none; } .ltx_indexlist { margin-left:1em; padding-left:1em;} /* Listings */ .ltx_listing { display:block; margin: 1ex 3em 1ex 0em; overflow-x:auto; text-align: left; } .ltx_float .ltx_listing { margin: 0; } .ltx_listingline { white-space:nowrap; min-height:1em; } .ltx_lst_numbers_left .ltx_listingline .ltx_tag { background-color:transparent; margin-left:-3em; width:2.5em; position:absolute; text-align:right; } .ltx_lst_numbers_right .ltx_listingline .ltx_tag { background-color:transparent; width:2.5em; position:absolute; right:3em; text-align:right; } /* position:absolute; left:0em; max-width:0em; text-align:right; } */ .ltx_parbox {text-indent:0em; } /* NOTE that it is CRITICAL to put position:relative outside & absolute inside!! I wish I understood why! Outer box establishes resulting size, neutralizes any outer positioning, etc; inner establishes position of stuff to be rotated */ .ltx_transformed_outer { position:relative; bottom:0pt;left:0pt; overflow:visible; } .ltx_transformed_inner { display:block; position:absolute;bottom:0pt;left:0pt; } .ltx_transformed_inner > .ltx_p {text-indent:0em; margin:0; padding:0; } /* If simulating a table (html5), try to get rowspan to work...sorta? */ span.ltx_rowspan { position:absolute; top:0; bottom:0; } /* by default, p doesn't indent */ .ltx_p { text-indent:0em; white-space:normal; } /* explicit control of indentation (on ltx_para) */ .ltx_indent > .ltx_p:first-child { text-indent:2em!important; } .ltx_noindent > .ltx_p:first-child { text-indent:0em!important; } /*====================================================================== Columns */ .ltx_page_column1 { width:44%; float:left; } /* IE uses % of wrong container*/ .ltx_page_column2 { width:44%; float:right; } .ltx_page_columns > .ltx_page_column1 { width:48%; float:left; } .ltx_page_columns > .ltx_page_column2 { width:48%; float:right; } .ltx_page_columns:after { content:"."; display:block; height:0; clear:both; visibility:hidden; } /*====================================================================== Borders and such */ .ltx_tabular { display:inline-table; border-collapse:collapse; } .ltx_tabular.ltx_centering { display:table; } .ltx_thead, .ltx_tfoot, .ltx_tbody { display:table-row-group; } .ltx_tr { display:table-row; } .ltx_td, .ltx_th { display:table-cell; } .ltx_framed { border:1px solid black;} .ltx_tabular .ltx_td, .ltx_tabular .ltx_th { padding:0.1em 0.5em; } /* regular lines */ .ltx_border_t { border-top:1px solid black; } .ltx_border_r { border-right:1px solid black; } .ltx_border_b { border-bottom:1px solid black; } .ltx_border_l { border-left:1px solid black; } /* double lines */ .ltx_border_tt { border-top:3px double black; } .ltx_border_rr { border-right:3px double black; } .ltx_border_bb { border-bottom:3px double black; } .ltx_border_ll { border-left:3px double black; } /* Light lines */ .ltx_border_T { border-top:1px solid gray; } .ltx_border_R { border-right:1px solid gray; } .ltx_border_B { border-bottom:1px solid gray; } .ltx_border_L { border-left:1px solid gray; } /* Framing */ .ltx_framed_rectangle { border-style:solid; border-width:1px; } .ltx_framed_top { border-top-style:solid; border-top-width:1px; } .ltx_framed_left { border-left-style:solid; border-left-width:1px; } .ltx_framed_right { border-right-style:solid; border-right-width:1px; } .ltx_framed_bottom, .ltx_framed_underline { border-bottom-style:solid; border-bottom-width:1px; } .ltx_framed_topbottom { border-top-style:solid; border-top-width:1px; border-bottom-style:solid; border-bottom-width:1px; } .ltx_framed_leftright { border-left-style:solid; border-left-width:1px; border-right-style:solid; border-right-width:1px; } /*====================================================================== Misc */ /* .ltx_verbatim*/ .ltx_verbatim { text-align:left; } /*====================================================================== Meta stuff, footnotes */ .ltx_note_content { display:none; } /*right:5%; */ .ltx_note_content { max-width: 70%; font-size:90%; left:15%; text-align:left; background-color: white; padding: 0.5em 1em 0.5em 1.5em; border: 1px solid black; border-radius: 0 5px 5px 5px; box-shadow: 5px 5px 10px gray; } .ltx_note_mark { color:blue; } .ltx_note_type { font-weight: bold; } .ltx_note { display:inline-block; text-indent:0; } /* So we establish containing block */ .ltx_note_content .ltx_note_mark { position:absolute; left:0.2em; top:-0.1em; } .ltx_note:hover .ltx_note_content, .ltx_note .ltx_note_content:hover { display:block; position:absolute; z-index:10; } .ltx_ERROR { color:red; } .ltx_rdf { display:none; } .ltx_missing { color:red;} .ltx_nounicode { color:red; } /*====================================================================== SVG (pgf/tikz ?) basics */ /* Stuff appearing in svg:foreignObject */ .ltx_svg_fog foreignObject { margin:0; padding:0; overflow:visible; } .ltx_svg_fog foreignObject > p { margin:0; padding:0; display:block; } /*.ltx_svg_fog foreignObject > p { margin:0; padding:0; display:block; white-space:nowrap; }*/ /*====================================================================== Low-level Basics */ /* Note that LaTeX(ML)'s font model doesn't map quite exactly to CSS's */ /* Font Families => font-family */ .ltx_font_serif { font-family: serif; } .ltx_font_sansserif { font-family: sans-serif; } .ltx_font_typewriter { font-family: monospace; } /* dingbats should be converted to unicode? */ /* Math font families handled within math: script, symbol, fraktur, blackboard ? */ /* Font Series => font-weight */ .ltx_font_bold { font-weight: bold; } .ltx_font_medium { font-weight: normal; } /* Font Shapes => font-style or font-variant */ .ltx_font_italic { font-style: italic; font-variant:normal; } .ltx_font_upright { font-style: normal; font-variant:normal; } .ltx_font_slanted { font-style: oblique; font-variant:normal; } .ltx_font_smallcaps { font-variant: small-caps; font-style:normal; } .ltx_font_oldstyle { font-variant: oldstyle-nums; /* experimental css3 ? Doesn't seem to work!*/ font-style:normal; -moz-font-feature-settings: "onum"; -ms-font-feature-settings: "onum"; -webkit-font-feature-settings: "onum"; font-variant-numeric: oldstyle-nums; } .ltx_font_mathcaligraphic { font-family: "Lucida Calligraphy", "Zapf Chancery","URW Chancery L"; } /* .ltx_font_mathscript { ? } */ cite { font-style: normal; } .ltx_red { color:red; } /*.ltx_centering { text-align:center; margin:auto; }*/ /*.ltx_inline-block.ltx_centering,*/ /* Hmm.... is this right in general? */ .ltx_centering { display:block; margin:auto; text-align:center; } /* Dubious stuff */ .ltx_hflipped { display:inline-block; -moz-transform: scaleX(-1); -o-transform: scaleX(-1); -webkit-transform: scaleX(-1); transform: scaleX(-1); filter: FlipH; -ms-fliter: "FlipH"; } .ltx_vflipped { display:inline-block; -moz-transform: scaleY(-1); -o-transform: scaleY(-1); -webkit-transform: scaleY(-1); transform: scaleY(-1); filter: FlipV; -ms-fliter: "FlipV"; } /* .ltx_phantom handled in xslt */ /* A somewhat amsart-like style */ .ltx_title_document { text-transform:uppercase; font-size: 150%; text-align:center; } .ltx_title_section { font-family:serif; font-weight:normal; font-variant:small-caps; font-size: 120%; text-align:center; } .ltx_authors { text-align:center; } .ltx_author_before,.ltx_author_after, .ltx_role_author .ltx_personname { text-transform:uppercase; } .ltx_abstract { margin-bottom:1em; } .ltx_title_abstract { font-family:serif; font-weight:normal; font-variant:small-caps;} .ltx_title_keywords, .ltx_title_classification { display:inline; font-size:100%; font-weight:normal; font-style:italic; } .ltx_role_dedicatory { font-size:100%; font-style:italic; text-align:center; margin:1em;} .ltx_figure .ltx_tag, .ltx_table .ltx_tag { font-family:serif; font-variant:small-caps;} .ltx_title_subsection, .ltx_title_subsubsection { font-family:serif; font-weight:bold;} .ltx_title_subsection:after, .ltx_title_subsubsection:after, .ltx_title_paragraph:after { content:"."; } /* abstract, subsection, subsubsection and paragraph titles are runin */ .ltx_abstract .ltx_title, .ltx_abstract .ltx_title + .ltx_p { display:inline; } .ltx_subsection .ltx_title, .ltx_subsection .ltx_title + .ltx_para, .ltx_subsection .ltx_title + .ltx_para > .ltx_p { display:inline; } .ltx_subsubsection .ltx_title, .ltx_subsubsection .ltx_title + .ltx_para, .ltx_subsubsection .ltx_title + .ltx_para > .ltx_p { display:inline; } .ltx_paragraph .ltx_title, .ltx_paragraph .ltx_title + .ltx_para, .ltx_paragraph .ltx_title + .ltx_para > .ltx_p { display:inline; } .ltx_title +.ltx_para > .ltx_p {text-indent:2em; } .ltx_itemize, .ltx_enumerate, .ltx_description { margin-left: clamp(1em, 3%, 2em); } .ltx_title_document { font-size:170%; text-align:center; margin:0.5em 0 0.5em 0; } .ltx_authors, .ltx_role_author { text-align:center; margin:0.5em 0 0.5em 0; } .ltx_role_author .ltx_personname { font-size: 120%; } .ltx_date { text-align:center; font-size: 120%; margin:0.5em 0 0.5em 0; } .ltx_subtitle { text-align:center; font-size: 120%; padding-left:0.2em; margin-left:-0.5em; } .ltx_title_abstract { text-align:center; font-size: 100%; font-weight:bold; } .ltx_abstract { margin-left:4em; margin-right:4em; } .ltx_title_acknowledgements, .ltx_title_keywords, .ltx_title_classification { text-align:left; font-size: 100%; font-weight:bold; margin:0.5 0 0 0; } .ltx_appendix, .ltx_section, .ltx_subsection, .ltx_subsubsection { margin-top:1.5em; } .ltx_paragraph, .ltx_subparagraph { margin-top:1.0em; } .ltx_title_appendix, .ltx_title_section, .ltx_title_bibliography { font-size:140%; font-weight:bold; margin-bottom:1em; } .ltx_title_subsection { font-size:120%; font-weight:bold; margin-bottom:1em; } .ltx_title_subsubsection { font-size:100%; font-weight:bold; margin-bottom:1em; } /* Paragraph & Subparagraph titles should be runin! */ .ltx_title_paragraph { font-size:100%; font-weight:bold; display:inline; margin-right:1em; } .ltx_paragraph .ltx_title, .ltx_paragraph .ltx_title + .ltx_para, .ltx_paragraph .ltx_title + .ltx_para > .ltx_p { display:inline; } .ltx_title_subparagraph { font-size:100%; font-weight:bold; display:inline; margin-left:2em; margin-right:1em; } .ltx_subparagraph .ltx_title, .ltx_subparagraph .ltx_title + .ltx_para, .ltx_subparagraph .ltx_title + .ltx_para > .ltx_p { display:inline; } .ltx_figure { text-align:center; margin:auto; margin:0.5em; } .ltx_table { text-align:center; margin:auto; margin:0.5em; } /* experimental: flex model, to center short captions, justify long ones */ .ltx_table, .ltx_figure { display:flex; flex-direction:column; align-items:center; } .ltx_table .ltx_caption, .ltx_figure .ltx_caption { text-align:justify; } /* first p in para gets indented */ .ltx_para > .ltx_p:first-child { text-indent:2em; } /* except the initial in a section */ section > .ltx_title +.ltx_para > .ltx_p, section > .ltx_title +.ltx_date +.ltx_para > .ltx_p {text-indent:0em; } .ltx_title_abstract + .ltx_p {text-indent:2em; } .ltx_itemize, .ltx_enumerate { margin-left:1em; } .ltx_theorem { margin-top:0.5em; margin-bottom:0.5em; } .ltx_theorem .ltx_title { margin-bottom:0.1em; } .ltx_theorem .ltx_title + .ltx_para, .ltx_theorem .ltx_title + .ltx_para .ltx_p, .ltx_theorem .ltx_title + .ltx_p { margin-top:0em; }

"Babeş-Bolyai" University

Faculty of Mathematics and Physics

Research Seminars

Seminar on Functional Analysis and Numerical Methods

Preprint Nr.1, 1983, pp. 175-182


On Optimal Iterative Methods

by
Ion Pavaloiu and Ioan Serb

In this work we study iterative methods for solving the equation:

(1) x=f⁢(x),

Orf:I→Xis a given application and(X,r)is a complete metric space.

A l’application fwe will attach an applicationF:Xk→X, which has the following property:

(2) F⁢(x,x,…,x)=f⁢(x),x∈X.

To solve equation ( 1 ) we consider the multi-step iterative method, of the form:

(3) xn+1=F⁢(xn,xn−1,…,xn−k+1),n=k−1,k,k+1,…

Andi0,i1,…,ik−1is a permutation of numbers0,1,…,k−1, SOi0+n−k+1,i1+n−k+1,…,ik−1+n−k+1will be a permutation that matches the numbersn−k+1,n−k+2,…,n.

For such a permutation we obtain a new multi-step iterative method:

(4) xn+1=F⁢(xi0+n−k+1,xi1+n−k+1,…,xik−1+n−k+1),

Orn=k−1,k,…

In this way, from the given applicationF, we obtain k!iterative methods, which are generally distinct.

If under certain conditions imposed on the applicationF, all these k!iterative methods converge to a solution of equation ( 1 ), the problem arises of choosing the method for which the best error bound is obtained.

In the following we will assume thatFverifies the following relationship:

(5) r⁢(F⁢(in1,in2,…,ink),F⁢(in1,in2,…,ink))≤(∑i=1kai⁢r⁢(ini,ini)b)1b,

whatever the elementsin1,in2,…,ink,in1,in2,…,ink∈X, Orai≥0Andb>0are given and0<∑i=1kai<1.

We note that ifFverifies conditions ( 2 ) and ( 5 ) thenfchecks the condition:

r⁢(f⁢(x),f⁢(and)) =r⁢(F⁢(x,x,…,x),F⁢(and,and,…,and))
≤(∑i=1kai⁢r⁢(x,and)b)1b=(∑i=1kai)1b⁢r⁢(x,and),

and the fact that0<∑i=1kai<1, it follows thatfis a contraction, that is, equation ( 1 ) has only one solution.

On the other hand, if we designate byℱb⁢(a1,a2,…,ak)all functionsFverifying conditions ( 2 ) and ( 5 ), then ifb<b′it resultsℱb⁢(a1,a2,…,ak)⊂ℱb′⁢(a1,a2,…⁢ak).

Forb<b′we deduce from Hölder's inequality:

(∑i=1kai⁢r⁢(ini,ini)b)1b =[∑i=1k(aibb′⁢r⁢(ini,ini)b)⁢ai1−bb′]1b
≤(∑i=1kai⁢r⁢(ini,ini)b′)1b′⁢[∑i=1k(ai1−bb′)b′b′−b]b′−bb⁢b′
=(∑i=1kai⁢r⁢(ini,ini)b′)1b′⁢(∑i=1kai)b′−bb⁢b′
<(∑i=1kai⁢r⁢(ini,ini)b′)1b′

SOℱb⁢(a1,a2,…,ak)⊂ℱb′⁢(a1,a2,…,ak).

Let's consider the numbersai∈ℝ, i=1,2,…,k, which verify the relations:

(6) a1≥a2≥…≥ak−1≥ak≥0,

And

(6′) 0<∑i=1kai<1.

We now consider the equations

(7) P⁢(in)=ink−a1⁢ink−1−…−ak−1⁢in−ak=0,
(8) Q⁢(in)=ink−ak⁢ink−1−…−a2⁢in−a1=0,

And

(9) R⁢(in)=ink−ai1⁢ink−1−…−aik−1⁢in−aik=0,

Ori1,i2,…,ikis any permutation of the numbers 1,2,…,k.

LEMME.

Anda1,a2,…,akverify the conditions ( 6 ) and ( 6 ′ ) then each equation of the form ( 9 ) has a single positive root. If we denote by a the positive root of equation ( 7 ) and bybthe positive root of equation ( 8 ), then the positive root cof equation ( 9 ) verifies the following inequalities:

(10) 0<a≤c≤b<1.
Demonstration.

Eitheri1,i2,…,ika permutation of numbers1,2,…,kand eithersthe biggest clue for whichais≠0.

We obviously haveais+1=ais+2=⋯=aik=0Andais≠0.

Consider the functionψ⁢(in)=R⁢(in)/ink−s. SOψ⁢(0)=−ais<0And ψ⁢(1)=R⁢(1)=1−ai1−…−aik−1−aik=1−a1−…−ak>0.

SOψ⁢(in)has a positive root in the interval (0,1). It follows thatR⁢(in)=ink−s⁢ψ⁢(in)has a root in the interval(0,1).

The uniqueness of this root follows immediately by considering the equationf⁢(in)=−ink⁢R⁢(1/in)=0, for whichf′⁢(in)>0andin>0Andf⁢(0)=−1.

To prove inequality ( 10 ) it suffices to show that R⁢(a)≤R⁢(c)=0AndR⁢(b)≥R⁢(c)=0. In fact, we have:

R⁢(a)= ak−ai1⁢ak−1−…−aik−1⁢a−aik
= (a1−ai1)⁢ak−1+(a2−ai2)⁢ak−2+…+ak−aik
= (a−1)⁢(a1−ai1)⁢ak−2+(a1+a2−ai1−ai2)⁢ak−3+…+
+(a1+a2+…+ak−2−ai1−ai2−…−aik−2)⁢a+
+(a1+a2+…+ak−1−ai1−ai2−…−aik−1)≤0,

because0<a<1and from ( 6 ) it follows that

a1+a2+…+as−ai1−ai2−…−ais≥0;

for eachs=1,2,…,k−1.

Similarly, we show thatR⁢(b)≥0,Q.AND.D.

Now we sort the coefficientsaifrom ( 5 ) in descending order, that is:

(11) ai1≥ai2≥…≥aik≥0,

and we write

(11′) as=ais,s=1,2,…,k.

We consider the applicationG:Xk→X, given by the relation

(12) G⁢(in1,in2,…,ink)=F⁢(ini1,ini2,…,inik),

Ori1,i2,…,ikis the permutation of numbers 1,2,…,kwhich corresponds to the relation ( 11 ). We then obtain from ( 5 ), forG, the condition

(13) r⁢(G⁢(in1,in2,…,ink),G⁢(in1,in2,…,ink))≤(∑i=1kai⁢r⁢(ini,ini)b)1b,

whatever the elementsin1,in2,…,ink,in1,in2,…,ink∈X,Or

a1≥a2≥…≥ak≥0,b>0⁢…⁢…⁢.0<∑i=1kai<1.

To solve equation ( 1 ) we consider the following iterative method:

(14) xn+1=G⁢(xn,xn−1,…,xn−k+1),n=k−1,k⁢…,

Orx0,x1,…,xk−1∈X.

If we writeqi=r⁢(xi,xi+1),i=0,1,…, then from ( 13 ) and ( 14 ) we obtain

(15) qi=(∑s=1kas⁢qi−sb)1b,i=k,k+1,….

If in equation ( 7 )a1,a2,…,ak are given by ( 11 ′ ) and ifais the positive root of equation ( 7 ), thena1b∈(0,1).

EitherCa positive constant such that:

qs≤C⁢as/b⁢ For⁢s=0,1,…,k−1.

From ( 15 ) we obtain:

qi ≤(∑s=1kas⁢qi−sb)1b≤(C⁢∑s=1kas⁢ai−s)ib
=C⁢a(i−k)/b⁢(∑s=1kas⁢ak−s)1b=C⁢a(i−k)/b⁢ak/b=C⁢ai/b,

i=k,k+1,…,. Sincea1/b∈(0,1)it follows that the following(xn)n=0∞is fundamental. Let us designate byx¯its limit. It is obvious that x¯is the unique solution of equation ( 1 ) and

(16) r⁢(x¯,xn)≤C⁢an/b1−a1/b.

∎

So, using the lemma we obtain the following theorem:

THEOREM .

Let's suppose thatF∈ℱb⁢(a1,a2,…,ak)and that we give thek!iterative methods of the form ( 4 ), attached toF.

The iterative method for which the best error bound is obtained is that given by ( 14 ), withGdefined by ( 12 ) and which corresponds to the decreasing order of the sequence of coefficientsai,i=1,2,…,k,from ( 5 ).

Remarks .

1) Because, from theF∈ℱb⁢(a1,a2,…,ak)it follows thatfis a contraction, with the Lipschitz constant(∑i=1kai)1/b, it is clear that the iterative method

(17) xn+1=f⁢(xn),x0∈X,n=0,1,…

converges to the solutionx¯from equation ( 1 ). Furthermore, we have

(18) r⁢(x¯,xn)≤C1⁢(∑i=1kai)n/b1−(∑i=1kai)1/b.

We will demonstrate that∑i=1kai≤a, where a is the solution to equation ( 7 ), withas,s=1,2,…,k given by the formula ( 11 ′ ).

On a:

P⁢(∑i=1kai)=
=(∑i=1kai)k−a1⁢(∑i=1kai)k−1−a2⁢(∑i=1kai)k−2−…−ak−1⁢(∑i=1kai)−ak
≤(∑i=1kai)k−ai⁢(∑i=1kai)k−1−a2⁢(∑i=1kai)k−1−…−ak−1⁢(∑i=1kai)k−1−ak⁢(∑i=1kai)k−1
=(∑i=1kai)k−(∑i=1kai)k=0,

from which it follows that∑i=1kai≤a.

Be it nowF∈ℱb⁢(a1,a2,…,ak)Andf⁢(x)=F⁢(x,x,…,x). Formulas ( 16 ) and ( 18 ) valid for any metric and for each pair(F,f),with the above properties, show us that for the iterative method ( 17 ) we obtain a better error delimitation than in the case of other methods of the form ( 4 ). This does not mean that for each choice of the initial points and each applicationF,method ( 17 ) converges faster than all iterative methods of the form ( 4 ).

For example,X=ℝ, r=|⋅|,f(x)=14sinx, F⁢(x,and)=12⁢sin⁡x2⁢cos⁡and2. On a

|F⁢(x1,and1)−F⁢(x2,and2)|≤14⁢|x1−x2|+14⁢|and1−and2|,

for eachx1,x2,and1,and2∈ℝ; F∈ℱ1⁢(14,14)AndF⁢(x,x)=f⁢(x).

We consider the following iterative methods:

(19) xn+1=14⁢sin⁡xn,n=0,1,…,x0=p2,
(20) xn+2=12⁢sin⁡xn+12⁢cos⁡xn2,n=0,1,…,x0=p2,x1=0,
(21) xn+2=12⁢cos⁡xn+12⁢sin⁡xn2,n=0,1,…,x0=p2,x1=0.

In this case the sequence ( 20 ) converges faster than the sequence ( 19 ) or the sequence ( 21 ).

2) There exist iterative methods of the form ( 4 ) respectively ( 17 ), for which the bounds ( 16 ) respectively ( 18 ) cannot improve, i.e.:

limn→∞¯⁢r⁢(x¯,xn)1/n=a1/b,

respectively

limn→∞¯⁢r⁢(x¯,xn)1/n=(∑i=1kai)1/b.
Example .

EitherX=ℝ,r=|⋅|,f(x)=56x,

F⁢(x,and) =12⁢x+13⁢and,
F⁢(x,x) =f⁢(x),
|F⁢(x1,and1)−F⁢(x2,and2)| ≤12⁢|x1−x2|+13⁢|and1−and2|,b=1.

Either

(22) xn+1=56⁢xn,n=0,1,…,x0=1.
(23) xn+2=12⁢xn+1+13⁢xn,n=0,1,…,x0=1,x1=1.
(24) xn+2=13xn+1+12xn,n=0,1,…<x0=1,x1=1.

In this case we obtain the formulas:

(22′) limn→∞(xn)1n=56=12+13,
(23′) limn→∞xn1/n=a,

where a is the positive root of the equationx2−12⁢x−13=0, that's to saya=3+5712. In the same way we obtain

(24′) limn→∞xn1/n=b,

orbis the positive root of the characteristic equation x2−13⁢x−12=0that's to sayb=1+196. It is clear that56<3+5712<1+196.

Bibliography

1983

Related Posts