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

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 16, средняя оценка - 4.75
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
#1

Нужен рабочий пример решения СЛАУ методом Гаусса - C++

21.05.2014, 05:25. Просмотров 2193. Ответов 18
Метки нет (Все метки)

Нужен рабочий (ну или более менее рабочий) код для решения СЛАУ методом Гауса на С++.
Уже замучился искать
0
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
21.05.2014, 05:25
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Нужен рабочий пример решения СЛАУ методом Гаусса (C++):

Решения СЛАУ методом Гаусса по шагу - C++
Добрый день! Очередной раз прошу советов! Имеется СЛАУ, которое требуется решить методом Гаусса на с++. Соответственно, на...

Разработать программу для решения СЛАУ методом Гаусса. - C++
Разработать программу для решения СЛАУ методом Гаусса. вот заданный массив 7x1-x2+2x3+x4+x5=7 2x1+5x2-x3-x5=3 ...

Помогите написать программу для решения слау методом Гаусса - C++
Осталось как обычно пару дней до сдачи курсовой, а она еще не готова! Помогите пожалуйста. Вот код, который чуть-чуть работает Главная...

Кто сможет? Написать функцию для решения СЛАУ методом Гаусса - C++
Кто может написать функцию для решения СЛАУ методом Гаусса используя базовые знания.

Метод Гаусса для решения СЛАУ - C++
Здравствуйте! У меня есть алгоритм решения СЛАУ методом Гаусса, но почему-то он не работает((( Я проверил промежуточные матрицы, все ОК! И...

Метод Гаусса (подстановки) решения СЛАУ - C++
Всем доброго времени суток, написал я одну недопрограмму, потому как сделал её наполовину. Не хватает мне опыта и знания математики и опыта...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 797
Записей в блоге: 10
21.05.2014, 05:40 #2
Зайди ко мне в блог ;
Там есть метод Йордана-Гаусса . Это почти тот же метод
только мы зануляем заодно и верхний треугольник
0
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
21.05.2014, 07:15  [ТС] #3
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
124
125
126
127
128
129
130
131
#include <stdlib.h>//malloc, system("pause")
#include <stdio.h>  //i/o
#include <conio.h>  //getch
#include <math.h>
 
void ShowVector(int n, double * vec)
{
        for(int i = 0; i < n; i++)
                printf("%.3f ",vec[i]);
        printf("\r\n");
}
 
void PryamoiHod(int n, double **a, double *b)
{
        double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                        v     = b[im];
                        b[im] = b[k];
                        b[k]  = v;
                }
                for(i = k + 1; i < n; i++)
                {
                        v               = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        b[i]    = b[i] - v*b[k];
                                                if(v != 0)
                        for(j = k + 1; j < n; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}
 
void ObratniHod(int n, double **a, double *b, double *x)
{
        double s = 0;
        x[n - 1] = 1.0*b[n - 1]/a[n - 1][n - 1];
        for(int i = n - 2, j; 0 <= i; i--)
        {
                s = 0;
                for(j = i + 1; j < n; j++)
                {
                        s = s+a[i][j]*x[j];
                }
                x[i] = 1.0*(b[i] - s)/a[i][i];
        }
}
 
int main()
{
        int i,j,n;
        double **a, *b, *x;
        do
        {
                printf("Enter NUM of equations: ");
                scanf("%d",&n);
                //Выделяем память под матрицу А и векторы В и Х
                a = (double **)malloc(n*sizeof(double));
                b = (double *)malloc(n*sizeof(double));
                x = (double *)malloc(n*sizeof(double));
                for(i = 0; i < n; i++)
                {
                        a[i] = (double *)malloc(n*sizeof(double));
                        //Ввод a
                        for(j = 0; j < n; j++)
                        {
                                printf("a[%d][%d] = ",i + 1,j + 1);
                                scanf("%lf",&a[i][j]);
                        }
                }
                //Ввод b
                for(i = 0; i < n; i++)
                {
                        printf("b[%d] = ",i + 1);
                        scanf("%lf",&b[i]);
                }
                
                printf("\tSee input\r\n");
                printf("Matrix A:\r\n");
                for(i = 0; i < n; i++)
                        ShowVector(n, a[i]);
                printf("Vector B:\r\n");
                ShowVector(n, b);
                
                printf("\tSolving on Gauss method\r\n");
                PryamoiHod(n, a, b);
                printf("Forvard Gauss course\r\n");//Прямой ход
                printf("Matrix A:\r\n");
                for(i = 0; i < n; i++)
                        ShowVector(n, a[i]);
                printf("Vector B:\r\n");
                ShowVector(n, b);
 
                ObratniHod(n, a, b, x);
                printf("Back Gauss course\r\n");//Обратный ход
                printf("Matrix A:\r\n");
                for(i = 0; i < n; i++)
                        ShowVector(n, a[i]);
                printf("Vector B:\r\n");
                ShowVector(n, b);
 
                printf("Results :\r\n");
                ShowVector(n, x);
 
                printf("Press Y for new input\r\n");
                //Чистим память
                free((void *)a);
                free((void *)b);
                free((void *)x);
        }
        while(toupper(getch()) == 'Y'); // ????????????????
        return 0;
}
Вот нашел, вроде норм работает.
Может кто то закомментировать?

Добавлено через 45 минут
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
void PryamoiHod(int n, double **a, double *b)
{
        double v;
        for(int k = 0,i,j,im; k < n - 1; k++)
        {
                im = k;
                for(i = k + 1; i < n; i++)
                {
                        if(fabs(a[im][k]) < fabs(a[i][k]))
                        {
                                im = i;
                        }
                }
                if(im != k)
                {
                        for(j = 0; j < n; j++)
                        {
                                v                = a[im][j];
                                a[im][j] = a[k][j];
                                a[k][j]  = v;
                        }
                        v     = b[im];
                        b[im] = b[k];
                        b[k]  = v;
                }
                for(i = k + 1; i < n; i++)
                {
                        v               = 1.0*a[i][k]/a[k][k];
                        a[i][k] = 0;
                        b[i]    = b[i] - v*b[k];
                                                if(v != 0)
                        for(j = k + 1; j < n; j++)
                        {
                                a[i][j] = a[i][j] - v*a[k][j];
                        }
                }
        }
}
Про комментируйте что тут происходит (с математической точки зрения), плс

Добавлено через 32 секунды
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
void ObratniHod(int n, double **a, double *b, double *x)
{
        double s = 0;
        x[n - 1] = 1.0*b[n - 1]/a[n - 1][n - 1];
        for(int i = n - 2, j; 0 <= i; i--)
        {
                s = 0;
                for(j = i + 1; j < n; j++)
                {
                        s = s+a[i][j]*x[j];
                }
                x[i] = 1.0*(b[i] - s)/a[i][i];
        }
}
и тут...

Добавлено через 2 минуты
Как прописать эти же процедуры с помощью new и delete
C++
1
2
3
                a = (double **)malloc(n*sizeof(double));
                b = (double *)malloc(n*sizeof(double));
                x = (double *)malloc(n*sizeof(double));
C++
1
2
3
4
                //Чистим память
                free((void *)a);
                free((void *)b);
                free((void *)x);
Добавлено через 21 секунду
help me, please =)
0
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 797
Записей в блоге: 10
21.05.2014, 15:17 #4
C++
1
2
3
a = new double[n];
b = new double[n];
x = new double[n];
Добавлено через 34 секунды
C++
1
2
3
delete [] a;
delete [] b;
delete [] x;
Добавлено через 2 минуты
Цитата Сообщение от andreyananas Посмотреть сообщение
с математической точки зрения
Алгоритм Гаусса сам знаете ?
Получив треугольную матрицу
мы выводим элементы уже (Сами значения)
снизу вверх сперва http://www.cyberforum.ru/cgi-bin/latex.cgi?x_{n} , потом http://www.cyberforum.ru/cgi-bin/latex.cgi?x_{n - 1} , и так далее ..
0
-=ЮрА=-
Заблокирован
Автор FAQ
21.05.2014, 15:42 #5
Цитата Сообщение от andreyananas Посмотреть сообщение
Вот нашел, вроде норм работает.
Может кто то закомментировать?
- когда ты брал код из этой ветки Методом Гаусса решить систему n линейных алгебраических уравнений, удосужился бы посмотреть пост 4 той же темы. Совет : почитай FAQ по математике прежде чем задавать вопросы http://www.cyberforum.ru/faq/thread436065.html

Добавлено через 5 минут
Цитата Сообщение от Qazan Посмотреть сообщение
a = new double[n];
- неверно т.к а предствляет собой матрицу а не вектор
Цитата Сообщение от andreyananas Посмотреть сообщение
a = (double **)malloc(n*sizeof(double));
int ** a = new double * [n];
2
Qazan
211 / 59 / 9
Регистрация: 30.04.2013
Сообщений: 797
Записей в блоге: 10
21.05.2014, 15:45 #6
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
Сообщение от Qazan
a = new double[n];
да , точно
0
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
22.05.2014, 11:47  [ТС] #7
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
- когда ты брал код из этой ветки Методом Гаусса решить систему n линейных алгебраических уравнений, удосужился бы посмотреть пост 4 той же темы.
Я код не там взял, а с твоего поста Метод Гаусса с визуализацией

Добавлено через 1 час 12 минут
Цитата Сообщение от -=ЮрА=- Посмотреть сообщение
int ** a = new double * [n];
Почему int???? а не double **a = new double *[n]??
0
-=ЮрА=-
Заблокирован
Автор FAQ
22.05.2014, 21:00 #8
Цитата Сообщение от andreyananas Посмотреть сообщение
Почему int???? а не double **a = new double *[n]??
- да да double ** ты всё верно заметил.
0
aiwprton805
70 / 70 / 26
Регистрация: 30.03.2013
Сообщений: 189
Завершенные тесты: 1
22.05.2014, 21:29 #9
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
double *gauss::solution()
{
    for(int i=0; i<n; i++){
        double tmp = matrix[i][i];
        for(int j=n; j>=i; j--){
            matrix[i][j] /= tmp;}
        for(int j=i+1; j<n; j++){
            tmp = matrix[j][i];
            for(int k=n; k>=i; k--){
                matrix[j][k] -= tmp * matrix[i][k];}
        }
    }
 
    double *vector = new double[n];
    vector[n-1] = matrix[n-1][n];
    for(int i=n-2; i>=0; i--){
        vector[i] = matrix[i][n];
        for(int j=i+1; j<n; j++){
            vector[i] -= matrix[i][j] * vector[j];}
    }
    return vector;
}
0
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
24.05.2014, 01:12  [ТС] #10
Можете помочь с модернизации кода.
Программа должна записывать в файл СЛАУ и считывать с него.

Добавлено через 11 часов 26 минут
ап хелп плс)
0
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
24.05.2014, 13:01  [ТС] #11
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
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
#include <iostream>
#include <stdio.h>
#include <windows.h>
#include <conio.h>
#include <math.h>
 
using namespace std;
 
void Show(int, double*);
void PryamojHod(int, double**, double*);
void ObratnyHod(int, double**, double*, double*);
 
int main()
{
    SetConsoleOutputCP(1251);
    int n;
    cout << "\n Введите количество переменных: ";
    while(!(cin >> n) || (cin.peek() != '\n') || (n <=0))
    {
        cin.clear();
        while(cin.get() != '\n');
        cout << "\n\a Ошибка! \n для новой попытки, нажмите любую кнопку...";
        getch();
        system("cls");
        cout << "\n Введите количество переменных: ";
    }
 
    // выдиляем память под матрицу А и векторы В и Х
    double **a = new double *[n];
    for (int i = 0; i < n; i++)
        a[i] = new double [n];
    double *b = new double [n];
    double *x = new double [n];
 
    // ввод данных
    for(int i = 0; i < n; i++)
    {
        for(int j = 0; j < n; j++)
        {
            cout << " Введите " << (j+1) << "й коэффициент " << (i+1) << "го уравнения: ";
            while(!(cin >> a[i][j]) || (cin.peek() != '\n'))
            {
                cin.clear();
                while(cin.get() != '\n');
                cout << "\n Ошибка ввода!\n для новой попытки, нажмите любую кнопку...";
                getch();
                system("cls");
                cout << " Введите " << (j+1) << "й коэффициент " << (i+1) << "го уравнения: ";
            }
        }
        cout << " Введите " << (i+1) << "й свободный член " << (i+1) << "го уравнения: ";
        while(!(cin >> b[i]) || (cin.peek() != '\n'))
        {
            cin.clear();
            while(cin.get() != '\n');
            cout << "\n Ошибка ввода!\n для новой попытки, нажмите любую кнопку...";
            getch();
            system("cls");
            cout << " Введите " << (i+1) << "й свободный член " << (i+1) << "го уравнения: ";
        }
    }
 
    // вывод СЛАУ
    cout << "\n Матрица А: " << endl;
    for(int i = 0; i < n; i++)
        Show(n, a[i]);
    cout << "\n Вектор B: " << endl;
    Show(n, b);
 
 
    PryamojHod(n, a, b);
    ObratnyHod(n, a, b, x);
 
    cout << "\n Результат: " << endl;
    Show(n, x);
 
    // освобождаем память
    for (int i = 0; i < n; i++)
        delete [] a[i];
    delete [] a;
    delete [] b;
    delete [] x;
 
    return 0;
}
 
void Show(int n, double *vec)
{
    for(int i = 0; i < n; i++)
    {
        cout << "\t" << vec[i];
    }
    cout << endl;
}
 
void PryamojHod(int n, double **a, double *b)
{
    double v = 0;
    for(int k(0), im; k < n; k++)
    {
        im = k;
        for(int i = k+1; i < n; i++)
        {
            if(fabs(a[im][k]) < fabs(a[i][k]))
                im = i;
        }
        if(im != k)
        {
            for(int j(0); j < n; j++)
            {
                v = a[im][j];
                a[im][j] = a[k][j];
                a[k][j] = v;
            }
            v = b[im];
            b[im] = b[k];
            b[k] = v;
        }
        for(int i = k + 1; i < n; i++)
        {
            v = 1.0*a[i][k]/a[k][k];
            a[i][k] = 0;
            b[i] = b[i]-v*b[k];
            if(v != 0)
            {
                for(int j = k + 1; j < n; j++)
                    a[i][j] = a[i][j]-v*a[k][j];
            }
        }
    }
}
 
void ObratnyHod(int n, double **a, double *b, double *x)
{
    double s = 0;
    x[n-1]=1.0*b[n-1]/a[n-1][n-1];
    for(int i = n - 2; i >= 0; i--)
    {
        s = 0;
        for(int j = i + 1; j < n; j++)
        {
            s=s+a[i][j]*x[j];
        }
        x[i]=1.0*(b[i]-s)/a[i][i];
    }
}
Вот программа, что то решает правильно, а что то нет(
(СЛАУ методом Гаусса)
Вот скрин, какой результат получается в некоторых случаях.
Почему так? Если такая СЛАУ не решается, тогда как прописать проверку: сможет программа решить СЛАУ или нет?
0
Миниатюры
Нужен рабочий пример решения СЛАУ методом Гаусса  
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
24.05.2014, 13:02  [ТС] #12
п.с. я ввел: 1 2 3 4 5 6 7 8 9 0 1 2 3 4 5 6 7 8 9 0
0
Vergil
42 / 36 / 14
Регистрация: 11.11.2009
Сообщений: 243
24.05.2014, 13:34 #13
У тебя две строчки одинаковых(не помню как называются), и 4 неизвестных. Должно быть 4 разные строки и 4 неизвестных.
В случае, если число уравнений в совместной системе получилось меньше числа неизвестных, то тогда ответ будет записываться в виде фундаментальной системы решений.
Фундаментальная система решений (ФСР) представляет собой набор линейно независимых решений однородной системы уравнений.
0
andreyananas
22 / 22 / 9
Регистрация: 15.10.2013
Сообщений: 862
Завершенные тесты: 2
24.05.2014, 14:15  [ТС] #14
Цитата Сообщение от Vergil Посмотреть сообщение
У тебя две строчки одинаковых(не помню как называются), и 4 неизвестных. Должно быть 4 разные строки и 4 неизвестных.
В случае, если число уравнений в совместной системе получилось меньше числа неизвестных, то тогда ответ будет записываться в виде фундаментальной системы решений.
Фундаментальная система решений (ФСР) представляет собой набор линейно независимых решений однородной системы уравнений.
Последний раз редактировалось Vergil; Сегодня в 12:38.
Можешь написать, как изобразить это в коде?
0
-=ЮрА=-
Заблокирован
Автор FAQ
24.05.2014, 17:18 #15
andreyananas, ты ввёл своим 1,2,34 систему у которой множество решений, я уже писал тебе прочти FAQ там оговрено что решение по Гауссу существует если определитель системы не равен нулю, а ты что ввёл?
Уже хотя бы систему составил бы!Хочешь результаты из файла и в файл почитай об freopen либо об перенаправлении вывода с помощью rdbuf тогда в коде надо будет всего лишь 2 строчки записать.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
24.05.2014, 17:18
Привет! Вот еще темы с ответами:

Метод Гаусса решения СЛАУ с полным выбором. C++ - C++
Написал код реализующий метод Гаусса с полным выбором ведущего элемента, но ответ выдает не правильный. Помогите найти ошибку #include...

СЛАУ методом Гаусса - C++
Всем доброго времени суток. По-сколько я новичек в программирование, то вот скачал код, для решения СЛАУ методом Гауса. Выдает...

Метод Гаусса для решения СЛАУ с использованием одномерного массива - C++
Доброго всем времени суток. У меня есть прога для двумерного массива #include &lt;stdio.h&gt; #include &lt;iostream&gt; #include &lt;conio.h&gt; ...

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


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

Или воспользуйтесь поиском по форуму:
Yandex
Объявления
24.05.2014, 17:18
Ответ Создать тему
Опции темы

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