Lineare Iterationsverfahren < Lin. Gleich.-systeme < Numerik < Hochschule < Mathe < Vorhilfe
 
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Frage) beantwortet    |    | Datum: |  17:10 Do 06.11.2008 |    | Autor: |  esel |   
	   
	  
 | Aufgabe |   Konvergiert das lineare Iterationsverfahren
 
[mm] x^{(i+1)} [/mm] = [mm] Ax^{(i)}+c, [/mm] i = 0,1,2,...
 
mit
 
[mm] A=\bruch{1}{4}\pmat{ 2 & 0 &0 \\ -4 & -2 & 0 \\ -2& -1 & 0}, [/mm] c = [mm] \vektor{1 \\ 0 \\ -1}
 [/mm] 
für jeden Startwert [mm] x^{(0)}? [/mm]  |  
  
Spektralradius der iterationsmatrix ist doch konvergenzkriterium, heiß: ich muss nur diesen ausrechnen, und wenn dieser <1 ist dann konvergiert das Iterationsverfahren, oder?
 
 
Liebe Grüße
 
Anna
 
 
 
 
 
 
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
 
 
      | 
     
    
   | 
  
 |          | 
 
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Antwort) fertig    |    | Datum: |  08:55 Fr 07.11.2008 |    | Autor: |  fred97 |   
	   
	   Ja.
 
 
Wie groß ist denn der Spektralradius ?
 
 
FRED
 
 
      | 
     
    
   | 
  
 
 |   
|                  | 
  
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Frage) beantwortet    |    | Datum: |  21:15 Fr 07.11.2008 |    | Autor: |  esel |   
	   
	   Beim Berechnen des Spektralradius bin ich jetzt auf ein Problem gestoßen. Den Radius davon erhalte ich ja doch durch Eigenwertberechnung der Iterationsmatrix. Nur ist die gleich die Nullmatrix??
 
 
[mm] x^{i+1}= -(L+D)^{-1}R [/mm] * [mm] x^{i} [/mm] +  [mm] (L+D)^{-1} [/mm] * c
 
[mm] -(L+D)^{-1}R= [/mm] - ( [mm] \pmat{ 0 & 0 & 0\\ -1 & 0 & 0 \\ -\bruch{1}{2} & -\bruch{1}{4} & 0} [/mm] +  [mm] \pmat{ \bruch{1}{2} & 0 & 0 \\ 0 & -\bruch{1}{2} & 0 \\ 0 & 0 & 0} ^{-1})*\pmat{ 0 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 }
 [/mm] 
 
 
      | 
     
    
   | 
  
 
 |   
|                          | 
   
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Antwort) fertig    |    | Datum: |  14:42 Sa 08.11.2008 |    | Autor: |  Denny22 |   
	   
	   Hallo,
 
 
es ist der Spektralradius Deiner Matrix $A$ gemeint und der ist
 
 
[mm] $\rho(A)=\frac{1}{2}<1$
 [/mm] 
 
da $A$ die Eigenwerte [mm] $-\frac{1}{2},0,\frac{1}{2}$ [/mm] besitzt. Damit konvergiert Dein Verfahren für alle Startwerte.
 
 
Gruß Denny
 
 
      | 
     
    
   | 
  
 
 |   
|          | 
 
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Antwort) fertig    |    | Datum: |  14:43 Sa 08.11.2008 |    | Autor: |  Denny22 |   
	   
	  
  
> Konvergiert das lineare Iterationsverfahren
 
>  [mm]x^{(i+1)}[/mm] = [mm]Ax^{(i)}+c,[/mm] i = 0,1,2,...
 
>  mit
 
>  [mm]A=\bruch{1}{4}\pmat{ 2 & 0 &0 \\ -4 & -2 & 0 \\ -2& -1 & 0},[/mm] 
 
> c = [mm]\vektor{1 \\ 0 \\ -1}[/mm]
 
>  für jeden Startwert [mm]x^{(0)}?[/mm]
 
 
Ja.
 
 
>  Spektralradius der iterationsmatrix ist doch 
 
> konvergenzkriterium, heiß: ich muss nur diesen ausrechnen, 
 
> und wenn dieser <1 ist dann konvergiert das 
 
> Iterationsverfahren, oder?
 
 
Das ist korrekt. Deine Iterationsmatrix ist $A$. Berechne den Spektralradius.
 
 
>  
 
> Liebe Grüße
 
>  Anna
 
 
Gruß Denny
 
 
      | 
     
    
   | 
  
 
 |   
|                  | 
  
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Frage) beantwortet    |    | Datum: |  14:50 Sa 08.11.2008 |    | Autor: |  Joan2 |   
	   
	   Aber wozu wurde dann der Vektor c gegeben, wenn ich den zum Berechnen nicht verwenden muss?
 
 
      | 
     
    
   | 
  
 
 |   
|                          | 
   
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Antwort) fertig    |    | Datum: |  15:13 Sa 08.11.2008 |    | Autor: |  Denny22 |   
	   
	   Ich bin zwar auch nicht so fit darin, aber siehe mal zum Beispiel bei Splitting-Verfahren auf der Seite:
 
 
 http://de.wikipedia.org/wiki/Spektralradius
 
 
Die Spektralradiusbedingung bezieht sich nur auf die Matrix vor dem [mm] $x^{(k)}$ [/mm] (bei Dir die Matrix $A$) und nicht auf die Inhomogenität. Du weist mit dem Spektralradius von $A$, dass das Iterationsverfahren
 
 
[mm] $x^{(k+1)}=Ax^{(k)}$
 [/mm] 
 
für jeden Startwert [mm] $x^{(0)}$ [/mm] konvergiert, d.h. es konvergiert gegen einen Vektor $x$ mit
 
 
$x=Ax$
 
 
Wenn Du von diesem Iterationsverfahren weist, dass es für jeden Startwert [mm] $x^{(0)}$ [/mm] konvergiert, dann weist Du auch von dem Iterationsverfahren
 
 
[mm] $x^{(k+1)}=Ax^{(k)}+c$
 [/mm] 
 
dass es für jeden Startwert [mm] $x^{(0)}$ [/mm] konvergiert. Denn wenn [mm] $Ax^{(k)}$ [/mm] gegen $x$ kovergiert, dann konvergiert [mm] $Ax^{(k)}+c$ [/mm] gegen $x+c$. Mache Dir mal darüber Gedanken.
 
 
Gruß Denny
 
 
      | 
     
    
   | 
  
 
 |   
|                                  | 
    
 
   | 
  
 
  
   
    
     
	   | Status: | 
	   		           				(Mitteilung) Reaktion unnötig    |    | Datum: |  15:17 Sa 08.11.2008 |    | Autor: |  Joan2 |   
	   
	   DANKE!!!! ^^ Das war super hilfreich für mich
 
 
 
 
      | 
     
    
   | 
  
 
 |   
  
   |