Sia f : A ⊆ R n → R f: A \subseteq \R^n \to \R f : A ⊆ R n → R e assumiamo che
∀ i = 1 , … , n , ∃ ∂ f ∂ x i \forall\ i = 1,\ldots,n,\ \exists\ \frac{\partial f}{\partial x_i} ∀ i = 1 , … , n , ∃ ∂ x i ∂ f . La
derivata seconda rispetto alla direzione j j j -esima della derivata i i i -esima è
definita come:
∂ 2 f ∂ x j ∂ x i ( x ) = ∂ ∂ x j ( ∂ f ∂ x i ) ( x ) \frac{\partial^2 f}{\partial x_j\ \partial x_i}(\mathbf{x}) = \frac{\partial}{\partial x_j} \left(\frac{\partial f}{\partial x_i}\right)(\mathbf{x}) ∂ x j ∂ x i ∂ 2 f ( x ) = ∂ x j ∂ ( ∂ x i ∂ f ) ( x )
Se j = i j = i j = i , la derivata parziale seconda è denominata pura ; altrimenti, è
detta mista .
La matrice Hessiana di f : R n → R f: \R^n \to \R f : R n → R è la matrice n × n n \times n n × n definita da:
D f 2 ( x ) = [ ∂ 2 f ∂ x 1 2 ( x ) … ∂ 2 f ∂ x 1 ∂ x n ( x ) ⋮ ⋮ ∂ 2 f ∂ x n ∂ x 1 ( x ) … ∂ 2 f ∂ x n 2 ( x ) ] D_f^2(\mathbf{x}) = \begin{bmatrix}
\frac{\partial^2 f}{\partial x_1^2}(\mathbf{x}) & \ldots & \frac{\partial^2 f}{\partial x_1\ \partial x_n}(\mathbf{x}) \\
\vdots & & \vdots \\
\frac{\partial^2 f}{\partial x_n\ \partial x_1}(\mathbf{x}) & \ldots & \frac{\partial^2 f}{\partial x_n^2}(\mathbf{x})
\end{bmatrix} D f 2 ( x ) = ∂ x 1 2 ∂ 2 f ( x ) ⋮ ∂ x n ∂ x 1 ∂ 2 f ( x ) … … ∂ x 1 ∂ x n ∂ 2 f ( x ) ⋮ ∂ x n 2 ∂ 2 f ( x )
In generale, le derivate parziali seconde miste opposte tra loro (rispetto alla
diagonale dell’Hessiana) possono essere diverse.
∂ 2 f ∂ x j ∂ x i ≠ ∂ 2 f ∂ x i ∂ x j \frac{\partial^2 f}{\partial x_j\ \partial x_i} \neq \frac{\partial^2 f}{\partial x_i\ \partial x_j} ∂ x j ∂ x i ∂ 2 f = ∂ x i ∂ x j ∂ 2 f
L’uguaglianza tra queste derivate in un punto x \mathbf{x} x è garantita se le due
derivate parziali seconde sono continue in un intorno di x \mathbf{x} x :
∂ 2 f ∂ x j ∂ x i , ∂ 2 f ∂ x i ∂ x j ∈ C 0 ( x ) ⟹ ∂ 2 f ∂ x j ∂ x i ( x ) = ∂ 2 f ∂ x i ∂ x j ( x ) \frac{\partial^2 f}{\partial x_j\ \partial x_i}, \frac{\partial^2 f}{\partial x_i\ \partial x_j} \in C^0(\mathbf{x})
\implies \frac{\partial^2 f}{\partial x_j\ \partial x_i}(\mathbf{x}) = \frac{\partial^2 f}{\partial x_i\ \partial x_j}(\mathbf{x}) ∂ x j ∂ x i ∂ 2 f , ∂ x i ∂ x j ∂ 2 f ∈ C 0 ( x ) ⟹ ∂ x j ∂ x i ∂ 2 f ( x ) = ∂ x i ∂ x j ∂ 2 f ( x )
Si dice che f f f sia di classe C m ( A ) C^m(A) C m ( A ) se tutte le sue derivate parziali
m m m -esime sono continue e f ∈ C m − 1 ( A ) f \in C^{m - 1}(A) f ∈ C m − 1 ( A ) .
Problema : Data una funzione f : B ( p 0 , r ) ⊆ R 2 → R f: B(\mathbf{p}_0, r) \subseteq \R^2 \to \R f : B ( p 0 , r ) ⊆ R 2 → R ,
di classe C m ( B ( p 0 , r ) ) C^m(B(\mathbf{p}_0, r)) C m ( B ( p 0 , r )) ; come la si può approssimare con un
polinomio di 2 variabili di ordine m m m nel modo migliore possibile?
Preso un punto p 0 = ( x 0 , y 0 ) ∈ R 2 \mathbf{p}_0 = (x_0, y_0) \in \R^2 p 0 = ( x 0 , y 0 ) ∈ R 2 , si definisce come polinomio
di ordine m m m di 2 variabili, centrato in p 0 \mathbf{p}_0 p 0 , una funzione
T : R 2 → R T: \R^2 \to \R T : R 2 → R del tipo:
T m , p 0 ( x , y ) = ∑ h = 0 m ∑ i = 0 h c i , h − i ( x − x 0 ) i ( y − y 0 ) h − i T_{m, \mathbf{p}_0}(x, y) = \sum_{h = 0}^m \sum_{i = 0}^h c_{i, h-i}\ (x - x_0)^i\ (y - y_0)^{h - i} T m , p 0 ( x , y ) = h = 0 ∑ m i = 0 ∑ h c i , h − i ( x − x 0 ) i ( y − y 0 ) h − i
dove i coefficienti c i , h − i ∈ R c_{i, h - i} \in \R c i , h − i ∈ R sono ( m + 1 ) ( m + 2 ) 2 \frac{(m + 1)(m + 2)}{2} 2 ( m + 1 ) ( m + 2 ) in
totale.
L’approssimazione di f ( p ) f(\mathbf{p}) f ( p ) in p 0 \mathbf{p}_0 p 0 è data dal valore del
polinomio più un errore di ordine superiore:
f ( p ) = T m , p 0 ( p ) + o ( ∣ ∣ p − p 0 ∣ ∣ m ) f(\mathbf{p}) = T_{m, \mathbf{p}_0}(\mathbf{p}) + o(||\mathbf{p} - \mathbf{p}_0||^m) f ( p ) = T m , p 0 ( p ) + o ( ∣∣ p − p 0 ∣ ∣ m )
Presi p 0 ∈ R 2 \mathbf{p}_0 \in \R^2 p 0 ∈ R 2 e f ∈ C 2 ( B ( p 0 , r ) ) f \in C^2(B(\mathbf{p}_0, r)) f ∈ C 2 ( B ( p 0 , r )) :
T 2 , p 0 ( p ) = f ( p 0 ) + ∇ f ( p 0 ) ⋅ ( p − p 0 ) + 1 2 ( D f 2 ( p 0 ) ( p − p 0 ) ) ⋅ ( p − p 0 ) T_{2, \mathbf{p}_0}(\mathbf{p}) = f(\mathbf{p}_0)
+ \nabla f(\mathbf{p}_0) \cdot (\mathbf{p} - \mathbf{p}_0)
+ \frac{1}{2} (D_f^2(\mathbf{p}_0)\ (\mathbf{p} - \mathbf{p}_0)) \cdot (\mathbf{p} - \mathbf{p}_0) T 2 , p 0 ( p ) = f ( p 0 ) + ∇ f ( p 0 ) ⋅ ( p − p 0 ) + 2 1 ( D f 2 ( p 0 ) ( p − p 0 )) ⋅ ( p − p 0 )