Witam wszystkich.
Mam problem z funkcja odwracajaca macierz metoda eliminacji Gaussa. Chodz o to ze gdy chce dopisac do macierzy odwracanej macierz jednostkową zgodnie z algorytmem prosta funkcja oparta na 2 petlach chrzani mi ta macierz wyjsciową i nie wiem dlaczego. Prosze o jakas wskazówke co robie tutaj zle albo jak to mozna by inaczej zrobic. Z góry dziekuje.
[URL=http://img863.imageshack.us/i/beztytuutw.jpg/][IMG]http://img863.imageshack.us/img863/8754/beztytuutw.jpg[/IMG][/URL]
Pytanie 1: czemu wyświetlasz 9 kolumn tablicy MATRIX, a nie 10?
Pytanie 2: czy liczba kolumn w tablicy MATRIX nie jest za mała? Po danych wyjściowych mogę się domyślać, że MATRIX ma tylko 8 kolumn, potrzemne zaś 10.
Witam
Wyswietlam 10 tylko cos wcieło przy robieniu screena
for (i=0;i<n;i++)
for (j=0;j<n;j++)
if (i==j)
MATRIX[i][n+j] = 1;
else
MATRIX[i][n+j] = 0;
fprintf(stream,"\n\nMacierz kontrolna5max \n");
for(i=0;i<n;i++){
fprintf(stream ,"\n");
for(j=0;j<n*2;j++){
fprintf(stream,"%5g ",MATRIX[i][j]);
printf("\n");
}
}
Dane do macierzy wejsciowej wczytywane są z pliku. Problem polega na tym dlaczego przy dodaniu tej macierzy jednostkowej z prawej strony pojawiają sie 0 i 1 w macierzy wyjsciowej z lewej strony. Cały program muszę napisać do soboty wiec bede widzieczny za kazdą wskazówke.
Pokombinowałem z tym co powiedziałes o ilosci kolumn.
malloc((2*n)* sizeof(float));
pomogło dzieki
Wita. Znów prosze o pomoc bo teraz utknałem w innym nmiejscu. Po dodaniu macierzy jednostkowej z prawej strony jest taka o to macierz:
Po Podzieleniu wierszy macierzy przez odpowiednie elementy otrzymalem 1 na przekotnej macierzy głównej i to jest ok, ale potem zgodnie z algorytmem Gaussa- Jordana trzeba wyzerowac pozostale elementy w kolumnach. I tu niby proste ale wynik niezadawalajacy.
for (i=0;i<n;i++)
{
MATRIXTEMP=MATRIX[i][i];
for (j=0;j<n*2;j++)
{
MATRIX[i][j]= (MATRIX[i][j]/ MATRIXTEMP);
}
for (int k=0;k<n;k++)
{
MATRIXTEMP=MATRIX[k][i];
if (k!=i)
{
for (j=0;j<n*2;j++);
MATRIX[k][j]=MATRIX[k][j]-(MATRIXTEMP*MATRIX[i][j]);
}
}
};
fprintf(stream,"\n\nMacierz wyjsciowa \n");
for(i=0;i<n;i++){
fprintf(stream ,"\n");
for(j=0;j<n*2;j++){
fprintf(stream,"%10f ",MATRIX[i][j]);
}
}
Z góry dziekuje za pomoc
Widze ze chyba nikt nie wie jak to zrobic wiec prosze zamknąc temat