Nummerierung der Unbekannten < DGL < Numerik < Hochschule < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 14:19 Sa 16.05.2009 | | Autor: | jumape |
| Aufgabe | Die 2D Poissongleichung mit Dirichlet Randbedingungen im EInheitsquadrat,
[mm] -\Deltau(x,y)=f(x,y) [/mm] im [mm] \Omega=(0,1)², [/mm] u(x,y)=g(x,y) auf [mm] \Gamma=[0,1]² [/mm] \ [mm] \Omega [/mm] sei auf folgende weise diskretisiert:
[mm] \Delta_hu(x,y)=\bruch{1}{h²}[u(x-h,y)+u(x+h,y)+u(x,y-h)+u(x,y+h)-4u(x,y)]
[/mm]
mit Maschenweite h=1/n in jeder Raumrichtung. Nach Elimination der Randbedingungen ergibt sich ein lineares Gleichungssystem mit N=(n-1)² Unbekannten. Bestimmen Sie die zugehörige Matrix bei schachbrettartiger Nummerierung für ungerades n.
|
Als Tip ist angegeben das Rechengebiet aufzuspalten:
[mm] \Omega_h=\Omega_h^{gg}\cup\Omega_h^{uu}\cup\Omega_h^{gu}\cup\Omega_h^{ug}
[/mm]
wobei uu für die felder steht die bezüglich x und y eine ungerade nummerierung haben ug für die die bezüglich x ungerade und bezüglich y gerade nummeriert sind,...
Die Menge der roten, bzw. schwarzen Felder ist dann:
[mm] \Omega_h^{rot}=\Omega_h^{gg}\cup\Omega_h^{uu} [/mm] und [mm] \Omega_h^{schwarz}=\Omega_h^{gu}\cup\Omega_h^{ug}
[/mm]
Schachbrettaartige Nummerierung heißt:
4 9 5
7 3 8
1 6 2
In dieser Reihenfolge werden die Gitterpunkte durchlaufen, also ist der Vektor u der an die Matrix multipliziert wird für ein Feld mit 9 Punkten:
[mm] u=(u_{11},u_{13},u_{22},u_{31},u_{33},u_{12},u_{21},u_{23},u_{32})
[/mm]
Da die Diskretisierung immer nur auf den gegenwärtigen Punkt und auf Felder zurückgreift die eine andere FArbe haben, und erst alle Felder in der einen und dann die in der anderen Farbe drankommen in u, ist die Matrix eine Blockmatrix mit 4 Blöcken die oben links und unten rechts jeweils die Diagonalmatrix mmit -4 auf der Diagonalen ist und von der ich aber leider noch nicht weiß was oben rechts und unten links passiert. Es sind wohl Bandmatrizen mit Bandbreite 5. Die Diagonale und die beiden außen haben nur einsen die beiden eingeschlossenen sind einmal eins dann dreimal 0... Jedenfalls ist das bei meinem Beispiel mit 16 Feldern so. Ich habe aber keine Ahnung wie ich das auf irgendwelche Kenntnisse fundieren kann und was das für ein System sein soll.
Es wäre nett wenn mir jemand helfen könnte.
Mfg jumape
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 14:20 Mi 20.05.2009 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
