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

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

Войти
Регистрация
Восстановить пароль
 
 
Рейтинг: Рейтинг темы: голосов - 17, средняя оценка - 4.82
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
#1

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

25.02.2014, 18:15. Просмотров 2538. Ответов 32
Метки нет (Все метки)

Всем привет. Кто знаком с методом простой итерации???
Вот есть код, он вроде работает, но в итоге ничего не дает, в коде описаны комментарии к работе, подскажите что не так??
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;   
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
25.02.2014, 18:15
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод простой итерации (C++):

Метод простой итерации - C++
Вот такая задача! Используя метод простой итерации, решить с точностью ε = 10 такие нелинейные системы уравнений. Помогите реализовать в...

Метод простой итерации - C++
Здраствуйте! Нужно реализовать в C++ метод простой итерации. #include &lt;iostream&gt; #include &lt;iomanip&gt; #include &lt;conio.h&gt; #include...

Метод простой итерации - C++
Найти приближенное решение уравнения x = f(x), где f(x) = α*exp(x), при α = 0, 2 методом простой итерации: xn+1 = f(xn). В качестве...

Метод простой итерации - C++
Помогите с рекурсией!!! Задача: Найти корень заданного уравнения методом простой итерации с заданной точностью. Напомним, что в этом...

Метод простой итерации - C++
Люди помогите пожалуйста очень нужно для курсачя, а в С++ не тум тум. Метод простой итерации найти все действительные корни...

Метод простой итерации - C++
решение системы линейных уравнений методом простой итерации

32
S_el
2113 / 1633 / 308
Регистрация: 15.12.2013
Сообщений: 6,577
27.02.2014, 19:10 #16
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
Ev_Hyper, дело в том что я вообще понимаю этого кода, и мне просто нужно сделать чтоб работало, простенько как нибудь, мне лишь бы сдать, помогите прошу, вообще нервоз
Почему вы решили,что этот метод работает неправильно?
0
Ev_Hyper
Заблокирован
27.02.2014, 19:11 #17
Alexandr1966, препод знает, что метод работает не всегда. Он скорее спросит, почему такое ограничение.

У вас те значения по заданию такие?
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
28.02.2014, 08:07  [ТС] #18
Ev_Hyper, преподу плевать если честно)
главное чтоб вот мой конкретный пример работал, который я выше указал)))
0
S_el
2113 / 1633 / 308
Регистрация: 15.12.2013
Сообщений: 6,577
28.02.2014, 09:19 #19
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
главное чтоб вот мой конкретный пример работал, который я выше указал)))
А если он специально подобран,что его не решить методом простой итерации?
0
Ev_Hyper
Заблокирован
28.02.2014, 09:31 #20
Alexandr1966, вы хотите, чтобы я за вас "исправил" программу, так, чтобы она свела конкретно ваш пример к корректному виду?
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
28.02.2014, 15:30  [ТС] #21
Ev_Hyper, S_el, парни, короче, я нашел в интернете код для метода простой итерации, судя по всему предыдущий код был код для метода Зейделя, вот для итерации, но он по прежнему не работает по примеру, который я привожу, может тут вы подскажите что не так, коды похожи, НО ЧЕМ-ТО ОТЛИЧАЮТСЯ.

Матрица А:
-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

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
#include <iostream>
#include <vector>
#include <cmath>
 
using namespace std;
 
int main()
{
    // Считываем размер вводимой матрицы
    int size;
    cin >> size;
    
    // Будем хранить матрицу в векторе, состоящем из 
    // векторов вещественных чисел
    vector <vector <long double> > matrix;
    
    // Матрица будет иметь размер (size) x (size + 1),
    // c учетом столбца свободных членов    
    matrix.resize (size);
    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;
    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]);
    }
 
    return 0;   
}
0
S_el
2113 / 1633 / 308
Регистрация: 15.12.2013
Сообщений: 6,577
28.02.2014, 21:27 #22
У вашей системы такое решение?
x1 = 1.1680182102717
x2 = 2.199459382558
x3 = 1.8651301749893
x4 = - 0.67292644757433
0
Ev_Hyper
Заблокирован
28.02.2014, 21:39 #23
S_el, на ваш вопрос уже был дан ответ.
Метод простой итерации
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
01.03.2014, 05:12  [ТС] #24
S_el, да да, вы как сделали?
1,1680
2,1995
1,8651
-0,6729
0
S_el
2113 / 1633 / 308
Регистрация: 15.12.2013
Сообщений: 6,577
01.03.2014, 10:40 #25
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
S_el, да да, вы как сделали?
1,1680
2,1995
1,8651
-0,6729
Я проверил обычным методом решения.
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
02.03.2014, 06:15  [ТС] #26
S_el, сможете пожалуйста код написать чтоб решал этот пример?
0
Ev_Hyper
Заблокирован
02.03.2014, 12:48 #27
Alexandr1966, я же вам уже написал, что нужно исправить, чтобы программа заработала.
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
02.03.2014, 14:15  [ТС] #28
Ev_Hyper, честно сказать я вообще не понимаЮ, помогите переделать(
0
Ev_Hyper
Заблокирован
02.03.2014, 16:22 #29
Цитата Сообщение от Alexandr1966 Посмотреть сообщение
Ev_Hyper, честно сказать я вообще не понимаЮ, помогите переделать(
Alexandr1966, помочь переделать или сделать за вас?
0
Alexandr1966
26 / 26 / 1
Регистрация: 17.12.2012
Сообщений: 438
02.03.2014, 17:17  [ТС] #30
Ev_Hyper, если есть возможность переделать этот код под тот который мне нужен, то помогите буду очень вам признателен. Если невозможно его переделать и нужно делать всё с нуля, то я не знаю, если вы захочете мне помочь сделать по новой всё, то цены вам не будет, вы будете просто моим спасителем)
0
02.03.2014, 17:17
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
02.03.2014, 17:17
Привет! Вот еще темы с ответами:

Метод простой итерации - C++
Добрый день! нужна помощь, помогите пожалуйста с решением задачи Решить методом простой итерации уравнение вида х=Ф(Х). Очевидное...

Метод простой итерации - C++
Найти решение уравнения x=f(x) методом Эйткена-Стеффенсона, в котором от заданного начального Xo три очередных приближения находятся по...

Метод простой итерации - C++
Помогите решить методом простой итерации уравнение вида x=f(x).Очередное приближение корня находится по формуле {x}_{n}=f({x}_{n-1});...

Метод простой итерации - C++
#include &lt;iostream&gt; #include &lt;cmath&gt; using namespace std; double f(double x) { return (-2.5*pow(x, 3) + 0.25*pow(x, 2) - 3); ...


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

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

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