Subjects numerical analysis

Diff Finite Scheme C91D13

Step-by-step solutions with LaTeX - clean, fast, and student-friendly.

Search Solutions

Diff Finite Scheme C91D13


1. **Énoncé du problème** : On considère l'équation différentielle $$-u''(x) + c(x)u(x) = f(x), \quad x \in [0,1],$$ avec les conditions aux limites $$u(0) = 0$$ et $$u(1) = 0$$, où $$c(x) > 0$$. 2. **Schéma aux différences finies** : On discrétise l'intervalle $$[0,1]$$ en $$N+1$$ points avec un pas $$h = \frac{1}{N+1}$$. Pour $$i = 1, \ldots, N$$, on approxime $$u''(x_i)$$ par $$u''(x_i) \approx \frac{u_{i-1} - 2u_i + u_{i+1}}{h^2}$$. Le schéma devient donc : $$-\frac{u_{i-1} - 2u_i + u_{i+1}}{h^2} + c_i u_i = f_i,$$ avec $$c_i = c(x_i)$$ et $$f_i = f(x_i)$$. 3. **Matrice du système linéaire** : On obtient un système $$A \mathbf{u} = \mathbf{f}$$ où $$A$$ est une matrice tridiagonale définie par : - Diagonale principale : $$\frac{2}{h^2} + c_i > 0$$ - Diagonales secondaires : $$-\frac{1}{h^2}$$ 4. **Inversibilité de la matrice** : La matrice $$A$$ est strictement diagonale dominante et symétrique définie positive car $$c(x) > 0$$. Donc, $$A$$ est inversible. 5. **Consistance** : Le schéma est une approximation d'ordre 2 de la dérivée seconde, donc l'erreur locale est $$O(h^2)$$. 6. **Stabilité** : La matrice $$A$$ est définie positive, donc le système est stable. 7. **Convergence** : Par le théorème de Lax, la consistance et la stabilité impliquent la convergence du schéma. 8. **Ordre du schéma** : L'ordre global de convergence est $$O(h^2)$$. **Réponse finale** : Le schéma aux différences finies est donné par $$-\frac{u_{i-1} - 2u_i + u_{i+1}}{h^2} + c_i u_i = f_i,$$ avec $$A$$ matrice tridiagonale définie positive et inversible. Le schéma est consistant, stable, convergent et d'ordre 2.