Форум программистов, компьютерный форум, киберфорум
C++ Builder
Войти
Регистрация
Восстановить пароль
 
Рейтинг 5.00/6: Рейтинг темы: голосов - 6, средняя оценка - 5.00
artem2a
1

компьютерный код интегрирующий систему уравнений Гамильтона для модели Хенона-Хейлеса методом Рунге-Кутта 4-го порядка

22.05.2013, 20:48. Просмотров 1077. Ответов 1
Метки нет (Все метки)


Помогите, пожалуйста, построить график.
Задание:
Написать компьютерный код интегрирующий систему уравнений Гамильтона для модели Хенона-Хейлеса
(dx/dt)=px, (dy/dt)=py,
(dpx/dt)=-x-2xy, (dpy/dt)=y2-y-x2
методом Рунге-Кутта 4-го порядка.
Построить сечения Пуанкаре(график зависимости импульса от координаты) для энергий E=1/12, 1/8, 1/6.


ФАЙЛ Unit1.cpp
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
49
#include <math.h>
#include <iostream.h>
#include "runge_kutta.h"
//---------------------------------------------------------------------------
#pragma hdrstop
//---------------------------------------------------------------------------
#pragma argsused
double m = 1.0;
void rhs(double *y,double *dy,double t)
{
    dy[0] = y[2];
    dy[1] = y[3];
    dy[2] = -(y[0]+2*y[0]*y[1]);
    dy[3] = -(y[1]+y[0]*y[0]-y[1]*y[1]);
}
double EDS(double t)
{
return 0;
}
 
int main(int argc, char* argv[])
{
    int n = 4;
    double ht = 0.05;
    RungeKutta *rk = new RungeKutta(n,ht,rhs);
 
    double
        Etot=1/8,
        V=(1/12),
        w=(Etot-V)*2*m,
        *in = new double [n];
    in[0] = 0.1;
    in[1] = 0.2;
    in[2] = m*sqrt(w);
    in[3] = 0.3;
 
    rk->setInitialValues(in);
 
    double t = 0.1;
    for(int it=0;it<100;it++)
    {
        rk->doOneRKStep(t);
        t += ht;
 
        cout << "x=" << rk->y1[1] << "  y="<< rk->y1[3] << "\n";
    }
    return 1;
}
//---------------------------------------------------------------------------

ФАЙЛ Unit4.cpp
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
49
50
//---------------------------------------------------------------------------
#include <vcl.h>
#pragma hdrstop
 
#include "Unit4.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm4 *Form4;
 
//---------------------------------------------------------------------------
__fastcall TForm4::TForm4(TComponent* Owner)
        : TForm(Owner)
{
    DrawingBoard=new Graphics::TBitmap;
    DrawingBoard->Width = this->PaintBox1->Width;
    DrawingBoard->Height = this->PaintBox1->Height;
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1MouseDown(TObject *Sender,
      TMouseButton Button, TShiftState Shift, int X, int Y)
{
    IsDrawing = True;
    StartX = X;
    StartY = Y;
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1MouseUp(TObject *Sender,
      TMouseButton Button, TShiftState Shift, int X, int Y)
{
    if( IsDrawing )
    {
        DrawingBoard->Canvas->Rectangle(StartX, StartY, X, Y);
        PaintBox1->Canvas->Draw(0, 0, DrawingBoard);
        IsDrawing = False;
    }
}
//---------------------------------------------------------------------------
void __fastcall TForm4::PaintBox1Paint(TObject *Sender)
{
    PaintBox1->Canvas->Draw(0, 0, DrawingBoard);
}
//---------------------------------------------------------------------------
 
 
void __fastcall TForm4::Button1Click(TObject *Sender)
{
}
 
//---------------------------------------------------------------------------
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
22.05.2013, 20:48
Ответы с готовыми решениями:

Решение уравнений методом Рунге-Кутта 4 порядка
напишите програму для решения уравнений методом Рунги-кута 4 порядка))

Решение дифференциальных уравнений четвертого порядка методом Эйлера и методом Рунге-Кутта
Форумчане прошу помочь решить уравнение f :=y+y'*x+y''+y'''*x методами Рунге-Кутта и Эйлера. Если...

Решение системы 6 уравнений методом Рунге-Кутта 4 порядка
Подскажите пожалуйста, как переделать мою программу для решения системы 6 уравнений методом рунге...

Решение системы уравнений методом Рунге-Кутта 4-го порядка
Всем привет! Помогите составить код для решения данной системы уравнений.

1
Эксперт С++
8354 / 6101 / 610
Регистрация: 10.12.2010
Сообщений: 28,387
Записей в блоге: 29
23.05.2013, 16:32 2
Метод Рунге-Кутта : https://www.cyberforum.ru/post4579895.html
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
23.05.2013, 16:32

Заказываю контрольные, курсовые, дипломные и любые другие студенческие работы здесь или здесь.

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

Решение системы дифференциальных уравнений Методом Рунге-Кутта 4 порядка
В соседней теме я уже выкладывал решение системы дифф. уравнений через ode23....

Метод Рунге-Кутта 4 порядка для уравнений n-го порядка
Здравствуйте. Задание: Написать метод Рунге-Кутта 4 порядка для уравнений n-го порядка. Я...

Коши для ОДУ 1-го порядка методом Рунге–Кутта 4-го порядка
==Решить задачу Коши для ОДУ 1-го порядка методом Рунге–Кутта 4-го порядка. Количество узлов – 6....


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

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

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