Описание тега ode45

Вопросы с тегом
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 ответ

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 ответ

Почему моя третья функция MATLAB выводит только нули при использовании ode45?

Мне нужно смоделировать отрицательную, положительную и простую регуляцию гена для моего урока системной биологии с использованием MATLAB. Проблема в том, что функции для отрицательного и простого регулирования работают, но функция положительного рег…
05 дек '18 в 19:37
0 ответов

Решаем систему ОД MATLAB с од45

У меня есть 4 переменные y1, y2,y3,y4, и я знаю их вторичные производные как функции от y1, y2, y3, y4: y1"= y1(y1,y2,y3,y4) y2"= y1(y1,y2,y3,y4) y3"= y1(y1,y2,y3,y4) y4"= y1(y1,y2,y3,y4) Я пытаюсь использовать ode45, чтобы получить y1, y2,y3,y4 Я п…
22 окт '17 в 04:15
0 ответов

Вмешательство в оду решатель, использующий для цикла

Я пытаюсь запустить систему уравнений для моделирования популяций пчел и клещей с инфекциями. Тем не менее, как только популяция пчел достигнет определенного порога размера, я хочу, чтобы половина пчел ушла, а также определенный процент клещей тоже.…
01 мар '18 в 20:59
1 ответ

Решение системы из 4 ODE второго порядка в Matlab с использованием ODE45

Мне нужно решить эту систему уравнений второго порядка, используя ODE45 в Matlab. Я знаком только с использованием ODE45, возможно, с одним или двумя уравнениями, но не так много. Вот что у меня есть, но я не уверен, как это исправить: function seco…
18 фев '16 в 05:26
0 ответов

MATLAB График зависимости переменной в оде

У меня проблема в MATLAB. Я пытаюсь построить переменную внутри функции Оде. Я называю свою функцию с этим кодом: %Above is the code for the IC and constants %Initial conditions [x0, vx0, y0, vy0] IC = [x_0; w_x; y_0; w_y]; %Simulationszeitvektor ts…
14 май '16 в 20:39
1 ответ

Как решить связанное дифференциальное уравнение в Matlab с помощью ode45

У меня есть два дифференциальных уравнения: da/dt=a(.3/a^3+.7)^1/2 и dτ/dt=1/a. Начальные условия: t=0; а = 1 и т = 0 соответственно. Как я могу решить уравнения в Matlab? Мне нужно рассчитать различные значения a, t и τ, а также построить график τ …
28 ноя '18 в 12:17
1 ответ

Почему я получаю только NaN и Inf при моделировании ODE45 в Октаве?

Я использую Octave с ODE45 для моделирования системы уравнений ODE. Но проблема в том, что моделирование ODE дает неправильные значения. Посмотрите на этот код Octave: function dx = dynamik(t, x) b1 = 1000; b2 = 2000; m1 = 10; m2 = 7; M = 2000; g = …
09 июн '17 в 16:22
0 ответов

Если операторы в функции ode45

Я использую следующий сигнал в качестве входных данных в коде Matlab, который я сделал. Я поместил его в файл функции ode45, но результаты, которые я получаю, не имеют смысла. Я попытался использовать решатель с фиксированным шагом ode4, и результат…
26 апр '16 в 17:34
1 ответ

Matlab ode45 против ode23, различные решения

Я использовал ode45 а также ode23 для системы не жестких дифференциальных уравнений. Однако оба метода сходятся к несколько другому решению. Как я могу узнать, какой из них правильный? Смотрите прикрепленный участок, где синий ode45красный ode23, Пу…
27 янв '16 в 19:12
1 ответ

Система дифференциальных уравнений с ode45 в matlab

У меня есть эта модель для глюкозы и инсулина, а также система дифференциальных уравнений: Куда: G(t) - концентрация глюкозы в плазме в момент времени t I(t) - концентрация инсулина в плазме в момент времени t X(t)- интерстициальный инсулин в момент…
26 май '16 в 07:46
0 ответов

Crank Nicolson MATLAB

Есть ли в MATLAB встроенный решатель, который решает данный ODE (порядок>1) с использованием метода Кранка-Николсона? (Обратите внимание, что у меня есть система ODE, поэтому я мог бы решить каждый из них, используя Crank-Nicolson. Я пробовал ode45,…
28 авг '15 в 16:30
0 ответов

ODE 3-го порядка в Simulink - Симуляция остановлена

Как мне запустить следующее? Я построил его в третий раз, и у меня до сих пор нет решения. Я пытаюсь запустить симуляцию в MATLAB/Simulink. Все это было сделано в этой статье, но без Simulink. Код MATLAB содержит параметры и начальные условия для ин…
21 мар '17 в 20:35
0 ответов

Matlab: ODE45+ метод съемки не работает

Я пытаюсь численно решить три одновременных функции ODE: dw/ds = sin(theta) dx1/ds = cos(theta) dtheta/ds = - a * cos(theta) * (w-w0) + a * sin(theta) * x1 Я решаю для столбца, который разделен на две половины, я хочу построить w против x1 и иметь п…
15 июл '16 в 02:15
3 ответа

Имитация функции ode45 из MATLAB в Python

Мне интересно, как экспортировать функцию MATLAB ode45 в Python. Согласно документации это должно быть следующим: MATLAB: [t,y]=ode45(@vdp1,[0 20],[2 0]); Python: import numpy as np def vdp1(t,y): dydt= np.array([y[1], (1-y[0]**2)*y[1]-y[0]]) return…
24 янв '18 в 17:14
2 ответа

ODE45 с очень большими числами в качестве ограничений

2-й ODE для решения в MATLAB: ( (a + f(t))·d²x/dt² + (b/2 + k(t))·dx/dt ) · dx/dt - g(t) = 0 Граничное условие: dx/dt(0) = v0 где t это время, x это позиция dx/dt это скорость d2x/dt2 это ускорение a, b, v0 являются постоянными f(t), k(t) а также h(…
03 янв '17 в 12:22
1 ответ

MATLAB - Как прочитать одно значение из вектора? (задействовано ode45)

Я пытаюсь работать с кодом, который включает в себя ode45, в том, что в уравнении, которое будет решено, я должен ввести одно значение из вектора, который будет меняться в зависимости от времени решения уравнения, я имею в виду, если это x=x(0), тог…
10 сен '15 в 03:09
0 ответов

ode45 с двумя степенями свободы MATLAB

Я пытаюсь решить дифференциальное уравнение с помощью функции ode45, которая включает две степени свободы, направления x и y. Когда в моей функции oem ниже я пытаюсь сохранить все четыре переменные, которые мне нужны, в матрицу для подключения к ode…
26 янв '16 в 19:41
0 ответов

Построение трех функций на одном выходе

Я определил функцию function dzdt=odefun2222(t,z) beta=1+exp(-t); delta=1+exp(-t); rho=1+exp(-t); gama=exp(-t); theta=exp(-t); f1=1/2+1/(t+1); f2=1/2+1/(t+1); dzdt=zeros(6,1); dzdt(1)=z(2)-f1*z(1); dzdt(2)=(1/4+1/(t+1)-beta)*z(1)+gama*z(3)-f1*z(2)-f…
26 июл '18 в 15:52