Описание тега runge-kutta
Методы Рунге – Кутты представляют собой важное семейство неявных и явных итерационных методов, которые используются во временной дискретизации для аппроксимации решений обыкновенных дифференциальных уравнений.
1
ответ
Численный метод Рунге-Кутты Плохое приближение
Я пытаюсь использовать метод Рунге-Кутты, чтобы сравнить его с функцией lsode. Но он работает довольно плохо, все остальные методы, которые я использовал (Фовардс и Эйлер, Хеун) для сравнения с lsode, делают работу лучше, до такой степени, что они п…
12 ноя '18 в 00:18
0
ответов
Ошибка конвергенции для Рунге-Кутты в Python
У меня есть этот кусок кода (я положил только часть этого) в Python: from __future__ import division import numpy as np from pylab import * from numpy import linalg as LA a = [[5/36,2/9-sqrt(15)/15,5/36-sqrt(15)/30],[5/36+sqrt(15)/24,2/9,5/36-sqrt(1…
06 дек '17 в 20:31
1
ответ
Оценка Runge Kutta не работает для размера шага, кроме 1
Я преобразовывал реализацию метода Эйлера из 4 связанных дифференциальных уравнений в реализацию Рунге Кутты 4-го порядка. Я вполне уверен, что у меня есть общий подход, и я понял, как применять RK4, но я не занимался наполовину серьезной математико…
17 мар '16 в 01:51
1
ответ
Ошибка в больших итерациях n-мерного решателя Рунге-Кутты 4-го порядка
Код, представленный здесь, работает нормально, но, скажем, чтобы найти точки бифуркации при изменении омега для (-30,30), а не (10,30), таким образом, изменяя "int o" с 2000 на 6000, на экране появляется следующее сообщение: Необработанное исключени…
15 июн '17 в 15:09
1
ответ
Адаптивный пошаговый метод для метода Рунге-Кутты 5-го порядка в Фортране
Я хочу решить систему уравнений, используя метод Рунге-Кутты 5-го порядка с адаптивным методом шага. Я нашел полезный код, написанный Танер Акгун. Вот код: c c Adaptive Size Method for 5th Order Runge-Kutta Method c (Based on Numerical Recipes.) c c…
07 июн '16 в 12:06
1
ответ
C-код для простого маятника не дает ожидаемых результатов
Я пытаюсь написать код для сети нейронов FitzHugh NAgumo в C, используя интегратор RK4. Поскольку это не сработало, я решил попробовать что-то более простое, простую маятниковую систему. Я использую отдельные функции - одну для возврата дифференциал…
22 фев '17 в 08:35
1
ответ
Matlab, поле склона, система ODE Эйлера и Рунге-Кутта 2-го порядка
У меня есть система дифференциальных уравнений х ' = топор - на y' = bx + ay Мне нужно найти приближенное решение, используя явный метод Эйлера и Рунге-Кутта второго порядка, условия a = 0, b = 1, x (0) = 0, y (0) = 1, и, более того, Используя Wolfr…
19 мар '18 в 22:19
0
ответов
Модель Pymc3 с ODE Solver с использованием Theano
Я использую модель, где средний ответ зависит от решения ODE. Я пытаюсь приспособить эту модель, используя pymc3, но у меня возникают проблемы (связанные с отсутствием тестовых значений) при присоединении решателя ODE к модели. модель y_t Логнормаль…
14 фев '18 в 11:08
3
ответа
Использование Рунге-Кутта для решения связанных дифференциальных уравнений
У меня есть система связанных уравнений: уравнение гидростатического равновесия, уравнение неразрывности массы и уравнение состояния идеального газа. Это, в математической грамматике, \frac{dP}{dr}=- \rho*g, где \rho это плотность и g это гравитацио…
28 июн '15 в 14:55
1
ответ
Интеграция Runge kutta для отклика подвеса
Я должен использовать интеграцию Runge Kutta, чтобы решить и построить реакцию ускорения подвески транспортного средства для определенного интервала времени. Когда на графике изображен отклик, он должен иметь несколько колебаний, но из-за демпфирова…
01 ноя '16 в 07:53
1
ответ
Индекс превышает погрешность размеров матрицы в методе Рунге-Кутты: Matlab
Я пытаюсь сделать код с пошаговым шагом времени, используя метод Рунге-Кутты 4-го порядка, но сталкиваюсь с проблемами при правильном индексировании одного из моих значений. Мой код: clc; clear all; L = 32; M = 32; N = 32; % No. of elements Lx = 2; …
01 май '15 в 14:01
1
ответ
Почему численный результат отличается (RK45)?
Это тест для дифференциального уравнения с использованием runge-kutta45: f(x,y)= (-5*x - y/5)^1/8 + 10 почему численный результат отличается? Я использовал: function Rk_JL() f(x,y)= (-5*x - y/5)^1/8 + 10 tspan = 0:0.001:n y0 = [0.0, 1.0] return ODE.…
11 июн '16 в 02:22
1
ответ
Может ли Эйлер быть лучше, чем Рунге-Кутта для некоторых функций?
Я пытаюсь решить упражнения из "Нелинейной динамики и хаоса" Стивена Строгаца. Ожидается, что в упражнении 2.8.3, 2.8.4 и 2.8.5 будет реализован метод Эйлера, улучшенный метод Эйлера и метод Рунге-Кутта (4-й порядок) соответственно для задачи началь…
31 окт '18 в 08:15
1
ответ
Рунге Кутта в Фортране
Я пытаюсь реализовать метод Рунге Кутты в Фортране и сталкиваюсь с проблемой сходимости. Я не знаю, сколько кода мне нужно показать, поэтому я опишу проблему подробно и, пожалуйста, подскажите, что мне следует добавить / удалить в / из поста, чтобы …
21 май '17 в 15:18
1
ответ
Правило средней точки для Matlab
Здравствуйте, меня попросили создать код Matlab для правила средней точки. У меня есть код для метода Эйлера, поэтому я должен внести некоторые изменения, но я изо всех сил пытаюсь это сделать, у меня есть следующее function H = heun(f,a,b,ya,M) h =…
04 окт '18 в 08:43
1
ответ
Решение двойного маятника по Рунге-Кутте Фортран четвертого порядка
Я пытаюсь численно решить двойной маятник Рунге-Кутты четвертого порядка. Решаемые дифференциальные уравнения находятся на следующей странице: http://www.myphysicslab.com/dbl_pendulum.html Есть даже анимация, показывающая углы 1 и 2. Я хочу решить у…
01 дек '14 в 18:38
1
ответ
Значения для моделирования хаотического рассеяния не совпадают с базовым случаем
Мой первый пост о переполнении стека, будьте осторожны. Я написал код для отслеживания положения на плоскости x,y частицы массы M на потенциале V(r), описываемой четырехмерной системой уравнений движения M(dv/dt)=-grad V(r), dr/dt=v, Которые решаютс…
21 ноя '16 в 03:45
2
ответа
Как я могу ввести функцию и использовать ее в функции дескриптора
В этом коде пользовательские коэффициенты ввода ODE, и когда все коэффициенты в числовом коде работают нормально, но я хочу ввести f как функцию типа sin(2*t), как я могу ввести f, чтобы использовать ее в функции z2perim? я пытаюсь ввести f в качест…
30 дек '18 в 21:32
0
ответов
О спектре мощности и верхнем пределе достоверности данных временного ряда
Пока у меня есть связанная система с 5 переменными, и я использую метод Рунге-Кутты для интеграции. sigma=9.95;k=28;b=8/3;C1=0.1;C2=1;Od=1;Om=10;Sm=10; Ss=1;Spd=10;Sigma=100;C3=0.01;C4=0.01;C5=1;C6=0.001; ntime=4000; %X=zeros(5,ntime); dt=0.001; T=d…
02 янв '19 в 08:33
0
ответов
Интегрирование дифференциальных уравнений движения для орбиты в Рунге Кутте в C
Я пытаюсь интегрировать дифференциальные уравнения для планетарной орбиты, используя runge-kutta четвертого порядка в C. Первое уравнение, которое я попытался интегрировать, было для позиции, dr/dt = sqrt((2/mu)*(E-(k/r)) - pow(1, 2)/(pow(mu, 2)*pow…
13 июл '16 в 17:17