24 / 6 / 3
Регистрация: 06.12.2015
Сообщений: 308
1

Фрагмент из метода Зейделя работает почему-то неправильно

09.12.2016, 12:09. Показов 489. Ответов 1
Метки нет (Все метки)

Значения на первой итерации высчитывает верно x[1] = -0.8, y[1] = -2.05, а последующие итерации выводит с ошибкой.
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
    
double A [2][3] = {{-5,-1,-4},{-4,4,-5}};//матрица с коэффициентами при двух неизвестных
double epsilon = 0.001;
double x[3], y[3]; // массивы для значений x и y 
    x[0] = 0, y[0] = 0;
do
    {
        int i = 0;
        i++;
        x[i + 1] = (A[i][i + 2] - A[i][i + 1] * y[i]) / A[i][i];
        y[i + 1] = (A[i + 1][i + 2] - A[i + 1][i] * x[i + 1]) / A[i + 1][i + 1];
        y[i + 1] = y[i];
        cout << "Решение на " << i << "-ой " << "итерации - "
            << "x = " << round(x[i + 1] * 10000) / 10000 << "y =" << round(y[i + 1] * 10000) / 10000 << "\n";
    } while (x[i + 1]<epsilon);
__________________
Помощь в написании контрольных, курсовых и дипломных работ, диссертаций здесь
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
09.12.2016, 12:09
Ответы с готовыми решениями:

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

почему не работает фрагмент?
Программа считывает данные из файла, заполняет ими массив bf static char bf; ReadFile(fl,...

Почему смена текста работает только один раз? Что неправильно в коде и почему?
Нужно чтобы при щелчке мыши на поверхности кнопки текст метки Label изменялся на текст окна...

Нужно перевести программу из паскаль в делфи . Сама программа этот реализация метода гауса и метода зейделя
Это все в паскале сделано. В программе идет в начале описание меня а аотом сама подпрограмма...

1
1481 / 1198 / 819
Регистрация: 29.02.2016
Сообщений: 3,579
09.12.2016, 20:02 2
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
#include <iostream>
#include <time.h>
 
using namespace std;
 
double GetNorm(int ndim, double *vector)
{
    double result = 0.0;
    for (int i = 0; i < ndim; i++)
        result += vector[i] * vector[i];
    return sqrt(result);
}
void  Seidel(int size, double **A, double *b, int MaxIterations, double tolerance , double *x)
{
    const double epsilon = 1.0e-500;
    double  *dx, *xOld;
 
    dx = new double[size];
    xOld = new double[size];
    for(int i = 0; i < size; ++i)
        x[i] = 0.;
    double entry, diagonal;
    do
        {
            for(int i = 0; i < size; ++i)
                xOld[i]=x[i];
            for (int i = 0; i < size; i++)
            {
                entry = b[i];
                diagonal = A[i][i];
                if (abs(diagonal) < epsilon){
                    cout << "Diagonal element is too small!" << endl;
                    system("pause");
                    exit(0);
                }
                for (int j = 0; j < i; j++)
                    entry -= A[i][j] * x[j];
                for (int j = i + 1; j < size; j++)
                    entry -= A[i][j] * xOld[j];
                x[i] = entry / diagonal;
            }
            for(int i = 0; i < size; ++i)
                dx[i] = x[i] - xOld[i];
        }
        while (GetNorm(size,dx) > tolerance);
    }
 
int main()
{
    double A[3][3] = {{4, 0, 1},{0, 3, 2},{1, 2, 4}};
    double b[3] = { 2, 1, 3 };
    int n=3;
    double **A1;
    A1 = new double *[n];
    double *b1;
    b1 = new double[n];
    double *x1;
    x1 = new double[n];
    for(int i = 0; i < n; ++i)
    {
        A1[i] = new double [n];
        b1[i]= b[i];
    }
    for(int i = 0; i < n; ++i)
        for(int j = 0; j < n; ++j)
            A1[i][j] = A[i][j];
 
    Seidel(n, A1, b1, 10, 1.0e-4, x1);
    cout << "Solution:"<<endl;
    cout << " x1[0] = " << x1[0] << endl;
    cout << " x1[1] = " << x1[1] << endl;
    cout << " x1[2] = " << x1[2] << endl;
    system("pause");
    return 0;
}
Добавлено через 17 минут
кое что исправил
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
#include <iostream>
#include <time.h>
 
using namespace std;
 
double GetNorm(int ndim, double *vector)
{
    double result = 0.0;
    for (int i = 0; i < ndim; i++)
        result += vector[i] * vector[i];
    return sqrt(result);
}
void  Seidel(int size, double **A, double *b, int &MaxIterations, double tolerance , double *x)
{
    const double epsilon = 1.0e-500;
    double  *dx, *xOld;
 
    dx = new double[size];
    xOld = new double[size];
    for(int i = 0; i < size; ++i)
        x[i] = 0.;
    double entry, diagonal;
    do
        {
            for(int i = 0; i < size; ++i)
                xOld[i]=x[i];
            for (int i = 0; i < size; i++)
            {
                entry = b[i];
                diagonal = A[i][i];
                if (abs(diagonal) < epsilon){
                    cout << "Diagonal element is too small!" << endl;
                    system("pause");
                    exit(0);
                }
                for (int j = 0; j < i; j++)
                    entry -= A[i][j] * x[j];
                for (int j = i + 1; j < size; j++)
                    entry -= A[i][j] * xOld[j];
                x[i] = entry / diagonal;
            }
            for(int i = 0; i < size; ++i)
                dx[i] = x[i] - xOld[i];
            MaxIterations++;
        }
        while (GetNorm(size,dx) > tolerance);
    }
 
int main()
{
    double A[3][3] = {{4, 0, 1},{0, 3, 2},{1, 2, 4}};
    double b[3] = { 2, 1, 3 };
    int n=3;
    int MaxIterations=0;
    double **A1;
    A1 = new double *[n];
    double *b1;
    b1 = new double[n];
    double *x1;
    x1 = new double[n];
    for(int i = 0; i < n; ++i)
    {
        A1[i] = new double [n];
        b1[i]= b[i];
    }
    for(int i = 0; i < n; ++i)
        for(int j = 0; j < n; ++j)
            A1[i][j] = A[i][j];
 
    Seidel(n, A1, b1, MaxIterations, 1.0e-4, x1);
    cout << "MaxIterations = " << MaxIterations << endl;
    cout << "Solution:"<<endl;
    cout << " x1[0] = " << x1[0] << endl;
    cout << " x1[1] = " << x1[1] << endl;
    cout << " x1[2] = " << x1[2] << endl;
    system("pause");
    return 0;
}
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
09.12.2016, 20:02
Помогаю со студенческими работами здесь

Объясните, пожалуйста, как работает фрагмент кода. Почему при i = 0 вводится лишь 2 строки?
Здравствуйте! Объясните, пожалуйста, как происходит ввод строк в этом фрагменте кода. string...

Почему неправильно работает сортировка?
сортировка: for i := 1 to n-1 do begin if gruppa1.grupa&gt;gruppa1.grupa then begin ...

Почему программа работает неправильно
не могу понять почему программа работает не правильно # include &lt;stdio.h&gt; # include &lt;conio.h&gt; #...

Почему программа неправильно работает
string IntToStr(int a) { if (a == 0) return &quot;0&quot;; string str = &quot;&quot;; int b = a; if (a &lt; 0)...


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

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

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