Rozważ IVP dla systemu ODE , . Najczęściej ten problem jest uważany za sztywny, gdy macierz Jacobiego ma zarówno wartości własne o bardzo dużej ujemnej części rzeczywistej, jak i wartości własne o bardzo małej ujemnej części rzeczywistej (uważam tylko za stabilne walizka).∂...
Mam ODE: u′=−1000u+sin(t)u′=−1000u+sin(t)u'=-1000u+sin(t) u(0)=−11000001u(0)=−11000001u(0)=-\frac{1}{1000001} Wiem, że ten konkretny ODE jest sztywny, analitycznie. Wiem również, że jeśli użyjemy jawnej (krokowej) metody krokowej (Euler, Runge-Kutta, Adams itp.), Metoda powinna zwrócić bardzo...
Rozwiązuję układ dyfuzyjny reakcji Turinga za pomocą następującego kodu C ++. Jest zbyt wolny: dla tekstury 128 x 128 pikseli dopuszczalna liczba iteracji wynosi 200 - co powoduje 2,5 sekundy opóźnienia. Potrzebuję 400 iteracji, aby uzyskać interesujący obraz - ale 5 sekund oczekiwania to za dużo....
W klasach MES zwykle przyjmuje się za pewnik, że macierz sztywności jest pozytywnie określona, ale po prostu nie rozumiem, dlaczego. Czy ktoś mógłby podać jakieś wyjaśnienie? Na przykład możemy rozważyć problem Poissona: -∇2)u = f,−∇2u=f, -\nabla^2 u = f, którego macierz sztywności jest: K.I...