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

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

Войти
Регистрация
Восстановить пароль
 
uZZZZer
0 / 0 / 0
Регистрация: 17.05.2015
Сообщений: 3
#1

Рунге-Кутта 4го порядка из книги (The Art of Scientific Computing) - C++

06.12.2015, 15:48. Просмотров 258. Ответов 1
Метки нет (Все метки)

Доброго всем дня!
Нашел в сети книгу с исходниками, если кому нужно могу скинуть.
Press W.H., Teukolsky S.A., Vetterling W.T., Flannery B.P. - Numerical Recipes The Art of Scientific Computing - 2007

Но вопрос в следующем, в общем нужно мне стало решить уравнение Ван-дер-Поля методом Рунге-Кутта 4го порядка, сам проект с примером легко без проблем компилировался,запустился под Visual C++ 2015, подставил свое уравнение, все работает, Но! решение очень быстро расходится с моделью в Mathcad, начиная с 6го шага. Правильно ли я все сделал? Подскажите.
Прикладываю архив с проектом и маткадовской моделькой.
 Комментарий модератора 
П.5.18.Правил
Запрещено размещать задания и решения в виде картинок и других файлов с их текстом.
0
Вложения
Тип файла: rar c.rar (753.1 Кб, 3 просмотров)
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
06.12.2015, 15:48
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Рунге-Кутта 4го порядка из книги (The Art of Scientific Computing) (C++):

Решение диф. уравнения для колебаний методом Рунге—Кутты 4го порядка.(С++) - C++
Ребят помогите пожалуйста решить такое на C++ уравнение x''=(-w^2)*x где w-const там замену нужно сделать что то вроде y=x' и y'=x'' и...

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

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

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

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

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

1
uZZZZer
0 / 0 / 0
Регистрация: 17.05.2015
Сообщений: 3
06.12.2015, 23:33  [ТС] #2
А в чем я был не прав, я про оформление?
 Комментарий модератора 
Выложите текст прикрепленного файла в сообщение


Добавлено через 53 минуты
Кликните здесь для просмотра всего текста
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
/* Driver for routine rk4 */
 
#include <stdio.h>
#define NRANSI
#include "nr.h"
#include "nrutil.h"
 
#define N 2
 
void derivs(float x,float y[],float dydx[])
{
    dydx[1] = y[2];
    dydx[2]=(1-y[1]* y[1])*y[2]-y[1];
    //dydx[3]=y[2]-(2.0/x)*y[3];
    //dydx[4]=y[3]-(3.0/x)*y[4];
}
 
int main(void)
{
    int i,j;
    float h, x = 0.00, *y, *dydx, *yout;
 
    y=vector(1,N);
    dydx=vector(1,N);
    yout=vector(1,N);
    y[1] = 0.01;//bessj0(x);
    y[2] = 0.00;//bessj1(x);
    //y[3]=bessj(2,x);
    //y[4]=bessj(3,x);
    derivs(x,y,dydx);
    //printf("\n%16s %5s %12s %12s %12s\n",
        //"Bessel function:","j0","j1","j3","j4");
    for (i=1;i<=130;i++) {
        h = 0.1*i;
        //x = x + h;
        rk4(y,dydx,N,x,h,yout,derivs);
        printf("\nfor a step size of: %6.2f\n",h);
        printf("%12s","rk4:");
        for (j=1;j<=2;j++) printf(" %12.6f",yout[j]);
        //printf("\n%12s %12.6f %12.6f %12.6f %12.6f\n","actual:",
            //bessj0(x+h),bessj1(x+h),bessj(2,x+h),bessj(3,x+h));
    }
    free_vector(yout,1,N);
    free_vector(dydx,1,N);
    free_vector(y,1,N);
    return 0;
}
#undef NRANSI


Добавлено через 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
#define NRANSI
#include "nrutil.h"
 
void rk4(float y[], float dydx[], int n, float x, float h, float yout[],
    void (*derivs)(float, float [], float []))
{
    int i;
    float xh,hh,h6,*dym,*dyt,*yt;
 
    dym=vector(1,n);
    dyt=vector(1,n);
    yt=vector(1,n);
    hh=h*0.5;
    h6=h/6.0;
    xh=x+hh;
    for (i=1;i<=n;i++) yt[i]=y[i]+hh*dydx[i];
    (*derivs)(xh,yt,dyt);
    for (i=1;i<=n;i++) yt[i]=y[i]+hh*dyt[i];
    (*derivs)(xh,yt,dym);
    for (i=1;i<=n;i++) {
        yt[i]=y[i]+h*dym[i];
        dym[i] += dyt[i];
    }
    (*derivs)(x+h,yt,dyt);
    for (i=1;i<=n;i++)
        yout[i]=y[i]+h6*(dydx[i]+dyt[i]+2.0*dym[i]);
    free_vector(yt,1,n);
    free_vector(dyt,1,n);
    free_vector(dym,1,n);
}
#undef NRANSI


Добавлено через 6 часов 19 минут
Помогите, а то я себе голову всю сломал!
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
06.12.2015, 23:33
Привет! Вот еще темы с ответами:

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

Метод Рунге-Кутта 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) { ...

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

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


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

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

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