Pytania dotyczące tego, czy sekwencja iteracji generowana metodą iteracyjną ma jeden lub więcej punktów granicznych i czy te punkty graniczne mają prawidłowe właściwości.
Jak rozumiem, istnieją dwie główne kategorie iteracyjnych metod rozwiązywania liniowych układów równań: Metody stacjonarne (Jacobi, Gauss-Seidel, SOR, Multigrid) Metody podprzestrzeni Kryłowa (Gradient sprzężony, GMRES itp.) Rozumiem, że większość metod stacjonarnych działa poprzez iteracyjne relaksowanie (wygładzanie) trybów błędu Fouriera. Jak rozumiem, metoda gradientu sprzężonego (metoda podprzestrzeni Kryłowa) działa poprzez „przechodzenie” przez …
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 nieprawidłowych wyników, mimo że są spójne i stabilne. Na przykład w tym artykule pokazano, w …
Jaki jest teoretyczny współczynnik zbieżności dla narzędzia FFT Poison? Rozwiązuję równanie Poissona: przy n ( x , y , z ) = 3∇2)V.H.( x , y, z) = - 4 πn ( x , y, z)∇2)V.H.(x,y,z)=-4πn(x,y,z)\nabla^2 V_H(x, y, z) = -4\pi n(x, y, z) w domenie[0,2]×[0,2]×[0,2]z okresowymi warunkami brzegowymi. Ta …
tło Rozwiązuję wariant równania Ornsteina-Zernike z teorii płynów. Abstrakcyjnie problem można przedstawić jako rozwiązanie problemu punktu stałego , gdzie A jest operatorem alegbraicznym integro, a c ( r ) jest funkcją rozwiązania (funkcja korelacji bezpośredniej OZ). Rozwiązuję za pomocą iteracji Picarda, w której zapewniam wstępne rozwiązanie próbne c 0 ( …
Załóżmy, że mam następującą interesującą funkcję: Ma pewne nieprzyjemne właściwości, takie jak jego pochodna, która nie jest ciągła przy racjonalnych wielokrotnościach . Podejrzewam, że zamknięty formularz nie istnieje.πfa( x ) = ∑k ≥ 1sałatak xk2)( 2 - cosk x ).f(x)=∑k≥1coskxk2(2−coskx). f(x) = \sum_{k\geq1} \frac{\cos k x}{k^2(2-\cos kx)}. ππ\pi Mogę to …
Według Wikipedii współczynnik konwergencji wyraża się jako konkretny stosunek norm wektorowych. Próbuję zrozumieć różnicę między szybkościami „liniowymi” i „kwadratowymi” w różnych punktach czasu (w zasadzie „na początku” iteracji i „na końcu”). Czy można stwierdzić, że: ek+1ek+1e_{k+1}xk+1xk+1x_{k+1}∥ek∥‖ek‖\|e_k\| z kwadratową konwergencją norma błędu iteracji x_ {k + 1} jest ograniczona przez \ …
Próbuję rozwiązać następujący układ równań dla zmiennych i (wszystkie pozostałe są stałymi):x 2P., x1P.,x1P,x_1x2)x2)x_2 A ( 1 - P)2)- k1x1= 0A P.2)- k2)x2)= 0( 1 - P) ( r1+ x1)4L.1- P( r1+ x2))4L.2)= 0ZA(1-P.)2)-k1x1=0ZAP.2)-k2)x2)=0(1-P.)(r1+x1)4L.1-P.(r1+x2))4L.2)=0\frac{A(1-P)}{2}-k_1x_1=0 \\ \frac{AP}{2}-k_2x_2=0 \\ \frac{(1-P)(r_1+x_1)^4}{L_1}-\frac{P(r_1+x_2)^4}{L_2}=0 Widzę, że mogę przekształcić ten układ równań w pojedyncze równanie jednej zmiennej …
W informatyce często spotykamy duże układy liniowe, które musimy rozwiązać za pomocą niektórych (skutecznych) środków, np. Metod bezpośrednich lub iteracyjnych. Jeśli skupimy się na tym drugim, jak możemy ustalić, że iteracyjna metoda rozwiązywania dużych układów liniowych jest w praktyce zbieżna? Oczywiste jest, że możemy przeprowadzić analizę prób i błędów (por. …
Zacząłem uczyć się OpenFOAM z samouczka Cavity, który jest dostępny na stronie internetowej . Podczas eksperymentowania z różnymi liczbami Reynoldsa w sekcji „2.1.8.2 Uruchamianie kodu” samouczek mówi, aby ponownie uruchomić solver, ponieważ „Rozsądne jest zwiększenie czasu rozwiązania”. Ale kiedy to zrobiłem, nie mogłem znaleźć żadnej różnicy między przepływem we wnęce …
Chcę wiedzieć, które z rozwiązują klasyczny liniowych (np Gauss-Seidel Jacobiego, SOR) gwarantowane są zbieżne do problemu , gdzie jest pozytywne pół definitywna i oczywiścieA x = bZAx=bAx=bZAZAAb ∈ i m ( A )b∈jam(ZA)b \in im(A) (Uwaga jest półokreślona i nieokreślona)ZAZAA
W metodzie samodzielnego pola Hartree-Focka rozwiązywania niezależnego od czasu elektronicznego równania Schroedingera dążymy do zminimalizowania energii stanu podstawowego, , układu elektronów w polu zewnętrznym w odniesieniu do wyboru orbitali spinowych, { χ i } .E0E0E_{0}{χi}{χi}\{\chi_{i}\} Czynimy to poprzez iteracyjne rozwiązywanie 1 elektronów Hartree-Focka gdzie x i jest wirowanie / przestrzenne …
Metoda Newtona rozwiązywania równań nieliniowych jest znana, że zbiega się kwadratowo, gdy domysły początkowe są „wystarczająco blisko” do rozwiązania. Co jest „wystarczająco blisko”? Czy istnieje literatura na temat struktury tego basenu przyciągania?
Wiem, że częściowe przybliżenie liniowe elementu skończonego z spełnia pod warunkiem, że U jest wystarczająco gładki i f \ w L ^ 2 (U) .uhuhu_hΔ u ( x ) = f( x )w Uu ( x ) = 0na ∂UΔu(x)=fa(x)w Uu(x)=0na ∂U \Delta u(x)=f(x)\quad\text{in }U\\ u(x)=0\quad\text{on }\partial U ∥ u …
Rozważmy, że masz problem z nieskończoną przestrzenią w przestrzeni Hilberta lub Banacha (pomyśl o PDE lub problemie optymalizacji w takiej przestrzeni) i masz algorytm, który słabo zbiega się z rozwiązaniem. Jeśli dyskretyzujesz problem i zastosujesz odpowiedni dyskretny algorytm do problemu, wówczas słaba zbieżność to zbieżność w każdej współrzędnej, a zatem …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.