Lineares Gleichungssystem

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 15:10:51 GMT

blurry333 — Sun, 08 Jun 2014 15:10:51 GMT

Hallo,

ich braechte einen C++ Code zum loesen von linearen Gleichungssystemen.
Ich hab das zwar im Studium gemacht, aber mein Wissen scheint wie erloschen.
Bevor ich da jetzt tagelang rumprogrammiere, waere es gut wenn es schon ein fertiges STueck Code irgendwo im Netz gaebe. Kennt da jemand eine Seite ?

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 15:12:49 GMT

Nathan — Sun, 08 Jun 2014 15:12:49 GMT

Ernsthaft? Lineare Gleichugnssysteme?

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 15:17:21 GMT

blurry333 — Sun, 08 Jun 2014 15:17:21 GMT

ja wieso ? Dachte auch nicht dass ich das jemals wieder brauche, aber muss da gerade was programmieren und dazu brauch ich es.

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 15:33:47 GMT

blurry333 — Sun, 08 Jun 2014 15:33:47 GMT

Ich stell mir das so vor, ich uebergebe der Funktion eine Koeffizientenmatrix und den Ergebnisvektor und sie spuckt mir die Loesung aus

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 15:47:07 GMT

Mechanics — Sun, 08 Jun 2014 15:47:07 GMT

LAPACK

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 19:56:02 GMT

blurry333 — Sun, 08 Jun 2014 19:56:02 GMT

Mich wunderts dass ich hier im Netz keinen Code dafür finde. Selbst ich habe den Algorithmus früher schon mal programmiert. Aber das war in C#. da muss es doch irgendwas im Netz geben. Wenn ich das jetzt selber mache , dann dauert es ewig.

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 20:16:44 GMT

blurry333 — Sun, 08 Jun 2014 20:16:44 GMT

//  Gauss'sches Eliminationsverfahren mit Zeilenpivotisierung
//    Argumente:
//    double a[N][N+1] erweiterte Koeffizientenmatrix   Read/Write
//    int n            Anzahl der Gleichungen           Read
//    double x[N]      Loesungen                        Write
//  Resultat:
//  int Fehlercode  0 fuer Fehler, 1 fuer Erfolg 

int Gaussalg (double a[][NMAX], int n, double x[]) 
{
  int   i, j;                    // Zeile, Spalte
  int   s;                       // Elimininationsschritt
  int   pzeile;                  // Pivotzeile
  int   fehler = 0;              // Fehlerflag
  double f;                      // Multiplikationsfaktor
  const double Epsilon = 0.01;   // Genauigkeit
  double Maximum;                // Zeilenpivotisierung
  extern FILE *fout;

  s = 0;
  do {             // die einzelnen Eliminationsschritte
    fprintf(fout, "Schritt %2i von %2i\n", s+1, n-1);
    Maximum = fabs(a[s][s]);   // groesstes Element
    pzeile = s ;               // suchen
    for (i = s+1; i < n; i++)
      if (fabs(a[i][s]) > Maximum) {
        Maximum = fabs(a[i][s]) ;
        pzeile = i;
      }
    fehler = (Maximum < Epsilon);
    if (fehler) break;           // nicht loesbar 

    if (pzeile != s)  // falls erforderlich, Zeilen tauschen
    { double h;
      for (j = s ; j <= n; j++) {
        h = a[s][j];
        a[s][j] = a[pzeile][j];
        a[pzeile][j]= h;
      }
    }

    // Elimination --> Nullen in Spalte s ab Zeile s+1
    for (i = s + 1; i < n; i++ ) {
      f = -(a[i][s]/a[s][s]);       // Multiplikationsfaktor
      a[i][s] = 0.0;
      for (j = s+1; j <= n ; j++)   // die einzelnen Spalten
        a[i][j] += f*a[s][j];       // Addition der Zeilen i, s
    }
    #if DEBUG
       MatOut (stdout, a, n, n+1);
    #endif
    s++;
  } while ( s < n-1 ) ;

  if (fehler) 
  {
    fprintf (fout, "gauss: Gleichungssystem nicht loesbar\n");
    return 0; 
  }
  else 
  {
    // Berechnen der Loesungen aus der entstandenen Dreiecksmatrix
    // letzte Zeile
    x[n-1] =  a[n-1][n] / a[n-1][n-1];       
    // restliche Zeilen
    for (i = n-2 ; i >= 0; i-- ) 
    {
      for (j = n-1 ; j > i ; j-- ) 
      {
          a[i][n] -= x[j]*a[i][j];    // rechte Seite berechnen
      } 
      x[i] = a[i][n] / a[i][i];       // Loesung
    }
    return 1;  
  }

}

Reply to Lineares Gleichungssystem on Sun, 08 Jun 2014 20:17:44 GMT

blurry333 — Sun, 08 Jun 2014 20:17:44 GMT

Was haltet ihr von dem hier.

Reply to Lineares Gleichungssystem on Mon, 09 Jun 2014 07:35:56 GMT

Marthog — Mon, 09 Jun 2014 07:35:56 GMT

Und jetzt mach den nochmal in schön. Halte die Variablen lokaler, nutze standard-Algorithmen (z.B. swap, min) und Container statt pointer. Außerdem gehört Ausgabe in Datei und Berechnung in getrennte Funktionen.

Reply to Lineares Gleichungssystem on Mon, 09 Jun 2014 04:23:50 GMT

hustbaer — Mon, 09 Jun 2014 04:23:50 GMT

Marthog schrieb:

Halte die Variablen globaler

Du meinst lokaler

Marthog schrieb:

Und jetzt mach den nochmal in schön.

Das ist nicht sein Code
http://www.zeiner.at/informatik/c/Matrix.html

Reply to Lineares Gleichungssystem on Mon, 09 Jun 2014 07:37:09 GMT

Marthog — Mon, 09 Jun 2014 07:37:09 GMT

hustbaer schrieb:

Du meinst lokaler

Ui. Das ist mir jetzt peinlich. Ich habe es mal editiert, nicht das noch jemand dadurch auf komische Gedanken kommt.