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

Метод простой итерации - C++

Восстановить пароль Регистрация
 
 
Рейтинг: Рейтинг темы: голосов - 17, средняя оценка - 4.82
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
25.02.2014, 18:15     Метод простой итерации #1
Всем привет. Кто знаком с методом простой итерации???
Вот есть код, он вроде работает, но в итоге ничего не дает, в коде описаны комментарии к работе, подскажите что не так??
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
#include "stdafx.h"
#include "iostream"
#include "vector"
#include "cmath"
#include "locale.h"
 
using namespace std;
 
int main()
{
    // Считываем размер вводимой матрицы
    setlocale(LC_ALL, "rus");
    int size;
    cout << "Введите размер вводимой матрицы: ";
    cin >> size;
    cout << endl;
    
    // Будем хранить матрицу в векторе, состоящем из 
    // векторов вещественных чисел
    vector <vector <long double> > matrix;
    
    // Матрица будет иметь размер (size) x (size + 1),
    // c учетом столбца свободных членов    
    matrix.resize (size);
    cout << "Введите элементы матрицы " << size << "x" << size+1 <<":" <<endl;
    for (int i = 0; i < size; i++)
    {
        matrix[i].resize (size + 1);
 
        for (int j = 0; j < size + 1; j++)
        {
            cin >> matrix[i][j];
        }
    }
 
    // Считываем необходимую точность решения
    long double eps;
    cout << endl;
    cout << "Введите точность решения: ";
    cin >> eps;
 
    // Введем вектор значений неизвестных на предыдущей итерации,
    // размер которого равен числу строк в матрице, т.е. size,
    // причем согласно методу изначально заполняем его нулями
    vector <long double> previousVariableValues (size, 0.0);
 
    // Будем выполнять итерационный процесс до тех пор, 
    // пока не будет достигнута необходимая точность    
    while (true)
    {
        // Введем вектор значений неизвестных на текущем шаге       
        vector <long double> currentVariableValues (size);
 
        // Посчитаем значения неизвестных на текущей итерации
        // в соответствии с теоретическими формулами
        for (int i = 0; i < size; i++)
        {
            // Инициализируем i-ую неизвестную значением 
            // свободного члена i-ой строки матрицы
            currentVariableValues[i] = matrix[i][size];
 
            // Вычитаем сумму по всем отличным от i-ой неизвестным
            for (int j = 0; j < size; j++)
            {
                if (i != j)
                {
                    currentVariableValues[i] -= matrix[i][j] * previousVariableValues[j];
                }
            }
 
            // Делим на коэффициент при i-ой неизвестной
            currentVariableValues[i] /= matrix[i][i];
        }
 
        // Посчитаем текущую погрешность относительно предыдущей итерации
        long double error = 0.0;
        
        for (int i = 0; i < size; i++)
        {
            error += abs (currentVariableValues[i] - previousVariableValues[i]);
        }
 
        // Если необходимая точность достигнута, то завершаем процесс
        if (error < eps)
        {
            break;
        }
 
        // Переходим к следующей итерации, так 
        // что текущие значения неизвестных 
        // становятся значениями на предыдущей итерации
        previousVariableValues = currentVariableValues;
    }
 
    // Выводим найденные значения неизвестных с 8 знаками точности
    for (int i = 0; i < size; i++)
    {
        printf ("%.8llf ", previousVariableValues[i]);
    }
    cout << endl;
    system("pause");
    return 0;   
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
25.02.2014, 18:15     Метод простой итерации
Посмотрите здесь:

Метод простой итерации C++
Метод простой итерации C++
метод простой итерации C++
[C++] Метод простой итерации C++
C++ Метод простой итерации
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
25.02.2014, 18:19     Метод простой итерации #2
Alexandr1966, тестовый пример выкладывайте
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
25.02.2014, 18:24  [ТС]     Метод простой итерации #3
Вот тестовый пример
Миниатюры
Метод простой итерации  
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
25.02.2014, 18:27     Метод простой итерации #4
Alexandr1966, мне нужен пример, который уже решен методом простой итерации. То есть условие и решение.

Не по теме:

А как вы программируете мат. методы без тестовых примеров?

Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
26.02.2014, 14:13  [ТС]     Метод простой итерации #5
Ev_Hyper, вот смотрите:

Матрица А:
-0,1 1 0,1 0,4
0,2 0,4 0,3 1
1 -0,1 0,1 0,2
0,1 0,1 1 0,3

Матрица b:
2
1
1
2

Ответы матрица Х:
1,1680
2,1995
1,8651
-0,6729

Добавлено через 2 минуты
Я не знаю, для чего в коде погрешность и размер матрицы (size * size+1), в общем что то здесь не так

Добавлено через 6 часов 43 минуты
Ребята помогите кто знает
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
26.02.2014, 14:43     Метод простой итерации #6
Alexandr1966, откуда вы этот код брали?
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
26.02.2014, 16:28  [ТС]     Метод простой итерации #7
Ev_Hyper, нашел в инете где-то или на компах в универе

Добавлено через 16 минут
Ev_Hyper, помоги сделать пожалуйста правильно это метод
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
26.02.2014, 19:09     Метод простой итерации #8
Alexandr1966, какие сроки?

Добавлено через 14 минут
Alexandr1966, метод работает правильно.
Код
Введите размер вводимой матрицы: 4

Введите элементы матрицы 4x5:
3
0
0
-3
21
0
11
1
2
21
1
0
5
0
48
-1
2
7
11
61
3 0 0 -3 21
0 11 1 2 21
1 0 5 0 48
-1 2 7 11 61

Введите точность решения: 0.001
8,00070948 0,99990793 7,99939336 1,00062094
Для продолжения нажмите любую клавишу . . .
Почитайте:
http://dxdy.ru/topic46434.html

То, есть тут просто неправильный выход из итераций. И нет проверки на корректность входящей системы.
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
27.02.2014, 03:45  [ТС]     Метод простой итерации #9
Ev_Hyper, а почему матрица 4 на 5 вводиться, мне надо чтоб просто 4 на 4 вводилась и матрица В
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
27.02.2014, 09:23     Метод простой итерации #10
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
Ev_Hyper, а почему матрица 4 на 5 вводиться, мне надо чтоб просто 4 на 4 вводилась и матрица В
Alexandr1966, так захотел автор кода. Переделывайте так, как вам удобно. Тут же сложность только в правильном программировании алгоритма.
S_el
1907 / 1502 / 296
Регистрация: 15.12.2013
Сообщений: 5,918
27.02.2014, 09:38     Метод простой итерации #11
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
а почему матрица 4 на 5 вводиться, мне надо чтоб просто 4 на 4 вводилась и матрица В
Сложнее работать с 2-мя массивами,чем с одним.Оставляйте как есть.
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
27.02.2014, 14:39  [ТС]     Метод простой итерации #12
Ev_Hyper, вот смотрите:
как я понял, почему размер матрицы size * size+1 - потому что программа предлагает вести сразу расширенную матрицу. Так вот, в том примере, который я вам скидывал, я так и ввел, потом ввел точность решения, нажал Enter, а в ответ ничего. Вопрос, в чём дело???

Матрица А:
-0,1 1 0,1 0,4
0,2 0,4 0,3 1
1 -0,1 0,1 0,2
0,1 0,1 1 0,3

Матрица b:
2
1
1
2

Вот как программа работает, ничего не выдает в итоге, а у вас почему то выдал ответ.
Миниатюры
Метод простой итерации  
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
27.02.2014, 14:49  [ТС]     Метод простой итерации #13
S_el, взгляните пожалуйста, что не так тут
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
27.02.2014, 16:59     Метод простой итерации #14
Alexandr1966, я вам еще раз говорю - метод специфический. Вы не все системы сможите им решить.
Добавьте выход из цикла после 200 итерации, или через погрешность. А еще лучше попробуйте сделать проверку на корректность входных данных. А может даже и правильное преобразование.
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
27.02.2014, 18:43  [ТС]     Метод простой итерации #15
Ev_Hyper, дело в том что я вообще понимаю этого кода, и мне просто нужно сделать чтоб работало, простенько как нибудь, мне лишь бы сдать, помогите прошу, вообще нервоз
S_el
1907 / 1502 / 296
Регистрация: 15.12.2013
Сообщений: 5,918
27.02.2014, 19:10     Метод простой итерации #16
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
Ev_Hyper, дело в том что я вообще понимаю этого кода, и мне просто нужно сделать чтоб работало, простенько как нибудь, мне лишь бы сдать, помогите прошу, вообще нервоз
Почему вы решили,что этот метод работает неправильно?
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
27.02.2014, 19:11     Метод простой итерации #17
Alexandr1966, препод знает, что метод работает не всегда. Он скорее спросит, почему такое ограничение.

У вас те значения по заданию такие?
Alexandr1966
25 / 25 / 1
Регистрация: 17.12.2012
Сообщений: 425
28.02.2014, 08:07  [ТС]     Метод простой итерации #18
Ev_Hyper, преподу плевать если честно)
главное чтоб вот мой конкретный пример работал, который я выше указал)))
S_el
1907 / 1502 / 296
Регистрация: 15.12.2013
Сообщений: 5,918
28.02.2014, 09:19     Метод простой итерации #19
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
главное чтоб вот мой конкретный пример работал, который я выше указал)))
А если он специально подобран,что его не решить методом простой итерации?
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
28.02.2014, 09:31     Метод простой итерации
Еще ссылки по теме:

Метод простой итерации C++
Метод простой итерации C++
Метод простой итерации C++

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

Или воспользуйтесь поиском по форуму:
Ev_Hyper
 Аватар для Ev_Hyper
1806 / 1627 / 435
Регистрация: 15.12.2013
Сообщений: 5,790
28.02.2014, 09:31     Метод простой итерации #20
Alexandr1966, вы хотите, чтобы я за вас "исправил" программу, так, чтобы она свела конкретно ваш пример к корректному виду?
Yandex
Объявления
28.02.2014, 09:31     Метод простой итерации
Ответ Создать тему
Опции темы

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