1 / 1 / 2
Регистрация: 07.07.2014
Сообщений: 39
1

Решение СЛАУ методом Гаусса

14.04.2016, 09:53. Показов 409. Ответов 0
Метки нет (Все метки)

Доброго времени суток! Помогите, пожалуйста исправить код

Нужно реализовать решение СЛАУ методом Гаусса, коэффициенты матрицы заданы рациональными дробями
программа приводит матрицу к треугольному виду, а вот сами иксы считает неправильно (выводит х=0/0)


C++ (Qt)
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
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
#include <iostream>
#include <fstream>
#include <cmath>
using namespace std;
 
void input(int **nom, int **denom, int m, int n, int*Bnom, int*Bdenom, int**cnom, int**cdenom)
    {
        int i,j;  char ch;
        {
            cout << "\nBведите коефициенты при неизвестных СЛАУ: \n";
            for ( i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    cin >> nom[i][j]; //числитель
                    cin >> ch;  
                    cin >> denom[i][j];  //знаменатель
                }
            }
        }
//----------------------------------------------------------------------------------------------------------------------------------------
        for (i = 0; i < m; i++)
        {
            for (int j = 0; j < n; j++)
            {
                cnom[i][j] = nom[i][j]; 
                cdenom[i][j] = denom[i][j]; 
            }
        }
//----------------------------------------------------------------------------------------------------------------------------------------
        cout << ("Введите столбец свободных членов: \n");
        for (i = 0; i < m; i++)
        {
            cin >> Bnom[i];
            cin >> ch;
            cin >> Bdenom[i];
        }
//----------------------------------------------------------------------------------------------------------------------------------------
        cout << "Cистема уравнений: \n\n";
        for (int i = 0; i<m; i++)
        {
            for (int j = 0; j<n; j++)
            {
                if (nom[i][j]<0 || denom[i][j]<0) cout << "(" << nom[i][j] << ch << denom[i][j] << ")" << "x" << j + 1;
                else cout << nom[i][j] << ch << denom[i][j] << "x" << j + 1 << " ";
                if (j != (n - 1)) cout << " + ";
            }
            cout << " = " << Bnom[i] << ch<< Bdenom[i] << "\n";
        }
//----------------------------------------------------------------------------------------------------------------------------------------
        cout << "Mатрица системы: \n" << endl;
        for ( i = 0; i < m; i++)
        {
            for (j = 0; j < n; j++)
                cout << nom[i][j] << ch << denom[i][j] << "\t";
                cout << endl;
        }
//----------------------------------------------------------------------------------------------------------------------------------------
        cout << "\n Столбец свободных членов: \n";
        for (i = 0; i<m; i++)
        {
            cout  << Bnom[i] << ch << Bdenom[i] << "\t";
            cout << "\n";
        }
    }
 
 
void Gauss(int**nom, int**denom, int**cnom, int**cdenom, int *Bnom, int *Bdenom, int*X, int*X2, int n)
{
    int i, j, k;
    int kf, kf2;
 
    for (int i = 0; i<n; i++)
        for (int j = 0; j < n; j++)
        {
            nom[i][j] = cnom[i][j];
            denom[i][j] = cdenom[i][j];
        }
//----------------------------------------------------------------------------------------------------------------------------------------
    for (i = 0; i<n - 1; i++)
        for (k = i + 1; k<n; k++)
        {
            kf = nom[k][i] * denom[i][i];
            kf2 = denom[k][i] * nom[i][i];
 
            nom[k][i] = 0;
            denom[k][i] = 0;
 
            for (j = i + 1; j<n; j++)
            {
                nom[k][j] = nom[k][j] * kf2*denom[i][j] - kf*nom[i][j] * denom[k][j];
                denom[k][j] = denom[k][j] * kf2*denom[i][j];
            }
            Bnom[k] = Bnom[k] * kf2*Bdenom[i] - kf*Bnom[i];
            Bdenom[k] = Bdenom[k] * kf2*Bdenom[i];
        }
//----------------------------------------------------------------------------------------------------------------------------------------
    cout << "\nМатрица преобразованная до треугольной: \n";
    for (i = 0; i<n; i++)
    {
        for (j = 0; j<n; j++)
            cout << nom[i][j] << '/' << denom[i][j] << '\t';
        cout << endl;
    }
//----------------------------------------------------------------------------------------------------------------------------------------
    cout << "\nПреобразованный столбец свободных членов: \n";
    for (i = 0; i<n; i++)
    {
        cout << Bnom[i] << '/' << Bdenom[i] << '\t';
        cout << endl;
    }
//----------------------------------------------------------------------------------------------------------------------------------------
    for (int m = n - 1; m >= 0; m--) //ошибка где-то в этом цикле
    {
        int sum = 0;
        int sum2 = 0;
 
        for (int i = n - 1; i > m; i--) 
        {
            sum += X[i] * nom[m][i];
            sum2 += X2[i] * denom[m][i];
        }
        X[m] = ((Bnom[m] * sum2) - (sum*Bdenom[m]))*denom[m][m];
        X2[m] = Bdenom[m] * sum2*nom[m][m];
    }
 
    for (int i = 0; i < n; i++)
        cout << "x" << i + 1 << " = " << X[i] << '/' << X2[i] << "\n";
 
}
 
int main()
{
    setlocale(LC_ALL, "Russian");
    int n, m;
    cout << "Введите количество уравнений СЛАУ:\n";
    cin >> m;
    cout << "\nВведите количество неизвестных СЛАУ:\n";
    cin >> n;
 
    int**nom = new int*[m];
    for (int i = 0; i<m; i++)
        nom[i] = new int[n];
 
    int**denom = new int*[m];
    for (int i = 0; i<m; i++)
        denom[i] = new int[n];
 
    int**cnom = new int*[m];
    for (int i = 0; i<m; i++)
        cnom[i] = new int[n];
 
    int**cdenom = new int*[m];
    for (int i = 0; i<m; i++)
        cdenom[i] = new int[n];
 
    int *Bnom = new int[m];
    int *Bdenom = new int[m];
 
    int *X = new int[n];
    int *X2 = new int[n];
 
    input(nom, denom, m, n, Bnom, Bdenom, cnom, cdenom);
    Gauss(nom, denom, cnom, cdenom, Bnom, Bdenom, X, X2, n);
}
__________________
Помощь в написании контрольных, курсовых и дипломных работ, диссертаций здесь
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
14.04.2016, 09:53
Ответы с готовыми решениями:

Решение СЛАУ методом Гаусса
У меня программа для решения слау методом Гаусса с выбором главного элемента по столбцам. Что-то...

Решение СЛАУ методом Гаусса
В институте задали написать программу решения СЛАУ методом Гаусса, но я абсолютный чайник (знаю,...

Решение СЛАУ методом Гаусса
Всем привет))) Друзья, сколько я уже лазию по инету и по этому форуму, я никак не могу найти...

Решение СЛАУ методом Гаусса
помогите разобраться!!ВЫдает 85 ошибок!!!! #include &lt;stdio.h&gt; #include &lt;conio.h&gt; #include...

0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
14.04.2016, 09:53
Помогаю со студенческими работами здесь

Реализовать решение СЛАУ методом Гаусса
Не могу разобраться в чём проблема. Пытаюсь реализовать метод гаусса, вроде бы всё правильно, но...

Решение СЛАУ методом Гаусса, индексы
Собственно программа решает слау методом гаусса. Как найти индексы главного элемента на второй...

Решение СЛАУ методом Гаусса с выбором ведущего элемента
1) Написать вычислительную программу на языке программирования С++ для решения СЛАУ методом Гаусса...

Найти решение СЛАУ методом Гаусса (без перестановки строк)
Итак, требуется найти решение СЛАУ методом Гаусса (без перестановки строк). Проще говоря, есть...

Почему программа прерывается с ошибкой? Реализовывал решение СЛАУ методом Гаусса
#include &lt;iostream&gt; #include &lt;cmath&gt; using namespace std; int main() { int n=0; int...

Численное решение СЛАУ методом Гаусса с организацией хранения матрицы в виде одномерного массива
Помогите кому не сложно.


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

Или воспользуйтесь поиском по форуму:
1
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2022, CyberForum.ru