Pytania oznaczone «stability»

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

Twierdzenie Lax o równoważności stwierdza, że ​​spójność i stabilność schematu numerycznego dla problemu liniowej wartości początkowej jest koniecznym i wystarczającym warunkiem konwergencji. Jednak w przypadku problemów nieliniowych metody numeryczne mogą być bardzo prawdopodobne w przypadku...

Próbuję zaimplementować następującą funkcję w zmiennoprzecinkowym podwójnej precyzji z niskim błędem względnym : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Jest to szeroko stosowane w aplikacjach statystycznych w celu dodania...

Chciałbym wiedzieć, czy istnieje szybki sposób na obliczenie odległości euklidesowej dwóch wektorów w oktawie. Wydaje się, że nie ma do tego żadnej specjalnej funkcji, więc czy powinienem po prostu użyć formuły z

Metody numeryczne rozwiązywania PDE (lub ODE) dzielą się na dwie szerokie kategorie: metody jawne i niejawne. Metody niejawne pozwalają na większe stabilne kroki czasowe, ale wymagają więcej pracy na krok. W przypadku hiperbolicznych PDE powszechnie uważa się, że metody niejawne zwykle się nie...

W gazecie natknąłem się na zagadkową uwagę PJ van der Houwen, Rozwój metod Runge-Kutty dla równań różniczkowych cząstkowych, Appl. Num. Matematyka 20: 261,1996 W wierszach 8ff na stronie 264 van der Houwen pisze: „W przypadku wielomianów Taylora oznacza to, że wyobrażony przedział stabilności...

Po wykonaniu około matematyki związane ze stabilnością elementów problemu 3D Stokesa I nieznacznie szoku sobie sprawę, że nie jest stabilna przez dowolną czworościenne siatki. Mówiąc dokładniej, jeśli masz element, w którym wszystkie węzły i trzy z czterech aspektów leżą na granicy domeny z...

Pracuję nad rozwiązaniem sprzężonych jednowymiarowych równań poroelastyczności (model Biota ), podanych jako: ∂−(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 w dziedzinieΩ=(0,1)i z warunkami brzegowymi:...

Zaimplementowałem solver z Eulerem wstecznym w Pythonie 3 (używając numpy). Dla własnej wygody i jako ćwiczenie napisałem również małą funkcję, która oblicza przybliżoną różnicę skończoną różnicy gradientu, aby nie zawsze musiałem określać analitycznie jakobian (jeśli to w ogóle...

Przepływy o dużej liczbie Reynoldsa wytwarzają bardzo cienkie warstwy przyścienne. Jeśli rozdzielczość ściany jest używana w symulacji dużego wiru, współczynnik kształtu może być rzędu . Wiele metod staje się niestabilnych w tym systemie, ponieważ stała inf-sup degraduje się jako pierwiastek...

Załóżmy, że mam funkcję o wartościach rzeczywistych niektórych zmiennych które chcę ocenić liczbowo. Zasadniczo wzór na może zawierać produkty, racjonalności, funkcje transcendentalne itp. I będzie musiał długo badać analityczną stabilność numeryczną. Albo będzie to przynajmniej czasochłonne, aby...

W tej chwili mam kod, który używa algorytmu Crank-Nicholson, ale myślę, że chciałbym przejść do algorytmu wyższego rzędu w celu pomiaru czasu. Wiem, że algorytm Crank-Nicholson jest stabilny w dziedzinie, w której chcę pracować, ale martwię się, że niektóre inne algorytmy mogą nie być. Wiem, jak...

Przeczytałem to (1) Czynności źle uwarunkowane należy wykonać przed czynnościami dobrze uwarunkowanymi. Jako przykład należy obliczyć x z- yzxz-yzxz-yz jako ( x - y) z(x-y)z(x-y)z ponieważ odejmowanie jest źle uwarunkowane, podczas gdy mnożenie nie. Jednak analiza błędów pierwszego rzędu obu...

Korzystałem z SVD Intela MKL ( dgesvdprzez SciPy) i zauważyłem, że wyniki są znacząco różne, kiedy zmieniam precyzję pomiędzy float32i float64kiedy moja matryca jest źle uwarunkowana / nie w pełni ustawiona. Czy istnieje przewodnik dotyczący minimalnej ilości regularyzacji, którą powinienem dodać,...

Czytam artykuł [1], w którym rozwiązują następujące równanie nieliniowe przy użyciu metod różnic skończonych. Analizują także stabilność schematów za pomocą analizy stabilności von Neumanna. Jednak, jak zdają sobie sprawę autorzy, dotyczy to tylko liniowych PDE. Autorzy więc to poprzez...