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

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

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 157, средняя оценка - 4.62
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
#1

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

23.01.2011, 23:51. Просмотров 22907. Ответов 9
Метки нет (Все метки)

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)
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.01.2011, 23:51
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Метод Рунге-Кутта 4-го порядка (C++):

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

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

Метод рунге-кутта 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 воспринимать не хочет. Программа просто его пропускает. Помогите...

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

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

9
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:00 #2
может потому что n - целочисленное? Там автоматически округляется в 23 строке
0
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
24.01.2011, 00:09  [ТС] #3
ну так я (2+1)/0.3 = 10 колличество разбиений целое число
0
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:31 #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 строчках

Я первокур, метода не знаю, потому фиг его знает, как правильно)
0
cap
1 / 1 / 0
Регистрация: 27.10.2010
Сообщений: 51
24.01.2011, 00:35  [ТС] #5
ну да ; пропустил а так ошибки нет просто алгоритмов рунге кутта 4го порядка несколько видов я использовал этот
0
ser4ega
27 / 27 / 3
Регистрация: 15.11.2009
Сообщений: 143
24.01.2011, 00:55 #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");
}
то расхожений опять почти нет, видимо дело в ом, что при малых икс погрешности получаются довольно весомыми по сравнению с иксом и игриком, как то так)
короче методы все виноваты
0
kNOpk@
2 / 2 / 0
Регистрация: 10.01.2011
Сообщений: 9
24.01.2011, 03:09 #7
писала лабу по этой теме, может посмотри, может поможет. сам пример

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

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

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

Спасибо!
 Комментарий модератора 
Обращайтесь в раздел фриланс http://www.cyberforum.ru/freelance/
0
20.05.2014, 18:56
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
20.05.2014, 18:56
Привет! Вот еще темы с ответами:

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

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

Метод Рунге-Кутта пятого порядка, модификация Мерсона - C++
помогите сделать этот метод при условии - автоматический выбор шага по контрольным членами((( заранее спасибо!!

Метод Рунге-Кутта 4 порядка для решения системы дифференциальных уравнений - C++
Имеется код программы, но по не очевидным для меня причинам работает не корректно, во время отладки выдает сообщения типа : Run-Time...


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

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

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