Als «banded-matrix» getaggte Fragen

2
Schreiben der Finite-Differenz-Matrix der Poisson-Gleichung mit Neumann-Randbedingungen
Ich bin daran interessiert, die Poisson-Gleichung mit dem Finite-Differenzen-Ansatz zu lösen. Ich möchte besser verstehen, wie man die Matrixgleichung mit Neumann-Randbedingungen schreibt. Würde jemand das Folgende überprüfen, ist es richtig? Die Finite-Differenz-Matrix Die Poisson-Gleichung, ∂2u(x)∂x2=d(x)∂2u(x)∂x2=d(x) \frac{\partial^2u(x)}{\partial x^2} = d(x) kann durch eine Finite-Differenz-Matrix-Gleichung angenähert werden, 1(Δx)2M∙u^=d^1(Δx)2M∙u^=d^ \frac{1}{(\Delta x)^2} \textbf{M}\bullet \hat …
4
Wie kann man Variablen neu anordnen, um eine gebänderte Matrix mit minimaler Bandbreite zu erzeugen?
Ich versuche, eine 2D-Poisson-Gleichung durch endliche Differenzen zu lösen. Dabei erhalte ich eine spärliche Matrix mit nur Variablen in jeder Gleichung. Wenn die Variablen beispielsweise wären, würde die Diskretisierung ergeben:555UUU Ui−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jUi−1,j+Ui+1,j−4Ui,j+Ui,j−1+Ui,j+1=fi,jU_{i-1,j} + U_{i+1,j} -4U_{i,j} + U_{i,j-1} + U_{i,j+1} = f_{i,j} Ich weiß, dass ich dieses System durch eine iterative Methode …
1
Wie löst LAPACK tridiagonale Systeme und warum?
In meinem Projekt muss ich bei jedem Zeitschritt ein paar tridiagonale Matrizen lösen, daher ist es wichtig, einen guten Löser für diese zu haben. Ich habe meine eigene Implementierung gemacht, nur die klassische Art, wie es auf Wikipedia beschrieben wird. Ich habe dann versucht, stattdessen Lapack zu verwenden, und zu …
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.