viele lineare Gleichungssysteme lösen

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Schrank

Forum-Century


	Beiträge: 203

	Anmeldedatum: 10.07.08

	Wohnort: Berlin

	Version: ---

Verfasst am: 07.08.2008, 16:03 Titel: viele lineare Gleichungssysteme lösen

Hallo,
ich bin mal wieder am Ende meines Lateins. Ich möchte lineare Gleichunssysteme lösen. Für ein einzelnes System funktioniert das ganze mit x = A \ b ja auch schon wunderbar.
Nur muss ich so ca. 5000 Gleichunssysteme lösen und möchte daher die Verwendung einer Schleife vermeiden.
Nur irgendwie ist mir nicht klar, wie ich Matlab dazu bringe. In welche Form muss ich meine Matrizen bringen, damit matlab tut was ich möchte.

Ich hoffe ich hab mein Problem verständlich machen können.
Danke im voraus, der Schrank.

nschlange

Ehrenmitglied


	Beiträge: 1.318

	Anmeldedatum: 06.09.07

	Wohnort: NRW

	Version: R2007b

Verfasst am: 07.08.2008, 16:29 Titel:

Hi,

hast Du auch 5000 unterschiedliche Matrizen und rechte Seiten?
Ich glaube, dann wäre nicht die Schleife der Flaschenhals.
_________________

Viele Grüße
nschlange

"Chuck Norris ejakuliert fluessigen Stahl!"

Schrank

Themenstarter

Forum-Century


	Beiträge: 203

	Anmeldedatum: 10.07.08

	Wohnort: Berlin

	Version: ---

Verfasst am: 07.08.2008, 18:22 Titel:

Hallo nschlange,
schon das Du dich mal wieder meiner Probleme annimmst. Smile

A und b aus dem Beispiel sind jedes mal unterschiedlich. Glücklicherweise sind es aber jedesmal nur 3 Gleichungen mit 3 Unbekannten. So das ich eigentlich Hoffnung auf eine flinke Methode hatte.

Dachte daran aus A eine 3*3*5000 Matrix zu machen und aus b eine 3*5000 und dann soetwas wie .\ zu benutzen.
Nur das ich ja nicht elementweise sondern 3*3 Matrizenweise gerechnet haben möchte...
naja wenns nicht anders geht, dann bastel ich halt ne Schleife.

Vielleicht denke ich aber auch nur wieder zu kompliziert. Hintergrund der Rechnung ist, das ich ca. 5000Frames einer Messung von Punkten im 3d-Raum von einem globalen zu einem durch Punkte definierten Koordinatensystem transformieren muss. Von Interesse sind dann letztlich nur die Koordianten zweier Punkte im neuen Koordinatensystem (pro Frame).
Dazu berechne ich für jedes Frame den Koordiantenursprung und die Koordiantenachsen und bestimme über die Gleichungssysteme die transformierten Koordinaten der beiden Punkte (genau genommen muss ich also 10000LGS lösen).

Also falls jemand noch eine Idee hat die mir das Warten verkürzen könnte, immer raus damit. Ansonsten Danke fürs lesen. Wink

steve

Ehrenmitglied


	Beiträge: 2.022

	Anmeldedatum: 03.09.07

	Wohnort: Wien

	Version: R2023b

Verfasst am: 08.08.2008, 08:26 Titel:

Moin,

ich hab leider im Moment keine Idee, ob man das über Indizierung lösen könnte, aber

Code:

tic
A = rand(3,3,5000);
b = rand(3,5000);
C = zeros(size(b));
for k=1:5000
C(:,k) = A(:,:,k)\b(:,k);
end
toc

Funktion ohne Link?

ergibt bei mir

Zitat:

Elapsed time is 0.084823 seconds.

Ich denke, dass ist noch im Rahmen der Möglichkeiten, oder?!

Gruß
Alex
_________________

>> I told me to.

____________________________________
Matlab Cheat Sheet
goMatlab-Knigge - dran gehalten?!
Schon in den FAQ gesucht?
Ist vielleicht bei den Skripten oder den Tutorials was für dich dabei?

Schrank

Themenstarter

Forum-Century


	Beiträge: 203

	Anmeldedatum: 10.07.08

	Wohnort: Berlin

	Version: ---

Verfasst am: 08.08.2008, 10:23 Titel:

Hallo steve,
das ist wesentlich schneller als ich erwartet habe! Mit der Zeit kann ich, auch wenn auf meinem Senioren-Pc wahrscheinlich etwas länger, sehr gut leben.

Bleibt nurnoch das ich es unelegant finde, vielleicht setze ich mich ja mal im Urlaub ran und versuche einen hübscheren Weg zu finden. Cool

... wahrscheinlich aber nicht Very Happy

Danke


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.