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

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

Войти
Регистрация
Восстановить пароль
 
DenProx
Техник
310 / 168 / 5
Регистрация: 09.10.2009
Сообщений: 3,057
#1

Рунге-Кутта из 4 порядка в 3-й - C++

20.06.2013, 14:55. Просмотров 357. Ответов 0
Метки нет (Все метки)

Доброго времени суток. У меня есть такой кусок кода, для численного метода Рунге-Кутта 4 порядка:

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
///////////////////////////////////////////////////////////////////////
    //                  - Метод Рунге-Кутта 4 порядка -                  //
    ///////////////////////////////////////////////////////////////////////
    for (j = 0; j < n; j++) y0[j] = y[j];   // Копируем его в y0
    f (y0, q0, xi);         // Заполняем q0, основываясь на значениях из y0
    for (j = 0; j < n; j++) q0[j] *= h; // Делаем q0
    xi += h;                // Следующий шаг
    // ... а остальные 3 добываем с помощью метода Рунге-Кутта 4 порядка.
    for (i = 0; i < 3; i++)     // i - КАКОЕ ЗНАЧЕНИЕ УЖЕ ЕСТЬ
    {                   // А ВЫЧИСЛЯЕМ ЗНАЧЕНИЯ Y[i+1]!!!!
        // Сначала нужны коэффициенты k1
        // Элемент y[i, j] = y0 + (i * n) + j = y0[i * n + j]
        f (&y0[i * n], k1, xi); // Вычислим f(xi, yi) = k1 / h
        // И для каждого дифференциального уравнения системы проделываем
        // операции вычисления k1, а также подготовки в ya аргумента для
        // вычисления k2
        for (j = 0; j < n; j++)
        {
            k1[j] *= h;                 // Вычислим наконец-то k1
            ya[j] = y0[i*n+j] + k1[j] / 2.;
            // И один из аргументов для функции 
        }                           // вычисления k2
        f (ya, k2, xi + (h / 2.));          // Вычислим f(xi,yi) = k2 / h
        for (j = 0; j < n; j++)
        {                           // Вычислим наконец-то k2
            k2[j] *= h;
            ya[j] = y0[i*n+j] + k2[j] / 2.; // И один из аргументов для функции 
        }                           // вычисления k3
        f (ya, k3, xi + h / 2.);            // Вычислим f(xi,yi) = k3 / h
        for (j = 0; j < n; j++)
        {
            k3[j] *= h;                 // Вычислим наконец-то k3
            ya[j] = y0[i*n+j] + k3[j];  // И один из аргументов для функции 
        }                       // вычисления k4
        f (ya, k4, xi + h);         // Вычислим f(xi,yi) = k4 / h
        for (j = 0; j < n; j++) k4[j] *= h; // Вычислим наконец-то k4
        // Надо вычислить приращение каждой функции из n
        for (j = 0; j < n; j++)         // Вычисляем следующее значение
                                // функции
                                // Y[i+1] = Yi + ...
    y0[(i+1)*n+j] = y0[i*n+j] + (k1[j] + 2. * k2[j] + 2 * k3[j] + k4[j]) / 6.;
                                // И новое значение q[i+1]
        f (&y0[(i+1)*n], &q0[(i+1)*n], xi); // qi = f (xi, yi);
        for (j = 0; j < n; j++) q0[((i+1)*n)+j] *= h;
        xi += h;                    // Следующий шаг    }
Помогите пожалуйста переделать его на 3-й порядок. В численных методах я не силен, разбираться времени нет ... сроки поджимают. Думаю тут разница в 2 строки примерно должна быть, одну убрать, другую немного подредактировать.

Заранее спасибо.

Добавлено через 19 часов 4 минуты
Тема закрыта. Изменения не потребовались 3-го порядка требовалось уравнение, а метод можно и 4-го )
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
20.06.2013, 14:55
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Рунге-Кутта из 4 порядка в 3-й (C++):

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

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

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

Метод Рунге-Кутта 4 порядка - C++
Помогите найти ошибку в методе рунге-Кутта 4 порядка System::System(double m, const Vector3D&amp; g, double C, double f, double p, double S,...

Метод Рунге-Кутта 3 порядка - C++
Начерикал что-то,вроде бы работает,но не уверен,да и не доходит,как под условия сделать задачи,тыкните,пожалуйста #include &quot;stdafx.h&quot; ...

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

0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.06.2013, 14:55
Привет! Вот еще темы с ответами:

Метод рунге-кутта 3 порядка - C++
дана функция d(y(x))/dx=e^x-2y(x) Нач. условия y(0)=e Код#include&lt;stdio.h&gt; #include&lt;conio.h&gt; #include&lt;math.h&gt; float fun(float...

Метод Рунге-Кутта второго порядка - C++
В общем есть задание. Задание к-е на скрине: Ток в электрической цепи описывается дифференциальным уравнением (которое на скрине),...

Метод Рунге-Кутта четвертого порядка - C++
Доброго времени суток. Задание по выч. мату. координаты х рассчитываются просто по шагу, но у присваивается значение начальной точки и не...

Метод Рунге-Кутта 4 порядка (исправить код) - C++
Всем привет. Нужно решить пример методом Рунге-Кутта 4-го порядка точности. Пример: y=1/2*x*y x0=0 - начало интервала. y0=1 -...


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

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

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