16 09 15

Page 1

Metodi Analitici e Numerici per l’Ingegneria Quinto appello 16 Settembre 2015 DOMANDA 1 Si consideri il seguente problema:    −u′′(x) + 2u(x) = sinx x ∈ (0,π ) u(0) = 0 u′(π) = 0. 1. (T) Si scriva la formulazione debole del problema, indicando chiaramente quali sono gli spazi funzionali e con quale motivazione vengono scelti. 2. (E) Si dimostri la coercivit` a della forma bilineare. 3. (T) Si disegnino le funzioni di base dello spazio di elementi finiti lineari definiti su una partizione unifome dell’intervallo [0 ,π ] di passo h = 1/3, che approssimano lo spazio funzionale scelto per la formulazione debole. 4. (T) Si scriva l’approssimazione di Galerkin del problema, Au = f ottenuta utilizzando lo spazio degli elementi finiti lineari descritto al punto 3, indicando come gli elementi della matrice A e del vettore f dipendono dagli elementi della base a elementi finiti. DOMANDA 2 Si consideri il seguente sistema di equazioni lineari Ax =b con A =   4 1 0 0 −1 4 1 0 0 −1 4 1 0 0 −1 4   e b = [1, 5, 3,−1]T . 1. (E) Si verifichi una condizione sufficiente per la fattorizzazione LU (senza pivoting). 2. (M) Si calcoli la fattorizzazione LU mediante opportuno comando Matlab. 3. (M) Si implementino gli algoritmi di sostituzione in avanti e indietro. 4. (M) Si risolva il sistema. DOMANDA 3 Si consideri il problema di Cauchy { y′(t) = 3y(t) + sint y(0) = 1 1. (E) Si risolva il problema utilizzando la trasformata di Laplace. 2. (T) Si ricavi l’algoritmo di Eulero all’indietro (ossia si scriva esplicitamente quale approssi- mazione si utilizza) prima nel caso generale di un problema di Cauchy del primo ordine, poi applicato a questo problema (giustificando i passaggi). In particolare, si sfrutti la linearit` a del problema per ricavare una formula esplicita per un+1 in funzione…