Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
fitochay
14 / 14 / 1
Регистрация: 14.03.2015
Сообщений: 113
1

Метод Гаусса. Выполняется не верно

19.09.2015, 21:16. Просмотров 141. Ответов 0
Метки нет (Все метки)

Здравствуйте. Захотелось написать решение слау методом Гаусса. Можете подсказать что я сделал не так?

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
#include <iostream>
 
using namespace std;
 
const int N=3;
 
class Matrix
{
private:
    float **m;
public:
    void CreateMatrix();
    void SetMatrix();
    void PrintMatrix();
    void RMatrix();
};
 
void Matrix :: CreateMatrix()
{
    m = new float * [N];
    for (int i(0); i<N; i++)
        m[i] = new float [N+1];
}
 
void Matrix :: SetMatrix()
{
    for (int i(0); i<N; i++)
        for (int j(0); j<N+1; j++)
        {
            cout << "a[" << i << "][" << j << "]: ";
            cin >> m[i][j];
        }
}
 
void Matrix :: PrintMatrix()
{
    for (int i(0); i<N; i++)
        for (int j(0); j<N+1; j++)
            cout << endl << "a[" << i << "][" << j << "]: " << m[i][j];
}
 
void Matrix :: RMatrix()
{
    if (m[0][0] != 1)
    {
        float k = m[0][0];
        for (int j(0); j< N+1; j++)
            m[0][j] /= k;
    }
    else
    {
        for (int col(0); col<N-1; col++)
        {
            for (int row(col+1); row < N; row++)
            {
                float k = m[row][col]/m[col][col];
                float l;
                for (int j(col);j< N+1; j++)
                {
                    m[row][j] /= k;
                    m[row][j] -= m[col][j];
                }
                if ((m[row][row] != 1)&&(m[col][col] != 0))
                {
                    l = m[row][row];
                    for (int j(row); j<N+1; j++)
                        m[row][j] /= l;
                }
            }
        }
    }
    for (int col(N-1); col > 0; col --)
    {
        for (int row(col -1); row >= 0; row --)
        {
            float k = m[row][col];
            for (int j(col); j < N+1; j++)
            {
                m[row][j] -= m[col][j]*k;
            }
        }
    }
}
 
int main()
{
    Matrix a;
    a.CreateMatrix();
    a.SetMatrix();
    a.RMatrix();
    a.PrintMatrix();
    return 0;
}
Ведет себя вполне адекватно с матрицами где нет отрицательных чисел, но в остальных случаях все печально.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
19.09.2015, 21:16
Ответы с готовыми решениями:

СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
Помогите ребят. Не могу построить алгоритмы для этих методов Язык C++

метод гаусса
найти метод гаусса n+1 значение определителя в n+1 точке, а затем решить...

метод Гаусса
дана система линейных уравнений, в графическом режиме надо написать программу....

Метод Гаусса
Добрый вечер, можете написать программу для решения СЛАУ методом гаусса на Си,...

Метод Гаусса
Здравствуйте, у меня такая проблема. Мне нужно написать программу, которая...

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
19.09.2015, 21:16

Метод Гаусса в С++
Подскажите please, как запустить эту программу, через C++??? Решение системы...

Метод Гаусса
Всем привет. Ребят, что скажите по поводу такой реализации метода Гаусса: ...

Метод Гаусса
Как решить квадратную матрицу методом Гаусса, я ее привел к треугольному виду,...


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

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

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru