Форум программистов, компьютерный форум, киберфорум
Наши страницы
С++ для начинающих
Войти
Регистрация
Восстановить пароль
 
Teapot1488
3 / 3 / 1
Регистрация: 10.03.2013
Сообщений: 125
#1

Решение методом Рунге - Кутта. Переполнение переменной - C++

01.12.2016, 00:57. Просмотров 231. Ответов 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
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <windows.h>
#include <fstream>
#include <cstdio>
#include <iostream>
#include <dos.h>
#include <Math.h>
 
using namespace std;
 
 
void set_max_size()
{
    HANDLE out_handle = GetStdHandle(STD_OUTPUT_HANDLE);
    COORD maxWindow = GetLargestConsoleWindowSize(out_handle); // размер самого большого возможного консольного окна
    SMALL_RECT srctWindow = { 0, 0, maxWindow.X - 5, maxWindow.Y - 5 };
    SMALL_RECT minWindow = { 0, 0, 0, 0 };
    SetConsoleWindowInfo(out_handle, true, &minWindow);
    SetConsoleScreenBufferSize(out_handle, maxWindow);
    SetConsoleWindowInfo(out_handle, true, &srctWindow);
    cout << "max console size X = " << maxWindow.X << "; " << "max console size Y = " << maxWindow.Y <<endl;
}
 
int main()
{
 
    set_max_size();     // размер окна
    system("color F0"); // цвет
 
    double tn, tk, yn, h, Fe, xe, ye, Fem, xem, yem, Fr, xr, yr, k1, k2, k3, k4;
    int i, j, k, n;
    cout << "Differential equitation: t' = y^2*t " <<endl;
    cout << "tn=";
    cin >> tn;
    cout << "tk=";
    cin >> tk;
    cout << "yn=";
    cin >> yn;
    cout << "h=";
    cin >> h;
    n=static_cast<int>((tk-tn)/h)+1;
    cout << "Quantity of ierations n = ";
    cout << n << endl;
    getch();
    //Eiler method
    xe=tn; ye=yn;
    cout << "Eiler method" <<endl;
    for (i = 0; i < n; i++)
    {
    Fe= Fe+h*(ye*ye*xe);
    xe=xe+h;
    cout << Fe << endl;
    }
    getch();
    //Eiler modified
    xem=tn; yem=yn;
    cout << "Eiler modified method" <<endl;
    for (j = 0; j < n; j++)
    {
    Fem= Fem+(h/2)*(yem*yem*xem+(yem+h)*(yem+h)*(xem+h));
    cout << Fem << endl;
    xem=xem+h;
    }
    getch();
    //Runge-Kuta
    xr=tn; yr=yn;
    cout << "Runge - Kuta method" <<endl;
    for (k = 0; k < n; k++)
    {
    k1=h*(xr*yr*yr);
    k2=h*((xr+h/2)*(yr+k1/2)*(yr+k1/2));
    k3=h*((xr+h/2)*(yr+k2/2)*(yr+k2/2));
    k4=h*((xr+h)*(yr+k3)*(yr+k3));
    yr=yr+(k1+2*k2+3*k3+k4)/6;
    cout << yr << endl;
    xr=xr+h;
    }
 
getch();
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
01.12.2016, 00:57
Я подобрал для вас темы с готовыми решениями и ответами на вопрос Решение методом Рунге - Кутта. Переполнение переменной (C++):

Решение задачи Коши методом Рунге-Кутта
Решение задачи Коши методом Рунге-Кутта и сравнить результат с MATLAB не...

Решение дифференциального уравнения методом Рунге-Кутта
решить уравнение v'=u/m*m'. обе производные меняются по времени масса...

Решение системы диффуров методом Рунге-Кутта
#include &lt;stdio.h&gt; #include &lt;stdafx.h&gt; #include &lt;math.h&gt; #include &lt;iostream&gt;...

Решение уравнений методами Рунге-Кутта и Эйлера
Решить уравнение при помощи методов Эйлера и Рунге Кутта, описывающее смещение...

Численно решить ОДУ методом Рунге-Кутта 4 порядка
Доброго времени суток. Мне нужно численно решить ОДУ методом Рунге-Кутта 4...

Решение обычных деференциальных уравнений методами Ейлера, Рунге-Кутта
Доброго времени суток!Сдаю академ разницу в университете и вот по информатике...

3
RJM
3 / 3 / 2
Регистрация: 21.09.2016
Сообщений: 83
01.12.2016, 02:30 #2
Это консольное приложение? По формулам ничего не скажу, но по коду:

- в конце int main () у тебя метод ничего не возвращает, хотя он описан как int. Соответственно и виснет все.
- в конце кода какой-то огрызок с getch();}, зачем?
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
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
#include <windows.h>
#include <fstream>
#include <cstdio>
#include <iostream>
#include <dos.h>
#include <Math.h>
 
using namespace std;
 
 
void set_max_size()
{
    HANDLE out_handle = GetStdHandle(STD_OUTPUT_HANDLE);
    COORD maxWindow = GetLargestConsoleWindowSize(out_handle); // размер самого большого возможного консольного окна
    SMALL_RECT srctWindow = { 0, 0, maxWindow.X - 5, maxWindow.Y - 5 };
    SMALL_RECT minWindow = { 0, 0, 0, 0 };
    SetConsoleWindowInfo(out_handle, true, &minWindow);
    SetConsoleScreenBufferSize(out_handle, maxWindow);
    SetConsoleWindowInfo(out_handle, true, &srctWindow);
    cout << "max console size X = " << maxWindow.X << "; " << "max console size Y = " << maxWindow.Y <<endl;
}
 
int main()
{
 
    set_max_size();     // размер окна
    system("color F0"); // цвет
 
    double tn, tk, yn, h, Fe, xe, ye, Fem, xem, yem, Fr, xr, yr, k1, k2, k3, k4;
    int i, j, k, n;
    cout << "Differential equitation: t' = y^2*t " <<endl;
    cout << "tn=";
    cin >> tn;
    cout << "tk=";
    cin >> tk;
    cout << "yn=";
    cin >> yn;
    cout << "h=";
    cin >> h;
    n=static_cast<int>((tk-tn)/h)+1;
    cout << "Quantity of ierations n = ";
    cout << n << endl;
    getch();
    //Eiler method
    xe=tn; ye=yn;
    cout << "Eiler method" <<endl;
    for (i = 0; i < n; i++)
    {
    Fe= Fe+h*(ye*ye*xe);
    xe=xe+h;
    cout << Fe << endl;
    }
    getch();
    //Eiler modified
    xem=tn; yem=yn;
    cout << "Eiler modified method" <<endl;
    for (j = 0; j < n; j++)
    {
    Fem= Fem+(h/2)*(yem*yem*xem+(yem+h)*(yem+h)*(xem+h));
    cout << Fem << endl;
    xem=xem+h;
    }
    getch();
    //Runge-Kuta
    xr=tn; yr=yn;
    cout << "Runge - Kuta method" <<endl;
    for (k = 0; k < n; k++)
    {
    k1=h*(xr*yr*yr);
    k2=h*((xr+h/2)*(yr+k1/2)*(yr+k1/2));
    k3=h*((xr+h/2)*(yr+k2/2)*(yr+k2/2));
    k4=h*((xr+h)*(yr+k3)*(yr+k3));
    yr=yr+(k1+2*k2+3*k3+k4)/6;
    cout << yr << endl;
    xr=xr+h;
    system("pause");
    return 0;
    }
Пробуй теперь.
0
Teapot1488
3 / 3 / 1
Регистрация: 10.03.2013
Сообщений: 125
01.12.2016, 03:31  [ТС] #3
в вашем коде всё также работает как и у меня. Один в один. С той лишь разницей, что вы перенесли в цикл паузу.
Программа вылетает из-за того, что переменная слишком быстро растет и вылазит за пределы 1,8*10^308
0
RJM
3 / 3 / 2
Регистрация: 21.09.2016
Сообщений: 83
01.12.2016, 05:36 #4
у тебя 2 переменные (Fem и Fe) не инициализированы. 53 и 63 строка.
Попробовал у себя в ВС запустить.

Добавлено через 3 минуты
Атакаи рабочий код скинул.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
01.12.2016, 05:36
Привет! Вот еще темы с решениями:

Численное интегрирование системы дифференциальных уравнений методом Рунге - Кутта
Доброго времени суток. Пытаюсь разработать функцию для численного...

Функция для численного интегрирования системы дифференциальных уравнений методом Рунге - Кутта
Доброго времени суток. У меня есть вот такая задачка: Разработать функцию...

Решение дифференциального уравнения методом Рунге-Кутты
Не могу понять почему не правильно выдает значения. #include &quot;stdafx.h&quot;...

Решение задачи Коши методом Эйлера и Рунге-Кутты
Доброго времени суток всем:) Писал прогу для решения задачи Коши методом...


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

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

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