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

Метод Рунге-Кутта 4-го порядка - C++

Восстановить пароль Регистрация
 
Рейтинг: Рейтинг темы: голосов - 157, средняя оценка - 4.62
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
23.01.2011, 23:51     Метод Рунге-Кутта 4-го порядка #1
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
#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f(double x, double y)
{
      return(2*y/x+x*x*pow(2.7,x));               
}
 
int main(void){
   int n,i;
   double x,y,h,end,k1,k2,k3,k4;
   
   x=-1;
   
   y=-2.33
   
   end=2
   
   h=0.3;   
 
   n=(end-x)/h;
   
   for(i=1;i<=n;i++){
      k1=h*f(x, y);
      k2=h*f(x+h/4, y+k1/4);
      k3=h*f(x+0.5*h, y+0.5*k2);
      k4=h*f(x+h, y+k1-2*k2+2*k3);
      x+=h;
      y+=(k1+4*k3+k4)/6;
      cout<<"x="<<x<<" "<<"y="<<y<<endl;
   }
 
}

почему когда я ставлю h=0,3 мы идем по x перескакивая через 0...значит деления на ноль не происходит но результат все равно неверен когда ставим интервал от 1 до 2 тоесть ноль вобще не затрагиваем все нормально

вот какие должны быть результаты при x=2 y=18.36



аналитическое решение для проверки Y(x)=x*x(pow(2.7,x)-2.7)
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.01.2011, 23:51     Метод Рунге-Кутта 4-го порядка
Посмотрите здесь:

C++ Метод Рунге-Кутта пятого порядка, модификация Мерсона
метод Рунге-Кутта C++
Метод Рунге-Кутта C++
C++ Метод Рунге-Кутта 3 порядка
C++ Метод Рунге-Кутта 4 порядка (исправить код)
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:00     Метод Рунге-Кутта 4-го порядка #2
может потому что n - целочисленное? Там автоматически округляется в 23 строке
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
24.01.2011, 00:09  [ТС]     Метод Рунге-Кутта 4-го порядка #3
ну так я (2+1)/0.3 = 10 колличество разбиений целое число
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:31     Метод Рунге-Кутта 4-го порядка #4
это, а в алгоритме в27, 28 строчках нет ошибки? Может там так
C++
1
2
 k2=h*f(x+h/2, y+(h/2)/k1);
      k3=h*f(x+0.5*h, y+0.5*k2*h);
ну и так далее. я на вики слазил)
да и еще пропущены ; в 17 и 19 строчках

Я первокур, метода не знаю, потому фиг его знает, как правильно)
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
24.01.2011, 00:35  [ТС]     Метод Рунге-Кутта 4-го порядка #5
ну да ; пропустил а так ошибки нет просто алгоритмов рунге кутта 4го порядка несколько видов я использовал этот
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:55     Метод Рунге-Кутта 4-го порядка #6
я хз. Там вначале около точки -1 расхождениЯ при использовании алгоритма программы и проверочного - расхождения небольшие, но постепенно становЯтся все больше. Наверно ошибка метода.Если ставить вот так
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
#include <stdio.h>
#include <iostream>
#include <math.h>
using namespace std;
 
double f(double x, double y)
{
          return(2*y/x+x*x*pow(2.7,x));               
}
 
int main(void){
   int n,i;
   double x,y,h,end,k1,k2,k3,k4;
   
   x=1.1;
   
   y=0.341;
   
   end=2;
   
   h=0.3;   
 
   //n=(end-x)/h;
   
   for(i=1;i<=4;i++){
      k1=h*f(x, y);
      k2=h*f(x+h/4, y+k1/4);
      k3=h*f(x+0.5*h, y+0.5*k2);
      k4=h*f(x+h, y+k1-2*k2+2*k3);
      x+=h;
      y+=(k1+4*k3+k4)/6;
      cout<<"x="<<x<<" "<<"y="<<y<<endl;
      cout<<x*x*(pow(2.7,x)-2.7)<<endl;
   }
 system("PAUSE");
}
то расхожений опять почти нет, видимо дело в ом, что при малых икс погрешности получаются довольно весомыми по сравнению с иксом и игриком, как то так)
короче методы все виноваты
kNOpk@
2 / 2 / 0
Регистрация: 10.01.2011
Сообщений: 9
24.01.2011, 03:09     Метод Рунге-Кутта 4-го порядка #7
писала лабу по этой теме, может посмотри, может поможет. сам пример

y'=y+(1+x)*pow(y,2). y(0)=1 h=0.1; a=0; b=0.5.

4 знака после запятой.
Вложения
Тип файла: doc ЛР3.doc (88.0 Кб, 1176 просмотров)
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
24.01.2011, 16:51  [ТС]     Метод Рунге-Кутта 4-го порядка #8
делол в том что у меня в моей основной функции (2*y/x+x*x*pow(2.7,x)) есть 2y/x тоесть если я выберу шаг 0.1 то будет деление на ноль ..выводит неверные ответы это естественно...но вот почему когда беру шаг 0.3 x то тогда не равен 0 .. но выводит теже ответ как и как будто происходит деление на ноль
Гульзана
Сообщений: n/a
12.04.2012, 07:15     Метод Рунге-Кутта 4-го порядка #9
можете показать в этом примере как найти погрешность или он его вообще не ищет? в примере откуда взята k, end, и h???

Добавлено через 42 минуты
пожалуйста объясните что такое end = 2 в этой программе! никак не могу понять!!!!!!!!
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.05.2014, 18:56     Метод Рунге-Кутта 4-го порядка
Еще ссылки по теме:

Метод рунге-кутта 3 порядка C++
Метод Рунге-Кутта 2-го порядка C++
C++ Метод Рунге-Кутта четвертого порядка

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

Или воспользуйтесь поиском по форуму:
Tamara14
0 / 0 / 0
Регистрация: 20.05.2014
Сообщений: 3
20.05.2014, 18:56     Метод Рунге-Кутта 4-го порядка #10
Здравствуйте.
Требуется помощь в написании программы.
нужно проинтегрировать КА с солнечным парусом методом Рунге Кутта 4го порядка
мат модель есть .
И еще одна программа уже написанная ,нужно только проверить))
Ответьте пожалуйста на e-mail!
Я там все вышлю ,и договоримся об оплате!

Спасибо!
 Комментарий модератора 
Обращайтесь в раздел фриланс http://www.cyberforum.ru/freelance/
Yandex
Объявления
20.05.2014, 18:56     Метод Рунге-Кутта 4-го порядка
Ответ Создать тему
Опции темы

Текущее время: 15:35. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru