12.02.2002 - 16:14:40
#0693
ALK LeBigMac
Commander
 *spammmmmer* new |
* Lösung linearer Gleichungssysteme: Gesamtschrittverfahren (JACOBI) */
/* 'A . x = B' */
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <math.h>
int main()
{
int p=5, q=5; /* Maximale Anzahl der Zeilen und Spalten festlegen */
float a[p][q], b[p][q];
int i=1, j=1;
int Za; /* Zeilen von A */
int Zb; /* Zeilen von B */
int Sa; /* Spalten von A */
int Sb; /* Spalten von B */
int n=1, t=0;
float v=0, w=0;
float sumZ[p], sumS[q];
int maxi=0, maxj=0;
float r;
float x[p][q]; /* Lösungsvektor */
char s[200], *z;
printf ( "Geben Sie die erste Matrize ein:\n");
while (fgets(s, 200, stdin) != 0){
j=1;
z = strtok(s, " ");
a[i][j] = atoi(z);
j++;
while((z = strtok(NULL, " ")) != 0){
a[i][j] = atoi(z);
j++;
}
i++;
}
Za=i-1;
Sa=j-1;
i=1;
printf ( "Geben Sie die zweite Matrize ein:\n");
while (fgets(s, 200, stdin) != 0){
j=1;
z = strtok(s, " ");
b[i][j] = atoi(z);
j++;
while((z = strtok(NULL, " ")) != 0){
b[i][j] = atoi(z);
j++;
}
i++;
}
Zb=i-1;
Sb=j-1;
printf("Matrize A:\n");
for (i=1; i<=Za; i++){
for (j=1; j<=Sa; j++){
printf("%4.0f", a[i][j]);
printf(" ");
}
printf("\n");
}
printf("Matrize B:\n");
for (i=1; i<=Zb; i++){
for (j=1; j<=Sb; j++){
printf("%4.0f", b[i][j]);
printf(" ");
}
printf("\n");
}
if ( Za!=Zb || Sb!=1){
printf ( "Fehler bei der Eingabe\n");
return 0;
}
/* Fälle in denen das Verfahren nicht konvergiert */
sumZ[0]=0;
for (i=1; i<=Za; i++){
for (j=1; j<=Sa; j++){
sumZ[i] = sumZ[i]+abs(a[i][j]);
}
/* sumZ = Summe der Werte einer Zeile */
if (sumZ[i] > v){
v=sumZ[i]; /* sumZ[i] = Summe der werte der Zeile i */
maxi=i; /* Zeile in der sumZ maximal ist */
}
}
sumS[0]=0;
for (j=1; j<=Sa; j++){
for (i=1; i<=Za; i++){
sumS[j] = sumS[j]+abs(a[i][j]);
}
if (sumS[j] > w){
w=sumS[j]; /* sumS[j] = Summe der Werte der Spalte j */
maxj=j; /* Spalte in der sumS maximal ist */
}
}
/*printf ("maxi=%d\tmaxj=%d\tsumZ=%f\tsumS=%f\n",maxi,maxj,v,w);*/
if (v >= 2*abs(a[maxi][maxi]) &&
w >= 2*abs(a[maxj][maxj])){
printf ("Konvergenzkriterien sind nicht erfüllt\n");
printf ("%4.0f >= 2*|%4.0f| und %4.0f >= 2*|%4.0f|\n", v, a[maxi][maxi],
w, a[maxj][maxj]);
}
for (i=1; i<=Sa; i++){
x[i][1]=0; /* Initialisierung vom Startvektor */
}
printf("Lösungsvektor\n");
for (n=1; n<=20 ; n++){
for (i=1; i<=Sa; i++){
r=0;
for (j=1 ; j<=Sa; j++){
if (j!=i){
r=r+a[i][j]*x[j][n];
}
}
t=n+1;
x[i][t] = (1/a[i][i]) * (b[i][1]-r);
}
/*r= (x[i][n]-x[i][n-1] );
printf ("%f\n", r);*/
}
for (i=1; i<=Sa; i++)
printf("x[%d][%d]=%f\n", i,n, x[i][n]);
return 0;
}
Nötig hab ich's nicht, aber haben will ich's trotzdem!
[ ∧ ]
|
