Resultados da pesquisa a pedido "ode"
Substituindo valores negativos em um modelo (sistema de ODEs) por zero
Atualmente, estou trabalhando na solução de um sistema de equações diferenciais ordinárias usandodeSolve, e estava pensando se há alguma maneira de impedir que os valores diferenciais das variáveis fiquem abaixo de zero. Eu já vi algumas outras ...
atualização do ode scipy set_f_params dentro da função definida como set_solout
Ao integrar um ode ao scipy, o ode aceita uma função com mais argumentos que t e y. Por exemplo: def fun(t, y, param1, param2): e o valor desses argumentos pode ser definido usandoset_f_params método. No entanto, ao usar tambémset_solout ...
Equações diferenciais em Java
Estou tentando criar um programa de simulação simples do modelo SIR-epidemias em java. Basicamente, o SIR é definido por um sistema de três equações diferenciais: S '(t) = - l (t) * S (t) I '(t) = l (t) * S (t) - g (t) * I (t) R '(t) = g (t) * I ...
evento aleatório matlab / oitava ode45
Tenho alguns problemas para entender como implementar eventos em oitava / matlab, na resolução de equações diferenciais. Considere, por exemplo, este código simples para resolver a equação diferencial y '= -y: function dy = odefun(t,y) dy = -y; ...
Erro absoluto dos métodos ODE45 e Runge-Kutta em comparação com solução analítica
Gostaria que alguém pudesse ajudar com o seguinte problema. Eu tenho o seguinte ODE: dr/dt = 4*exp(0.8*t) - 0.5*r ,r(0)=2, t[0,1] (1)Eu resolvi (1) de duas maneiras diferentes. Por meio doMétodo Runge-Kutta (4ª ordem) e por meio deode45 ...
Scipy.integrate.ode.set_solout funciona?
oscipy.integrate.ode a interface para rotinas de integração fornece um método para interromper a integração se uma restrição for violada a qualquer etapa,set_solout. No entanto, não consigo fazer esse método funcionar, mesmo nos exemplos mais ...
Solução de pêndulo duplo usando GSL
Estou aprendendo a usar o GSL para resolver a ODE. Eu queria resolver o problema do pêndulo duplo usando as funções GSL ODE. Eu decidi usar estas ...
Erro ODE Runge Kutta MATLAB
então estou tentando criar uma função Runge Kutta e este é o meu código: function [t,U] = RK(f, n, eta, interv) h = (interv(2)-interv(1))/n; t = interv(1):h:interv(2); v(1) = eta(1); w(1) = eta(2); for i=1:n k1 = f([v(i),w(i)]); k2 = ...
Resolvendo um sistema de odes (com constante de alteração!) Usando scipy.integrate.odeint?
Atualmente, tenho um sistema de odes com uma constante dependente do tempo. Por exemplo. 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_dt = -x*y+c*y-z return [dx_dt, dy_dt, dz_dt]As constantes são ...
precisa entender melhor como o rtol, atol funciona em scipy.integrate.odeint
Aquiscipy.integrate.odeint é chamado com seis problemas diferentes de ode padrão comrtol = atol de1E-06 para1E-13. Eu observei a diferença máxima entre os resultados em todas as tolerâncias maiores, menos as menores, para obter algum tipo de ...