ode

0зной

1ответ

Я хочу передать простой аргумент в моей функции ode45. Моя функция заключается в следующем: function dxdt = state(t,x,vgth,vgval) p=1; k=10^0.7; window1=1-((2*x)-1).^(2*p); dxdt=k*(vgval-vgth+1.2)

0зной

1ответ

Ошибка в `ode45` matlab

Основная функция ode приведена ниже. function dxdt = state(t,x,vgth,vgval1,vgval2) vgval=vgval1+vgval2; p=1; k=10^0.7; window1=1-((2*x)-1).^(2*p); dxdt=k*(vgval-vgth+1.2)*window1; end Сценарий

0зной

1ответ

Расчет временных шагов в Matlab ODE solver

Я попытался выяснить, как MATLAB вычисляет размер шага (а не начальный) для решения ODE с помощью, например, решения ode45. Исходный код действительно сложный, так кто-нибудь знает, что он работает?

3зной

1ответ

нужно, чтобы лучше понять, как RTOL, Атол работа в scipy.integrate.odeint

Здесь scipy.integrate.odeint вызывается с шестью различными стандартными проблемами оды с rtol = atol от 1E-06 до 1E-13. Я посмотрел на максимальную разницу между результатами при всех больших допуска

2зной

2ответ

Решение системы одов (с изменением константы!) С использованием scipy.integrate.odeint?

В настоящее время у меня есть система одов с постоянной времени. Например. def fun(u, t, a, b, c): x = u[0] y = u[1] z = u[2] dx_dt = a * x + y * z dy_dt = b * (y-z) dz_d

0зной

1ответ

Интеграция уравнения свободы множественной степени свободы

Я хотел бы найти отклик простой динамической системы с 2 степенями свободы, определяемой ODE: m ddx(t) + k ddx(t) = f(t). Система возбуждается с силой f1 = f01 * sin(omega*t) в точке 1: Матрицы массы