Static Wikipedia February 2008 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - en - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu

Web Analytics
Cookie Policy Terms and Conditions Satz von Kantorowitsch - Wikipedia

Satz von Kantorowitsch

aus Wikipedia, der freien Enzyklopädie

Der Satz von Kantorowitsch ist eine Aussage der angewandten Mathematik und garantiert die Konvergenz des Newton-Verfahrens unter minimalen Voraussetzungen. Er wurde von Leonid Witaljewitsch Kantorowitsch 1940 erstmals veröffentlicht.

Inhaltsverzeichnis

[Bearbeiten] Voraussetzungen

Es seien X\subset\R^n eine offene Teilmenge und F:\R^n\supset X\to\R^n eine differenzierbare Funktion, deren Ableitung lokal Lipschitz-stetig ist.

D.h. für jedes \mathbf x\in X existiere die Jacobi-Matrix F'(x) der partiellen Ableitungen und es gebe für jede beschränkte Teilmenge U\subset X eine Konstante L > 0 mit

\|F'(\mathbf x)-F'(\mathbf y)\|\le L\;\|\mathbf x-\mathbf y\| für beliebige \mathbf{x,y}\in U.

Die Norm der Differenz der Jacobi-Matrizen ist die induzierte Matrixnorm. Diese in die Vektornorm aufgelöst ergibt die Bedingung

\|F'(\mathbf x)(v)-F'(\mathbf y)(v)\|\le L\;\|\mathbf x-\mathbf y\|\,\|v\|

für beliebige Punkte x,y\in U und Tangentialvektoren v\in\R^n.


In X sei ein Punkt \mathbf x_0\in U bekannt, so dass die Jacobi-Matrix F'(\mathbf x_0) invertierbar ist. Sei h_0=F'(\mathbf x_0)^{-1}F(\mathbf x_0) der Newtonschritt und \mathbf x_1=\mathbf x_0-h_0 das nächste Glied der Newton-Iteration.

Es bezeichne \|h_0\|=\|\mathbf x_1-\mathbf x_0\| die Länge des Newtonschritts.

[Bearbeiten] Aussage

Liegt die Kugel B=\bar K(\mathbf x_1,\|h_0\|) um den Punkt \mathbf x_1 mit der Länge des ersten Newtonschritts als Radius noch vollständig in U und ist die Ungleichung

\alpha_0=M\,\|F'(\mathbf x_0)^{-1}\|\,\|h_0\|\le \tfrac12

erfüllt, wobei M die Lipschitz-Konstante auf B ist, dann

  1. gibt es eine eindeutige Lösung der Vektorgleichung F(\mathbf x)=0 innerhalb der abgeschlossenen Kugel B=\bar K(\mathbf x_1,\|h_0\|) und
  2. konvergiert die Newton-Iteration \mathbf x_{k+1}=\mathbf x_k-F'(\mathbf x_k)^{-1}F(\mathbf x_k) mit Startpunkt \mathbf x_0 mit wenigstens linearer Konvergenzgeschwindigkeit zu dieser Lösung.

[Bearbeiten] Verallgemeinerung

Der normierte Raum (\R^n,\|\cdot\|) kann durch einen beliebigen Banach-Raum in Definitions- und Wertebereich ersetzt werden, die Differenzierbarkeit ist dann durch die Frechet-Ableitung definiert.

Auch im endlichdimensionalen Fall kann man die Normen in Definitions- und Wertebereich unterschiedlich wählen. Mit

\|x\|_{\text{Def}}=\|F'(x_0)(\,x\,)\|_{\text{Wert}}

ergibt sich z.B., dass

\|F'(x_0)^{-1}\|=1

gilt. Die einfachere Form der Konvergenzbedingung ist jedoch aufzuwiegen gegen die komplexere Form der Absätzung zur Lipschitz-Konstanten.


[Bearbeiten] Beweisskizze

Man kann zeigen, dass für ein konvexes Gebiet U mit Lipschitz-Konstante M der ersten Ableitung immer die Ungleichung

\|F(x+h)\|\le \|F(x+h)-F(x)-F'(x)(h)\|+\tfrac12 M\,\|h\|^2

gilt, falls x und x+h in U enthalten sind. Für x0 und x1 = x0 + h0 mit dem Newtonschritt h0 folgt insbesondere

\|F(x_1)\|\le\tfrac12 M\,\|h_0\|^2\le \tfrac12 M\,\|F'(x_0)^{-1}\|\,\|F(x)\|\,\|h_0\|=\tfrac12 \alpha_0 \,\|F(x)\|.

Wegen

\|F'(x_1)-F'(x_0)\|\le M\|h_0\|\le\alpha_0\,\|F'(x_0)^{-1}\|^{-1}

ist F'(x1) nach dem Satz zur Neumann-Reihe ebenfalls invertierbar und es gilt

\|F'(x_1)^{-1}\| \le\tfrac1{1-\alpha_0}\,\|F'(x_0)^{-1}\|\le 2\,\|F'(x_0)^{-1}\|

Diese beiden Abschätzungen kann man zusammenfassen zu einer Abschätzung des nächsten Newtonschrittes h1 = − F'(x1) − 1F(x1):

\|h_1\|\le\alpha_0\,\|h_0\|\le\tfrac12\,\|h_0\|

und der die Konvergenz kontrollierenden Kenngröße

\alpha_1=M\|F'(x_1)^{-1}\|\,\|h_1\|\le2\alpha_0^2\le\tfrac12.

Die Kugel um x2 = x1 + h1 mit Radius \|h_1\| ist vollständig in B und damit in X enthalten, die Lipschitz-Konstante der kleineren Kugel kann nur kleiner sein als M. Es sind also alle Voraussetzungen für den nächsten Schritt hergestellt. Per Induktion wird dies auf die gesamte Newton-Iteration fortgesetzt. Es ergibt sich eine Folge von ineinander enthaltenen Kugeln, deren Radius sich in jedem Schritt mindestens halbiert. Der gemeinsame Durchschnitt aller Kugeln ist also genau ein Punkt, der auch Grenzwert der Newton-Iteration ist. Die Funktionswerte der Newton-Iteration reduzieren sich in jedem Schritt auf ein Viertel des vorhergehenden Funktionswertes, bilden also eine Nullfolge. Der Grenzwert der Newton-Iteration löst also die Vektorgleichung F(x)=0.

[Bearbeiten] Quellen

  • John H. Hubbard und Barbara Burke Hubbard: Vector Calculus, Linear Algebra, and Differential Forms: A Unified Approach,

Matrix Editions, Lesebeispiel der dritten Auflage (vorauss. März 2007)


[Bearbeiten] Literatur

  • Kantorowitsch, L. (1948): Funktionalanalysis und angewandte Mathematik (russ.). UMN3, 6 (28), 89–185.
  • Kantorowitsch, L. W.; Akilow, G. P. (1964): Funktionalanalysis in normierten Räumen. Berlin .
Static Wikipedia 2008 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - en - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu -

Static Wikipedia 2007 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - en - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu -

Static Wikipedia 2006 (no images)

aa - ab - af - ak - als - am - an - ang - ar - arc - as - ast - av - ay - az - ba - bar - bat_smg - bcl - be - be_x_old - bg - bh - bi - bm - bn - bo - bpy - br - bs - bug - bxr - ca - cbk_zam - cdo - ce - ceb - ch - cho - chr - chy - co - cr - crh - cs - csb - cu - cv - cy - da - de - diq - dsb - dv - dz - ee - el - eml - eo - es - et - eu - ext - fa - ff - fi - fiu_vro - fj - fo - fr - frp - fur - fy - ga - gan - gd - gl - glk - gn - got - gu - gv - ha - hak - haw - he - hi - hif - ho - hr - hsb - ht - hu - hy - hz - ia - id - ie - ig - ii - ik - ilo - io - is - it - iu - ja - jbo - jv - ka - kaa - kab - kg - ki - kj - kk - kl - km - kn - ko - kr - ks - ksh - ku - kv - kw - ky - la - lad - lb - lbe - lg - li - lij - lmo - ln - lo - lt - lv - map_bms - mdf - mg - mh - mi - mk - ml - mn - mo - mr - mt - mus - my - myv - mzn - na - nah - nap - nds - nds_nl - ne - new - ng - nl - nn - no - nov - nrm - nv - ny - oc - om - or - os - pa - pag - pam - pap - pdc - pi - pih - pl - pms - ps - pt - qu - quality - rm - rmy - rn - ro - roa_rup - roa_tara - ru - rw - sa - sah - sc - scn - sco - sd - se - sg - sh - si - simple - sk - sl - sm - sn - so - sr - srn - ss - st - stq - su - sv - sw - szl - ta - te - tet - tg - th - ti - tk - tl - tlh - tn - to - tpi - tr - ts - tt - tum - tw - ty - udm - ug - uk - ur - uz - ve - vec - vi - vls - vo - wa - war - wo - wuu - xal - xh - yi - yo - za - zea - zh - zh_classical - zh_min_nan - zh_yue - zu