Fragen zu den algorithmischen / rechnerischen Aspekten der linearen Algebra, einschließlich der Lösung linearer Systeme, Probleme der kleinsten Quadrate, Eigenprobleme und anderer solcher Fragen.
Angenommen, ich habe das ursprüngliche große, spärliche lineare System: . Nun, ich habe nicht A - 1 als ein zu groß , um Faktor oder jede Art von Zersetzung ist A , aber davon ausgehen , dass ich habe die Lösung x 0 mit einem gefunden iterativen lösen.Ax0=b0Ax0=b0A\textbf{x}_0=\textbf{b}_0A−1A−1A^{-1}AAAx0x0\textbf{x}_0 Jetzt möchte …
Sei ein endlichdimensionaler Vektorraum mit der Normund sei F: V \ rightarrow \ mathbb R eine begrenzte lineare Funktion. Es wird nur als Blackbox angegeben.VVV∥⋅∥‖⋅‖\|\cdot\|F:V→RF:V→RF : V \rightarrow \mathbb R Ich möchte die Norm von F schätzen FFF(von oben und unten). Da FFF eine Blackbox ist, besteht die einzige Möglichkeit …
Ich recherchiere über die Struktur in den Schur-Ergänzungen und finde ein interessantes Phänomen: Angenommen, A stammt von 5 - Punkt Laplace. Wenn ich zur Berechnung der LU-Faktorisierung eine verschachtelte Dissektionsreihenfolge und eine multifrontale Methode verwende und dann den letzten Schur-Komplementblock überprüfe, hat er einen niedrigen Rang für die nicht diagonalen …
Ich weiß, dass wir zur Lösung des symmetrischen Eigenwertproblems das Sylvester-Trägheitsgesetz verwenden können, dh die Anzahl der Eigenwerte von A kleiner als a entspricht der Anzahl der negativen Einträge von D, wobei die Diagonalmatrix D von der stammt LDL - Faktorisierung von A - a I = L D L …
Bei einer generischen dünn besetzten Matrix mit m << n (Korrektur: m ≪ n 2 ) Nicht-Null-Elementen (typischerweise m ∈ O ( n ) ). A ist generisch in dem Sinne, dass es keine spezifischen Eigenschaften (z. B. positive Bestimmtheit) aufweist und keine Struktur (z. B. Streifenbildung) angenommen wird.A ∈ …
Es wird gezeigt (Yousef Saad, Iterative Methoden für spärliche lineare Systeme , S. 260), dassc o n d( A.'A ) ≈ c o n d( A )2cond(A′A)≈cond(A)2cond(A'A) \approx cond(A)^2 Gilt das auch für ?A A.'AA′AA' Im Fall ist mit , beachten , dass IEINAAN.× M.N×MN\times MN.≪ M.N≪MN \ll Mc o …
Angenommen , das folgende lineare System gegeben wobei die Laplace - gewichteten ist bekannt , positiv sein definit ein eindimensionaler Nullraum von spannte und die Übersetzungsvarianz von , dh , ändert den Funktionswert nicht (dessen Ableitung ). Die einzigen positiven Einträge von befinden sich in seiner Diagonale, was eine Summe …
Ich habe einige Fragen zu folgenden Themen: Ich versuche, die Schrödinger-Gleichung in 1D unter Verwendung der Kurbel-Nicolson-Diskretisierung zu lösen, gefolgt von der Invertierung der resultierenden tridiagonalen Matrix. Mein Problem hat sich nun zu einem Problem mit periodischen Randbedingungen entwickelt, und deshalb habe ich meinen Code geändert, um den Sherman Morrison-Algorithmus …
Unbestimmte Matrizensysteme treten beispielsweise bei der Diskretisierung von Sattelpunktproblemen durch gemischte finite Elemente auf. Die Systemmatrix kann dann in das Formular eingefügt werden ( A.B.B.tC.)(ABtBC)\begin{pmatrix} A & B^t \\ B & C\end{pmatrix} wobei negativ (semi) -definit ist, positiv (semi) definit ist und willkürlich ist. Natürlich können Sie je nach Konvention …
Bei der Lösung spärlicher linearer Systeme mit direkten Faktorisierungsmethoden wirkt sich die verwendete Ordnungsstrategie erheblich auf den Füllfaktor von Nicht-Null-Elementen in den Faktoren aus. Eine solche Bestellstrategie ist die verschachtelte Dissektion. Ich frage mich, ob es möglich ist, die verschachtelte Dissektionsreihenfolge vorzeitig zu erstellen, wenn nur die Gitterparameter gegeben sind …
Der Titel sagt das meiste davon. Ich suche eine leichte und benutzerfreundliche Bibliothek, die ich für Android (NDK) -Projekte verwenden kann. Für dichtes Material verwende ich gerne Eigen, aber ich habe nicht viele umfassende (und dokumentierte!) Bibliotheken für spärliches Material gefunden, die in einem Projekt "nur funktionieren". PETSc scheint mit …
Ich bin ein Neuling in Eigenwertalgorithmen, aber etwas macht mich auf sich aufmerksam. Der QR-Algorithmus arbeitet mit reellen / komplexen Matrizen und erzeugt reelle / komplexe Eigenwerte. Es können jedoch keine komplexen Eigenwerte aus einer realen Matrix erzeugt werden . Hier ein vereinfachtes Beispiel, geschrieben in Julia und abgeleitet von …
Ich habe die folgende Situation: Ich habe eine Folge von Vektoren und für jeden möchte ich das Produkt berechnen, wobei zu Beginn festgelegt ist. Obwohl es keine Informationen über die Struktur von , hat normalerweise ein bestimmtes Muster, in dem viele Werte wiederholt werden, und ich möchte diese Produkte so …
Ich versuche, ein APSP-Problem (All-Pair Shortest Path) in einem gewichteten Diagramm zu lösen. Dieser Graph ist eigentlich ein 1, 2 oder 3-dimensionales Gitter, und die Gewichte an jeder Kante repräsentieren den Abstand zwischen seinen beiden Eckpunkten. Was ich haben möchte, ist die Entfernung des geodätischen Graphen (kürzester Weg durch den …
Ich schreibe eine kleine numerische Bibliothek für 2x2-, 3x3- und 4x4-Matrizen (real, unsymmetrisch). Viele numerische Analysetexte empfehlen dringend, die Wurzeln des charakteristischen Polynoms nicht zu berechnen, und empfehlen die Verwendung des doppelt verschobenen QR-Algorithmus. Die Größe der Matrizen lässt mich jedoch fragen, ob es einfach ausreichen kann, einfach das charakteristische …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.