Вопросы по теме 'runge-kutta'
Рунге Кутта в C для уравнения Лоренца
я пытаюсь вычислить систему Лоренца, используя метод Рунге Кутты, но я не могу найти, где в моем коде есть ошибка. Когда я запускаю его, система переходит в статическую точку, и я должен получить бабочку (аттрактор Лоренца). Я думаю, что это что-то...
3128 просмотров
schedule
24.12.2023
Абсолютная погрешность методов ODE45 и Рунге-Кутта по сравнению с аналитическим решением
Буду признателен, если кто-то поможет со следующей проблемой. У меня есть следующий ODE:
dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1)
Я решил (1) двумя разными способами. С помощью метода Рунге-Кутты (4-й порядок) и с помощью...
6659 просмотров
schedule
20.12.2022
Рунге-Кутта (RK4) в Haskell, проблема с системой типов
Я пытаюсь реализовать Рунге-Кутта 4-го порядка в Haskell, но мне сложно использовать систему типов Haskell для этой задачи. Может кто поможет? Я хочу изменить типы State и DState на классы типов в следующем коде:
data State = State Double...
358 просмотров
schedule
19.02.2022
Ошибка MATLAB Runge Kutta с использованием inline/subsref
Кто-нибудь может помочь мне с функцией А, почему бы не запустить? Не могу понять ошибку, может плохо пишет? потому что в B все в порядке с графикой, проверьте мой метод функций Рунге Кутта 2 (f, g) и исходную задачу из mathcad извините за мой плохой...
289 просмотров
schedule
25.02.2022
Метод Рунге-Кутты в Python
Я написал код метода рунге-кутта на питоне, но каждый раз, когда программа реализует какое-либо исчисление, программе требуется дифференциальное уравнение.
это мой код:
from math import *
import numpy as np
#Initial Values
n=input("Enter the...
3057 просмотров
schedule
05.06.2024
Оптимизация производительности RK4 в Python
Я разрабатываю классическую схему Рунге-Кутты (RK4) для большого количества связанных уравнений в Python 2.7. Поскольку будет более сотни связанных уравнений 1-го порядка, циклы for будут чертовски большими, и я ищу некоторые подсказки по...
534 просмотров
schedule
27.06.2023
Использование лямбда-функций в алгоритме RK4
Есть два способа реализации классической схемы Рунге-Кутта в Python, показанные здесь . Первый использует лямбда-функции, второй - без них.
Какой из них будет быстрее и почему?
352 просмотров
schedule
17.11.2023
Mathematica NDSolve выдает ошибку
Мне нужно решить дифференциальное уравнение вида w'=g(t,w(t)) , где g определяется следующим образом
g[t_, w_] := {f1[t, {w[[3]], w[[4]]}], f2[t, {w[[3]], w[[4]]}], w[[1]],w[[2]]};
и f1, f2
f1[t_, y_] := Sum[\[Mu][[i]] (s[[i]] -...
117 просмотров
schedule
15.01.2023
Решение уравнения Лоренца с помощью РК-4 на C ++
Я написал код для решения уравнений Лоренца методом RK-4 на C ++. Мне нужно построить график аттрактора, и мне трудно решить связанное дифференциальное уравнение 3 первого порядка с использованием метода RK-4. Вот мой код:
/* Solving 3 coupled...
810 просмотров
schedule
09.07.2023
C ++: Простой калькулятор Рунге-Кутты четвертого порядка
Я делаю калькулятор для решения нескольких уравнений Рунге-Кутты четвертого порядка, которые мы делали в классе, и хотя мне удалось заставить калькулятор работать и работать, значения, которые он дает, не совсем правильные, и я не могу понять, почему....
112 просмотров
schedule
17.08.2022
Как установить фиксированный размер шага с помощью scipy.integrate?
Я ищу способ установить фиксированный размер шага для решения моей проблемы начального значения методом Рунге-Кутта в Python. Соответственно, как я могу сказать scipy.integrate.RK45 поддерживать постоянное обновление (размер шага) для своей...
3379 просмотров
schedule
14.09.2022
Как решить следующий вопрос, используя предоставленный метод Рунге-Кутта на python
Тело вопроса:
Парашютист массой м при вертикальном свободном падении испытывает силу аэродинамического сопротивления F = cy'² (' c раз y prime square '), где y измеряется вниз от начала падения, а y является функцией времени ( y'...
189 просмотров
schedule
03.08.2022
Есть ли лучший способ привести в порядок этот кусок кода? Это 4-й порядок Рунге-Кутты с 4 ODE.
def xpos(x,y,t): #dx/dt = v_x
return vx
def ypos(x,y,t): #dy/dt = v_y
return vy
def xvel(x,y,t): #dv_x/dt = -GMx/r^3
return -G*M*x/((x)**2 + (y)**2)**1.5
def yvel(x,y,t): # dv_y/dt = -GMy/r^3
return -G*M*y/((x)**2 +...
108 просмотров
schedule
04.12.2022
Трудности с функцией RK4 VBA. Функция возвращает только #ССЫЛКА! ценность
Н/Д. Уже не вопрос......
133 просмотров
schedule
10.12.2022
Рунге Кутта прекращается при достижении состояния
Пример прыгающего мяча, показанный в решателе Matlab ODE, имеет способ прекратить интеграцию, когда срабатывает четное число ( https://www.mathworks.com/help/matlab/math/ode-event-location.html ). В этом примере он завершается, когда высота равна 0....
42 просмотров
schedule
30.03.2024
Как создать единый вектор с двумя измерениями?
Я использовал уравнение движения (закон Ньютона) для простого сценария пружины и массы, включив его в данное 2-е уравнение ОДУ y + (k / m) x = 0; у (0) = 3; у '(0) = 0.
Затем я смог запустить код, который вычисляет и сравнивает точное решение с...
60 просмотров
schedule
11.07.2022
Проблема по алгоритму Рунге Кутта Фельберга
Я написал код для Рунге-Кутты 4-го порядка, который отлично работает для системы дифференциальных уравнений:
import numpy as np
import matplotlib.pyplot as plt
import numba
import time
start_time = time.clock()
@numba.jit()
def V(u,t):...
272 просмотров
schedule
07.05.2022
TypeError: объект «модуль» не может быть вызван с помощью rungekutte
у меня есть домашнее задание, и я использовал свой код здесь, чтобы решить симметричное уравнение с использованием метода Рунге-Кутты
import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as integrate
from...
67 просмотров
schedule
29.06.2023