Kiedy chce się obliczyć pochodne numeryczne, metoda przedstawiona przez Bengt Fornberg tutaj (i zgłaszane tutaj ) jest bardzo wygodne (zarówno precyzyjne i proste do wykonania). Jako, że oryginalny artykuł pochodzi z 1988 roku, chciałbym wiedzieć, czy istnieje dziś lepsza alternatywa (jako (lub prawie tak) prosta i bardziej precyzyjna)?
Mam pytanie dotyczące kodowania warunków brzegowych dla mechaniki bryły (elastyczność liniowa). W szczególnym przypadku muszę użyć różnic skończonych (3D). Jestem bardzo nowy w tym temacie, więc być może niektóre z poniższych pytań mogą być bardzo podstawowe. Aby doprowadzić do mojego konkretnego problemu, przede wszystkim chcę pokazać, co już zaimplementowałem (aby …
Próbuję rozwiązać równanie typu: (−∂2∂x2−f(x))ψ(x)=λψ(x)(−∂2∂x2−f(x))ψ(x)=λψ(x) \left( -\tfrac{\partial^2}{\partial x^2} - f\left(x\right) \right) \psi(x) = \lambda \psi(x) Gdzie f(x)f(x)f(x) ma prosty biegun o wartości 000 , dla najmniejszych NNN wartości własnych i wektorów własnych. Warunki brzegowe to: ψ(0)=0ψ(0)=0\psi(0) = 0 i ψ(R)=0ψ(R)=0\psi(R)=0 , a ja patrzę tylko na funkcję ponad (0,R](0,R](0,R] . …
Czy ktoś może mi pomóc znaleźć książki na temat rozwiązań numerycznych (różnic skończonych i metod Cranka – Nicolsona) równań Poissona i dyfuzji, w tym przykładów na nieregularnej geometrii, takich jak dziedzina składająca się z obszaru między prostokątem a okręgiem (zwłaszcza książek lub linków na przykładach kodu MATLAB w tym przypadku)?
Wdrażam artykuł „ Optymalny transport masowy do rejestracji i wypaczania ”, moim celem jest umieszczenie go w Internecie, ponieważ po prostu nie mogę znaleźć żadnego kodu eulerowskiego transportu masowego w Internecie, co byłoby interesujące, przynajmniej dla społeczności naukowej zajmującej się przetwarzaniem obrazu. Artykuł można podsumować w następujący sposób: - znajdź …
Używam schematu różnic skończonych Cranka-Nicolsona do rozwiązania równania cieplnego 1D. Zastanawiam się, czy zasada maksimum / minimum równania ciepła (tj. Że maksimum / minimum występuje w stanie początkowym lub na granicach) również obowiązuje dla rozwiązania dyskretnego. Prawdopodobnie wynika to z faktu, że Crank-Nicolson jest stabilnym i zbieżnym schematem. Wygląda jednak …
Rozważ następujący problem Wuv=FWuv=F W_{uv} = F którym warunek wymuszający może zależeć od ( formuła znajduje się poniżej w Edycji 1 ) oraz i jego pierwszych pochodnych. Jest to równanie fali wymiarowej 1 + 1. Mamy wstępne dane określone w .u,vu,vu,vWWW{u+v=0}{u+v=0}\{u+v = 0\} Interesuje mnie rozwiązanie w dziedzinie zależności przedziału …
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 „zamrożenie” nieliniowego terminu, tzn. termin , gdzie „uważa się za reprezentujące …
Tło: Zbudowałem tylko jedno działające rozwiązanie numeryczne dla 2d Navier-Stokes na kurs. Było to rozwiązanie dla przepływu wnęki napędzanej pokrywką. Kurs omawiał jednak szereg schematów dyskretyzacji przestrzennych i dyskretyzacji czasowych. Wziąłem też więcej zajęć z manipulowaniem symbolami w stosunku do NS. Niektóre z podejść numerycznych do obsługi konwersji równania analitycznego …
Stosuję metodę różnic skończonych do układu 3 równań sprzężonych. Dwa równania nie są sprzężone, jednak trzecie równanie łączy się z dwoma pozostałymi. Zauważyłem, że zmieniając porządek równań, powiedzmy od do że macierz współczynników staje się symetryczna.( x , y, z)(x,y,z)(x, y, z)( x , z, y)(x,z,y)(x, z, y) Czy jest …
Jak aproksymować liczbę warunków dużej matrycy GGG, gdyby GGG jest kombinacją transformacji Fouriera FFF (niejednorodne lub jednolite), różnice skończone RRRi macierze diagonalne SSS? Matryce są bardzo duże i nie są przechowywane w pamięci i są dostępne tylko jako funkcje. W szczególności mam następującą macierz: Gμ=SHFHFS+μRHRGμ=SHFHFS+μRHRG_\mu=S^HF^HFS+\mu R^HR Chcę zbadać związek między …
Gdy preferowane jest użycie wielomianów Bernsteina do przybliżenia funkcji ciągłej zamiast stosowania tylko następujących wstępnych metod analizy numerycznej: „Wielomiany Lagrange'a”, „Proste operatory różnic skończonych”. Pytanie dotyczy porównania tych metod.
Załóżmy, że miałem następujący okresowy problem z poradą 1D: ∂u∂t+c∂u∂x=0∂u∂t+do∂u∂x=0\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0 w Ω=[0,1]Ω=[0,1]\Omega=[0,1] u(0,t)=u(1,t)u(0,t)=u(1,t)u(0,t)=u(1,t) u(x,0)=g(x)u(x,0)=g(x)u(x,0)=g(x) gdzie ma nieciągłość skoku przy . g(x)g(x)g(x)x∗∈(0,1)x∗∈(0,1)x^*\in (0,1) Rozumiem, że w przypadku schematów liniowych różnic skończonych wyższego niż pierwszego rzędu, fałszywe oscylacje występują w pobliżu nieciągłości, jak jest …
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.