Вопросы по теме 'runge-kutta'

Рунге-Кутта (RK4) для системы дифференциальных уравнений на Java
Это предложение в основном является результатом этой темы: Дифференциальные уравнения в Java . В основном, Я попытался последовать совету Джейсона С. и реализовать численные решения дифференциальных уравнений с помощью метода Рунге-Кутта (RK4)....
7226 просмотров

Нужна помощь в исправлении моей реализации RK4
Я был бы признателен, если бы кто-нибудь, более опытный в реализации, помог бы мне обнаружить мою логическую ошибку в моем текущем коде. В течение последних нескольких часов я застрял с реализацией и тестированием различных размеров шага для...
402 просмотров

Задача Джулии - решатель Рунге-Кутты модели ФитцХью – Нагумо в частных производных
Я новичок в языке программирования Julia, поэтому не знаю, как оптимизировать код. Я слышал, что Julia должна быть быстрее по сравнению с Python, но я написал простой код Julia для решения задачи ФитцХью– Модель Нагумо , и, похоже, она не быстрее...
617 просмотров
schedule 07.07.2023

Аппроксимация ошибки Рунге-Кутты 4-го порядка для различных временных шагов
from __future__ import division import numpy as np import matplotlib.pyplot as plt def f(x, t): #function return -x def exact(t): #exact solution return np.exp(-t) def Rk4(x0, t0, dt): #Runge-Kutta Fourth Order Error t...
321 просмотров
schedule 27.03.2023

Рунге – Кутта 4-го порядка
Давайте рассмотрим, что у меня есть система из 4 ОДУ: dX / dt = F (X), где X - вектор (4-размерность), а F: R ^ 4 -> R ^ 4. F называется vectorDE_total_function, и я пытаюсь вычислить решение с помощью RK-4. def solvingDES(): previous_vector...
409 просмотров

Решите систему уравнений с помощью Runge Kutta 4: Matlab
Я хочу решить систему ТРЕХ дифференциальных уравнений с помощью метода Рунге Кутта 4 в Matlab ( Ode45 не допускается). После долгих поисков все, что я смог найти в Интернете, - это либо непонятные примеры, либо общие объяснения, которые вообще...
9824 просмотров
schedule 20.10.2022

Пример кода адвекции частиц 4-го порядка Рунге-Кутты
Я пытался внедрить интеграцию RK4 в симуляцию, которую я делаю. Приведенная ниже функция - это моя лучшая попытка использовать RK4 для интегрирования по силовому полю в трех измерениях на основе уравнений на странице 12 по адресу этот сайт. В...
568 просмотров
schedule 19.01.2023

Реализация Рунге-Кутты для системы двух дифференциальных уравнений
У меня есть следующая система дифференциальных уравнений: И согласно бумаге, они сказали, что я могу решить это численно, используя РК 4-го порядка. Как видите, два последних уравнения связаны, и я построил матрицу (Вероятность), которая...
1677 просмотров

Детерминированная и стохастическая часть уравнения
Я ищу численный метод, который может решить как детерминированное, так и стохастическое уравнение. Я знаю, что в детерминированном случае метод РК четвертого порядка ценен и очень эффективен. К сожалению, они не были успешно применены к...
221 просмотров

Реализация метода Ньютона для поиска начальных значений с Дорманом Принцем для решения дифференциальных уравнений на языке C
Следующий код работает как шарм для решения в нем системы дифференциальных уравнений (функция fcn в коде) с правильными начальными значениями. Однако суть задачи состоит в том, чтобы заменить начальные значения y_1 (0) и y_2 (0) некоторыми случайными...
74 просмотров

Метод Рунге-Кутта четвертого порядка для орбит вокруг Солнца, задающий элемент массива с последовательностью
Я реализовал метод runge-kutta четвертого порядка, как всегда, но на этот раз я получил ошибку Value, я полагаю, это для решения двух дифференциальных уравнений с использованием функции, которая возвращает массив из четырех элементов. Может...
311 просмотров

Метод Рунге-Кутты для численного решения нелинейного дифференциального уравнения второго порядка - Octave
Итак, у меня есть это аналитически неразрешимое нелинейное дифференциальное уравнение второго порядка: x '' = (- 1 / x) ((x ') ² + gx-gh + ax') к которому я некоторое время пытался применить метод Рунге-Кутта в Octave, но не получил хороших...
163 просмотров

Точно ли SciPy RK45 (resolve_ivp) подсчитывает количество вычислений функции?
Я хочу измерить производительность моего собственного интегратора ODE по сравнению с SciPy RK45 . Таким образом, мне нужно точно знать количество вычислений функций правой части, которые использует RK45. Кто-нибудь знает, точное ли число сол.нфев,...
88 просмотров
schedule 19.06.2022

Валидация толерантности в Рунге Кутта Фельберг
Я написал код для решения системы ОДУ с использованием метода Рунге Кутта Фельберга: import numpy as np import matplotlib.pyplot as plt import numba import time start_time = time.clock() @numba.jit() def S1(x,u): alpha = 0.01 C1 =...
69 просмотров
schedule 23.06.2022