14 / 14 / 1
Регистрация: 14.03.2015
Сообщений: 113
1

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

19.09.2015, 21:16. Показов 432. Ответов 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
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
19.09.2015, 21:16
Ответы с готовыми решениями:

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

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

метод гаусса
Помогите,пожалуйста Преподаватель попросил расписать программу по блокам #include &quot;stdio.h&quot;...

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

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

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

Метод Гаусса
Всем привет. Ребят, что скажите по поводу такой реализации метода Гаусса: #include &lt;stdio.h&gt;...

Метод Гаусса!!!
Всем привет! )) У меня вот такая проблемка! Нужно решить СЛАУ методом Гаусса с выбором главного...

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


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

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

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