Форум программистов, компьютерный форум CyberForum.ru

Решение системы уравнений методом Гаусса. Исправить код - C++

Восстановить пароль Регистрация
 
Katan0073
0 / 0 / 0
Регистрация: 26.02.2014
Сообщений: 9
19.03.2014, 23:43     Решение системы уравнений методом Гаусса. Исправить код #1
Добрый день!
Столкнулся с небольшой проблемой, но не удается от неё избавиться. Программа решает систему уравнений методом гаусса.
Суть проблемы: при вводе конкретных коэффициентов уравнений:
C++
1
2
3
0,3  0,1  2 | 2
1   3   0 |1
0.9  0.3  6 |0
при приведении матрицы к стандартному виду последний коэффициент (a33) не обращается в ноль, хотя должен. Причину нашел, но я не знаю, как с ней справиться. Попытался показать проблему на скрине. Как я понял, в этом последнем случае коэффициент k не округляется до 3, а остается 2,99999976. Почему?
Проблемный участок помечен: "приводим к стандартному виду"
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
#include "stdafx.h"
#include <iostream>
#include <iomanip>
#include <stdio.h>
#include <stdlib.h>
#include <conio.h>
#pragma warning(disable: 4996)
using namespace std;
void vvodsohr(); void chtreshsohr(); void nevyazka();
int n, i, j;
 
FILE *ish, *rez;
 
int main()
{
    while (true){
        cout << "MENU:" << endl << "1. Zapis' matricy v file." << endl << "2. Reshenie matricy." << endl << "3. Nevyazka." << endl << "0. Vyhod." << endl;
        cin >> i;
        switch (i){
        case 1: vvodsohr(); break;
        case 2: chtreshsohr(); break;
        case 3: nevyazka(); break;
        case 0: return 0;
        }
    }
}
 
void vvodsohr(){
 
    ish = fopen("ishodnye_dannye.txt", "wb");
    cout << "Vvedite razmer matricy (nxn): ";
    cin >> n;
    fwrite(&n, sizeof(int), 1, ish);
    float **a = new float*[n], *b = new float[n]; //инициализация динамических массивов
    for (i = 0; i < n; i++) a[i] = new float[n];
    for (i = 0; i<n; i++){//заполнение массива и запись в файл
        cout << "Vvodite elementy i svobodny chlen stroki " << i + 1 << ":";
        for (j = 0; j < n; j++){
            cin >> a[i][j];
            fwrite(&a[i][j], sizeof(float), 1, ish);
        }
        cin >> b[i];
        fwrite(&b[i], sizeof(float), 1, ish);
    }
 
    cout << endl << "File sozdan, vvodimye dannye sohraneny." << endl << endl;
    fclose(ish);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
}
 
void chtreshsohr(){
    ish = fopen("ishodnye_dannye.txt", "rb");
    fread(&n, sizeof(float), 1, ish);
    float **a = new float*[n], *b = new float[n];
    for (i = 0; i < n; i++) a[i] = new float[n];
 
    for (i = 0; i < n; i++){//считывание из файла в массив.
        for (j = 0; j < n; j++)
            fread(&a[i][j], sizeof(float), 1, ish);
        fread(&b[i], sizeof(float), 1, ish);
    }
 
 
 
 
    int f, jmax; float k, max, *p = NULL, l; //приводим к стандартному виду
    for (i = 0; i < n; i++){
        max = a[i][i];
        jmax = i;
        for (j = i; j < n; j++)
        if ((a[j][i])>max) { max = a[j][i]; jmax = j; }
 
        p = a[i];
        a[i] = a[jmax];
        a[jmax] = p;
        l = b[i];
        b[i] = b[jmax];
        b[jmax] = l;
 
        for (f = i + 1; f < n; f++){
            k = a[f][i] / a[i][i];
            b[f] = b[f] - k*b[i];
            for (j = i; j < n; j++)
                a[f][j] = a[f][j] - a[i][j] * k;
        }
    }
 
    cout << "Standartnaya:" << endl;
    for (i = 0; i < n; i++){//вывод приведенной матрицы
        for (j = 0; j < n; j++)cout << setw(9) << a[i][j] << " ";
        cout << "|" << b[i] << endl;
    }
 
    for (i = 0; i < n; i++)
    if (a[i][i] == 0){
        cout << "Sistema ne imeet resheniy" << endl;
        cout << "Press any key to exit...";
        getch();
        exit(1);
    }
 
    float s = 0, *x = new float[n];//находим x
    for (i = 0; i<n; i++) x[i] = 0;
    for (i = n - 1; i >= 0; i--){
        s = b[i];
        for (j = n - 1; j>i; j--) s -= x[j] * a[i][j];
        x[i] = s / a[i][i];
    }
 
    for (i = 0; i < n; i++) cout << "x" << i + 1 << "=" << x[i] << " ";
    cout << endl << endl;
    fclose(ish);
    rez = fopen("rezult.txt", "wb");//запись иксов в файл
    for (i = 0; i < n; i++)
        fwrite(&x[i], sizeof(float), 1, rez);
    fclose(rez);
    for (i = 0; i < n; i++) delete[]a[i];
    delete[]a;
    delete[]b;
    delete[]x;
}
Миниатюры
Решение системы уравнений методом Гаусса. Исправить код  
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
19.03.2014, 23:43     Решение системы уравнений методом Гаусса. Исправить код
Посмотрите здесь:

C++ Системы линейных уравнений методом Гаусса
C++ Решение системы СЛАУ одним из методов Гаусса: исправить код
Найти решение системы линейных уравнений методом Гаусса с выбором главного элемента в столбце C++
Системы линейных уравнений методом Гаусса C++
C++ Решение уравнений методом Гаусса
Решение систем линейных уравнений методом Гаусса C++
C++ Решение системы линейных алгебраических уравнений методом Гаусса
C++ Написать программу для решения системы уравнений методом Гаусса (c++)

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 00:58. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru