Pytania otagowane jako ode

Zwykłe równania różniczkowe (ODE) zawierają funkcje tylko jednej zmiennej niezależnej i jednej lub więcej ich pochodnych w odniesieniu do tej zmiennej. Ten znacznik jest przeznaczony do zadawania pytań dotyczących modelowania zjawisk za pomocą ODE, rozwiązywania problemów z różniczkami różniczkowymi i innych powiązanych aspektów.

17
Czy istnieje wysokiej jakości nieliniowy solver programowania dla Pythona?
Mam kilka trudnych, niewypukłych problemów globalnej optymalizacji do rozwiązania. Obecnie używam MATLAB's Optimization Toolbox (konkretnie fmincon()z algorytmem = 'sqp'), co jest dość skuteczne . Jednak większość mojego kodu znajduje się w języku Python i chciałbym również przeprowadzić optymalizację w języku Python. Czy istnieje solver NLP z powiązaniami Pythona, z którym …

3
Jaki jest stan techniki w równoległych metodach ODE?
Obecnie szukam równoległych metod integracji ODE. Istnieje wiele nowej i starej literatury opisującej szeroki wachlarz podejść, ale nie znalazłem żadnych ostatnich badań ani artykułów przeglądowych opisujących ogólnie ten temat. Jest książka Burrage [1], ale ma ona prawie 20 lat i dlatego nie obejmuje wielu bardziej nowoczesnych pomysłów, takich jak algorytm …

3
Jak można przetestować numeryczną implementację solvera ODE?
To pytanie zostało przeniesione z Mathematics Stack Exchange, ponieważ można na nie odpowiedzieć w Computational Science Stack Exchange. Migrował 6 lat temu . Zaraz zacznę pracować nad biblioteką oprogramowania numerycznych solverów ODE i mam problemy z opracowaniem testów dla implementacji solvera. Moją ambicją jest, że biblioteka ostatecznie obejmie rozwiązują zarówno …
26 ode 


2
Co oznacza „symplektyczny” w odniesieniu do integratorów numerycznych i czy używa ich odeint SciPy?
W tym komentarzu napisałem: ... domyślny integrator SciPy, który, jak zakładam, używa tylko metod symplektycznych. w którym mam na myśli SciPy odeint, który używa albo „metody niesztywnej (Adamsa)”, albo „metody sztywnej (BDF)”. Według źródła : def odeint(func, y0, t, args=(), Dfun=None, col_deriv=0, full_output=0, ml=None, mu=None, rtol=None, atol=None, tcrit=None, h0=0.0, hmax=0.0, …


4
Definicja sztywnego układu ODE
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).∂ fy′=f(x,y)y′=f(x,y)y'=f(x,y)y(x0)=y0y(x0)=y0y(x_0)=y_0∂f∂y(x0,y0)∂f∂y(x0,y0)\frac{\partial f}{\partial y}(x_0,y_0) Z drugiej strony, w przypadku …
17 ode  stiffness 

2
Opcje rozwiązywania systemów ODE na GPU?
Chciałbym wdrożyć systemy rozwiązywania ODE na GPU w „trywialnie równoległym” ustawieniu. Na przykład, wykonując analizę wrażliwości z 512 różnymi zestawami parametrów. Idealnie chciałbym wykonać rozwiązywanie ODE za pomocą inteligentnego adaptacyjnego solvera pomiaru czasu, takiego jak CVODE, zamiast stałego pomiaru czasu, takiego jak Forward Euler, ale działającego na GPU NVIDIA zamiast …
16 ode  gpu 

5
Dlaczego numeryczne rozwiązanie ODE odchodzi od niestabilnej równowagi?
Chciałbym zasymulować zachowanie systemu podwójnie wahadłowego. System jest manipulatorem robota o 2 stopniach swobody, który nie jest uruchamiany i dlatego będzie zachowywał się głównie jak podwójne wahadło pod wpływem grawitacji. Jedyną główną różnicą w przypadku podwójnego wahadła jest to, że składa się on z dwóch sztywnych ciał o właściwościach masy …

3
Metody numeryczne dla nieciągłych ODE rs
jakie są najnowocześniejsze metody numerycznego rozwiązywania ODE z nieciągłą prawą stroną? Najbardziej interesują mnie częściowe gładkie funkcje prawej strony, np. Znak. Próbuję rozwiązać równanie następującego typu: x˙v˙=v={(|Fexternal|−|Ffriction|)sign(Fexternal)0:|Fexternal|&lt;|Ffriction|:otherwisex˙=vv˙={(|Fexternal|−|Ffriction|)sign⁡(Fexternal):|Fexternal|&lt;|Ffriction|0:Inaczej\begin{align*} \dot x &= v\\ \dot v &= \begin{cases} (|F_\text{external}| - |F_\text{friction}|) \mathop{\rm sign} (F_\text{external}) & :|F_\text{external}| < |F_\text{friction}|\\ 0 & : \text{otherwise} \end{cases} …
15 ode 

4
Optymalna metoda ODE dla stałej liczby ocen RHS
W praktyce czas wykonywania rozwiązywania liczbowego IVP x˙(t)=f(t,x(t)) for t∈[t0,t1]x˙(t)=f(t,x(t)) for t∈[t0,t1] \dot{x}(t) = f(t, x(t)) \quad \text{ for } t \in [t_0, t_1] x(t0)=x0x(t0)=x0 x(t_0) = x_0 jest często zdominowane przez czas oceny prawej strony (RHS)fff . Załóżmy zatem, że wszystkie inne operacje są natychmiastowe (tj. Bez kosztów obliczeniowych). …

3
Najlepsze praktyki dotyczące opisywania modeli opartych na agentach
Pracuję dość intensywnie w biologii matematycznej / epidemiologii, gdzie większość prac związanych z modelowaniem / naukami obliczeniowymi jest nadal zdominowana przez zestawy ODE, co prawda czasami dość skomplikowane zestawy. Jedną z zalet tych modeli jest to, że można je łatwo opisać i powielić. Tabela wartości parametrów i same równania, a …

3
Wybór wielkości kroku za pomocą ODE w Matlabie
Cześć i dzięki za poświęcenie czasu na spojrzenie na moje pytanie. To jest zaktualizowana wersja mojego pytania, które opublikowałem wcześniej w physics.stackexchange.com Obecnie badam wirnik ekscytonowy 2D Kondensat Bosego-Einsteina i jestem ciekawy stanu podstawowego tego układu. Matematyczna metoda dojścia do stanu podstawowego nazywa się metodą wyobrażoną czasu . Metoda jest …

1
Algorytmy liniowego układu ODE
Zastanawiam się: jaki jest najlepszy algorytm do rozwiązania gdzie jest prawdziwą macierzą . A nie jest wyraźnie zależne od czasu, zwykle rzadkie, ale niekoniecznie pasmowe. Jego wartości własne mają nie dodatnie części rzeczywiste. A jest również diagonalizowalny, ale może być zbyt duży, aby pełna diagonalizacja była wydajna obliczeniowo.reuret= A ureuret=ZAu\begin{equation} …


Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.