Pytania oznaczone «numerical-analysis»

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

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

Jaki jest stan techniki w przybliżeniu całek silnie oscylacyjnych zarówno w jednym wymiarze, jak i w wyższych wymiarach z dowolną

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

Nauczyłem się o Metodzie Elementów Skończonych (także trochę o innych metodach numerycznych), ale nie wiem, jaka jest dokładnie definicja tych dwóch błędów i różnic między

Zgodnie z moim poprzednim pytaniem próbuję zastosować warunki brzegowe do tej niejednorodnej siatki o skończonej objętości, Chciałbym zastosować warunek brzegowy typu Robin do lhs domeny ( , tak aby:x=xL)x=xL)x=x_L) σL=(dux+au)∣∣∣x=xLσL=(dux+au)|x=xL \sigma_L = \left( d u_x + a u \right)...

Cichą wiedzę można uzyskać z doświadczenia, zastanawiałem się tylko, czy ktoś wcześniej widział coś podobnego do tego. Wykres pokazuje początkowy warunek (zielony) dla równania rada-dyfuzja, następnie roztwór przy iteracji 200 (niebieski), a następnie ponownie przy iteracji 400...

Próbuję zrozumieć niektóre wyniki i doceniłbym kilka ogólnych uwag na temat rozwiązywania problemów nieliniowych. Równanie Fishera (nieliniowe PDE z dyfuzją reakcyjną), ut= dux x+ βu ( 1 - u ) = F.( u )ut=reuxx+βu(1-u)=fa(u) u_t = du_{xx} + \beta u (1 - u) = F(u) w formie dyskretnej, u′jot=...

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

Regularnie biorę udział w tak zwanych „konkursach programistycznych”, w których rozwiązujesz trudne problemy algorytmiczne za pomocą własnego kodu i umiejętności rozwiązywania problemów w ograniczonym czasie. Aby zapoznać się z przykładowymi przykładami tego, jak mogą one wyglądać, wyszukaj...

Częstym problemem w statystyce jest obliczanie pierwiastka kwadratowego odwrotnego symetrycznej dodatniej macierzy określonej. Jaki byłby najbardziej efektywny sposób obliczenia tego? Natknąłem pewnym literaturze (które nie zostały jeszcze przeczytane) oraz jakiegoś przypadkowego kodu R tutaj ,...

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

Mam problem z implementacją funkcji numerycznie. Cierpi na tym fakt, że przy dużych wartościach wejściowych wynik jest bardzo dużą liczbą razy bardzo małą liczbą. Nie jestem pewien, czy katastrofalne anulowanie jest właściwym terminem, więc proszę mnie poprawić, jeśli tak jest. Dowody na coś nie...

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

Niech a być funkcją w tych zmiennych.f( → x ):[0,1]n→Cx⃗ = ( x1, x2), … , Xn) ∈ [ 0 , 1 ]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Czy istnieje schemat rekurencyjny dla tej iterowanej...

Podczas dyskretyzacji MES i rozwiązania problemu dyfuzji reakcji, np. przy 0 < ε ≪ 1 (pojedyncze zaburzenie), rozwiązanie problemu dyskretnego zwykle wykazuje warstwy oscylacyjne blisko granicy. Z Ohm = ( 0 , 1 ) , ε = 10 - 5 i liniowych elementów skończonych, roztwór u h wygląda...

Próbuję odwrócić liczbowo następującą transformację całkowitą: F(y)=∫∞0yexp[−12(y2+x2)]I0(xy)f(x)d xF(y)=∫0∞yexp⁡[−12(y2+x2)]I0(xy)f(x)dxF(y) = \int_{0}^{\infty} y\exp{\left[-\frac{1}{2}(y^2 + x^2)\right]} I_0\left(xy\right)f(x)\;\mathrm{d}x Więc dla danego F(y)F(y)F(y) muszę aproksymować f( x...

Począwszy od równania doradczego w formie konserwatorskiej. ut=(a(x)u)xut=(a(x)u)x u_t = (a(x)u)_x gdzie a(x)a(x)a(x) jest prędkością zależną od przestrzeni, uuu jest stężeniem gatunku, który jest zachowany. Dyskretyzacja strumienia (gdzie strumień f=a(x)uf=a(x)uf=a(x)u jest zdefiniowany na...

W książce Nocedal & Wright o optymalizacji numerycznej znajduje się stwierdzenie w sekcji 2.2 (strona 27): „Ogólnie rzecz biorąc, łatwiej jest zachować niezmienność skali dla algorytmów wyszukiwania linii niż dla algorytmów regionu zaufania”. W tej samej sekcji mówią o posiadaniu nowych...

Jak mogę zastąpić metodę Eulera czwartym rzędem Runge-Kutty, aby określić ruch swobodnego spadku o nie stałej wielkości grawitacyjnej (np. Swobodny spadek z 10 000 km nad ziemią)? Do tej pory napisałem prostą integrację metodą Eulera: while() { v += getMagnitude(x) * dt; x += v * dt; time +=...