Описание тега ode
An ODE (ordinary differential equation, in contrast to partial differential equation) is a differential equation containing a function or functions of exactly one independent variable.
1
ответ
Ошибка при передаче аргументов в дескрипторы функций
Предположим, у нас есть следующая функция: function f=lorenz(t,x,a,b,c) % solve differential equation like this %dx/dt=a*(y-x) %dy/dt=-x*z+b*x-y %dz/dt=xy-c*z/3 f=zeros(3,1);% preallocate result f(1)=a*(x(2)-x(1)); f(2)=-x(1)*x(3)+b*x(1)-x(2); f(3)=…
19 фев '17 в 21:08
1
ответ
Чтение файла, соответствующего временному шагу в решении ODE
Я пытаюсь решить ODE в Matlab, и у меня есть следующая проблема: мой код как показано ниже: xinit=[0.19;25;0;7]; t=0:768:76800; %% 101 cells [t1,y]=ode45(@Model_Bio,t,xinit); В функции @Model_Bio у меня есть параметр, который мне нужен для чтения ег…
03 июл '15 в 10:53
1
ответ
Интеграция ODE с дискретизированными значениями
Я хочу использовать scipy.integrate.ode решатель. Я могу определить вызываемую функцию f только как массив дискретных точек (потому что это зависит от результатов интегрирования из предыдущих итераций). Но из документации кажется, что интегратор ожи…
01 янв '13 в 21:56
2
ответа
Ошибка при использовании *- внутренняя размерная матрица должна совпадать
У меня возникла проблема с задачей, где я должен сравнить методы решения ОДУ: ode45, методы Эйлера и Гаусса-Лежандра. Здесь я должен рассчитать ошибки для разных шагов. h=[0.01 0.05 0.1 0.5]; func = @(t, y) -2*y+t*sin(t); opts = odeset('Reltol',1e-1…
01 июн '18 в 21:12
0
ответов
C++ integrate_adaptive().
Прежде всего мой код: #include <iostream> #include <fstream> #include <boost/array.hpp> #include <boost/numeric/odeint.hpp> using namespace std; using namespace boost::numeric::odeint; typedef boost::array< double , 5 >…
21 апр '18 в 15:54
2
ответа
Не удалось решить простой ODE с Octave
Я новичок в Octave, поэтому я пытаюсь заставить работать несколько простых примеров, прежде чем переходить к более сложным проектам. Я пытаюсь решить ODE dy/dx = a*x+b, но безуспешно. Вот код: %Funzione retta y = a*x + b. Ingressi: vettore valori t;…
18 ноя '17 в 07:03
1
ответ
Решение динамической системы ODE с помощью Python
Я пытаюсь решить динамическую систему с тремя переменными состояния V1,V2,I3 и затем построить их на 3D-графике. Мой код пока выглядит следующим образом: from scipy.integrate import ode import numpy as np import matplotlib.pyplot as plt from mpl_too…
18 ноя '18 в 14:19
1
ответ
Вектор<double> быстрее, чем double*: почему?
Вот цикл, который я пробовал с std::vector<double> и с простой старой double*, Для 10 миллионов элементов векторная версия постоянно работает примерно в 80% случаев, когда double* версия берет; для практически любой стоимости N, вектор заметно…
14 июл '11 в 19:26
2
ответа
Несколько экземпляров scipy.integrate.ode
Я хотел бы использовать экземпляры scipy.integrate.ode (или scipy.integrate.odeint) в нескольких потоках (по одному на каждое ядро ЦП) для решения нескольких IVP одновременно. Однако в документации сказано: "Этот интегратор не является повторно вх…
15 дек '15 в 14:23
1
ответ
Python ODE45 IndexError: назначение списка вне диапазона
Я пытаюсь скопировать скрипт ODE, который я запускаю в Matlab, на Python. Вот скрипт Matlab: t0 = 0; tfinal = 25; q1 = 1; q2 = 1; q1dot = 0; q2dot = 0; % ODE variables times = [t0 tfinal]; stateVars=[q1 q1dot q2 q2dot]; % ODE45 options = odeset('rel…
29 авг '18 в 00:19
1
ответ
Должен ли я перейти на MTL4 с Eigen, если мне также нужно использовать ODE?
Я ищу матричные библиотеки C++ для работы в системе Linux (Fedora). Намерение состоит в том, чтобы реализовать непрерывные аттракторные нейронные сети и тому подобное для вычислительной нейробиологии. Я уже начал использовать Eigen. Тем не менее, я …
08 янв '13 в 03:17
1
ответ
Python ода первого порядка, как решить эту проблему с помощью Sympy
Когда я пытаюсь решить эту первую оду, используя Sympy, как показано ниже: import sympy y = sympy.Function('y') t = sympy.Symbol('t') ode = sympy.Eq(y(t).diff(t),(1/y(t))*sympy.sin(t)) sol = sympy.dsolve(ode,y(t)) csol=sol.subs([(t,0),(y(0),-4)]) # …
04 фев '17 в 08:43
1
ответ
Остановка ODE решателя с условием
Я хотел бы прекратить ODE решатель, когда зависимая переменная превысила определенное значение. Рассмотрим простую модель ODE library("deSolve") dModel<- function(t, y, parms) {list(c(y))} out<-ode(c(1),seq(0,100,1),dModel, parms=NULL) Я хочу …
09 июл '18 в 22:30
1
ответ
Численный метод Рунге-Кутты Плохое приближение
Я пытаюсь использовать метод Рунге-Кутты, чтобы сравнить его с функцией lsode. Но он работает довольно плохо, все остальные методы, которые я использовал (Фовардс и Эйлер, Хеун) для сравнения с lsode, делают работу лучше, до такой степени, что они п…
12 ноя '18 в 00:18
1
ответ
Решение DDE в матричной форме с использованием Matlab
Мне нужно решить следующее уравнение с задержкой: Xdot = A*x(t) + B*U(t) + E(t) U(t) = K*S*x(t-0.02) or U(t) = K*S*x(t-0.02)/alpha в котором: A: 6*6 matrix K and alpha are scalar S: 1*6 vector E and U and x and B: 6*1 vectors. Я использую следующий …
30 июн '17 в 07:25
1
ответ
Scipy.ode "Ошибка проверки неоднократно"
Я пытаюсь решить систему связанных сложных ODE в Python, используя scipy.ode с интегратором zvode. Но это сообщение об ошибке появляется, когда я запускаю код. ZVODE-- At T(=R1) and step size H(=R2), the error test failed repeatedly or with abs(H) =…
19 июн '18 в 15:30
1
ответ
Вектор возврата из ODE45 имеет меньшую длину, чем вектор начальных условий, при попытке решить связанный diffEQ
Я пытаюсь использовать ODE45, чтобы найти решение для 2 вращающихся стержней, вращающихся в вертикальной плоскости, которые имеют пружину кручения, которая создает момент на стержнях только тогда, когда угол между ними отличается от 90 градусов. Я п…
07 дек '16 в 15:36
1
ответ
Оценка Runge Kutta не работает для размера шага, кроме 1
Я преобразовывал реализацию метода Эйлера из 4 связанных дифференциальных уравнений в реализацию Рунге Кутты 4-го порядка. Я вполне уверен, что у меня есть общий подход, и я понял, как применять RK4, но я не занимался наполовину серьезной математико…
17 мар '16 в 01:51
1
ответ
Как исправить ошибку TypeError: в setindex! в DifferentialEquations.jl
Недавно я начал с Джулии (v1.0.3) DifferentialEquations.jl пакет. Я попытался решить простую систему ODE, с той же структурой, что и моя реальная модель, но гораздо меньше. В зависимости от решателя, который я использую, пример либо решает, либо выд…
29 янв '19 в 13:21
1
ответ
Что было бы эквивалентно "MaxSteps" с использованием решателя ODE GSL?
Я хочу воспроизвести решатель ODE, созданный с использованием Mathematica с GSL. Вот код Mathematica, который использует NDSolve: result[r_] := NDSolve[{ s'[t] == theta - (mu*s[t]) - ((betaA1*IA1[t] + betaA2*IA2[t] + betaB1*IB1[t] + betaB2*IB2[t]) +…
06 окт '13 в 14:04