Вопросы по теме 'ode'

Как разгадать жесткую оду с Python?
Я новичок в Python. Я пытаюсь переключить некоторые программы, которые у меня есть в Matlab. Мне нужно решить жесткое уравнение оды, все входные данные которого - матрицы. В Matlab я использую [ttT,uT] = ode23s('SST',t,fT);
3514 просмотров
schedule 18.12.2023

Установка условия на выходе MATLAB ode45 (во время работы)
Я решаю систему ОДУ с ode45 для многих различных параметров (эти параметры являются коэффициентами дифференциальных уравнений), и я хочу найти параметры, для которых решение меньше (не больше), чем заданное значение. Как я могу установить условие...
1023 просмотров
schedule 19.03.2024

Matlab: решение связанных ОДУ с изменяющимися во времени параметрами
Я пытаюсь решить набор ODE в Matlab (я думаю, ODE45 - правильный инструмент), которые имеют изменяющиеся во времени параметры. Вот иллюстрация: x' = x/y+z*a(t) y' = y*x+a(t)/y+b(t) z' = x*y + c(t)/z + d У меня есть данные для a(t), b(t) и...
675 просмотров
schedule 18.10.2022

Использование dopri5 для построения системы ОДУ в матричной форме
Система уравнений, которые меня интересуют, выглядит следующим образом: Я смог построить их, изменив пример, который кто-то опубликовал, выполнив следующие действия: import scipy as sp import pylab as plt import numpy as np import...
72 просмотров
schedule 21.11.2023

Выполнение одного шага ODE
Можно ли сделать один шаг с интегратором scipy ODE? Мне известен аргумент step для scipy.integrate.ode.integrate , но я не могу понять, как установить аргумент t таким образом, чтобы он гарантированно выполнял ровно один полный шаг. Очевидный...
542 просмотров
schedule 02.04.2024

Использование цикла в ODE для графического сравнения различных параметров R
Я использую пакет deSolve для построения пары дифференциальных уравнений (если интересно, прочтите http://www.maa.org/press/periodicals/loci/joma/the-sir-model-for-spread-модели-болезни-дифференциального-уравнения ). Моя конечная цель - создать...
1357 просмотров
r ode
schedule 02.02.2024

R: Решение ODE с пакетом deSolve с использованием матриц в качестве входных данных
Я пытаюсь решить обыкновенные дифференциальные уравнения в R. У меня есть матрицы в качестве начальных значений и параметров, которые я уже определил. Но когда я пытаюсь ее решить, я получаю следующую ошибку, которая не появляется, когда я ввожу...
1463 просмотров
schedule 23.10.2023

Как решить систему ODE с scipy.integrate.odeint с переменной, зависящей от времени
Чтобы узнать о решениях систем ODE на Python. В рамках этой модели существует уравнение, которое определяет численность населения на каждый день на основе коэффициента рождаемости, смертности и первоначальной численности населения. затем на...
3486 просмотров
schedule 24.04.2024

одовый решатель Matlab
Я пытаюсь понять синтаксис ode45 и не понимаю, почему иногда это ode45 (odefun, tspan, y0), а иногда ode45 (@odefun, tspan, y0). Я хотел бы понять значение символа @ перед функцией odefun. В частности, мне непонятен следующий синтаксис: (он должен...
78 просмотров
schedule 13.11.2022

Перенос дифференциального уравнения в Matlab ODE45 - была ли моя ошибка?
В течение нескольких часов я размышляю над проблемой, которую, как мне кажется, легко решить. Я просто не понимаю. В задании для Coursera (Aerial Robotics) я попросил запрограммировать контроллер PD с помощью функции ODE45 из Matlab. Целью является...
98 просмотров

Python ODE45 IndexError: назначение списка вне допустимого диапазона
Я пытаюсь продублировать сценарий ODE, который я запускаю в Matlab, на Python. Вот скрипт Матлаба: t0 = 0; tfinal = 25; q1 = 1; q2 = 1; q1dot = 0; q2dot = 0; % ODE variables times = [t0 tfinal]; stateVars=[q1 q1dot q2 q2dot]; % ODE45 options =...
302 просмотров
schedule 11.04.2024

установка переменной состояния на заданное значение с помощью ode15s
Я использую ode15 для имитации/решения набора ОДУ. Я хотел бы реализовать функцию, при которой при достижении заданного условия во время моделирования число в модели изменяется программно (например, константа индикатора) на фиксированное время, а...
76 просмотров
schedule 01.06.2024

Построение одного решения на фазовой плоскости
Я пытаюсь построить решение своей оды с помощью solver integration.odeint, однако я не получаю решения, когда пытаюсь построить его. Я не вижу, где неправильная формулировка моего кода. Пожалуйста, найдите ниже: import numpy as np import...
1220 просмотров
schedule 14.11.2023

Решите уравнение в питоне за заданный интервал времени с заданным начальным условием
Я хочу решить уравнение в python за интервал времени I = [0,10] с начальным условием (x_0, y_0) = (1,0) и значениями параметров μ ∈ {-2, -1, 0, 1, 2} с помощью функции scipy.integrate.odeint Затем я хочу построить решения (x (t; x_0, y_0),...
544 просмотров
schedule 09.11.2023

Как ввести временные шаги в scipy.integrate.RK45
Реализация scipy.integrate.RK45 не указывает, где указывать время, в которое должна выполняться интеграция. Параметр ввода "t_bound" кажется последней эпохой. Для интегрирования необходимо использовать опцию «шаг», которая не имеет ввода,...
2269 просмотров

Как пройти через точные точки с помощьюsolve_ivp?
У меня есть система ОДУ, и я хочу изменить значение переменных, когда решатель достигает точного момента времени. То, что я пытаюсь сделать, похоже на этот пример в Джулии То, что я пытался сделать, это использовать else и проверить, достигнуто...
574 просмотров
schedule 11.12.2022

Размеры входной матрицы для работы в Octave/Matlab
Я пытаюсь использовать свою собственную функцию Эйлера для решения небольших систем ОДУ с помощью обычного метода Эйлера в Octave. Я заставил его работать для матричного ввода 1x1, но я хочу использовать его и для матричного ввода 2x2. Мой код...
201 просмотров
schedule 08.12.2023

Решение дифференциальных уравнений в Matlab — In-Vitro Dissolution
Я пытаюсь решить проблему, аналогичную этой: Решение дифференциальных уравнений в Matlab Однако на этот раз речь идет не о введении препарата в подкожную клетчатку и его последующем растворении, а о более простой ситуации, когда суспензии дают...
33 просмотров
schedule 27.03.2024

Как правильно реализовать решатель системы принудительной массовой пружины с переменной, зависящей от силы, с использованием scipy.integrate.odeint
Я пытаюсь для каждого из 5 случаев численно интегрировать через функцию odeint, функцию массы пружины, причем параметр F изменяется со временем. Однако он представляет ошибку: строка 244, в odeint int(bool(tfirst))) ValueError: установка...
172 просмотров
schedule 02.12.2022

Настройка вектора u при использовании lsim
У меня есть неприятная ошибка при использовании Matlab, где я пытаюсь смоделировать систему непрерывного времени в дискретном времени. Ts = 0.01; A=[-0.313 0 56.7; 0 56.7 0; -0.0139 0 0.426]; B = [0.232; 0; 0.0203]; C = [0 1 0]; D = 0;...
39 просмотров
schedule 10.10.2022