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

Метод Гаусса-Зейделя - C++

Восстановить пароль Регистрация
 
Yuregg
Сообщений: n/a
06.01.2014, 17:07     Метод Гаусса-Зейделя #1
Ребят, нужна помощь. Программа вроде бы написана правильно, корни выдает подходящие. По идее, должно выполняться около 10-20 итераций, для нахождения корней. Но, вместо этого, выполняется от 40 до 600 итераций. Если выполняется более 100 итераций, корни выдает в виде -1.IND. Что в программе написано не так?
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
#include "stdafx.h"
#include "math.h"
#include <iostream>
#include <iomanip>
#include <time.h>
using namespace std;
 
bool converge(double *xk, double *xkp, double eps, int n)
{
    for (int i = 0; i < n; i++) 
    {
        if (fabs(xk[i] - xkp[i]) >= eps) 
            return false;
    }
    return true;
}
 
int _tmain(int argc, _TCHAR* argv[])
{
    srand((unsigned)time(NULL));
    double *p, *x, **a, *b, eps;
    int n,i,j,k=0;
    cin >> n;
    cin >> eps,
    p = new double[n];
    x = new double[n];
    a = new double*[n];
    b = new double[n];
    for (i = 0; i < n; i++)
    {
        a[i] = new double[n];
        b[i] = ((double)rand()/RAND_MAX)*100;
        x[i] = 0.0;
        for (j = 0; j < n; j++)
        a[i][j] = ((double)rand()/RAND_MAX)*100;
    
    }
    cout << setw (13) << "A" << setw (16) << "B";
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout<< endl;
        for (j = 0; j < n; j++)
            cout << a[i][j] << " ";
        cout << " " << b[i];
    }
    do
    {
        k++;
        for (int i = 0; i < n; i++)
         p[i] = x[i];
 
        for (int i = 0; i < n; i++)
        {
            double var = 0;
            for (int j = 0; j < i; j++)
                var += (a[i][j] * x[j]);
            for (int j = i + 1; j < n; j++)
                var += (a[i][j] * p[j]);
            x[i] = (b[i] - var) / a[i][i];
        }
    }
    while (!converge(x, p, eps, n));
    cout<< endl;
    for (i = 0; i < n; i++)
    {
        cout << x[i] << endl;
    }
    cout << k << endl;
    system ("pause");
 
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.01.2014, 17:07     Метод Гаусса-Зейделя
Посмотрите здесь:

параллельный алгоритм решенный методом Гаусса-Зейделя C++
Метод Зейделя C++
C++ Метод итераций и метод Зейделя
метод зейделя C++
Метод Зейделя C++
C++ СЛАУ. Метод обратной матрицы, метод Гаусса, метод Крамера, метод Зейделя
C++ Метод Зейделя-Гаусса С++
Метод Гаусса-Зейделя C++

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

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

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