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

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

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 40, средняя оценка - 4.85
Slayerz
2 / 2 / 1
Регистрация: 14.07.2012
Сообщений: 26
Завершенные тесты: 1
#1

Метод Рунге-Кутта - C++

14.07.2012, 11:57. Просмотров 5267. Ответов 9
Метки нет (Все метки)

Всем добрый день!
Задали практическую работу, в которой необходимо проинтегрировать функцию методом Рунге-Кутта 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
    float func= 0, F = 0, P = 0, U = 0, t = 0, e = 0, lambda = 0, gamma = 0; // Переменные для функции
    float Runge, K1 = 0, K2 = 0, K3 = 0, K4 = 0; //Переменные для Рунге-Кутта 4го порядка
 
    cout << "Insert p: ";
    cin >> P;
    cout << "\n";
    cout << "Insert U: ";
    cin >> U;
    cout << "\n";
    cout << "Insert F: ";
    cin >> F;
    cout << "\n";
    cout << "Insert e: ";
    cin >> e;
    cout << "\n";
    cout << "Insert gamma: ";
    cin >> gamma;
    cout << "\n";
    cout << "Insert lmbda: ";
    cin >> lambda;
    cout << "\n";
    cout << "Insert t: ";
    cin >> t;
    cout << "\n";
 
    func = (F*(sqrt(P/U))) * (sin(gamma)*cos(lambda) + ((e*pow(cos(gamma), 2) + 2*cos(gamma) + e)/(1 + cos(gamma)))*sin(lambda)); // Функция
 
    K1 = func + t*func;
    cout << "Answer for K1: " << K1 << endl; //Элементы из формулы Рунге-Кутта, их вычисление
    K2 = func + 0.5*K1;
    cout << "Answer for K2: " << K2 << endl;
    K3 = func + 0.5*K2;
    cout << "Answer for K3: " << K3 << endl;
    K4 = func + K3;
    cout << "Answer for K4: " << K4 << endl;
 
    Runge = func + (K1 + K2 + K3 + K4)/6; // Вычисление Рунге-Кутта
    cout << "Answer for Runge-Kutt(4th): ";
    cout << Runge << endl;
 
    system("pause");
    return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
14.07.2012, 11:57
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод Рунге-Кутта (C++):

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

Метод Рунге-Кутта - C++
как реализовать программу для уравнения вида dv/dt=const*dm/dt методом рунге-кутта. как именно заменить производную справа

метод Рунге-Кутта - C++
Всем доброго времени суток, меня интересует, пользуясь методом Рунге -Кута для диф уравнений, я получу в ответе точки, графика диф...

Метод Рунге-Кутта - C++
Нужно решить диффур второго порядка. При решении выводятся нуля. Где я оплошал? #include &lt;math.h&gt; using namespace std; ...

Метод Рунге-Кутта - C++
Помогите написать метод Рунге Кутта для решение систем диф. уравнений из 4 уравнений. dx/dt=Px=f1(t,x,y,Px,Py) dy/dt=Py=f2(t,x,y,Px,Py)...

метод Рунге-Кутта - C++
Помогите написать программу на хаотичное движение бильярдных шаров методом Рунге-Кутта. При столкновении шары должны отталкиваться, ну и...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Catstail
Модератор
22540 / 10945 / 1775
Регистрация: 12.02.2012
Сообщений: 18,083
14.07.2012, 14:25 #2
На глаз - все верно, кроме одного: метод Рунге-Кутта - это не метод интегрирования функции, а метод решения обыкновенного дифференциального уравнения.

Да, написано только тело цикла, а где переход к следующей точке?
Slayerz
2 / 2 / 1
Регистрация: 14.07.2012
Сообщений: 26
Завершенные тесты: 1
14.07.2012, 14:41  [ТС] #3
До этого я еще не успел дойти, так как запутался, как корректнее это реализовать.
Как указано в задании, что нужно выбрать интервал интервал точек от 60 до 600.
То есть я должен запустить цикл для t(это переменная для интервала, как я понял), к примеру, с шагом в 70 и поместить в этот цикл все формулы Рунге-Кутта, чтобы цикл пробежался и выдал все возможные ответы. Вот пример как я написал цикл для точек.

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
...........................................................
cout << "Insert t: ";
    cin >> t;
    for(int i = 0; i < 10; i++)
    {       
        t = t + 70;
        cout << "\n";
 
        func = (F*(sqrt(P/U))) * (sin(gamma)*cos(lambda) + ((e*pow(cos(gamma), 2) + 2*cos(gamma) + e)/(1 + cos(gamma)))*sin(lambda)); // Функция
 
        K1 = func + t*func;
        cout << "Answer for K1: " << K1 << endl; //Элементы из формулы Рунге-Кутта, их вычисление
        K2 = func + 0.5*K1;
        cout << "Answer for K2: " << K2 << endl;
        K3 = func + 0.5*K2;
        cout << "Answer for K3: " << K3 << endl;
        K4 = func + K3;
        cout << "Answer for K4: " << K4 << endl;
 
        Runge = func + (K1 + K2 + K3 + K4)/6; // Вычисление Рунге-Кутта
        cout << "Answer for Runge-Kutt(4th): ";
        cout << Runge << endl;
    }
.........................................................
Catstail
Модератор
22540 / 10945 / 1775
Регистрация: 12.02.2012
Сообщений: 18,083
14.07.2012, 14:51 #4
Странно то, что твоя "функция" не зависит от t... Или я что-то проглядел?
Slayerz
2 / 2 / 1
Регистрация: 14.07.2012
Сообщений: 26
Завершенные тесты: 1
14.07.2012, 15:28  [ТС] #5
В самой функции t нет, благодаря Рунге-Кутта я ищу df/dt.
Сама t используется в 1ом шаге Рунге-Кутта: K1 = func + t*func

Добавлено через 25 минут
Если, я конечно правильно понял принцип работы данного метода*
Catstail
Модератор
22540 / 10945 / 1775
Регистрация: 12.02.2012
Сообщений: 18,083
14.07.2012, 15:51 #6
"Если, я конечно правильно понял принцип работы данного метода" - похоже, что не совсем понял...
Slayerz
2 / 2 / 1
Регистрация: 14.07.2012
Сообщений: 26
Завершенные тесты: 1
14.07.2012, 22:17  [ТС] #7
Почитав еще теории и поискав примеры в интернете, пришел к выводу, что t является h(шагом интегрирования) просто в моем случае, на оси x откладывается время. То есть h = t, потому h я заменил на t. Вот, вроде вот так.
Catstail
Модератор
22540 / 10945 / 1775
Регистрация: 12.02.2012
Сообщений: 18,083
15.07.2012, 08:04 #8
Дело не в том, как обозначать нез. переменную, а в том, что Ваша "функция" от аргумента не зависит, т.е. F=const
Slayerz
2 / 2 / 1
Регистрация: 14.07.2012
Сообщений: 26
Завершенные тесты: 1
22.07.2012, 17:24  [ТС] #9
Я попробовал переписать код через перегрузку оператора, по совету преподавателя, ибо в уравнении будет 6 переменных, и метод Рунге тоже, но теперь происходит ошибка компиляции. Решил начать с простенького уравнения, чтобы как следует понять метод.
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
float step = 0, t;
    float x1;
    float x;
 
    cout << "Insert Integrator's step: ";
    cin >> t;
    cout << endl;
    cout << "Insert x: ";
    cin >> x;
    cout << endl;
 
    for(int i = 0; i < 10; i++)
    {       
        step = step + t;
        cout << "\n";
        float f(float x);
        {
            x1 = 2*x + 4;
        }
 
        K1 = step*f(x1);
        K2 = step*f(x1+K1/2);
        K3 = step*f(x1+K2/2);
        K4 = step*f(x1 + K3);
 
        Runge = f(x) + (1/6)*(K1 + 2*K2 + 2*K3 + K4); // Вычисление Рунге-Кутта
        cout << "Answer for Runge-Kutt(4th): ";
        cout << Runge << endl;
    }
    system("pause");
    return 0;
Jupiter
22.07.2012, 17:35     Метод Рунге-Кутта
  #10

Не по теме:

Цитата Сообщение от Slayerz Посмотреть сообщение
(1/6)
целочисленное деление

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

Метод рунге-кутта 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...

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

ДУ. Метод Рунге-Кутта. Погрешность. - C++
Сам Рунге-Кутта III порядка. - &quot;реализовано&quot;. Огорчает погрешность, не могу понять как улучшить точность (и не достаточно ли такой?). ...

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


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

Или воспользуйтесь поиском по форуму:
Yandex
Объявления
22.07.2012, 17:35
Ответ Создать тему
Опции темы

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