Considere un IVP para el sistema ODE , y ( x 0 ) = y 0 . Más comúnmente, este problema se considera rígido cuando la matriz de Jacobi ∂ fy′=f(x,y)y′=f(x,y)y'=f(x,y)y(x0)=y0y(x0)=y0y(x_0)=y_0tienedosvalores propios con gran parte real negativa y valores propios con parte real negativa muy pequeña...
Tengo una EDO: u ( 0 ) = - 1tu′= - 1000 u + s i n ( t )u′=−1000u+sin(t)u'=-1000u+sin(t) u ( 0 ) = - 11000001tu(0 0)=-11000001u(0)=-\frac{1}{1000001} Sé que este ODE particular es rígido, analíticamente. También sé que si usamos un método de paso de tiempo explícito (hacia adelante) (Euler,...
Estoy resolviendo el sistema de reacción-difusión de Turing con el siguiente código C ++. Es demasiado lento: para una textura de 128x128 píxeles, el número aceptable de iteraciones es 200, lo que resulta en 2.5 segundos de retraso. Necesito 400 iteraciones para obtener una imagen interesante, pero...
En las clases de FEM, generalmente se da por sentado que la matriz de rigidez es positiva definida, pero no puedo entender por qué. ¿Alguien podría dar alguna explicación? Por ejemplo, podemos considerar el problema de Poisson: cuya matriz de rigidez es: que es simétrico y positivo definido....