Introduzione ai metodi matematici e numerici
La modellistica di un fenomeno si suddivide in quattro fasi:
- Analisi del fenomeno: attraverso semplificazioni, analogie ed astrazioni si effettua un'analisi quantitativa.
- Analisi del modello: si ricerca esistenza, regolarità e unicità e quindi si effettua un'analisi qualitativa.
- Individuazione del metodo
- Implementazione
Analisi dell'errore
Quando i risultati di un algoritmo possono considerarsi attendibili allora l'algoritmo si dice stabile. Un problema si dice problema mal condizionato se piccole variazioni sui dati inducono grandi variazioni sul risultato.
Complessità computazionale: numero di operazioni elementari che caratterizzano l'algoritmo. Per la risoluzione di un problema si possono applicare algoritmi diversi più o meno complessi.
Ad esempio, per il calcolo del polinomio P(x) = a0 + a1x + a2x2 + ... + amxm in un punto x, si possono seguire i seguenti algoritmi:
Algoritmo 1
- P = a0
- y = x
- PER i = 1, ..., m-1
- P = P + aiy
- y = yx
- P = P + amy
Algoritmo 2 (Ruffini-Horner)
- P = am
- PER i = 1, ..., m
- P = P x + am-i
Rappresentazioni Floating Point
Con x ∈ ℝ x = sgn(x) 0.d1d2...dt ● βp = 1.mβp dove β = base, p = esponente, m = d1d2...dt mantissa, 0 ≤ di ≤ β - 1.
L'unicità della rappresentazione è garantita da d1 ≠ 0, β = in rappresentazione normalizzata, 1/β ≤ m < 1.
32 Bit (Singola Precisione)
- | s | - 8 p - | 23 m |
- -127 ≤ P ≤ 128
- P ∈ [L, U]
64 Bit (Doppia Precisione)
- | s | - 11 p - | 52 m |
- -1023 ≤ P ≤ 1024
- P ∈ [L, U]
Numeri di macchina
Sono tutti i numeri rappresentabili esattamente in floating point. ℑ(β, t, L, U) = { x ∈ ℝ | x = sgn(x) 0, d1...dt βe ∧ 0 ≤ di ≤ β-1 ∧ L ≤ e ≤ U }
Se x ∉ ℑ bisogna effettuare un'approssimazione per cui ẋ = m·βe ∈ ℑ
Troncamento (Chopping)
m = 0.d1d2...dtdt+1...ds ⇒ ẋ = m(m) = 0.d1d2...dt
Errore assoluto EA = ẋ - x
Errore relativo ER = ẋ - x/x
Esempio
| m | ẋ | |ẋ - m| |
|---|---|---|
| 0.14 | 0.1 | 0.04 |
| 0.12 | 0.1 | 0.02 |
| 0.15 | 0.1 | 0.05 |
| 0.19 | 0.1 | 0.09 |
|EA| = |ẋ - x| = |ẋβ - mβe| = |ẋ - m| βe 1-t βe = βe+t
|ER| = |ẋ - x|/|ẋ| = |ẋβ - mβe|/|mβe| = |ẋ - m|/|m| 1-t
POICHÉ 1/β ≤ m 1/m ≤ β
Arrotondamento (Rounding)
(x) = βpm̃ = (m) se dt+1 (m + 1/2 βt) se dt+1 ≥ β/2
Esempio
| m | m̃ | |m̃ - m| |
|---|---|---|
| 0.11 | 0.1 | 0.01 |
| 0.12 | 0.1 | 0.02 |
| 0.15 | 0.2 | 0.05 |
Altri argomenti
- Analisi dell'errore
- Sistemi lineari
- Ricerca degli autovalori
- Approssimazioni di funzioni per interpolazioni
- Spline
- B-Spline
- Minimi quadrati
- Approssimazione di funzioni nel continuo
- Approssimazione trigonometrica
- Integrazione numerica
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
Scarica il documento per vederlo tutto.
-
Metodi numerici e modelli matematici - 139 di 139
-
Metodi matematici e numerici - parte 2
-
Metodi matematici e numerici - parte 3
-
Metodi matematici e numerici - parte 4