6 / 6 / 5
Регистрация: 23.12.2012
Сообщений: 112

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

23.12.2012, 23:49. Показов 1592. Ответов 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;
}
Так же я выложил сам исходник
Вложения
Тип файла: rar RungeKutta.rar (466.6 Кб, 42 просмотров)
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
23.12.2012, 23:49
Ответы с готовыми решениями:

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

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

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

1
6 / 6 / 5
Регистрация: 23.12.2012
Сообщений: 112
24.12.2012, 01:37  [ТС]
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
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
24.12.2012, 01:37
Помогаю со студенческими работами здесь

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

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

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

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

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


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

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

Новые блоги и статьи
Как украсить новогоднюю елку с Q# и Qiskit
EggHead 24.06.2025
Что может быть необычнее, чем применить законы квантовой механики для украшения новогодней елки? Пока другие развешивают обычные гирлянды, я решил объединить свою страсть к квантовым вычислениям с. . .
Системы нулевого доверия на C#
UnmanagedCoder 24.06.2025
Традиционная архитектура безопасности работает по принципу средневекового замка: создаём высокие стены вокруг корпоративной сети, укрепляем ворота межсетевыми экранами и системами обнаружения. . .
Снова не мой путь. Циклическое среднее, я обеими руками за проверку условия, в ракурсе данной задачи - циклическое среднее в топку.
Hrethgir 24.06.2025
Привет. Такой вопрос - нужно выводить среднее математическое между двумя направлениями, интервал значений которых может лежать в диапазоне одного оборота по кругу. Проблема заключается в том, что. . .
Деплой Flask приложения
py-thonny 23.06.2025
За годы работы с Flask я натыкался на одни и те же грабли достаточно часто, чтобы наконец научится их обходить. И сегодня хочу поделится опытом, который сбережет вам немало нервных клеток. Начнем с. . .
WebAssembly и контейнеры в .NET Aspire для оркестрации распределенных архитектур
ArchitectMsa 23.06.2025
Я наблюдаю, как WebAssembly (или просто WASM) постепенно выходит за рамки своего первоначального предназначения — исполнения кода на стороне браузера. Теперь эта технология проникает в серверную. . .
Непрерывная интеграция для пакета Python
Mr. Docker 22.06.2025
Было 4 часа утра пятницы, когда я выпустил новую версию нашей внутренней библиотеки для обработки данных. Релиз 0. 5. 2 содержал небольшой фикс для обработки дат в ISO формате, что может пойти не так?. . .
Продвинутый ETL на C# из OLTP БД в хранилище
stackOverflow 22.06.2025
Работая в сфере корпоративной аналитики, я постоянно сталкиваюсь с одним и тем же - нужны чистые, структурированные и, главное, свежие данные. Без них современные аналитические системы, машинное. . .
Мастер-класс по микросервисам на Node.js
Reangularity 21.06.2025
Node. js стал одной из самых популярных платформ для микросервисной архитектуры не случайно. Его неблокирующая однопоточная модель и событийно-ориентированный подход делают его идеальным для. . .
Управление Arduino из WPF приложения
Wired 21.06.2025
Зачем вообще связывать Arduino с WPF-приложением? Казалось бы, у Arduino есть собственная среда разработки, своя экосистема, свои способы управления. Однако при создании серьезных проектов. . .
Звёздная пыль
kumehtar 20.06.2025
Я просто это себе представляю: как создавался этот мир. Как энергия слипалась в маленькие частички. Как они собирались в первые звёзды, как во вселенной впервые появился Свет. Как эти звёзды. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2025, CyberForum.ru