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; ...

Solucionadores de ode do Matlab: mudança de estado e tempo especificado

Estou resolvendo um conjunto de ODEs (dy / dt) em t = 0, todas as condições iniciais t = 0 y_0 = (0,0,0). Posso adicionar algum número aos valores de y em mo...

Como resolver um sistema de equações diferenciais ordinárias com parâmetros dependentes do tempo

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 = ...

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 ...

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 ...

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 ...

resolvendo duas equações diferenciais de dimensão em python com scipy

Eu sou um novato em python. Eu tenho um sistema diferencial simples, que consiste em duas variáveis e duas equações diferenciais e condições iniciaisx0=1, y0=2: dx/dt=6*y dy/dt=(2t-3x)/4yagora estou tentando resolver essas duas ...

istemas ODE complexos em scipy

Estou tendo problemas para calcular a equação óptica de bloch, que é um sistema ODE de primeira ordem com valores complexos. Eu descobri que o scipy pode resolver esse sistema, mas a página deles oferece muito pouca informação e eu mal consigo ...

Várias instâncias scipy.integrate.ode

Eu gostaria de usarscipy.integrate.ode(ou scipy.integrate.odeint) emvários threads(um para cada núcleo da CPU) para resolver vários IVPs por vez. No entanto, a documentação diz: "Este integrador não é reentrante. Você não pode ter duas instâncias ...