Lineares Gleichungssystem lösen?

[axelf98]

Zitat von paresy:

könntest du mir vielleicht nen tipp geben ? bzw vllt eine kleine demo?

Kein Problem.. Ich hab aber gesehen, dass diese Parameter gar nicht gebraucht werden.. Wohl aber die Backup-Funktion, die man vorher füllen muss, damit am Ende das System noch mal geprüft werden kann... Unten mal ein Anwendungsbeispiel mit Editfeldern...

Die Unit ist nebenbei gesagt noch nicht perfekt, aber löst so gut wie alles, was lösbar ist und hat keine Probleme wie Gauß!

Delphi-Quellcode:

var LoesungSys : Loesungssystem;
    DieMatrix : TMatrix;
    Loesungen : TLoesung;
....

  setlength(DieMatrix,0,0);
  setlength(DieMatrix,Gleichungen,Unbekannte+2);  // vorbereiten

 // Daten aus Edits holen...
 for i := 0 to Gleichungen-1 do
  for j := 0 to Unbekannte do
    begin
    try
     if Felder[i,j].Text <> '' then
     DieMatrix[i,j] := StrtoFloat(Felder[i,j].Text) ELSE
     DieMatrix[i,j] := 0;
     if j = unbekannte then
     begin
      DieMatrix[i,j] :=  - DieMatrix[i,j];
      if DieMatrix[i,j] = 0 then
      begin
       DieMatrix[i,j] := 1E-1000;  // Rundungsfehler wegen der Null
       genau := false;
      end;
     end;
    except
     MessageDLG('Fehler bei der Eingabe in Zeile '+inttostr(i+1)+', Spalte '+ inttostr(j+1),mtinformation,[mbOK],0);
      Felder[i,j].SetFocus;
     Felder[i,j].SelectAll;
     exit;
    end;
   end;
   self.LoesungSys.backup(DieMatrix,Gleichungen,unbekannte+1);  // WICHTIG ZUM ÜBERPRÜFEN SPÄTER
   setlength(loesungen,unbekannte);
                // Lösen...
   try
  Loesungen := LoesungSys.Systemloesen(DieMatrix,Gleichungen,unbekannte+1,unbekannte,1,true);       // Lösen...
   except
   MessageDLG('Das Gleichungssystem ist nicht lösbar!',mtinformation,[mbOK],0);
   self.StatusBar.Panels[1].Text := 'System nicht lösbar.';
   abort;
   end;

   self.meldungen := self.LoesungSys.Meldungen; // Fehler und Meldungen