numerical-integration

0зной

1ответ

Решение ОДУ с ode45 после интерполирования

Я работаю со следующим кодом и не может найти то, что это не так: xx = 0:1/50:1; v = 3.*exp(-xx)-0.4*xx; xq = xx; vq = @(xq) interp1(xx,v,xq); tspan = 0:1/50:1; x0 = 3; [~, y2] = ode45(@(t,x)vq

8зной

2ответ

Абсолютная ошибка методов ODE45 и Runge-Kutta по сравнению с аналитическим решением

Я был бы признателен, если кто-то может помочь в следующем вопросе. У меня есть следующие: ODE dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1) я решил (1) двумя различными способами. С помощью мет

1зной

1ответ

с использованием опции «Refine» в ode45 Matlab

Я пытаюсь использовать ode45 в MAtlab и хочу исправить количество точек, которые использует MAtlab (количество шагов времени). Использование опции 'refine' в ode45 кажется не поможет. Например, если я

2зной

2ответ

Трапецеидальное правило в Python

Я пытаюсь реализовать трапецеидальное правило в Python 2.7.2. Я написал следующую функцию: def trapezoidal(f, a, b, n): h = float(b - a)/n s = 0.0 s += h * f(a) for i in range(1, n

2зной

1ответ

Проблемы с квад в SymPy

Может кто-нибудь объяснить, почему: from sympy.mpmath import quad x, y = symbols('x y') f, g = symbols('f g', cls=Function) f = x g = x+1 u_1 = lambda x: f + g quad(u_1,[-1,1]) дает ошибку и f

1зной

1ответ

Использование локатора событий в Matlab при решении ODE с несколькими выходными аргументами

У меня возникла проблема с локатором событий в Matlab. Я смотрю на связанный ODE, который представляет два интерфейса, которые нестабильны. У меня нет проблем с этим, все, что я пытаюсь сделать, это н

2зной

1ответ

matlab/октава случайное событие ode45

У меня есть некоторые проблемы в понимании того, как реализовать события в октаве/матлабе в решении дифференциальных уравнений. Рассмотрим, например, этот простой код для решения дифференциального ура