Гауссово исключение и решение приведенной системы - PullRequest
Новая
       3

Гауссово исключение и решение приведенной системы

0 голосов
/ 28 мая 2018

Я пытаюсь реализовать исключение Гаусса с поворотом в C, но решения выходят из строя даже при самом простом вводе.Вот код, который я написал 

#include <stdio.h>
#include <math.h>

#define N 15

void gauss(float a[][N], float *b, int n);
void solution(float a[][N], float *b, float *x, int n);

int main(void)
{
    int n, i, j;

    printf("Enter an integer n<=15: ");
    scanf("%d", &n);

    float a[n][N], b[n], x[n];

    for(i=0; i<n; i++)
        for(j=0; j<n; j++) {
            printf("A[%d][%d] = ", i+1, j+1);
            scanf(" %f", &a[i][j]);
        }

    printf("\n");

    for(i=0; i<n; i++) {
        printf("b[%d] = ", i+1);
        scanf(" %f", b+i);
    }

    gauss(a, b, n);

    solution(a, b, x, n);

    printf("The solution is ");

    for(i=0; i<n; i++)
        printf("%g ", b[i]);

    printf("\n");

    return 0;

}

void gauss(float a[][N], float *b, int n)
{
    int i, j, r, l;

    float pivot, temp, m;

    for(r=0; r<n-1; r++) {
        for(l=i=r+1, pivot=fabs(a[l][r]); i<n; i++)
            if(fabs(a[i][r]) > pivot) {
                l = i;
                pivot = fabs(a[l][r]);
            }

        if(pivot > fabs(a[r][r])) {
            for(j=0; j<n; j++) {
                temp = a[r][j];
                a[r][j] = a[l][j];
                a[l][j] = temp;
            }

            temp = b[r];
            b[r] = b[l];
            b[l] = temp;
        }

        for(i=r+1; i<n; i++) {
            m = a[i][r]/a[r][r];

            for(j=0; j<n; j++)
                a[i][j] = a[i][j] - m*a[r][j];

            b[i] = b[i] - m*b[r];
        }
    }
}

void solutions(float a[][N], float *b, float *x, int n)
{
    int i, j;
    float sum = 0;

    for(i=n-1; i>=0; i--) {
        for(j=i; j<n; j++)
            sum += a[i][j]*x[j];
        x[i] = (b[i] - sum)/a[i][i];
        sum = 0;
    }
}

Я проверил его с помощью тривиальной системы, и получился следующий результат: 

Enter an integer n<=15: 2

A[1][1] = 1
A[1][2] = 1
A[2][1] = 1
A[2][2] = -1

b[1] = 1
b[2] = 0

The solution is 1 -1

Я попытался настроить индексы, но этоне сработало, тогда я записал на бумаге приведенный выше пример, но не смог найти ни одной ошибки.Что не так?

1 Ответ

0 голосов
/ 28 мая 2018

Первая проблема заключается в том, что вы печатаете неправильные значения.Вы печатаете b, но вы сохранили решение в x

for(i=0; i<n; i++)
    printf("%g ", b[i]);
                 ^^^
                  should be x

Кроме того, ваша функция solutions кажется неправильной.

1) Вам следует сбросить sum междуциклы.

2) Вы используете x неинициализированный

Попробуйте что-то вроде: 

void solutions(float a[][N], float *b, float *x, int n)
{
    int i, j;
    float sum;

    for(i=n-1; i>=0; i--) {
      sum = 0;                       // Notice this line
      for(j=i+1; j<n; j++)           // Notice this line, i.e. j=i+1
      {
        sum += a[i][j]*x[j];
      }
      x[i] = (b[i] - sum)/a[i][i];
      sum = 0;
    }
}
...