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

Метод Рунге-Кутта и многое другое

23.12.2012, 23:49. Просмотров 1239. Ответов 1
Метки нет (Все метки)


У меня есть задание:
 Комментарий модератора 
Правила:
5.18 Запрещено размещать задания в виде картинок и других файлов с их текстом.
надо проверить правильно ли я сделал ее?
Также нужна помощь с графиками, которые написанны в задании
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
51
52
53
54
55
56
57
58
59
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
 
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
 double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i;
double h = 0.1, k1, k2, k3, k4, m1, m2, m3, m4;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series1->AddXY(u,x);
   }
      Edit8->Text = x;
for(i = 0; i<N-1; i++)
{
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   Series2->AddXY(u,v);
}
      Edit9->Text = v;
}
Так же я выложил сам исходник
0
Вложения
Тип файла: rar RungeKutta.rar (466.6 Кб, 39 просмотров)
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
23.12.2012, 23:49
Ответы с готовыми решениями:

Метод Рунге-Кутта четвертого порядка (метод прогноза и коррекции)
Дано ОДУ x'''-x''x+(x')^2=0. Пытаюсь решить ОДУ методом Рунге-Кутта по данному примеру:...

Метод Эйлера и метод Рунге-Кутта: проверить код
Доброго времени суток. Хотел бы обратится к вам за помощью. Я написал программку, которая решает...

Графики - Метод Рунге-Кутта и Метод Адамса
Вот задание:

Графики - Метод Рунге-Кутта и Метод Адамса
Пожалуйста помогите!!! очень срочно нужно сделать графики (два в одном): метод Рунге-Кутта и метод...

1
6 / 6 / 5
Регистрация: 23.12.2012
Сообщений: 112
24.12.2012, 01:37  [ТС] 2
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
51
52
53
54
55
56
57
58
59
//---------------------------------------------------------------------------
 
#include <vcl.h>
#include <iostream.h>
#include <conio.h>
#include <math.h>
#pragma hdrstop
#include "Unit1.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TForm1 *Form1;
//---------------------------------------------------------------------------
__fastcall TForm1::TForm1(TComponent* Owner)
        : TForm(Owner)
{
}
//---------------------------------------------------------------------------
double f(double u, double a, double s, double x, double w, double d)
{
  return a*cos(s)+2*s*u+w*w*x;
}
 double f1(double u, double v)
{
  return u;
}
 
void __fastcall TForm1::Button1Click(TObject *Sender)
{
double a = Edit1->Text.ToDouble();
double u = Edit2->Text.ToDouble();
double s = Edit3->Text.ToDouble();
double x = Edit4->Text.ToDouble();
double w = Edit5->Text.ToDouble();
double d = Edit6->Text.ToDouble();
double v = Edit4->Text.ToDouble();
int N = Edit7->Text.ToDouble(), i, E=5;
double k1, k2, k3, k4, m1, m2, m3, m4;
double h = 0.1;
for(i = 0; i<N-1; i++)
   {
                   k1 = h*f(a, u, s, x, w, d);
                   k2 = h*f(u + h/2, x + (k1)/ 2, a, w, s, d);
                   k3 = h*f(u + h/2, x + (k2)/ 2, a, w, s, d);
                   k4 = h*f(u + h, x + k3, a, w, s, d);
                   x+= (k1 + 2*k2 + 2*k3 + k4)/6;
                   u+= h;
                   m1 = h*f1(u,v);
                   m2 = h*f1(u + h/2, v + (k1)/ 2);
                   m3 = h*f1(u + h/2, v + (k2)/ 2);
                   m4 = h*f1(u + h, v + k3);
                   v+= (m1 + 2*m2 + 2*m3 + m4)/6;
                   u+= h;
                   double T = f(u,a,s,x,w,d);
                   Series1->AddXY(T,x+v);
   }
      Edit8->Text=printf("%1.2f \n", x);
      Edit9->Text = v;
}
Программу улучшил, теперь нужно как то построить гармоническое и затухающие колебание, помогите пожалуйста
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
24.12.2012, 01:37

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

Метод Рунге-Кутта
Помогите, пожалуйста, с задачей! :( Разработать схему алгоритмов и программу численного решения...

Метод Рунге-Кутта
Здравствуйте, помогите пожалуйста решить проблему. Решаю простой дифур dy=-2*x*y Создаю m-файл...

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

Метод Рунге-Кутта
Сдесь искал, недопонял.


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

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

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