0 / 0 / 0
Регистрация: 12.06.2015
Сообщений: 10
1

Метод Рунге-Кутта 4 порядка для решения системы дифференциальных уравнений

30.10.2015, 05:59. Показов 12174. Ответов 2
Метки нет (Все метки)

Имеется код программы, но по не очевидным для меня причинам работает не корректно, во время отладки выдает сообщения типа :

Run-Time Check Failure #2 - Stack around the variable 'k4' was corrupted.

Run-Time Check Failure #2 - Stack around the variable 'k3' was corrupted.


Подскажите в чем проблема.

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
104
#include <stdio.h>
 
#include <math.h>
#include <iostream>
 
using namespace std;
double f(int i, double x, double y[4]) {
    switch (i) {
    case 1:return y[4]; break;
    case 2: return y[1]; break;
    case 3:return y[2]; break;
    case 4:return y[3]; break;
    default: break;
    }
}
/*************************************************************************
Один шаг метода Рунге-Кутта четвертого порядка для решения
системы дифферециальных уравнений.
 
procedure SystemRungeKuttStep(
const x : Real;
const h : Real;
const n : Integer;
var y : array of Real);
 
Алгоритм совершает один шаг метода для системы
диффуров y[i]'=F(i,x,y) для i=1..n
 
Начальная точка имеет кординаты (x,y[1], ..., y[n])
 
После выполнения алгоритма в переменной y содержится состояние
системы в точке x+h
*************************************************************************/
void step(double x, double h, int n, double y[4])
{
    int i = 0;
    double yt[4];
    double k1[4];
    double k2[4];
    double k3[4];
    double k4[4];
 
    for (i = 1; i <= n; i++)
    {
        k1[i] = h*f(i, x, y);
    }
    for (i = 1; i <= n; i++)
    {
        yt[i] = y[i] + 0.5*k1[i];
    }
    for (i = 1; i <= n; i++)
    {
        k2[i] = h*f(i, x + h*0.5, yt);
    }
    for (i = 1; i <= n; i++)
    {
        yt[i] = y[i] + 0.5*k2[i];
    }
    for (i = 1; i <= n; i++)
    {
        k3[i] = h*f(i, x + h*0.5, yt);
    }
    for (i = 1; i <= n; i++)
    {
        yt[i] = y[i] + k3[i];
    }
    for (i = 1; i <= n; i++)
    {
        k4[i] = h*f(i, x + h, yt);
    }
    for (i = 1; i <= n; i++)
    {
        y[i] = y[i] + (k1[i] + 2.0*k2[i] + 2.0*k3[i] + k4[i]) / 6;
    }
}
/*************************************************************************
Алгоритм решает систему диффуров y[i]'=F(i,x,y) для i=1..n
методом Рунге-Кутта 4 порядка.
 
Начальная точка имеет кординаты (x,y[1], ..., y[n])
 
До конечной точки мы добираемся через n промежуточных
с постоянным шагом h=(x1-x)/m
 
Результат помещается в переменную result[4]
*************************************************************************/
void solvesystemrungekutta(double x, double x1, int steps, double result[4]) {
 
    for (int i = 1; i <= steps - 1; i++)
    {
        step(x + i*(x1 - x) / steps, (x1 - x) / steps, 4, result);
    }
}
 
int main() {
    //первая координата 
    double temporaryresult[4];
    solvesystemrungekutta(2, 3, 3, temporaryresult);
    for (int i = 1; i <= 4; i++) {
        printf("Htpekmnf", temporaryresult[i], '\n');
    }
    char c = getchar();
    return 0;
}
__________________
Помощь в написании контрольных, курсовых и дипломных работ, диссертаций здесь
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
30.10.2015, 05:59
Ответы с готовыми решениями:

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

Метод Рунге-Кутта для системы из 5 уравнений
Нужно написать программу для решения системы дифуров. Задано 5 уравнений, необходимо вывести...

Метод Рунге-Кутта для системы уравнений
Написал метод рунге-кутта для динамо риккитаке, но проблема в том, что данные мне формулы не...

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

2
Модератор
Эксперт С++
12107 / 9778 / 5910
Регистрация: 18.12.2011
Сообщений: 26,247
30.10.2015, 10:47 2
См. ссылки внизу страницы
0
0 / 0 / 1
Регистрация: 12.10.2015
Сообщений: 207
15.04.2016, 10:36 3
Какие ссылки?
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
15.04.2016, 10:36
Помогаю со студенческими работами здесь

Какое условие поставить на выходе из цикла в методе Рунге Кутта 4 порядка для системы ДУ 1 порядка?
Пока имеется вот такой код для решения системы ДУ 1 порядка, если имеется ошибка исправьте...

Метод Рунге-Кутта 2-го порядка
Написал код программы. С компиляцией вроде бы нет проблем. А цикл for воспринимать не хочет....

Метод Рунге-Кутта 4-го порядка
#include &lt;stdio.h&gt; #include &lt;iostream&gt; #include &lt;math.h&gt; using namespace std; double f(double...

Метод Рунге-Кутта 4-го порядка
Правильно ли я записал выражения для коэффициентов моего уравнения(система Лоренца) #include...


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

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

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