Springen naar inhoud

[wiskunde] gauss-seidel


  • Log in om te kunnen reageren

#1

dirkwb

    dirkwb


  • >1k berichten
  • 4173 berichten
  • Moderator

Geplaatst op 01 november 2008 - 18:29

[attachment=2637:1.PNG]

De jacobi- iteratie snap ik wel en kan ik met niet veel moeite op mijn GR uitrekenen. Maar bij Gauss-seidel moet je de geupdate component gebruiken: ik snap niet hoe je dat moet uitvoeren. Kan iemand me vertellen hoe dat in elkaar steekt?
Quitters never win and winners never quit.

Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#2

TD

    TD


  • >5k berichten
  • 24052 berichten
  • VIP

Geplaatst op 02 november 2008 - 02:15

Bij Gauss-Seidel gebruik je voor de berekening van de i-de variabele (n-de iteratie) de reeds berekende waarden van de variabelen vůůr i - uiteraard niet die vanaf i want die heb je nog niet!
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

#3

dirkwb

    dirkwb


  • >1k berichten
  • 4173 berichten
  • Moderator

Geplaatst op 02 november 2008 - 11:17

Kan ik dan gewoon in die formule van Gauss-Seidel die L naar links toe halen en oplossen naar u(i+1)?
Quitters never win and winners never quit.

#4

TD

    TD


  • >5k berichten
  • 24052 berichten
  • VIP

Geplaatst op 02 november 2008 - 11:30

Ik snap nu (denk ik) waarom de formule zal kloppen, gisteren was het laat. Er staat inderdaad al ui+1 rechts bij iteratie i+1 maar door de L-matrix die ervoor staat, zullen enkel onbekenden uj met j<k als je uk berekent in iteratie i+1 met een niet-nulle coŽfficiŽnt bijdragen.
Wanneer je deze formule echt implementeert schrijf je de vergelijking uit voor je verschillende variabelen, je zal dan zien dat u1 enkel de oude u1 en u2 gebruikt terwijl je voor u2 al de nieuwe u1 gebruikt, maar de oude u2.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)

#5

dirkwb

    dirkwb


  • >1k berichten
  • 4173 berichten
  • Moderator

Geplaatst op 02 november 2008 - 11:40

Wanneer je deze formule echt implementeert schrijf je de vergelijking uit voor je verschillende variabelen, je zal dan zien dat u1 enkel de oude u1 en u2 gebruikt terwijl je voor u2 al de nieuwe u1 gebruikt, maar de oude u2.

Dus je berekent eerst de oplossing met u=[0 0]T en dan krijg je de waarde van u1 in iteratie 1. Vervolgens stop je die nieuwe u1 met u2 =0 in het algoritme om de (nieuwe) u2 voor iteratie 1 te vinden, toch?
Quitters never win and winners never quit.

#6

TD

    TD


  • >5k berichten
  • 24052 berichten
  • VIP

Geplaatst op 02 november 2008 - 11:53

Inderdaad, maar misschien wordt het eenvoudiger om het per variabele te implementeren.
"Malgré moi, l'infini me tourmente." (Alfred de Musset)





0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures