next up previous
Next: Variační metody Up: Okrajová úloha Previous: Metoda střelby

Metoda sítí (konečných diferencí)

Položme $x_0 \equiv a$ a $x_{N+1} \equiv b$. Vložíme mezi $a$ a $b$ body $x_1, \dots, x_n$. Budeme hledat aproximaci řešení v uvedených bodech. Nejjednodušší je ekvidistantní krok $\Delta x_k = x_{k+1} - x_k = (b - a)/(N + 1) = h$. Derivace lze nahradit konečnými diferencemi různě. Například

\begin{displaymath}
v'(x) \approx \frac{v (x+h) - v(x)}{h}
\end{displaymath}

a protože platí

\begin{displaymath}
v'(x) - \left[ \frac{v(x + h) - v(x)}{h} \right] = \frac{1}{2}
v''(x) h + O(h^2) ,
\end{displaymath}

máme metodu prvního řádu přesnosti. Pokud ovšem nahradíme derivaci vztahem

\begin{displaymath}
v'(x) \approx \frac{v (x+h) - v(x-h)}{2h},
\end{displaymath}

dostaneme metodu druhého řádu, protože platí

\begin{displaymath}
v'(x) - \left[ \frac{v(x + h) - v(x - h)}{2h} \right] = -
\frac{h^2}{3!} v'''(x) h + O(h^4).
\end{displaymath}

Pro druhou derivaci platí vztah druhého řádu přesnosti

\begin{displaymath}
v''(x) = \frac{v(x+h) - 2 v(x) + v(x-h)}{h^2} - \frac{h^2}{12}
v^{(4)}(x) + O(h^4).
\end{displaymath}

Ukázka pro lineární diferenciální rovnici

\begin{displaymath}\!\!\!\!\!\!\!\!
a(x) v'' + b(x) v' + c(x) v~= d(x) \qquad x \in
\langle 0, 1 \rangle \qquad v(0) = \alpha \quad v(1) = \beta
\end{displaymath}

Tuto rovnici lze pro všechna $i = 1, \dots,
N$ aproximovat diferenční rovnicí

\begin{displaymath}
a_i \frac{v_{i+1} - 2 v_i + v_{i-1}}{h^2} + b_i \frac{v_{i+1} -
v_{i-1}}{2h} + c_i v_i = d_i.
\end{displaymath}

Tuto rovnici můžeme pro $i = 1, \dots,
N$ upravit na tvar

\begin{displaymath}
\underbrace{(- a_i + b_i \frac{h}{2})}_{r_i} v_{i-1} +
\unde...
...nderbrace{(-a_i - b_i
\frac{h}{2})}_{q_i} v_{i+1} = - h^2 d_i.
\end{displaymath}

Potom řešíme soustavu lineárních rovnic pro $v_i$

\begin{displaymath}
\left[
\begin{array}{ccccccc}
p_1 & q_1 & 0 & & \dots & & 0 ...
...\\ d_{n-1} \\
d_n + q_n \frac{\beta}{h^2}
\end{array}\right].
\end{displaymath}

Lze dokázat, že pokud $h \to 0$, pak pro $\forall \ i = 1, \dots, N$ platí $\vert v_i - v(x_i)\vert \to 0$.

Pozn. Pokud okrajové podmínky obsahují derivace, musíme je též aproximovat. Ke zvýšení řádu této aproximace často užíváme virtuálních bodů $x_{-1}$ a $x_{N+2}$.

Pozn. Při řešení okrajových úloh pro obyčejné diferenciální rovnice metodou sítí je nutno řešit systémy lineárních rovnic s pásovou maticí.


next up previous
Next: Variační metody Up: Okrajová úloha Previous: Metoda střelby
Jiri Limpouch
2000-05-25