Prova scritta metodi matematici per l’ingegneria 17/01/2013
Modulo di calcolo numerico
D. Durante
Compito 1
- L’errore di troncamento dello schema è:
2 2 3 4 ! Δt ∂ u Δx ∂ u ν ∂ u 2 4 − − τ (Δt, Δx) = + a + O(Δt2, Δx4) 2 3 4 2 ∂t 6 ∂x 2 ∂x
Lo studio della stabilità comporta la soluzione della disequazione:
1|G| ≤ 1−1 + 2β(1 cos θ) + i c sin θ 2 ≥ − ≥ ≥
con c = aΔt/Δx e β = νΔt/Δx. Osservando che 2 1 cos θ 0 e β 0 per definizione, e che la parte immaginaria è sempre minore dell’unità, la disequazione risulta soddisfatta per ogni θ senza particolari limitazioni, ovvero lo schema risulta incondizionatamente stabile.
- | ∈ La serie di Fourier di f (x) = sin(x)| per x [0, π] si scrive:
∞ 2 4 X f (x) = + cos(2kx) 2−π π(1 4k2) k=1
Ovviamente in x = 0 la serie converge come 1/k2, dal momento che la funzione è pari ma la sua derivata prima è discontinua in x = kπ per k = 0, . . ..